Dynamic connectivity analyses to inform management of the invasive American mink and its native competitor, the European mink

Forecasting habitat suitability and connectivity can be central to both controlling range expansion of invasive species and promoting native species conservation, especially under changing climate conditions. This study aimed to identify and prioritize areas in Spain to control the expansion of one of the most harmful invasive species in Europe, the American mink, while conserving its counterpart, the endangered European mink, under current and future conditions. We used ensemble habitat suitability and dynamic connectivity models to predict species ranges and movement routes considering likely climate change under three emission scenarios. Then, using habitat availability metrics, we prioritized areas for invasive mink control and native mink conservation and classified them into different management zones that reflected the overlap between species and threat from American to European minks. Results suggest that both species are likely to experience declines in habitat and connectivity under climate change scenarios with significantly larger declines by the end of the century for European minks (72 and 80% respectively) than for American minks (41 and 32%). Priority areas for management of both species varied over time and across emission scenarios, with a general shift in priority habitat towards the North-East of the study area. Our findings demonstrate how habitat suitability and dynamic connectivity approaches can guide long-term management strategies to control invasive species and conserve native species while accounting for likely landscape changes. The simultaneous study of both invasive and native species can support prioritized management action and inform management planning of the intensity, extent, and techniques of intervention depending on the overlap between species.


Introduction
Climate change is projected to alter the threats posed by invasive species to native species conservation (Vilà et al. 2007;Hellmann et al. 2008;Polaina et al. 2021).Invasive species endanger native species through predation, competition for resources, and disease transmission (Genovesi et al. 2012).
Vol:. ( 1234567890) Because invasive species are frequently more generalist in their habitat selection and adapt more easily to changing landscape conditions than native specialist species, climate change has been found to exacerbate the detrimental effects posed by invasive species (Ravi et al. 2022).Specifically, shifts in climate have the potential to alter the introduction, establishment, distribution, and spread of invasive species (Polaina et al. 2021).These changes may also affect invasive species' impact on native species, altering the distribution and degree of vulnerability of the native species and areas of overlap of both species (Hellmann et al. 2008).To respond to these potential shifts, both invasive species control measures and native species conservation may require dynamic approaches that can support responsive management actions across changing conditions to prevent the establishment, spread, and threat of invasive species and preserve native species over time.
The impacts of climatic changes on native and invasive species' ability to thrive depends largely on the availability, distribution, suitability, and connectivity of their habitats.Habitat suitability and connectivity modeling are common approaches used to locate and quantify a species' potential distribution and movements as a function of multiple biotic and abiotic environmental factors such as climate and land use or cover (Correa Ayram et al. 2015;Guisan et al. 2017).As these factors change over time (Mora et al. 2013;Song et al. 2018), so do habitat suitability and landscape connectivity (Saura et al. 2011;Rubio et al. 2012;Carvalho et al. 2021;Goicolea and Mateo-Sánchez 2022).Despite the widespread recognition that habitat suitability and associated connectivity are dynamic processes, connectivity has been more commonly modeled as a static, immutable feature (Zeller et al. 2020).However, several studies have found that static approaches to connectivity modeling fail to effectively quantify changes in landscape connectedness or identify priority areas for management measures under changing conditions (Wimberly 2006;Bishop-Taylor et al. 2018;Jennings et al. 2020;Goicolea and Mateo-Sánchez 2022).Recent advances in dynamic connectivity modeling approaches account for changes in landscape suitability and connectivity stability and can be used to identify persistent areas that are likely to support species expansion to potential future ranges (Keeley et al. 2018;Zeller et al. 2020;Goicolea and Mateo-Sánchez 2022).Some of these dynamic approaches generate predictions over multiple time periods and use future climate predictions to forecast potential changes in species expansion and adaptation to climate and related landscape changes (Beltrán et al. 2014;Bishop-Taylor et al. 2018).These dynamic models can be powerful tools to identify areas likely to persist important for focal species connectivity and thus guide relevant and sustainable conservation measures targeting long-term management (Jennings et al. 2020;Goicolea and Mateo-Sánchez 2022).
Despite the well-documented importance of maintaining or enhancing habitat integrity and connectivity in the face of habitat loss and fragmentation (Rudnick et al. 2012;Correa Ayram et al. 2015), there can also be negative consequences of landscape connectivity.Connectivity can facilitate the spread and expansion of unwanted disturbances such as catastrophic fires, pests, or invasive species, increasing the threats to species of conservation concern (Drake et al. 2017;de la Fuente et al. 2018).Previous studies have mainly focused on enhancing the distribution and movements of specialist, native, or threatened species yet rarely focus on limiting expansion of detrimental processes.However, concurrent modeling of species of conservation concern and the threats they face can be used to better identify management and conservation actions to mitigate threats and protect at-risk species.This concurrent modeling has the potential to be particularly beneficial in managing native species and non-native invasive competitors.To delineate priority conservation action areas (Wilson et al. 2009) where management can promote the persistence of native species and limit the impacts and expansion of invasive competitors, dynamic connectivity models can be combined with habitat availability metrics (Saura and Pascual-Hortal 2007;Saura and Rubio 2010).Habitat availability metrics are particularly useful for informing conservation and management efforts because they are flexible and can assess adequacy, permanence, and irreplaceability of habitat patches and linkages, three important criteria in conservation prioritization.Adequacy denotes how suitable and connected areas are to favor species populations' survival and growth.These characteristics can arise from habitat suitability and connectivity models and depend on the environmental conditions and location of the area.The second criterion is the continuity or permanence of the adequacy of these conservation areas over time.In this sense, dynamic connectivity approaches are essential to identify areas likely to support persistently suitable and connected habitat even under projected changes.Lastly, habitat availability metrics can also be used to consider the irreplaceability of areas, i.e., whether alternative landscape elements could contribute similarly to the conservation goal.Therefore, these habitat availability metrics can be useful to guide conservation and management efforts to enhance habitat suitability and connectivity for native species and reduce that of invasive species under changing conditions.
Here, we identified and compared current and future priority areas to conserve a critically endangered native species, the European mink (Mustela lutreola), and control the expansion of its invasive competitor, the American mink (Neovison vison), in Spain.Specifically, we studied the current and potential future distribution and functional connectivity of both species with habitat suitability and dynamic connectivity models accounting for the adequacy, permanence, and irreplaceability of the areas.To address the uncertainty associated with future predictions, we employed multiple statistical techniques and evaluated suitability and connectivity over four time under three emissions scenarios.Additionally, we integrated habitat availability and connectivity results for both species to identify the most important areas to control American minks and to preserve European minks over this century.

