# Soil temperature effects on the structure and diversity of plant and invertebrate communities in a natural warming experiment

Sinikka I. Robinson, Orla B. Mclaughlin, Bryndis Marteinsdottir, Eoin J. O'Gorman

<figure><img src="image_1.png" /></figure>

## To cite this version:

Sinikka I. Robinson, Orla B. Mclaughlin, Bryndis Marteinsdottir, Eoin J. O'Gorman. Soil temperature effects on the structure and diversity of plant and invertebrate communities in a natural warming experiment. Journal of Animal Ecology, 2018, 87 (3), pp.634-646. 10.1111/1365-2656.12798. hal-02625664

HAL Id: hal-02625664

https://hal.inrae.fr/hal-02625664

Submitted on 26 May 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

Journal of Animal Ecology

<figure><img src="image_3.png" /></figure>

BRITISH
ECOLOGICAL
SOCIETY

# Soil temperature effects on the structure and diversity of plant and invertebrate communities in a natural warming experiment

Sinikka I. Robinson<fnref n="1" />,<fnref n="2" /> ID
Eoin J. O'Gorman<fnref n="1" /> ID

Órla B. McLaughlin<fnref n="3" />

Bryndís Marteinsdóttir<fnref n="4" />,<fnref n="5" /> ID

**Correspondence**

Eoin J. O'Gorman
Email: e.ogorman@imperial.ac.uk

**Funding information**

NERC, Grant/Award Number: NE/ L011840/1 and NE/M020843/1; British Ecological Society, Grant/Award Number: 4009-4884 and SR16/1152; Royal Society, Grant/Award Number: RG140601; Finnish Cultural Foundation

Handling Editor: Mariano Rodriguez-Cabal

**Abstract**

1. Global warming is predicted to significantly alter species physiology, biotic interactions and thus ecosystem functioning, as a consequence of coexisting species exhibiting a wide range of thermal sensitivities. There is, however, a dearth of research examining warming impacts on natural communities.

2. Here, we used a natural warming experiment in Iceland to investigate the changes in above-ground terrestrial plant and invertebrate communities along a soil temperature gradient (10°C-30°C).

3. The α-diversity of plants and invertebrates decreased with increasing soil temperature, driven by decreasing plant species richness and increasing dominance of certain invertebrate species in warmer habitats. There was also greater species turnover in both plant and invertebrate communities with increasing pairwise temperature difference between sites. There was no effect of temperature on percentage cover of vegetation at the community level, driven by contrasting effects at the population level.

4. There was a reduction in the mean body mass and an increase in the total abundance of the invertebrate community, resulting in no overall change in community biomass. There were contrasting effects of temperature on the population abundance of various invertebrate species, which could be explained by differential thermal tolerances and metabolic requirements, or may have been mediated by changes in plant community composition.

5. Our study provides an important baseline from which the effect of changing environmental conditions on terrestrial communities can be tracked. It also contributes to our understanding of why community-level studies of warming impacts are imperative if we are to disentangle the contrasting thermal responses of individual populations.

**KEYWORDS**

Arctic, climate change, Hengill, invertebrate community, natural experiment, pitfall, soil temperature, vegetation

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

© 2018 The Authors. *Journal of Animal Ecology* published by John Wiley & Sons Ltd on behalf of British Ecological Society.

<footnote n="1">¹Department of Life Sciences, Imperial College London, Ascot, UK</footnote>

<footnote n="2">²Faculty of Biological and Environmental Sciences, University of Helsinki, Lahti, Finland</footnote>

<footnote n="3">³Agroécologie, AgroSup Dijon, INRA, Université Bourgogne Franche-Comté, Dijon, France</footnote>

<footnote n="4">⁴Institute of Life and Environmental Sciences, University of Iceland, Reykjavík, Iceland</footnote>

<footnote n="5">⁵The Soil Conservation Service of Iceland, Hella, Iceland</footnote>
# 1 | INTRODUCTION

The average global surface temperature has increased by 0.8°C since 1880 and is predicted to rise by at least 1.5°C during the next century (IPCC, 2014). Evidence for the ecological impacts of global warming at the population level is already substantial (see Parmesan & Yohe, 2003; Walther, 2010 for review), for example altered geographical distributions (Chen, Hill, Ohlemuller, Roy, & Thomas, 2011), shifts in phenology (Menzel, Sparks, & Estrella, 2006) and decreasing body size (Daufresne, Lengfellner, & Sommer, 2009). Species do not exist in isolation, however, and complex networks of trophic interactions make it difficult to extrapolate these population-level impacts to the community or ecosystem levels (Ings et al., 2009; Walther, 2010; Woodward, Perkins, & Brown, 2010).

Understanding and quantifying the impacts of warming across multiple levels of biological organisation is important for modelling the ecological and evolutionary dynamics of ecosystem change (Montoya & Raffaelli, 2010). Climate change is already generating new communities and restructuring assemblages of species as a consequence of shifts in geographical distributions and/or adaptation to novel climatic conditions (Lurgi, López, & Montoya, 2012). It is likely that coexisting species will have a wide range of thermal sensitivities, as well as physiological and behavioural adaptations, inducing differential responses to warming and precipitating effects on the physiology of individuals, biotic interactions and ecosystem functioning. Warming should thus have effects across all levels of biological organisation, rooted in the relationship between temperature, metabolism and body mass (Brose et al., 2012; Brown, Gillooly, Allen, Savage, & West, 2004).

Cold-adapted populations are likely to decline dramatically in abundance or become locally extinct with warming (Hering et al., 2009; Somero, 2010; Thomas et al., 2004). This is expected to result in decreasing α-diversity, either through a decline in species richness and/or a decline in evenness, as more resilient populations with higher thermal optima begin to dominate the community (O'Gorman et al., 2012; Sharp et al., 2014). Such changes may be complicated if new, warm-adapted species invade to offset the loss of cold-adapted species (Krajick, 2004; Lejeusne, Chevaldonne, Pergent-Martini, Bourdouresque, & Perez, 2010; Walther et al., 2002). This could maintain overall species richness, but the implications for evenness of the community are more difficult to predict (Andrew & Hughes, 2005; Friberg et al., 2009; Hillebrand, Soininen, & Snoeijns, 2010; Woodward, Dybkjaer, et al., 2010). For example, the decline of rove beetles with warming of an agroecosystem was balanced by increasing dominance of ground beetles, resulting in no overall effect of temperature on species richness (Berthe, Derocles, Lunt, Kimball, & Evans, 2015).

Consumer responses to temperature will be driven by resource availability. There is evidence for increasing primary production as a consequence of warming which may be particularly apparent in colder, high latitude areas, where earlier snowmelt and warmer conditions stimulate growth (Klanderud & Totland, 2005; Rustad et al., 2001). This increased productivity may offset the predicted

reduction in plant biomass with warming due to stronger grazing pressure, thus supporting the higher metabolic demands of grazing species in a warmer environment. For example, warming and thus earlier snowmelt increased insect herbivore species richness and damage to several plant species in a long-term manipulation experiment (Roy, Güsewell, & Harte, 2004). Plant diversity and community composition are likely to change as a consequence of warming, however (Klein, Harte, & Zhao, 2004; Walker et al., 2006), with evidence for a shift in plant dominance hierarchies and species evenness as a consequence of different thermal tolerances (Klanderud & Totland, 2005). Such changes in vegetation could alter habitat complexity and thus the associated invertebrate community, with plant and insect diversity shown to be positively correlated (Muren, Hoffmann, & Kwak, 2003).

Warming impacts on community composition will be mediated by the physiology of populations. In particular, body size underlies community- and population-level responses to warming because it is directly shaped by temperature-dependent metabolic processes (Brose et al., 2012; Brown et al., 2004). Communities are generally predicted to have a smaller mean body size in warmer regions either because they are made up of smaller sized species, or because the body size of species' populations within the community are smaller (Daufresne et al., 2009). There is widespread empirical support for declining body mass with warming at the community and population levels (see Daufresne et al., 2009; Kingsolver & Huey, 2008 for reviews). The mechanisms underlying these patterns remain contentious (see Gardner, Peters, Kearney, Joseph, & Heinsohn, 2011; Klok & Harrison, 2013), with thermoregulation in endotherms (Porter & Kearney, 2009), "rate-size" trade-offs in ectotherms (DeLong, 2012) and competition for limiting nutrients in unicellular algae (Reuman, Holt, & Yvon-Durocher, 2014) among the explanations. These models are not universally accepted and exceptions exist, for example in Baltic phytoplankton (Rüger & Sommer, 2012), centipedes (Vedel, Chipman, Akam, & Arthur, 2008), grasshoppers (Walters & Hassall, 2006), stream invertebrates (O'Gorman et al., 2012) and freshwater diatoms (Adams et al., 2013). Nevertheless, warming is predicted to alter the size-structure of communities, with far-reaching impacts on biotic interactions and ecosystem functioning (Montoya & Raffaelli, 2010).

