Nonlinear processes in seagrass colonisation explained by simple clonal growth rules

Understanding the depth limit of the seagrass Cymodocea nodosa as a critical transition: Field and modeling evidence

Changes in light and sediment conditions can sometimes trigger abrupt regime shifts in seagrass meadows resulting in dramatic and unexpected die-offs of seagrass. Light attenuates rapidly with depth, and in seagrass systems with non-linear behaviours, can serve as a sharp boundary beyond which the meadow transitions to bare sand. Determining system behaviour is therefore essential to ensuring resilience is maintained and to prevent stubborn critical ecosystem transitions caused by declines in water quality. Here we combined field and modelling studies to explore the transition from meadow to bare sand in the seagrass Cymodocea nodosa at the limit of its depth distribution in a shallow, light-limited bay. We first describe the relationship between light availability and seagrass density along a depth gradient in an extensive unfragmented meadow (Alfacs bay, NE Spain). We then develop a simple mechanistic model to characterise system behaviour. In the field, we identified sharp decline in shoot density beyond a threshold of ∼1.9 m depth, shifting from a vegetated state to bare sand. The dynamic population model we developed assumes light-dependent growth and an inverse density-dependent mortality due to facilitation between shoots (mortality rate decreases as shoot density increases). The model closely tracked our empirical observations, and both the model and the field data showed signs of bistability. This strongly suggests that the depth limit of C. nodosa is a critical transition driven by photosynthetic light requirements. While the mechanisms still need to be confirmed with experimental evidence, recognizing the non-linear behaviour of C. nodosa meadows is vital not only in improving our understanding of light effects on seagrass dynamics, but also in managing shallow-water meadows. Given the shallow threshold (<2m), light-limited systems may experience significant and recalcitrant meadow retractions with even small changes in sediment and light conditions. Understanding the processes underlying meadow resilience can inform the maintenance and restoration of meadows worldwide.

Spatial patterns in ecological systems: from microbial colonies to landscapes

Self-organized spatial patterns are ubiquitous in ecological systems and allow populations to adopt non-trivial spatial distributions starting from disordered configurations. These patterns form due to diverse nonlinear interactions among organisms and between organisms and their environment, and lead to the emergence of new (eco)system-level properties unique to self-organized systems. Such pattern consequences include higher resilience and resistance to environmental changes, abrupt ecosystem collapse, hysteresis loops, and reversal of competitive exclusion. Here, we review ecological systems exhibiting self-organized patterns. We establish two broad pattern categories depending on whether the self-organizing process is primarily driven by nonlinear density-dependent demographic rates or by nonlinear density-dependent movement. Using this organization, we examine a wide range of observational scales, from microbial colonies to whole ecosystems, and discuss the mechanisms hypothesized to underlie observed patterns and their system-level consequences. For each example, we review both the empirical evidence and the existing theoretical frameworks developed to identify the causes and consequences of patterning. Finally, we trace qualitative similarities across systems and propose possible ways of developing a more quantitative understanding of how self-organization operates across systems and observational scales in ecology.

Letter to the editor regarding the article "Taking advantage of seagrass recovery potential to develop novel and effective meadow rehabilitation methods" by Alagna et al., published in Marine Pollution Bulletin, 149: 2019 (110578)

Alagna et al. (2019) suggest new transplantation methods for Posidonia oceanica (Linnaeus) Delile, inspired by its natural recovery process after disturbance due to dredging operations for gas-pipelines. They observe that P. oceanica vegetative fragments naturally settled only on loose calcareous stones deployed to fill the trenches of the gas-pipeline. No recovery was noted on dead matte, sand and large calcarenitic boulders. Following a new pilot restoration project currently ongoing in the same area, we demonstrate that natural recovery also occurs on dead matte. After examining other alternative transplantation methods for P. oceanica, the Authors suggest using their "habitat enhancement units" method for the restoration of seagrasses, not only on rocky bottom but also on sand and other bare substrate requiring general environmental restoration. Here we express disagreement on certain issues reported in the paper.

Distribution of growth directions in meadows of clonal plants