Species and study area
The study area is Spain (excluding insular areas) in southwest Europe (Fig. 1).It occupies 492,175 km 2 in area, and is composed of the Mediterranean forests, woodlands and scrubs, and temperate broadleaf and mixed forests biomes (Olson et al. 2001).Spain is expected to experience large climate and landuse changes (Loarie et al. 2009;Song et al. 2018) and these changes may affect population dynamics of both native and invasive species.In fact, Spain is considered one of the regions with the highest concentration of vulnerable terrestrial species to climate change (Pacifici et al. 2015).
The European mink (Mustela lutreola) has been recognized as one of the most endangered mammals in the world (Maran et al. 2016).In the mid 1800's, European mink were found throughout most of continental Europe (Youngman 1990).However, in the last years, over 85% of their range has disappeared, and is currently restricted to small, isolated, and decreasing populations in eastern Europe, southern France, and northern Spain.The causes of this decline are a combination of invasive species competition and disease transmission, habitat loss and fragmentation, water pollution, and overhunting (Maran and Henttonen 1995;Mañas et al. 2001;Zuberogoitia et al. 2013).European minks in Spain (Fig. 1) occupy less than 2300 km of watercourses and have a population size of about 500 individuals with low genetic variability (Michaux et al. 2005;Palazón and Melero 2014;Cabria et al. 2015).
Another mink species, the American mink (Neovison vison) is one of the most harmful exotic and invasive species in Europe (Genovesi et al. 2012).Both mink species are small, semi-aquatic, and carnivorous mustelids with considerable niche overlap in food consumption (Sidorovich et al. 2010;Zuberogoitia et al. 2013) and preferred habitat, i.e. habitat adjacent to water with high, dense vegetation cover (Fuller et al. 2016;Maran et al. 2016).However, American minks are more adaptable to poor habitat quality, in addition to being bigger, more aggressive, and more successful breeders than European minks (Palazón and Melero 2014;Zuberogoitia et al. 2014).American minks were introduced in Europe from North America in the first half of the twentieth century for fur farming.Escaped or released individuals formed wild populations, and rapidly spread to be widely distributed around Europe.In Spain, American minks have established feral widespread populations across 12,530 km of rivers with a population of over 30,000 individuals (Fig. 1., MAGRAMA 2014).The establishment and spread of the American mink is considered one of the main causes of decline of European mink (Sidorovich et al. 2010;Maran et al. 2016), with similar negative impacts on other native species of conservation concern such as European polecats (Mustela putorius), stoats (Mustela erminea), and sea birds (MacDonald and Harrington 2003;Ahola et al. 2006;Banks et al. 2008;Schüttler et al. 2009;Melero et al. 2012;García-díaz et al. 2013).American mink has also been linked to economic damage

