## New Research In

### Physical Sciences

### Social Sciences

#### Featured Portals

#### Articles by Topic

### Biological Sciences

#### Featured Portals

#### Articles by Topic

- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology

# Ecological processes can synchronize marine population dynamics over continental scales

Edited* by Simon A. Levin, Princeton University, Princeton, NJ, and approved April 2, 2010 (received for review December 17, 2009)

## Abstract

Determining the relative importance of local and regional processes for the distribution of population abundance is a fundamental but contentious issue in ecology. In marine systems, classical theory holds that the influence of demographic processes and dispersal is confined to local populations whereas the environment controls regional patterns of abundance. Here, we use spatial synchrony to compare the distribution of population abundance of the dominant mussel *Mytilus californianus* observed along the West Coast of the United States to that predicted by dynamical models undergoing different dispersal and environmental treatments to infer the relative influence of local and regional processes. We reveal synchronized fluctuations in the abundance of mussel populations across a whole continent despite limited larval dispersal and strong environmental forcing. We show that dispersal among neighboring populations interacts with local demographic processes to generate characteristic patterns of spatial synchrony that can govern the dynamic distribution of mussel abundance over 1,800 km of coastline. Our study emphasizes the importance of dispersal and local dynamics for the distribution of abundance at the continental scale. It further highlights potential limits to the use of “climate envelope” models for predicting the response of large-scale ecosystems to global climate change.

Synchronized fluctuations in abundance among spatially segregated populations are common in nature and can be used to quantify and understand the distribution of abundance in space and time (1). Synchrony can be induced by local intrinsic processes such as dispersal among populations and strong interactions with mobile predators or regional extrinsic processes such as spatially correlated environmental variability (1). Although these processes are well known, identifying their relative contribution to patterns of synchrony remains a challenge (1). Recent work has shown that when the processes that contribute to synchrony can be studied in isolation, be it via natural barriers to dispersal among populations (2, 3) or experimental manipulation (4), synchrony patterns can be ascribed to their underlying cause. However, when intrinsic and extrinsic causes of synchrony co-occur, as is the case in most systems, assigning synchrony patterns to any specific causal process becomes onerous (1). Here, we show that in marine populations experiencing both intrinsic and extrinsic sources of synchrony, the shape of spatial synchrony patterns can be used to infer the cause of synchrony and explain the regional distribution of abundance.

Marine population theory has relied mostly on the environment to explain the regional (>1,000 km) dynamics of populations. This focus is motivated by the lengthy pelagic larval stage commonly found in marine organisms, during which the larvae can be transported over large distances by strong nearshore currents (5). The potential for large-scale transport, along with the difficulties associated with measuring larval dispersal, has prompted many studies to assume either completely closed (no exportation of larvae to other populations) or completely open (no coupling between larval production and recruitment) demography (5, 6, but see refs. 7, 8). This assumption, typically associated with equilibrium dynamics at the local scale, has emphasized the effect of large-scale heterogeneity in nearshore environmental conditions on recruitment (i.e., supply-side theory) to explain the regional dynamics of marine populations (6, 9). However, recent progress on the quantification of larval dispersal distance has motivated the relaxation of demographic openness in several marine species (10–12). In light of these recent developments, we relax the assumptions of demographic openness and local equilibrium dynamics and compare the distribution of population abundance, predicted by dynamical metapopulation models undergoing different dispersal and environmental treatments, to that of the dominant mussel *Mytilus californianus*, observed along the West Coast of the United States, to assess the relative importance of nearshore environmental heterogeneity and dispersal.

## Results and Discussion

### Role of the Environment in Natural Mussel Populations.

