# Network Analysis Shows Asymmetrical Flows within a Bird Metapopulation

Emilio R. Rojas, Cédric Sueur, Pierre-Yves Henry, Blandine Doligez, Gérard

Wey, Olivier Dehorter, Sylvie Massemin

► To cite this version:

Emilio R. Rojas, Cédric Sueur, Pierre-Yves Henry, Blandine Doligez, Gérard Wey, et al.. Network Analysis Shows Asymmetrical Flows within a Bird Metapopulation. PLoS ONE, 2016, 11 (11), pp.e0166701. 10.1371/journal.pone.0166701 . hal-01429872

HAL Id: hal-01429872

https://hal.sorbonne-universite.fr/hal-01429872

Submitted on 9 Jan 2017

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Network Analysis Shows Asymmetrical Flows within a Bird Metapopulation

Emilio R. Rojas¹, Cédric Sueur¹, Pierre-Yves Henry²,³, Blandine Doligez⁴, Gérard Wey⁵, Olivier Dehorter³, Sylvie Massemin¹*, Groupe Cigognes France⁵¶

1 IPHC, UNISTRA, CNRS, Strasbourg, France, 2 Mécanismes adaptatifs et Evolution (MECADEV UMR7179), Sorbonne Universités, MNHN, CNRS, Muséum National d'Histoire Naturelle, 1 avenue du Petit Château, Brunoy, France, 3 Centre de Recherches sur la Biologie des Populations d'Oiseaux (CRBPO), Centre d'Ecologie et des Sciences de la Conservation (CESCO UMR7204), Sorbonne Universités, MNHN, CNRS, UPMC, Muséum National d'Histoire Naturelle, Paris, France, 4 UMR CNRS 5558-LBBE, "Biométrie et Biologie Évolutive", Villeurbanne cedex, France, 5 Groupe Cigognes France & APRECIAL, Colmar, France

¶ Membership of the Group Cigognes France is provided in the Acknowledgments.

* sylvie.massemin@iphc.cnrs.fr

## Abstract

How the spatial expansion of a species changes at a human time scale is a process difficult to determine. We studied the dispersal pattern of the French white stork population, using a 21-year ringing/resighting dataset. We used the graph-theory to investigate the strength of links between 5 populations (North-East, North-West, Centre, West, and South) and to determine factors important for the birds' movements. Two clusters of populations were identified within the metapopulation, with most frequent movements of individuals between North-Eastern and Centre populations, and between North-Western and Western populations. Exchanges of individuals between populations were asymmetrical, where North-Eastern and North-Western populations provided more emigrants than they received immigrants. Neither the geographical distance between populations, nor the difference in densities influenced the number of individuals exchanging between populations. The graph-theory approach provides a dynamic view of individual movements within a metapopulation and might be useful for future population studies in the context of conservation.

## Introduction

Theory of metapopulation focuses on the dynamics of individual exchanges between populations [1,2]. It is based on the accumulation of field records concerning the dispersal behaviour of individuals and/or observed gene flow between populations [3]. The majority of past studies focused on the impact of habitat fragmentation, since animal populations were patchily distributed, forming metapopulations. Another important aspect that has been studied concerns the risk of species extinction in dependence of population size and the degree of connectivity between populations. The role of individual dispersal between populations is crucial for the persistence of a metapopulation, in particular for threatened species (e.g. the Cougar, *Puma concolor* or the Florida scrub-jay, *Aphelocoma coeruleescens* [4,5]). Ecological factors also need

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Rojas ER, Sueur C, Henry P-Y, Doligez B, Wey G, Dehorter O, et al. (2016) Network Analysis Shows Asymmetrical Flows within a Bird Metapopulation. PLoS ONE 11(11): e0166701. doi:10.1371/journal.pone.0166701

**Editor:** Antoni Margalida, University of Lleida, SPAIN

**Received:** June 20, 2016

**Accepted:** November 2, 2016

**Published:** November 28, 2016

**Copyright:** © 2016 Rojas et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** Data are from a dataset coordinated by the CRBPO. We attached a file with the dataset.

**Funding:** ER was supported by the Ministère français de l'Enseignement et de la Recherche. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. The other authors received no specific funding for this work.