Clonal growth plants are abundant in both terrestrial and marine ecosystems. Some marine species are particularly important since they provide essential ecosystem services in the shores of all continents except Antarctica. For the appropriate modelization of clonal growth the discretization of the growth direction angle has to be carefully treated to correctly describe the dynamics of the meadow. Specifically, determining the minimum number of growth directions is important to maximize the efficiency of numerical simulations. We show that the presence of neutral modes in the growth direction angle tends to make the distribution of the growth directions uniform in the presence of noise, allowing us to choose the minimal number of discrete angles compatible with the branching. We also show that the formation of spatial patterns induces small differences in the population density within these different growth directions.

A model for the biomass-density dynamics of seagrasses developed and calibrated on global data

BACKGROUND

Seagrasses are foundation species in estuarine and lagoon systems, providing a wide array of services for the ecosystem and the human population. Understanding the dynamics of their stands is essential in order to better assess natural and anthropogenic impacts. It is usually considered that healthy seagrasses aim to maximize their stand biomass (g DW m-2) which may be constrained by resource availability i.e., the local environment sets a carrying capacity. Recently, this paradigm has been tested and reassessed, and it is believed that seagrasses actually maximize their efficiency of space occupation-i.e., aim to reach an interspecific boundary line (IBL)-as quick as possible. This requires that they simultaneously grow in biomass and iterate new shoots to increase density. However, this strategy depresses their biomass potential.

RESULTS

to comply with this new paradigm, we developed a seagrass growth model that updates the carrying capacities for biomass and shoot density from the seagrass IBL at each time step. The use of a joint biomass and density growth rates enabled parameter estimation with twice the sample sizes and made the model less sensitive to episodic error in either of the variables. The use of instantaneous growth rates enabled the model to be calibrated with data sampled at widely different time intervals. We used data from 24 studies of six seagrass species scattered worldwide. The forecasted allometric biomass-density growth trajectories fit these observations well. Maximum growth and decay rates were found consistently for each species. The growth rates varied seasonally, matching previous observations.

CONCLUSIONS

State-of-art models predicting both biomass and shoot density in seagrass have not previously incorporated our observation across many seagrass species that dynamics depend on current state relative to IBL. Our model better simulates the biomass-density dynamics of seagrass stands while shedding light on its intricacies. However, it is only valid for established patches where dynamics involve space-filling, not for colonization of new areas.

A functional-structural model of ephemeral seagrass growth influenced by environment

Background and Aims

Ephemeral seagrasses that respond rapidly to environmental changes are important marine habitats. However, they are under threat due to human activity and are logistically difficult and expensive to study. This study aimed to develop a new functional-structural environmentally dependent model of ephemeral seagrass, able to integrate our understanding of ephemeral seagrass growth dynamics and assess options for potential management interventions, such as seagrass transplantation.

Methods

A functional-structural plant model was developed in which growth and senescence rates are mechanistically linked to environmental variables. The model was parameterized and validated for a population of Halophila stipulacea in the Persian Gulf.

Key Results

There was a good match between empirical and simulated results for the number of apices, net rhizome length or net number of internodes using a 330 d simulation. Simulated data were more variable than empirical data. Simulated structural patterns of seagrass rhizome growth qualitatively matched empirical observations.

Conclusions

This new model successfully simulates the environmentally dependent growth and senescence rates of our case-study ephemeral seagrass species. It produces numerical and visual outputs that help synthesize our understanding of how the influence of environmental variables on plant functional processes affects overall growth patterns. The model can also be used to assess the potential outcomes of management interventions like seagrass transplantation, thus providing a useful management tool. It is freely available and easily adapted for new species and locations, although validation with more species and environments is required.

Fairy circle landscapes under the sea

Short-scale interactions yield large-scale vegetation patterns that, in turn, shape ecosystem function across landscapes. Fairy circles, which are circular patches bare of vegetation within otherwise continuous landscapes, are characteristic features of semiarid grasslands. We report the occurrence of submarine fairy circle seascapes in seagrass meadows and propose a simple model that reproduces the diversity of seascapes observed in these ecosystems as emerging from plant interactions within the meadow. These seascapes include two extreme cases, a continuous meadow and a bare landscape, along with intermediate states that range from the occurrence of persistent but isolated fairy circles, or solitons, to seascapes with multiple fairy circles, banded vegetation, and "leopard skin" patterns consisting of bare seascapes dotted with plant patches. The model predicts that these intermediate seascapes extending across kilometers emerge as a consequence of local demographic imbalances along with facilitative and competitive interactions among the plants with a characteristic spatial scale of 20 to 30 m, consistent with known drivers of seagrass performance. The model, which can be extended to clonal growth plants in other landscapes showing fairy rings, reveals that the different seascapes observed hold diagnostic power as to the proximity of seagrass meadows to extinction points that can be used to identify ecosystems at risks.

