# Identifying refugia and corridors under climate change conditions for the Sichuan snub-nosed monkey (Rhinopithecus roxellana) in Hubei Province, China

Yu Zhang, Céline Clauzel, Jia Li, Yadong Xue, Yuguang Zhang, Gongsheng

Wu, Patrick Giraudoux, Li Li, Diqiang Li

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

► To cite this version:

Yu Zhang, Céline Clauzel, Jia Li, Yadong Xue, Yuguang Zhang, et al.. Identifying refugia and corridors under climate change conditions for the Sichuan snub-nosed monkey (Rhinopithecus roxellana) in Hubei Province, China. Ecology and Evolution, 2019, 9 (4), pp.1680-1690. 10.1002/ece3.4815. hal-02013080

HAL Id: hal-02013080

https://hal.science/hal-02013080

Submitted on 24 May 2023

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

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

<figure><img src="image_2.png" /></figure>
# Identifying refugia and corridors under climate change conditions for the Sichuan snub-nosed monkey (*Rhinopithecus roxellana*) in Hubei Province, China

Yu Zhang<sup>1,2 ID</sup> | Céline Clauzel<sup>2,3</sup> | Jia Li<sup>1 ID</sup> | Yadong Xue<sup>1</sup> | Yuguang Zhang<sup>1</sup> |
Gongsheng Wu<sup>2,4</sup> | Patrick Giraudoux<sup>2,5</sup> | Li Li<sup>2,4</sup> | Diqiang Li<sup>1</sup>

## Correspondence

Céline Clauzel and Li Li, Key Lab of Hazard Risk Management and Wildlife Management and Ecosystem Health, Yunnan University of Finance and Economics, Kunming, China.
Email: celine.clauzel@univ-paris-diderot.fr; lilyzsu@126.com
and
Diqiang Li, Chinese Academy of Forestry/ Key Laboratory of Biodiversity of National Forestry and Grassland Administration, Research Institute of Forest Ecology, Environment and Protection, Beijing, China.
Email: lidiqiang_caf@163.com

## Funding information

National Key Technology R&D Program of China, Grant/Award Number: 2013BAD03B02 and 2013BAD03B03

## Abstract

Using a case study of an isolated management unit of Sichuan snub-nosed monkey (*Rhinopithecus roxellana*), we assess the extent that climate change will impact the species' habitat distribution in the current period and projected into the 2050s. We identify refugia that could maintain the population under climate change and determine dispersal paths for movement of the population to future suitable habitats. Hubei Province, China. We identified climate refugia and potential movements by integrating bioclimatic models with circuit theory and least-cost model for the current period (1960–1990) and the 2050s (2041–2060). We coupled a maximum entropy algorithm to predict suitable habitat for the current and projected future periods. Suitable habitat areas that were identified during both time periods and that also satisfied home range and dispersal distance conditions were delineated as refugia. We mapped potential movements measured as current flow and linked current and future habitats using least-cost corridors. Our results indicate up to 1,119 km<fnref n="2" /> of currently suitable habitat within the study range. Based on our projections, a habitat loss of 67.2% due to climate change may occur by the 2050s, resulting in a reduced suitable habitat area of 406 km<fnref n="2" /> and very little new habitat. The refugia areas amounted to 286 km<fnref n="2" /> and were located in Shennongjia National Park and Badong Natural Reserve. Several connecting corridors between the current and future habitats, which are important for potential movements, were identified. Our assessment of the species predicted a trajectory of habitat loss following anticipated future climate change. We believe conservation efforts should focus on refugia and corridors when planning for future species management. This study will assist conservationists in determining high-priority regions for effective maintenance of the endangered population under climate change and will encourage increased habitat connectivity.

## KEYWORDS

circuit theory, climate change adaptation, connectivity, golden snub-nosed monkey, habitat suitability, refuge

<footnote n="1"><sup>1</sup>Chinese Academy of Forestry/
Key Laboratory of Biodiversity of
National Forestry and Grassland
Administration, Research Institute of Forest
Ecology, Environment and Protection,
Beijing, China</footnote>

<footnote n="2"><sup>2</sup>Key Lab of Hazard Risk Management and Wildlife Management and Ecosystem Health, Yunnan University of Finance and Economics, Kunming, China</footnote>

<footnote n="3"><sup>3</sup>LADYSS, UMR7533-CNRS, University Paris Diderot, Sorbonne Paris Cité, Paris, France</footnote>

<footnote n="4"><sup>4</sup>School of Urban Management and Resource Environment, Yunnan University of Finance and Economics, Kunming, China</footnote>

<footnote n="5"><sup>5</sup>Chrono-Environnement, UMR 6249 CNRS, University of Bourgogne Franche-Comté, Besançon, France</footnote>
# 1 | INTRODUCTION

Outcomes of climate change challenge the conservation and persistence of many species (Gouveia et al., 2016), especially those considered to be endemic and threatened with higher extinction risks (Lambers, 2015). Climate change can lead to changes in characteristics of movement and species distribution (Lister, Brocki, & Ament, 2015), including range shifts, habitat contractions and expansions and fragmentation (Parmesan, 2006; Struebig et al., 2015). Projections of spatial variations in species distributions and the identification of regions where species could persist and later expand during climate change are essential for targeted conservation efforts (Lambers, 2015; Struebig et al., 2015). Generally, refugia have relatively stable climates and thus facilitate the survival of many taxa under quickly changing environmental conditions (Ashcroft, 2010; Keppel et al., 2015; Li et al., 2016).

The identification of refugia relies on recognition of the current and future distributions of species (Keppel et al., 2015). Habitat suitability assessments based on species distribution modeling (SDM) have been widely used to understand the species response to environmental changes (Gouveia et al., 2016) and to identify refugia (Keppel et al., 2012). Alternatively, climate change velocity can also be used to identify refugia, as these areas characterized by low climate change velocity (Keppel et al., 2015; Sandel et al., 2011). In addition, potential refugia may be recognized by investigating analogues set to current climatic conditions and disturbance distributions (Keppel et al., 2015). However, determining if a species can adjust to shifting climatic conditions by altering its current distribution to locations of climatic refugia or optimal habitat is challenging (Lambers, 2015) because the accessibility of analogue climate habitats is constrained by landscape permeability from anthropogenic modifications and species dispersal capabilities (Littlefield, McRae, Michalak, Lawler, & Carroll, 2017). Thus, maintaining habitat connectivity, the degree to which organisms can move through a landscape (Taylor, Fahrig, Henein, & Merriam, 1993), is the most commonly recommended strategy for conserving species diversity in a changing climate (Heller & Zavaleta, 2009; Nuñez et al., 2013).