**Competing Interests:** The authors have declared that no competing interests exist.
to be studied to understand the mechanisms that, for instance lead to a wider geographical dis-
tribution of a metapopulation [6]. While biogeography studied the expansion of species at the
level of islands or continents and over periods of several thousand years, determining the role
of individual dispersal behaviour may help to define the dynamics of metapopulations on a
finer geographical and temporal scale. Harrison [7] proposed that individual qualities deter-
mined how individuals disperse. For instance, young individuals have a propensity to leave the
natal site, to travel long distances over unfavourable habitat, have a higher ability to locate new
habitats and, therefore, a higher capacity to establish new populations from a small numbers of
founders. How the spatial expansion of a species changes at the level of a human time scale (i.e.
less than 40 years) is difficult to determine. To study the intensity of individual exchanges
between populations and the factors influencing dispersal patterns, we used the white stork
(*Ciconia ciconia*) as our model species. Until the 1970's, the distribution of white storks was
restricted to the North-Eastern region of France (i.e. Alsace) but has expended in recent
decades [8]. In 1961, the Alsatian population consisted of 145 pairs but declined to only nine
pairs in 1974 [9]. Such a fast, drastic decline was also observed in other white stork populations
of Western Europe [10–12]. After the implementation of several conservation programs in
Switzerland and the Alsatian region, these populations increased again [8,13,14]. In France,
the North-Eastern population size increased to an almost five fold level in 2015, when com-
pared with the situation in 1961 (Groupe Cigognes France). Other populations started to
establish themselves in different regions of France during the 1980's [15], for example in Cha-
rente-Maritime [16,17] or in the North-Western regions of France [8]. Moreover, the coloniz-
ing white storks have the capacity to install in new habitats. Dispersal in juveniles is high (natal
dispersal: 70% of individuals) and is also observed in adults (breeding dispersal: 6%) [16–18].
Hence, juveniles make up good candidates to disperse and found new populations. The build
up of the French white stork metapopulation, after its near extinction in the 1970's, has not
been analysed yet. Information is only available for the Charente-Maritime population and
only for the period between 1991 and 1996 [16]. We studied the dispersal patterns of white
storks on the national scale (i.e. France) and investigated the importance of the North-Eastern
population as a source population for the establishment of new populations within other
regions of France. For this purpose, we used the graph-theory to study the strength and
dynamics of links between populations. Graph-theory is a mathematical approach dealing
with the connectivity inside a network [19]. In biology, it has been used in particularly for the
analysis of social networks [20–23]. Its use widened to study landscape networks [19,24–26]
and it has recently emerged as a powerful tool to study the exchanges of individuals between
populations [27,28].

In our analysis, we used a ringing/resighting dataset concerning mainland France (~544,000 km²), covering a 21-year period (from 1988 to 2008). Such a large spatial and temporal coverage is essential to investigate the movement patterns within a metapopulation for far-ranging species, such as white storks. The main aim of our study was to investigate the potential asymmetry of individual exchanges between populations to diagnose ‘source-sink’ dynamics during the build-up of a metapopulation. We studied several factors that can affect the dispersal patterns of white storks: 1) attractive or repulsive effects of conspecifics, which have been reported in the literature [29–31]. We tested the effect of the relative density of con-specifics to investigate the trade-off between competition and gregarious behaviour [30,32]. 2) The effect of inter-population distances was evaluated, with the prediction that the majority of individuals would disperse between populations that are closer to each other. 3) The presence of population sub-groups with preferential exchanges of individuals, called clusters, was also investigated in order to assess dispersal patterns inside the metapopulation.
# Materials and Methods

## Study species and data collection

The white stork is a large bird that lives in close association with humans [33] and is present in numerous cultivated habitats in Central and Western Europe [34]. After a critical drop of the Western European stork populations in the 1970's [11,35], conservation programs succeeded to increase bird numbers to their former levels or even beyond. In France, the current white stork population size is higher than before the drop [16,35,36].

White storks were ringed as nestlings (~1.5-month old) using uniquely marked individual rings that could be read up to a distance of 50-70 m using spotting scopes. Ringing and resighting data (i.e. individual identity, date and location) were collected from 1988 to 2008 by volunteer ornithologists and bird ringers that monitor stork populations in France (coordinated by CRBPO until 2008 [37]). Ringing effort (i.e. the proportion of hatched nestlings in a given year that were ringed) decreased relative to the observed population growth. While almost all nestlings were ringed at the beginning of the study period, this was not the case later on, because of the large number of nestlings. However, we are confident in the robustness of our conclusions relative to this drop in marking probability because while the ringing effort varied along the study period, the resighting effort remained similar for each population (i.e. almost all breeders were followed).

We defined the breeding season as the period from the latest laying date recorded for the species [38] to the earliest fledging date [12,39,40], extending from mid-April to mid-July. Immature, non-breeding individuals are known to have erratic movements during the first years of their life [41-43] and were therefore excluded from the analysis. We only considered individuals that were resighted as adults (≥ 4 years old) and for which the probability to be breeders at the time of resighting was high [12,16]. We obtained a total of 1720 dispersal events (for 1087 individuals), containing 239 dispersal events between populations and 1481 dispersal events within populations.

Dispersers were defined as breeding individuals observed at a distance greater than 5km from where they had been observed during the previous breeding season (breeding dispersal) or from their nest of birth (natal dispersal; [12,16]). This 5-km threshold distance was set according to the foraging distance of breeders during the breeding season around their nest [44,45]. Dispersal distance and direction were determined from geographical coordinates of ringing and/or resighting sites. The distance between two populations was defined as the Euclidian distance between the centroids of populations. The mean distance between the centroids of populations in our study was 351.8 ± 30.3 km, a value higher than the mean dispersing distance typically found in this species (18 ±41.7 km in Germany [46], 94 ±132.2 km in Poland [43], and 67.3 ±136.2 km in France, [18]).

