An Individual-Based Model of the Red Alga Agarophyton chilense Unravels the Complex Demography of Its Intertidal Stands

Abstract

Algal demographic models have been developed mainly to study their life cycle evolution or optimize their commercial exploitation. Most commonly, structured-aggregated population models simulate the main life cycle stages considering their fertility, growth and survival. Their coarse resolution results in weak predictive abilities since neglected details may still impact the whole. In our case, we need a model of Agarophyton chilense natural intertidal populations that unravels the complex demography of isomorphic biphasic life cycles and be further used for: (i) introduction of genetics, aimed at studying the evolutionary stability of life cycles, (ii) optimizing commercial exploitation, and (iii) adaptation for other species. Long-term monitoring yield 6,066 individual observations and 40 population observations. For a holistic perspective, we developed an Individual-Based Model (IBM) considering ploidy stage, sex stage, holdfast age and survival, frond size, growth and breakage, fecundity, spore survival, stand biomass, location and season. The IBM was calibrated and validated comparing observed and estimated sizes and abundances of gametophyte males, gametophyte females and tetrasporophytes, stand biomass, haploid:dipoid ratio (known as H:D or G:T), fecundity and recruitment. The IBM replicated well the respective individual and population properties, and processes such as winter competition for light, self-thinning, summer stress from desiccation, frond breakage and re-growth, and different niche occupation by haploids and diploids. Its success depended on simulating with precision details such as the holdfasts’ dynamics. Because “details” often occur for a reduced number of individuals, inferring about them required going beyond statistically significant evidences and integrating these with parameter calibration aimed at maximized model fit. On average, the population was haploid-dominated (H:D > 1). In locations stressed by desiccation, the population was slightly biased toward the diploids and younger individuals due to the superior germination and survival of the diploid sporelings. In permanently submerged rock pools the population was biased toward the haploids and older individuals due to the superior growth and survival of the haploid adults. The IBM application demonstrated that conditional differentiation among ploidy stages was responsible for their differential niche occupation, which, in its turn, has been argued as the driver of the evolutionary stability of isomorphic biphasic life cycles.