Increasing connectivity is critical for species' ability to track rapidly changing climates and for reconfiguring habitats to facilitate access to more suitable habitats (McGuire, Lawler, McRae, Nuñez, & Theobald, 2016). Previous studies of connectivity analysis are used to identify areas that promote species movements between their current habitats (Brodie et al., 2015; Larue & Nielsen, 2008; Wang et al., 2014). However, this analysis does not consider that suitable areas will shift with climate change; thus, it does not determine key areas of importance under changing climate conditions. Although some papers have addressed the shifting of habitats by modeling connectivity for climate change, for example, Nuñez et al., (2013), Littlefield et al. (2017) and Brost and Beier (2012), the majority of these studies have not considered specific species, but instead addressed “coarse-filter” conservation that ignores species-specific traits. Therefore, assessments of the impacts of climate change on

connectivity and distribution of specific species are needed, particularly those with high conservation values (Schmitz et al., 2015).

The Sichuan snub-nosed monkey (*Rhinopithecus roxellana*) is a primate species endemic to China (Li, 2004) and is listed as Endangered by the International Union for Conservation of Nature (Long & Richardson, 2008). *R. roxellana* currently occurs in three isolated temperate montane forest regions in China with a total of approximately 22,000 individuals, comprised of approximately 16,000 individuals in the Sichuan-Gansu population, approximately 5,500 individuals in the Shanxi population, and approximately 1,000 individuals in the Hubei population (Chang, Luo, et al., 2012). *R. roxellana* has suffered a population decline of more than 50% over the past half-century (Li, Pan, & Oxnard, 2002). Currently, the major threats to its population are tourism-related activities and continued habitat loss (Long & Richardson, 2008). Hubei Province is at the easternmost edge of the species distribution and harbors over 1,000 individuals (Liu et al., 2015). This population is considered as a stand-alone management unit (Chang, Luo, et al., 2012), determined to be a group of conspecific individuals among which the degree of connectivity is sufficiently low and thus should be monitored and managed separately (Palsbøll, Bérubé, & Allendorf, 2007; Taylor & Dizon, 1999). The relatively lower genetic diversity, genetically distinct status and small population size of this population makes it more vulnerable to environmental change than the other two populations (Li et al., 2007; Luo, Pan, Liu, & Li, 2012; Pan et al., 2009). This arboreal species lives in temperate broadleaf and coniferous forests (Chang, Liu, Yang, Li, & Vigilant, 2012), and the impact of climate change on vegetation is expected to decrease the availability of its suitable habitat (Luo et al., 2015; Xiang et al., 2011). These predicted range reductions will force the Hubei population monkeys to migrate to higher elevations over time (Lou et al., 2015). In addition, there is little research identifying the refugia of *R. roxellana* and how it may disperse from its currently suitable habitats to future habitats. Identifying the extent to which the *R. roxellana* could be affected by climate change and the dispersal paths between its current and future habitats are therefore important for effective conservation management.

Here, we built a *R. roxellana* distribution model based on monkey occurrences associated with bioclimatic and environmental variables. We aimed to (a) assess the extent that climate change will impact the *R. roxellana* habitat distribution from the current period to the 2050s, (b) identify refugia locations that could maintain the population under climate change and dispersal paths that could facilitate the movement of the population to suitable future habitats, and (c) provide suggestions for the conservation of *R. roxellana* in Hubei Province under environmental change.

# 2 | METHODS

## 2.1 | Study area and population of *R. roxellana*

This study was conducted at the main distribution area of *R. roxellana* in Hubei Province (Figure 1). The study area was expanded appropriately to cover larger areas in the surrounding region, including the
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** Map of the study area in Hubei Province and the occurrences and locations of subpopulations of *Rhinopithecus roxellana* (DLT, JHL, and QJP refer to the Dalongtan, Jinhoulin, and Qianjiaping subpopulations, respectively). The population in Hubei Province was found in the Shennongjia (SNJ) National Park, Shennongjia Forestry District, and Badong Natural Reserve</figcaption></figure>

entire Shennongjia (SNJ) Forestry District, the Badong Nature Reserve, and other areas within the rectangular range in an effort to prevent omissions of *R. roxellana* likely habitat. The SNJ National Park is located in the SNJ Forestry District; however, this paper uses the term “SNJ Forestry District” to refer to the area not in the SNJ National Park. This paper mainly involves three areas, that is, The SNJ Forestry District, SNJ National Park, and Badong Nature Reserve (Figure 1). There are three monkey subpopulations, Dalongtan (DLT), Jinhoulin (JHL), and Qianjiaping (QJP), that inhabit the SNJ National Park (Yang, Liao, Yu, & Yao, 2008). A substantial area of the SNJ National Park suffered deforestation between the 1950s and the early 1980s (Zhu, 1992). In recent years, tourism at the SNJ has developed rapidly (Xiang et al., 2011), resulting in human disturbances near scenic locations and habitat fragmentation caused by the increased traffic. Human disturbances in the SNJ Forestry District are mainly concentrated in the spring, when villagers collect Chinese medicinal herbs.

## 2.2 | Occurrence records

The occurrence records of *R. roxellana* (N = 1,225) comprised data from standard monitoring activities by the staff of the SNJ National

Park from May to December 2013, camera-trap monitoring data and the published literature (Luo et al., 2015; Su et al., 2004). To diminish the impact of spatial autocorrelation, we filtered multiple records by randomly selecting one record in a 1 km × 1 km grid cell. A total of 234 occurrence records were used in this study (Figure 1).

## 2.3 | Environmental variables