We defined five geographical populations, where individuals originated from (i.e. where individuals were born and ringed). The five main populations (West, North-East, North-West, South, and Centre) contained 83% of all ringed individuals of our dataset and 86% of all dispersal events.

In our analysis, we defined 5 temporal classes, with one class spanning the period of 9 years (1988-1996; due to the low number of ringed birds) and four classes spanning a period of 3 years each (1997-1999, 2000-2002, 2003-2005, and 2006-2008). We also considered the density differences observed between populations to investigate if the number of individuals within a population influenced the movements from one population to another. As we used correlation matrices in our analysis, which required one value per pair (dyad) of populations, we used density differences as an index for the effects of conspecifics density. Density was estimated based on the counts of breeding pairs made by the French national network of storks
ornithologists (“Groupe Cigognes France”). In the following, the number of individuals was divided by the surface area of the territory occupied by the population (i.e. number of individuals per 100km²).

## Population network indices

We calculated the following parameters for the entire observation period (21 years) and for each temporal class. Our metapopulation consisted of 5 “nodes” inside the network, the five populations. Links between these nodes represent the numbers of between-population dispersal events. Links can be oriented so as to distinguish between symmetrical and asymmetrical exchanges (for instance, source-sink dispersal) between populations. We calculated the strength of nodes as the sum of dispersal events between a given population and all others [47].

We performed a cluster analysis that indicated if and how populations were connected to each other. To this end, the 1481 dispersal events within their population were excluded from these following analyses. We extracted the maximum modularity (Max. Mod.), which is the fraction of internal connections in each cluster minus the expected fraction, if connections were distributed randomly [21], and also determined the cophenetic correlation coefficients CCC [47]. Modularity indicates to what extent a group is clustered, with a value that ranges between 0 (indicating a perfectly random distribution of exchanges) and 1 (indicating a strong hierarchical clustering [21]). A modularity greater than 0.3 is usually considered to indicate useful divisions of the data [48–50]. The CCC indicates how well the dendrogram (the graphic representation of data clustering) matches the matrix of dispersal events between populations (nodes). A CCC value greater than 0.8 indicates a good match [47,50]. A cluster out/in ratio was also calculated for populations included in clusters. For a given population, the number of interactions outside its cluster was divided by the number of interactions inside its own cluster. This ratio indicated the population role in the cluster. A low out/in ratio indicated populations that interacted mostly within their cluster, while a high out/in ratio suggested external exchanges with other clusters/outside populations.

## Statistical analyses

We tested for correlations between 1) the number of individuals dispersing between two populations and the density differences between these populations; 2) the number of individuals exchanged in relation to the distance between their centroids. We obtained a matrix for each of these parameters and correlations were conducted using a Mantel Z-test [51] with 10,000 permutations.

We further investigated the attractiveness (more individuals arrive in a population than leave) or repulsiveness of a population for dispersing events. We designed three matrices: 1) a matrix containing the number of dispersing individuals within their own populations for each temporal class; 2) a matrix containing the number of individuals arriving in each population for each temporal class; and 3) a matrix with the number of individuals leaving the populations for each temporal class. Thus, we were able to perform matrix correlations using Mantel Z-tests with 10,000 permutations, to test attractiveness or repulsiveness of each population.

The asymmetry between the number of emigrating and immigrating individuals for a given population was established by means of permutation-based chi-squared tests for asymmetry (with 10,000 permutations). Network measurements and analyses were performed using SOC-PROG v2.5 [47] and R (v3.1.1; R Development Core Team) softwares. Gephi [47] was used for drawing figures. Means are presented with SE.
# Results

## Dispersal at metapopulation level

The geographical distance between populations did not influence the number of individuals that they exchanged (Mantel Z-test, Z = -0.245; p = 0.851). The cluster analysis identified two distinct clusters in the global metapopulation network: cluster A, containing the North-East/ Centre dyad and cluster B, containing the North-West/West dyad (Total: Max. Mod. = 0.507, CCC = 0.985; Fig 1, Table 1).

A detailed examination of the temporal dynamics revealed that clusters A and B only appeared during the last years considered, from 2004 onward, with a clear increase in modularity (Fig 1, Table 1). Hence, both clusters stabilised only recently. The dispersal events occurring inside (In) and outside (Out) of the clusters (out/in ratio) varied between populations (permutation-based chi-square test of asymmetry: $\chi^2 = 13.94, p < 0.01$; Table 2), regardless of the temporal class considered. Individuals from the Centre population in Cluster A and from the Western population in Cluster B dispersed in equal proportions in and out of the clusters. On the other hand, the North-Eastern population in Cluster A and the North-Western population in Cluster B contained mainly individuals that dispersed within their respective clusters (i.e. within-cluster dispersal). In other words, North-Eastern and North-Western populations mainly provided individuals to the other population of their cluster (respectively Centre and West populations), whereas these last populations contributed with relatively more between-cluster exchanges of individuals (i.e. they equally provided individuals to their cluster as to the other cluster).