We first focus on the role of environmental heterogeneity by quantifying the strength and the consistency of the relationship between nearshore environmental conditions and the abundance (% cover) of *M. californianus*. Although nearshore environmental conditions have a strong effect on patterns of recruitment (13, 14), that effect seems lost on the regional distribution of *M. californianus* cover (15) (Fig. 1*A*). Indeed, the spatial correlation between the environmental conditions and the mean annual *M. californianus* cover is relatively weak and inconsistent through time (Fig. 1*A*), regardless of the temporal lag used (Fig. S1), the temporal scale over which the environment is averaged (Fig. S2) or whether the analysis explicitly accounts for the spatial heterogeneity observed in nearshore conditions (Fig. S3). This weak spatial correlation reflects the mismatch between the persistent spatial gradient in the environment and the spatiotemporal variability in the cover of *M. californianus* (Fig. S3). This mismatch leads to local correlations whose strength and sign vary in space and time (Fig. S3). Once spatially averaged, these highly variable local correlations lead to weak and inconsistent spatial correlations at the regional scale (Fig. 1*A* and Fig. S2). Spatial synchrony, which measures the correlation in the time series of pairs of sites as a function of the lag distance that separates them (*Materials and Methods*), more succinctly reveals the same discrepancy between the persistent spatial gradient in nearshore environmental conditions and the more complex spatiotemporal patterns in the cover of *M. californianus* (Fig. 1*B*). All environmental variables undergo a slow and statistically significant linear decrease in synchrony with increasing lag distance (Fig. 1*B*, see *Materials and Methods* and Table 1 for statistical details). However, the cover of *M. californianus* shows a statistically significant nonlinear pattern, oscillating between synchrony and asynchrony with increasing lag distance (Fig. 1*B* and Table 1). This discrepancy suggests that intrinsic processes (i.e., dispersal and species interactions) rather than local nearshore environmental conditions may control the spatial synchrony patterns exhibited by *M. californianus*.

### Effect of Dispersal and Species Interactions on Natural and Model Mussel Populations.

To elucidate how dispersal and species interactions can generate the complex spatial synchrony patterns observed in natural *M. californianus* populations, we develop metapopulation models that describe disturbance–succession (16, 17) and predator–prey (18–20) dynamics in a network of mussel populations connected by dispersal. Dispersal was implemented as a symmetrical kernel in the successional model (*Materials and Methods*) and as a symmetrical and uniform nearest-neighbor process in the predator–prey model (*SI Text*). We subjected these metapopulation models to environmental variability treatments based on the nearshore environmental conditions observed along the West Coast of the United States by varying mussel fecundity (21) according to a linear spatial gradient (i.e., spatial environmental variability) and a linear spatial gradient with normally distributed white noise (i.e., spatiotemporal environmental variability; *Materials and Methods*). We now perform a factorial experiment on the metapopulation models by using different dispersal and environmental treatments to assess the relative importance of dispersal and environmental variability for generating nonlinear spatial synchrony patterns of abundance that are compatible with those observed in natural mussel populations.

When there is no dispersal (i.e., full local retention of larvae), regional mussel dynamics are strictly controlled by environmental heterogeneity. Under this scenario, the model metapopulations experiencing either spatial or spatiotemporal environmental variability predict a weak and inconsistent spatial correlation between the mussel cover and the 1-year lagged environmental conditions (Fig. 2 *A* and *B*). This weak and inconsistent relationship is compatible with the results from our survey data (Fig. 1*A*) and previously published accounts (15). However, model metapopulations predict that environmental variability induces a rapid decay in the spatial synchrony pattern of the mussel cover, which is incompatible with the nonlinear spatial synchrony pattern observed in natural populations of *M. californianus* (Fig. 2 *C* and *D*). We now introduce dispersal among populations and vary the environment and the scale of dispersal to determine their relative contribution to spatial synchrony.

When dispersal is limited to neighboring populations (herein “limited dispersal”; 8.6% of the spatial domain for results in Figs. 3 and 4), model metapopulations undergoing either spatial or spatiotemporal environmental variability predict a weak and inconsistent relationship between mussel cover and 1-year lagged environmental conditions (Fig. 3 *A* and *C*) that is similar to the one observed in the survey data (Fig. 1*A*). However, regional dispersal (44% of the spatial domain in Figs. 3 and 4) leads to a strong and dynamical relationship between the mussel cover and the environment (Fig. 3 *B* and *D*) that is inconsistent with the weak spatial correlation observed between populations of *M. californianus* and nearshore environmental conditions (Fig. 1*A*). This suggests that limited dispersal between fluctuating populations might be an important driver of abundance in natural populations of *M. californianus*. Indeed, this role of limited dispersal is made more evident through the analysis of spatial synchrony.