We obtained a set of climatic variables from the WorldClim database (http://www.worldclim.org) at a 30" resolution to define the species distribution models for current climate (average for 1960–1990) and future climatic (average for 2041–2060) scenarios. The projection was based on the general circulation models of HadGEM₂-AO of the Fifth Assessment Report of the International Panel of Climate Change. The simulation used HadGEM₂-AO which performed well in the East Asia region (Baek et al., 2013). The scenario was described in terms of the Representative Concentration Pathways 4.5, which predict an average global temperature increase of 0.9–2.0°C by the 2050s (UNFCCC, 2015). The time horizon of the 2050s was selected because it is a period far enough in the future for significant climate changes to have occurred (Young et al., 2012).
Other environmental variables (Supporting Information Table S1) used to construct the *R. roxellana* distribution models included density of rivers, roads and settlements and vegetation types, which were obtained from a 1:1,000,000 map of China (National Geomatics Center of China). The elevation was derived from a digital elevation model with a resolution of 30" from the WorldClim database. Non-climate variables are not available for the 2050s, and vegetation changes occur slowly, so we kept these variables static in our projections (Li, Liu, Xue, Zhang, & Li, 2017).

All variables were resampled at a resolution of 1 km × 1 km and put into the same projection using ArcGIS 10.1 (ESRI Inc., Redlands, CA, USA). To reduce multicollinearity, a reduced version was produced by eliminating Pearson's correlation coefficients of $|r| > 0.8$ (Cord, Klein, Mora, & Dech, 2014; Supporting Information Table S2), leaving seven variables to construct the *R. roxellana* distribution model (Temperature Seasonality (Bio4), Min Temperature of Coldest Month (Bio6), Precipitation Seasonality (Bio15), Settlement density, Road density, River density, and Vegetation type; Table S1).

## 2.4 | Species distribution model

We employed the maximum entropy algorithm (MaxEnt 3.3.3k), one of the best performing approaches in modeling species distribution, with presence-only data (Elith & Yates, 2015) to construct the habitat suitability for *R. roxellana* under the current day scenario and then projected the spatial information into the 2050s. We used the default settings for the MaxEnt model (Phillips, Anderson, & Schapire, 2006), with the exception of dividing the percentage at 75% of the occurrence data into training set for model construction and the percentage at 25% into a random test set for evaluation of model performance. We conducted a subsampling procedure with 15 replicates (Khatchikian, Sangermano, Kendell, & Livdahl, 2011).

To evaluate the model performance, we used the threshold-independent area under the receiver operating characteristic curve (AUC) with value ranges from 0 to 1. AUC values close to 1 indicate perfect model agreement (Phillips et al., 2006). Variable importance was estimated by the permutation importance method (Searcy & Shaffer, 2016). The logistic results of the MaxEnt model were considered to represent the probabilities of species occurrence (Phillips & Dudik, 2008). We then reclassified the results into presence and absence using the average of the maximum training sensitivity plus specificity (Liu, Berry, Dawson, & Pearson, 2005; Songer, Delion, Biggs, & Huang, 2012). Areas with probability values above the threshold were regarded as suitable habitats.

## 2.5 | Assessing habitat vulnerability and climate refugia

We assessed the impact of climate change on *R. roxellana* based on the suitable habitat changes between the current period and the 2050s. The aim was to identify vulnerable habitats, that is, current suitable habitat that will be lost by 2050, and climatic refugia, that is, areas where suitable habitat was present in the current period and

in the 2050s projection. The home range of the species is 18.3 km², of which 7.4 km² is considered core area (Tan, Guo, & Li, 2007). The daily path length varies from 0.75 to 5 km with a mean of 2.1 km (Tan et al., 2007). *R. roxellana* has a male-based dispersal system (Chang et al., 2014; Huang et al., 2017). We parameterized our models of refugia according to those values; for example, patch areas had to be >7.4 km², and the distance to the nearest patch had to be <2.1 km and located in a protected area for more effective management and protection. Three indicators were employed to assess the habitat vulnerability, including the percentage of suitable habitat area change (AC), the percentage of currently suitable area that was lost by the 2050s ($SH_L$), and the percent increase of future suitable area by the 2050s ($SH_I$). Indicators were calculated as follows:

$$AC = \frac{A_F - A_C}{A_C} \times 100\%$$

$$SH_L = \frac{A_C - A_{FC}}{A_C} \times 100\%$$

$$SH_I = \frac{A_F - A_{FC}}{A_C} \times 100\%$$

where $A_F$ is the area of the projected suitable habitat for *R. roxellana* under the 2050s climatic scenario; $A_C$ is the area of the modeled current suitable habitat; $A_{FC}$ is the constant area of the suitable habitat in both the current period and the 2050s (Irina, Flemming, Jenschristian, & Carsten, 2007; Li et al., 2017).

The difference in the average elevation of the suitable habitat between the time periods was examined using a Mann-Whitney U test. Statistical analyses were carried out using SPSS 19.0 software.

## 2.6 | Habitat connectivity analysis

We used Circuitscape v4.0 (McRae, Shah, & Mohapatra, 2013) to quantify potential species movement between suitable habitat patches in the current period and the 2050s. The Circuitscape model connectivity was based on the circuit theory, predicting the movement patterns of random walkers between source and target cells across a landscape. High densities of current flow indicate important movement between habitat patches (McRae, Dickson, Keitt, & Shah, 2008). We ran Circuitscape using the pairwise mode for the current period and the 2050s. The suitable habitat patches in each period were treated as nodes (the source and target). We used the inverse of the logistic output from our MaxEnt model as a measure of movement resistance for *R. roxellana* and rescaled from 1 to 100 to construct a resistance layer (Li et al., 2016). For cells with a value larger than the threshold (maximum training sensitivity plus specificity was given by the MaxEnt model), resistance was set to 1. For cells with a value smaller than
the threshold, resistance was set to (threshold - "value") × 100/
threshold. We mapped the current flow between suitable habitats
for two states (current and the future climate scenario) to visual-
ize habitat connectivity.

We linked suitable habitats in the 2050s scenario to model po-
tential dispersal paths from current to future habitats. We mapped
the least-cost path and least-cost corridors between suitable hab-
itats whose distance to the nearest patch was <2.1 km. The resis-
tance layer identical to the input used in Circuitscape. The least-cost
path and corridors were implemented in LinkageMapper (McRae &
Kavanagh, 2011), and truncated at a cost distance of 200,000 cost
unit for visualization.

3 | RESULTS

3.1 | Habitat suitability model

The MaxEnt model for *R. roxellana* provided satisfactory results, with an AUC value of 0.955 (±0.005). Temperature seasonality (bio04, 57.8%, Supporting Information Table S1) contributed the most to the model, according to the permutation importance, followed by the precipitation seasonality (bio15, 21.9%), settlement density (10.4%), vegetation type (4.6%), road density (4.1%), river density (0.7%), and min temperature of the coldest month (bio06, 0.6%). The average threshold value for the measure of suitable habitat was 0.185.

3.2 | Habitat suitability for the current period and the 2050s

For the current period, the area with a habitat suitability value
higher than 0.185 was 1,119 km², predominately concentrated
in SNJ National Park (722 km²), the southwest region of the SNJ
Forestry District (247 km²), and the northern region of the Badong
Nature Reserve (61 km²) (Figure 2a). In the 2050s, the suitable
habitat area was dramatically reduced to 406 km², representing
a decrease of -63.7%. The most suitable habitat retreated to SNJ
National Park (293 km², AC = -59.4%) and the northern region of
Badong Nature Reserve (60 km², AC = -1.6%). The SNJ Forestry
District lost the majority of its suitable habitat, with only 11 km² of
existing suitable habitat area remaining in the 2050s (AC = -95.5%)
(Table 1, Figure 2b).

Suitable habitat scarcely increased in the 2050s (SH<sub>I</sub> = 3.5%, Table 1, Supporting Information Figure S1). Most of the habitat was shown to be vulnerable to climate change (SH<sub>L</sub> = 67.2%). Among the main areas, the habitat of the Badong Nature Reserve had the lowest vulnerability to climate change (SH<sub>L</sub> = 9.8%), followed by the habitat of the SNJ National Park (SH<sub>L</sub> = 62.2%) and the SNJ Forestry District (SH<sub>L</sub> = 96.0%).

Climate change would result in a shift of the suitable habitat of *R. roxellana* to higher elevations (Figure 3). The mean elevation of the suitable habitat in the 2050s was projected to be 2,183.66 ± 325.35 m, which is significantly higher (Z = -3.554, *p* = 0.000) than the elevation of the currently suitable habitat (2,118.90 ± 315.94 m).

<figure><img src="image_4.png" /><figcaption>**FIGURE 2** Habitat suitability for *Rhinopithecus roxellana* in Hubei Province. (a) depicts the models under the current climate. (b) is the projection of habitat suitability by the 2050s</figcaption></figure>

3.3 | Climate refugia and potential movement under climate change

By intersecting areas suitable for *R. roxellana* in the current and future scenarios and filtering areas according to species-specific parameters, a total area of 286 km² was identified as climate refugia (Figure 4b). In accordance with the parameters, the patch areas had to be >7.4 km², and the distance to the nearest patch had to be <2.1 km and located in a protected area.

The model of circuit theory highlighted "pinch points" of high
level movement. In the current situation, several areas exhibited
potentially high current flow; therefore, they likely indicate critical
pathways for possible movement across suitable habitats (Figure 5a).
As habitats shrink and narrow under the conditions modeled for the
2050s, movement will be facilitated by narrowed habitats and se-
verely restricted (Figure 5b). From the least-cost path and least-cost
corridors, we modeled linkages between current and future habitats
<table><caption>**TABLE 1** Estimates of the suitable habitat area (km²) for *Rhinopithecus roxellana* in Hubei Province in the current period and the 2050s, with the percentage change values</caption><thead><tr><th></th><th>A<sub>C</sub> (km²)</th><th>A<sub>F</sub> (km²)</th><th>A<sub>FC</sub> (km²)</th><th>AC (%)</th><th>SH<sub>L</sub> (%)</th><th>SH<sub>I</sub> (%)</th></tr></thead><tbody><tr><td>Study area</td><td>1,119</td><td>406</td><td>367</td><td>-63.7</td><td>67.2</td><td>3.5</td></tr><tr><td>SNJ Forestry District</td><td>247</td><td>11</td><td>10</td><td>-95.5</td><td>96.0</td><td>0.4</td></tr><tr><td>SNJ National Park</td><td>722</td><td>293</td><td>273</td><td>-59.4</td><td>62.2</td><td>2.8</td></tr><tr><td>Badong Nature Reserve</td><td>61</td><td>60</td><td>54</td><td>-1.6</td><td>11.5</td><td>9.8</td></tr></tbody></table>

Notes. A<sub>C</sub> is the area of the modeled current suitable habitat, A<sub>F</sub> is the area of the projected suitable habitat under the climate scenario of the 2050s, and A<sub>FC</sub> is the constant area of suitable habitat both in the current and the future (2050s) periods. AC, SH<sub>L</sub>, and SH<sub>I</sub> refer to the percentage of suitable habitat area change, currently suitable area that was lost by 2050s and increased in future suitable area by the 2050s, respectively.

(Figure 4a). These linkages emphasized routeways in which *R. roxellana* can navigate from current to future habitats.

## 4 | DISCUSSION

In this paper, we assessed the impact of climate change on the *R. roxellana* habitat range and identified refugia and corridors under climate change conditions. Our assessment of the species predicted a trajectory of habitat loss following anticipated future climate change. Additionally, we identified strategic areas that should be prioritized during species preservation efforts. Our analysis provides a perspective to evaluate the impact of climate change on habitat connectivity for a species of conservation concern, and proposes targeted actions that consider species characteristics.

### 4.1 | Projected effects of climate change

There have been few published articles about the effects of the changing climate on the *Rhinopithecus* genus to date. A case study about *Rhinopithecus bieti*, which used climate and socio-economic scenarios to model land cover changes until 2050, predicted that the area of suitable habitats will be reduced by 15% by 2050 (Li, Wu, Xue, He, & Giraudoux, 2011). The potential distribution area of *Rhinopithecus avunculus* will be 20% smaller in 2020 than it is under current climate conditions (Van, Manh, & Hoang, 2010). According to our results, *R. roxellana* in Hubei will be the most seriously impacted species of the *Rhinopithecus* genus, with a habitat decrease of 67.2% and a meager 3.5% addition of newly formed habitat by the 2050s. This suggests that the effects of climate change will greatly exacerbate the vulnerability of *R. roxellana*. These results are consistent with the research of Luo et al. (2015), which predicted a habitat reduction of almost 70% by 2050 in the SNJ area and an increase in the mean elevation of suitable habitat. Compared with Luo's analysis, which omitted the Badong Nature Reserve, we considered the entire distribution area of the Hubei population. Among the three subpopulations, the DLT and JHL subpopulations may be forced to recolonize southward due to substantial habitat loss, likely resulting in intraspecific resource competition that will cause population decline of *R. roxellana*.

### 4.2 | Refugia and potential movement under climate change

According to the ecological requirements of the species, we further refined the area of refugia that can sustain subpopulations and the corridors that can serve as priority sanctuaries. Such refugia were primarily located in the SNJ National Park and the Badong Natural Reserve, which suggests the high potential of these regions for use in preserving the population under climate change (Figure 4). Preservation efforts between the two protected areas should focus on building trans-boundary cooperation systems to ensure habitat connectivity.

Under climate change conditions, *R. roxellana* will be forced to shift their range away from the current habitat areas that will no longer be suitable and arrive at future habitats. The dispersal path of the least-cost model can improve movement into the 2050s' optimal habitat (Figure 4a). Moreover, the map connecting the suitable habitats under the current and future scenarios using circuit theory highlights different corridor types that are important for potential movements and emphasizes that corridors are likely

<figure><img src="image_5.png" /><figcaption>**FIGURE 3** Impacts of climate change on the potential distribution of *Rhinopithecus roxellana* in Hubei Province over elevation. The y-axis represents the frequencies of suitable habitat cell numbers, with 100 m elevation intervals, for the current (light gray) period and the 2050s (dark gray)</figcaption></figure>
<figure><img src="image_6.png" /></figure>

**FIGURE 4** (a) The least-cost path and corridors between current and future suitable habitats, least-cost corridors truncated at cost distance of 200,000 cost units. (b) Climate refugia and corridors where priority conservation efforts could be focused. The potential corridor areas numbered 1, 2, and 3 are the more critical areas for current movements and tracking climate change. The corridor areas numbered 4, 5, 6, and 7 are more critical under future conditions

to vary in importance with time (Figure 5). The potential corridor areas numbered one, two, and three were the more vital areas for current movements and tracking climate change. Corridors four, five, six, and seven will become important under future conditions. Thus, a staged process of conservation actions could be implemented as follows: Conservation actions could be immediately focused on corridor areas 1–3 through assessment of habitat quality and restoration of degraded habitats; the future corridor areas 4–7 should also be protected to ensure they remain ecologically intact for their future use. Rehabilitation of unsuitable habitats could be taken into consideration around these areas. Priority conservation efforts should be focused on these strategic areas. Likewise, long-term efforts are also needed to monitor not only habitat quality, including its phenology, food availability, and plant community dynamics, but also the status and movement trends of each subpopulation.

## 4.3 | Implication for species conservation

Because climate change is anticipated to alter both the movements and distributions of species, there is a need to take climate change into account during the promotion of adaptive conservation strategies (Schmitz et al., 2015). Investigating connectivity without considering shifts in future habitats would omit areas that are crucial for species movement from current to future habitats. Understanding the connection from current habitats to future habitats can allow for prediction of significant areas necessary for the climate-induced movement and where interventions for conservation efforts could be focused. For example, Littlefield et al. (2017) identified key areas likely to facilitate climate-induced species movement across western North America. For habitat reduction at the small study scale, like that investigated in this study, the changes in suitable habitat locations also influence the pathways that animals need to track climate change and
<figure><img src="image_7.png" /><figcaption>**FIGURE 5** Potential *Rhinopithecus roxellana* movements in Hubei Province, based on circuit theory, between suitable habitats of (a) the current climate and (b) the climate scenario of the 2050s. The color ramp reflects the absolute potential movement values for each situation. Thus, colors are not directly comparable, and the results should be interpreted regarding the relative importance of the potential movement area</figcaption></figure>

reach emerging suitable habitats. This also highlights the importance of protecting the current habitat for such situations. The prediction of the potential movement between future habitats will allow for advanced preparation against the negative impacts of climate change.

Connectivity analysis that considers climate change typically focuses on two scales—regional and individual species. The regional scale takes several species into account simultaneously; therefore, the ecological niche needs of each species, such as specific climatic condition, elevations, slopes, and special vegetation types, cannot be considered in detail. Instead, it uses a coarsely approximate climatic condition to determine current habitat distribution and to predict future conditions. For example, Littlefield et al. (2017) used the multivariate similarity of climates through a principal component analysis to identify climate analogue habitats. The analysis based on the individual species scale, as a refinement of the analysis based on the regional scale, usually requires a species distribution model to simulate future habitats, which can focus on specific requirements. For example, Kang, Minor, Lee, and Park (2016) used a species-based MaxEnt model to valuate changes in the extent and connectivity of castor aralia (*Kalopanax septemlobus*) habitat, taking bioclimatic and topographic variables into account. The individual species scale analysis can make a complete corridor prediction involving species dispersal capability.

Identifying resistance surface is challenging in connectivity analysis (Milanesi et al., 2017). Ideally, resistance values would be parameterized with empirical data, but because of a shortage of such information, expert knowledge is often used (Stevenson-Holt, Watts, Bellamy, Nevin, & Ramsey, 2014). However, the use of expert opinion is seen as subjective, human-centric and unreliable (Puyravaud, Cushman, Davidar, & Madappa, 2017). Thus, the resistance surface converted from the habitat suitability index is superior to the subjective method of expert opinion to some extent. It takes into account a variety of variables affecting the distribution of species, and facilitates the rapid assessment of the connectivity of specific species.

## 4.4 | Contributions and limits of the method

Our model identified climate-induced shifts in the habitat distribution of *R. roxellana* and determined the most effective regions enabling *R. roxellana* to shift its distribution range. However, we made several simplifying assumptions. First, our approach used a relatively coarse resolution of climate data and did not take fine-scale characteristics into account, which might be central to long-term survival (Struebig et al., 2015) and provide climate microrefugia for *R. roxellana*. We also used a static map of vegetation and human disturbances, because changes in the vegetation community lag
in response to climate change (Barbet-Massin, Thuiller, & Jiguet, 2012), and the main distribution area was in protected areas where construction projects need strict approval. However, we do not know precisely what kind of changes in land use under climate change might occur with an increasing human population and global economy shifts. Here, we assume that the area will be relatively protected. In addition, our analysis was a species-specific approach that cannot be directly extended to other species, although the preservation of *R. roxellana* could benefit other symmetric animal species. Nonetheless, we note that *R. roxellana* in Hubei Province is a stand-alone management unit (Chang, Luo, et al., 2012), and that use a relatively small spatial extent was valuable for comparing the effects of climate change.

The identification of refugia areas and analysis of climate connectivity will allow conservationists to determine the most effective regions for maintaining a population of *R. roxellana* and increasing habitat connectivity, in the context of the climate changes that are predicted for this century.

## ACKNOWLEDGEMENTS

This research was supported by the National Key Technology R&D Program of China (2013BAD03B02 and 2013BAD03B03). We thank all the staff of the Shennongjia National Park for helping us collect the occurrence data of the Sichuan snub-nosed monkey. The authors would like to thank Charlotte Hacker for assistance in editing this paper. This study benefited from inputs of the GDRI EHEDE network (http://gdri-ehede.univ-fcomte.fr).

## CONFLICT OF INTEREST

None declared.

## AUTHOR CONTRIBUTIONS

Yu Zhang and Diqiang Li conceived and designed the work. Yadong Xue, Yuguang Zhang, and Gongsheng Wu collected the data. Yu Zhang, Jia Li analyzed and interpreted the data. Yu Zhang drafted the article. Celine Clauzel, Patrick Giraudoux, and Li Li critically revised the article. Yu Zhang, Céline Clauzel, Jia Li, Yadong Xue, Yuguang Zhang, Gongsheng Wu, Patrick Giraudoux, Li Li, and Diqiang Li gave final approval of the version to be published.

## DATA ACCESSIBILITY

Climate data and MaxEnt input files for this study will be available at the Dryad after the paper publishing.

## ORCID

Yu Zhang [iD](https://orcid.org/0000-0003-1676-4147)

Jia Li [iD](https://orcid.org/0000-0002-0794-4376)

## REFERENCES

Ashcroft, M. B. (2010). Identifying refugia from climate change. *Journal of Biogeography*, *37*, 1407–1413. https://doi.org/10.1111/j.1365-2699.2010.02300.x

Baek, H. J., Lee, J., Lee, H. S., Hyun, Y. K., Cho, C. H., Kwon, W. T., ... Byun, Y.-H. (2013). Climate change in the 21st century simulated by HadGEM2-AO under representative concentration pathways. *Asia-Pacific Journal of Atmospheric Sciences*, *49*, 603–618. https://doi.org/10.1007/s13143-013-0053-7

Barbet-Massin, M., Thuiller, W., & Jiguet, F. (2012). The fate of European breeding birds under climate, land-use and dispersal scenarios. *Global Change Biology*, *18*, 881–890. https://doi.org/10.1111/j.1365-2486.2011.02552.x

Brodie, J. F., Giordano, A. J., Dickson, B., Hebblewhite, M., Bernard, H., Mohd-Azlan, J., ... Ambu, L. (2015). Evaluating multispecies landscape connectivity in a threatened tropical mammal community. *Conservation Biology*, *29*, 122–132. https://doi.org/10.1111/cobi.12337

Brost, B. M., & Beier, P. (2012). Use of land facets to design linkages for climate change. *Ecological Applications*, *22*, 87–103. https://doi.org/10.1890/11-0213.1

Chang, Z., Liu, Z., Yang, J., Li, M., & Vigilant, L. (2012). Noninvasive genetic assessment of the population trend and sex ratio of the Shennongjia population of Sichuan snub-nosed monkeys (*Rhinopithecus roxellana*). *Chinese Science Bulletin*, *57*, 1135–1141. https://doi.org/10.1007/s11434-011-4872-6

Chang, Z. F., Luo, M. F., Liu, Z. J., Yang, J. Y., Xiang, Z. F., Li, M., & Vigilant, L. (2012). Human influence on the population decline and loss of genetic diversity in a small and isolated population of Sichuan snub-nosed monkeys (*Rhinopithecus roxellana*). *Genetica*, *140*, 105–114. https://doi.org/10.1007/s10709-012-9662-9

Chang, Z., Yang, B., Vigilant, L., Liu, Z., Ren, B., Yang, J., ... Li, M. (2014). Evidence of male-biased dispersal in the endangered Sichuan snub-nosed monkey (*Rhinopithecus roxellana*). *American Journal of Primatology*, *76*, 72–83.

Cord, A. F., Klein, D., Mora, F., & Dech, S. (2014). Comparing the suitability of classified land cover data and remote sensing variables for modeling distribution patterns of plants. *Ecological Modelling*, *272*, 129–140. https://doi.org/10.1016/j.ecolmodel.2013.09.011

Elith, J., & Yates, C. J. (2015). A statistical explanation of MaxEnt for ecologists. *Diversity and Distributions*, *17*, 43–57. https://doi.org/10.1111/j.1472-4642.2015.00725.x

Gouveia, S. F., Souza-Alves, J. P., Rattis, L., Dobrovolski, R., Jerusalinsky, L., Beltrão-Mendes, R., & Ferrari, S. F. (2016). Climate and land use changes will degrade the configuration of the landscape for titi monkeys in eastern Brazil. *Global Change Biology*, *22*, 2003–2012. https://doi.org/10.1111/gcb.13162

Heller, N. E., & Zavaleta, E. S. (2009). Biodiversity management in the face of climate change: A review of 22 years of recommendations. *Biological Conservation*, *142*, 14–32. https://doi.org/10.1016/j.biocon.2008.10.006

Huang, Z. P., Bian, K., Liu, Y., Pan, R. L., Qi, X. G., & Li, B. G. (2017). Male dispersal pattern in Golden snub-nosed monkey (*Rhinopithecus roxellana*) in Qinling Mountains and its conservation implication. *Scientific Reports*, *7*. https://doi.org/10.1038/srep46217

Irina, L., Flemming, S., Jenschristian, S., & Carsten, R. (2007). Potential impacts of climate change on the distributions and diversity patterns of European mammals. *Biodiversity and Conservation*, *16*, 3803–3816. https://doi.org/10.1007/s10531-007-9181-7

Kang, W., Minor, E. S., Lee, D., & Park, C. R. (2016). Predicting impacts of climate change on habitat connectivity of *Kalopanax septemlobus* in South Korea. *Acta Oecologica*, *71*, 31–38. https://doi.org/10.1016/j.actao.2016.01.005

Keppel, G., Mokany, K., Wardell-Johnson, G. W., Phillips, B. L., Welbergen, J. A., & Reside, A. E. (2015). The capacity of refugia for
conservation planning under climate change. *Frontiers in Ecology and the Environment*, 13, 106-112. https://doi.org/10.1890/140055

Keppel, G., Van Niel, K. P., Wardell-Johnson, G. W., Yates, C. J., Byrne, M., Mucina, L., ... Franklin, S. E. (2012). Refugia: Identifying and understanding safe havens for biodiversity under climate change. *Global Ecology and Biogeography*, 21, 393-404. https://doi.org/10.1111/j.1466-8238.2011.00686.x

Khatchikian, C., Sangermano, F., Kendell, D., & Livdahl, T. (2011). Evaluation of species distribution model algorithms for fine-scale container breeding mosquito risk prediction. *Medical and Veterinary Entomology*, 25, 268-275. https://doi.org/10.1111/j.1365-2915.2010.00935.x

Lambers, J. H. R. (2015). Extinction risks from climate change. *Science*, 348, 501-502. https://doi.org/10.1126/science.aab2057

Larue, M. A., & Nielsen, C. K. (2008). Modelling potential dispersal corridors for cougars in midwestern North America using least-cost path methods. *Ecological Modelling*, 212, 372-381. https://doi.org/10.1016/j.ecolmodel.2007.10.036

Li, B., Pan, R., & Oxnard, C. E. (2002). Extinction of Snub-nosed monkeys in China during the past 400 years. *International Journal of Primatology*, 23, 1227-1244.

Li, J., Liu, F., Xue, Y., Zhang, Y., & Li, D. (2017). Assessing vulnerability of giant pandas to climate change in the Qinling Mountains of China. *Ecology and Evolution*, 7, 4003-4015. https://doi.org/10.1002/ece3.2981

Li, J., McCarthy, T. M., Wang, H., Weckworth, B. V., Schaller, G. B., Mishra, C., ... Beissinger, S. R. (2016). Climate refugia of snow leopards in High Asia. *Biological Conservation*, 203, 188-196. https://doi.org/10.1016/j.biocon.2016.09.026

Li, L., Wu, G., Xue, Y., He, Y., & Giraudoux, P. (2011). Habitat and connectivity of the black snub-nosed monkey (*Rhinopithecus bieti*) populations in Yunnan: A climate change scenario. In *Vlth European Congress of Mammalogy*, Paris, France.

Li, M., Liu, Z., Gou, J., Ren, B., Pan, R., Su, Y., ... Wei, F. (2007). Phylogeography and population structure of the golden monkeys (*Rhinopithecus roxellana*): Inferred from mitochondrial DNA sequences. *American Journal of Primatology*, 69, 1195-1209. https://doi.org/10.1002/ajp.20425

Li, Y. (2004). The effect of forest clear-cutting on habitat use in Sichuan snub-nosed monkey (*Rhinopithecus roxellana*) in Shennongjia Nature Reserve, China. *Primates*, 45, 69-72. https://doi.org/10.1007/s10329-003-0059-3

Lister, N. M., Brocki, M., & Ament, R. (2015). Integrated adaptive design for wildlife movement under climate change. *Frontiers in Ecology and the Environment*, 13, 493-502. https://doi.org/10.1890/150080

Littlefield, C. E., McRae, B. H., Michalak, J. L., Lawler, J. J., & Carroll, C. (2017). Connecting today's climates to future climate analogs to facilitate movement of species under climate change. *Conservation Biology*, 31, 1397-1408. https://doi.org/10.1111/cobi.12938

Liu, C., Berry, P. M., Dawson, T. P., & Pearson, R. G. (2005). Selecting thresholds of occurrence in the prediction of species distributions. *Ecography*, 28, 385-393. https://doi.org/10.1111/j.0906-7590.2005.03957.x

Liu, Z., Liu, G., Roos, C., Wang, Z., Xiang, Z., Zhu, P., ... Li, M. (2015). Implications of genetics and current protected areas for conservation of 5 endangered primates in China. *Conservation Biology*, 29, 1508-1517. https://doi.org/10.1111/cobi.12581

Long, Y., & Richardson, M. (2008). *Rhinopithecus roxellana*. The IUCN Red List of Threatened Species 2008: e.T19596A8985735. https://doi.org/10.2305/IUCN.UK.2008.RLTS.T19596A8985735.en

Luo, M., Pan, H., Liu, Z., & Li, M. (2012). Balancing selection and genetic drift at major histocompatibility complex class II genes in isolated populations of golden snub-nosed monkey (*Rhinopithecus roxellana*). *BMC Evolutionary Biology*, 12, 207. https://doi.org/10.1186/1471-2148-12-207

Luo, Z., Zhou, S., Yu, W., Yu, H., Yang, J., Tian, Y., ... Wu, H. (2015). Impacts of climate change on the distribution of Sichuan snub-nosed monkeys

(*Rhinopithecus roxellana*) in Shennongjia area, China: Climate change impacts on snub-nosed monkey. *American Journal of Primatology*, 77, 135-151.

McGuire, J. L., Lawler, J. J., McRae, B. H., Nuñez, T. A., & Theobald, D. M. (2016). Achieving climate connectivity in a fragmented landscape. *Proceedings of the National Academy of Sciences of the United States of America*, 113, 7195-7200. https://doi.org/10.1073/pnas.1602817113

McRae, B. H., Dickson, B. G., Keitt, T. H., & Shah, V. B. (2008). Using circuit theory to model connectivity in ecology, evolution, and conservation. *Ecology*, 89, 2712-2724. https://doi.org/10.1890/07-1861.1

McRae, B. H., & Kavanagh, D. M. (2011). Linkage mapper connectivity analysis software. The Nature Conservancy. Retrieved from http://www.circuitscape.org/linkagemapper.

McRae, B. H., Shah, V. B., & Mohapatra, T. K. (2013). *Circuitscape 4 user guide*. The Nature Conservancy. Retrieved from: http://www.circuitscape.org.

Milanesi, P., Holderegger, R., Caniglia, R., Fabbri, E., Galaverni, M., & Randi, E. (2017). Expert-based versus habitat-suitability models to develop resistance surfaces in landscape genetics. *Oecologia*, 183, 67-79. https://doi.org/10.1007/s00442-016-3751-x

Nuñez, T. A., Lawler, J. J., Mcrae, B. H., Pierce, D. J., Krosby, M. B., Kavanagh, D. M., ... Tewksbury, J. J. (2013). Connectivity planning to address climate change. *Conservation Biology*, 27, 407-416. https://doi.org/10.1111/cobi.12014

Palsbøll, P. J., Bérubé, M., & Allendorf, F. W. (2007). Identification of management units using population genetic data. *Trends in Ecology & Evolution*, 22, 11-16. https://doi.org/10.1016/j.tree.2006.09.003

Pan, D., Hu, H. X., Meng, S. J., Men, Z. M., Fu, Y. X., & Zhang, Y. (2009). A high polymorphism level in *Rhinopithecus roxellana*. *International Journal of Primatology*, 30, 337-351. https://doi.org/10.1007/s10764-009-9347-0

Parmesan, C. (2006). Ecological and evolutionary responses to recent climate change. *Annual Review of Ecology Evolution and Systematics*, 37, 637-669. https://doi.org/10.1146/annurev.ecolsys.37.091305.110100

Phillips, S. J., Anderson, R. P., & Schapire, R. E. (2006). Maximum entropy modeling of species geographic distributions. *Ecological Modelling*, 190, 231-259. https://doi.org/10.1016/j.ecolmodel.2005.03.026

Phillips, S. J., & Dudík, M. (2008). Modeling of species distributions with maxent: New extensions and a comprehensive evaluation. *Ecography*, 31, 161-175. https://doi.org/10.1111/j.0906-7590.2008.5203.x

Puyravaud, J. P., Cushman, S. A., Davidar, P., & Madappa, D. (2017). Predicting landscape connectivity for the Asian elephant in its largest remaining subpopulation. *Animal Conservation*, 20, 225-234. https://doi.org/10.1111/acv.12314

Sandel, B., Arge, L., Dalsgaard, B., Davies, R. G., Gaston, K. J., Sutherland, W. J., & Svenning, J. (2011). The influence of late quaternary climate-change velocity on species endemism. *Science*, 334, 660-664. https://doi.org/10.1126/science.1210173

Schmitz, O. J., Lawler, J. J., Beier, P., Groves, C., Knight, G., Boyce, D. A., ... Trainor, A. M. (2015). Conserving biodiversity: Practical guidance about climate change adaptation approaches in support of land-use planning. *Natural Areas Journal*, 35, 190-203. https://doi.org/10.3375/043.035.0120

Searcy, C. A., & Shaffer, H. B. (2016). Do ecological niche models accurately identify climatic determinants of species ranges? *American Naturalist*, 187, 423-435. https://doi.org/10.1086/685387

Songer, M., Delion, M., Biggs, A., & Huang, Q. (2012). Modeling impacts of climate change on giant panda habitat. *International Journal of Ecology*, 2012, 1-12. https://doi.org/10.1155/2012/108752

Stevenson-Holt, C. D., Watts, K., Bellamy, C. C., Nevin, O. T., & Ramsey, A. D. (2014). Defining landscape resistance values in least-cost connectivity models for the invasive grey squirrel: A comparison of approaches using expert-opinion and habitat suitability modelling. *PLoS ONE*, 9, e112119. https://doi.org/10.1371/journal.pone.0112119
Struebig, M. J., Fischer, M., Gaveau, D. L. A., Meijaard, E., Wich, S. A., Gonner, C., ... Kramer-Schadt, S. (2015). Anticipated climate and land- cover changes reveal refuge areas for Borneo's orang-utans. *Global Change Biology*, *21*, 2891–2904. https://doi.org/10.1111/gcb.12814

Su, H., Ma, Q., Lin, Y., Chen, Q., Qiao, Y., & Tan, M. (2004). Survey on population of the golden monkey (*Rhinopithecus roxellana*) at Badong County, Hubei. *Acta Theriologica Sinica*, *24*, 84–87, 70. (in Chinese).

Tan, C. L., Guo, S., & Li, B. (2007). Population structure and ranging patterns of *Rhinopithecus roxellana* in Zhouzhi National Nature Reserve, Shaanxi, China. *International Journal of Primatology*, *28*, 577–591. https://doi.org/10.1007/s10764-007-9147-3

Taylor, B. L., & Dizon, A. E. (1999). First policy then science: Why a management unit based solely on genetic criteria cannot work. *Molecular Ecology*, *8*, S11–S16. https://doi.org/10.1046/j.1365-294X.1999.00797.x

Taylor, P. D., Fahrig, L., Henein, K., & Merriam, G. (1993). Connectivity is a vital element of landscape structure. *Oikos*, *68*, 571–573. https://doi.org/10.2307/3544927

UNFCCC, (2015). Adoption of the Paris agreement. Report No. FCCC/CP/2015/L.9/Rev.1. Paris, France. Retrieved from http://unfccc.int/resource/docs/2015/cop21/eng/l09r01.pdf.

Van, P. T., Manh, V. V., & Hoang, T. M. (2010). Using environmental niche model to study the distribution of Tonkin snub-nosed monkey (*Rhinopithecus avunculus*) in the Northeastern Vietnam under some climate change scenarios. In 24th International Conference on Informatics for Environmental Protection in cooperation with Intergeo2010: "Integration of Environmental Information in Europe". Bonn, Germany.

Wang, F., Mcshea, W. J., Wang, D., Li, S., Zhao, Q., Wang, H., & Lu, Z. (2014). Evaluating landscape options for corridor restoration between giant panda reserves. *PLoS ONE*, *9*, e105086. https://doi.org/10.1371/journal.pone.0105086

Xiang, Z., Yu, Y., Yang, M., Yang, J., Niao, M., & Li, M. (2011). Does flagship species tourism benefit conservation? A case study of the golden snub-nosed monkey in Shennongjia National Nature Reserve. *Chinese Science Bulletin*, *56*, 2553–2558. https://doi.org/10.1007/s11434-011-4613-x

Yang, J., Liao, M., Yu, H., & Yao, H. (2008). Protection and research status of golden monkeys in Shennongjia. *World Sci-Tech R & D*, *30*, 418–421. (in Chinese).

Young, B. E., Hall, K. R., Byers, E. A., Gravuer, K., Hammerson, G., Redder, A., & Szabo, K. (2012). *Rapid assessment of plant and animal vulnerability to climate change* (pp. 129–152). Chicago, IL: University of Chicago Press.

Zhu, Z. (1992). Status and conservation of snub-nosed monkey at Shennongjia Nature Reserve. *Chinese Journal of Wildlife*, *3*, 16–17. (in Chinese).

## SUPPORTING INFORMATION

Additional supporting information may be found online in the Supporting Information section at the end of the article.

**How to cite this article:** Zhang Y, Clauzel C, Li J, et al. Identifying refugia and corridors under climate change conditions for the Sichuan snub-nosed monkey (*Rhinopithecus roxellana*) in Hubei Province, China. *Ecol Evol*. 2019;9:1680–1690. https://doi.org/10.1002/ece3.4815