Dispersal links between populations (i.e. edges of the network) were asymmetrical, particularly between the North-Eastern and Centre populations (permutation-based chi-square test of asymmetry: $\chi^2 = 46.43, p < 0.01$), the North-Eastern and Southern populations ($\chi^2 = 11.84, p < 0.01$), and, finally, between the North-Western and Western populations ($\chi^2 = 14.22, p < 0.01$; Fig 1). This asymmetry appeared to be time-dependant and no evidence for asymmetry was found before 2003. Between 2003 and 2005, asymmetry existed only between the North-Eastern and Centre populations ($\chi^2 = 25.33, p < 0.01$) and the North-Western and Western populations ($\chi^2 = 6.37, p = 0.01$; Table 1; the former providing more dispersers than the latter). In the most recent years, between 2006 and 2008, the pattern of asymmetry remained between the North-Eastern and Centre populations ($\chi^2 = 16.89, p < 0.01$), and between the North-Western and Western populations ($\chi^2 = 17.29, p < 0.01$), but also became significant between the North-Eastern and the Southern populations ($\chi^2 = 5.00, p = 0.03$; Fig 1). This suggests that an asymmetry of the edges in our metapopulation appeared only very recently. In addition, we calculated nodal indices to test the strength of links between populations or nodes for each temporal class. This showed that the strength of nodes increased with time, and the strongest nodes were observed during the most recent years for the North-Eastern and Centre populations (permutation-based F-test, 10,000 permutations: F = 32.2, p = 0.030, $r^2 = 0.829$; Fig 1).

## Dispersal at population level

The North-Eastern population was always larger than the Centre population, regardless of the period considered (Fig 1), especially when the two clusters were significantly different (mean ±SE number of individuals/year: 283±8 vs. 25±10 for the North-Eastern and Centre population, respectively in 2000–2002; 390±21 vs. 78±5 in 2003–2005, and 444±27 vs. 115±18 in 2006–2008). The North-Western population, on the other hand, was always smaller than the
1988 - 1996

<figure><img src="image_4.png" /></figure>

1997 - 1999

<figure><img src="image_5.png" /></figure>

2000 - 2002

<figure><img src="image_6.png" /></figure>

2003 - 2005

<figure><img src="image_7.png" /></figure>

2006 - 2008

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

50 Ind.

<figure><img src="image_11.png" /></figure>

25 Ind.

<figure><img src="image_12.png" /></figure>

5 Ind.

<figure><img src="image_13.png" /></figure>
**Fig 1. Visual representation of the metapopulation structure for white storks in France during the five periods considered.**

The curves thickness is proportional to the number of individuals dispersing between the populations. Each node represents one population and is positioned in its center. The size of the nodes is proportional to the number of individuals. Codes for the populations are: W for West, S for South, C for Centre, NE for North-East and NW for North-West.

doi:10.1371/journal.pone.0166701.g001

Western population (48±6 vs 218±63 for the North-Western and Western population, respectively in 2000–2002, 89±8 vs. 397±63 in 2003–2005, and 120±11 vs. 553±57 in 2006–2008).

The difference in densities did not influence the number of individuals exchanging between populations (Mantel Z-test, Z = 0.188; p = 0.250). However, the number of dispersing events within a population can influence its attractiveness/repulsiveness. We found no evidence for an attractiveness effect (Mantel Z-test, Z = 0.197; p = 0.208), while we did find significant evidence for a repulsiveness effect (Mantel Z-test, Z = 0.871; p = 0.008). This indicates that a high level of internal dispersal activity leads to more individuals leaving the population, thus, increasing the contribution of the population to metapopulation dynamics.

Discussion