Species occurrences and environmental variables data
Fig. 1 The study area (Spain) and estimated distribution of the invasive (American mink) and the endangered native (European mink) species as of 2012 according to the Spanish inventory of terrestrial species (MITECO 2012) Data for both mink species were obtained from the yearly updated Spanish inventory of terrestrial species (MITECO 2012).We gathered 1396 American and 197 European mink occurrences at 10-km resolution (Fig. 1).The bioclimatic variables (Table 1) were obtained from Climatologies at High resolution for the Earth's Surface Areas (CHELSA) data (Karger et al. 2017) with 30 arcsec resolution but resampled to species occurrences resolution (10 km).We downloaded these bioclimatic variables for the historic period of 1985-2010 and three future time periods 2011-2040, 2041-2070, 2071-2100 (hereafter referred to as t1, t2, t3, and t4).The variables of the future periods were derived from the GFDL-ESM4 (National Oceanic and Atmospheric Administration, Geophysical Fluid Dynamics Laboratory, Princeton, NJ 08540, USA) Global Climatic Model under three different shared socio-economic and emissions pathways (SSP): 126, 375 and 585 (hereafter referred to as e1, e2, and e3) corresponding to the Representative Concentration Pathways (RCP) 2.6 (an optimistic and sustainable scenario), 7 (medium-high emissions scenario), and 8.5 (high fossil-fuel development) (Meinshausen et al. 2020).In total, we developed models for 10 scenarios (one historic and nine future).
Although climate variables are likely to influence species habitat selection at a coarse scale (Wiens 1989;Pearson et al. 2002;Thuiller et al. 2004), prior studies have shown that both mink species respond to finer scale habitat structure (Yamaguchi et al. 2003;Harrington and Macdonald 2008;Crego et al. 2018).Therefore, we refined our models with land use and cover variables previously identified as affecting mink habitat suitability at a finer resolution (Yamaguchi et al. 2003;Fuller et al. 2016).Specifically, we incorporated 2018 CORINE land cover maps (European Union, Copernicus Land Monitoring Service 2018) at 100-m resolution complemented with roads (ESRI and Garmin International 2021), rivers (Ministerio para la Transición Ecológica y el Reto Demográfico 2018), and tree cover density (European Environment Agency 2018).We lacked fine-scaled future forecasts for landscape variables and assumed static land use and cover over time (Costanza and Terando 2019) even though adding accurate information on land use and human pressure is of utmost importance in species' habitat selection (Parks et al. 2020).