Evidences of adaptive traits to rocky substrates undermine paradigm of habitat preference of the Mediterranean seagrass Posidonia oceanica

Posidonia oceanica meadows are acknowledged as one of the most valuable ecosystems of the Mediterranean Sea. P. oceanica has been historically described as a species typically growing on mobile substrates whose development requires precursor communities. Here we document for the first time the extensive presence of sticky hairs covering P. oceanica seedling roots. Adhesive root hairs allow the seedlings to firmly anchor to rocky substrates with anchorage strength values up to 5.23 N, regardless of the presence of algal cover and to colonise bare rock without the need for precursor assemblages to facilitate settlement. Adhesive root hairs are a morphological trait common on plants living on rocks in high-energy habitats, such as the riverweed Podostemaceae and the seagrass Phyllospadix scouleri. The presence of adhesive root hairs in P. oceanica juveniles suggests a preference of this species for hard substrates. Such an adaptation leads to hypothesize a new microsite driven bottleneck in P. oceanica seedling survival linked to substrate features. The mechanism described can favour plant establishment on rocky substrates, in contrast with traditional paradigms. This feature may have strongly influenced P. oceanica pattern of colonisation through sexual propagules in both the past and present.

Contemporary connectivity is sustained by wind- and current-driven seed dispersal among seagrass meadows

BACKGROUND

Seagrasses are clonal marine plants that form important biotic habitats in many tropical and temperate coastal ecosystems. While there is a reasonable understanding of the dynamics of asexual (vegetative) growth in seagrasses, sexual reproduction and the dispersal pathways of the seeds remain poorly studied. Here we address the potential for a predominantly clonal seagrass, P. australis, to disperse over long distances by movement of floating fruit via wind and surface currents within the coastal waters of Perth, Western Australia. We first simulated the dominant atmospheric and ocean forcing conditions that are known to disperse these seagrass seeds using a three-dimensional numerical ocean circulation model. Field observations obtained at 8 sites across the study area were used to validate the model performance over ~2 months in summer when buoyant P. australis fruit are released into the water column. P. australis fruit dispersal trajectories were then quantified throughout the region by incorporating key physical properties of the fruit within the transport model. The time taken for the floating fruit to release their seed (dehiscence) was incorporated into the model based on laboratory measurements, and was used to predict the settlement probability distributions across the model domain.

RESULTS

The results revealed that high rates of local and regional demographic connectivity among P. australis meadows are achieved via contemporary seed dispersal. Dispersal of seeds via floating fruit has the potential to regularly connect meadows at distances of 10s of kilometres (50% of seeds produced) and infrequently for meadows at distances 100 s km (3% of seeds produced).

CONCLUSIONS

The spatial patterns of seed dispersal were heavily influenced by atmospheric and oceanographic conditions, which generally drove a northward pattern of connectivity on a regional scale, but with geographical barriers influencing finer-scale connectivity pathways at some locations. Such levels of seed dispersal infer greater levels of ecological and genetic connectivity and suggest that seagrasses are not just strongly clonal.

Toxic effects of increased sediment nutrient and organic matter loading on the seagrass Zostera noltii

As a result of anthropogenic disturbances and natural stressors, seagrass beds are often patchy and heterogeneous. The effects of high loads of nutrients and organic matter in patch development and expansion in heterogeneous seagrass beds have, however, poorly been studied. We experimentally assessed the in situ effects of sediment quality on seagrass (Zostera noltii) patch dynamics by studying patch (0.35 m diameter) development and expansion for 4 sediment treatments: control, nutrient addition (NPK), organic matter addition (OM) and a combination (NPK+OM). OM addition strongly increased porewater sulfide concentrations whereas NPK increased porewater ammonium, nitrate and phosphate concentrations. As high nitrate concentrations suppressed sulfide production in NPK+OM, this treatment was biogeochemically comparable to NPK. Sulfide and ammonium concentrations differed within treatments, but over a 77 days period, seagrass patch survival and expansion were impaired by all additions compared to the control treatment. Expansion decreased at porewater ammonium concentrations >2,000 μmol L(-1). Mother patch biomass was not affected by high porewater ammonium concentrations as a result of its detoxification by higher seagrass densities. Sulfide concentrations >1,000 μmol L(-1) were toxic to both patch expansion and mother patch. We conclude that patch survival and expansion are constrained at high loads of nutrients or organic matter as a result of porewater ammonium or sulfide toxicity.

