- Open Access
The dual role of rivers in facilitating or hindering movements of the false heath fritillary butterfly
Movement Ecologyvolume 3, Article number: 4 (2015)
Species movement responses to landscape structures have been studied using a variety of methods, but movement research is still in need of simple methods that help predicting and comparing movements across structurally different landscapes. We demonstrate how habitat-specific movement models can be used to disentangle causes of differentiated movement patterns in structurally different landscapes and to predict movement patterns in altered and artificial landscapes. In our case study, we studied the role of riparian landscapes to the persistence of the endangered false heath fritillary butterfly (Melitaea diamina) in its newly discovered coastal distribution region in Finland. We compared the movement parameters of the riparian population to two reference populations by using capture-recapture data and habitat-specific diffusion modelling, and analysed the role of the river and riverbank buffer zones in facilitating or hindering false heath fritillary movement with movement simulations.
The riparian population of the false heath fritillary did not show major differences to reference populations in terms of movement parameters within breeding habitat, high-quality matrix and low-quality matrix. However, movement simulations showed that the habitat-specific movement parameters estimated for the false heath fritillary can lead into markedly different movement patterns in structurally different landscapes. An artificial riparian landscape mimicking those of the coastal distribution resulted into more directional, longitudinal movements both parallel and perpendicular to the river than a more mosaic-like landscape, but the existence of the river in the landscape reduced movements across the river.
Our study demonstrates how habitat-specific movement models enable comparisons of movement patterns across structurally different real, altered and artificial landscapes. As such, they can be used to compare movement parameters across populations, to study the effects of management interventions to endangered species and to identify areas that have high sensitivity to individual movement. In our case study, the river is shown to perform a dual role for the movements of the riparian false heath fritillary population. Whereas the river acts as a moderate movement barrier for the false heath fritillary, the longitudinal configuration of riverbank habitats provides a means especially for the male false heath fritillaries to move across the landscape.
Studies on species movement responses to landscape structures increase our understanding on how habitat fragmentation, human-generated landscape structures or management interventions affect individual movement , landscape connectivity or population viability. An example of a landscape structure in which movements have been studied by many researchers are riparian corridors, which provide habitats and migration pathways across human-modified landscapes for a large group of species and thus contribute to the conservation of terrestrial wildlife [2-4]. Riverbank buffer zones, which have historically been left astray because of susceptibility to flooding , are nowadays often protected by or managed according to environmental regulations [5,6]. Riverbanks are also influenced by high levels of soil moisture  and flood-based natural disturbance dynamics that can maintain a continuum of successional habitats [8-11]. Thus they may provide a variety of habitats ranging from riparian forest stripes  to wet open meadows .
Species occurrence and movements in riverbank habitats have been studied in several taxa, including moths , spiders , forest birds [12,15] and mammalian predators . Among butterflies, many species inhabiting moist open habitats have been found to occur or move along riverbanks, like the Monarch butterflies , the Clouded Apollo  and the St. Francis satyr butterfly . Study methods have included direct observations of movement [12,16,20,21], capture-recapture methods , spatiotemporal patterns in species sightings during migration periods , densities of species sightings at riverbanks versus other habitats [13,15,18] and population genetic studies [14,19]. Other studies have shown that the ability and willingness to cross rivers vary between species and taxa [22-25] and that the hostility of the river to the species has an effect on the likelihood of river crossings [26-28]. Relocation and mobbing call experiments have been used to study movements parallel and perpendicular to rivers in birds [29-32], demonstrating differences in species movement along rivers versus across them  and indicating the dual role of rivers in the formation of both movement barriers and functional movement corridors [34-37] for many terrestrial species.
Experimental data on species movement responses to tested landscape structures can be analysed by straightforward statistical analyses and can give reliable estimates on the effects of such structures on species movement. In contrast, making inferences on species’ movement in alternative landscapes is often not possible, since either the data lack the overall landscape context or the analysis model does not estimate the overall decision-making logic of the individuals studied. Movement researchers have expressed a need for a larger variety of methods that enable predicting species’ movement patterns in structurally different landscapes based on observed data [32,36] besides the existing potential methods and studies [36-39].
Habitat-specific movement models can be used to infer species’ movement parameters indirectly from capture-recapture data and habitat type maps of the study landscape [40-42]. In this study, we demonstrate how they can be used to disentangle causes of differentiated movements in structurally different landscapes and to predict movements in altered and artificial landscapes. Habitat-specific diffusion models have been successfully applied to non-territorial species with a presumably relatively simple movement logic, e.g. to some species of butterflies, the movement of which can be approximated by random walk at sufficiently large spatial scales [40-44]. Movement rate, measured by the diffusion parameter, is often expected to differ between at least three landscape types: breeding habitat (BH), high-quality matrix (HQM) and a more hostile low-quality matrix (LQM) [41,45,46]. Additionally, butterflies typically show edge-mediated behaviour [40,47-51], i.e. habitat selection at the edges between any two habitat types, which behaviour can be implemented in diffusion models by assuming that the probability of the butterfly being on preferred side of the edge is k times higher than the probability of it being in the other side of the edge, where the parameter k measures relative habitat preference .
The aim of our case study was to find out whether the structural configuration of riparian landscapes could be a key factor that explains the persistence of the endangered false heath fritillary butterfly Melitaea diamina (Lang, 1789) in its newly discovered coastal distribution region in Finland. The false heath fritillary is a moist meadow specialist that has, like many butterfly species, suffered from agricultural modernization and drainage of moist soils within the recent decades . In Finland, sightings of the false heath fritillary have been made lately only in two distinctive regions of a narrow geographical range (Figure 1). The inland population has been subjected to various conservation measures, whereas the distribution along the West Coast of Finland was only properly revealed by a monitoring project in 2009–2012. It has remained unclear what has provided the means of persistence for the species at its newly discovered distribution.
We hypothesised that (H1) the coastal populations of the false heath fritillary are similar to the inland populations of the species in terms of habitat-specific movement parameters, but that (H2) the longitudinally structured riparian landscapes of the coastal distribution result into more directional, longitudinal movements than the mosaic-like landscapes of the inland populations. We assumed that the resulting connectivity across larger spatial scales could decrease metapopulation susceptibility to the effects of local climate variations, thus enforcing metapopulation persistence in the riparian landscapes [54,55]. To test the first part of our hypothesis (H1), we compared the habitat-specific movement parameters of the riparian coastal population to those of two inland reference populations (Figure 1A-C, Table 1) via posterior comparisons and by contrasting the riparian capture-recapture data against data simulated for the riparian landscape based on the reference movement models. Having verified the overall similarity of movement in the three populations, we tested the second part of our hypothesis (H2) by creating a generalised movement model for the false heath fritillary in Finland, using it to simulate movements in the riparian landscape and in altered and artificial landscapes and testing whether riparian landscapes would result into more directional, longitudinal movements than more mosaic-like landscapes.
Results and discussion
Comparison of movement parameters across study areas
Posterior comparisons of the model parameters (Table 2, Figure 2) did not show major differences in habitat-specific movement parameters between the riparian population and the two inland reference populations. In comparison to the reference populations, the riparian population shows a trend towards lower movement rate in the low-quality matrix (D LQM ), higher preference for the high-quality matrix (k HQM ), higher female mortality (m) and higher capture probability of females (p). However, the statistical support for these differences remained weak (Table 3), which is partly due to the large amount of posterior uncertainty in the riparian movement model (Figure 2). This was in turn caused by low population density during the year of the capture-recapture experiment in the riparian population, resulting in a low number of individuals and recaptures relative to the number of days sampled (Table 1). Thus, direct comparisons of movement parameters among the three populations were consistent, without any strong evidence for population specific movement parameters.
False heath fritillary movement parameters proved similar across populations also when the riparian capture-recapture data was contrasted against data simulated for the riparian landscape based on the reference movement models (Figures 3 and 4, Table 4). Among the test statistics we used, the distribution of times between marking and last recapture (Figure 3) acts as a proxy for life-span, whereas the observed total movement distance (Figure 4) acts as a proxy for overall movement activity. The exception for the good fit between the riparian data and the predictions of the models parameterised by the reference populations is the number of females that were observed only during one day or which moved less than 50 m in total between marking and last recapture, which number was underestimated by the model predictions (Figures 3 and 4, Table 4).
The above results suggest that false heath fritillary populations are not differentiated in a major way in terms of their movement behaviour across the two distribution regions. The relatively small observed differences across distribution regions are more likely to be caused by variation in the habitat types (which were classified here in three very broad categories), or differences in weather conditions during the capture-recapture experiments, rather than e.g. genetic or other biologically important differences.
Movement parameters in the two sexes and across habitat types
A generalised movement model for the false heath fritillary in Finland, created by pooling the riparian and REF1 data sets, provided relatively narrow posterior distributions especially for mortality, capture probability and habitat type preference (Figure 2). It differed from the model parameters of REF1 (Figure 2) mostly for the movement rate in the breeding habitat (D BH ), for which the original REF1 data did not contain information due to its capture-recapture study design . According to this model, males moved faster than females in the breeding habitat (D BH ) and high-quality matrix (D HQM ), and had higher mortality rates (m) and capture probabilities (p; Figure 2 and Table 3), suggesting a median expected life-time (1/m) of ten days for females and seven days for males. Both sexes showed higher preference for and had lower movement rates in the breeding habitat in comparison to the high-quality matrix and showed higher preference for the high-quality matrix in comparison to the low-quality matrix . Females also followed the expectation  that the rate of movement is faster in the low-quality matrix than in the breeding habitat . These differences are in line with the previous view of the false heath fritillary as a rather sedentary butterfly, with males patrolling for females for mating and thus being slightly more mobile than females [56-60].
In the generalised movement model and in the riparian movement model, estimated habitat preference (relative to breeding habitat) for the river (k R ) were markedly lower in comparison to the preference of high-quality matrix (k HQM ); there was only a small overlap between the posteriors in both models and a nearly hundred-fold difference between the medians in the riparian model (Figure 2). Estimated habitat preference for the river (k R ) also had slightly lower medians and narrower posterior distributions in comparison to those of the low-quality matrix (k LQM ; Figure 2), but these differences did not gain much statistical support in posterior comparisons . In line with this, the number of river crossings in the riparian data was below the medians predicted by the reference movement models, in which movements across the river were predicted as if the river was part of the low-quality matrix (Table 4), but nevertheless within the 95% posterior credibility intervals. The results thus suggest that the river fits well within the low-quality matrix habitat type in terms of movement rates (D) and boundary responses (k).
A movement simulation based on the posterior medians of the generalised movement model (Figure 5) illustrates the sedentary nature of the false heath fritillary butterflies, individuals being likely to spend most of their time along a network of connected habitat patches (Figure 5, dark colours). Individuals starting flight at a well-connected patch (panels C-D) would be more likely to spend time in the habitat patch network than individuals starting movement at a more peripheral patch (panels A-B), who would end up spending some of their time in an area of the low-quality matrix that did not contain habitat patches (upper right). Female movement is somewhat more confined to the habitat patches than that of the males. These patterns highlight the importance of the well-connected habitat patch network provided by the riverbank buffer zones in the riparian landscape.
False heath fritillary movement patterns in structurally different landscapes
Movement simulations carried out using three artificial landscapes (Figure 6) demonstrate how different landscape elements interplay to inhibit or facilitate false heath fritillary movements in a riparian landscape. First, we compare false heath fritillary movement patterns between a riparian landscape consisting of a low-quality matrix river and symmetric riverbank habitats (Figure 6, panels D-F), and a control landscape where an equal amount of habitat and low-quality matrix is distributed evenly across the landscape (panels A-C). For male false heath fritillaries, the probabilities of hitting a target patch both from the North (parallel to the river) and from the East and West (perpendicular to the river) are higher in the riparian landscape than in a control landscape (Figure 7, panels C-D). Long directional movements perpendicular to the river occur due to the absence of habitat patches in the surrounding landscape, causing the males to fly fast across the low-quality matrix (D LQM ) while searching for habitat patches. If habitat patches are distributed evenly across the landscape, males are likely to stop at habitat patches and may change direction when leaving a patch. Long directional movements parallel to the river are co-products of relatively high movement rates in, and a high preference of, the breeding habitat (D BH and k BH ). For female false heath fritillaries, differences in movement patterns are otherwise similar to those of males, except that probabilities of hitting a target patch from the North (parallel to the river) are not higher in the riparian landscape than in the control landscape due to lower estimated movement rates in the breeding habitat (D BH , Figure 8). Hitting probabilities diagonal to the river are likely to be higher than those parallel to the river for females (Figure 6, panel F).
We next compare false heath fritillary movement patterns between the riparian landscape (Figure 6, panels D-F) and a similar landscape where the river has been replaced with breeding habitat (panels G-I). The probabilities of hitting a target patch from the East (perpendicular to the river) are lower for both male and female false heath fritillaries in the riparian landscape (Figures 7 and 8, panel C) due to the low preference of entering the low-quality matrix environment (k LQM ). For male false heath fritillaries, the probabilities of hitting the target patch from the North and from the West are also lower (Figure 7, panel D), probably because of faster movement in (D LQM ) and lower preference of (k LQM ) the river environment; individuals that enter the river environment are either likely to fly fast past the target patch or are not likely to cross the river for another time.
Movement simulation carried out in the riparian landscape of the riparian study area (Figure 9) visualises regions where nuanced changes in the species’ movement parameter estimates would have largest effects to estimated hitting probabilities. While differences are small between a landscape where the river was considered a distinct habitat type with its own (median) movement rates and boundary responses (panels A-B) and a landscape where the river is modelled as part of low-quality matrix (panels C-D), a more detailed comparison (panels E-F) demonstrates that they are likely to be strongest for individuals that start movement on the opposite side of the river from the target patch (decreasing their probability of hitting the target patch) and for individuals starting from the same side as the target patch but close to the opposite riverbank (increasing their possibility to hit the target patch). This demonstrates the role of the riverbends of the riparian landscapes in increasing and decreasing connectivity between habitat patches.
The above simulations demonstrate that a longitudinally structured riparian landscape performs a dual role for the movements of the riparian false heath fritillary population. First, as hypothesized, longitudinally structured riparian landscapes result into more directional, longitudinal movements than the mosaic-like landscapes, thus generating habitat connectivity for the riparian false heath fritillaries. However, the river itself acts as a structural movement barrier for the false heath fritillaries, reducing movements across the river. Such dual effect has implications for conservation planning: low probabilities of river crossings decrease connectivity between the two sides of the river but riverbank habitats can increase connectivity within each of the sides .
Habitat-specific movement models as disentanglers of movement patterns from landscape structure
Our analyses show that habitat-specific movement parameters can lead into markedly different movement patterns in structurally different landscapes for some species. This is likely to be the case especially for habitat specialists, like the false heath fritillary, that spend most of their time in their breeding habitat. Habitat-specific movement models make it possible to predict the movements of such species in structurally different landscapes, and thus can be used to study the effects of management interventions, e.g. the placement of conservation sites, to the movements of endangered species and to identify areas in the landscape that have high sensitivity to their movement. The possibility to use also artificial landscapes for studying movements enable studying movements in structurally simple, symmetrical landscapes, which clarify the often complex effects of multiple landscape elements to movements, and series of slightly altered landscapes, which enables quantifying the effects landscape structure to individual movements.
Habitat-specific movement models can be parameterised using spatially explicit capture-recapture data, which kind of data are difficult to analyse with many other methods because it tells about movement behaviour only indirectly. Such data are often the only possibility for movement analysis for threatened species, for which relocation experiments may not be possible due to legislative restrictions. The indirect nature of capture-recapture data is an obvious drawback, both in terms of limited informativeness of the data as well as the need for more sophisticated analysis methods. For instance, cultivated fields and closed forests may be quite different environments for butterfly movement, but we have pooled them in the broader category of low-quality matrix based on prior analysis results , as the capture-recapture data lack resolution to disentangle if movements would differ among these habitat types. Similarly, it might be the case that the river environment differs from closed forests and cultivated fields in terms of butterfly movement, even though we could not prove such a difference in this study. The existence of such a difference would be most efficiently studied by relocation experiments, where butterflies would be followed to directly study their behaviour at habitat type edges, which unfortunately were not possible for our study species.
Our study demonstrates how habitat-specific movement models and already a simple representation of landscape structure enable comparisons of movement patterns across structurally different real, altered and artificial landscapes. As such, they can be used to compare movement parameters across populations, to study the effects of management interventions and to identify areas that have high sensitivity to individual movement.
In our case study, the river is shown to perform a dual role for the movements of the riparian false heath fritillary population. As habitat specialists , the false heath fritillaries have a high tendency to stay at habitat patches and thus their movement is very much dictated by the spatial configuration of their habitat patch network. Whereas rivers act as moderate movement barriers for the false heath fritillary, the longitudinal configurations of riverbank meadows provide a means especially for the male false heath fritillaries to move across the landscape.
The study sites, capture-recapture data and habitat classifications
The reference populations REF1 (located in Siitama; 61.6°N, 24.2°E) and REF2 (Sorila; 61.55°N, 23.9°E) belong to the same metapopulation system [56-58] and are located 15 km from each other near the city of Tampere. The REF1 study site contains a dense cluster of 14 habitat patches and other meadows classified as high-quality matrix within a landscape of forests and cultivated fields (Figure 1A). The REF2 study site is characterized by more sparsely located stripes of habitat patches and high-quality matrix, including open areas such as a powerline right-of-way at the Eastern edge (Figure 1B). The riparian study site at Merikarvia (61.86°N, 21.56°E) is characterized by the ~30-50 m wide Merikarvia River that twists across the landscape (Figure 1C), with habitat patches and high-quality matrix along riverbanks and some powerlines. At all study sites, monthly mean temperatures vary from −8°C to 16.5°C, the mean annual rainfall is approximately 650 mm and a permanent snow cover lasts for 3–5 months .
Habitat-specific movement analyses were based on habitat classification maps (Figure 1) that categorized the landscape into breeding habitat (open or semi-open meadows with the host plant Valeriana sambucifolia), high-quality matrix (open areas with nectar plants), low-quality matrix (e.g. closed forests and cultivated fields) and the river. Butterfly searching took place within all three terrestrial landscape classes based on a discrete set of search sites, so that the search effort could be described as a table of sites searched each day.
The capture-recapture data for REF1 has been described by Wahlberg  and Wahlberg et al. [56,60] and the original habitat classifications by Ovaskainen [40; model B]. In REF2, capture-recapture data was originally collected and analysed by Ovaskainen and Cabeza for a non-refereed conservation study . In REF2, 73, 12 and 27 polygons within the breeding habitat, the high-quality matrix and the low-quality matrix respectively with an average size of 0.57 ha were used as search sites. At this time point, the open areas (O) and forests (F) of REF1 were reclassified into areas of high-quality matrix and low-quality matrix (Figure 1). In the riparian population, we collected data using the EarthCape software . We split the previously delineated 21 large habitat patches into search polygons, each with a diameter of 62 m in minimum, and placed 20 longitudinal search areas of width 15 m into the remaining landscape for the capture-recapture study. Six of these search areas were categorised as high-quality matrix and 12 as low-quality matrix during habitat classification. We estimated population sizes during the sampling periods using the open population model of Rcapture . To account for births and deaths occurring during the sampling periods, and to find a model with smallest residuals, we grouped the capture-recapture data into two-day (REF2, RIPARIAN) or four-day (REF1) primary periods and, as instructed , removed outliers with standardized residuals greater than four before population size estimation. The details of the capture-recapture studies and the resulting data sets  are summarised in Table 1.
Estimation of movement parameters
We applied habitat-specific diffusion modelling with Bayesian inference [40,41] to estimate posterior distributions of the model parameters separately for males and females. We triangulated the habitat classification and capture-recapture site maps with the Mapper software  and analysed the triangulated maps, the search effort matrix and the observation matrix with the Disperse software . Disperse computes the likelihood of the data using the finite element method, and estimates posterior distributions via adaptive MCMC methods.
The biological parameters estimated were the habitat-specific diffusion coefficients D measuring the rates of individual movement, the habitat preference k for each habitat type relative to that in the breeding habitat (k BH = 1 by definition), and the mortality rate m. The observation model involves the parameter capture probability p, i.e. the probability of capturing an individual from the search site conditional on the individual being present. We used lognormal priors for k, D and m and a logit-normal prior for p. The prior medians were derived from a long-term study of a closely related species Melitaea cinxia , and we assumed wide credibility intervals to account for interspecies differences. We assumed the same priors for the river as for the low-quality matrix. To create the generalised movement model (GEN) for the false heath fritillary across distribution regions, we estimated the movement parameters for the riparian data using the posteriors of the most data-rich REF1 population (Figure 2, (Additional file 1)) as priors, which corresponds to the estimation of the model parameters from the combined data sets.
Comparison and cross-validation of parameter estimates
To assess the similarity of posterior distributions among the three populations and between the two sexes we calculated the posterior probability of difference by drawing 10000 pairs of random samples from each pair of the posterior distributions to be compared, as in Ovaskainen et al. . Posterior comparisons between the two sexes were based on the generalised movement model.
To test whether the movement models of the reference populations would correctly predict false heath fritillary movements in the riparian landscape, we generated posterior predictive data by drawing 1000 samples from the reference model posterior distribution and simulated movements for each sample, assuming for the river (which was absent from the reference populations) the parameters of the low-quality matrix. We set the simulations to start at the locations of real butterfly marking events at the riparian landscape, generated a movement track for each individual, and generated capture-recapture data assuming the same spatiotemporal search effort as in the field study. We adjusted the capture probability p for both males and females so that the proportions of recaptured individuals matched with the real proportions in the riparian data. We compared simulated data against the real data with respect to the following parameters: the number of river crossings, the distribution of total distance moved by recaptured individuals, the distribution of times from marking until last recapture, and the ratio of migrated individuals. To retrieve parameter values for all individuals, we generated 1000 sums of random pairs of male and female results.
To analyse the roles of the high-quality matrix and the river for false heath fritillary movements, we created movement simulations for the riparian landscape based on the generalised movement model by using both the original, altered and artificial landscape maps. Simulated individuals either started flight from a pre-defined natal patch, after which we calculated the occupancy time density u(y), the time that the individual was expected to spend at any location y of the study area during its lifetime . Alternatively, we calculated the hitting probability q(y), the probability that the butterfly starting its flight from any location y would reach a target patch during its lifetime . We visualised the effect of the landscape alterations by plotting the difference of q(y) between the original and the altered landscape at each map location y. We created credibility intervals for the movement statistics by generating them for 100 parameter combinations sampled from the posterior distribution.
Availability of supporting data
The data sets supporting the results of this article are available in the Dryad repository, doi:10.5061/dryad.j54vv, [http://doi.org/10.5061/dryad.j54vv].
Nathan R, Getz WM, Revilla E, Holyoak M, Kadmon R, Saltz D, et al. A movement ecology paradigm for unifying organismal movement research. PNAS. 2008;105:19052–9.
Bentrup G, Kellerman T. Where should buffers go? Modeling riparian habitat connectivity in northeast Kansas. J Soil Water Conserv. 2004;59:209–15.
Ives CD, Hose GC, Nipperess DA, Taylor MP. Environmental and landscape factors influencing ant and plant diversity in suburban riparian corridors. Landsc Urban Plan. 2011;103:372–82.
Dallimer M, Rouquette JR, Skinner AMJ, Armsworth PR, Maltby LM, Warren PH, et al. Contrasting patterns in species richness of birds, butterflies and plants along riparian corridors in an urban landscape. Divers Distrib. 2012;18:742–53.
European Commission. The European Union’s Biodiversity Action Plan: halting the loss of biodiversity by 2010 – and beyond. Luxembourg: Office for Official Publications of the European Communities; 2008 [http://ec.europa.eu/environment/nature/info/pubs/docs/brochures/bio_brochure_en.pdf]
Gerard M, El Kahloun M, Rymen J, Beauchard O, Meire P. Importance of mowing and flood frequency in promoting species richness in restored floodplains. J Appl Ecol. 2008;45:1780–9.
Henry AC, Hosack DA, Johnson CW, Rol D, Bentrup G. Conservation corridors in the United States: benefits and planning guidelines. J Soil Water Conserv. 1999;54:645–50.
Naiman RJ, Melillo JM, Lock MA, Ford TE, Reice SR. Longitudinal patterns of ecosystem processes and community structure in a subarctic river continuum. Ecology. 2014;68:1139–56.
Jäkäläniemi A, Tuomi J, Siikamäki P. Conservation of species in dynamic landscapes: divergent fates of Silene tatarica populations in riparian habitats. Conserv Biol. 2006;20:844–52.
Tabacchi E, Correll DL, Hauer R, Pinay G, Planty-Tabacchi AM, Wissmar RC. Development, maintenance and role of riparian vegetation in the river landscape. Freshw Biol. 1998;40:497–516.
Januschke K, Brunzel S, Haase P, Hering D. Effects of stream restorations on riparian mesohabitats, vegetation and carabid beetles. Biodivers Conserv. 2011;20:3147–64.
Gillies CS, St Clair CC. Riparian corridors enhance movement of a forest specialist bird in fragmented tropical forest. PNAS. 2008;105:19774–9.
Mönkkönen M, Mutanen M. Occurrence of moths in boreal forest corridors. Conserv Biol. 2003;17:468–75.
Lambeets K, Breyne P, Bonte D. Spatial genetic variation of a riparian wolf spider Pardosa agricola (Thorell, 1856) on lowland river banks: the importance of functional connectivity in linear spatial systems. Biol Conserv. 2010;143:660–8.
Machtans CS, Villard M, Hannon SJ. Use of riparian buffer strips as movement corridors by forest birds. Conserv Biol. 1996;10:1366–79.
Hilty JA, Merenlender AM. Use of riparian corridors and vineyards by mammalian predators in Northern California. Conserv Biol. 2004;18:126–35.
Dingle H, Zalucki MP, Rochester WA, Armijo-Prewitt T. Distribution of the monarch butterfly, Danaus plexippus (L.) (Lepidoptera: Nymphalidae), in western North America. Biol J Linn Soc. 2005;85:491–500.
Meier K, Kuusemets V, Luig J, Mander Ü. Riparian buffer zones as elements of ecological networks: case study on Parnassius mnemosyne distribution in Estonia. Ecol Eng. 2005;24:531–7.
Milko LV, Haddad NM, Lance SL. Dispersal via stream corridors structures populations of the endangered St. Francis’ satyr butterfly (Neonympha mitchellii francisci). J Insect Conserv. 2012;16:263–73.
Shirley SM. Movement of forest birds across river and clearcut edges of varying riparian buffer strip widths. For Ecol Manag. 2006;223:190–9.
Kuefler D, Hudgens B, Haddad NM, Morris WF, Thurgate N. The conflicting role of matrix habitats as conduits and barriers for dispersal. Ecology. 2010;91:944–50.
Lamborot M, Eaton L. The Maipo River as a biogeographical barrier to Liolaemus monticola (Tropiduridae) in the mountain ranges of central Chile. J Zool Syst Evol Res. 1997;35:105–11.
Puth LM, Wilson KA. Boundaries and corridors as a continuum of ecological flow control: lessons from rivers and streams. Conserv Biol. 2001;15:21–30.
Díaz-Muñoz SL. Role of recent and old riverine barriers in fine-scale population genetic structure of Geoffroy’s tamarin (Saguinus geoffroyi) in the Panama Canal watershed. Ecol Evol. 2012;2:298–309.
Juen L, Marco P. Dragonfly endemism in the Brazilian Amazon: competing hypotheses for biogeographical patterns. Biodivers Conserv. 2012;21:3507–21.
Hayes FE, Sewlal JN. The Amazon River as a dispersal barrier to passerine birds: effects of river width, habitat and taxonomy. J Biogeogr. 2004;31:1809–18.
Matsuda I, Tuuga A, Akiyama Y, Higashi S. Selection of river crossing location and sleeping site by proboscis monkeys (Nasalis larvatus) in Sabah, Malaysia. Am J Primatol. 2008;70:1097–101.
Chilena R, Natural DEH. River crossings by Ateles geoffroyi and Alouatta pigra in southern Mexico: a preliminary report. Rev Chil Hist Nat. 2010;83:435–42.
Tremblay MA, St. Clair CC. Factors affecting the permeability of transportation and riparian corridors to the movements of songbirds in an urban landscape. J Appl Ecol. 2009;46:1314–22.
Tremblay MA, St. Clair CC. Permeability of a heterogeneous urban landscape to the movements of forest songbirds. J Appl Ecol. 2011;48:679–88.
St. Clair CC. Comparative Permeability of Roads, rivers, and meadows to songbirds in Banff National Park. Conserv Biol. 2003;17:1151–60.
Bélisle M, St. Clair CC. Cumulative effects of barriers on the movement of forest birds. Conserv Ecol. 2002;5:9.
Kuefler D, Haddad NM. Local versus landscape determinants of butterfly movement behaviors. Ecography. 2006;29:549–60.
Dennis RLH, Dapporto L, Dover JW, Shreeve TG. Corridors and barriers in biodiversity conservation: a novel resource-based habitat perspective for butterflies. Biodivers Conserv. 2013;22:2709–34.
Hess G, Fischer R. Communicating clearly about conservation corridors. Landsc Urban Plan. 2001;55:195–208.
Chetkiewicz CB, St. Clair CC, Boyce MS. Corridors for conservation: integrating pattern and process. Annu Rev Ecol Evol Syst. 2006;37:317–42.
Levey DJ, Bolker BM, Tewksbury JJ, Sargent S, Haddad N. Effects of landscape corridors on seed dispersal by birds. Science. 2005;309:146–8.
Adriaensen F, Chardon JP, De Blust G, Swinnen E, Villalba S, Gulinck H, et al. The application of ‘least-cost’ modelling as a functional landscape model. Landsc Urban Plan. 2003;64:233–47.
Goulart FF, Salles P, Machado RB. How may agricultural matrix intensification affect understory birds in an Atlantic Forest landscape? A qualitative model on stochasticity and immigration. Ecol Inform. 2013;18:93–106.
Ovaskainen O. Habitat-specific movement parameters estimated using mark-recapture data and a diffusion model. Ecology. 2004;85:242–57.
Ovaskainen O, Luoto M, Ikonen I, Rekola H, Meyke E, Kuussaari M. An empirical test of a diffusion model: predicting clouded Apollo movements in a novel environment. Am Nat. 2008;171:610–9.
Ovaskainen I, Rekola H, Meyke E, Arjas E. Bayesian methods for analysing movements in heterogeneous landscapes from mark-recapture data. Ecology. 2008;89:542–54.
Patlak C. Random walk with persistence and external bias. Bull Math Biophys. 1953;15:311–38.
Turchin P. Quantitative analysis of movement: measuring and modeling population redistribution in animals and plants. Sunderland, MA: Sinauer; 1998.
Bennett VJ, Pack SM, Smith WP, Betts MG. Sex-biased dispersal in a rare butterfly and the implications for its conservation. J Insect Conserv. 2013;17:949–58.
Severns PM, McIntire EJB, Schultz CB. Evaluating functional connectivity with matrix behavior uncertainty for an endangered butterfly. Landsc Ecol. 2013;28:559–69.
Kuussaari M, Nieminen M, Hanski I. An experimental study of migration in the Glanville fritillary butterfly Melitaea cinxia. J Anim Ecol. 1996;65:791–801.
Ries L, Debinski DM. Butterfly responses to habitat edges in the highly fragmented prairies of Central Iowa. J Anim Ecol. 2001;70:840–52.
Schultz CB, Crone EE. Edge-mediated dispersal behavior in a prairie butterfly. Ecology. 2001;82:1879–92.
Schtickzelle N, Baguette M. Behavioural responses to habitat patch boundaries restrict dispersal and generate emigration-patch area relationships in fragmented landscapes. J Anim Ecol. 2003;72:533–45.
Maciel GA, Lutscher F. How individual movement response to habitat edges affects population persistence and spatial spread. Am Nat. 2013;182:42–52.
Ovaskainen O, Cornell SJ. Biased movement at a boundary and conditional occupancy times for diffusion processes. J Appl Probab. 2003;40:557–80.
van Swaay C, Warren M, Lois G. Biotope use and trends of European butterflies. J Insect Conserv. 2006;10:189–209.
Hanski I. Metapopulation dynamics. Nature. 1998;396:41–9.
Poos MS, Jackson DA. Impact of species-specific dispersal and regional stochasticity on estimates of population viability in stream metapopulations. Landsc Ecol. 2012;27:405–16.
Wahlberg N, Moilanen A, Hanski I. Predicting the occurrence of endangered species in fragmented landscapes. Science. 1996;273:1536–8.
Moilanen A, Cabeza M. Single-species dynamic site selection. Ecol Appl. 2002;12:913–26.
Cabeza M, Arponen A, Jäättelä L, Kujala H, van Teeffelen A, Hanski I. Conservation planning with insects at three different spatial scales. Ecography. 2010;33:54–63.
Wahlberg N. The life history and ecology of Melitaea diamina (Nymphalidae) in Finland. Nota Lepidoptera. 1997;20:70–81.
Wahlberg N, Klemetti T, Selonen V, Hanski I. Metapopulation structure and movements in five species of checkerspot butterflies. Oecologia. 2002;130:33–43.
Ovaskainen O, Cabeza M. Tummaverkkoperhonen Sorilassa: Näkökulmia alueen kaavoitusprosessiin. 2007 [http://www.tampere.fi/material/attachments/t/5m2jqu1Cg/tummaverkkoperhonensorilassa290507.pdf]
Krauss J, Steffan-dewenter I, Tscharntke T. How does landscape context contribute to effects of habitat fragmentation on diversity and population density of butterflies? J Biogeogr. 2003;30:889–900.
Ilmatieteen laitos. Suomen maakuntien ilmasto. In: Kersalo J, Pirinen P, editors. [Ilmatieteen laitos: Raportteja, vol. 2009:8.]. Helsinki: Yliopistopaino; 2009. p. 1–185.
Baillargeon S, Rivest L-P. Rcapture: loglinear models for capture-recapture in R. J Stat Softw. 2007;19:1–31.
Cabeza M, Fabritius H, Lehtomäki J, Ovaskainen O, Rönkä K, Wahlberg N. Data from: The dual role of rivers in facilitating or hindering movements of the false heath fritillary butterfly. Movement Ecol. 2015. [http://doi.org/10.5061/dryad.j54vv]
Ojanen SP, Nieminen M, Meyke E, Pöyry J, Hanski I. Long-term metapopulation study of the Glanville fritillary butterfly (Melitaea cinxia): survey methods, data management, and long-term population trends. Ecol Evol. 2013;3:3713–37.
Ovaskainen O. Analytical and numerical tools for diffusion based movement models. Theor Popul Biol. 2008;73:198–211.
We thank three anonymous referees for comments to the manuscript, Mar Cabeza for the REF2 data and comments to the manuscript, Niklas Wahlberg for the REF1 data, Joona Lehtomäki, Sanna Mäkeläinen, Tinto Aaltonen, Jussi Iso-Tuisku, Marko Schrader, Tiina Avomaa, Veera Piironen and Anttu Lipponen for fieldwork, and Evgeniy Meyke for help with EarthCape and the Mapper software. The Academy of Finland (Grant no. 250444 to OO), the Doctoral Programme in Wildlife Biology Research, Helsingin yliopiston Professorien rouvat ry., the Kuopio Naturalists’ society, Societas pro Flora et Fauna Fennica and the Entomological Society of Helsinki are thanked for funding.
The authors declare that they have no competing interests.
HF took part in the collection of RIPARIAN data, took part in parameter estimation, model validation, movement simulations and the creation of figures, created both the movement simulations and figures for the artificial landscapes and wrote the manuscript. KR took part in the collection of RIPARIAN data, prepared habitat classifications for the riparian landscape and took part in the estimation of movement parameters for the riparian and generalised movement models, model validation and in the preparation of the movement simulations. OO participated in the design of fieldwork in the REF2 landscape, estimated movement parameters for REF2 and supervised the RIPARIAN study. All authors have contributed to and approved the final manuscript.
Prior and posterior distributions of the false heath fritillary. An Excel file. The medians (and 95% credibility intervals) of the posterior distributions for false heath fritillary males and females across the three study areas. RIPARIAN refers to the riparian population, REF1 and REF2 to the two reference populations and GEN to the generalised movement model (see “Methods”). The parameters k (dimensionless), D (m2 day−1), m (day−1 ), and p (dimensionless) refer to relative habitat preference, movement rate, mortality and capture probability, respectively. The Subscripts BH, HQM, LQM and R stand for breeding habitat, high-quality matrix, low-quality matrix and the river, respectively (see “Methods”).