Warming is predicted to reduce the standing stock of primary producers and herbivores due to the greater metabolic demands of consumer species in warmer environments (O'Connor, Gilbert, & Brown, 2011), although this may be offset by greater primary production. Consumer losses may be further compounded by ingestion inefficiency due to the weaker temperature dependence of feeding relative to metabolism (Lemoine & Burkepile, 2012; Vucic-Pestic, Ehnes, Rall, & Brose, 2011). Larger organisms at the highest trophic levels are likely to be the most susceptible as a consequence of their greater metabolic demands and lower population densities (Brown et al., 2004; Petchey, McPhearson, Casey, & Morin, 1999). Loss of the largest species and the increased prevalence of smaller organisms at lower trophic levels should lead to an overall decline in community biomass (Yvon-Durocher et al., 2011). The heterogeneity of
species thermal responses or physiological characteristics that determine susceptibility to warming could, however, lead to variable effects on individual populations (Dollery, Hodkinson, & Jónsdóttir, 2006; Dong, Hou, Ouyang, & Zhang, 2013; O'Gorman et al., 2012).

Here, we make use of a high latitude natural soil warming experiment to investigate the impact of environmental temperature on the diversity, size-structure, and biomass of terrestrial plant and invertebrate communities. We tested the following four hypotheses: (1) α-diversity of plants and invertebrates decreases with increasing soil temperature; (2) species turnover of plants and invertebrates increases with increasing pairwise temperature difference between sites; (3) the mean body mass of invertebrates declines with increasing soil temperature at the community- and population levels; and (4) percentage cover of vegetation and the total abundance and biomass of invertebrates decrease with increasing soil temperature at the community level, with variable effects at the population level.

## 2 | MATERIALS AND METHODS

### 2.1 | Study site description

Geothermal regions that are not confounded by extreme physical or chemical variables have been identified as ideal systems for studying the impacts of environmental temperature on naturally occurring communities (O'Gorman et al., 2014). The Hengill valley in Iceland

(64.03°N 21.18°W) contains a catchment of 16 geothermally heated streams that have been extensively studied and monitored over the past decade. Freshwater research at Hengill revealed a decline in macroinvertebrate evenness with increasing stream temperature and increased species turnover along the stream temperature gradient, but no effect of temperature on species richness (Friberg et al., 2009; Woodward, Perkins, et al., 2010). Studies have also shown mixed effects of stream temperature on population-level body mass (Adams et al., 2013; O'Gorman et al., 2012) and increased community-level abundance and biomass with increasing temperature (Hannesdottir, Gislason, Olafsson, Olafsson, & O'Gorman, 2013; O'Gorman et al., 2012, 2016).

The various stream temperatures arise from differential geo-thermal heating of the landscape, which also creates a soil temperature gradient. Here, we chose three locations on each bank of the 16 streams previously studied at Hengill for a total of 96 habitat patches. Note that we consider the three locations on each bank as pseudoreplicates and so we present mean values from these three habitat patches for all environmental and biotic data for a total of 32 independent sites (Figure 1). GPS coordinates and spot measurements of soil temperature at 10 cm depth were taken between 12.00 hr and 18.00 hr on the same day at each habitat patch during pilot work in August 2012 using a digital soil thermometer (Digitron FM15). Soil temperature was monitored at each habitat patch during the main study period in July 2013 using Maxim Integrated DS1921G

<figure><img src="image_4.png" /></figure>

**FIGURE 1** The 32 terrestrial sites in the Hengill geothermal region. Site names are annotated on the map in bold, consisting of a number (corresponding to the stream names in previous publications from the system) and a letter (corresponding to the bank of the stream: L = left; R = right). Each site name is followed by the mean soil temperature (°C) as calculated from 48-hr temperature logger data in July 2013. The "=" symbol, used at two points along the main river, indicates a much greater distance than shown on the map. Three habitat patches were selected along the stream bank at each terrestrial site. Note that habitat patches were established directly opposite each other on the left and right bank of each stream (see inset). Five pitfall traps were established at each habitat patch in a 1 m² area, located c. 30 cm from the stream bank, as shown in the inset
Thermocron iButton temperature loggers. One iButton was buried 10 cm below the soil surface at each site and the mean (±1 SD) temperature for each bank was calculated from measurements recorded every 10 min over a 48-hr time period, that is including both night and day (Table S1).

## 2.2 | Soil analyses

Soil moisture was measured at each habitat patch using a soil moisture probe (Imko GmbH TRIME-TDR) in July 2013. Five soil cores of 1.5 × 5 cm were taken from c. 10 cm below the surface using a soil auger during the same time period. The soil was homogenised, dried at 60°C for 96 hr and finely milled before analysis. Soil pH was measured from 10 g of the dry, milled soil, with a 1:5 ratio of soil to distilled water, using a pH probe (Dr.Meter® 0.01pH Resolution pH Meter). A further 5 g was tested for total carbon and total nitrogen content at the Forest Research Centre for Ecosystems, Society, and Biosecurity (Farnham, UK), using a combustion method with a Carlo Erba CN analyser (Flash 1112 series).

## 2.3 | Vegetation survey

A vegetation survey was conducted in July 2013 to determine the above-ground plant community composition. A 0.5 × 0.5 m quadrat was placed randomly on the ground at each habitat patch and photographed using a Nikon camera. The photographs were analysed according to the methods for sampling vascular plants suggested by Peet, Wenthworth, and White (1998), using cover classes to describe the percentage cover of each species: Class 1 (0%-1%), Class 2 (1%-2%), Class 3 (2%-5%), Class 4 (5%-10%), Class 5 (10%-25%), Class 6 (25%-50%), Class 7 (50%-75%), Class 8 (75%-95%) and Class 9 (95%-99%). The mid-points of the cover classes were used for further analysis of vegetation cover (e.g. if a plant species had a cover in Class 2 [1%-2%], we gave it the value 1.5%). All vascular plants, except for graminoids, were identified to species level with nomenclature following Kristinsson and Sigurdsson (2010), while the percentage cover of grasses, mosses, lichens and litter was also noted.

## 2.4 | Invertebrate sampling

Above-ground terrestrial invertebrate communities were sampled using pitfall traps in August 2012 and July 2013. Five pitfall traps were established at each habitat patch in a 1 m² area, with one trap at each of the four corners and one in the centre (for a more comprehensive trapping area). White plastic cups of 7 cm diameter and 8.5 cm depth were filled with 10 ml of ethylene glycol and 30 ml of stream water, and left for 48 hr before collection (after Woodcock, 2007). Ethylene glycol was used for pitfall trapping as it prevents the escape of invertebrates from the traps and acts as a preservative. During collection, samples from the five traps at each habitat patch were combined into a 250-µm sieve and stored in 70% ethanol. Note that a fine mesh size was used to prevent loss of small organisms such as mites and springtails. Terrestrial invertebrates were

identified to species level where possible. We will refer to all taxa as species from henceforth.

## 2.5 | Statistical analysis

All statistical analyses were carried out in R 3.2.0. Pearson correlations were used to test for relationships between temperature and pH, moisture, total carbon and total nitrogen (using the *cor* function in the *stats* package). Pairwise distances between sites were computed from the GPS coordinates, using the *earth.dist* function in the *fossil* package. A Mantel test was then used to test for spatial dependency of the temperature gradient, by comparing pairwise temperature difference to the pairwise distance between sites, using the *mantel* function in the *vegan* package.

The α-diversity of the plant and invertebrate assemblages was assessed as Shannon-Weiner diversity (*H'*, using the *diversity* function in the *vegan* package), species richness (*S*) and Pielou's evenness (*J'* = *H'*/In[*S*]) at each site. To avoid assumptions about linearity, the effect of temperature on the α-diversity of invertebrates and plants in July 2013 was explored with generalised additive models (GAM: *gam* function in the *mgcv* package, with *k* = 5 and *bs* = "cr", after Keele, 2008). Analyses incorporating the main and interactive effects of temperature and year on the α-diversity of invertebrates are presented in Supporting Information. Species turnover of the vegetation and invertebrate assemblages was quantified as the Sørensen similarity between pairwise combinations of all sites: $\beta = (2S_T/(S_1 + S_2)) - 1$, where $S_1$ and $S_2$ are the species richness for sites 1 and 2, respectively, and $S_T$ is the total species richness of the two sites. Note that greater similarity corresponds to lower species turnover between sites. A Mantel test was used to analyse the relationship between community similarity and pairwise temperature difference between sites for both plants and invertebrates.

Invertebrate population abundance in each habitat patch was estimated as the number of individuals of a species found in the five traps, with abundance at the site level taken as the mean value of the three habitat patches at each site. The total abundance of the invertebrate community was the sum of all population abundances per site. Individual body masses were estimated by measuring one linear dimension (body length) and converting to dry weight (mg) using length-weight relationships (Table S2, Figure S1). All individuals were measured apart from mites and springtails, where 30 individuals were measured for each habitat patch. From these data, the mean body mass of each species and the abundance-weighted arithmetic mean body mass (henceforth "mean body mass") of the invertebrate community were estimated for each site. Invertebrate population biomass was estimated as the abundance of each species per habitat patch multiplied by its mean body mass, while the total biomass of the invertebrate community was the sum of all population biomasses.

The effect of temperature on the percentage cover of the plant community and the mean body mass, total abundance and total biomass of the invertebrate community in July 2013 were explored with GAM. Analyses incorporating the main and interactive effects
of temperature and year on the mean body mass, total abundance and total biomass of the invertebrate community are presented in Supporting Information. For the population-level analyses, we explored the effect of temperature on the percentage cover of vegetation taxa and the mean body mass, abundance and biomass of each invertebrate species with GAM. We excluded species that were found at fewer than 10 sites to avoid biases associated with small sample size. All abundance, mean body mass and biomass data were transformed with $log_{10}(x)$, or $log_{10}(x + 1)$ if zeros were present, to meet the assumptions of normality and homogeneity of variance and given our a priori expectation that these metrics should approximate a log-normal distribution (Jonsson, Cohen, & Carpenter, 2005). The false discovery rate was used to avoid issues with multiple comparisons in the population-level analyses (Benjamini & Hochberg, 1995),

using the *p-adjust* function in the *stats* package, where *n* was the number of population-level comparisons for percentage cover of vegetation taxa or the mean body mass, abundance or biomass of invertebrate species.

To determine whether the plant community may be indirectly mediating changes in the invertebrate community, we repeated the GAM analyses of invertebrate community metrics with plant community metrics included as covariates with temperature. Here, either the species richness, evenness, diversity, mean body mass, total abundance or total biomass of invertebrates was the response variable, and either the species richness, evenness, diversity or percentage cover of plants was a covariate with soil temperature. Note that structural equation modelling could not be implemented as a tool to determine relationships between our measured variables because

<figure><img src="image_5.png" /></figure>

**FIGURE 2** Relationships between soil temperature and various metrics of α-diversity in July 2013: (a) plant species richness; (b) invertebrate species richness; (c) Pielou's evenness for plants; (d) Pielou's evenness for invertebrates; (e) Shannon diversity for plants; (f) Shannon diversity for invertebrates. Solid and dashed lines are the predicted fitting and 95% confidence intervals, respectively, from significant GAM models
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** Declining Sørensen similarity in community composition with increasing pairwise temperature difference between sites in July 2013: (a) plants (y = 0.756–0.006x); (b) invertebrates (y = 0.618–0.004x)</figcaption></figure>

of small sample size (i.e. 32 independent sites; Kline, 2015) and poor fittings between models and our data (i.e. root mean square error of approximation >0.2 in all analyses; MacCallum, Browne, & Sugawara, 1996).

All models were tested for spatial autocorrelation using the spline.correlog and mantel.test functions in the ncf package, with x coordinates corresponding to distance between sites along a North-South axis, y coordinates corresponding to distance between sites along a West-East axis, and model residuals as the observations at each site. There was little evidence for spatial autocorrelation of residuals in any analysis (Figures S2-S4), so no spatial structure was incorporated into any of the models.

## 3 | RESULTS

### 3.1 | Environmental variables

No significant correlation was found between temperature and any of the environmental variables measured in July 2013: pH (Pearson correlation: *r* = -.20, *t*<sub>28</sub> = -1.085, *p* = .287), moisture (Pearson correlation: *r* = -.03, *t*<sub>30</sub> = -0.166, *p* = .869), total carbon (Pearson correlation: *r* = -.32, *t*<sub>30</sub> = -1.860, *p* = .073) or total nitrogen (Pearson correlation: *r* = -.28, *t*<sub>30</sub> = -1.628, *p* = .114). Temperature was also not correlated with spatial distance between sites in either 2012 (Mantel test: *r* = .07, *p* = .256) or 2013 (Mantel test: *r* = .01, *p* = .373).

### 3.2 | Vegetation

There was a significant decrease in plant species richness with increasing soil temperature (GAM: *r*<sup>2</sup> = .12, *F* = 5.183, *p* = .030; Figure 2a), partially supporting our first hypothesis. There was, however, no significant relationship between soil temperature and the Pielou's evenness (GAM: *r*<sup>2</sup> = .04, *F* = 2.099, *p* = .155; Figure 2c) or Shannon diversity (GAM: *r*<sup>2</sup> = .07, *F* = 3.294, *p* = .080; Figure 2e) of the plant community at Hengill. There was a significant decrease in community similarity (i.e. increased species turnover) with increasing pairwise temperature difference between sites (Mantel test: *r* = .327, *p* = .006; Figure 3a), supporting our second hypothesis.

In contrast to our fourth hypothesis, there was no significant effect of temperature on the percentage cover of plants at the community level (GAM: *r*<sup>2</sup> < .01, *F* = 0.099, *p* = .755). There was limited support for our fourth hypothesis at the population level, with a significant decrease in the percentage cover of the Alpine bistort *Persicaria vivipara* with increasing soil temperature (Table S3; Figure S5a). There was an increase in the percentage cover of the willow-herb *Epilobium* sp. from 10°C to 20°C and a decrease in percentage cover at higher soil temperatures (Table S3; Figure S5b). There was also a significant increase in the percentage cover of lichens, the thymeleaf speedwell *Veronica serpyllifolia* and the marsh violet *Viola palustris* with increasing soil temperature (Table S3; Figure S5c-e), in direct contrast to our fourth hypothesis. There were no significant effects of temperature on the percentage cover of mosses, grasses, litter or any other plant species (Table S3).

### 3.3 | Invertebrate community

There was no effect of temperature on invertebrate species richness (GAM: *r*<sup>2</sup> < .01, *F* = 0.133, *p* = .718; Figure 2b), but there was a significant reduction in Pielou's evenness (GAM: *r*<sup>2</sup> = .38, *F* = 9.420, *p* < .001; Figure 2d) and Shannon diversity (GAM: *r*<sup>2</sup> = .34, *F* = 7.984, *p* = .002; Figure 2f) of the invertebrate community as soil temperature increased, partially supporting our first hypothesis. There was also a significant decrease in community similarity (i.e. increased species turnover) with increasing pairwise temperature difference between sites (Mantel *r* = .218, *p* = .022; Figure 3b), supporting our second hypothesis. Similar results were found in the August 2012 pilot study (Figure S6).

There was a significant decrease in the mean body mass of the invertebrate community with increasing soil temperature (GAM: *r*<sup>2</sup> = .19, *F* = 8.273, *p* = .007, Figure 4a), in line with our third hypothesis. A similar response was detected in August 2012 (Figure S7a). There was limited support for our third hypothesis at the population level, with a significant decrease in the mean body mass of parasite wasps (Platygastridae spp.; Table S4; Figure S8a). However, there was an increase in the mean body mass of springtails (Collembola spp.) and true bug nymphs (Hemiptera spp.) from 10°C to 15°C
<figure><img src="image_7.png" /><figcaption>**FIGURE 4** Relationships between soil temperature and (a) mean body mass; (b) total abundance; (c) total biomass of the invertebrate community in July 2013. Solid and dashed lines are the predicted fitting and 95% confidence intervals, respectively, from significant GAM models</figcaption></figure>

before declining at higher temperatures (Table S4; Figure S8b,c) and no significant effects on the remaining species (Table S4).

In contrast to our fourth hypothesis, there was a significant increase in the total abundance of terrestrial invertebrates with increasing soil temperature (GAM: $r^2$ = .22, *F* = 5.879, *p* = .014, Figure 4b). A similar response was detected in August 2012 (Figure S7b). There were contrasting effects of temperature on population-level abundances, partially supporting our fourth hypothesis. Here, adult flies (Diptera spp.), the harvestman *Mitopus morio*, the wolf

spider *Pardosa palustris* and the carabid beetle *Patrobus septentrionis* decreased in abundance with increasing soil temperature (Table S4; Figure S9a-d). Conversely, earthworms (Annelida spp.), springtails, and the carabid beetles *Pterostichus diligens* and *Pterostichus nigrita* increased in abundance with increasing soil temperature (Table S4; Figure S9e-h). There were no significant changes in population abundance of the remaining invertebrate species across the soil temperature gradient (Table S4).

In contrast to our fourth hypothesis, there was no effect of temperature on the total biomass of terrestrial invertebrates with increasing soil temperature in July 2013 (GAM: $r^2$ < .01, *F* = 0.731, *p* = .402, Figure 4c) or August 2012 (Figure S7c). There were contrasting effects of temperature on population-level biomasses, partially supporting our fourth hypothesis. Here, *P. palustris* and *P. septentrionis* decreased in biomass with increasing soil temperature (Table S4; Figure S10a,b), while there was an increase in the biomass of the wolf spider *Pirata piraticus* and both *P. diligens* and *P. nigrita* with increasing soil temperature (Table S4; Figure S10c-e). There were no significant changes in biomass of the remaining invertebrate species across the soil temperature gradient (Table S4).

## 3.4 | Plant-mediated effects on invertebrates

There was a significant effect of Pielou's evenness of the plant community on both Pielou's evenness (GAM: $r^2$ = .45, *F* = 5.968, *p* = .021) and Shannon diversity (GAM: $r^2$ = .41, *F* = 5.024, *p* = .033) of the invertebrate community when it was included as a covariate with soil temperature in the analysis. There were no other significant effects of plant community metrics on invertebrate community metrics when the former were included as covariates in the analyses of soil temperature effects on the latter (Table S5).

# 4 | DISCUSSION

## 4.1 | Diversity and percentage cover of plants

There were fewer plant species found as soil temperature increased (Figure 2a), which likely reflects the dearth of species with higher thermal optima in the regional species pool (Kristinsson & Sigurdsson, 2010). Similar effects have been found along elevational gradients, where invasion of warm-adapted species could not offset the loss of native plant species at higher temperatures (De Sassi, Lewis, & Tylianakis, 2012). The early developmental stages of plants are likely to be more sensitive to environmental constraints, with warming shown to decrease the species richness of emerging seedlings (Lloret, Penuelas, & Estiarte, 2004). This may be particularly true for *P. vivipara*, which was rarely found at warmer soil temperatures (Figure S5a) and is most susceptible to heat stress following germination (Marcante, Erschbamer, Buchner, & Neuner, 2014), with a lower thermal tolerance than many of its alpine competitors, such as *Epilobium* sp. (Schwienbacher, Navarro-Cano, Neuner, & Erschbamer, 2012). It was perhaps surprising that there was a greater percentage cover of *V. palustris* in warmer soil (Figure S5e)
given that it typically experiences reduced emergence of seedlings at higher temperatures (Klanderud, Meineri, Töpper, Michel, & Vandvik, 2017), although this may be driven by its increased performance at higher temperatures after loss of dominant competitors (Olsen, Töpper, Skarpaas, Vandvik, & Klanderud, 2016).

While competition for space and soil nutrients should have been lower as a result of fewer species, there were no changes in evenness and thus dominance of plants across the temperature gradient. However, large-scale patterns in declining diversity are often masked by variation at regional or local scales (Walker et al., 2006). There was also an increased turnover of plant species as soil temperature increased (Figure 3a), which is predicted for many plant communities under a warming climate (Thuiller, 2004) and is most likely related to temperature niche separation amongst species (Bertrand et al., 2011). Thus, plant species with a low thermal optimum (e.g. *P. vivipara*) are likely to be excluded in the warmest sites, where they may be replaced by more warm-adapted plants (e.g. *Epilobium* sp.). This phenomenon has been demonstrated in highland areas of France in response to just a 1°C increase in temperature over recent decades (Bertrand et al., 2011). Note that increases in the percentage cover of lichens and *V. serpyllifolia* with increasing soil temperature were driven by a single influential data point in each case (Cook's distance of 9.2 and 9.8 respectively). As such, these trends are unlikely to be representative of general warming impacts on these taxa, especially given the widespread negative effect of temperature on terrestrial lichens (Lang et al., 2012; Walker et al., 2006).

## 4.2 | Invertebrate diversity

Invertebrate species evenness, although not richness, decreased in the warmer sites, leading to an overall decline in Shannon diversity (Figure 2). This suggests that warming had a marked effect on community composition, with the wolf spider *P. piraticus* and beetles from the genus *Pterostichus* dominating the warmer habitat patches (Figure S10). These species have a high thermal tolerance (Berthe et al., 2015; Nørgaard, 1951) and may be able to disperse to their optimum temperature due to their mobility (e.g. Juliano, 1983). These changes in evenness are supported by evidence for a reduction in invertebrate diversity following experimental warming due to shifts in species dominance (Berthe et al., 2015; Villalpando, Williams, & Norby, 2008). There was also high invertebrate species turnover along the temperature gradient, as seen in thermally influenced marine benthic communities (Hillebrand et al., 2010). Here, the harvestman *M. morio* and the beetle *P. septentrionis* in particular were absent from sites >17°C, where they were replaced by the *Pterostichus* beetles (see Figure S9). Similar effects have also been observed in the freshwater streams at Hengill, where community composition changed dramatically in the warmer streams, with a few dominant species reducing the evenness of the invertebrate community (Woodward, Dybkjaer, et al., 2010).

Changes in invertebrate community composition may also have been indirectly mediated by the vegetation, rather than just direct effects of temperature on invertebrate physiology and behaviour

(De Sassi et al., 2012). For example, changes in the evenness of the plant community across the soil temperature gradient were shown to independently contribute to the reduction in evenness and diversity of the invertebrate community, as suggested by the stronger temperature effect on invertebrate community composition after taking plant evenness into account (Table S5). Changes in plant community composition have previously been shown to alter the abundance and species richness of herbivorous invertebrates (Zhou et al., 2015), which might help to explain the changes we observed in the abundance of springtails in our study and even some of the predatory beetles whose diet frequently consists of plant material (Dawson, 1965). Similarly, less diverse vegetation provides fewer refugia for herbivorous invertebrates, resulting in greater top-down control by predators (Sanders, Nickel, Grützner, & Platner, 2008). For example, the predatory *Pterostichus* beetles may be aggregating at warmer patches due to the availability of their springtail prey there (Figure S9; Dawson, 1965; Mundy, Allen-Williams, Underwood, & Warrington, 2000). Indeed, declining plant species richness has stronger negative effects on progressively lower trophic levels in above-ground food webs (Scherber et al., 2010), which could help to explain the loss of invertebrate diversity seen here (Figure 2f). Given that the sampling period for the study took place under the near-equilibrium conditions of late summer, and not during the community assembly phase, it would be interesting to examine plant and invertebrate species composition as a function of time following the spring snowmelt and early plant growth, right through to winter snowfall. Transient differences in plant community assembly as a result of soil temperature decrease could result in legacy effects on the associated invertebrate community (Bale et al., 2002; Dollery et al., 2006; Liu, Reich, Li, & Sun, 2011), especially considering how plant identity may play a more important role in insect herbivore colonisation than temperature alone (Andrew & Hughes, 2007).

## 4.3 | Invertebrate body mass

The mean body mass of the invertebrate community decreased with increasing soil temperature (Figure 4a), in line with our third hypothesis and echoing recent meta-analyses on freshwater communities (Daufresne et al., 2009) and oceanic phytoplankton (Moran, Lopez-Urrutia, Calvo-Diaz, & Li, 2010). This suggests that temperature effects on the behaviour and physiology of invertebrates were strong despite the relatively interconnected nature of the warm and cold sites that we studied and indications that temperature effects are much weaker on the body size of terrestrial than aquatic organisms (Forster, Hirst, & Atkinson, 2012). The general lack of temperature effects on mean body mass at the population level (Table S4) suggests that this finding was driven by the observed species turnover in the invertebrate community (Figure 4b,c), with smaller species at higher soil temperatures rather than smaller individuals within each population in the community (except for parasitoid wasps; see Figure S8a). This finding is of particular note given the short distances between sites of different temperature in the Hengill system
and indicates a clear preferential use of warmer habitat by smaller species (e.g. springtails; see Figure S9f).

In contrast to the population-level component of our third hypothesis, we only observed a decrease in the mean body mass of parasitoid wasps with increasing soil temperature. Furthermore, there was an increase in the mean body mass of springtails and true bug nymphs from 10 to 15°C, before their body mass decreased at higher temperatures (Figure S8b,c), with no effects on the remaining invertebrate species. Variable population-level responses to warming have been regularly reported (Gardner et al., 2011; Watt, Mitchell, & Salewski, 2010), including idiosyncratic trends among insects (Mousseau, 1997; Shelomi, 2012; Walters & Hassall, 2006) and a decrease in the body size of some Carabidae at cooler, high latitudes (Ikeda, Tsuchiya, Nagata, Ito, & Sota, 2012). It is possible that population-level variability in temperature-size responses are a consequence of developmental and functional constraints on body mass and individual plasticity (Angilletta, Steury, & Sears, 2004; Forster, Hirst, & Woodward, 2011; Moya-Larano et al., 2012). Voltinism (i.e. the number of generations per year) may also help to explain the contrasting temperature-size responses of terrestrial species, with long-lived univoltine organisms able to achieve large size over a single growing season, while multivoltine organisms often maximise their fitness by maturing earlier at smaller size (Horne, Hirst, & Atkinson, 2015). Contradictions of temperature-size rules have also been observed in the Hengill streams, at both the community and population levels, potentially as a consequence of temperature-mediated increase in resource availability, or enhanced rates of growth and reproduction (e.g. Adams et al., 2013; O'Gorman et al., 2012). Clearly, more research is needed on this topic to determine the idiosyncratic differences in temperature-size responses for various taxonomic groups.

## 4.4 | Invertebrate abundance and biomass

There was a positive effect of temperature on the total abundance of invertebrates (Figure 4b), in contrast to our fourth hypothesis. Similar increases in the abundance of soil invertebrates with warming are also predicted for Arctic areas in response to changes in vegetation type (Nielsen & Wall, 2013). Shifts in plant community composition, such as those observed here (Figure 3a), have also been shown to mediate a greater abundance of herbivores at higher temperatures (De Sassi et al., 2012). While there were more invertebrates at the warmer sites, the concurrent reduction in mean body size ensured that overall community biomass was unaffected by temperature (Figure 4c). Resource usage by the consumer community thus seemingly kept pace with resource provisioning, with no overall change in percentage cover of the vegetation community as soil temperature increased. Faster rates of vegetation regeneration at higher temperatures (Nishar et al., 2017) may play a key role in meeting the higher metabolic demands of herbivores and subsequently their predators. Arctic and subarctic organisms are also predicted to exhibit relatively small absolute shifts in metabolic rate compared to those at lower latitudes due to their cooler baseline

temperatures and the exponential relationship between metabolism and temperature (Dillon, Wang, & Huey, 2010).

There were contrasting population-level responses to temperature (Figures S9 and S10), which is most likely due to variation in the thermal tolerances of individual species. For example, *M. morio* and *P. septentrionis* are most commonly associated with cold alpine environments (Bakken, 1985; Hein, Pape, Finch, & Löffler, 2014), reflecting their occurrence at only the coldest sites here (Figure S9). Indeed, *M. morio* prefers cool, moist habitats because it is susceptible to desiccation at higher temperatures (Hein et al., 2014). In contrast, *Pterostichus* beetles have been shown to thrive after experimental warming (Berthe et al., 2015), mirroring the increased abundance and biomass of *P. diligens* and *P. nigrita* here (Figures S9 and S10). Contrasting effects of temperature have also been shown on invertebrate populations in German forests (Müller-Kroehling, Jantsch, Fischer, & Fischer, 2014).

The abundance and biomass of species in riparian ecosystems may also be influenced by aquatic subsidies. For example, adult insect emergence can alter terrestrial riparian communities, with the diets of dominant beetle and spider species in some studies shown to consist largely of aquatic insects (Paetzold, Schubert, & Tockner, 2005). Greater insect emergence has been demonstrated from warmer waters (Greig et al., 2012), including at our study system (Hannesdottir et al., 2013), which may provide supplementary resources to predatory invertebrates on the surrounding banks, leading to their aggregation at warmer soil patches. Due to the natal site fidelity observed for many insect species (e.g. Krosch, Baker, Mather, & Cranston, 2011), terrestrial invertebrates may return to these warmer and resource-rich patches across generations, resulting in consistently higher abundances of terrestrial invertebrates at these sites. This trend may be amplified by the effect of temperature on voltinism, with warming shown to increase the number of generations per year in European butterflies and moths (Altermatt, 2010). Thus, temperature effects on stream subsidies to the terrestrial invertebrate community should be explored and quantified to better understand their contribution to the observed patterns.

## 5 | CONCLUSIONS

Warming is widely predicted to negatively affect the diversity, body size and biomass of natural communities (Daufresne et al., 2009; Millennium Ecosystem Assessment 2005; Pimm, 2009). We found support for reductions in the diversity of plant and ground-dwelling invertebrate communities at higher temperatures. We observed a smaller mean body size of invertebrates at the community level in warmer environments, although generally not at the population level, offering mixed support for temperature-size rules. The surprising increase in total abundance of the invertebrate community at warmer sites was offset by their smaller size, resulting in no net change in invertebrate biomass across the temperature gradient. This suggests that consumers had sufficient resources, despite their higher energy demands at warmer temperatures. Inconsistent population-level
effects of temperature on the percentage cover of vegetation and the abundance and biomass of invertebrate species in the system highlight the need for community-level research to understand the ecosystem-level consequences of these context-dependent effects. An analysis of changes in trophic network structure in response to temperature would be a logical next step, to determine how the observed changes in plant and invertebrate community structure might alter the flow of energy through the food web. Understanding the underlying mechanisms of population- and community-level change will help to identify the emergent properties of ecosystems and mitigate the consequences of climate change.

## ACKNOWLEDGEMENTS

We thank Gísli Már Gíslason and Jón S. Ólafsson for providing research support and facilities, and Francois Bochereau for help with the soil chemistry analyses. We acknowledge the funding support of NERC (NE/L011840/1, NE/M020843/1), the British Ecological Society (4009-4884, SR16/1152, 7283/5350), the Royal Society (RG140601), the recruitment fund of the University of Iceland and the Finnish Cultural Foundation.

## AUTHORS' CONTRIBUTIONS

O.B.M. and E.J.O. conceived the ideas and designed methodology; O.B.M., E.J.O. and B.M. collected the data; S.I.R., E.J.O., and B.M. analysed the data; S.I.R. and E.J.O. led the writing of the manuscript. All authors contributed critically to the drafts and gave final approval for publication.

## DATA ACCESSIBILITY

All data is archived with the NERC Environmental Information Data Centre https://doi.org/10.5285/0f074839-1630-4ccd-aa63-84d0da16b28a (O'Gorman, Robinson, McLaughlin, & Marteinsdóttir, 2017).

## ORCID

Sinikka I. Robinson ID http://orcid.org/0000-0002-8831-0497

Bryndís Marteinsdóttir ID http://orcid.org/0000-0003-3779-7327

Eoin J. O'Gorman ID http://orcid.org/0000-0003-4507-5690

## REFERENCES

Adams, G. L., Pichler, D. E., Cox, E. J., O'Gorman, E. J., Seeney, A., Woodward, G., & Reuman, D. C. (2013). Diatoms can be an important exception to temperature-size rules at species and community levels of organization. *Global Change Biology*, *19*, 3540–3552.

Altermatt, F. (2010). Climatic warming increases voltinism in European butterflies and moths. *Proceedings of the Royal Society B: Biological Sciences*, *1685*, 1281–1287. https://doi.org/10.1098/rspb.2009.1910

Andrew, N. R., & Hughes, L. (2005). Diversity and assemblage structure of phytophagous Hemiptera along a latitudinal gradient predicting the

potential impacts of climate change. *Global Ecology and Biogeography*, *14*, 249–262. https://doi.org/10.1111/j.1466-822x.2005.00149.x

Andrew, N. R., & Hughes, L. (2007). Potential host colonization by insect herbivores in a warmer climate: A transplant experiment. *Global Change Biology*, *13*, 1539–1549. https://doi.org/10.1111/j.1365-2486.2007.01393.x

Angilletta, M. J., Steury, T. D., & Sears, M. W. (2004). Temperature, growth rate, and body size in ectotherms: Fitting pieces of a life-history puzzle. *Integrative & Comparative Biology*, *44*, 498–509. https://doi.org/10.1093/icb/44.6.498

Bakken, H. (1985). Cold hardiness in the alpine beetles *Patrobus septentrionis* and *Calathus melanocephalus*. *Journal of Insect Physiology*, *31*, 447–453. https://doi.org/10.1016/0022-1910(85)90091-5

Bale, J. S., Masters, G. H., Hodkinson, I. D., Awmack, C., Martijn Bezemer, T., Brown, V. K., ... Whittaker, J. B. (2002). Herbivory in global climate change research: Direct effects of rising temperature on insect herbivores. *Global Change Biology*, *8*, 1–16. https://doi.org/10.1046/j.1365-2486.2002.00451.x

Benjamini, Y., & Hochberg, Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. *Journal of the Royal Statistical Society Series B*, *57*, 289–300.

Berthe, S. C. F., Derocles, S. A. P., Lunt, D. H., Kimball, B. A., & Evans, D. M. (2015). Simulated climate-warming increases Coleoptera activity-densities and reduces community diversity in a cereal crop. *Agriculture, Ecosystems & Environment*, *2010*, 11–14. https://doi.org/10.1016/j.agee.2015.05.001

Bertrand, R., Lenoir, J., Piedallu, C., Riofrío-Dillon, G., De Ruffray, P., Vidal, C., ... Gegout, J. C. (2011). Changes in plant community composition lag behind climate warming in lowland forests. *Nature*, *479*, 517–520. https://doi.org/10.1038/nature10548

Brose, U., Dunne, A., Montoya, J. M., Petchey, O. L., Schneider, F. D., & Jacob, U. (2012). Climate change in size-structured ecosystems. *Philosophical Transactions of the Royal Society B: Biological Sciences*, *367*, 2903–2912. https://doi.org/10.1098/rstb.2012.0232

Brown, J. H., Gillooly, J. F., Allen, A. P., Savage, V. M., & West, G. B. (2004). Toward a metabolic theory of ecology. *Ecology*, *85*, 1771–1789. https://doi.org/10.1890/03-9000

Chen, I., Hill, J. K., Ohlemuller, R., Roy, D. B., & Thomas, C. D. (2011). Rapid range shifts of species associated with high levels of climate warming. *Science*, *333*, 1024–1026. https://doi.org/10.1126/science.1206432

Daufresne, M., Lengfellner, K., & Sommer, U. (2009). Global warming benefits the small in aquatic ecosystems. *Proceedings of the National Academy of Sciences of the United States of America*, *106*, 12788–12793. https://doi.org/10.1073/pnas.0902080106

Dawson, N. (1965). A comparative study of the ecology of eight species of fenland Carabidae (Coleoptera). *Journal of Animal Ecology*, *34*, 299–314. https://doi.org/10.2307/2652

De Sassi, C., Lewis, O. T., & Tylianakis, J. M. (2012). Plant-mediated and nonadditive effects of two global change drivers on an insect herbivore community. *Ecology*, *93*, 1892–1901. https://doi.org/10.1890/11-1839.1

DeLong, J. P. (2012). Experimental demonstration of a 'rate-size' trade-off governing body size optimization. *Evolutionary Ecology Research*, *14*, 343–352.

Dillon, M. E., Wang, G., & Huey, R. B. (2010). Global metabolic impacts of recent climate warming. *Nature*, *467*, 704–706. https://doi.org/10.1038/nature09407

Dollery, R., Hodkinson, I. D., & Jónsdóttir, I. S. (2006). Impact of warming and timing of snow melt on soil microarthropod assemblages association with Dryas-dominated plant communities on Svalbard. *Ecography*, *29*, 111–119. https://doi.org/10.1111/j.2006-0906-7590.04366.x

Dong, Z., Hou, R., Ouyang, Z., & Zhang, R. (2013). Tritrophic interaction influenced by warming and tillage: A field study on winter wheat,
aphids and parasitoids. *Agriculture, Ecosystems & Environment*, **181**, 144-148. https://doi.org/10.1016/j.agee.2013.09.009

Forster, J., Hirst, A. G., & Atkinson, D. (2012). Warming-induced reductions in body size are greater in aquatic than terrestrial species. *Proceedings of the National Academy of Sciences of the United States of America*, **109**, 19310-19314. https://doi.org/10.1073/pnas.1210460109

Forster, J., Hirst, A. G., & Woodward, G. (2011). Growth and development rates have different thermal responses. *The American Naturalist*, **178**, 668-678. https://doi.org/10.1086/662174

Friberg, N., Dybkjaer, J. B., Olafsson, J. S., Gislason, G. M., Larsen, S. E., & Lauridsen, T. L. (2009). Relationship between structure and function in streams contrasting in temperature. *Freshwater Biology*, **54**, 2051-2068. https://doi.org/10.1111/j.1365-2427.2009.02234.x

Gardner, J. L., Peters, A., Kearney, M. R., Joseph, L., & Heinsohn, R. (2011). Declining body size: A third universal response to warming? *Trends in Ecology & Evolution*, **26**, 285-291. https://doi.org/10.1016/j.tree.2011.03.005

Greig, H. S., Kratina, P., Thompson, P. L., Palen, W. J., Richardson, J. S., & Shurin, J. B. (2012). Warming, eutrophication, and predator loss amplify subsidies between aquatic and terrestrial ecosystems. *Global Change Biology*, **18**, 504-514. https://doi.org/10.1111/j.1365-2486.2011.02540.x

Hannesdottir, E. R., Gislason, G. M., Olafsson, J. S., Olafsson, O. P., & O'Gorman, E. J. (2013). Increased stream productivity with warming supports higher trophic levels. *Advances in Ecological Research*, **48**, 285-333. https://doi.org/10.1016/B978-0-12-417199-2.00005-7

Hein, N., Pape, R., Finch, O. D., & Löffler, J. (2014). Alpine activity patterns of Mitopus morio (Fabricius, 1779) are induced by variations in temperature and humidity at different scales in central Norway. *Journal of Mountain Science*, **11**, 644-655.

Hering, D., Schmidt-Kloiber, A., Murphy, J., Lücke, S., Zamora-Muñoz, C., López-Rodríguez, M. J., ... Graf, W. (2009). Potential impact of climate change on aquatic insects: A sensitivity analysis for European caddisflies (Trichoptera) based on distribution patterns and ecological preferences. *Aquatic Sciences*, **71**, 3-14. https://doi.org/10.1007/s00027-009-9159-5

Hillebrand, H., Soininen, J., & Snoeijs, P. (2010). Warming leads to higher species turnover in a coastal ecosystem. *Global Change Biology*, **16**, 1181-1193. https://doi.org/10.1111/j.1365-2486.2009.02045.x

Horne, C. R., Hirst, A., & Atkinson, D. (2015). Temperature-size responses match latitudinal-size clines in arthropods, revealing critical differences between aquatic and terrestrial species. *Ecology Letters*, **18**, 327-335. https://doi.org/10.1111/ele.12413

Ikeda, H., Tsuchiya, Y., Nagata, N., Ito, M. T., & Sota, T. (2012). Altitudinal life-cycle and body-size variation in ground beetles of the genus *Carabus* (subgenus *Ohomopterus*) in relation to the temperature conditions and prey earworms. *Pedobiologia*, **55**, 67-73. https://doi.org/10.1016/j.pedobi.2011.10.008

Ings, T. C., Montoya, J. M., Bascompte, J., Blüthgen, N., Brown, L., Dormann, C. F., ... Woodward, G. (2009). Review: Ecological networks - Beyond food webs. *Journal of Animal Ecology*, **78**, 253-269. https://doi.org/10.1111/j.1365-2656.2008.01460.x

IPCC. (2014). *Climate Change 2014: Synthesis Report*. In R. K. Pachauri & L. A. Meyer (Eds.), *Contribution of working groups I, II and III to the fifth assessment report of the intergovernmental panel on climate change*. Geneva, Switzerland: IPCC, 151 pp.

Jonsson, T., Cohen, J. E., & Carpenter, S. R. (2005). Food webs, body size, and species abundance in ecological community description. *Advances in Ecological Research*, **36**, 1-84.

Juliano, S. A. (1983). Body size, dispersal ability, and range size in North American species of Brachinus (Coleoptera: Carabidae). *The Coleopterists Bulletin*, **37**, 232-238.

Keele, L. J. (2008). *Semiparametric regression for the social sciences* (p. 230). Chichester, UK: John Wiley & Sons.

Kingsolver, J. G., & Huey, R. B. (2008). Size, temperature, and fitness: Three rules. *Evolutionary Ecology Research*, **10**, 251-268.

Klanderud, K., Meineri, E., Töpper, J., Michel, P., & Vandvik, V. (2017). Biotic interaction effects on seedling recruitment along bioclimatic gradients: Testing the stress-gradient hypothesis. *Journal of Vegetation Science*, **28**, 347-356. https://doi.org/10.1111/jvs.12495

Klanderud, K., & Totland, O. (2005). The relative importance of neighbours and abiotic environmental conditions for population dynamic parameters of two alpine species. *Journal of Ecology*, **93**, 493-501. https://doi.org/10.1111/j.1365-2745.2005.01000.x

Klein, J. A., Harte, J., & Zhao, X. (2004). Experimental warming causes large and rapid species loss, dampened by stimulated grazing on the Tibetan Plateau. *Ecology Letters*, **7**, 1170-1179. https://doi.org/10.1111/j.1461-0248.2004.00677.x

Kline, R. B. (2015). *Principles and practice of structural equation modeling*. New York, NY: Guilford Publications, 534 pp.

Klok, C. J., & Harrison, J. F. (2013). The temperature size rule in arthropods: Independence of macro-environmental variables but size dependent. *Integrative & Comparative Biology*, **53**, 557-570. https://doi.org/10.1093/icb/ict075

Krajick, K. (2004). All downhill from here? *Science*, **303**, 1601-1602.

Kristinsson, H., & Sigurdsson, S. V. (2010). A guide to the flowering plants and ferns of Iceland (3rd ed., p. 311). Reykjavik, Iceland: Mál og Menning.

Krosch, M. N., Baker, A. M., Mather, P. B., & Cranston, P. S. (2011). Spatial population genetic structure reveals strong natal site fidelity in *Echinocladius martini* (Diptera: Chironomidae) in northeast Queensland, Australia. *Freshwater Biology*, **56**, 1328-1341. https://doi.org/10.1111/j.1365-2427.2010.02571.x

Lang, S. I., Cornelissen, J. H., Shaver, G. R., Ahrens, M., Callaghan, T. V., Molau, U., ... Aerts, R. (2012). Arctic warming on two continents has consistent negative effects on lichen diversity and mixed effects on bryophyte diversity. *Global Change Biology*, **18**, 1096-1107. https://doi.org/10.1111/j.1365-2486.2011.02570.x

Lejeusne, C., Chevaldonne, P., Pergent-Martini, C., Bourdouresque, C. F., & Perez, T. (2010). Climate change effects on a miniature ocean: The highly diverse, highly impacted Mediterranean Sea. *Trends in Ecology and Evolution*, **25**, 250-260. https://doi.org/10.1016/j.tree.2009.10.009

Lemoine, N. P., & Burkepile, D. E. (2012). Temperature-induced mismatches between consumption and metabolism reduce consumer fitness. *Ecology*, **93**, 2483-2489. https://doi.org/10.1890/12-0375.1

Liu, Y., Reich, P. B., Li, G., & Sun, S. (2011). Shifting phenology and abundance under experimental warming alters trophic relationships and plant reproductive capacity. *Ecology*, **92**, 1201-1207. https://doi.org/10.1890/10-2060.1

Lloret, F., Penuelas, J., & Estiarte, M. (2004). Experimental evidence of reduced diversity of seedlings due to climate modification in a Mediterranean-type community. *Global Change Biology*, **10**, 248-258. https://doi.org/10.1111/j.1365-2486.2004.00725.x

Lurgi, M., López, B. C., & Montoya, J. M. (2012). Climate change impacts on body size and food web structure on mountain ecosystems. *Philosophical Transactions of the Royal Society B: Biological Sciences*, **367**, 3050-3057. https://doi.org/10.1098/rstb.2012.0239

MacCallum, R. C., Browne, M. W., & Sugawara, H. M. (1996). Power analysis and determination of sample size for covariance structure modeling. *Psychological Methods*, **1**, 130. https://doi.org/10.1037/1082-989X.1.2.130

Marcante, S., Erschbamer, B., Buchner, O., & Neuner, G. (2014). Heat tolerance of early developmental stages of glacier foreland species in the growth chamber and in the field. *Plant Ecology*, **215**, 747-758. https://doi.org/10.1007/s11258-014-0361-8

Menzel, A., Sparks, T. H., & Estrella, N. (2006). European phenological response to climate change matches the warming pattern.
Global Change Biology, 12, 1969-1976. https://doi.org/10.1111/j.1365-2486.2006.01193.x

Millennium Ecosystem Assessment. (2005). *Millennium ecosystem assessment ecosystems and human well-being: Biodiversity synthesis*. Washington, DC: World Resources Institute.

Montoya, J. M., & Raffaelli, D. (2010). Climate change, biotic interactions and ecosystem services. *Philosophical Transactions of the Royal Society B: Biological Sciences*, 365, 2013-2018. https://doi.org/10.1098/rstb.2010.0114

Moran, X. A. G., Lopez-Urrutia, A., Calvo-Diaz, A., & Li, W. K. W. (2010). Increasing importance of small phytoplankton in a warmer ocean. *Global Change Biology*, 16, 1137-1144. https://doi.org/10.1111/j.1365-2486.2009.01960.x

Mousseau, T.A. (1997). Ectotherms follow the converse to Bergmann's rule. *Evolution*, 51, 630-632. https://doi.org/10.1111/j.1558-5646.1997.tb02453.x

Moya-Larano, J., Verdeny-Vilalta, O., Rowntree, J., Melguizo-Ruiz, N., Montserrat, M., & Laiolo, P. (2012). Climate change and evolutionary dynamics in food webs. *Advances in Ecological Research*, 47, 1-36.

Müller-Kroehling, S., Jantsch, M. C., Fischer, H. S., & Fischer, A. (2014). Modelling the effects of global warming on the ground beetle (Coleoptera: Carabidae) fauna on beech forest in Bavaria, Germany. *European Journal of Entomology*, 111, 35-49. https://doi.org/10.14411/eje.2014.005

Mundy, C., Allen-Williams, L., Underwood, N., & Warrington, S. (2000). Prey selection and foraging behaviour by *Pterostichus cupreus* L. (Col., Carabidae) under laboratory conditions. *Journal of Applied Entomology*, 124, 349-358. https://doi.org/10.1046/j.1439-0418.2000.00491.x

Muren, C., Hoffmann, F., & Kwak, M. M. (2003). Insect diversity on yellow Asteraceae in road verges in the Netherlands. *Proceedings of the Section Experimental and Applied Entomology - Netherlands Entomological Society*, 14, 115-118.

Nielsen, U. N., & Wall, D. H. (2013). The future of soil invertebrate communities in polar regions: Different climate change responses in the Arctic and Antarctic? *Ecology Letters*, 16, 409-419. https://doi.org/10.1111/ele.12058

Nishar, A., Bader, M. K. F., O'Gorman, E. J., Deng, J., Breen, B., & Leuzinger, S. (2017). Temperature effects on biomass and regeneration of vegetation in a geothermal area. *Frontiers in Plant Science*, 8, 249.

Nørgaard, E. (1951). On the ecology of two lycosid spiders (*Pirata piraticus* and *Lycosa pullata*) from a Danish sphagnum bog. *Oikos*, 3, 1-21. https://doi.org/10.2307/3565128

O'Connor, M. I., Gilbert, B., & Brown, C. J. (2011). Theoretical predictions for how temperature affects the dynamics of interacting herbivores and plants. *The American Naturalist*, 178, 626-638. https://doi.org/10.1086/662171

O'Gorman, E. J., Benstead, J. P., Cross, W. F., Friberg, N., Hood, J. M., Johnson, P. W., ... Woodward, G. (2014). Climate change and geo-thermal ecosystems: Natural laboratories, sentinel systems, and future refugia. *Global Change Biology*, 20, 3291-3299. https://doi.org/10.1111/gcb.12602

O'Gorman, E. J., Olafsson, O. P., Demars, B. O., Friberg, N., Gudbergsson, G., Hannesdottir, E. R., ... Gislason, G. M. (2016). Temperature effects on fish production across a natural thermal gradient. *Global Change Biology*, 9, 3206-3220. https://doi.org/10.1111/gcb.13233

O'Gorman, E. J., Pichler, D. E., Adams, G., et al. (2012). Impacts of warming on the structure and functioning of aquatic communities: Individual- to ecosystem-level responses. *Advances in Ecological Research*, 47, 81-176. https://doi.org/10.1016/B978-0-12-398315-2.00002-8

O'Gorman, E. J., Robinson, S. I., McLaughlin, Ó. B., & Marteinsdóttir, B. (2017). Invertebrate biomass, vegetation cover, and environmental data from Hengill, Iceland.

NERC Environmental Information Data Centre. https://doi.org/10.5285/0f074839-1630-4ccd-aa63-84d0da16b28a

Olsen, S. L., Töpper, J. P., Skarpaas, O., Vandvik, V., & Klanderud, K. (2016). From facilitation to competition: Temperature-driven shift in dominant plant interactions affects population dynamics in semi-natural grasslands. *Global Change Biology*, 22, 1915-1926. https://doi.org/10.1111/gcb.13241

Paetzold, A., Schubert, C. J., & Tockner, C. (2005). Aquatic terrestrial linkages along a braided-river: Riparian arthropods feeding on aquatic insects. *Ecosystems*, 8, 748-759. https://doi.org/10.1007/s10021-005-0004-y

Parmesan, C., & Yohe, G. (2003). A globally coherent fingerprint of climate change impacts across natural systems. *Nature*, 421, 37-42. https://doi.org/10.1038/nature01286

Peet, R. K., Wenthworth, T. R., & White, P. S. (1998). A flexible, multipurpose method for recording vegetation composition and structure. *Castanea*, 63, 262-274.

Petchey, O. L., McPhearson, P. T., Casey, T. M., & Morin, P. J. (1999). Environmental warming alters food-web structure and ecosystem function. *Nature*, 402, 69-72. https://doi.org/10.1038/47023

Pimm, S. L. (2009). Climate disruption and biodiversity. *Current Biology*, 19, 595-601. https://doi.org/10.1016/j.cub.2009.05.055

Porter, W. P., & Kearney, M. (2009). Size, shape and the thermal niche of endotherms. *Proceedings of the National Academy of Sciences of the United States of America*, 106, 666-672.

Reuman, D. C., Holt, R. D., & Yvon-Durocher, G. (2014). A metabolic perspective on competition and body size reductions with warming. *Journal of Animal Ecology*, 83, 59-69. https://doi.org/10.1111/1365-2656.12064

Roy, B. A., Güsewell, S., & Harte, J. (2004). Response of plant pathogens and herbivores to a warming experiment. *Ecology*, 85, 2570-2581. https://doi.org/10.1890/03-0182

Rüger, T., & Sommer, U. (2012). Warming does not always benefit the small - Results from a plankton experiment. *Aquatic Botany*, 97, 64-68. https://doi.org/10.1016/j.aquabot.2011.12.001

Rustad, L. E., Campbell, J. L., Marion, G. M., Norby, R. J., Mitchell, M. J., Hartley, A. E., ... GCTE-NEWS (2001). A meta-analysis of the response of soil respiration, net nitrogen mineralization, and aboveground plant growth to experimental warming. *Oecologia*, 126, 543-562. https://doi.org/10.1007/s004420000544

Sanders, D., Nickel, H., Grützner, T., & Platner, C. (2008). Habitat structure mediates top-down effects of spiders and ants on herbivores. *Basic and Applied Ecology*, 9, 152-160. https://doi.org/10.1016/j.baae.2007.01.003

Scherber, C., Eisenhauer, N., Weisser, W. W., Schmid, B., Voigt, W., Fischer, M., ... Allan, E. (2010). Bottom-up effects of plant diversity on multitrophic interactions in a biodiversity experiment. *Nature*, 468, 553-556. https://doi.org/10.1038/nature09492

Schwienbacher, E., Navarro-Cano, J. A., Neuner, G., & Erschbamer, B. (2012). Correspondence of seed traits with niche position in glacial foreland succession. *Plant Ecology*, 213, 371-382. https://doi.org/10.1007/s11258-011-9981-4

Sharp, C. E., Brady, A. L., Sharp, G. H., Grasby, S. E., Stott, M. B., & Dunfield, P. F. (2014). Humboldt's spa: Microbial diversity is controlled by temperature in geothermal environments. *ISME Journal*, 8, 1166-1174. https://doi.org/10.1038/ismej.2013.237

Shelomi, M. (2012). Where are we now? Bergmann's rule sensu lato in insects. *The American Naturalist*, 180, 511-519. https://doi.org/10.1086/667595

Somero, G. N. (2010). The physiology of climate change: How potentials for acclimatization and genetic adaptation will determine 'winners' and 'losers'. *The Journal of Experimental Biology*, 213, 912-920. https://doi.org/10.1242/jeb.037473

Thomas, C. D., Cameron, A., Green, R. E., Bakkenes, M., Beaumont, L. J., Collingham, Y. C., ... Williams, S. E. (2004). Extinction risk from
climate change. *Nature*, 427, 145-148. https://doi.org/10.1038/
nature02121

Thuiller, W. (2004). Patterns and uncertainties of species' range shifts under climate change. *Global Change Biology*, *10*, 2020-2027. https://doi.org/10.1111/j.1365-2486.2004.00859.x

Vedel, V., Chipman, A. D., Akam, M., & Arthur, W. (2008). Temperature-dependent plasticity of segment number in an arthropod species: The centipede *Strigamia maritime*. *Evolution & Development*, *10*, 487-492. https://doi.org/10.1111/j.1525-142X.2008.00259.x

Villalpando, S. N., Williams, R. S., & Norby, R. J. (2008). Elevated air temperature alters an old-field insect community in a multi factor climate change experiment. *Global Change Biology*, *15*, 930-942.

Vucic-Pestic, O., Ehnes, R. B., Rall, B. C., & Brose, U. (2011). Warming up the system: Higher predator feeding rates but lower energetic efficiencies. *Global Change Biology*, *17*, 1301-1310. https://doi.org/10.1111/j.1365-2486.2010.02329.x

Walker, M. D., Wahren, C. H., Hollister, R. D., Henry, G. H., Ahlquist, L. E., Alatalo, J. M., ... Wookey, P. A. (2006). Plant community responses to experimental warming across the tundra biome. *Proceedings of the National Academy of Sciences of the United States of America*, *103*, 1342-1346. https://doi.org/10.1073/pnas.0503198103

Walters, J., & Hassall, M. (2006). The temperature-size rule in ectotherms: May general explanations exist after all? *The American Naturalist*, *167*, 510-523.

Walther, G. R. (2010). Community and ecosystem responses to recent climate change. *Philosophical Transactions of the Royal Society B: Biological Sciences*, **365**, 2019-2024. https://doi.org/10.1098/rstb.2010.0021

Walther, G., Post, E., Convey, P., Menzel, A., Parmesan, C., Beebee, T. J. C., ... Bairlein, F. (2002). Ecological responses to recent climate change. *Nature*, **416**, 389-395. https://doi.org/10.1038/416389a

Watt, C., Mitchell, S., & Salewski, V. (2010). Bergmann's rule; a concept cluster? *Oikos*, *119*, 89-100. https://doi.org/10.1111/j.1600-0706.2009.17959.x

Woodcock, B. A. (2007). Chapter 3: Pitfall trapping in ecological studies. In S. R. Leather (Ed.), *Insect sampling in forest ecosystems* (pp. 37-53). Oxford, UK: Blackwell Science Ltd.

Woodward, G., Dybkjaer, J. B., Olafsson, J. S., Gislason, G. M., Hannesdóttir, E. R., & Friberg, N. (2010). Sentinel systems on a razor's edge: Effects of warming on Arctic geothermal stream ecosystems. *Global Change Biology*, **16**, 1979-1991. https://doi.org/10.1111/j.1365-2486.2009.02052.x

Woodward, G., Perkins, D. M., & Brown, L. E. (2010). Climate change and freshwater ecosystems: Impacts across multiple levels of organization. *Philosophical Transactions of the Royal Society B: Biological Sciences*, **365**, 2093-2106. https://doi.org/10.1098/rstb.2010.0055

Yvon-Durocher, G., Montoya, J. M., Trimmer, M., & Woodward, G. (2011). Warming alters the size spectrum and shifts the distribution of biomass in freshwater ecosystems. *Global Change Biology*, **17**, 1681-1694. https://doi.org/10.1111/j.1365-2486.2010.02321.x

Zhou, H., Zou, X., Wang, D., Wan, S., Wang, L., & Gou, J. (2015). Responses of community-level plant-insect interactions to climate warming in a meadow steppe. *Scientific Reports*, **5**, 18654.

## SUPPORTING INFORMATION

Additional Supporting Information may be found online in the
supporting information tab for this article.

**How to cite this article:** Robinson SI, McLaughlin ÓB, Marteinsdóttir B, O'Gorman EJ. Soil temperature effects on the structure and diversity of plant and invertebrate communities in a natural warming experiment. *J Anim Ecol.* 2018;87:634-646. https://doi.org/10.1111/1365-2656.12798