DOI: 10.5061/dryad.7n474b1
title: Data from: Warming impacts on early life stages increase the vulnerability and delay the population recovery of a long-lived habitat-forming macroalga

description: 1. Understanding the combined effects of global and local stressors is
 crucial for conservation and management, yet challenging due to the
 different scales at which these stressors operate. Here we examine the
 effects of one of the most pervasive threats to marine biodiversity, ocean
 warming, on the early life stages of the habitat-forming macroalga
 Cystoseira zosteroides, its long-term consequences for population
 resilience and its combined effect with physical stressors. 2. First, we
 performed a controlled laboratory experiment exploring the impacts of
 warming on early life stages. Settlement and survival of germlings were
 measured at 16ºC (control), 20ºC and 24ºC and both processes were affected
 by increased temperatures. Then, we integrated this information into
 stochastic, density-dependent integral projection models (IPM). 3.
 Recovery time after a minor disturbance significantly increased in warmer
 scenarios. The stochastic population growth rate (λs) was not strongly
 affected by warming alone, as high adult survival compensated for
 thermal-induced recruitment failure. Nevertheless, warming coupled with
 recurrent physical disturbances had a strong impact on λs and population
 viability. 4. Synthesis: The impact of warming effects on early stages may
 significantly decrease the natural ability of habitat-forming algae to
 rebound after major disturbances. These findings highlight that, in a
 global warming context, populations of deep-water macroalgae will become
 more vulnerable to further disturbances, and stress the need to
 incorporate abiotic interactions into demographic models.
description: CzIPMR code to estimate the
 recovery time for Cystoseira zosteroides populations after a major
 disturbance at different temperature scenarios treatments. In addition,
 stochastic population growth rate (λs) and quasi-extinction probability at
 increasing frequency of two major disturbances at increasing temperature
 scenarios. These analyses correspond to the figures 4 and 5 of Capdevila
 et al. 2018 JEcol.MixedEffectsparamsParameter values needed for the Integral Projection Models used to model the life cycle and population dynamics of Cystoseira zosteroides. This includes seven demographic processes: 1.survival (σ), 2.growth (γ), 3.fertility (φ), 4.recruits per capita (δ(N)), 5.probability of settlement of recruits (ε), 6.early survival of recruits (σs) and 7.recruits size probability distribution.IPMFunctionsFunctions required to run the CzIPM.R script. This script contains the description of the growth, survival and fecundity functions used to build the IPMs.1. The best-fitted model for survival (σ) was a logistic mixed effect model including size as fixed factors and population nested in years as a random factor. 2. For growth (γ), the best-fitted model was a linear mixed effect model, with size as fixed factor and population nested in year as random factor. 3. Fertility (φ(z)), was estimated as the relation between reproductive status (reproductive vs. non-reproductive) and size with a binomial regression. 4. Recruitment per capita (δ(N)) is density-dependent in C. zosteroides (Capdevila et al., 2015), so a generalized linear model with Poisson error distribution and a log-link function was fitted, correlating the recruit:adult ratio as a function of the adult density. 5. To model the effect of temperature on the probability of settlement (ε) we used a generalized linear mixed models (GLMM), with a Poisson error distribution and a logit link function, the independent variable was the number of zygotes, temperature was treated as a fixed variable and we used the ID of each quadrat of the Petri dishes as a random variable. 6. To model the effect of temperature and time (fixed factors) on germling survival (σs), we used a GLMM with a binomial error distribution and a logit link function, with the ID of each quadrat of each Petri dish as a random variable to deal with the lack of independence between observations repeated at different times and a binomial error distribution was assumed to deal with the binary response variable (survive vs. die). 7. The size distribution of recruits was estimated as a normal probability function. In addition, the function required to project the density-dependent and stochastic IPMs is provided.modsumDensity-dependent function, relating the number of Cystoseira zosteroides recruits with the number of adults. It is a generalized linear model (GLM) with Poisson error distribution and a log-link function, correlating the recruit:adult ratio with the adult density. This file is needed to run the code CzIPM.R.settData on the impacts of temperature (16ºC, 20ºC and 24ºC) on the settlement of Cystoseira zosteroides early stages. This file is needed to perform the projections in CzIPM.R code.survrecData on the impacts of the temperature treatments (16ºC, 20ºC and 24ºC) to early survival of Cystoseira zosteroides. This file is required to run the code CzIPM.R.