<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Ecol. Evol.</journal-id>
<journal-title>Frontiers in Ecology and Evolution</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Ecol. Evol.</abbrev-journal-title>
<issn pub-type="epub">2296-701X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fevo.2023.1078649</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Ecology and Evolution</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>An objective approach to select surrogate species for connectivity conservation</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Dutta</surname><given-names>Trishna</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>*</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1985190"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>De Barba</surname><given-names>Marta</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<xref ref-type="aff" rid="aff5"><sup>5</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Selva</surname><given-names>Nuria</given-names>
</name>
<xref ref-type="aff" rid="aff6"><sup>6</sup></xref>
<xref ref-type="aff" rid="aff7"><sup>7</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/408291"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Fedorca</surname><given-names>Ancuta Cotovelea</given-names>
</name>
<xref ref-type="aff" rid="aff8"><sup>8</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/313981"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Maiorano</surname><given-names>Luigi</given-names>
</name>
<xref ref-type="aff" rid="aff9"><sup>9</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/107332"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Thuiller</surname><given-names>Wilfried</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1052212"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zedrosser</surname><given-names>Andreas</given-names>
</name>
<xref ref-type="aff" rid="aff10"><sup>10</sup></xref>
<xref ref-type="aff" rid="aff11"><sup>11</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Signer</surname><given-names>Johannes</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Pfl&#xfc;ger</surname><given-names>Femke</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff12"><sup>12</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2067835"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Frank</surname><given-names>Shane</given-names>
</name>
<xref ref-type="aff" rid="aff10"><sup>10</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Lucas</surname><given-names>Pablo M.</given-names>
</name>
<xref ref-type="aff" rid="aff6"><sup>6</sup></xref>
<xref ref-type="aff" rid="aff9"><sup>9</sup></xref>
<xref ref-type="aff" rid="aff13"><sup>13</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2067489"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Balkenhol</surname><given-names>Niko</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/210328"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Wildlife Sciences, Faculty of Forest Sciences and Forest Ecology, University of Goettingen</institution>, <addr-line>G&#xf6;ttingen</addr-line>, <country>Germany</country></aff>
<aff id="aff2"><sup>2</sup><institution>Resilience Programme, European Forest Institute, Platz der Vereinten Nationen</institution>, <addr-line>Bonn</addr-line>, <country>Germany</country></aff>
<aff id="aff3"><sup>3</sup><institution>Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, Laboratoire d&#x2019;Ecologie Alpine (LECA)</institution>, <addr-line>Grenoble</addr-line>, <country>France</country></aff>
<aff id="aff4"><sup>4</sup><institution>Department of Biology, Biotechnical Faculty, University of Ljubljana, Jamnikarjeva</institution>, <addr-line>Ljubljana</addr-line>, <country>Slovenia</country></aff>
<aff id="aff5"><sup>5</sup><institution>DivjaLabs Ltd., Alja&#x17e;eva ulica</institution>, <addr-line>Ljubljana</addr-line>, <country>Slovenia</country></aff>
<aff id="aff6"><sup>6</sup><institution>Institute of Nature Conservation, Polish Academy of Sciences, Adama Mickiewicza</institution>, <addr-line>Krak&#xf3;w</addr-line>, <country>Poland</country></aff>
<aff id="aff7"><sup>7</sup><institution>Departamento de Ciencias Integradas, Facultad de Ciencias Experimentales, Centro de Estudios Avanzados en F&#xed;sica, Matem&#xe1;ticas y Computaci&#xf3;n, Universidad de Huelva</institution>, <addr-line>Huelva</addr-line>, <country>Spain</country></aff>
<aff id="aff8"><sup>8</sup><institution>Wildlife Department, National Institute for Research and Development in Forestry &#x201c;Marin Dracea&#x201d;</institution>, <addr-line>Brasov</addr-line>, <country>Romania</country></aff>
<aff id="aff9"><sup>9</sup><institution>Department of Biology and Biotechnologies &#x201c;Charles Darwin&#x201d;, Sapienza University of Rome</institution>, <addr-line>Rome</addr-line>, <country>Italy</country></aff>
<aff id="aff10"><sup>10</sup><institution>Department of Natural Science and Environmental Health, University of South-Eastern Norway</institution>, <addr-line>B&#xf8;</addr-line>, <country>Norway</country></aff>
<aff id="aff11"><sup>11</sup><institution>Institute for Wildlife Biology and Game Management, University for Natural Resources and Life Sciences</institution>, <addr-line>Vienna</addr-line>, <country>Austria</country></aff>
<aff id="aff12"><sup>12</sup><institution>Department of Conservation Biology, University of Goettingen</institution>, <addr-line>G&#xf6;ttingen</addr-line>, <country>Germany</country></aff>
<aff id="aff13"><sup>13</sup><institution>Departamento de Biolog&#xed;a Vegetal y Ecolog&#xed;a, Facultad de Biolog&#xed;a, Universidad de Sevilla</institution>, <addr-line>Sevilla</addr-line>, <country>Spain</country></aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Janet Rachlow, University of Idaho, United States</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Erin Poor, The Nature Conservancy, United States; Aaron Adams, Bonefish and Tarpon Trust, United States</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Trishna Dutta, <email xlink:href="mailto:trishnad@gmail.com">trishnad@gmail.com</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>12</day>
<month>07</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>11</volume>
<elocation-id>1078649</elocation-id>
<history>
<date date-type="received">
<day>24</day>
<month>10</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>19</day>
<month>06</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Dutta, De Barba, Selva, Fedorca, Maiorano, Thuiller, Zedrosser, Signer, Pfl&#xfc;ger, Frank, Lucas and Balkenhol</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Dutta, De Barba, Selva, Fedorca, Maiorano, Thuiller, Zedrosser, Signer, Pfl&#xfc;ger, Frank, Lucas and Balkenhol</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<sec>
<title>Introduction</title>
<p>Connected landscapes can increase the effectiveness of protected areas by facilitating individual movement and gene flow between populations, thereby increasing the persistence of species even in fragmented habitats. Connectivity planning is often based on modeling connectivity for a limited number of species, i.e., &#x201c;connectivity umbrellas&#x201d;, which serve as surrogates for co-occurring species. Connectivity umbrellas are usually selected <italic>a priori</italic>, based on a few life history traits and often without evaluating other species.</p>
</sec>
<sec>
<title>Methods</title>
<p>We developed a quantitative method to identify connectivity umbrellas at multiple scales. We demonstrate the approach on the terrestrial large mammal community (24 species) in continental Europe at two scales: 13 geographic biomes and 36 ecoregions, and evaluate the interaction of landscape characteristics on the selection of connectivity umbrellas.</p>
</sec>
<sec>
<title>Results</title>
<p>We show that the number, identity, and attributes of connectivity umbrellas are sensitive to spatial scale and human influence on the landscape. Multiple species were selected as connectivity umbrellas in 92% of the geographic biomes (average of 4.15 species) and 83% of the ecoregions (average of 3.16 species). None of the 24 species evaluated is by itself an effective connectivity umbrella across its entire range. We identified significant interactions between species and landscape attributes. Species selected as connectivity umbrellas in regions with low human influence have higher mean body mass, larger home ranges, longer dispersal distances, smaller geographic ranges, occur at lower population densities, and are of higher conservation concern than connectivity umbrellas in more human-influenced regions. More species are required to meet connectivity targets in regions with high human influence (average of three species) in comparison to regions with low human influence (average of 1.67 species).</p>
</sec>
<sec>
<title>Discussion</title>
<p>We conclude that multiple species selected in relation to landscape scale and characteristics are essential to meet connectivity goals. Our approach enhances objectivity in selecting which and how many species are required for connectivity conservation and fosters well-informed decisions, that in turn benefit entire communities and ecosystems.</p>
</sec>
</abstract>
<kwd-group>
<kwd>connectivity</kwd>
<kwd>umbrella species</kwd>
<kwd>Europe</kwd>
<kwd>objectivity</kwd>
<kwd>surrogate species</kwd>
<kwd>landscape</kwd>
<kwd>human influence</kwd>
</kwd-group>
<contract-num rid="cn001">01LC1614A</contract-num>
<contract-sponsor id="cn001">BiodivERsA<named-content content-type="fundref-id">10.13039/100019184</named-content>
</contract-sponsor>
<counts>
<fig-count count="5"/>
<table-count count="0"/>
<equation-count count="2"/>
<ref-count count="46"/>
<page-count count="11"/>
<word-count count="5567"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Conservation and Restoration Ecology</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>    <p>Biodiversity is highly threatened due to increasing human influence on the planet. We are currently witnessing the so-called sixth mass extinction, primarily due to rapid habitat loss and fragmentation, which is further exacerbated by climate change (<xref ref-type="bibr" rid="B15">Ceballos et&#xa0;al., 2017</xref>). Habitat alteration restricts species to small populations in isolated patches with a high rate of ecosystem decay (<xref ref-type="bibr" rid="B16">Chase et&#xa0;al., 2020</xref>). This reduces movement and gene flow between populations (<xref ref-type="bibr" rid="B17">Crooks et&#xa0;al., 2017</xref>), and ultimately elevates the risk of species extinction (<xref ref-type="bibr" rid="B25">Gilpin and Soul&#xe9;, 1986</xref>). Halting biodiversity loss therefore requires urgent, concerted and sustained efforts to protect high-quality sites capable of sustaining viable populations, while simultaneously supporting connectivity among populations and movement of individuals across the landscape (<xref ref-type="bibr" rid="B6">Boyd et&#xa0;al., 2008</xref>).</p>
<p>Landscape connectivity buffers the effects of local extinction processes by facilitating the movement and effective dispersal of individuals, thus maintaining gene flow between populations, and supporting species range-shifts in response to changing climate and land-use regimes (<xref ref-type="bibr" rid="B36">Robillard et&#xa0;al., 2015</xref>). By increasing the potential for species dispersal into climatically suitable areas (<xref ref-type="bibr" rid="B4">Belote et&#xa0;al., 2017</xref>), ecological connectivity fosters resilient ecosystems that are effective in delivering ecosystem services (<xref ref-type="bibr" rid="B32">Mitchell et&#xa0;al., 2013</xref>). Consequently, connectivity has been recognized as an integral component of biodiversity conservation in international agreements, e.g. Aichi Target 11 of the United Nations Convention on Biological Diversity (<xref ref-type="bibr" rid="B12">CBD, 2010</xref>), the current draft of the post-2020 CBD framework (<xref ref-type="bibr" rid="B13">CBD, 2020</xref>), and in the European Union&#x2019;s Habitats Directive and Biodiversity Strategy for 2030 (<xref ref-type="bibr" rid="B21">European Commission, 2020</xref>).</p>
<p>Surrogate approaches are frequently used in conservation planning to compensate for incomplete ecological knowledge on all species in any ecosystem (<xref ref-type="bibr" rid="B43">Wiens et&#xa0;al., 2008</xref>). Two types of surrogates used in connectivity planning are species-based (or fine-filter) approaches where one or few species are used as surrogates for several co-occurring species (<xref ref-type="bibr" rid="B11">Caro and O&#x2019;Doherty, 1999</xref>), and habitat-based (or coarse-filter) approaches wherein landscape naturalness (<xref ref-type="bibr" rid="B39">Theobald et&#xa0;al., 2012</xref>) or land-facets (<xref ref-type="bibr" rid="B8">Brost and Beier, 2012</xref>) are used as surrogates for the species that inhabit them. Evaluations of coarse vs fine-filter approaches suggest significant trade-offs (<xref ref-type="bibr" rid="B27">Krosby et&#xa0;al., 2015</xref>), but overall the fine-filter approach has been found to be more effective (<xref ref-type="bibr" rid="B31">Meurant et&#xa0;al., 2018</xref>).</p>
<p>The species-based surrogate approach should aim to identify a suite of surrogate species, which we refer to as a &#x201c;connectivity umbrella&#x201d; in this paper, that encompass areas most likely to be used by several co-occurring species. Connectivity umbrellas are often chosen based on a few traits such as body mass or home range size, their charisma, conservation status (<xref ref-type="bibr" rid="B3">Beier et&#xa0;al., 2009</xref>; <xref ref-type="bibr" rid="B40">Wang et&#xa0;al., 2018</xref>), or because they have previously been identified for other purposes such as flagships for fund-raising, indicators of ecosystem health, or umbrellas for habitat management (<xref ref-type="bibr" rid="B10">Caro and Girling, 2010</xref>). Pre-selected connectivity umbrellas based on few criteria such as charismatic large species, have been shown to be ineffective at representing other species in multiple studies (<xref ref-type="bibr" rid="B3">Beier et&#xa0;al., 2009</xref>; <xref ref-type="bibr" rid="B18">Cushman and Landguth, 2012</xref>; <xref ref-type="bibr" rid="B40">Wang et&#xa0;al., 2018</xref>).</p>
<p>In addition to species-based surrogates being better performing than habitat based surrogates, research suggests that multiple species are better than any single surrogate species (<xref ref-type="bibr" rid="B31">Meurant et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B40">Wang et&#xa0;al., 2018</xref>), surrogate species are effective within their own taxonomic group (<xref ref-type="bibr" rid="B7">Brodie et&#xa0;al., 2015</xref>), and a diverse set of species reflects the needs of other co-occurring species (<xref ref-type="bibr" rid="B18">Cushman and Landguth, 2012</xref>). Despite the emphasis on multi-species connectivity models (<xref ref-type="bibr" rid="B45">Wood et&#xa0;al., 2022</xref>), there are only a few methods that use data-driven approaches to identify which and how many species could represent the connectivity needs of all species in any given region.</p>
<p>In this study, we address this research gap in connectivity science. Our goal is to develop an approach that increases the objectivity in selecting connectivity umbrella species from a pool of candidate species. We demonstrate the approach at multiple scales and evaluate the interaction of landscape characteristics and connectivity umbrella selection using the terrestrial large-mammal community in continental Europe as a case study.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>Methods</title>
<p>For the large mammal community in Europe, we modelled potential connectivity between protected areas (PAs) and created a connectivity target representing high-probability connectivity areas for multiple species present in the region. We then ranked each species in the region using a cumulative measure that indicates how well a species represents the target and connectivity network of co-occurring species. Starting with the highest-ranked species, we sequentially added the connectivity network of all species, and plotted the contribution, i.e., percent increase of the connectivity target for each addition (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1</bold></xref>). The final connectivity umbrella comprises the combination of species that are required to reach a certain coverage of the target (set at 95% in this study), as adding more species contributed only marginally towards reaching the connectivity target. All steps of the analyses are presented in greater detail in <xref ref-type="supplementary-material" rid="SM1"><bold>Figure S1</bold></xref>. All quantitative analysis was performed using R statistical software (<uri xlink:href="https://www.r-project.org/">https://www.r-project.org/</uri>) and all maps were created in ArcGIS 10.3.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Conceptual figure of the selection process for connectivity umbrellas. The oval represents the region of analysis; the dark gray polygon represents the connectivity target. Five hypothetical species are shown in different colors and arranged according to their overall rank. Species are then added sequentially while quantifying their contribution towards the connectivity target <bold>(A&#x2013;C)</bold> until an asymptote is reached. In this example, four species (blue, green, orange, maroon) are required to cover 95% of the target [dashed line in <bold>(C)</bold>].</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-11-1078649-g001.tif"/>
</fig>
<p>The four main steps involved in this approach are: (1) Modelling connectivity for each species, (2) Generating connectivity targets per region, (3) Ranking species in each region, and finally (4) Determining the suite of umbrella species.</p>
<sec id="s2_1">
<label>2.1</label>
<title>Step 1. Modelling connectivity for each species</title>
<p>We modelled connectivity between Protected Areas (PAs) within the suitable habitat for all large terrestrial mammals in Europe. We first selected species for the analysis, identified PAs to be used as nodes for connectivity analysis for each species, generated species-specific resistance surfaces that represented the potential difficulty of movement between nodes, and estimated potential connectivity for each species (<xref ref-type="supplementary-material" rid="SM1"><bold>Figure S1</bold></xref>, Step 1).</p>
<sec id="s2_1_1">
<label>2.1.1</label>
<title>Step 1a. Selecting species and traits for the study</title>
<p>We selected large terrestrial mammals, species with a mean body mass larger than 3 kg (<xref ref-type="bibr" rid="B9">Cardillo et&#xa0;al., 2005</xref>), that are native to Europe. We used the PanTHERIA database to obtain data on species traits (<xref ref-type="bibr" rid="B26">Jones et&#xa0;al., 2009</xref>) and putative maximum dispersal distances based on allometric scaling equations from previous studies (<xref ref-type="bibr" rid="B37">Santini et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B42">Whitmee and Orme, 2013</xref>). Details of the species used in this study are presented in <xref ref-type="supplementary-material" rid="SM1"><bold>Table S1</bold></xref>.</p>
</sec>
<sec id="s2_1_2">
<label>2.1.2</label>
<title>Step 1b. Selecting nodes</title>
<p>Within the extent of occurrence and suitable habitat of each species (<xref ref-type="bibr" rid="B29">Maiorano et&#xa0;al., 2013</xref>), we selected PAs (IUCN Level I&#x2013;IV) and Natura 2000 sites that reported presence of the species (<xref ref-type="bibr" rid="B23">European Environment Agency, 2020</xref>). We calculated the Euclidean distance between each pair of nodes and only retained pairs that were equal to or less than the putative maximum dispersal distance presented in <xref ref-type="supplementary-material" rid="SM1"><bold>Table S1</bold></xref>. As a result, each species had a unique set of pairwise nodes.</p>
</sec>
<sec id="s2_1_3">
<label>2.1.3</label>
<title>Step 1c. Generating the resistance surface</title>
<p>We used previously published habitat suitability maps of mammals in Europe (<xref ref-type="bibr" rid="B29">Maiorano et&#xa0;al., 2013</xref>), which were generated from three environmental variables &#x2013; land cover, elevation, and distance to water. The habitat suitability models contained pixel values of 0, 1, 2 corresponding to non-habitat, marginal habitat, suitable habitat, respectively, for each species. Because animals are more tolerant of suboptimal habitat during dispersal (<xref ref-type="bibr" rid="B30">Mateo-S&#xe1;nchez et&#xa0;al., 2015</xref>), we coded the suitability values as 0 (unsuitable) and 1 (suitable, where we pooled suitable and marginal habitats) and resampled each raster using a 3&#xd7;3 moving window. For each of the 24 species, this step produced a continuous habitat suitability raster at 1&#xa0;km resolution. We then converted this suitability surface to a landscape resistance raster, which represents the ability of an organism to cross a particular environment (<xref ref-type="bibr" rid="B46">Zeller et&#xa0;al., 2012</xref>).</p>
<p>Previous studies have found that resistance is often a non-linear negative exponential function of habitat suitability (<xref ref-type="bibr" rid="B20">Dudaniec et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B30">Mateo-S&#xe1;nchez et&#xa0;al., 2015</xref>). Therefore, we converted habitat suitability surfaces using a negative exponential function (<xref ref-type="bibr" rid="B2">Balkenhol et&#xa0;al., 2020</xref>) derived from the following formula:</p>
<disp-formula>
<mml:math display="block" id="M2">
<mml:mrow>
<mml:mtext>R</mml:mtext>
<mml:mo>=</mml:mo>
<mml:mtext>CS</mml:mtext>
<mml:mo>&#xd7;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mn>100</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>HS)</mml:mtext>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where R is the resistance assigned to a specific cell, CS is the cell size in meters (here 1000&#xa0;m) and HS is the habitat suitability associated with that cell. This conversion assumes that animals are more tolerant of suboptimal habitat and respond less strongly to low habitat suitability during dispersal. The final resistance values ranged from 1000&#x2013;10,000 resistance units.</p>
</sec>
<sec id="s2_1_4">
<label>2.1.4</label>
<title>Step 1d: Corridor simulation</title>
<p>We mapped potential connectivity for each species using the stochastic modelling program LSCorridors (<xref ref-type="bibr" rid="B35">Ribeiro et&#xa0;al., 2017</xref>). LSCorridors uses the resistance surface to produce potential corridors between each pair of nodes by simulating multiple dispersal events. It then adds all pairwise maps to produce the overall connectivity map for the species, quantified as the Route Selection Frequency Index (RSFI), a measure of the probability that a particular pixel will be used by any dispersing individual. We conducted 50 simulations between each pair of nodes, without landscape influence (Measures by Pixel, MP method, scale = 1000m, variability = 2). With these settings, stochasticity between simulations is introduced by different starting and ending points within the nodes for each simulation and by the variability factor in the resistance map, which represents uncertainty in the values used to produce the resistance surface. We estimated potential connectivity between a total of 53023 pairs of nodes, which amounts to a total of 2,651,150 pairwise simulations.</p>
<p>At the end of all simulations for each species, we further pruned the pairwise simulations to retain only those PA-pairs that had a mean corridor length equal to or less than the maximum putative dispersal distance. Since the maximum RSFI values were different for each of the species, we normalized the final connectivity maps of the 24 mammals for further analyses. We conducted further analyses at two spatial scales as described below.</p>
</sec>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Step 2: Generating connectivity targets per region</title>
<p>Defining regions: We conducted analyses at two hierarchical levels based on ecoregion classifications within broad geographical regions in continental Europe. We combined the six biomes present in Europe (<xref ref-type="bibr" rid="B34">Olson et&#xa0;al., 2001</xref>) with the six broad geographic regions (excluding islands) which produced 13 geographic biomes, henceforth referred simply as biomes in this study. Nested within these 13 biomes there are a total of 36 ecoregions (eight biomes have 31 nested ecoregions while five biomes do not have ecoregions) (<xref ref-type="supplementary-material" rid="SM1"><bold>Figure S2</bold></xref>).</p>
<p>Connectivity target: Our approach requires users to pre-define a spatially explicit connectivity target. We tested the impact of different connectivity targets based on three different conservation goals. The three targets are aimed at protecting connectivity networks for all species in the region (Target EQ), only high-quality connectivity networks (Target Q), or connectivity networks of conservation-concern species only (Target T) (<xref ref-type="supplementary-material" rid="SM1"><bold>Figure S1</bold></xref>, Step 2). Details of how these targets were derived are explained here:</p>
<list list-type="bullet">
<list-item>
<p>Target EQ: This target represents the best quality corridor cells for all species equally. It is produced by first extracting high-quality corridor cells (top 25 percentile values for each species), and then summing the high-quality corridor cells for all species.</p>
</list-item>
<list-item>
<p>Target Q: This target represents only the best quality corridor cells irrespective of which species are represented in the target. It is produced by first adding all the normalized species corridor maps into one composite raster, and then extracting the cells that contain the top 25 percentile values.</p>
</list-item>
<list-item>
<p>Target T: This target represents all corridor cells used by species of conservation concern, which we defined as species endemic to Europe or with an IUCN red list status of Near Threatened or higher (<xref ref-type="supplementary-material" rid="SM1"><bold>Table S1</bold></xref>). In Europe, these species are <italic>Lynx pardinus, Bison bonasus, Gulo gulo, Marmota marmota, Rupicapra pyrenaica, Capra ibex</italic>, and <italic>Capra pyrenaica.</italic> This target is produced by adding corridor cells of all conservation-concern species in any given region. Because nine of the 13 biomes contain species of conservation concern, this target existed only in nine biomes.</p>
</list-item>
</list>
<p>We chose target EQ to demonstrate our data-driven approach to select connectivity umbrella species as we believe most conservation actions would be concerned about several species in the region, rather than just high-quality regions (target Q) or only species of conservation concern (target T). We compared connectivity umbrellas selected for target EQ with target Q and T to test the sensitivity of our approach and results are presented in <xref ref-type="supplementary-material" rid="SM1"><bold>Figure S3</bold></xref>.</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Step 3: Ranking species</title>
<p>To rank the species in each region, we calculated an umbrella score &#x201c;<italic>Uir&#x201d;</italic>, which quantifies the suitability of each species as a connectivity umbrella. For each region, we used the connectivity maps of each species, wherein the normalized RSFI value represents the probability of connectivity of each cell for that species and a set of five criteria using the calculation:</p>
<disp-formula>
<mml:math display="block" id="M1">
<mml:mrow>
<mml:mi>U</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mo>&#x223c;</mml:mo>
<mml:mi>N</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>A</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>F</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>S</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:math>
</disp-formula>
<p><italic>U</italic> is the overall umbrella score for species <italic>i</italic> in region <italic>r</italic>, calculated as the product of five criteria listed below:</p>
<p><italic>N</italic> is the total number of species present in the region <italic>r</italic> whose corridors overlap with that of species <italic>i</italic>, calculated as the total number of overlapping species (<xref ref-type="supplementary-material" rid="SM1"><bold>Figure S1</bold></xref>, Step 3a).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Connectivity umbrella species selected to reach 95% of the connectivity target in the thirteen biomes within major geographic regions in Europe. Species of conservation concern in Europe-wide assessments are represented by black silhouettes, all other species are represented by gray silhouettes.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-11-1078649-g002.tif"/>
</fig>
<p><italic>A</italic> is the mean area of overlapping corridors of species <italic>i</italic> with all species corridors present in region <italic>r</italic>, calculated as the mean area (km<sup>2</sup>) of pairwise overlap of species <italic>i</italic> with all other species (<xref ref-type="supplementary-material" rid="SM1"><bold>Figure S1</bold></xref>, Step 3b).</p>
<p><italic>F</italic> is the fraction of the connectivity target in region <italic>r</italic> covered by species <italic>i</italic>, calculated as the percent area of intersection between the connectivity target and the corridor network of species <italic>i</italic> (<xref ref-type="supplementary-material" rid="SM1"><bold>Figure S1</bold></xref>, Step 3c).</p>
<p><italic>S</italic> is the mean strength of pairwise overlapping corridors, calculated as the mean of Bhattacharyya&#x2019;s affinity index (<xref ref-type="bibr" rid="B5">Bhattacharya, 1943</xref>) that measures volumetric intersection of each pair of species (<xref ref-type="bibr" rid="B24">Fieberg and Kochanny, 2005</xref>). Bhattacharyya&#x2019;s affinity index is a statistical measure of affinity between two populations that assumes they use space independently of one another. Values range from zero (no overlap) to 1 (complete overlap) (<xref ref-type="supplementary-material" rid="SM1"><bold>Figures S1</bold></xref>, Step 3d).</p>
<p><italic>P</italic> is the proportion of the total connectivity network of species <italic>i</italic> that lies within the connectivity target of region r, calculated as the sum of all pixel values within region <italic>r</italic> (<xref ref-type="supplementary-material" rid="SM1"><bold>Figure S1</bold></xref>, Step 3e).</p>
<p>Criteria <italic>N</italic> (number of overlapping species)<italic>, A</italic> (overlapping area with other species)<italic>, F</italic> (fraction of target covered by the species) characterize the numerical suitability of a species at representing other species and the target. Criteria <italic>S</italic> (strength of overlapping corridors) and <italic>P</italic> (proportion of the species corridors within target) characterize the quality of a species at representing co-occurring species corridors and the target. Combining these criteria effectively select species that are representative of co-occurring species and the region under consideration (<xref ref-type="supplementary-material" rid="SM1"><bold>Figure S1</bold></xref>, Step 3e).</p>
</sec>
<sec id="s2_4">
<label>2.4</label>
<title>Step 4: Determining the suite of umbrella species</title>
<p>Within each region, we ranked each species by its <italic>Uir</italic> value so that the species with the highest value received the highest rank. Effectively, this means that species with the highest ranks are those that overlapped with many other species, and did so over a larger area, with the highest strength of overlap, represented a large part of the connectivity network, and had a large part of their connectivity network within the region of interest. We plotted the species by their rank vs the additional contribution they made to the connectivity target. We first plotted the proportion of the connectivity target achieved by adding the corridors of the highest ranked species. We continued adding species according to their ranks, calculating the additional coverage of the connectivity target (that was not already covered by the previous species) for each species addition. We defined the suite of umbrella species as the number and combination of species required to cover at least 95% area of the connectivity target as adding more species beyond this asymptote did not make much improvement towards achieving the connectivity target (<xref ref-type="supplementary-material" rid="SM1"><bold>Figure S1</bold></xref>, Step 4).</p>
</sec>
<sec id="s2_5">
<label>2.5</label>
<title>Evaluating the role of underlying landscape</title>
<p>Finally, to evaluate the role of landscape characteristics on the selection of connectivity umbrellas, we conducted a cluster analysis of all 36 ecoregions based on their mean human footprint, mean landscape fragmentation and the mean latitude. The human footprint index (Wildlife Conservation Society &#x2013; WCS and Center for International Earth Science Information Network &#x2013; CIESIN, Columbia University, 2005), generated from human population pressure, land use and infrastructure, and access, represents the cumulative anthropogenic impact on the environment, whereas the landscape fragmentation index (<xref ref-type="bibr" rid="B22">European Environment Agency, 2016</xref>) measured as the effective mesh density, i.e. landscape patches per 1000 km<sup>2</sup>, represents a measure of the degree to which the landscape is fragmented. We performed agglomerative hierarchical cluster analysis on ecoregions to identify distinct clusters in our dataset. We used Ward&#x2019;s method of agglomeration, which produces clusters of more equal size by keeping distances within the clusters as small as possible (<xref ref-type="bibr" rid="B41">Ward, 1963</xref>). We identified five unique groups using a dendrogram on the cluster solutions that we had obtained using the hclust function available with R library fastcluster (<xref ref-type="bibr" rid="B33">M&#xfc;llner, 2013</xref>).</p>
<p>Within the clusters, we evaluated six attributes of connectivity umbrella species (mean body mass, home range size, population density, putative maximum dispersal distance, geographic range in Europe and conservation status (<xref ref-type="supplementary-material" rid="SM1"><bold>Table S1</bold></xref>). We conducted Kruskal-Wallis tests to identify significant differences of attributes of connectivity umbrella species in the different clusters and used a <italic>post-hoc</italic> Dunns test to identify which clusters were significantly different using Holms-adjusted p-values to account for multiple comparisons.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<p>Using the approach presented in <xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1</bold></xref>, we demonstrate that a systematic and objective assessment of which and how many species are required to meet connectivity needs of multiple co-occurring species is indeed possible. In this process, we discovered some important patterns.</p>
<p>Multiple species are required to meet connectivity goals at both scales considered. Ninety-two percent of the biomes (<xref ref-type="fig" rid="f2"><bold>Figure&#xa0;2</bold></xref>) and 83% of the ecoregions (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3</bold></xref>) had more than one connectivity umbrella species. Both the number and identity of species are scale dependent. More species are selected at the biome level (average 4.15, range 1&#x2013;9 species; <xref ref-type="fig" rid="f2"><bold>Figure&#xa0;2</bold></xref>) than at the ecoregion level (average 3.16, range 1&#x2013;9 species; <xref ref-type="fig" rid="f3"><bold>Figures&#xa0;3A</bold></xref>, <xref ref-type="supplementary-material" rid="SM1"><bold>S4</bold></xref>). The identity of the species selected also varies across scales. For example, six species are selected in the temperate broadleaf mixed forest biome (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3A</bold></xref>), whereas several of the species selected in the eleven nested ecoregions (ranging from 1 to 6) are exclusive to the ecoregions, i.e., not selected in the biome-level analysis (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3B</bold></xref>).</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>The number <bold>(A)</bold> and identity <bold>(B)</bold> of connectivity umbrella species vary across spatial scales. More species are selected at the broad biome level [represented by large dots in <bold>(A)</bold>] than at the finer ecoregion level [represented by smaller dots in <bold>(A)</bold>]. Biomes are represented by different colors. Connectivity umbrellas selected at finer scale often include unique and conservation-concern species [grey boxes in <bold>(B)</bold>], not selected at the coarse scale white box in <bold>(B)</bold>. Species identified to be of conservation concern in ecoregional assessments are represented by black silhouettes.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-11-1078649-g003.tif"/>
</fig>
<p>Connectivity coverage achieved using single species is highly variable (<xref ref-type="supplementary-material" rid="SM1"><bold>Figure S5</bold></xref>). For example, if only the highest-ranking species were to be selected, an average of 75% of the regional targets would be covered, ranging from ~14% by the Iberian lynx (<italic>Lynx pardinus</italic>) in Southwest Iberian Mediterranean sclerophyllous and mixed forests to ~97% by the wolf (<italic>Canis lupus</italic>) in South Apennine mixed montane forests. If only the highest-ranking species of conservation concern is selected, an average of 25% of the individual regional targets would be covered, ranging from ca. 1% by the bison in the Carpathian Mountains to 95% by the wolverine in the Scandinavian Tundra. Some ecoregions in the Scandinavian and Apennine peninsulas are exceptions to this general pattern.</p>
<p>Our results indicate that even when species of conservation concern are present in a region, they may not necessarily be the most suitable connectivity umbrellas. A total of seven conservation-concern species are present in 17 ecoregions (47% of ecoregions, average 1.3 conservation-concern species/region), but are selected as connectivity umbrellas in only 9 ecoregions (~53% where they are present). The number or proportion of conservation-concern species present in any region is not significantly correlated to the number or proportion of conservation-concern species in the suite of connectivity umbrella species (Pearson&#x2019;s r = 0.04 and 0.26, respectively; p &gt; 0.05 in both cases). We used Red List status at the European level because it was consistent with the European extent of our study, but threat assessments at regional scales may also be relevant. For example, none of the species selected as connectivity umbrellas in the temperate broadleaf mixed forest biome are of conservation concern at the European scale (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3B</bold></xref>), but the brown bear (<italic>Ursus arctos</italic>) and Eurasian lynx (<italic>Lynx lynx</italic>), selected at the ecoregion scale are categorized as vulnerable and endangered in the Rhodope and Dinaric Mountain ecoregions (<xref ref-type="bibr" rid="B38">Temple and Terry, 2007</xref>).</p>
<p>Of the 24 species evaluated, not even one species is an effective connectivity umbrella across its entire range (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4</bold></xref>). Even the species selected most often, the wolverine (<italic>Gulo gulo</italic>) and moose (<italic>Alces alces</italic>), were only selected as connectivity umbrellas in 60% of the regions in which they occur. While there is a significantly higher chance of being selected as an umbrella species for wider ranging species, i.e., species that occur in many regions (Pearson&#x2019;s r = 0.73, p&lt;.05), this effect is lost when comparing the occurrence of species across multiple ecoregions with the proportion of the times that it is selected as a connectivity umbrella (Pearson&#x2019;s r = 0.13, p = 0.66).</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>No species is selected as a connectivity umbrella everywhere it occurs. For each species in the x-axis, the y-axis represents the number of times it was selected as a connectivity umbrella everywhere it is present. For example, the distribution ranges of the moose (<italic>Alces alces</italic>), wildcat (<italic>Felis silvestris</italic>), and chamois (<italic>Rupicapra rupicapra</italic>) are shown in pink, and the ecoregions where they are selected as connectivity umbrellas are shown in hashed lines.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-11-1078649-g004.tif"/>
</fig>
<p>We identified five unique clusters of ecoregions that ranged from very low to very high human influence (<xref ref-type="supplementary-material" rid="SM1"><bold>Figure S6</bold></xref>). Mountain ranges and regions at higher latitudes are generally less fragmented and impacted by humans than regions in the plains or at lower latitudes. The maximum dispersal distance (Kruskal-Wallis chi-squared = 10.942, df = 4, p-value = 0.027), home range (chi-squared = 12.401, df = 4, p-value = 0.014), proportion of range in Europe (chi-squared = 21.014, df = 4, p-value = 0.0003), and conservation status (chi-squared = 9.9392, df = 4, p-value = 0.041) of connectivity umbrella species are significantly different between the clusters. On average, species selected as connectivity umbrellas in regions with low human influence have higher mean body mass, larger home ranges, longer dispersal distances, smaller geographic ranges, occur at lower population densities, and are of higher conservation concern than connectivity umbrellas in more human-influenced regions (<xref ref-type="fig" rid="f5"><bold>Figures&#xa0;5A&#x2013;F</bold></xref>). Proportionally more species are required to meet connectivity targets in regions with high human influence (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5G</bold></xref>).</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Attributes of species across clusters arranged from lowest (1) to highest (5) human impact. The means and standard errors of attributes of connectivity umbrellas (dark blue) and non-umbrellas (light blue) are shown in <bold>(A&#x2013;F)</bold>, along with the mean attribute value in the cluster as a grey bar. The proportion of species selected as connectivity umbrellas in each cluster is shown in <bold>(G)</bold>. Significant difference of mean values between clusters are shown on the top in solid lines and between umbrella and non-umbrella species within clusters are shown at the bottom in panels <bold>(A&#x2013;F)</bold>. * indicates a p-value &#x2264; 0.05, ** a p-value &#x2264; 0.01, and *** a p-value &#x2264; 0.001.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-11-1078649-g005.tif"/>
</fig>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<p>Our approach increases objectivity when selecting umbrella species for connectivity conservation. It can be applied to multiple-taxonomic groups and to varying spatial scales. Unlike static target-setting in conservation prioritization exercises, our approach is flexible and can be adapted to test which species would suit different conservation goals and scenarios.</p>
<p>Two methods have been developed to specifically select species-based surrogates for connectivity planning. <xref ref-type="bibr" rid="B28">Lechner et&#xa0;al. (2017)</xref> use a cluster analysis on habitat requirements, connectivity elements, and gap-crossing distance thresholds to group twelve candidate mammal species into five distinct dispersal guilds. Connectivity for each guild is then modelled separately to identify which guild of species are more sensitive to habitat fragmentation. <xref ref-type="bibr" rid="B31">Meurant et&#xa0;al. (2018)</xref> present an approach to select multi-taxa surrogates from a pool of candidate species that represent a diversity of habitat needs and movement abilities. From a set of fourteen candidate species, the authors create a reference spatial prioritization map and run multiple permutations and combinations of the fourteen candidate species to identify the species subset that best retains different thresholds of the reference spatial priority map. Our approach is different from both methods. First, we include the entire terrestrial large mammal community in Europe without going through a pre-selection process as in <xref ref-type="bibr" rid="B31">Meurant et&#xa0;al. (2018)</xref>, where they selected fourteen candidate species from a pool of 48 mammals, 216 birds, and 32 amphibians and reptiles based on a previous study (<xref ref-type="bibr" rid="B1">Albert et&#xa0;al., 2017</xref>). Second, we do not group species into dispersal guilds, diversity of habitat, or movement requirements. We treat each species independently. Further, instead of using a conservation concern filter in selecting the species to include in the candidate pool, we use the results to guide and test how frequently species of conservation concern are suitable connectivity umbrellas. Finally, in addition to addressing which and how many species are required for connectivity planning, we also evaluate the species-surrogates approach across two different geographic scales and assess the characteristics of connectivity umbrellas in landscapes with varying levels of human influence.</p>
<p>Our results re-confirm previous findings on the fine-filter surrogate approach that multiple species are better than any single surrogate species (<xref ref-type="bibr" rid="B31">Meurant et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B40">Wang et&#xa0;al., 2018</xref>) and <italic>a priori</italic> selection of certain groups, e.g., carnivores, may not be suitable connectivity umbrellas (<xref ref-type="bibr" rid="B3">Beier et&#xa0;al., 2009</xref>). Most suites of species selected in our analyses consist of a combination of herbivores and carnivores. Such mixed suites of species have diverse ecological attributes and dispersal abilities and are therefore more likely to accommodate the ecological requirements of diverse terrestrial mammal communities and their differing sensitivities to environmental conditions. In addition to corroborating previous findings, we show that both the number and the identity of species differ by the scale of analysis, and that species selected at finer scales are not subsets of those selected at broader scales (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3</bold></xref>).</p>
<p>Our approach is flexible in two ways. First, the result is a suite of species that are needed to cover 95% of the connectivity target. Depending on the regional context, it is possible to sub-select from within this suite of species, especially when species provide very little percent improvement towards the overall target. Second, the evaluation of sensitivity to different connectivity targets (<xref ref-type="supplementary-material" rid="SM1"><bold>Figure S3</bold></xref>) provides a way to select species that are suitable for a diversity of conservation goals. For example, moose and wolverine in the northern regions, red deer (<italic>Cervus elaphus</italic>) and Iberian Ibex (<italic>Capra pyrenaica</italic>) in the Mediterranean forests, woodland and scrub of the Iberian Peninsula and wolf in the Apennine peninsula, red deer and red fox (<italic>Vulpes vulpes</italic>) in temperate broadleaf and mixed forest in central Europe, and the Alpine ibex in the Alps, are species selected in all three connectivity targets which were created with very different conservation goals. In general, target EQ required more species than target Q or T. This is intuitive because target EQ represents all species, target Q represents only high-quality corridors irrespective of species, and target T only represents species of conservation concern. The identity of species selected across the targets changed in their rank, but in general the top-ranking species were consistent across the different targets, and the pool of species that was most highly ranked remained similar. This redundancy in the selection of species across disparate conservation goals can be a valuable tool to incorporate more certainty about including a particular species to be an umbrella for connectivity conservation.</p>
<p>An important finding is the interaction between the attributes of landscape and species selected as connectivity umbrellas. The interaction between species and landscape indicates that species selected to be surrogates for connectivity conservation need to be tailored to the level of human influence in the landscape under consideration. To be effective in highly humanized landscapes, numerically more species with a certain set of attributes (widely distributed but smaller-bodied species that occur at higher population densities and disperse over shorter distances) may be required. These attributes are generally not associated with umbrella species but may be the key to developing region-specific effective connectivity plans.</p>
<p>The aim of our study is not to prescribe which species should be used for connectivity conservation in Europe. Rather, we use continental Europe as an example to demonstrate an approach that is quite flexible and can accommodate many different conservation and modelling choices. Due to the continental scale of analyses, we use several simplifications in our modelling to optimize computational time. We used coarse-scale (1 km) existing habitat suitability maps to generate resistance surfaces and model connectivity between PAs, but it is noteworthy that source populations may also be present in habitat patches outside of the PA network, and individual species are likely to be impacted by more and different sets of variables. Therefore, we caution the use of the species presented here to be effective and &#x201c;correct&#x201d; connectivity umbrellas. The species presented here can be starting points, but we recommend users to create resistance surfaces using variables and resolutions that are relevant to their region of interest. It will most likely be computationally plausible to run regional models using non-protected habitat patches as source populations at a finer resolution and include more variables to develop more realistic resistance surfaces.</p>
<p>Our approach can easily be extended to other systems not presented here. For example, we analyze only terrestrial mammals, but other taxonomic groups of interest can easily be used to conduct similar analyses. We define three different connectivity targets based on different conservation goals, but connectivity targets may be defined in many other ways. Finally, we use a threshold of 95% of the target as our goal (guided by our observation that an asymptote is reached around 95%), but users may decide to use a different threshold. Our approach should be applicable in any focal landscape where it is possible to generate a connectivity network for multiple species.</p>
<p>Overall, we demonstrate that an objective selection of surrogate species for connectivity planning is possible. Considering multiple species, not just in relation to species attributes but in relation to landscape scale and characteristics may be the key towards identifying connectivity umbrellas that better represent the needs of entire species communities. In practice, decision-making in conservation is a complex process, involving many stakeholders, and is highly contextual. When species are tied in rank and contribution to the connectivity target, managers or conservationists could select species with a higher social acceptance or other issues that are relevant locally or regionally. Rather than advocating for one species over another, we present an approach to remove the guesswork over which species should be selected as connectivity umbrellas and increase objectivity in the selection of multispecies conservation strategies.</p>
<p>Conservation targets, whether modest at 30% of terrestrial land (<xref ref-type="bibr" rid="B14">CBD, 2023</xref>), or ambitious at 50% (<xref ref-type="bibr" rid="B19">Dinerstein et&#xa0;al., 2017</xref>), are unachievable without integrating connectivity conservation, especially in rapidly changing ecosystems in the Anthropocene. The key to selecting umbrella species for connectivity conservation is to move beyond arbitrary selection of species, for the selection process to be objective, and in relation to the scale and attributes of the landscape. Only then can the overall goal of connectivity conservation be realized and be effective in restoring ecosystems and decreasing biodiversity loss.</p>
</sec>
<sec id="s5" sec-type="data-availability">
<title>Data availability statement</title>
<p>The original contributions presented in the study are publicly available. This data can be found here: <ext-link ext-link-type="uri" xlink:href="https://figshare.com/s/23831fdc52a9dcc60a36">https://figshare.com/s/23831fdc52a9dcc60a36</ext-link>.</p>
</sec>
<sec id="s6" sec-type="author-contributions">
<title>Author contributions</title>
<p>TD, NB, and MDB contributed to conception and design of the study. JS contributed towards the analyses and SF&#x2019;s contributions helped improve the figures. TD wrote the first draft of the manuscript. All authors contributed to the article and approved the submitted version.</p>
</sec>
</body>
<back>
<sec id="s7" sec-type="funding-information">
<title>Funding</title>
<p>This research is part of project BearConnect, funded through the 2015&#x2013;2016 BiodivERsA COFUND call, with national funders Agence Nationale de la Recherche (ANR), France (Grant number: ANR-16-EBI3-0003), National Science Center (NCN), Poland (Grant number: 2016/22/Z/NZ8/00121), Federal Ministry of Education and Research (BMBF) and DLR Project Management Agency (DLR-PT), Germany (Grant number: 01LC1614A), Romanian National Authority for Scientific Research and Innovation, CCCDI &#x2013; UEFISCDI, Romania (Grant number: BiodivERsA3-2015-147-BearConnect (96/2016) within PNCDI III, and Norwegian Research Council (RCN), Norway (Grant number: 269863).</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>We thank Phylopic and all the creators who provided animal silhouettes in the public domain which we used in the different figures in this paper. We would like to thank Drs. Fanie Pelletier, Moreno Di Marco, Thomas Mueller, Ninon Meyer, and Nina Gerber for their reviews on previous versions of this manuscript. We acknowledge support by the Open Access Publication Funds of the University of Goettingen.</p>
</ack>
<sec id="s8" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>Author MDB is a founder of DivjaLabs Ltd.</p>
<p>The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s9" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s10" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fevo.2023.1078649/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fevo.2023.1078649/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="DataSheet_1.pdf" id="SM1" mimetype="application/pdf"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Albert</surname> <given-names>C. H.</given-names>
</name>
<name>
<surname>Rayfield</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Dumitru</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Gonzalez</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Applying network theory to prioritize multi-species habitat networks that are robust to climate and land-use change</article-title>. <source>Conserv. Biol.</source> <volume>31</volume>, <fpage>1383</fpage>&#x2013;<lpage>1396</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/cobi.12943</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Balkenhol</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Schwartz</surname> <given-names>M. K.</given-names>
</name>
<name>
<surname>Inman</surname> <given-names>R. M.</given-names>
</name>
<name>
<surname>Copeland</surname> <given-names>J. P.</given-names>
</name>
<name>
<surname>Squires</surname> <given-names>J. S.</given-names>
</name>
<name>
<surname>Anderson</surname> <given-names>N. J.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Landscape genetics of wolverines <italic>(Gulo gulo)</italic>: scale-dependent effects of bioclimatic, topographic, and anthropogenic variables</article-title>. <source>J. Mammal.</source> <volume>101</volume>, <fpage>790</fpage>&#x2013;<lpage>803</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/jmammal/gyaa037</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Beier</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Majka</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Newell</surname> <given-names>S. L.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Uncertainty analysis of least-cost modeling for designing wildlife linkages</article-title>. <source>Ecol. Appl.</source> <volume>19</volume>, <fpage>2067</fpage>&#x2013;<lpage>2077</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1890/08-1898.1</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Belote</surname> <given-names>R. T.</given-names>
</name>
<name>
<surname>Dietz</surname> <given-names>M. S.</given-names>
</name>
<name>
<surname>Jenkins</surname> <given-names>C. N.</given-names>
</name>
<name>
<surname>McKinley</surname> <given-names>P. S.</given-names>
</name>
<name>
<surname>Irwin</surname> <given-names>G. H.</given-names>
</name>
<name>
<surname>Fullman</surname> <given-names>T. J.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>). <article-title>Wild, connected, and diverse: building a more resilient system of protected areas</article-title>. <source>Ecol. Appl.</source> <volume>27</volume>, <fpage>1050</fpage>&#x2013;<lpage>1056</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/eap.1527</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bhattacharya</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>1943</year>). <article-title>On a measure of divergence between two statistical populations defined by their probability distributions</article-title>. <source>Bull. Calcutta Math. Soc</source> <volume>35</volume>, <fpage>99</fpage>&#x2013;<lpage>109</lpage>.</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Boyd</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Brooks</surname> <given-names>T. M.</given-names>
</name>
<name>
<surname>Butchart</surname> <given-names>S. H. M.</given-names>
</name>
<name>
<surname>Edgar</surname> <given-names>G. J.</given-names>
</name>
<name>
<surname>Fonseca</surname> <given-names>G. A. B. D.</given-names>
</name>
<name>
<surname>Hawkins</surname> <given-names>F.</given-names>
</name>
<etal/>
</person-group>. (<year>2008</year>). <article-title>Spatial scale and the conservation of threatened species</article-title>. <source>Conserv. Lett.</source> <volume>1</volume>, <fpage>37</fpage>&#x2013;<lpage>43</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1755-263X.2008.00002.x</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brodie</surname> <given-names>J. F.</given-names>
</name>
<name>
<surname>Giordano</surname> <given-names>A. J.</given-names>
</name>
<name>
<surname>Dickson</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Hebblewhite</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Bernard</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Mohd-Azlan</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>Evaluating multispecies landscape connectivity in a threatened tropical mammal community</article-title>. <source>Conserv. Biol.</source> <volume>29</volume>, <fpage>122</fpage>&#x2013;<lpage>132</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/cobi.12337</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brost</surname> <given-names>B. M.</given-names>
</name>
<name>
<surname>Beier</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Use of land facets to design linkages for climate change</article-title>. <source>Ecol. Appl.</source> <volume>22</volume>, <fpage>87</fpage>&#x2013;<lpage>103</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1890/11-0213.1</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cardillo</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Mace</surname> <given-names>G. M.</given-names>
</name>
<name>
<surname>Jones</surname> <given-names>K. E.</given-names>
</name>
<name>
<surname>Bielby</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Bininda-Emonds</surname> <given-names>O. R. P.</given-names>
</name>
<name>
<surname>Sechrest</surname> <given-names>W.</given-names>
</name>
<etal/>
</person-group>. (<year>2005</year>). <article-title>Multiple causes of high extinction risk in large mammal species</article-title>. <source>Science</source> <volume>309</volume>, <fpage>1239</fpage>&#x2013;<lpage>1241</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1126/science.1116030</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Caro</surname> <given-names>T. M.</given-names>
</name>
<name>
<surname>Girling</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2010</year>). <source>Conservation by proxy: indicator, umbrella, keystone, flagship, and other surrogate species</source> (<publisher-loc>Washington, DC</publisher-loc>: <publisher-name>Island Press</publisher-name>).</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Caro</surname> <given-names>T. M.</given-names>
</name>
<name>
<surname>O&#x2019;Doherty</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>On the use of surrogate species in conservation biology</article-title>. <source>Conserv. Biol.</source> <volume>13</volume>, <fpage>805</fpage>&#x2013;<lpage>814</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1046/j.1523-1739.1999.98338.x</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="confproc">
<person-group person-group-type="author">
<collab>CBD</collab>
</person-group> (<year>2010</year>). &#x201c;<article-title>Strategic plan for biodiversity 2011&#x2013;2020 and the aichi targets</article-title>,&#x201d; in <conf-name>Report of the Tenth Meeting of the Conference of the Parties to the Convention on Biological Diversity</conf-name>. <publisher-loc>Nagoya, Japan</publisher-loc>. Available at: <uri xlink:href="https://wedocs.unep.org/20.500.11822/27533">https://wedocs.unep.org/20.500.11822/27533</uri>.</citation>
</ref>
<ref id="B13">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>CBD</collab>
</person-group> (<year>2020</year>). <source>Update of the zero draft of the post-2020 global biodiversity framework (Convention on biological diversity)</source>. Available at: <uri xlink:href="https://www.cbd.int/doc/c/3064/749a/0f65ac7f9def86707f4eaefa/post2020-prep-02-01-en.pdf">https://www.cbd.int/doc/c/3064/749a/0f65ac7f9def86707f4eaefa/post2020-prep-02-01-en.pdf</uri> (Accessed June 14, 2023)</citation>
</ref>
<ref id="B14">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>CBD</collab>
</person-group> (<year>2023</year>) <source>2030 targets and guidance notes. secretariat of the convention on biological diversity</source>. Available at: <uri xlink:href="https://www.cbd.int/gbf/targets/">https://www.cbd.int/gbf/targets/</uri> (Accessed <access-date>June 14, 2023</access-date>).</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ceballos</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Ehrlich</surname> <given-names>P. R.</given-names>
</name>
<name>
<surname>Dirzo</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Biological annihilation via the ongoing sixth mass extinction signaled by vertebrate population losses and declines</article-title>. <source>Proc. Natl. Acad. Sci.</source> <volume>114</volume>, <fpage>E6089</fpage>&#x2013;<lpage>E6096</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1073/pnas.1704949114</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chase</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Blowes</surname> <given-names>S. A.</given-names>
</name>
<name>
<surname>Knight</surname> <given-names>T. M.</given-names>
</name>
<name>
<surname>Gerstner</surname> <given-names>K.</given-names>
</name>
<name>
<surname>May</surname> <given-names>F.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Ecosystem decay exacerbates biodiversity loss with habitat loss</article-title>. <source>Nature</source>. <volume>584</volume>, <fpage>238</fpage>&#x2013;<lpage>243</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41586-020-2531-2</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Crooks</surname> <given-names>K. R.</given-names>
</name>
<name>
<surname>Burdett</surname> <given-names>C. L.</given-names>
</name>
<name>
<surname>Theobald</surname> <given-names>D. M.</given-names>
</name>
<name>
<surname>King</surname> <given-names>S. R. B.</given-names>
</name>
<name>
<surname>Di Marco</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Rondinini</surname> <given-names>C.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>). <article-title>Quantification of habitat fragmentation reveals extinction risk in terrestrial mammals</article-title>. <source>Proc. Natl. Acad. Sci.</source> <volume>114</volume>, <fpage>7635</fpage>&#x2013;<lpage>7640</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1073/pnas.1705769114</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cushman</surname> <given-names>S. A.</given-names>
</name>
<name>
<surname>Landguth</surname> <given-names>E. L.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Multi-taxa population connectivity in the northern rocky mountains</article-title>. <source>Ecol. Model.</source> <volume>231</volume>, <fpage>101</fpage>&#x2013;<lpage>112</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ecolmodel.2012.02.011</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dinerstein</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Olson</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Joshi</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Vynne</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Burgess</surname> <given-names>N. D.</given-names>
</name>
<name>
<surname>Wikramanayake</surname> <given-names>E.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>). <article-title>An ecoregion-based approach to protecting half the terrestrial realm</article-title>. <source>BioScience</source> <volume>67</volume>, <fpage>534</fpage>&#x2013;<lpage>545</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/biosci/bix014</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dudaniec</surname> <given-names>R. Y.</given-names>
</name>
<name>
<surname>Rhodes</surname> <given-names>J. R.</given-names>
</name>
<name>
<surname>Worthington Wilmer</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Lyons</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Lee</surname> <given-names>K. E.</given-names>
</name>
<name>
<surname>McAlpine</surname> <given-names>C. A.</given-names>
</name>
<etal/>
</person-group>. (<year>2013</year>). <article-title>Using multilevel models to identify drivers of landscape-genetic structure among management areas</article-title>. <source>Mol. Ecol.</source> <volume>22</volume>, <fpage>3752</fpage>&#x2013;<lpage>3765</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/mec.12359</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>European Commission</collab>
</person-group> (<year>2020</year>) <source>EU Biodiversity strategy for 2030. bringing nature back into our lives. communication from the commission to the European parliament, the council, the European economic and social committee and the committee of the regions</source>. Available at: <uri xlink:href="https://ec.europa.eu/environment/nature/biodiversity/strategy/index_en.htm">https://ec.europa.eu/environment/nature/biodiversity/strategy/index_en.htm</uri> (Accessed <access-date>May 23, 2020</access-date>).</citation>
</ref>
<ref id="B22">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>European Environment Agency</collab>
</person-group> (<year>2016</year>) <source>Landscape fragmentation indicator effective mesh density (seff) - major and medium anthropogenic fragmentation (FGA2_S_2016)</source>. Available at: <uri xlink:href="https://www.eea.europa.eu/data-and-maps/data/landscape-fragmentation-indicator-effective-mesh">https://www.eea.europa.eu/data-and-maps/data/landscape-fragmentation-indicator-effective-mesh</uri> (Accessed <access-date>July 21, 2020</access-date>).</citation>
</ref>
<ref id="B23">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>European Environment Agency</collab>
</person-group> (<year>2020</year>) <source>Natura 2000 data - the European network of protected sites [[/amp]]mdash; European environment agency</source>. Available at: <uri xlink:href="https://www.eea.europa.eu/data-and-maps/data/natura-11">https://www.eea.europa.eu/data-and-maps/data/natura-11</uri> (Accessed <access-date>January 25, 2020</access-date>).</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fieberg</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Kochanny</surname> <given-names>C. O.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Quantifying home-range overlap: the importance of the utilization distribution</article-title>. <source>J. Wildl. Manage.</source> <volume>69</volume>, <fpage>1346</fpage>&#x2013;<lpage>1359</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2193/0022-541X(2005)69[1346:QHOTIO]2.0.CO;2</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Gilpin</surname> <given-names>M. E.</given-names>
</name>
<name>
<surname>Soul&#xe9;</surname> <given-names>M. E.</given-names>
</name>
</person-group> (<year>1986</year>). &#x201c;<article-title>Minimum viable populations: processes of species extinction</article-title>,&#x201d; in <source>Conservation biology: the science of scarcity and diversity</source>. Ed. <person-group person-group-type="editor">
<name>
<surname>Soul&#xe9;</surname> <given-names>M. E.</given-names>
</name>
</person-group> (<publisher-loc>Sunderland, MA</publisher-loc>: <publisher-name>Sinauer Associates</publisher-name>), <fpage>19</fpage>&#x2013;<lpage>34</lpage>.</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jones</surname> <given-names>K. E.</given-names>
</name>
<name>
<surname>Bielby</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Cardillo</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Fritz</surname> <given-names>S. A.</given-names>
</name>
<name>
<surname>O&#x2019;Dell</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Orme</surname> <given-names>C. D. L.</given-names>
</name>
<etal/>
</person-group>. (<year>2009</year>). <article-title>PanTHERIA: a species-level database of life history, ecology, and geography of extant and recently extinct mammals</article-title>. <source>Ecology</source> <volume>90</volume>, <fpage>2648</fpage>&#x2013;<lpage>2648</lpage>. doi: <pub-id pub-id-type="doi">10.1890/08-1494.1</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Krosby</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Breckheimer</surname> <given-names>I.</given-names>
</name>
<name>
<surname>John Pierce</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Singleton</surname> <given-names>P. H.</given-names>
</name>
<name>
<surname>Hall</surname> <given-names>S. A.</given-names>
</name>
<name>
<surname>Halupka</surname> <given-names>K. C.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>Focal species and landscape &#x201c;naturalness&#x201d; corridor models offer complementary approaches for connectivity conservation planning</article-title>. <source>Landsc. Ecol.</source> <volume>30</volume>, <fpage>2121</fpage>&#x2013;<lpage>2132</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10980-015-0235-z</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lechner</surname> <given-names>A. M.</given-names>
</name>
<name>
<surname>Sprod</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Carter</surname> <given-names>O.</given-names>
</name>
<name>
<surname>Lefroy</surname> <given-names>E. C</given-names>
</name>
</person-group>. (<year>2017</year>). <article-title>Characterising landscape connectivity for conservation planning using a dispersal guild approach</article-title>. <source>Landsc. Ecol.</source> <volume>32</volume>, <fpage>99</fpage>&#x2013;<lpage>113</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10980-016-0431-5</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Maiorano</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Amori</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Capula</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Falcucci</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Masi</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Montemaggiori</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2013</year>). <article-title>Threats from climate change to terrestrial vertebrate hotspots in Europe</article-title>. <source>PloS One</source> <volume>8</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.pone.0074989</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mateo-S&#xe1;nchez</surname> <given-names>M. C.</given-names>
</name>
<name>
<surname>Balkenhol</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Cushman</surname> <given-names>S.</given-names>
</name>
<name>
<surname>P&#xe9;rez</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Dom&#xed;nguez</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Saura</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Estimating effective landscape distances and movement corridors: comparison of habitat and genetic data</article-title>. <source>Ecosphere</source> <volume>6</volume>, <fpage>1</fpage>&#x2013;<lpage>16</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1890/ES14-00387.1</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Meurant</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Gonzalez</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Doxa</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Albert</surname> <given-names>C. H.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Selecting surrogate species for connectivity conservation</article-title>. <source>Biol. Conserv.</source> <volume>227</volume>, <fpage>326</fpage>&#x2013;<lpage>334</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.biocon.2018.09.028</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mitchell</surname> <given-names>M. G.</given-names>
</name>
<name>
<surname>Bennett</surname> <given-names>E. M.</given-names>
</name>
<name>
<surname>Gonzalez</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Linking landscape connectivity and ecosystem service provision: current knowledge and research gaps</article-title>. <source>Ecosystems</source> <volume>16</volume>, <fpage>894</fpage>&#x2013;<lpage>908</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s10021-013-9647-2</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>M&#xfc;llner</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Fastcluster: fast hierarchical, agglomerative clustering routines for r and Python</article-title>. <source>J. Stat. Software</source> <volume>53</volume>, <fpage>1</fpage>&#x2013;<lpage>18</lpage>. doi: <pub-id pub-id-type="doi">10.18637/jss.v053.i09</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Olson</surname> <given-names>D. M.</given-names>
</name>
<name>
<surname>Dinerstein</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Wikramanayake</surname> <given-names>E. D.</given-names>
</name>
<name>
<surname>Burgess</surname> <given-names>N. D.</given-names>
</name>
<name>
<surname>Powell</surname> <given-names>G. V. N.</given-names>
</name>
<name>
<surname>Underwood</surname> <given-names>E. C.</given-names>
</name>
<etal/>
</person-group>. (<year>2001</year>). <article-title>Terrestrial ecoregions of the world: a new map of life on earth</article-title>. <source>BioScience</source> <volume>51</volume>, <fpage>933</fpage>&#x2013;<lpage>938</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1641/0006-3568(2001)051[0933:TEOTWA]2.0.CO;2</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ribeiro</surname> <given-names>J. W.</given-names>
</name>
<name>
<surname>Silveira dos Santos</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Dodonov</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Martello</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Brand&#xe3;o Niebuhr</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Ribeiro</surname> <given-names>M. C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>LandScape corridors (<sc>lscorridors</sc>): a new software package for modelling ecological corridors based on landscape patterns and species requirements</article-title>. <source>Methods Ecol. Evol</source>. <volume>8</volume>, <fpage>1425</fpage>&#x2013;<lpage>1432</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/2041-210X.12750</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Robillard</surname> <given-names>C. M.</given-names>
</name>
<name>
<surname>Coristine</surname> <given-names>L. E.</given-names>
</name>
<name>
<surname>Soares</surname> <given-names>R. N.</given-names>
</name>
<name>
<surname>Kerr</surname> <given-names>J. T.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Facilitating climate-change-induced range shifts across continental land-use barriers</article-title>. <source>Conserv. Biol.</source> <volume>29</volume>, <fpage>1586</fpage>&#x2013;<lpage>1595</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/cobi.12556</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Santini</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Marco</surname> <given-names>M. D.</given-names>
</name>
<name>
<surname>Visconti</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Baisero</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Boitani</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Rondinini</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Ecological correlates of dispersal distance in terrestrial mammals</article-title>. <source>Hystrix Ital. J. Mammal.</source> <volume>24</volume>, <fpage>181</fpage>&#x2013;<lpage>186</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4404/hystrix-24.2-8746</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Temple</surname> <given-names>H. J.</given-names>
</name>
<name>
<surname>Terry</surname> <given-names>A.</given-names>
</name>
</person-group> (Compilers) (<year>2007</year>). <article-title>The status and distribution of European mammals. (Luxembourg: Office for Official Publications of the European Communities).</article-title>. Available at: https://portals.iucn.org/library/sites/library/files/documents/RL-4-013.pdf (Accessed July 18, 2020).</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Theobald</surname> <given-names>D. M.</given-names>
</name>
<name>
<surname>Reed</surname> <given-names>S. E.</given-names>
</name>
<name>
<surname>Fields</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Soul&#xe9;</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Connecting natural landscapes using a landscape permeability model to prioritize conservation activities in the united states: connecting natural landscapes</article-title>. <source>Conserv. Lett.</source> <volume>5</volume>, <fpage>123</fpage>&#x2013;<lpage>133</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1755-263X.2011.00218.x</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>F.</given-names>
</name>
<name>
<surname>McShea</surname> <given-names>W. J.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Does one size fit all? a multispecies approach to regional landscape corridor planning</article-title>. <source>Divers. Distrib.</source> <volume>24</volume>, <fpage>415</fpage>&#x2013;<lpage>425</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/ddi.12692</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ward</surname> <given-names>J. H.</given-names>
</name>
</person-group> (<year>1963</year>). <article-title>Hierarchical grouping to optimize an objective function</article-title>. <source>J. Am. Stat. Assoc.</source> <volume>58</volume>, <fpage>236</fpage>&#x2013;<lpage>244</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/01621459.1963.10500845</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Whitmee</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Orme</surname> <given-names>C. D. L.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Predicting dispersal distance in mammals: a trait-based approach</article-title>. <source>J. Anim. Ecol.</source> <volume>82</volume>, <fpage>211</fpage>&#x2013;<lpage>221</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1365-2656.2012.02030.x</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wiens</surname> <given-names>J. A.</given-names>
</name>
<name>
<surname>Hayward</surname> <given-names>G. D.</given-names>
</name>
<name>
<surname>Holthausen</surname> <given-names>R. S.</given-names>
</name>
<name>
<surname>Wisdom</surname> <given-names>M. J.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Using surrogate species and groups for conservation planning and management</article-title>. <source>BioScience</source> <volume>58</volume>, <fpage>241</fpage>&#x2013;<lpage>252</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1641/B580310</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>Wildlife Conservation Society-WCS</collab>
<collab>Center for International Earth Science Information Network-CIESIN-Columbia University</collab>
</person-group> (<year>2005</year>). <source>Last of the wild project, version 2, 2005 (LWP-2): global human footprint dataset (Geographic)</source>. <publisher-loc>Palisades, New York</publisher-loc>: <publisher-name>NASA Socioeconomic Data and Applications Center (SEDAC)</publisher-name>. doi:&#xa0;<pub-id pub-id-type="doi">10.7927/H4M61H5F</pub-id>. Accessed May 23, 2020.</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wood</surname> <given-names>S. L. R.</given-names>
</name>
<name>
<surname>Martins</surname> <given-names>K. T.</given-names>
</name>
<name>
<surname>Dumais-Lalonde</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Tanguy</surname> <given-names>O.</given-names>
</name>
<name>
<surname>Maure</surname> <given-names>F.</given-names>
</name>
<name>
<surname>St-Denis</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Missing interactions: the current state of multispecies connectivity analysis</article-title>. <source>Front. Ecol. Evol.</source> <volume>10</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fevo.2022.830822</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zeller</surname> <given-names>K. A.</given-names>
</name>
<name>
<surname>McGarigal</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Whiteley</surname> <given-names>A. R.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Estimating landscape resistance to movement: a review</article-title>. <source>Landscape Ecology</source> <volume>27</volume>, <fpage>777</fpage>&#x2013;<lpage>797</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s10980-012-9737-0</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>