We conducted a network analysis concerning the strength of links between 5 populations within a white stork metapopulation that increased in number of individuals between 1988 and 2008. We found that the metapopulation is divided into two clusters, opposing North-Eastern/Centre populations to North-Western/Western populations. Clustering was detected starting in 2000, whereas the establishment of a metapopulation was already observed at the beginning of the study period (i.e. 1988). Clustering seems to be sensible to an increasing number of dispersing events, associated with an increase of the different white stork populations in France [16,35,36]. From 1999 onward, the number of dispersal events between populations further increased. The asymmetry analysis conducted in our study suggests a source-sink relationship between populations [52]. Thus, emigration from the North-Eastern and North-Western populations was greater than immigration, so that these populations, took on the role of “source” populations, while the other populations, Western, Centre, and Southern, were “receiver” populations, this immigration likely contributing to their maintenance or growth. However, while the North-Eastern population is larger than the Centre population, the North-Western population is smaller than the Western population. Consequently, there is no support for a density effect that has been observed in Polish and German white stork populations [16,17]. In Germany, the study of Itonaga and colleagues showed that the proportion of dispersers in white storks increased in parallel with the population size [46]. These findings support indirectly the influence of density on the dispersal rate. There are several explanations for this seemingly contradictory pattern of density-dependence in white storks. First, the U-shaped relationship between density and dispersal distance might be responsible for the observed pattern [53]. In colonial species, the dispersal probability/distance is high in populations with a moderate to low density (with birds moving to higher density populations, cf.

<table><caption>**Table 1. Temporal evolution of the asymmetry of links, cophenetic correlation coefficients (CCC), maximum modularities (Max. Mod.), and consecutive cluster determinations for the metapopulation.** Statistically significant results are presented in bold. For the asymmetry tests, “»” indicate the direction of the asymmetry. Codes for the populations are: W for West, S for South, C for Centre, NE for North-East and NW for North-West.</caption><thead><tr><th>Periods</th><th>1988–1996</th><th>1997–1999</th><th>2000–2002</th><th>2003–2005</th><th>2006–2008</th><th>Total</th></tr></thead><tbody><tr><td>CCC</td><td>0.948</td><td>0.954</td><td><strong>0.863</strong></td><td><strong>0.997</strong></td><td><strong>0.976</strong></td><td><strong>0.985</strong></td></tr><tr><td>Max. Mod.</td><td>0.272</td><td>0.285</td><td><strong>0.503</strong></td><td><strong>0.564</strong></td><td><strong>0.520</strong></td><td><strong>0.507</strong></td></tr><tr><td>Clusters</td><td>None</td><td>None</td><td>NE/C NW/W</td><td>NE/C NW/W</td><td>NE/C NW/W</td><td>NE/C NW/W</td></tr><tr><td>Asymmetry</td><td>None</td><td>None</td><td>None</td><td>NE»C NW»W</td><td>NE»C NE»SNW»W</td><td>NE»C NE»SNW»W</td></tr></tbody></table>

doi:10.1371/journal.pone.0166701.t001
<table><caption>**Table 2.** Number of dispersal events (links) occurring inside (In) and outside (Out) of the clusters for each population and the corresponding ratio between them, "Out/In". Codes for the populations are: W for West, C for Centre, NE for North-East and NW for North-West.</caption><thead><tr><th>Clusters</th><th>Populations</th><th>In</th><th>Out</th><th>Out/In</th></tr></thead><tbody><tr><td rowspan="2">Cluster A</td><td>NE</td><td>78</td><td>34</td><td>0.436</td></tr><tr><td>C</td><td>13</td><td>12</td><td>0.923</td></tr><tr><td rowspan="2">Cluster B</td><td>NW</td><td>52</td><td>8</td><td>0.159</td></tr><tr><td>W</td><td>20</td><td>14</td><td>0.700</td></tr></tbody></table>

doi:10.1371/journal.pone.0166701.t002

Allee effect) and also in populations with a moderate to high density (with birds moving to lower density populations, cf. negative density-dependence effect). Consequently, the dispersal probability could be low in the populations with a moderate density, underlining the U-shaped relationship. Secondly, the observed dispersal rate might not depend so much on the distance or the proximity between two populations but rather might be linked with the migratory behaviour of this species. Thus, the pattern observed in our study, contradictory to previous studies, might largely depend on the migratory axis, which could potentially mask a density effect. White storks coming back from their wintering grounds disperse and settle along the North/South migration axis [16,18]. This might explain why the Northern "source" populations (North-East and North-West) provided more individuals than they received, as the asymmetry and cluster out/in analyses have shown. However, the Western population cannot be considered as a sink population, since its population size increased at a higher rate than the North-Western population, and its population increase might be more likely related to high breeding success and habitat quality. In future analysis, a more biologically realistic picture of the 'French' white stork metapopulation dynamics should be obtained by taking into account international dispersal events: the North-Western population indeed largely extends in western Germany [9], the Eastern cluster is also connected to the reintroduced Swiss population [54], whereas the Southern and Western populations have exchanged some individuals with Spanish populations [17]. Finally, we have shown that the dispersal within a population can have an apparent repulsive effect. Populations with a high internal dispersal were also the ones also had the highest rate of between-populations dispersal. This suggests that the factors influencing dispersal initiation (departure), either within or between populations, most likely are the same [55]. We did not find a density dependence of dispersal events, but one should keep in mind that cost/benefit trade-offs for dispersal are influenced by different factors, not all tested in our study. These factors could be abiotic, such as climate variation [56,57], biotic, such as density, predation pressure or prey distribution [30,53,58], or genetic, like kin competition and inbreeding avoidance [3,59]. Interplay of all these factors can directly affect the metapopulation dynamics [31,55], especially the source/sink role of populations.

In conclusion, the aim of our study was to highlight the asymmetric contribution of different populations in the buildup of a metapopulation, after near extinction. Few studies have simultaneously considered dispersal movements by using gene flow and marked individuals [27,60] that reflect past and more recent dispersal between patches/populations [5]. It was possible with our time-series database to consider the dispersal behavior at both levels (intra- and interpopulational) and thus built the metapopulation network for the different time periods. For conservation purposes, this type of analysis is extremely useful, as it allows identifying the sources from which individuals emigrate to found and build-up new populations. It focuses on populations that receive less and less individuals, requiring special attention in order to guarantee the existence of endangered species populations.
**Supporting Information**

S1 Dataset. Supporting information-metapopulation-white storks.xlsx”
(XLSX)