Implications of extreme life span in clonal organisms: millenary clones in meadows of the threatened seagrass Posidonia oceanica

The maximum size and age that clonal organisms can reach remains poorly known, although we do know that the largest natural clones can extend over hundreds or thousands of metres and potentially live for centuries. We made a review of findings to date, which reveal that the maximum clone age and size estimates reported in the literature are typically limited by the scale of sampling, and may grossly underestimate the maximum age and size of clonal organisms. A case study presented here shows the occurrence of clones of slow-growing marine angiosperm Posidonia oceanica at spatial scales ranging from metres to hundreds of kilometres, using microsatellites on 1544 sampling units from a total of 40 locations across the Mediterranean Sea. This analysis revealed the presence, with a prevalence of 3.5 to 8.9%, of very large clones spreading over one to several (up to 15) kilometres at the different locations. Using estimates from field studies and models of the clonal growth of P. oceanica, we estimated these large clones to be hundreds to thousands of years old, suggesting the evolution of general purpose genotypes with large phenotypic plasticity in this species. These results, obtained combining genetics, demography and model-based calculations, question present knowledge and understanding of the spreading capacity and life span of plant clones. These findings call for further research on these life history traits associated with clonality, considering their possible ecological and evolutionary implications.

Modelling seagrass growth and development to evaluate transplanting strategies for restoration

BACKGROUND AND AIMS

Seagrasses are important marine plants that are under threat globally. Restoration by transplanting vegetative fragments or seedlings into areas where seagrasses have been lost is possible, but long-term trial data are limited. The goal of this study is to use available short-term data to predict long-term outcomes of transplanting seagrass.

METHODS

A functional-structural plant model of seagrass growth that integrates data collected from short-term trials and experiments is presented. The model was parameterized for the species Posidonia australis, a limited validation of the model against independent data and a sensitivity analysis were conducted and the model was used to conduct a preliminary evaluation of different transplanting strategies.

KEY RESULTS

The limited validation was successful, and reasonable long-term outcomes could be predicted, based only on short-term data.

CONCLUSIONS

This approach for modelling seagrass growth and development enables long-term predictions of the outcomes to be made from different strategies for transplanting seagrass, even when empirical long-term data are difficult or impossible to collect. More validation is required to improve confidence in the model's predictions, and inclusion of more mechanism will extend the model's usefulness. Marine restoration represents a novel application of functional-structural plant modelling.

Standardizing methods to address clonality in population studies

Although clonal species are dominant in many habitats, from unicellular organisms to plants and animals, ecological and particularly evolutionary studies on clonal species have been strongly limited by the difficulty in assessing the number, size and longevity of genetic individuals within a population. The development of molecular markers has allowed progress in this area, and although allozymes remain of limited use due to their typically low level of polymorphism, more polymorphic markers have been discovered during the last decades, supplying powerful tools to overcome the problem of clonality assessment. However, population genetics studies on clonal organisms lack a standardized framework to assess clonality, and to adapt conventional data analyses to account for the potential bias due to the possible replication of the same individuals in the sampling. Moreover, existing studies used a variety of indices to describe clonal diversity and structure such that comparison among studies is difficult at best. We emphasize the need for standardizing studies on clonal organisms, and particularly on clonal plants, in order to clarify the way clonality is taken into account in sampling designs and data analysis, and to allow further comparison of results reported in distinct studies. In order to provide a first step towards a standardized framework to address clonality in population studies, we review, on the basis of a thorough revision of the literature on population structure of clonal plants and of a complementary revision on other clonal organisms, the indices and statistics used so far to estimate genotypic or clonal diversity and to describe clonal structure in plants. We examine their advantages and weaknesses as well as various conceptual issues associated with statistical analyses of population genetics data on clonal organisms. We do so by testing them on results from simulations, as well as on two empirical data sets of microsatellites of the seagrasses Posidonia oceanica and Cymodocea nodosa. Finally, we also propose a selection of new indices and methods to estimate clonal diversity and describe clonal structure in a way that should facilitate comparison between future studies on clonal plants, most of which may be of interest for clonal organisms in general.