Habitat suitability modeling
To assess connectivity dynamics over time, we started by modeling snapshots of habitat suitability at different time periods and under our different emission scenarios.Habitat suitability was calculated as the product of climate and landscape suitability measurements (Fig. 2).We calculated climate suitability by modeling species distribution based on the relationship between species occurrences and bioclimatic variables.Landscape suitability was calculated combining land use, tree cover density, and distance to water information.We determined the relationships to these landscape variables based on published literature (Yamaguchi et al. 2003;Zschille et al. 2012;Fuller et al. 2016;Halbrook and Petach 2018) as the resolution of the species occurrences (10 km) was considered too coarse to be relevant in species land use and cover selection modeling.
We calculated climatic suitability for each species following an ensemble species distribution modeling approach (Araújo and New 2006) considering relevant bioclimatic variables (Table 1 and Appendix S1) with the biomod2 R package (Thuiller et al. 2021).Ensemble models combine a set of individual models to account for algorithm biases (Grenouillet et al. 2011) and limit the influence of uncertainties associated with a single statistical technique for future forecasts (Araújo and New 2006;Beaumont et al. 2016).In this study, we combined six different models representing a range of statistical techniques: generalized linear models (GLM) (McCullagh and Nelder 1989), generalized additive models (GAM) (Hastie and Tibshirani 1986), generalized boosting models (GBM) (Friedman 2001), maximum entropy (MAXENT) (Phillips et al. 2006), multivariate adaptive regression splines (MARS) (Friedman 1991), and random forests (RF) (Breiman 2001).Background points were randomly selected across the study area 10-km grid at a 3:1 ratio with occurrence points (Barbet-Massin et al. 2012).We assumed there was no bias in the species occurrence data, as they were reported at a coarse, grid-based scale across all of Spain.Models were calibrated with the historic bioclimatic variables.We assessed the performance of the individual models with tenfold cross-validation (Pearce 2000) using the area under the curve (AUC) performance metric (Fan;et al. 2006).To assemble the final ensemble model for each species, we retained cross-validation runs of individual models with an AUC greater than 0.75 and averaged them in proportion to their AUC (weighted average).We also evaluated the ensemble model performance using the same cross-validation procedure, and the AUC, TSS, and Kappa evaluation statistics.Lastly, we projected the climate suitability model with clamping extrapolation to every time period (t1, t2, t3, t4) and SSP scenario (e1, e2, and e3) producing ten different climate suitability surfaces.To measure the risk of extrapolation of these projections (Owens et al. 2013), we developed a mobility-oriented parity analysis (MOP) to compare the calibration of climatic conditions (the variables at the historic scenario) with the transferred climate conditions (climate variables at future scenarios) with the kuenm R package (Cobos et al. 2019).
To calculate landscape suitability, we considered three factors within each 100 m pixel across our study area: type of land use, tree cover density, and distance to rivers and water bodies (Fuller et al. 2016).We assigned land use suitability according to the degree of human pressure for each land cover type (Table S2 in supplementary material): artificial land uses were considered unsuitable (0), agricultural uses moderately suitable (0.75), and natural areas assigned maximum suitability (1) (Fuller et al. 2016).For the tree cover component of landscape suitability, we assigned suitability values according to quartiles of tree cover density (Yamaguchi et al. 2003): a coverage of 100-75% was associated with the highest suitability (1), 75-50% with 0.75 suitability, 50-25% with 0.5 suitability, and 0-25% as unsuitable (0).For the final component of landscape suitability, we transformed distance to water into a suitability percentage with a negative exponential function assigning pixels adjacent to rivers or water bodies with a maximum suitability value (1) and with values decreasing exponentially for pixels farther from water (Zschille et al. 2012;Fuller et al. 2016;Halbrook and Petach 2018).To fit the exponential function, we assigned a suitability value of 0.05 to 900 m, the maximum distance of American minks to water sources reported in the literature (Crego et al. 2018).To calculate the final habitat suitability value of each cell, we multiplied the climate suitability and the three components of landscape suitability.We obtained ten different habitat suitability surfaces at 100 m resolution, one for each period and emissions scenario considered.We also transformed the habitat suitability surfaces to binary surfaces of habitat or non-habitat using the true skill statistic (TSS) threshold (Somodi et al. 2017) of the climate suitability models.We used this binary surface to delineate habitat patches (Fig. 3) in ArcGIS software (Esri, Redlands, California, USA) by aggregating suitable habitat cells within 1 km and assembling polygons of habitat with a minimum area of 70 ha.This 70 ha threshold is the estimated minimum area required for populations' long-term survival, reproduction, and proliferation, and was calculated as the longitudinal home range (7 km) multiplied by the pixel width (100 m).We verified that previous studies obtained similar home ranges (Yamaguchi and Mac-Donald 2003;Zabala et al. 2007;Fournier et al. 2008;Harrington and Macdonald 2008;Peters et al. 2009;Zschille et al. 2012;Palomares et al. 2017;Halbrook and Petach 2018).Lastly, we used the area of these habitat patches to estimate species abundance assuming a proportional relationship to patch size (Drake et al. 2021).Aside from the ten sets of habitat patches for each species (one per period and scenario), we calculated the overall potential habitat by aggregating the habitat patches of all periods and scenarios to represent the total area that could be colonized at any time by each species.

Habitat connectivity modeling
We used the delineated habitat patches and resistance surfaces derived from the habitat suitability surfaces as the basis for least-cost analyses and connectivity metrics (Fig. 3).We calculated the resistance surface, which represents the difficulty for the species to move across each cell of the landscape matrix (Spear et al. 2010;Zeller et al. 2012) for each species, period, and SSP scenario using the negative exponential transformation described in Keeley et al. (2016).We considered that both species would likely be more tolerant of moderately suitable areas during dispersal than when selecting habitat; therefore, we set an intermediate value for the parameter c (c = 2) to transform from habitat suitability to resistance to movement surfaces as the inverse exponential of the habitat suitability surfaces (Mateo-Sánchez et al. 2015;Keeley et al. 2016;Jennings et al. 2020) with the formula: We then calculated potential links or corridors (least-cost paths) (Adriaensen et al. 2003), as well as their Euclidean and effective (cost weighted) distances between patch edges and across the resistance surface with Linkage Mapper software (McRae and Kavanagh 2011).We set a cutoff maximum bounding circle buffer distance of 50 km (Harrington et al. 2014;Oliver et al. 2016) around every habitat patch to increase computability.We also obtained the ratio of the effective to Euclidean distances of all links and calculated the mean of all the ratios in the historic scenario as an estimate to transform from Euclidean to effective distance.
To evaluate connectivity over time for both species, we calculated two connectivity metrics: the overall connectivity of the ecological network and the importance (adequacy and irreplaceability) of each patch and link to maintain it.To calculate these metrics, we used graph-based habitat availability metrics (Appendix S3 in supplementary material, Pascual-Hortal and Saura 2006; Saura and Pascual-Hortal 2007; Saura and Rubio 2010) with the command line version of Conefor software (Saura and Torné 2009, www. conef or. org).These connectivity metrics were calculated separately for each species, time period, and SSP scenario.They account for both the habitat available within a patch (intrapatch connectivity) and the habitat available through the connections with other patches (interpatch connectivity).