In metapopulations experiencing limited dispersal, both the successional model (Fig. 4 *A*, *C*, and *E*) and the predator–prey model (Fig. S4 *A*, *C*, and *E*) accurately predict the nonlinear spatial synchrony pattern observed in natural populations of *M. californianus*, regardless of environmental variability. In contrast, in metapopulations undergoing regional dispersal, both the successional model (Fig. 4 *B*, *D*, and *F*) and the predator–prey model (Fig. S4 *A*, *C*, and *E*) predict that the spatial synchrony pattern displayed by mussels will match the quasilinear decay displayed by the environment, despite the strong but inconsistent spatial correlation between the environment and mussel cover (Fig. 3 *B* and *D*). Hence, in successional and predator–prey metapopulation models, limited dispersal is critical for the emergence of nonlinear spatial synchrony patterns that are compatible with those observed in natural populations of *M. californianus*. Nonlinear spatial synchrony patterns arise because limited dispersal (Fig. S5) couples neighboring populations and thus allows local fluctuations (Fig. S6) to scale up and generate complex, nonstationary spatiotemporal patterns at the regional scale (Figs. S7 and S8) that are robust to environmental forcing (Fig. 4 *B*, *D*, and *F* and Fig. S4 *B*, *D*, and *F*). This cross-scale interaction between local population dynamics and limited dispersal is a general property of metapopulations (Fig. 4 and Fig. S4) that merely requires that local populations undergo sustained fluctuations (Fig. S6) and that the average dispersal distance represents 3–10% of the spatial domain (Fig. S5). Regional dispersal prevents these cross-scale interactions by spatially synchronizing population fluctuations across the entire metapopulation and thus generating regular regional oscillations characterized by stationary and quasilinear spatial synchrony patterns (Fig. 4 *B*, *D*, and *F* and Figs. S4 *B*, *D*, and *F*, S9, and S10). The shape of spatial synchrony can thus be used to determine the relative influence of cross-scale interactions and environmental forcing on the distribution of abundance in space and time. This can be achieved by fitting nonlinear and linear statistical models representing, respectively, cross-scale interactions and environmental forcing to the observed spatial synchrony patterns and then comparing their performance using model selection (*Data Analysis* and Fig. S5).

### Using Spatial Synchrony to Quantify the Scale of Dispersal in Natural Populations.

Our metapopulation models show that the shape of nonlinear spatial synchrony patterns can be used to quantify the scale of dispersal. Indeed, the spatial range, defined as the lag distance at which synchrony first reaches zero (22), is systematically associated with the average dispersal distance in both our successional model (Fig. 4 *A*, *C*, and *E*) and our predator–prey model (Fig. S4 *A*, *C*, and *E*) when dispersal is limited. Applying this result to our survey data, we estimate that the scale of dispersal of *M. californianus* is ≈100 km (i.e., 6% of the domain). This estimate falls within the 95% confidence interval of dispersal distances documented for other bivalve species (72–220 km, see ref. 11) and is very similar to the scale of dispersal (97–115 km) of more closely related *Mytilus* species in other systems (11). Our estimate is based on the number of recruits that survive to the adult stage and is smaller than the 250-km estimate derived empirically by measuring the density of settlers in the same system (23). This is because spatial synchrony estimates integrate postsettlement processes that can limit the effective scale of dispersal.

### Spatial Synchrony Reveals the Relative Influence of Intrinsic and Extrinsic Processes on the Distribution of Population Abundance.

Theory has shown that local dynamics and limited dispersal can lead to complex spatial (24, 25) or spatiotemporal (26, 27) patterns at the regional scale that can promote coexistence (28), stability (27, 29), persistence (26, 29, 30), and functioning (24, 25) in spatial ecological systems. Similar complex spatiotemporal patterns have also been used to describe insect (31) and epidemic (32) outbreaks at regional and continental scales. However, linking local dynamics and limited dispersal to regional patterns has typically required extensive time series and characteristic spatial signatures such as Turing structures (24, 25), traveling waves (26, 31), or power laws (33). Here, we extend these theories by showing that in systems lacking these characteristic spatial signatures, the shape of spatial synchrony patterns can be used to infer the joint effect of local dynamics and dispersal on the regional distribution of abundance. Hence, our framework moves away from the use of correlations between abundance and the environment to infer the causal effect of abiotic processes on biological patterns. Instead, by partitioning long nonstationary time series into smaller quasistationary time series and applying spatial synchrony analysis, we show that the spatial and temporal properties of spatial synchrony patterns can be used to determine the relative influence of intrinsic and extrinsic factors on the regional distribution of abundance in natural systems: limited dispersal interacts with local intrinsic fluctuations to generate nonlinear and nonstationary spatial synchrony patterns that are robust to environmental forcing, whereas regional dispersal facilitates environmental forcing and leads to stationary and linear spatial synchrony patterns.