Clonal extent, apical dominance and networking features in the phalanx angiosperm Zostera noltii Hornem

Disaggregating seagrass meadows and studying its components separately (clones, ramets, shoots) can provide us insights on meadow dynamics and growth patterns. The clonal growth, dependent upon clonal rules may regulate and impose constraints to plant architecture and, therefore, determine how individual clones evolve into the environment. In order to investigate the relationship between clonal growth rules and clone architecture, the belowground network architecture of single-clones of the seagrass Zostera noltii was studied. Networks were traced in situ after washing out the overlying sediment, and network characteristics were measured using digital analysis: area covered by clone, total rhizome length, type of rhizomatic axes (main, secondary, tertiary, quaternary), number and length of the internodes, branching angles and branching frequencies. This approach revealed that Z. noltii is able to develop into large clones integrating up to 300 internodes, 676 cm of rhizome, 208 shoots and 4,300 cm2 of plant area. Internodal length depended on both, the distance to the apical shoot (time effect) and the axes type (apical dominance effect). However, average branching angle was independent of axis type (average 58.3 ± 0.75), but varied significantly depending on the distance from the apical shoot. This average branching angle allows Z. noltii maximize the rate of centrifugal expansion, maintaining a high density in colonized areas to produce close stands but also minimizing the investment in belowground biomass and ramets overlapping. The clonal architecture of Z. noltii seems to be regulated by the interaction of both, apical dominance strength and clonal integration distance. Moreover, clonal growth rules and growth pattern seem to constrain clonality through (clonal) plant architecture regulations (i.e. branching is restricted in secondary axes, similar average branching angles regardless the axes, the higher the distance to the apex the higher the number of internodes in secondary axes, shorter internodes in secondary and tertiary axes). Future research efforts should focus on how these complex relationships between apical dominance and clonal integration interact to elucidate the temporal (seasonal) and spatial scales of both processes and the outcome at the plant architectural level.

Fisher waves and front roughening in a two-species invasion model with preemptive competition

We study front propagation when an invading species competes with a resident; we assume nearest-neighbor preemptive competition for resources in an individual-based, two-dimensional lattice model. The asymptotic front velocity exhibits an effective power-law dependence on the difference between the two species' clonal propagation rates (key ecological parameters). The mean-field approximation behaves similarly, but the power law's exponent slightly differs from the individual-based model's result. We also study roughening of the front, using the framework of nonequilibrium interface growth. Our analysis indicates that initially flat, linear invading fronts exhibit Kardar-Parisi-Zhang (KPZ) roughening in one transverse dimension. Further, this finding implies, and is also confirmed by simulations, that the temporal correction to the asymptotic front velocity is of O(t(-2/3)).

Within-population spatial genetic structure, neighbourhood size and clonal subrange in the seagrass Cymodocea nodosa

Abstract The extent of clonality within populations strongly influences their spatial genetic structure (SGS), yet this is hardly ever thoroughly analysed. We employed spatial autocorrelation analysis to study effects of sexual and clonal reproduction on dispersal of the dioecious seagrass Cymodocea nodosa. Analyses were performed both at genet level (i.e. excluding clonal repeats) and at ramet level. Clonal structure was characterized by the clonal subrange, a spatial measure of the linear limits where clonality still affects SGS. We show that the clonal subrange is equivalent to the distance where the probability of clonal identity approaches zero. This combined approach was applied to two meadows with different levels of disturbance, Cadiz (stable) and Alfacs (disturbed). Genotypic richness, the proportion of the sample representing distinct genotypes, was moderate (0.38 Cadiz, 0.46 Alfacs) mostly due to dominance of a few clones. Expected heterozygosities were comparable to those found in other clonal plants. SGS analyses at the genet level revealed extremely restricted gene dispersal in Cadiz (Sp = 0.052, a statistic reflecting the decrease of pairwise kinship with distance), the strongest SGS found for seagrass species, comparable only to values for selfing herbaceous land plants. At Cadiz the clonal subrange extended across shorter distances (20-25 m) than in Alfacs (30-35 m). Comparisons of sexual and vegetative components of gene dispersal suggest that, as a dispersal vector within meadows, clonal spread is at least as important as sexual reproduction. The restricted dispersal and SGS pattern in both meadows indicates that the species follows a repeated seedling recruitment strategy.