Overall landscape connectivity
To calculate overall landscape connectivity, we used the Equivalent Connected Area index (ECA), which is defined as the area of a single habitat patch that would provide the same level of connectivity as the actual pattern of habitat patches in the whole landscape (Saura et al. 2011).This metric is a function of the dispersal capacity of the focal species (given in this case by its maximum dispersal distance), the landscape configuration (represented by the effective distance between patches), and patch attributes (in this case patch area) and is directly comparable among scenarios.We used a maximum dispersal distance of 40 km (Mitchell 1961;Harrington et al. 2014;Oliver et al. 2016) converted into an effective distance (i.e., multiplied by the mean effective to Euclidean ratio obtained previously).We matched this distance to a probability of 0.05 to fit a negative exponential function that represented the probability of connection between each pair of patches in relation to the effective distance between them.

Priority conservation and control areas
To assess patch importance, we calculated the change in the overall connectivity when systematically removing each patch (see appendix S3 in the supplementary material and Saura and Pascual-Hortal 2007).The most important (adequate and irreplaceable) patches are those whose removal would imply a bigger loss of overall connectivity.We calculated patch importance for each of the ten scenarios, but also calculated a generalized patch importance value by min-max normalizing and summing the importance across all the scenarios for each species to find the areas that remained stable in their importance despite global changes.
To most effectively manage areas of overlap between competing native and invasive species, it is important to evaluate changes in habitat suitability and connectivity relative to one another.To do that, we integrated the results for both species and classified potential habitat into three different management zones: the overlap zone, the expansion zone, and the low-risk zone.The overlap zone corresponds to the areas where overall potential habitat of both species co-occur and therefore management of both species should be considered as top priority.The expansion zone represents areas of overall potential American mink habitat that do not directly overlap with potential habitat of European minks, but is within a buffer of 40 km, the maximum dispersal distance.The expansion zone is where controlling the spread of American minks is most important.Finally, we defined the low-risk zone as the area where American mink overall potential habitat exceeds the maximum dispersal distance (> 40 km) from European mink habitat, or European mink overall potential habitat does not overlap with that of American minks.Within these three management zones, we ranked the importance of the areas based on generalized patch importance.Lastly, we determined the percentage of overlap of each management zone within the network for existing protected areas (UNEP-WCMC and IUCN 2022).

Habitat suitability
Every individual and ensemble climate suitability model achieved an AUC greater than 0.85 (Table S4).In the historic period (t1), the amount of suitable area (total patch area) was 19,418 km 2 for the American mink and 6171 km 2 for the European mink.Suitable habitat was predicted to decrease for both species across all future time periods and emission scenarios (Fig. 4 and Table S6).However, the pattern of decline was different for the two species.The suitable area for the American mink was predicted to decrease more drastically in the two first future periods (t2 and t3) while the biggest decrease in habitat suitability for the European mink was predicted at the end of the century (t4).Overall, results showed a greater proportional decrease in habitat amount for the European minks; the mean decrease in suitable area across the  , t2, t3, and t4 (1985-2010, 2011-2040, 2041-2070, 2071-2100) and species (American and European mink).Overall connectivity is measured with the metric Equivalent Connected Area (ECA).The three future periods (t2, t3, and t4) have three emissions scenarios.In A.1 and B.1, central lines show the mean across the three emission scenarios and the shaded areas represent the range between the minimum and the maximum habitat area (A) and ECA (B) according to the three scenarios.A.2 and B.2 show the proportion of suitable area and ECA of future periods (t2, t3, and t4) in relation to the historic period (t1) Vol:. ( 1234567890) three emission scenarios by t4 was 41% for American minks and 72% for European minks.For both species, the number of habitat patches generally decreased over time and with the increase of expected emissions.The size of these habitat patches also declined with the time and emissions for the European mink while the size of American mink patches was more variable (Table S6 of supplementary material).The overlap in suitable habitat remained constant across all periods and emission scenarios (Fig. S7).The overall potential habitat aggregated across all time periods and emission scenarios was 2,756,673 ha and 828,606 ha (141% and 134% of the amount of suitable area in the historic period) for American and European mink respectively.The risk of extrapolation increased with the time period and emissions scenario (Fig. S5).