Our results have important implications for identifying the processes that control the distribution of species abundance in natural systems and for predicting their response to global climate change. A popular approach for understanding and predicting species abundance distributions is to build species’ “climate envelopes” by either mapping the current species distribution to climate variables via correlation techniques or by determining the physiological tolerances of individual species (see review in ref. 34). This climate envelope can then be used to predict the future distribution of species under various climate change scenarios (34). However, climate envelope models have been criticized because they do not integrate the effects of species interactions and dispersal on the distribution of population abundance (35–38). A recent synthesis proposes integrating the effects of dispersal and species interactions into climate envelopes by adopting a hierarchical modeling framework (34). According to this hierarchical framework, species interactions and dispersal control the distribution of species abundance at smaller spatial scales (<200 km) whereas climate dominates at larger spatial scales (>200 km) (34). Here, we have shown that despite strong regional environmental forcing, local dynamics interact with limited dispersal to control the distribution of population abundance at scales that are much larger (>1,000 km) than that of dispersal (∼100 km). Hence, our work suggests that processes occurring at small spatial scales can interact synergistically to control the distribution of population abundance at large spatial scales. Such cross-scale interactions demonstrate the limitations of adopting climate envelope models based on hierarchical frameworks to understand the distribution of species abundance and predict the effects of global climate change.

Overall, by applying spatial synchrony analysis to a large data set of mussel populations along the West Coast of the United States, our study provides unique evidence that limited connectivity among local populations affects the dynamic distribution of abundance over >1,000 km. Our work supports the suggested shift toward a more dynamical approach to regional conservation, one that emphasizes patterns and processes across scales instead of those limited to the scale of the environment or dispersal.

## Materials and Methods

### Data Collection.

Abundance (percent cover) of *M. californianus* was quantified annually from 1999 to 2004 at 48 sites located along the West Coast of the United States and stretching from southern California to northern Washington (39) (32.7 °N to 48.4 °N). For each site, the cover of *M. californianus* was surveyed in 10 randomly placed 0.25 m^{2} quadrats for each of three 50-m transects located within the midintertidal zone. Mean annual sea surface temperature (SST, in °C), chlorophyll-a concentration (chl-a, in mg/m^{3}), and upwelling index (in m^{3}/s/100 m of coastline) data from 1997 to 2003 occurring within a 0.2 degree radius (1 degree radius for upwelling) of each of the 48 sites were obtained, respectively, from the advanced very high resolution radiometer (NOAA), the sea-viewing wide field-of-view sensor (NASA), and sea level pressure maps (Pacific Fisheries Environmental Laboratory). These environmental data series were validated by comparing them to in situ buoy measurements (see *Data Validation* in *SI Text*).

### Data Analysis.

Before conducting spatial synchrony analysis, all variables were detrended by subtracting the global mean time series from each site's time series to remove any bias caused by common large-scale trends (40). The distances between all pairs of sites were then computed and used to group the detrended time series data into equally spaced distance bins (66 km wide). The coefficient of synchrony for each bin was calculated by computing the correlation coefficient between the time series of all pairs of sites within the bin. The extent of the spatial synchrony analysis was restricted to half of the spatial domain to avoid large discrepancies in the number of pairs of sites within each bin (22). Statistical significance was determined by using a one-tailed test (α = 0.05) on 10,000 Monte Carlo randomizations (22). Specifically, for each bin, the *P* value was calculated by shuffling the data pairs within the bin 10,000 times, computing the coefficient of synchrony for each randomization, and then calculating the proportion of randomizations with a coefficient of synchrony greater than or equal to that obtained with the original data. The same one-tailed randomization technique was used to assess the statistical significance of the correlation between the annual cover of *M. californianus* at year *i* and each environmental variable at year *i-1* (1-year lag). We used a 1-year lag because it corresponds to the temporal scale at which the correlation between the annual cover of *M. californianus* and each nearshore environmental variable is statistically significant (Fig. S1).