**Acknowledgments**

Croupe Cigognes France is the consortium of French ornithologists involved in Stork moni-
toring in France, led by G. Wey (collectifcigognesfrance@yahoo.fr), namely G. Balanca, J-P
Baudet, J-Y Brie, V. Burot, P. Carruette, P. Chable, A. Chartier, D. Chevallier, T. Clauss, P.
Delaporte, H. Dugue, J-Y. Fournier, C. Goulevant, S. Merle, F. Noel, J. Pitois, N. Pointecou-
teau, G. Ranvier, T. Roi, N. Vincent-Martin, B. Vollot. All ringing data used in this study have
been collected by qualified, licensed volunteer bird ringers under the coordination of Centre
de Recherches sur la Biologie des Populations d'Oiseaux (CRBPO), Muséum National d'His-
toire Naturelle. All stork ringers and observers are warmly thanked for their contribution to
the long-term monitoring of this species. We also warmly thank J.-C. and C. Barbraud, A.
Fleury, F. Ibanez, J. Pourreau, P. Provost, P. Dagnas, S. Samtmann, G. Paulus, G. Massez and
A. Schierer for their data contribution. No individual affiliation available because all are affili-
ated to GCF (or CRBPO). E. Rojas was supported by the Ministère français de l'Enseignement
et de la Recherche. We also thank M. Enstipp and A. Nesterova for the help provided during
the preparation of the manuscript.

**Author Contributions**

**Conceptualization:** ER SM CS.

**Data curation:** P-YH GW OD.

**Formal analysis:** ER.

**Funding acquisition:** SM.

**Investigation:** ER SM CS P-YH.

**Methodology:** ER CS.

**Resources:** SM.

**Supervision:** SM.

**Visualization:** ER.

**Writing – original draft:** ER SM.

**Writing – review & editing:** ER SM CS PY-H BD.

**References**

1. Hanski I, Gilpin M. Metapopulation dynamics: brief history and conceptual domain. Biol J Linn Soc. 1991; 42: 3–16.

2. Muneepeerakul R, Azaele S, Levin S a., Rinaldo A, Rodriguez-Iturbe Ignacio I. Evolution of dispersal in explicitly spatial metacommunities. J Theor Biol. Elsevier; 2011; 269: 256–265. doi: 10.1016/j.jtbi.2010. 11.003 PMID: 21056581

3. Gandon S. Kin competition, the cost of inbreeding and the evolution of dispersal. J Theor Biol. 1999; 200: 345–64. doi: 10.1006/jtbi.1999.0994 PMID: 10525395

4. Sweanor LL, Logan K a., Hornocker MG. Cougar dispersal patterns, metapopulation dynamics, and conservation. Conserv Biol. 2000; 14: 798–808.
5. Coulon A, Fitzpatrick JW, Bowman R, Stith BM, Makarewich CA, Stenzler LM, et al. Congruent population structure inferred from dispersal behaviour and intensive genetic surveys of the threatened Florida scrub-jay (Aphelocoma caerulescens). Mol Ecol. 2008; 17: 1685–1701. doi: 10.1111/j.1365-294X.2008.03705.x PMID: 18371014

6. Kokko H, López-Sepulcre A. From individual dispersal to species range: perspectives for a changing world. Science (80-). 2006; 313: 789–790. doi: 10.1126/science.1128566 PMID: 16902127

7. Harrison S. Long-Distance Dispersal and Colonization in the Bay Checkerspot Butterfly, Euphydryas Editha Bayensis. Ecology. 1989; 70: 1236–1243.

8. Duquet M. Cigogne blanche. In: Yeatman-Berthelot D, Jarry G, editors. Nouvel atlas des oiseaux nicheurs de France 1985–1989. Société Or. Paris; 1994. pp. 108–111.

9. Schierer A. Quarante années d'observation et de recherche sur la Cigogne Blanche (Ciconia ciconia) en Alsace. Ciconia. 1986; 10: 1–12.

10. Schüz E. Results of the 111 International Census (1974) of the White Stork. Int Counc Bird Preserv Bull. 1979; 13: 173–179.

11. Bairlein F. Population studies of White storks (Ciconia ciconia) in Europe. In: Perrins CM, Lebreton J-D, Hirons GJM, editors. Bird Population Studies: Relevance to Conservation and Management. Oxford, UK: Oxford University Press; 1991. pp. 207–229.

12. Schaub M, Kania W, Köppen U. Variation of primary production during winter induces synchrony in survival rates in migratory white storks Ciconia ciconia. J Anim Ecol. 2005; 74: 656–666.

13. Schierer A. Cigogne blanche. In: Yeatman-Berthelot D, editor. Atlas des Oiseaux de France en Hiver. Societe Or. Paris; 1991. pp. 533–534.

14. Schierer A. Réintroduction de la Cigogne Blanche et Politique des Enclos en Alsace. In: Mériaux J-L, Schierer A, Tombal C, Tombal J-C, editors. Les Cigognes d'Europe, Actes du Colloque International. Metz; 1992. pp. 341–343.