Connectivity analyses
The effective to Euclidean distance ratio was 78 for the American and 71.6 for the European mink, and the maximum effective dispersal distances were 2864 km and 3120 km for each species respectively.In mapping our connectivity results, we observed that links connecting habitat patches for one species often crossed habitat patches or links of the other species (Fig. S8).In fact, in the historic period t1 26.3% of European mink links passed through American mink habitat patches or links.The predicted overall habitat connectivity (ECA) decreased for both species over time (Fig. 4).American mink's greatest loss of habitat connectivity was from time t2 to time t3 while for the European mink was from time t3 to time t4.By the end of the century, the impacts of climate change are expected to be worse for the European mink: connectivity decreased on average across the three emission scenarios 32% for American minks and 80% for European minks by t4.Also, increasing emissions were generally associated with a greater decline in habitat connectivity for both species.

Prioritization of action areas
We identified clear separation in the priority areas to maintain the overall species connectivity (i.e., areas with the greatest generalized patch importance) for the two species, with the northwest of Spain being most important for American minks and north central for European minks (Fig. 6).These priority areas for both species are predicted to shrink towards the north and east over time (Fig. S9) and with increasing emissions (Fig. S10).
Figure 6 shows the three zones of conservation and the order of priority within each zone to control American minks and to preserve European minks.Of the overall potential habitat for American mink, the overlap zone represented 18%, the expansion zone 25%, and the low threat zone 57%.17.4% of this overall potential habitat is inside Protected areas (19% in the overlap zone, 27% in the expansion zone, and 13% in low threat zone).For the European mink, 60% of their overall potential habitat belongs to the overlap zone and 40% to the low exposure zone.19.5% of the potential European mink habitat is protected (19% in the case of the overlap zone and 20% in the low exposure zone).

Discussion
Climate change, loss of connectivity through habitat fragmentation, and establishment and expansion of invasive species are interrelated phenomena that pose a threat to biodiversity and challenge conservation efforts.Effectively planning for and responding to these threats requires consideration of the complex and dynamic nature of the interactions among them.For example, enhancing connectivity is one of the most recommended strategies to support species' adaptation to climate change (Welch 2005;Heller and Zavaleta 2009;Keeley et al. 2018), but it also has the potential to facilitate the spread of invasive species.Additionally, climate change may affect both the threat posed by invasive species on ecosystems and native species (Hellmann et al. 2008;MITECO 2020) and species' habitat availability and connectivity (Keeley et al. 2018;Zeller et al. 2020;Goicolea and Mateo-Sánchez 2022).In this study, we presented a spatially explicit and dynamic approach to integrate all three phenomena to guide long-term and climatewise measures to preserve a native species and control its invasive competitor in a changing landscape.The combination of the habitat suitability and connectivity results from multiple time periods and scenarios allowed us to identify consistently important areas likely to enhance current and future habitat connectivity between contemporary patches, but also areas and routes that promote species' expansion to new areas (Goicolea and Mateo-Sánchez 2022).By focusing on one of the most worrying invasive species in Europe, the American mink, and one endangered native species affected by its invasion, the European mink, our findings demonstrated the importance of dynamic, integrated analyses.This approach also serves as a transferable framework to support to climate-wise conservation planning for other invasive species and their threatened native counterparts.
Decreasing habitat suitability and connectivity for both species under climate projections Our results confirmed that habitat availability and connectivity for both native and invasive species are dynamic phenomena likely to be affected by projected climatic changes over time (Figs. 4 and 5).
By 2071-2100 European mink habitat is expected to experience more substantial negative impacts under changing climate conditions, probably due to their higher degree of habitat selection specialization (Pacifici et al. 2015).American minks were also projected to experience habitat reduction, even higher than that of European minks between 2041 and 2070, which is notable given their ecological plasticity (Palazón and Melero 2014;Zuberogoitia et al. 2014) and the observed rapid range expansion in recent years (Põdra and Gómez 2018).This result suggests that even generalist species currently undergoing range expansion might be substantially affected by global changes.In this sense, climate change is likely to lead to a gradual decline in the amount of suitable area and potentially the distribution of American minks in Spain.Other studies have also predicted a contraction of the American mink range in Europe in future decades (Polaina et al. 2021).The period with the highest predicted loss of habitat and connectivity for the American mink (2041-2070) may offer an opportunity to intensify the control measures during a time when populations contract and control measures might be more effective.Additionally, minimizing the impact of the invasive species during this time could reduce other stressors and help European minks adapt to the greater loss of habitat availability and connectivity predicted by the end of the century (2071-2100).