We used model selection to detect linear and nonlinear spatial synchrony patterns in the environmental (SST, chl-a, and upwelling index) and *M*. *californianus* data sets. For each data set, we fit a linear statistical model and a nonlinear statistical model to the spatial synchrony pattern:where is the lag distance vector, represent fitted coefficients and is the normalized lag distance vector scaled to the domain . This scaling of the lag distance vector allows the cosine function in the nonlinear model to fit modal patterns of spatial synchrony over the spatial domain. For each statistical model, we calculated the Akaike information criterion corrected for small samples (41):where represents the maximum likelihood, represents the number of parameters in each statistical model, and represents the number of samples. The statistical model with the smallest value was selected for each data set (41).

### Successional Model.

The successional model describes local disturbance and recovery dynamics in a network of mussel populations that are connected by dispersal (17). Within populations, the successional dynamics observed in natural intertidal systems (16, 42) are represented as a mean-field implementation of a spatial process affecting the proportional abundance of (*i*) the dominant mussel , (*ii*) the wave disturbance , and (*iii*) the empty substrate (43). A maximum fraction of the proportional abundance of the dominant mussel species can be displaced by wave disturbances . A proportion of disturbances displaces mussels through a density-dependent contact process with aggregation (Moore neighborhood, ), whereas a proportion of disturbances is density independent. This disturbance dynamic is based on the assumption that wave disturbances destroy the byssal thread attachments of mussels around the edges of disturbed areas, thus making them temporarily more susceptible to further disturbance (43, 44). Hence, newly disturbed areas allow the local propagation of wave disturbances to adjacent mussel beds. Once the disturbance has propagated away from the newly disturbed area, the area transitions from the “wave disturbed” state to the “empty substrate” state. Similarly to disturbance, a maximum fraction of the empty substrate can be colonized by mussels. A proportion of colonization occurs through a density-independent process, whereas the remaining colonization is density dependent . Mussel colonization also depends on the production and recruitment of larvae. Within populations *x*, larval production is a function of local mussel proportional abundance and fecundity . The recruitment rate is described by a Poisson process (45) , where integrates (*i*) the total number of larvae produced and retained in populations at time and (*ii*) the total number of larvae produced in other populations *y* and dispersed to populations at time . The dynamics of the model are represented by the following integro-difference equation system for populations in a metapopulation consisting of populations:with:where is the symmetrical mussel dispersal kernel resulting from larval transport at a constant speed and with a time-dependent settlement rate (double Weibull distribution) (46), represents the proportion of larvae being dispersed, represents the scale of dispersal, and **L** represents a zero-based vector of population locations. We vary dispersal by manipulating , the proportion of larvae being dispersed ( means that all larvae are retained locally, whereas means that all larvae are dispersed), and , the scale of dispersal ( corresponds to regional dispersal and corresponds to limited dispersal). We assumed periodic boundary conditions for all simulations and set the dispersal rate to (i.e., no local retention) unless otherwise specified (i.e., Fig. 2 where ). All successional model simulations were performed for 256 populations and the results were analyzed over 2,000 posttransient time steps. Because our goal was to test the importance of local ecological processes and dispersal, all parameter values detailed above were selected to be representative of the broad parameter space characterized by spatiotemporal heterogeneity (17, 47). Here, we further assess the model's sensitivity to dispersal distance and to fecundity *f* as a means to determine the role of environmental variability in marine metapopulations.

Environmental variability can have a significant impact on the productivity of intertidal populations (21, 48). Here, we implement this effect by varying mussel fecundity spatially and spatiotemporally. Specifically, the spatial environmental treatment consists of varying the mussel fecundity *f* linearly from 3 to 7.5 over the entire spatial range . We generate the spatiotemporal environmental treatment by adding normally distributed white noise with zero mean and variance to the previously described spatial variation in fecundity . The spatial and spatiotemporal environmental treatments thus preserve the same mean fecundity as the successional model undergoing no environmental variability. These spatial and spatiotemporal treatments were chosen to roughly mimic the spatial and spatiotemporal properties of the environment along the West Coast of the United States. Applying the same spatial and spatiotemporal treatments to the mussel growth rate yields qualitatively similar results.