15. Dubois PJ, Le Maréchal P, Olioso G, Yésou P. Inventaires des Oiseaux de France. Paris: Nathan/VUEF; 2001.

16. Barbraud C, Barbraud J-C, Barbraud M. Population dynamics of the White Stork Ciconia ciconia in western France. Ibis. 1999; 141: 469–479.

17. Barbraud J-C, Barbraud C. La nidification de la Cigogne blanche (Ciconia ciconia) en Charente-Maritime de 1978 à 1987. Trajhasse. 1987; 16: 3–7.

18. Rojas ER. Déterminants proximaux et contexte évolutif de la dispersion chez deux espèces d'oiseaux, la Cigogne blanche (Ciconia ciconia) et le Cincle plongeur (Cinclus cinclus). PhD thesis, Université de Strasbourg. 2015.

19. Bunn AG, Urban DL, Keitt TH. Landscape connectivity: A conservation application of graph theory. J Environ Manage. 2000; 59: 265–278.

20. Bret C, Sueur C, Ngoubangoye B, Verrier D, Deneubourg J-L, Petit O. Social structure of a semi-free ranging group of mandrills (Mandrillus sphinx): a social network analysis. PLoS One. 2013; 8: e83015. doi: 10.1371/journal.pone.0083015 PMID: 24340074

21. Pasquaretta C, Levé M, Claidière N, van de Waal E, Whiten A, MacIntosh AJJ, et al. Social networks in primates: smart and tolerant species have more efficient networks. Sci Rep. 2014; 4: 7600. doi: 10.1038/srep07600 PMID: 25534964

22. Shimada M, Sueur C. The importance of social play network for infant or juvenile wild chimpanzees at Mahale Mountains National Park, Tanzania. Am J Primatol. 2014; 76: 1025–1036. doi: 10.1002/ajp.22289 PMID: 24990324

23. Kanngiesser P, Sueur C, Riedl K, Grossmann J, Call J. Grooming network cohesion and the role of individuals in a captive chimpanzee group. Am J Primatol. 2011; 73: 758–767. doi: 10.1002/ajp.20914 PMID: 21698658

24. Bradley B, Doran-Sheehy D, Lukas D. Dispersed male networks in western gorillas. Curr Biol. 2004; 14: 510–513. doi: 10.1016/j.cub.2004.02.062 PMID: 15043817

25. Estrada-Peña A. The relationships between habitat topology, critical scales of connectivity and tick abundance Ixodes ricinus in a heterogeneous landscape in northern Spain. Ecography (Cop). 2003; 26: 661–671.

26. O'Brien D, Manseau M, Fall A, Fortin M-J. Testing the importance of spatial configuration of winter habitat for woodland caribou: An application of graph theory. Biol Conserv. 2006; 130: 70–83.

27. Rozenfeld AF, Arnaud-Haond S, Hernández-García E, Eguíluz VM, Serrão EA, Duarte CM. Network analysis identifies weak and strong links in a metapopulation system. Proc Natl Acad Sci U S A. 2008; 105: 18824–18829. doi: 10.1073/pnas.0805571105 PMID: 19022909
28. Espinoza M, Lédée EJI, Simpfendorfer CA, Tobin, Andrew J, Heupel MR. Contrasting movements and connectivity of reef-associated sharks using acoustic telemetry: implications for management. Ecol Appl. 2015; 25: 2101–2118. PMID: 26910942

29. Gilroy JJ, Lockwood JL. Mate-finding as an overlooked critical determinant of dispersal variation in sex-
ually-reproducing animals. PLoS One. 2012; 7: e38091. doi: 10.1371/journal.pone.0038091 PMID:
22662269

30. Matthysen E. Density-dependent dispersal in birds and mammals. Ecography. 2005. pp. 403–416.

31. Bowler DE, Benton TG. Causes and consequences of animal dispersal strategies: relating individual
behaviour to spatial dynamics. Biol Rev Camb Philos Soc. 2005; 80: 205–25. PMID: 15921049

32. Itonaga N, Köppen U, Plath M, Wallschläger D. Declines in breeding site fidelity in an increasing popula-
tion of White Storks Ciconia ciconia. Ibis. 2011; 153: 636–639.

33. Tobolka M, Kuźniak S, Zolnierowicz KM, Sparks TH, Tryjanowski P. New is not always better: low breeding success and different occupancy patterns in newly built nests of a long-lived species, the white stork Ciconia ciconia. Bird Study. 2013; 60: 399–403.

34. Czajkowski M, Giergiczny M, Kronenberg J, Tryjanowski P. The economic recreational value of a white stork nesting colony: A case of "stork village" in Poland. Tour Manag. 2014; 40: 352-360.

35. Kanyamibwa S, Bairlein F, Schierer A. Comparison of survival rates between populations of the White Stork Ciconia ciconia in central Europe. Ornis Scand. 1993; 24: 297–302.

