Larger but younger fish when growth outpaces mortality in heated ecosystem
Abstract
Ectotherms are predicted to ‘shrink’ with global warming, in line with general growth models and the temperaturesize rule (TSR), both predicting smaller adult sizes with warming. However, they also predict faster juvenile growth rates and thus larger sizeatage of young organisms. Hence, the result of warming on the sizestructure of a population depends on the interplay between how mortality rate, juvenile and adult growth rates are affected by warming. Here, we use twodecade long time series of biological samples from a unique enclosed bay heated by cooling water from a nearby nuclear power plant to become 5–10 °C warmer than its reference area. We used growthincrement biochronologies (12,658 reconstructed lengthatage estimates from 2426 individuals) to quantify how >20 years of warming has affected body growth, sizeatage, and catch to quantify mortality rates and population size and age structure of Eurasian perch (Perca fluviatilis). In the heated area, growth rates were faster for all sizes, and hence sizeatage was larger for all ages, compared to the reference area. While mortality rates were also higher (lowering mean age by 0.4 years), the faster growth rates lead to a 2 cm larger mean size in the heated area. Differences in the sizespectrum exponent (describing how the abundance declines with size) were less clear statistically. Our analyses reveal that mortality, in addition to plastic growth and sizeresponses, is a key factor determining the size structure of populations exposed to warming. Understanding the mechanisms by which warming affects the size and the age structure of populations is critical for predicting the impacts of climate change on ecological functions, interactions, and dynamics.
Editor's evaluation
This work provides convincing evidence to refute a general tenet in biology, that warming induces smaller maximum body sizes in adult ectoterm individuals. Using a seminatural experiment in an exceptional manmade ecosystem, the authors demostrate that fish in waters warmed by a nearby nuclear plant grew faster but died younger, causing little effect on the size distribution of fish in the area. This work will be of interest to ecologists and physiologists interested in the impacts of global warming on natural communities.
https://doi.org/10.7554/eLife.82996.sa0Introduction
Ectotherm species, constituting 99% of species globally (Atkinson and Sibly, 1997; Wilson, 1992), are commonly predicted to shrink in a warming world (Gardner et al., 2011; Sheridan and Bickford, 2011). However, as the size distribution of many species spans several orders of magnitude, and temperature effects on size may depend on size or age, it is important to be specific about which sizes or life stages are predicted to shrink (usually mean or adult is meant). For instance, warming can shift size distributions without altering mean size if increases in juvenile sizeatage outweigh the decline in sizeatage in adults, which is consistent with the temperature size rule, TSR (Atkinson, 1994). Resolving how warming induces changes in population’s size distribution may thus be more instructive (Fritschie and Olden, 2016), especially for inferring warming effects on species’ ecological role, biomass production, or energy fluxes (Gårdmark and Huss, 2020; Yvondurocher et al., 2011). This is because key processes such as metabolism, feeding, growth, and mortality scale with body size (Andersen and Link, 2020; Blanchard et al., 2017; Brown et al., 2004; Pauly, 1980; Thorson et al., 2017; Ursin, 1967). Hence, as the value of these traits at mean body size is not the same as the mean population trait value (Bernhardt et al., 2018), the size distribution within a population matters for its dynamics and for how it changes under warming.
The population size distribution can be represented as a sizespectrum, which generally is the frequency distribution of individual body sizes (Edwards et al., 2017). It is often described in terms of the sizespectrum slope (slope of individuals or biomass of a size class over the mean size of that class on a loglog scale [Edwards et al., 2017; Sheldon et al., 1973; White et al., 2007]) or simply the exponent of the power law individual size distribution (Edwards et al., 2017). The sizespectrum thus results from temperaturedependent ecological processes such as body growth, mortality, and recruitment (Blanchard et al., 2017; Heneghan et al., 2019). Despite its rich theoretical foundation (Andersen, 2019) and usefulness as an ecological indicator (Blanchard et al., 2005), few studies have evaluated warming effects on the species sizespectrum in largerbodied species (but see Blanchard et al., 2005), and none in large scale experimental setups. There are numerous paths by which a species’ sizespectrum could change with warming (Heneghan et al., 2019). For instance, in line with TSR predictions, warming may lead to a smaller sizespectrum exponents (steeper slope) if the maximum size declines. However, changes in sizeatage and the relative abundances of juveniles and adults may alter this decline in the sizespectrum slope. Warming can also lead to elevated mortality (Barnett et al., 2020; Berggren et al., 2022; Biro et al., 2007; Pauly, 1980), partly because a faster pace of life with higher metabolic rates is associated with a shorter lifespan (Brown et al., 2004; Munch and Salinas, 2009) or due to direct lethal effects of extreme temperature events. This truncates the age distribution towards younger individuals (Barnett et al., 2017), which may reduce density dependence and potentially increase growth rates, thus countering the effects of mortality on the sizespectrum exponent. However, not all sizes may benefit from warming, as e.g. the optimum temperature for growth declines with size (Lindmark et al., 2022). Hence, the effect of warming on the sizespectrum depends on several interlinked processes affecting abundanceatsize and sizeatage.
Sizeatage is generally predicted to increase with warming for small individuals, but decrease for large individuals according to the mentioned TSR (Atkinson, 1994; Ohlberger, 2013). Several factors likely contribute to this pattern, such as increased allocation to reproduction (Wootton et al., 2022) and larger individuals in fish populations having optimum growth rates at lower temperatures (Lindmark et al., 2022). Empirical support in fishes for this pattern seems to be more consistent for increases in sizeatage of juveniles (Huss et al., 2019; Rindorf et al., 2008; Thresher et al., 2007) than declines in adult sizeatage (but see Baudron et al., 2014; Oke et al., 2022; Smoliński et al., 2020), for which a larger diversity in responses is observed among species (Barneche et al., 2019; e.g., Huss et al., 2019). However, most studies have been done on commercially exploited species, since longtime series are more common in such species. This may confound or interact with the effects of temperature because fishing mortality can affect densitydependent growth (van Gemert et al., 2018), but also select for slowgrowing individuals and changes in maturation processes, which also influences growth trajectories (Audzijonyte et al., 2016).
The effect of temperature on mortality rates of wild populations is often studied using amongspecies analyses (Pauly, 1980; Thorson et al., 2017). These relationships based on thermal gradients in space may not necessarily be the same as the effects of warming on mortality in single populations. Hence, the effects of warming on growth and sizeatage, and mortality within natural populations constitute a key knowledge gap for predicting the consequences of climate change on population size spectra.
Here, we used data from a unique, largescale 23yearlong heating experiment of a coastal ecosystem to quantify how warming changed fish body growth, mortality, and the size structure in an unexploited population of Eurasian perch (Perca fluviatilis, ‘perch’). We compare fish from this enclosed bay exposed to temperatures approximately 5–10 °C above normal (‘heated area’) with fish from a reference area in the adjacent archipelago (Figure 1). Using hierarchical Bayesian models, we quantify differences in key individual and populationlevel parameters, such as body growth, asymptotic size, mortality rates, and size spectra, between the heated and reference coastal areas.
Results
Analysis of perch sizeatage using the von Bertalanffy growth equation (VBGE) revealed that fish cohorts (year classes) in the heated area both grew faster initially (larger sizeatage) and reached larger predicted asymptotic sizes than those in the unheated reference area (Figure 2). The model with areaspecific VBGE parameters (${L}_{\infty}$ , $K$, and ${t}_{0}$) had the best outofsample predictive accuracy (the largest expected log pointwise predictive density for a new observation; Supplementary file 1a). Models where both ${L}_{\infty}$ and $K$ were shared did not converge (Supplementary file 1a). Both the estimated values for fish asymptotic length (${L}_{\infty}$) and growth coefficient ($K$) were larger in the heated compared to the reference area (Figure 2—figure supplement 8). We estimated that the asymptotic length of fish in the heated area was 16% larger than in the reference area (calculated as $\frac{{L}_{\infty \mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}{L}_{\infty \mathrm{r}\mathrm{e}\mathrm{f}}}{{L}_{\infty \mathrm{r}\mathrm{e}\mathrm{f}}}$) (${L}_{\infty \mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}=45.7[36.8,56.3]$, ${L}_{\infty \mathrm{r}\mathrm{e}\mathrm{f}}=39.4\left[35.4,43.9\right]$ , where the point estimate is the posterior median and values in brackets correspond to the 95% credible interval). The growth coefficient was 27% larger in the heated area (${K}_{\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}=0.19[0.15,\text{}0.23]$, ${K}_{\mathrm{r}\mathrm{e}\mathrm{f}}=0.15[0.12,\text{}0.17]$). These differences in growth parameters lead to fish being approximately 7%–11% larger in the heated area at any age relative to the reference area (Figure 2—figure supplement 4). Due to the last three cohorts (1995–1997) having large estimates of ${L}_{\infty \mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}$ and low $K$ (potentially due to their negative correlation and because of the young age with data far from the asymptote, Figure 2—figure supplements 3 and 5–6), we fit the same model with these cohorts omitted to evaluate the importance of those for the predicted difference between the areas. Without these, the predicted difference in sizeatage was still clear, but smaller (between 4%–7%, Figure 2—figure supplements 9–10).
In addition, we found that growth rates in the reference area were both slower and declined faster with size compared to the heated area (Figure 3). The best model for growth ($G=\alpha {L}^{\theta}$) had areaspecific $\alpha $ and $\theta $ parameters (Supplementary file 1b). Initial growth ($\alpha $) was estimated to be 18% faster in the heated than in the reference area (${\alpha}_{\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}=512\left[\mathrm{462,565}\right]$, ${\alpha}_{\mathrm{r}\mathrm{e}\mathrm{f}}=433\left[\mathrm{413,454}\right]$), and the growth of fish in the heated area declines more slowly with length than in the reference area (${\theta}_{\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}=1.13[1.16,1.11]$, ${\theta}_{\mathrm{r}\mathrm{e}\mathrm{f}}=1.18[1.19,1.16]$). The distribution of differences of the posterior samples for α and only had 0.3% and 0.2% of the density below 0, respectively (Figure 3C and E), indicating a high probability that lengthbased growth rates are faster in the heated area.
By analyzing the decline in catchperuniteffort over age, we found that the instantaneous mortality rate $Z$ (the rate at which log abundance declines with age) is higher in the heated area (Figure 4). $Z$ was estimated as a fixed effect, as the model where only intercepts varied among years had the best outofsample predictive ability. The overlap with zero is 0.07% for the distribution of differences between posterior samples of ${Z}_{\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}$ and ${Z}_{\mathrm{r}\mathrm{e}\mathrm{f}}$ (Figure 4C). We estimated ${Z}_{\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}$ to be 0.73 [0.66,0.79] and ${Z}_{\mathrm{r}\mathrm{e}\mathrm{f}}$ to be 0.62 [0.58,0.67], which corresponds to annual mortality rates (calculated as $1{e}^{Z}$) of 52% in the heated area and 46% in the reference area.
Lastly, analysis of the size and agestructure in the two areas revealed that, despite the faster growth rates, higher mortality, and larger maximum sizes in the heated area Figure 5A, the sizespectrum exponents were largely similar Figure 5B, C. In fact, the sizespectrum exponent was only slightly larger in the heated area (Figure 5B), and their 95% confidence intervals largely overlap. However, results from the lognormal model fitted to the size and agedistributions revealed that the average size was two centimeters longer and the average age 0.4 years younger in the heated compared to the reference area (Figure 6).
Discussion
Our study provides strong evidence for warminginduced differentiation in growth and mortality in a natural population of an unexploited, temperate fish species exposed to an ecosystemscale experiment with 5–10°C above normal temperatures for more than two decades. Interestingly, these effects largely, but not completely, counteract each other when it comes to population sizestructure—while the fish are younger, they are also larger on average. However, differences in the rate of decline in abundance with size are less pronounced between the areas. It is difficult to generalize these findings since it is a study of only a single species. It is, however, a unique climate change experiment, as experimental studies on fish to date are much shorter and often on scales much smaller than whole ecosystems, and longtime series of biological samples exist mainly for commercially exploited fish species (Baudron et al., 2014; Smoliński et al., 2020; Thresher et al., 2007) (in which fisheries exploitation affects sizestructure both directly and indirectly by selecting for fastgrowing individuals). While factors other than temperature could have contributed to the observed elevated growth and mortality, the temperature contrast is unusually large for natural systems (i.e. 5–10°C, which can be compared to the 1.35°C change in the Baltic Sea between 1982 and 2006 [Belkin, 2009]). Moreover, heating occurred at the scale of a whole ecosystem, which makes the findings highly relevant in the context of global warming.
Interestingly, our findings contrast with both broader predictions about declining mean or adult body sizes based on the GOLT hypothesis (Cheung et al., 2013; Pauly, 2021), and with intraspecific patterns such as the TSR (temperaturesize rule, Atkinson, 1994). The contrasts lie in that both asymptotic size and sizeatage of mature individuals, as well as the proportion of larger individuals, were slightly larger and higher in the heated area—despite the elevated mortality rates. This result was unexpected for two reasons: optimum growth temperatures generally decline with body size within species under food satiation in experimental studies (Lindmark et al., 2022), and fish tend to mature at smaller body sizes and allocate more energy into reproduction as it gets warmer (Niu et al., 2023; Wootton et al., 2022). Both patterns have been used to explain how growth can increase for small and young fish, while large and old fish typically do not benefit from warming. Our study species is no exception to these rules (Huss et al., 2019; Karås and Thoresson, 1992; Niu et al., 2023; Sandstrom et al., 1995). This suggests that growth dynamics under food satiation may not be directly proportional to those under natural feeding conditions (Railsback, 2022). It could also mean that while temperatures is near optimum for growth in the warmest months of the year for a 15 cm individual (and above optimum for larger fish as the optimum declines with size) (Huss et al., 2019; Lindmark et al., 2022), the exposure to such high temperatures is not enough to cause strong reductions in growth and eventually sizeatage. Our results highlight that we need to focus on understanding to what extent the commonly observed increase in sizeatage for juveniles in warm environments can be maintained as they grow older.
Our finding that mortality rates were higher in the heated area was expected—warming leads to faster metabolic rates (faster ‘pace of life’), which in turn is associated with a shorter life span (Brown et al., 2004; McCoy and Gillooly, 2008; Munch and Salinas, 2009). Extreme temperatures, which may be more common in warmed systems under natural variability, can also be lethal if e.g., acute oxygen demands cannot be met (Sandblom et al., 2016). Warming may further increase predation mortality, as predators’ feeding rates increase in order to meet the higher demand for food (Biro et al., 2007; Pauly, 1980; Ursin, 1967). However, most evidence to date of the temperature dependence of mortality rates in natural populations stems from acrossspecies studies (Gislason et al., 2010; Pauly, 1980; Thorson et al., 2017, but see Berggren et al., 2022; Biro et al., 2007). Acrossspecies relationships are not necessarily determined by the same processes as withinspecies relationships; thus, our finding of warminginduced mortality in a heated vs control environment in two nearby conspecific populations is important.
Since a key question for understanding the implications of warming on ectotherm populations is if larger individuals in a population become rarer or smaller (Ohlberger, 2013; Ohlberger et al., 2018), withinspecies mortality and growth responses to warming need to be further studied. Importantly, this requires accounting also for the effects of warming on growth, and how responses in growth and mortality depend on each other. For instance, higher mortality (predation or natural, physiological mortality) can release intraspecific competition and thus increase growth. While e.g., benthic invertebrate density was not affected by the initial warming of the heated area (Sandstrom et al., 1995), warminginduced mortality may have led to higher benthic prey availability per capita for the studied perch. Conversely, altered growth and body sizes can lead to changes in sizespecific mortality, such as predation or starvation, both of which are expected to change with warming (Thunell, 2023). In conclusion, individuallevel patterns such as the TSR can only be used to predict changes in the populationlevel size structure in limited cases, as it does not concern changes in abundanceatsize via mortality. Mortality may, however, be an important driver of the observed shrinking of ectotherms (PeraltaMaraver and Rezende, 2021). Understanding the mechanisms by which the size and agedistribution change with warming is critical for predicting how warming changes species functions and ecological roles (Audzijonyte et al., 2020; Fritschie and Olden, 2016; Gårdmark and Huss, 2020). Our findings demonstrate that a key to do this is to acknowledge temperature effects on both growth and mortality and how they interact.
Materials and methods
Data
We use sizeatage data from perch sampled annually from an artificially heated enclosed bay (‘the Biotest basin’) and its reference area, both in the western Baltic Sea (Figure 1). Heating started in 1980, the first analyzed cohort is 1981, and the first and last catch year is 1987 and 2003, respectively, to omit transient dynamics and acute responses, and to ensure we use cohorts that only experienced one of the thermal environments during its life. A grid at the outlet of the heated area (Figure 1) prevented fish larger than 10 cm from migrating between the areas (Adill et al., 2013; Huss et al., 2019), and genetic studies confirm the reproductive isolation between the two populations during this time period (Björklund et al., 2015). However, the grid was removed in 2004. Since then, fish growing up in the heated Biotest basin can easily swim out and fish caught in the reference area can no longer be assumed to be born there. Hence, we use data only up until 2003. This resulted in 12,658 lengthatage measurements from 2,426 individuals (i.e. multiple measurements per individual) from 256 net deployments.
We use data from fishing events using surveygillnets that took place in October in the heated Biotest basin and in August in the reference area when temperatures are most comparable between the two areas (Huss et al., 2019), because temperature affects catchability in static gears. The catch was recorded by 2.5 cm length classes during 1987–2000, and into 1 cm length groups in years 2001–2003. To express lengths in a common length standard, 1 cm intervals were converted into 2.5 cm intervals. The unit of catch data is hence the number of fish caught by 2.5 cm size class per net per night (i.e. a catchperuniteffort (CPUE) variable). All data from fishing events with disturbance affecting the catch (e.g. seal damage, strong algal growth on the gears, clogging by drifting algae) were removed (years 1996 and 1999 from the heated area in the catch data).
Lengthatage throughout an individuals' life was reconstructed for a random or lengthstratified subset of caught individuals each year (depending on which year, and in some cases, the number of fish caught). This was done using growthincrement biochronologies derived from annuli rings on the operculum bones, with control counts done on otoliths. Such analyses have become increasingly used to analyze changes in the growth and sizeatage of fishes (Essington et al., 2022; Morrongiello and Thresher, 2015). Specifically, an established powerlaw relationship between the distance of annual rings and fish length was used: $L=\kappa {R}^{s}$, where $L$ is the length of the fish, $R$ the operculum radius, $\kappa $ the intercept, and $s$ the slope of the line for the regression of logfish length on logoperculum radius from a large reference data set for perch (Thoresson, 1996). Backcalculated lengthatage was obtained from the relationship ${L}_{a}={L}_{s}{\left(\frac{{r}_{a}}{R}\right)}^{s}$, where ${L}_{a}$ is the backcalculated body length at age $a$, ${L}_{s}$ is the final body length (body length at catch), ${r}_{a}$ is the distance from the center to the annual ring corresponding to age $a$ and $s=0.861$ for perch (Thoresson, 1996). Since perch exhibits sexual sizedimorphism, and age determination together with backcalculation of growth was not done for males in all years, we only used females for our analyses.
Statistical analysis
Request a detailed protocolThe differences in sizeatage, growth, mortality, and size structure between perch in the heated and the reference area were quantified using hierarchical linear and nonlinear models fitted in a Bayesian framework. First, we describe each statistical model and then provide details of model fitting, model diagnostics, and comparison.
To describe individual growth throughout life, we fit the von Bertalanffy growth equation (VBGE) (Beverton and Holt, 1957; von Bertalanffy, 1938) on a $\mathrm{l}\mathrm{o}\mathrm{g}$ scale, describing length as a function of age to evaluate differences in sizeatage and asymptotic size: $\mathrm{l}\mathrm{o}\mathrm{g}\left({L}_{t}\right)=\mathrm{l}\mathrm{o}\mathrm{g}\left({L}_{\mathrm{\infty}}\left(1{e}^{\left(K\left(t{t}_{0}\right)\right)}\right)\right)$, where ${L}_{t}$ is the lengthatage $t$ (years), ${L}_{\infty}$ is the asymptotic size, $K$ is the Brody growth coefficient (${yr}^{1}$) and ${t}_{0}$ is the age when the average length was zero. Here and henceforth, log refers to natural logarithms. We used only age and size at catch, i.e. not backcalculated lengthatage. This was to have a simpler model and not have to account for parameters varying within individuals as well as cohorts, as mean sample size per individual was only ~5. We let parameters vary among cohorts rather than year of catch, because individuals within cohorts share similar environmental conditions and density dependence (Morrongiello and Thresher, 2015). Eight models in total were fitted with area dummycoded, with different combinations of shared and areaspecific parameters. We evaluated if models with areaspecific parameters led to better fit and quantified the differences in areaspecific parameters (indexed by subscripts $\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}$ and $\mathrm{r}\mathrm{e}\mathrm{f}$). The model with all areaspecific parameters can be written as:
where log lengths are $\mathrm{S}\mathrm{t}\mathrm{u}\mathrm{d}\mathrm{e}\mathrm{n}\mathrm{t}\text{}\mathrm{t}$ distributed to account for extreme observations, $\upsilon $, $\mu$, and $\sigma$ represent the degrees of freedom, mean, and the scale parameter, respectively. ${A}_{\mathrm{r}\mathrm{e}\mathrm{f}}$ and ${A}_{\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}$ are dummy variables such that ${A}_{\mathrm{r}\mathrm{e}\mathrm{f}}=1$ and ${A}_{\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}=0$ if it is the reference area, and vice versa for the heated area. The multivariate normal distribution in Equation 3 is the prior for the cohortvarying parameters ${L}_{\infty \mathrm{r}\mathrm{e}\mathrm{f}j}$, ${L}_{\infty \mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}j}$, ${K}_{\mathrm{r}\mathrm{e}\mathrm{f}j}$, and ${K}_{\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}j}$ (for cohorts $j=1981,$…,1997) (note that cohorts extend further back in time than the catch data), with hyperparameters ${\mu}_{{L}_{\infty \mathrm{r}\mathrm{e}\mathrm{f}}}$, ${\mu}_{{L}_{\infty \mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}}$, ${\mu}_{{K}_{\mathrm{r}\mathrm{e}\mathrm{f}}}$, ${\mu}_{{K}_{\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}}$ describing the population means and a covariance matrix with the betweencohort variation along the diagonal. We did not model a correlation between the parameters, hence offdiagonals are 0. The other seven models include some or all parameters as parameters common for the two areas, e.g., substituting ${L}_{\infty \mathrm{r}\mathrm{e}\mathrm{f}j}$ and ${L}_{\infty \mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}j}$ with ${L}_{\infty j}$. To aid the convergence of this nonlinear model, we used informative priors chosen after visualizing draws from prior predictive distributions (Wesner and Pomeranz, 2021) using probable parameter values (Figure 2—figure supplement 7; Figure 3—figure supplement 3). We used the same prior distribution for each parameter class for both areas to not introduce any other sources of differences in parameter estimates between areas. We used the following priors for the VBGE model: ${\mu}_{{L}_{\infty \mathrm{r}\mathrm{e}\mathrm{f},\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}}~\mathrm{N}\left(\mathrm{45,20}\right)$ , ${\mu}_{{K}_{\mathrm{r}\mathrm{e}\mathrm{f},\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}}~\mathrm{N}\left(0.2,0.1\right)$ , ${t}_{0\mathrm{r}\mathrm{e}\mathrm{f},\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}~\mathrm{N}\left(0.5,1\right)$, and $\upsilon ~\mathrm{g}\mathrm{a}\mathrm{m}\mathrm{m}\mathrm{a}\left(\mathrm{2,0.1}\right)$. $\sigma $ parameters, $\sigma}_{{L}_{\mathrm{\infty}\mathrm{r}\mathrm{e}\mathrm{f}}$, $\sigma}_{{L}_{\mathrm{\infty}\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}$, $\sigma}_{{K}_{\mathrm{r}\mathrm{e}\mathrm{f}}$, $\sigma}_{{K}_{\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}$ were given a $\mathrm{S}\mathrm{t}\mathrm{u}\mathrm{d}\mathrm{e}\mathrm{n}\mathrm{t}\text{}t\left(3,0,2.5\right)$ prior.
We also compared how body growth scales with body size (in contrast to length vs age). This is because sizeatage reflects lifetime growth history rather than current growth and may thus be large because growth was fast early in life, not because current growth rates are fast (Lorenzen, 2016). We therefore fit allometric growth models describing how specific growth rate scales with length: $G=\alpha {L}^{\theta}$ , where $G$, the annual specific growth between year $t$ and $t+1$, is defined as: $G=100\times \left(\mathrm{l}\mathrm{o}\mathrm{g}\left({L}_{t+1}\right)\mathrm{l}\mathrm{o}\mathrm{g}\left({L}_{t}\right)\right)$ and $L$ is the geometric mean length: $L={\left({L}_{t+1}\times {L}_{t}\right)}^{0.5}$. Here we use backcalculated lengthatage, resulting in multiple observations per individual. As with the VBGE model, we dummycoded area to compare models with different combinations of common and shared parameters. We assumed growth rates were $\mathrm{S}\mathrm{t}\mathrm{u}\mathrm{d}\mathrm{e}\mathrm{n}\mathrm{t}\text{}\mathrm{t}$ distributed, and the full model can be written as:
We assumed only $\alpha $ varied across individuals $j$ within cohorts $k$ and compared two models: one with $\theta $ common for the heated and reference area, and one with an areaspecific $\theta $. We used the following priors: ${\mu}_{{\alpha}_{\mathrm{r}\mathrm{e}\mathrm{f},\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}}\sim \mathrm{N}\left(500,\text{}100\right)$, ${\theta}_{\mathrm{r}\mathrm{e}\mathrm{f},\mathrm{h}\mathrm{e}\mathrm{a}\mathrm{t}}~\mathrm{N}(1.2,0.3)$ and $\upsilon ~\mathrm{g}\mathrm{a}\mathrm{m}\mathrm{m}\mathrm{a}\left(\mathrm{2,0.1}\right)$. $\sigma $, ${\sigma}_{\mathrm{i}\mathrm{d}:\mathrm{c}\mathrm{o}\mathrm{h}\mathrm{o}\mathrm{r}\mathrm{t}}$ and ${\sigma}_{\mathrm{c}\mathrm{o}\mathrm{h}\mathrm{o}\mathrm{r}\mathrm{t}}$ were all given a $\mathrm{S}\mathrm{t}\mathrm{u}\mathrm{d}\mathrm{e}\mathrm{n}\mathrm{t}\text{}t\left(3,0,13.3\right)$ prior.
We estimated total mortality by fitting linear models to the natural log of catch (CPUE) as a function of age (catch curve regression), under the assumption that in a closed population, the exponential decline can be described as ${N}_{t}={N}_{0}{e}^{Zt}$ , where ${N}_{t}$ is the population at time $t$, ${N}_{0}$ is the initial population size and $Z$ is the instantaneous mortality rate. This equation can be rewritten as a linear equation: $\mathrm{l}\mathrm{o}\mathrm{g}\left({C}_{t}\right)=\mathrm{l}\mathrm{o}\mathrm{g}\left(v{N}_{0}\right)Zt$, where ${C}_{t}$ is a catch at age $t$, if a catch is assumed proportional to the number of fish (i.e. ${C}_{t}=v{N}_{t}$). Hence, the negative of the age slope is the mortality rate, $Z$. To get catchatage data, we constructed areaspecific agelength keys using the subsample of the total (female) catch that was agedetermined. Age length keys describe the age proportions of each length category (i.e. a matrix with length category as rows, and ages as columns). The age composition is then estimated for the total catch based on the ‘probability’ of fish in each length category being a certain age. Due to the smallest and youngest fish not being representatively caught with the gillnet, the catch is domeshaped over size and age. We therefore followed the practice of selecting only ages on the descending right limb (Dunn et al., 2002; Figure 4—figure supplement 1). We fit this model with and without an $age\times area$interaction, and the former can be written as:
where ${\beta}_{0j}$ and ${\beta}_{1j}$ are the intercepts for the reference and heated areas, respectively, ${\beta}_{2j}$ is the age slope for the reference area and ${\beta}_{3j}$ is the difference between the age slope in the reference area and in the heated area. All parameters vary by cohort (for cohort $j=1981,\dots ,2000$). We use the default brms priors for these models, i.e., flat priors for the regression coefficients (Bürkner, 2017) and $\upsilon ~\mathrm{g}\mathrm{a}\mathrm{m}\mathrm{m}\mathrm{a}\left(\mathrm{2,0.1}\right)$. $\sigma $ and $\sigma}_{{\beta}_{0j,\dots ,3j}$ were given a $\mathrm{S}\mathrm{t}\mathrm{u}\mathrm{d}\mathrm{e}\mathrm{n}\mathrm{t}\text{}t\left(3,0,2.5\right)$ prior.
Lastly, we quantified differences in the average age and size distributions between the areas. We estimate the biomass sizespectrum exponent $\gamma $ directly, using the likelihood approach for binned data, i.e., the MLEbin method in the R package sizeSpectra (Edwards, 2020; Edwards et al., 2020; Edwards et al., 2017). This method explicitly accounts for uncertainty in body masses within size classes (bins) in the data and has been shown to be less biased than regressionbased methods or the likelihood method based on bin midpoints (Edwards et al., 2020; Edwards et al., 2017). We pooled all years to ensure negative relationships between biomass and size in the size classes (as the sign of the relationship varied between years). We also fitted lognormal models as data are positive and tailed to length and ageresolved catch data. Here, we assume that the catchability with respect to size does not differ between the areas, and, therefore, use the entire catch (Figure 6—figure supplement 1). In contrast to the catch curve regression, we do not need to filter representatively caught size or age classes. The lognormal models fitted to age or size (denoted ${y}_{age,length,i}$) model can be written as:
where ${\beta}_{0j}$ is the intercept for the reference area and ${\beta}_{1j}$ is the intercept for the heated area. These intercepts vary by year (for years $j=1987,\dots ,2003$). We use flat priors for the regression coefficients, and $\sigma $ was given a $\mathrm{S}\mathrm{t}\mathrm{u}\mathrm{d}\mathrm{e}\mathrm{n}\mathrm{t}\text{}t\left(3,0,2.5\right)$ prior, and compared models with and without random slopes.
All analyses were done using R (R Development Core Team, 2020) version 4.0.2 with R Studio (2021.09.1). The packages within the tidyverse (Wickham et al., 2019) collection was used to process and visualize data. Models were fit using the R package brms (Bürkner, 2018). For the nonlinear von Bertalanffy growth equation and the allometric growth model, we used informative priors to facilitate convergence. These were chosen by defining vague priors, and then progressively tightening these until convergence was achieved (Bürkner, 2017; Gesmann and Morris, 2020). We used prior predictive checks to ensure the priors were suitable (vague enough to include also unlikely predictions, but informative enough to ensure convergence), and the final prior predictive checks are shown in Figure 2—figure supplement 1 and Figure 3—figure supplement 1. We also explored priors vs posteriors to evaluate the influence of our informative priors visually (Figure 2—figure supplement 7; Figure 3—figure supplement 3). For the linear models (catch curve and mean size), which do not require the same procedure to achieve convergence typically, we used the default priors from brms as written above. We used three chains and 4000 iterations in total per chain. Models were compared by evaluating their expected predictive accuracy (expected log pointwise predictive density) using leaveoneout crossvalidation (LOOCV) (Vehtari et al., 2017) while ensuring Pareto $k$ values <0.7, in the R package loo (Vehtari et al., 2020). Results of the model comparison can be found in the Supplementary file 1. We used bayesplot (Gabry et al., 2019) and tidybayes (Kay, 2019) to process and visualize model diagnostics and posteriors. Model convergence and the fit were assessed by ensuring potential scale reduction factors ($\hat{R}$) were less than 1.1, suggesting all three chains converged to a common distribution (Gelman et al., 2003), and by visually inspecting trace plots, residuals QQplots, and with posterior predictive checks (Figures 2—4, Figure 6—figure supplement 2).
Data availability
All raw data and R code to clean and reproduce the results reported in this paper are available on GitHub (copy archived at Lindmark, 2023) and have been deposited on Zenodo.
References

ReportBiologisk recipientkontroll vid forsmarks kärnkraftverk (in Swedish)Öregrund.

Sizebased theory for fisheries adviceICES Journal of Marine Science 77:2445–2455.https://doi.org/10.1093/icesjms/fsaa157

Temperature and organism size—A biological law for ectotherms?Advances in Ecological Research 25:1–58.https://doi.org/10.1016/S00652504(08)602123

Why are organisms usually bigger in colder environments? Making sense of a life history puzzleTrends in Ecology & Evolution 12:235–239.https://doi.org/10.1016/s01695347(97)010586

Fish body sizes change with temperature but not all species shrink with warmingNature Ecology & Evolution 4:809–814.https://doi.org/10.1038/s4155902011710

Warming increases the cost of growth in a model vertebrateFunctional Ecology 33:1256–1266.https://doi.org/10.1111/13652435.13348

Oldgrowth fishes become scarce under fishingCurrent Biology 27:2843–2848.https://doi.org/10.1016/j.cub.2017.07.069

Warming temperatures and smaller body sizes: synchronous changes in growth of North Sea fishesGlobal Change Biology 20:1023–1031.https://doi.org/10.1111/gcb.12514

Rapid warming of large marine ecosystemsProgress in Oceanography 81:207–213.https://doi.org/10.1016/j.pocean.2009.04.011

Warmer water increases early body growth of Northern pike ( esox lucius ), but mortality has larger impact on decreasing body sizesCanadian Journal of Fisheries and Aquatic Sciences 79:771–781.https://doi.org/10.1139/cjfas20200386

Nonlinear averaging of thermal experience predicts population growth rates in a thermally variable environmentProceedings. Biological Sciences 285:20181076.https://doi.org/10.1098/rspb.2018.1076

On the dynamics of exploited fish populationsFishery Investigations London Series 19:21064.https://doi.org/10.1007/9789401121064

Do climate and fishing influence sizebased indicators of Celtic Sea fish community structure?ICES Journal of Marine Science 62:405–411.https://doi.org/10.1016/j.icesjms.2005.01.006

From bacteria to whales: using functional size spectra to model marine ecosystemsTrends in Ecology & Evolution 32:174–186.https://doi.org/10.1016/j.tree.2016.12.003

Brms: an R package for bayesian multilevel models using stanJournal of Statistical Software 80:i01.https://doi.org/10.18637/jss.v080.i01

Shrinking of fishes exacerbates impacts of global ocean changes on marine ecosystemsNature Climate Change 3:254–258.https://doi.org/10.1038/nclimate1691

Testing and recommending methods for fitting size spectra to dataMethods in Ecology and Evolution 8:57–67.https://doi.org/10.1111/2041210X.12641

Accounting for the bin structure of data removes bias when fitting size spectraMarine Ecology Progress Series 636:19–33.https://doi.org/10.3354/meps13230

Fitting growth models to otolith increments to reveal timevarying growthCanadian Journal of Fisheries and Aquatic Sciences 79:159–167.https://doi.org/10.1139/cjfas20210046

Visualization in Bayesian workflowJournal of the Royal Statistical Society Series A 182:389–402.https://doi.org/10.1111/rssa.12378

Individual variation and interactions explain food web responses to global warmingPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 375:20190449.https://doi.org/10.1098/rstb.2019.0449

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

Size, growth, temperature and the natural mortality of marine fishFish and Fisheries 11:149–158.https://doi.org/10.1111/j.14672979.2009.00350.x

Climate change impacts on marine ecosystems through the lens of the size spectrumEmerging Topics in Life Sciences 3:233–243.https://doi.org/10.1042/ETLS20190042

An application of a bioenergetics model to Eurasian perch (Perca fluviatilis L)Journal of Fish Biology 41:217–230.https://doi.org/10.1111/j.10958649.1992.tb02652.x

Optimum growth temperature declines with body size within fish speciesGlobal Change Biology 28:2259–2271.https://doi.org/10.1111/gcb.16067

Predicting natural mortality rates of plants and animalsEcology Letters 11:710–716.https://doi.org/10.1111/j.14610248.2008.01190.x

Demographic changes in Chinook salmon across the Northeast Pacific OceanFish and Fisheries 19:533–546.https://doi.org/10.1111/faf.12272

Warming leads to opposite patterns in weightatage for young versus old age classes of Bering sea walleye pollockCanadian Journal of Fisheries and Aquatic Sciences 79:1655–1666.https://doi.org/10.1139/cjfas20210315

On the interrelationships between natural mortality, growth parameters, and mean environmental temperature in 175 fish stocksICES Journal of Marine Science 39:175–192.https://doi.org/10.1093/icesjms/39.2.175

The gilloxygen limitation theory (GOLT) and its criticsScience Advances 7:eabc6050.https://doi.org/10.1126/sciadv.abc6050

Heat tolerance in ectotherms scales predictably with body sizeNature Climate Change 11:58–63.https://doi.org/10.1038/s4155802000938y

What we don't know about the effects of temperature on salmonid growthTransactions of the American Fisheries Society 151:3–12.https://doi.org/10.1002/tafs.10338

SoftwareR: A language and environment for statistical computingR Foundation for Statistical Computing, Vienna, Austria.

Growth, temperature, and density relationships of North Sea cod (Gadus morhua)Canadian Journal of Fisheries and Aquatic Sciences 65:456–470.https://doi.org/10.1139/f07150

Effects of temperature on life history variables in perchJournal of Fish Biology 47:652–670.https://doi.org/10.1111/j.10958649.1995.tb01932.x

The production of particles in the surface waters of the ocean with particular reference to the Sargasso SeaLimnology and Oceanography 18:719–733.https://doi.org/10.4319/lo.1973.18.5.0719

Shrinking body size as an ecological response to climate changeNature Climate Change 1:401–406.https://doi.org/10.1038/nclimate1259

BookMetoder För Övervakning Av Kustfiskbestånd in SwedishKustlaboratoriet, Fiskeriverket.

Predicting life history parameters for all fishes worldwideEcological Applications 27:2262–2276.https://doi.org/10.1002/eap.1606

ThesisFish life histories in a warming climate: a mechanistic basis of change and a community contextSwedish University of Agricultural Sciences, Uppsala, Sweden.https://doi.org/10.54612/a.61rases4mb

A mathematical model of some aspects of fish growth, respiration, and mortalityJournal of the Fisheries Research Board of Canada 24:2355–2453.https://doi.org/10.1139/f67190

Implications of lateinlife densitydependent growth for fishery sizeatentry leading to maximum sustainable yieldICES Journal of Marine Science 75:1296–1305.https://doi.org/10.1093/icesjms/fsx236

Practical Bayesian model evaluation using leaveoneout crossvalidation and WAICStatistics and Computing 27:1413–1432.https://doi.org/10.1007/s1122201696964

A quantitative theory of organic growth (inquiries on growth lawsHuman Biology 10:181–213.

Relationships between body size and abundance in ecologyTrends in Ecology & Evolution 22:323–330.https://doi.org/10.1016/j.tree.2007.03.007

Welcome to the tidyverseJournal of Open Source Software 4:1686.https://doi.org/10.21105/joss.01686

Smaller adult fish size in warmer water is not explained by elevated metabolismEcology Letters 25:1177–1188.https://doi.org/10.1111/ele.13989
Decision letter

David A DonosoReviewing Editor; Escuela Politécnica Nacional, Ecuador

Meredith C SchumanSenior Editor; University of Zurich, Switzerland

Christine StawitzReviewer
Our editorial process produces two outputs: (i) public reviews designed to be posted alongside the preprint for the benefit of readers; (ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.
Decision letter after peer review:
Thank you for submitting your article "Faster growth rates and higher mortality but similar sizespectrum in heated, largescale natural experiment" for consideration by eLife. Your article has been reviewed by 2 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Meredith Schuman as the Senior Editor. The following individual involved in the review of your submission has agreed to reveal their identity: Christine Stawitz (Reviewer #1).
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
Essential revisions:
Both reviewers raised questions about the validity/meaning of CPUE values for different treatments. Please revise in response to these comments and amend the text when possible.
Reviewer #1 (Recommendations for the authors):
Congratulations on an excellent paper – I really enjoyed reading it.
58: Is this a common way of measuring changes in size? Looking strictly at mean body size (w/o age structure) seems simplistic – but a citation would convince me this is a common practice
240 – 244: There was a bit more information in your Github repository – it looks like you started with the default priors then had to tighten them to achieve convergence in some of the VBGF models? I agree this is too much detail to include in the main body of the paper but for future scientists' use, I would recommend adding a supplementary section describing how you did this and a table with the priors you used so readers can find them without going to Github.
376: please add the full citation.
Figure 4A: Does it mean anything that CPUE has a lower rate at the beginning of the time series in the heated area than the unheated area? Or is the slope the only thing that matters?
Equations 1 12: I find it confusing to reuse mu, eta, and σ throughout the three different modeling approaches and would prefer to see a unique letter applied for each of the three models to improve clarity.
Reviewer #2 (Recommendations for the authors):
I'd like to see more consideration of faster life history vs higher mortality in this paper. Have these fish adapted a faster life history in response to warmed conditions that happen to result in higher mortality because of the demands of the faster pace of life? Or, given that adults did not shrink, is the increased mortality due to some other factor, like increased predation or stress from the heat? Perhaps this cannot be answered in this system, but it seems like a worthwhile consideration, especially given the unexpected results.
I did have some more substantial questions on the CPUE section, which I'll try to put together here. It seems that CPUE is always lower for the heated population. Could the unexpected increased size observed in heated adults be due to a tradeoff between slower growth due to temperature but higher growth due to a release from density dependence? I think the authors suggest this might be the case, but it would be helpful to state it more explicitly.
Can the authors explain how the initial population size is known and whether it differs between the heated and reference population? It seems there are many reasons why reproductive potential could differ between the heated and reference populations, especially if warming shifts reproductive investment (as would be expected). Why does the CPUE figure start at age 3? It seems younger fish were caught and the fishbase entry on perch says they can mature earlier than 3. Is the mortality rate they measure only applied to fish that have been 'recruited' to the survey? And if fish can mature younger than 3, are some being missed by looking at CPUE for 3+ only?
Finally, the grammar could use some editing, although the errors do not impact readability. I've made a few suggestions below but did not point out all places where grammar could be improved.
https://doi.org/10.7554/eLife.82996.sa1Author response
Reviewer #1 (Recommendations for the authors):
Congratulations on an excellent paper – I really enjoyed reading it.
58: Is this a common way of measuring changes in size? Looking strictly at mean body size (w/o age structure) seems simplistic – but a citation would convince me this is a common practice
Thank you. We think statements/predictions about shrinking (without specific mentioning of the age and/or sizestructure) often is made when talking about “large scale” or general shrinking of organisms as a universal rule. The temperaturesize rule (TSR) on the other hand, which we describe in more detail further down (line 68) describes changes in size over ontogeny, and hence does not make predictions about the sizedistribution of organisms or populations, which is also affected by other processes such as mortality. The reason we open the introduction with the universal prediction about shrinking (and the papers we cite here do focus on average shrinking of species or taxa, to address your specific question) is to make it clear that the predictions about organism shrinking with warming is based on various, sometimes related, biogeographical pattens across space and time, and on the other hand, individuallevel experimental rules. Therefore, we think there is a need to distinguish these patterns and rules because there are situations where they do not conform. The example we give is that the TSR predicts larger sizeatage for young specimens, meaning that mortality (and/or shrinking of adults) must increase to offset that sizeincrease for populationlevel shrinking to occur.
We have edited the second sentence (lines 64–66) to make it clearer.
240 – 244: There was a bit more information in your Github repository – it looks like you started with the default priors then had to tighten them to achieve convergence in some of the VBGF models? I agree this is too much detail to include in the main body of the paper but for future scientists' use, I would recommend adding a supplementary section describing how you did this and a table with the priors you used so readers can find them without going to Github.
Thank you for pointing out that need. Perhaps a detail, but for the record, brms forces the user to specify priors in nonlinear models, so we never started with the standard default priors that are used for general regression parameters. Instead, we first defined broad priors centered on what we believed were reasonable values for this fish species, and then progressively tightened them until convergence was achieved. We have added text in the model fitting paragraph (lines 330–338) describing the iterative approach to finding priors for the nonlinear model and removed the single sentence describing this approach under each model description and a reference. All priors are given in the main text under each model with inline equations, so perhaps it is not needed to put them also in a table.
376: please add the full citation.
Done
Figure 4A: Does it mean anything that CPUE has a lower rate at the beginning of the time series in the heated area than the unheated area? Or is the slope the only thing that matters?
We are not sure what is meant by “a lower rate of CPUE (catch per unit effort) in the beginning of the time series”, but we suspect the reviewer means that the catch rate is lower in the heated area on average (since Figure 4 depicted data and the global prediction without the random cohort effect). For mortality estimates, the slope is what matters, because it relates to the instantaneous mortality parameter in a population with exponential decay when log transformed, whereas the intercept does not matter directly for the mortality estimate. Note also that we have added a model for the average age, which is complementary to the analysis of mortality but is conceptually a different approach to showing the same thing: that the heated population consists of younger fish. As also Reviewer#2 pointed out, we do observe a difference in CPUE between the areas (i.e., the intercept), but we can only speculate why this is (possibly linked to size of the habitat).
Equations 1 12: I find it confusing to reuse mu, eta, and σ throughout the three different modeling approaches and would prefer to see a unique letter applied for each of the three models to improve clarity.
We thought about this but decided that it would probably be clearer to have a unique letter for the response variable, and a general definition for the parameters of the distribution, e.g., L_{i} for length and G_{i} for growth and mu for the mean of the distribution in question. The alternative would be to add more subscripts, but we have already 3 levels of subscripts in some models, and adding another level make it difficult to read. Instead, we hope the structure of the methods (one paragraph per model) the variable spelled out or abbreviated makes it clear that these are separate models.
Reviewer #2 (Recommendations for the authors):
I'd like to see more consideration of faster life history vs higher mortality in this paper. Have these fish adapted a faster life history in response to warmed conditions that happen to result in higher mortality because of the demands of the faster pace of life? Or, given that adults did not shrink, is the increased mortality due to some other factor, like increased predation or stress from the heat? Perhaps this cannot be answered in this system, but it seems like a worthwhile consideration, especially given the unexpected results.
It is difficult to say with the data we have available, and it becomes speculative. Since there are other fish predators (mainly pike and larger perch), which also need to feed at higher rates in the warmer water, it seems inevitable that higher predation mortality to some degree contributes to the higher mortality in perch. However, we do not know nor have any way of testing with the data we have available how important this is compared to e.g., physiological stress, ageing and paceoflife responses in the heated area experience. We do agree however, that it is an important point, and now bring forth this in the introduction, on lines 112–113, and in the discussion, line 497.
I did have some more substantial questions on the CPUE section, which I'll try to put together here. It seems that CPUE is always lower for the heated population. Could the unexpected increased size observed in heated adults be due to a tradeoff between slower growth due to temperature but higher growth due to a release from density dependence? I think the authors suggest this might be the case, but it would be helpful to state it more explicitly.
Yes, we agree that could be part of the explanation. As we replied to Reviewer#1, the issue is that the effect of living in a smaller area (essentially a small lake compared to the open coast, which could affect density), is largely confounded with the temperature effect. That said, there exists data from a few years prior to the warm water pollution in the heated basin (but after its construction), that indicate the CPUE was always lower in the basin (Huss et al., 2019) (doi: 10.1111/gcb.14637). The lower CPUE could therefore be an artefact of the enclosing of the basin, possibly related to its size. However, in those first years of the enclosing, but before the warm water pollution started, growth rates had not yet differentiated (Huss et al. 2019). However, the effect of density also depends on food availability. While prey data are scarce, Sandström et al. (1995) (doi: 10.1111/j.10958649.1995.tb01932.x) show that the density of benthic prey was relatively constant over time (including 4 years prior to warming), suggesting that it is the warming itself that drives the difference in size. We have expanded on this in the discussion (lines 516–518).
Can the authors explain how the initial population size is known and whether it differs between the heated and reference population? It seems there are many reasons why reproductive potential could differ between the heated and reference populations, especially if warming shifts reproductive investment (as would be expected). Why does the CPUE figure start at age 3? It seems younger fish were caught and the fishbase entry on perch says they can mature earlier than 3. Is the mortality rate they measure only applied to fish that have been 'recruited' to the survey? And if fish can mature younger than 3, are some being missed by looking at CPUE for 3+ only?
First, we just want to clarify that the initial size here refers to N_{0}, the recruitment or abundance at age 0, not the first population size in time. This is not directly observed but is the intercept in the catch curve regression. It likely differs between the areas for many reasons: e.g., the size of the system (potentially affecting the overall abundance and density of the heated population), and differences in reproductive investment due to life history optimization in warmer waters. However, with the latter it is less known if this also translates to a difference in actual “recruitment” to the fishable population, i.e., the abundance at age 3, due to unknown mortality rates early in life.
It is true that fish younger than age 3 were caught. But age 3 is the age when the catches start to descend, and it is thus assumed that only fish above this age are caught representatively (i.e., that catches of age 2 fish are lower is not because these individuals are rarer in the population, but because they are not caught as effectively in the sizeselective gill nets). Therefore, using fish older than 2 years is done to ensure that the catchability of certain ages does not affect the slope estimate, it should only be the true catch rate that does that. Ideally, we would want to avoid filtering the catch data like this, but this is the only way to handle the age/sizebased catchability of the fishing gear, and it is standard practice in catch curve regression (Dunn et al., 2002) (doi: https://doi.org/10.1016/S01657836(01)004076).
We have added a figure illustrating the domeshaped relationship between CPUE and age that is due to both catchability and declines in abundance by age (Figure 4—figure supplement 1), and we describe this procedure in the main text now (lines 280–283 and see also lines 310–312).
Finally, the grammar could use some editing, although the errors do not impact readability. I've made a few suggestions below but did not point out all places where grammar could be improved.
Thank you for pointing that out. We have now checked and corrected errors throughout the manuscript.
https://doi.org/10.7554/eLife.82996.sa2Article and author information
Author details
Funding
No external funding was received for this work.
Acknowledgements
We thank all staff involved in data collection, Jens Olsson and Göran Sundblad for discussions, Christine Stawitz and an anonymous reviewer for feedback that greatly improved the manuscript, and Forsmark Kraftgrupp AB for making data publicly available. This study was supported by SLU Quantitative Fish and Fisheries Ecology.
Senior Editor
 Meredith C Schuman, University of Zurich, Switzerland
Reviewing Editor
 David A Donoso, Escuela Politécnica Nacional, Ecuador
Reviewer
 Christine Stawitz
Version history
 Preprint posted: April 14, 2022 (view preprint)
 Received: August 25, 2022
 Accepted: April 10, 2023
 Version of Record published: May 9, 2023 (version 1)
Copyright
© 2023, Lindmark et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 695
 Page views

 79
 Downloads

 2
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading

 Ecology
 Evolutionary Biology
Groups of animals inhabit vastly different sensory worlds, or umwelten, which shape fundamental aspects of their behaviour. Yet the sensory ecology of species is rarely incorporated into the emerging field of collective behaviour, which studies the movements, populationlevel behaviours, and emergent properties of animal groups. Here, we review the contributions of sensory ecology and collective behaviour to understanding how animals move and interact within the context of their social and physical environments. Our goal is to advance and bridge these two areas of inquiry and highlight the potential for their creative integration. To achieve this goal, we organise our review around the following themes: (1) identifying the promise of integrating collective behaviour and sensory ecology; (2) defining and exploring the concept of a ‘sensory collective’; (3) considering the potential for sensory collectives to shape the evolution of sensory systems; (4) exploring examples from diverse taxa to illustrate neural circuits involved in sensing and collective behaviour; and (5) suggesting the need for creative conceptual and methodological advances to quantify ‘sensescapes’. In the final section, (6) applications to biological conservation, we argue that these topics are timely, given the ongoing anthropogenic changes to sensory stimuli (e.g. via light, sound, and chemical pollution) which are anticipated to impact animal collectives and grouplevel behaviour and, in turn, ecosystem composition and function. Our synthesis seeks to provide a forwardlooking perspective on how sensory ecologists and collective behaviourists can both learn from and inspire one another to advance our understanding of animal behaviour, ecology, adaptation, and evolution.

 Ecology
 Plant Biology
Global agrobiodiversity has resulted from processes of plant migration and agricultural adoption. Although critically affecting current diversity, crop diffusion from Classical antiquity to the Middle Ages is poorly researched, overshadowed by studies on that of prehistoric periods. A new archaeobotanical dataset from three Negev Highland desert sites demonstrates the first millennium CE&'s significance for longterm agricultural change in southwest Asia. This enables evaluation of the 'Islamic Green Revolution' (IGR) thesis compared to 'Roman Agricultural Diffusion' (RAD), and both versus crop diffusion during and since the Neolithic. Among the finds, some of the earliest aubergine (Solanum melongena) seeds in the Levant represent the proposed IGR. Several other identified economic plants, including two unprecedented in Levantine archaeobotanyjujube (Ziziphus jujuba/mauritiana) and white lupine (Lupinus albus)implicate RAD as the greater force for crop migrations. Altogether the evidence supports a gradualist model for Holocenewide crop diffusion, within which the first millennium CE contributed more to global agricultural diversity than any earlier period.