### Model Analysis.

We applied the spatial synchrony and model selection methods (*Data Analysis*) used on the survey data to the posttransient time series of the successional and the predator–prey (*SI Text*) models to assess their ability to generate nonlinear spatial synchrony patterns that are compatible with those observed in natural populations of *M. californianus*. Specifically, the model time series (2,000 posttransient time steps) was split into 10-time step windows and spatial synchrony analysis was conducted over each window. We chose 10-time step windows to approximate the temporal extent of our intertidal survey data. However, our results are robust to window size (Fig. S5 *B*, *D*, and *F*). For each time window, we fit the same linear and nonlinear statistical models described in *Data Analysis* to the spatial synchrony patterns generated by the metapopulation models. For each time window, the nonlinear statistical model was selected if (41). The model spatial synchrony patterns presented in all figures were computed by averaging the spatial synchrony patterns from all 10-time step windows in which the nonlinear statistical model was selected. When no compatible synchrony pattern exists for the entire model time series (i.e., for the regional dispersal treatment), the spatial synchrony patterns from *n* randomly selected time windows are computed and averaged (where *n* corresponds to the number of compatible synchrony patterns for the limited dispersal treatment). This allows for unbiased comparisons across dispersal treatments.

## Acknowledgments

We thank G. Fussmann, B. Cazelles, and J. van de Koppel for providing insightful comments on a previous version of the manuscript. T.C.G. was supported by a McGill Majors fellowship; F.G. was supported by a grant from the James S. McDonnell Foundation; B.A.M was supported by grants from the David and Lucile Packard Foundation, the Gordon and Betty Moore Foundation, the Wayne and Gladys Valley Foundation, the Andrew W. Mellon Foundation, and the National Science Foundation. This is publication number 360 from the Partnership for Interdisciplinary Studies of Coastal Oceans (PISCO), a large-scale, long-term consortium primarily funded by the David and Lucile Packard Foundation and the Gordon and Betty Moore Foundation.

## Footnotes

^{1}To whom correspondence should be addressed. E-mail: tarik.gouhier{at}gmail.com.Author contributions: T.C.G. and F.G. designed research; T.C.G. performed research; T.C.G. analyzed data; and T.C.G., F.G., and B.A.M. wrote the paper.

The authors declare no conflict of interest.

↵*This Direct Submission article had a prearranged editor.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.0914588107/-/DCSupplemental.

## References

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Connolly SR,
- Roughgarden J

- ↵
- Hastings A,
- Higgins K

- ↵
- ↵
- Roughgarden J,
- Gaines SD,
- Pacala SW

- ↵
- Swearer SE,
- et al.

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Guichard F

- ↵
- ↵
- ↵
- ↵
- Leslie HM,
- Breck EN,
- Chan F,
- Lubchenco J,
- Menge BA

- ↵
- Fortin M,
- Dale M

- ↵
- ↵
- ↵
- van de Koppel J,
- et al.

- ↵
- ↵
- Dieckmann U,
- Law R,
- Metz JAJ

- Jansen V,
- de Roos A

- ↵
- ↵
- ↵
- Hassell MP,
- Comins HN,
- May RM

- ↵
- Bjornstad ON,
- Peltonen M,
- Liebhold AM,
- Baltensweiler W

- ↵
- Viboud C,
- et al.

- ↵
- Solé R,
- Bascompte J

- ↵
- ↵
- Suttle KB,
- Thomsen MA,
- Power ME

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Burnham K,
- Anderson D

- ↵
- ↵
- ↵
- ↵
- ↵
- Neubert MG,
- Kot M,
- Lewis MA

- ↵
- ↵
- Menge BA,
- Daley BA,
- Wheeler PA,
- Strub PT

## Citation Manager Formats

### More Articles of This Classification

### Biological Sciences

### Related Content

- No related articles found.

### Cited by...

- No citing articles found.