36. Saether B- E, Grotan V, Tryjanowski P, Barbraud C, Engen S, Fulin M. Climate and spatio-temporal var-
iation in the population dynamics of a long distance migrant, the white stork. J Anim Ecol. 2006; 75: 80–
90. PMID: 16903045

37. Dehorter O, CRBPO. Bird ringing and movement database for France. Centre de Recherches sur la Biologie des Populations d'Oiseaux, Muséum National d'Histoire Naturelle, Paris, France. Accessed: http://crbpo.mnhn.fr/. 2013.

38. Massemin S, Gendner J- P, Samtmann S, Pichegru L, Wulgué A, Le Maho Y. The effect of migration strategy and food availability on White Stork Ciconia ciconia breeding success. Ibis. 2006; 148: 503–508.

39. Tryjanowski P, Kuzniak S. Population size and productivity of the White Stork Ciconia ciconia in relation to common vole Microtus arvalis density. Ardea. 2002; 90: 213–218.

40. Tryjanowski P, Jerzak L, Radkiewicz J. Effect of Water Level and Livestock on the Productivity and Numbers of Breeding White Storks. Waterbirds. 2005; 28: 378–382.

41. Vergara P, I. Aguirre J, Fernández-Cruz M. Arrival date, age and breeding success in white stork Cico-
nia ciconia. J Avian Biol. 2007; 38: 573–579.

42. Aguirre J, Vergara P. Younger, weaker white stork (Ciconia ciconia) nestlings become the best breed-
ers. Evol Ecol Res. 2007; 9: 355–364.

43. Chernetsov N, Chromik W, Dolata P, Profus P, Tryjanowski P. Sex-related natal dispersal of White Storks (Ciconia ciconia) in Poland: how far and where to? Auk. 2006; 123: 1103–1110.

44. Vergara P, Aguirre JI, Fargallo JA, Dávila JA. Nest-site fidelity and breeding success in White Stork Ciconia ciconia. Ibis. 2006; 148: 672–677.

45. Johst K, Brandl R, Pfeifer R. Foraging in a patchy and dynamic landscape: human land use and the white stork. Ecol Appl. 2001; 11: 60–69.

46. Itonaga N, Köppen U, Plath M, Wallschläger D. Breeding dispersal directions in the white stork (Ciconia ciconia) are affected by spring migration routes. J Ethol. 2010; 28: 393–397.

47. Whitehead H. SOCPROG programs: analysing animal social structures. Behav Ecol Sociobiol. 2009; 63: 765–778.

48. Sueur C, Petit O. Organization of group members at departure is driven by social structure in Macaca. Int J Primatol. 2008; 29: 1085–1098.

49. Newman MEJ. Analysis of weighted networks. Phys Rev E—Stat Nonlinear, Soft Matter Phys. 2004; 70.

50. Sueur C, Jacobs A, Amblard F, Petit O, King AJ. How can social network analysis improve the study of primate behavior? Am J Primatol. 2011; 73: 703–719. doi: 10.1002/ajp.20915 PMID: 21181869

51. Mantel N. The Detection of Disease Clustering and a Generalized Regression Approach. Cancer Res. 1967; 27: 209–220. PMID: 6018555

52. Mouquet N, Loreau M. Community patterns in source-sink metacommunities. Am Nat. 2003; 162: 544–557. doi: 10.1086/378857 PMID: 14618534

53. Kim S, Torres R, Drummond H. Simultaneous positive and negative density-dependent dispersal in a colonial bird species. Ecology. 2009; 90: 230–239. PMID: 19294928
54. Schaub M, Pradel R, Lebreton J- D. Is the reintroduced white stork (Ciconia ciconia) population in Swit-zerland self-sustainable? Biol Conserv. 2004; 119: 105-114.

55. Bonte D, Van Dyck H, Bullock JM, Coulon A, Delgado M, Gibbs M, et al. Costs of dispersal. Biol Rev Camb Philos Soc. 2012; 87: 290-312. doi: 10.1111/j.1469-185X.2011.00201.x PMID: 21929715

56. Møller AP, Flensted-Jensen E, Mardal W. Dispersal and climate change: a case study of the Arctic tern Sterna paradisaea. Glob Chang Biol. 2006; 12: 2005-2013.

57. Berg MP, Toby Kiers E, Driessen G, van der Heijden M, Kooi BW, Kuenen F, et al. Adapt or disperse: Understanding species persistence in a changing world. Glob Chang Biol. 2010; 16: 587-598.

58. Leturque H, Rousset F. Dispersal, Kin Competition, and the Ideal Free Distribution in a Spatially Hetero-geneous Population. Theor Popul Biol. 2002; 62: 169-180. PMID: 12167355

59. Bach LA, Thomsen R, Pertoldi C, Loeschcke V. Kin competition and the evolution of dispersal in an indi-vidual-based model. Ecol Modell. 2006; 192: 658-666.

60. Ortego J, García-Navas V, Ferrer ES, Sanz JJ. Genetic structure reflects natal dispersal movements at different spatial scales in the blue tit, Cyanistes caeruleus. Anim Behav. Elsevier Ltd; 2011; 82: 131-137.