Addressing future uncertainties
Climate suitability and connectivity models projected to future climate conditions can be associated with several sources of uncertainty (Heikkinen et al. 2006;Pacifici et al. 2015;Keeley et al. 2021).First, these models produce very different prediction and performance values depending on the statistical method used (Grenouillet et al. 2011;Pacifici et al. 2015;Beaumont et al. 2016).Given this variability, we applied multiple statistical methods and assembled their forecasts to reduce the uncertainty associated with any single modeling approach (Araújo and New 2006).Another common source of uncertainty in ecological models is the existence of future environmental conditions outside the training or calibration range (Zurell et al. 2012).Model projections to such novel future conditions can lead to erroneous extrapolations.Making cross-validation model evaluations is a recommended strategy to test how well the model performs when extrapolating outside the training range (Pearce and Ferrier 2000;Guisan et al. 2017).It is also recommended to check for the existence of novel environmental conditions to identify the areas where predictions could be spurious (Zurell et al. 2012;Owens et al. 2013).In this study we identified these areas with an MOP analysis (Fig. S5 in supplementary material).
Habitat and connectivity models projected to future scenarios are also subject to climatic uncertainties as different general circulation models, time periods, and emission scenarios forecast fluctuating climate conditions (Fernández et al. 2017).In fact, we obtained varying results depending on the time and emission scenario for both species (Figs. 4 and 5).To address this source of uncertainty, we separately considered the different periods and emission scenarios to obtain a range of potential outcomes that could be used to inform management strategies.The climatic uncertainty was also accounted for when prioritizing conservation areas.We found that the key areas to preserve the native species and control the invasive one did not remain static over time (Fig. S9) and the emissions scenario considered (Fig. S10).We found a trend of habitat reduction towards the east of the study area for both species over time and with higher emissions (Fig. S10).The similarities in the patterns we observed across time and emission scenarios suggest that higher greenhouse gas emissions may lead to an acceleration of climatic changes and associated effects.The independent study of habitat suitability at multiple periods let us identify both emerging and disappearing habitat and corridors over time that may inform improved management strategies that account for habitat and connectivity persistence.In the case of the invasive species, the early detection of novel habitats and corridors connecting them may guide more efficient measures that prevent invasive species' spread and establishment in advance instead of controlling and eradicating them from already invaded areas.It can also be useful to identify ephemerally suitable areas for the invasive or the native species to avoid focusing control and conservation resources and efforts on only temporarily suitable areas.
While the individual results from each time period and scenario informed our understanding of habitat suitability and connectivity dynamics for each species, the integration of these results into a generalized adequacy and irreplaceability value allowed us to account for the likelihood of permanence of these parameters over time for both species.Those areas that appeared to be highly adequate (suitable and connected) and irreplaceable in multiple periods and scenarios are more likely to keep being important for the species in the long-term.Additionally, priority areas that enhance connectivity at multiple time periods may also promote connectivity of species to their future suitable ranges which could inform management of expansion routes of both invasive and native species (Goicolea and Mateo-Sánchez 2022).
Integrating native and invasive species modeling to prioritize management actions Since the invasive and the native species in our study have similar habitat requirements (Zuberogoitia et al. 2013(Zuberogoitia et al. , 2014;;Zschille et al. 2014;Fuller et al. 2016;Halbrook and Petach 2018), their distribution may overlap and they may be affected similarly by the upcoming global changes and by management measures.The results from this and similarly designed studies can be used in planning for control of invasive species to avoid negatively affecting native species populations.Equally, these results can guide preservation measures that avoid promoting the expansion and improvement of invasive species populations.One main feature of this study is the prioritization of conservation and control areas considering the invasive and the threatened species together, to avoid misplacing resources that would have unintentional results on the other species.Despite the different rates of habitat loss for the two species and the past increasing overlap distributions between the two species found in previous studies (Põdra and Gómez 2018), here we found nearly constant overlapping habitat patches (nearly one-third) between the two species over time (Fig. S7), consistent with previously observed degrees of overlap (Põdra and Gómez 2018).However, when we consider the overall potential habitat aggregated across all time periods and emission scenarios, a much higher amount of habitat for European mink, nearly 60%, overlaps with that of American minks.Therefore, some areas may be suitable for both species in different periods.
We classified different management zones to control American minks and to preserve the European mink depending on the threat posed to the native species (Fig. 6).These categories divide the study area Vol:.( 1234567890) in localized areas that could be subjected to different timing, intensities of intervention, and types of measures to manage both the invasive and the native species (MAGRAMA 2014).The overlap zone comprises areas with likely direct competition between both species that may lead to a displacement or population decline of the native species.Managers should prioritize locations and management actions with higher intensities in these overlap zone to avoid further deleterious consequences for European minks.Managers can also target more holistic approaches in overlapping zones involving both measures to preserve and enhance European mink populations and to remove American minks (MacDonald and Harrington 2003).Measures to preserve the European mink in these areas will likely fail to protect the species if they are not combined with measures to reduce the threat of American minks.Additionally, control measures such as live-trapping and selective control procedures could be employed more judiciously and carried out by trained and experienced staff in the overlap zone due to the risk of harming European minks (MacDonald and Harrington 2003;Reynolds et al. 2013).These intensive and selective measures are usually expensive and difficult to implement and therefore are more likely to be successful when planned for localized and manageable areas such as this overlap management zone.More practical measures at broad scales could be used on the expansion, and low-risk zones where there is no threat of overlap between species.Control measures in the expansion and low-risk zones could include removal traps, castration measures, the use of dogs (Fuller et al. 2016), scent repellents (Baker and Macdonald 1999), or fences to avoid their expansion or access to specific areas (e.g., farms).However, careful consideration of the potential effects of management actions on other mustelid or native species remains a critical step in evaluating tools for controlling invasive species.The overall potential habitat of European minks inside the low-risk zone could undergo conservation measures that directly benefit European minks, such as habitat restoration, and recovery and artificial reintroduction of the species (Skorupski 2020).Other measures however should be applied throughout the study area to avoid further increasing the competition and disease transmission from American to European minks, such as the prevention of American mink escapes (MAGRAMA 2014) and eradication of diseases in fur farms (Mañas et al. 2001).
Our categorization of the priority areas aligns with the goals of the Spanish strategy to control American minks (MAGRAMA 2014) that differentiates two intensity levels and types of management depending on the distance to populations of European minks: first, eliminating the current American mink populations from the overlapping and peripheral areas of European mink distribution; and second, reducing their populations in the other areas of the country to minimize their colonizing potential.However, the Spanish strategy only considers the current distribution of both species, while this study considers both current and future distributions and connectivity.With integration of these future projections, conservation and control efforts can incorporate more proactive strategies.
Considering the overlap between management zones and protected areas also gives valuable information for management of both species as conservation plans to counteract climate change and invasive species effects can be implemented more easily there.The protected areas that overlap with American mink overall potential habitat can aid in control efforts by developing special programs to prevent species arrival, monitor their presence, remove individuals, or lessen the impact of their populations depending on the state of their distribution.Lastly, to enhance conservation efforts, the priority areas to preserve European minks could be evaluated for inclusion in the net of protected areas to avoid their exposure to land-use changes or anthropogenic pressure that would further threaten the fragile status of the species.Although future projections of land use change are not currently available for Spain at a scale relevant to our analyses, if produced in the future, this type of information could be overlaid and evaluated in the context of the management zones to determine where land-use decisions could affect European mink conservation.Such landscape scale decisions are particularly important given that 80% of European minks' potential habitat lies outside protected areas.The overall potential habitat and the generalized importance (Fig. 5) of the species should be considered in future protected area delineation to preserve highly adequate (suitable and connected), irreplaceable, and persistent areas for the European minks.Dynamic connectivity analyses to inform management of the invasive American mink and its native… 1 3 Vol.: (0123456789)

Conclusions
Climate change is likely to exacerbate the negative effects of invasive species while concurrently posing a threat to species of conservation concern.Habitat suitability and dynamic connectivity approaches can be useful tools to forecast species distributions and movements to anticipate management measures within the context of climate change and support adaptive planning that accounts for a changing landscape.Specifically, these tools can improve our understanding of how the landscape facilitates the functionality of ecological networks over time.In this study, climate change is expected to negatively affect American minks and European minks by 2100, with poorer outcomes for the native mink.However, results from the different emissions scenarios highlight the uncertainty surrounding habitat suitability projections.This uncertainty can be addressed by integrating multiple statistical techniques, periods, and emission scenarios to identify a range of potential outcomes for the species of interest and develop flexible and responsive management measures to adapt management strategies to new knowledge.Lastly, habitat suitability and connectivity models can guide both conservation and control measures, and integrating the results of invasive and threatened species to differentiate management zones can support the development of more effective management measures.Although our study was focused on the case of the American and European mink in Spain, our analytical framework can readily be applied to similar studies in different locales or with different focal species.

Fig. 2
Fig. 2 Diagram showing the workflow of the habitat suitability maps by combining A climate and B landscape suitability for each species.Climate suitability (A) was calculated by assembling several individual suitability models: Generalized Linear Model (GLM), generalized additive models (GAM), generalized boosting models (GBM), maximum entropy (MAXENT), multivariate adaptive regression splines (MARS),

Fig. 3
Fig. 3 Connectivity analysis workflow.Asterisks denote 10 different scenarios: one historic and three future scenarios, each of the future ones with three emissions scenarios

Fig. 5
Fig. 5 Generalized importance of habitat patches for A American mink and B European mink.The generalized patch importance was calculated as the sum of min-max normalized patch

Fig. 6
Fig. 6 Importance of areas to A control American minks and B conserve European mink in terms of connectivity separated in different management areas: overlap zone (potential habitat for both species), expansion zone (potential American mink habitat that within a buffer of 40 km of potential European