Integrating Climate Change Adaptation and Mitigation into Land Use Optimization: A Case Study in Huailai County, China

: Land use plays a crucial role in climate change adaptation and mitigation, as the reasonable design of land use distribution can positively impact these things. Therefore, research interest in climate change adaptation and mitigation strategies in land use and management has been growing. However, the adaptation and mitigation strategies have been handled separately at different dimensions and spatial levels. In this study, we presented a modeling framework for land use optimization that integrates climate change adaptation and mitigation, developed the model, and then applied it to Huailai County, wherein environmental and socioeconomic conditions are sensitive to climate change. The regional land use optimization model was combined with a linear programming model and a modiﬁed cellular automata model. Subsequently, the climate change adaptation and mitigation constraints, including ecological water demand, spatial suitability, and carbon sequestration, were incorporated into the model. The results indicate that most regions in the study area could adapt to and mitigate climate change with a constant land use pattern, and the land use conversion region under different climate change scenarios was primarily located in the topography transition region. The optimization results also reveal trade-offs between climate change adaptation and mitigation that were manifested with an increase in carbon sequestration and ecological water demand accompanied by decreases in the net income of agricultural production. Thus, it is necessary to simultaneously incorporate climate change adaptation and mitigation into land use optimization and management, and the proposed model provides a feasible method to incorporate them and balance their trade-offs in land use pattern optimization at a regional scale. pattern optimization integrated climate change adaptation and mitigation.


Introduction
Climate change is widely considered one of the greatest challenges to modern human civilization that has profound socio-economic and environmental impacts [1][2][3]. The general solution to minimize the negative impacts of climate change on humans and ecosystems can be categorized as climate change adaptation and mitigation [4]. Adaptation involves reducing costs by changing conditions and exploiting opportunities, mitigation involves reducing greenhouse gas (GHG) emission and mitigating future climate change [5,6].
Reasonable land use/cover management can play a positive role in climate change adaptation and mitigation [7]. Climate change adaptation and mitigation strategies for land use constitute key scientific foci for climate change sciences and land use sciences [8,9]. For adaptation, effective land use/cover management strategies can enhance the adaptive capacity and reduce vulnerability to the negative effects due to climate change (e.g., drought stress and water scarcity) [10,11]. For mitigation, land use and land cover change (LUCC) can be treated as a major source of GHG emissions [12,13]. Therefore, enhancing the function of land-based carbon sinks would be an effective climate change mitigation strategy.
Integrated climate change mitigation and adaptation studies which handle the topics of regional land use/cover management are scarce, because the corresponding measures,

Modelling Framework
In this study, we proposed a regional land use optimization model that could represent the optimized land use structure and spatial pattern required to adapt and mitigate different climate change scenarios. This model incorporated two types of constraint factors, that is, the climate change adaptation and mitigation constraints ( Figure 1). We adopted the linear programming (LP) and modified cellular automata (CA) algorithm for the optimization of land use. The LP algorithm was used to optimize the quantitative structure of different land use types, and the modified CA algorithm was applied to rearrange the spatial locations (pattern) of land use patches in the study area. In the model, the constraints for climate change adaptation primarily involved the suitability of land use, available water resources and social adaptation; the mitigation constraints specifically included the carbon sequestration of regional land use.

Study Area
The study area (Huailai County) is located in the southeast of APTZ of Northern China (115 • 16 -115 • 58 E, 40 • 04 -40 • 35 N) and covers an area of 178,701 ha ( Figure 2). It is dominated by a temperate semi-arid continental monsoonal climate, with an average annual precipitation of 409.61 mm in growth season (April to October) and mean annual temperature of 19.43 • C. The study area includes cultivated land, garden land, woodland, shrub/grassland, water, and construction land. Huailai County has recently been experiencing drastic LUCC, which along with climate change, has seriously affected the local ecosystem services of this area, such as grain yield and reservoir capacity [18]. Therefore, conducting this study in this region is important.

Study Area
The study area (Huailai County) is located in the southeast of APTZ of Northern China (115°16′-115°58′ E, 40°04′-40°35′ N) and covers an area of 178,701 ha ( Figure 2). It is dominated by a temperate semi-arid continental monsoonal climate, with an average annual precipitation of 409.61 mm in growth season (April to October) and mean annual temperature of 19.43 °C. The study area includes cultivated land, garden land, woodland, shrub/grassland, water, and construction land. Huailai County has recently been experiencing drastic LUCC, which along with climate change, has seriously affected the local ecosystem services of this area, such as grain yield and reservoir capacity [18]. Therefore, conducting this study in this region is important.

Climate Change Scenarios
According to the emission scenarios proposed by IPCC, our study designed three different climate change scenarios (A2, A1B, and B1) for land use optimization ( Table 1). The 2050 climate change prediction results for different emission scenarios were obtained from the Climate Wizard Webtool (http://www.climatewizard.org/, accessed on 25 February 2021). Scenario A2 represented an extreme hot and drought climate condition, with an average temperature of 23.06 °C, and total precipitation of 373.57 mm in the growing season; scenario A1B represented moderate high temperature and drought climate condition, with an average temperature of 22.94 °C, and total precipitation of 390.0 mm in the growing season; scenario B1 represented a warm and humid climate condition, with an average temperature of 20.94 °C. and total precipitation of 458.33 mm in the growing season.

Climate Change Scenarios
According to the emission scenarios proposed by IPCC, our study designed three different climate change scenarios (A2, A1B, and B1) for land use optimization ( Table 1). The 2050 climate change prediction results for different emission scenarios were obtained from the Climate Wizard Webtool (http://www.climatewizard.org/, accessed on 25 February 2021). Scenario A2 represented an extreme hot and drought climate condition, with an average temperature of 23.06 • C, and total precipitation of 373.57 mm in the growing season; scenario A1B represented moderate high temperature and drought climate condition, with an average temperature of 22.94 • C, and total precipitation of 390.0 mm in the growing season; scenario B1 represented a warm and humid climate condition, with an average temperature of 20.94 • C. and total precipitation of 458.33 mm in the growing season.

Data Collection
The data used in the study includes the land use map and meteorological data of 1990 and 2010, vegetation index and net primary productivity (NPP) dataset, digital elevation model (DEM), soil types, and properties map.
The land use maps of Huailai County in 1990 and 2010 were obtained from the Landsat images (30 m spatial resolution). The MODIS-NDVI and NPP datasets in 1990 and 2010 were obtained to identify ecological function hot areas. Additional details about remote sensing image data are listed in Table 2. The meteorological datasets in 1990 and 2010 were obtained for seven weather stations located in and around Huailai County, from the China Meteorological Data Sharing Service System (http://www.geodata.cn/, accessed on 20 May 2021). The datasets included the observation data of temperature, precipitation, and solar radiation. Then, the kriging interpolation method was used to generate the spatial data of these meteorological data.
The DEM data (30 m spatial resolution) were obtained from the Geospatial Data Website (http://www.gscloud.cn/, accessed on 20 May 2021). Spatial soil data were derived from the second General Soil Survey of Huailai County, and it comprised soil information, such as soil type, texture, nitrogen content, and organic matter content. All the spatial data were resampled to 30 m × 30 m resolution.
In addition, socio-economic statistical data, such as population, agricultural product yield and price, were primarily obtained from Huailai Statistics Yearbook of Social and Economics.  [13], and the IPCC pointed out that enhancing the land carbon sinks through adjusting land use patterns and management seems to be an effective plan [28]. Therefore, we used the regional carbon storage of land use as a proxy for carbon sequestration and climate change mitigation. It is the product of carbon density and area of different land use types. The carbon sequestration due to land use is contributed by vegetation carbon density and soil carbon density. The vegetation carbon density was estimated from the net primary productivity (NPP) by the Carnegie Ames Stanford Approach (CASA) model [29]. The soil carbon density of different land use types was surveyed and determined based on the soil organic carbon (SOC) content of surface soil  (40 cm) [30,31]. The calculated vegetation and soil carbon density of different land use in the study area are listed in Table 3. The regional carbon sequestration constraint was calculated using the following equation: where max 1 is the maximum value of regional total carbon storage of land use, X 1 , X 2 , X 3 and X 4 are the area (ha) of cultivated land, garden land, woodland, and shrub/grassland, respectively; 19.84, 27.53, 51.29, and 42.41 denote the carbon density (sum of vegetation and soil carbon density) of each type of land use (t/ha).

Social Adaptation
The important aspect of social adaptation to climate change involves improving farmers' incomes and livelihoods to reduce social vulnerability [9]. Thus, we established a net income target function as the regional social adaptation function: max 2 = 11,403.47X 1 + 29,177.68X 2 + 906.74X 3 + 4201.00X 4 (2) where max 2 is maximum value of the net income of regional cultivated land, garden land, woodland, and shrub/grassland; and X 1 , X 2 , X 3 , and X 4 are the same as in Formula (1). The constant represents the per unit area economic benefit of different land use types which were obtained from social and economic statistical data.

Available Water Resources
To minimize the negative impacts of limited water resources, and to adapt to drought, we calculated the ecological water demand of cultivated land, garden land, woodland, and shrub/grassland using the ecological modeling approach. The available water resources constraints are as follows: where X 1 , X 2 , X 3 , and X 4 are the same as those in Formula (1); the constant represents the water demand per unit area of four land use types (mm/ha). The quantitative results correspond to mean water demand statistical results of different land use types in 1990 and 2010. According to the results of Zhao [32] Huang, Yu [33], and Allen, Walter [34] the available ecological water of the study area under A2, A1B, and B1 scenarios were set as 4.64 × 10 8 m 3 , 4.84 × 10 8 m 3 , and 5.53 × 10 8 m 3 , respectively.

Land Use Suitability
In this study, the land use suitability evaluation was conducted to identify the areas that were suitable for specific land use in future climate change scenarios. The suitability of cultivated land, garden land, woodland, and shrub/grassland was evaluated based on climate, topography, soil, and social factors (Table 4). Construction land was evaluated based on topography and distance factors (Table 5).
To determine the climatic factor suitability score of agricultural land, the study established the grain and fruit yield estimation equation to assign cultivated land and garden land score and examined the relationship between vegetation NPP and climate index (ombrothermic) for woodland and grassland. Based on the collected meteorological and yield statistical data, the grain and fruit yield were fitted via the least square method, the estimation equation was y grain = 1.316 + e 0.011x and y fruit = 0.002x 2 − 0.478x + 28.351, where x denotes the ombrothermic index from June to September (calculation equation shown in Table 4), and R 2 values of the equations correspond to 0.49 and 0.62, respectively. We then predicted the grain and fruit yield of different climate change scenarios in 2050 based on the yield established suitability score, and the results are presented in Table 4. Based on the corresponding relationship between NPP and ombrothermic index, we determined the climate factor suitability score of woodland and grassland. The zones with an ombrothermic index in the range of 220-240 exhibited the highest NPP and corresponded to the most suitable zone for forests and grass growth, and the suitability score was 3. The zones with an ombrothermic index less than 210 exhibited the lowest NPP and were marginally suitable for woodland and grassland, and the suitability score was 1. The other zone was generally suitable for woodland and grassland, and the suitability score was 2.
The suitability scores of all factors were determined by consulting with a team of experts. In the evaluation process, each evaluation factor was divided into different classes, and a scale from 0 to 3 (0 = not suitable, 1 = marginally suitable, 2 = moderately suitable, and 3 = highly suitable) was used to assess the suitability of different classes (Tables 4 and 5). The weights of four agricultural land evaluation factors were determined via the analytic hierarchy process (AHP) [35,36]. The weight results of all factors are listed in Table 6.    The land use suitability evaluation results of cultivated land, garden land, woodland, and shrub/grassland are shown in Figure 3. The suitability results suggested that from scenarios A2 to A1B and B1 with the climate becoming wetter, the cultivated land gradually occupied the plain and expanded to the hills, the garden land extended to the low altitude plain, and the woodland and shrub/grassland suitability region expanded to the hills.   The land use structure optimization formula comprised the target function and constraint conditions. The target function was formed by maximum carbon sequestration (Formula (1)) and social adaptation (Formula (2)). We normalized the coefficients of each formula and established diverse weights for the two formulas in different scenarios. In A2 and A1B scenarios (extreme and moderate high temperature and drought scenarios), the weight of the normalized carbon sequestration formula was 0.6, and the social adaptation formula was 0.4. In scenario B1 (the warm and humid scenario), these normalized formulas were given equal weight.
The constraint conditions primarily depended on three aspects: available water resources (Formula (3)), land use constraints which are consistent with CO 2 emission background (Table 1), and land use optimization suggestions. The study established constraint conditions for six variables as follows: The first constraint condition was established as follows: where X 1 , X 2 , X 3 , and X 4 are the same as in Formula (1), and X 5 denotes area of water, and X 6 denotes the area of construction land. Subsequently, the following constraint conditions were established based on the development background of different emission and climate change scenarios (Table 1).
(1) We established the constraint condition of cultivated land based on the population and food demand. In the A2 and A1B scenarios (high population growth), population was predicted as 427,574 in 2050; in the B1 scenario (slower growth rate), the population was predicted as 397,400 in 2050.
(2) Based on the development description of different scenarios (Table 1), we proposed that the increase rate of construction land for A2, A1B, and B1 scenarios were 2%, 1.5%, and 1.0%, respectively. The constraint conditions are then established as follows: A2 : X 6 = 20,272 A1B : X 6 = 16,654 B1 : X 6 = 13,669 (3) The environmental protection policies were different between the three scenarios; B1 scenarios increasingly focused on environmental protection although A2 and A1B scenarios ignored these. Based on the above environmental protection background, we established the woodland, shrub/grassland, and water constraint conditions.
We extracted woodland (29,764.05 ha) and shrub/grassland (12,509.7 ha) preserved under A2 and A1B scenarios that are located in the high ecosystem function sub-basins; woodland (37,718.86 ha) and shrub/grassland (29,903.71 ha) preserved under B1 scenarios were located in the high and medium ecosystem function sub-basins. Hence, the constraint conditions can be established as follows: A2 : X 3 ≥ 29,764.05 X 4 ≥ 12,509.7 A1B : X 3 ≥ 29,764.05 X 4 ≥ 12,509.7 B1 : X 3 ≥ 37,718.86 X 4 ≥ 29,903.71 Subsequently, we proposed the water protection constraint conditions as follows: A2 and A1B scenarios preserved the water area in 2010 (8539 ha), and B1 scenario recovered to the range observed in 1990 (13,346 ha). The constraint conditions can be expressed as follows: A2 : X 5 = 8539 A1B : X 5 = 8539 B1 : X 5 = 13,346 Finally, the following two constraint conditions were constructed considering the land use optimization and management suggestions of study area.
(1) Based on numerous agroecological parameters, Zhang and Tang developed landscape pattern optimization research in the study area and proposed that cultivated land, garden land, and shrub/grassland should be located at the valley plain and low hilly areas [41]. Statistical results indicated that the total area of the region is 132,843 ha. Considering that construction land and water area were also primarily distributed in this region, we deducted these areas and obtained the area of cultivated land, garden land, and shrub/grassland located at the valley plain and low hilly areas. The constraint conditions can be expressed as follows: A2 : X 1 + X 2 + X 4 ≤ 104,032 A1B : (2) Zhong proposed that the garden land of the study area should not exceed 43,335 ha based on the water consumption analysis [42], and thus, the constraint condition can be established as follows: The optimized land use structure was determined using LINGO 10.0 software.

Land Use Pattern Optimization
To integrate with adaptation and mitigation strategies in climate change, the general land use pattern optimization procedure was determined and shown in Figure 4. The land use pattern optimization processes considered three basic criteria: (i) the constraint of land use area; (ii) land use suitability of different climate scenarios; (iii) judgment of neighborhood conditions. The land use area constraint was obtained from the land use structure optimization result (Section 2.6.1), and the land use suitability was from the land use suitability evaluation (Section 2.5.4). These data were used as input data to implement land use area judgment and the highest land use suitability cell selection in optimization processes. The land use optimization procedures were executed in terms of the land use in the order of woodland, construction land, cultivated land, garden land and shrub/grassland. The optimization procedure of the next one begins when the area of former land use exceeds the land use structure optimization result. The awaited optimized cell was selected based on the land use suitability result of different scenarios, with highest land use suitability.
In this study, the land use pattern optimization method is based on a modified CA algorithm, and it is also semi-automatic, which allows users to manually identify some optimization areas [43]. The CA algorithm is defined as a discrete spatial-temporal dynamic systems based on local rules and often includes four elements [44]. The elements of the modified CA algorithm are as follows: (i) the extent or grid data of land use/cover, (ii) a 5 × 5 kernel for neighborhood analysis, (iii) the manually identified optimization areas include construction land in 2010 and the spatial data of environmental protection areas (woodland, shrub/grassland and water environmental protection area in Section 2.6.1) as the initial states and input data, and (iv) neighborhood condition judgment rules. The land use allocation should consider the neighborhood conditions which represent the interrelationship among surrounding land use types. We used compatibility matrix and cell statistics to express the neighborhood conditions.
(1) Compatibility matrix. The value of spatial compatibility corresponds to 0 or 1 where 1 denotes the compatibility between the neighboring two land use types (generally adjacent to each other), and 0 denotes no compatibility (generally isolated from each other). (2) Cell number statistics. The rule of surrounding land use requires that the number of the cell belonging to optimized land use should exceed that of the other land use types.
When the cell met the requirements of neighborhood conditions, the cell value 0 was changed into explicit land use types. The land use allocation should consider the neighborhood conditions which represent the interrelationship among surrounding land use types. We used compatibility matrix and cell statistics to express the neighborhood conditions.
(1) Compatibility matrix. The value of spatial compatibility corresponds to 0 or 1 where 1 denotes the compatibility between the neighboring two land use types (generally adjacent to each other), and 0 denotes no compatibility (generally isolated from each other). (2) Cell number statistics. The rule of surrounding land use requires that the number of the cell belonging to optimized land use should exceed that of the other land use types.
When the cell met the requirements of neighborhood conditions, the cell value 0 was changed into explicit land use types.     Then, we compared the land use structure optimization results with those of 2010 and observed general land use optimization as follows: the area of shrub/grassland decreased while the cultivated land, construction land, garden land and woodland increased. Woodland is an important carbon sink land use type, and afforestation is beneficial to mitigate future climate change. The garden land exhibited higher profit, and the increased cultivation area aided in increasing the social adaptation capacity of the study area. In addition, under the warm and humid scenario (B1), shrub/grassland decreased the least, the extreme, drier and hotter scenarios (A2) exhibited a medium decrease, and the medium drier and hotter scenario (A1B) exhibited a high decrease. Table 7 also presents the optimization efficiency of the proposed model in future scenarios. The results indicate that the adaptation and mitigation efficiency of optimized land use significantly changed. The net income of land use optimization results increased by approximately 0.15 × 10 5 million compared with 2010. The carbon sequestration and ecological water demand of the optimization result decreased when compared with 2010. Furthermore, we compared the optimization results of future scenarios and observed only a low difference in adaptation and mitigation efficiency between the different climate change scenarios. This also demonstrated that rational adjustment of land use structure and its spatial pattern effectively decrease the effects of future climate change and maximize the climate change adaptation and mitigation efficiency of the land system.

Land Use Pattern Optimization Results
The land use pattern optimization maps under different climate change scenarios are shown in Figure 5. Evidently, future climate change scenarios exhibited the general land use pattern optimization contributions as follows: cultivated land, construction land and most of garden land located in the plain of the study area, a small part of the garden plots and most of the shrub/grassland located in the hills, and the woodland primarily located in the mountains of the study area. This was because the factors that affect the land use pattern included climatic factors, topography, and soil factors. Specifically, in the mountain-basin system (similar to the study area), the topography factors (altitude and slope) become increasingly important to identify land use suitability and spatial pattern. area, that between scenarios A1B and B1 corresponded to 74.2%, and that between scenarios A2 and B1 corresponded to 73.7%. The overlay results indicated that the changed land use regions among the different scenarios were primarily located in the transition region of the three topography regions, and the major land use types corresponded to shrub/grassland, garden land, and woodland. The results also imply that most regions in the study area could adapt to climate change with a constant land use pattern.

Incorporating Climate Change Adaptation and Mitigation into Land Use Planning
The results indicate that land use optimization improved the climate change adaptation and mitigation efficiency at a regional scale. Furthermore, in the study area, most regions with a constant land use pattern adapted to and mitigated climate change through rational land use planning. Therefore, decision makers should consider the role of land use planning when formulating climate change adaptation and mitigation strategies.
Additionally, the results indicate that the increase in ecological land use water demands was followed by an increase in climate carbon sequestration. Being the important land use type for carbon sequestration, the water demand for woodland was higher than the others. In the growing season, the water demands for shrub/grassland, garden plots, and cultivated land per unit area (1 ha) corresponded to 0.28 mm, 0.26 mm, and 0.24 mm, Furthermore, we compared the optimization results between each two scenarios and observed that the land use of most regions in the study area was unchanged ( Figure 6). The unaltered region between scenarios A2 and A1B corresponded to 84.1% of that study area, that between scenarios A1B and B1 corresponded to 74.2%, and that between scenarios A2 and B1 corresponded to 73.7%. The overlay results indicated that the changed land use regions among the different scenarios were primarily located in the transition region of the three topography regions, and the major land use types corresponded to shrub/grassland, garden land, and woodland. The results also imply that most regions in the study area could adapt to climate change with a constant land use pattern.
Land 2021, 10, x FOR PEER REVIEW 14 of 20 area, that between scenarios A1B and B1 corresponded to 74.2%, and that between scenarios A2 and B1 corresponded to 73.7%. The overlay results indicated that the changed land use regions among the different scenarios were primarily located in the transition region of the three topography regions, and the major land use types corresponded to shrub/grassland, garden land, and woodland. The results also imply that most regions in the study area could adapt to climate change with a constant land use pattern.

Incorporating Climate Change Adaptation and Mitigation into Land Use Planning
The results indicate that land use optimization improved the climate change adaptation and mitigation efficiency at a regional scale. Furthermore, in the study area, most regions with a constant land use pattern adapted to and mitigated climate change through rational land use planning. Therefore, decision makers should consider the role of land use planning when formulating climate change adaptation and mitigation strategies.
Additionally, the results indicate that the increase in ecological land use water demands was followed by an increase in climate carbon sequestration. Being the important land use type for carbon sequestration, the water demand for woodland was higher than the others. In the growing season, the water demands for shrub/grassland, garden plots, and cultivated land per unit area (1 ha) corresponded to 0.28 mm, 0.26 mm, and 0.24 mm,

Incorporating Climate Change Adaptation and Mitigation into Land Use Planning
The results indicate that land use optimization improved the climate change adaptation and mitigation efficiency at a regional scale. Furthermore, in the study area, most regions with a constant land use pattern adapted to and mitigated climate change through rational land use planning. Therefore, decision makers should consider the role of land use planning when formulating climate change adaptation and mitigation strategies.
Additionally, the results indicate that the increase in ecological land use water demands was followed by an increase in climate carbon sequestration. Being the important land use type for carbon sequestration, the water demand for woodland was higher than the others. In the growing season, the water demands for shrub/grassland, garden plots, and cultivated land per unit area (1 ha) corresponded to 0.28 mm, 0.26 mm, and 0.24 mm, respectively, whereas that of the woodland corresponded to 0.42 mm (this was approximately 1.5 times that of others). The efficiency of land use optimization results indicates that trade-offs and synergies existed between climate change adaptation and mitigation [45,46]. For example, woodland created more residual vegetation and exhibited increased carbon sequestration while demanding relatively higher ecological water, thereby leading to lower economic profits for farmers. Unlimited expansion of the carbon sequestration vegetation for climate change mitigation in water-limited and rural areas could induce potential conflict with demands for water and the improvement of human livelihoods [6]. Therefore, trade-offs exist between climate change mitigation and adaptation, and similar results were also found in other previous studies [11,47,48]. Wu, Wang suggested that water provision and water sustainability are among climate change adaptation policies and correspond to crucial issues of future land use management under the background of climate change, that is, a hot and dry climate [47]. In addition, the rural human livelihoods are dependent on environmental resources, and the success of mitigation measures typically depends on how well the community adapt to climate changes (e.g., drought and erratic rainfall) [26]. Thus, the efforts that can adapt to the climate changes in the areas and integrate natural resource management can only enhance or address livelihood issues and thus increase resilience in response to climate change in a sustainable manner [49].
The complex relationships also indicate that land use managers focus on one aspect of climate change adaptation or mitigation, which would be unsustainable land use and regional development states. Thus, strategies to minimize the trade-offs between them are equally necessary, and climate change adaptation and mitigation should be simultaneously considered to improve land use planning and decision making. The modeling methodology emphasized an integrated approach to optimally achieve the simultaneous multiple benefits of climate change adaptation and mitigation. Furthermore, by minimizing trade-offs between adaptation and mitigation, maximum efficiency results can be achieved. Moreover, a comprehensive analysis of climate change adaptation and mitigation in the land use pattern or land use management research can effectively resolve trade-offs between them and can promote sustainable development.
These are some analysis methods, such as the Pareto-optimal solutions, which could help decision makers to reveal the complicated relationship between the climate change adaptation and mitigation strategies, such as the [50,51], the objective function can be solved over a range of relative weight.

Modeling Study of Land Use Optimization Integrated Climate Change Adaptation and Mitigation in the Future
The proposed methodology optimized the land use pattern by integrating climate change adaptation and mitigation strategies on a multiscale. We developed a regional scale model using land use structure optimization results that integrated regional climate change optimization and mitigation constraints as the pattern optimization condition. The neighborhood condition and land use suitability map under different scenarios were used at a pixel scale. In addition, the land use optimization model proposed a simple model structure and exhibited user-defined organizational characteristics, such that researchers or decision makers could add any other adaptation and mitigation constraints and strategies.
There are a few limitations in this study that can be considered in further model research. First, the climate change adaptation and mitigation strategies or constraints incorporated into the model are limited, and a broader scope of adaptation and mitigation should be investigated. For example, the model considered limited adaptation strategies, and more constraints should be integrated into the model such as improved water use, adapted crop varieties and practices, soil and water conservation practices, and public awareness of climate change and variability [49]. Furthermore, the model did not consider the GHG emissions of land use while calculating climate change mitigation. Studies have emphasized that the shrub/grassland used for livestock is an important contributor to GHG emissions despite constituting the livelihoods of certain communities [14]. Second, the land use types modelled in the study area were approximate due to the limitation of source data and model complexity. To improve this, the vegetation or crop distribution data and the distribution model that can value the spatial suitability of species can be applied to concretize model results. The general species model is the maximum entropy (Maxent) model that is widely used for modeling species distributions and can obtain a suitability map of specific crops [11]. Finally, climate feedback from adaptation and mitigation efforts was not integrated into the model. The system dynamics (SD) model can be used to expand and model the feedback of a land use system, and it is used to solve economic, social, environment, and other system problems, such as nonlinear, multivariable, higher order, and multiple feedback [52]. The advantages of the SD model can be applied without sufficient data to analyze the features of the system, and effectively track and model the real system.
Finally, the spatial model established in this study ignored some key features of land use change, such as previous land use patterns and the spatiotemporal autocorrelation in land use patterns [53]. To further optimize the model proposed, we should also focus on the complex relationship between the socioeconomic system and land use system [54].

Climate Change Adaptation and Mitigation Strategies in Land Use Management
The modelling results indicate that the land use pattern change regions among different climate change scenarios are primarily located in the topography transition regions between plain, hills, and mountain. There were studies indicating that transitional zones among biomes or ecosystems are particularly sensitive to changing climatic conditions [10]. Specifically, some studies on montane regions exhibit a general movement of species upward with elevation. Luckman and Kavanagh [55] observed the upward movement of tree lines in Siberia and in the Canadian Rocky Mountains where the temperature had increased by 1.5 • C. A study across 26 mountains in Switzerland documented that alpine flora had expanded toward the summits since the 1940s [56]. All the studies indicated that an elevation shift in biomes or ecosystems occurs in the mountain area. Furthermore, considering the terrain constraints and climatic conditions, the spatial distribution of species changes with increasing elevation, such that the boundaries of biomes or ecosystems are always located in the topography transition region. The study area in this study was a typical mountain-basin system exhibiting significant differences in topography conditions, and the topography transition region of the study area corresponded to the key area for developing land use optimization and management to adapt to future climate change.
The climate change analysis of the study area indicated that the annual mean temperature has gradually increased at a rate of 0.37 • C/10a over the past 60 years, which exceeds the national average. The annual precipitation gradually decreased, and the warm and dry tendency of summer was more evident than other seasons [25]. Furthermore, to reveal the effect of dryness on land use optimization results, we proposed the continuous drought hypothesis (the available water resources are drier than the A2 scenario). The land use structure optimization result indicated that the area percentage of cultivated land increased to 44%, garden plots decreased to 8%, woodland decreased to 26%, and shrub/grassland decreased to 7%. The results also illustrated that the agricultural structure of the study area tended to be simple to adapt to a dry climate and continued to decrease the regional available water resources. Similar results were also observed in previous studies [47]. Many species suffered a reduced habitable area due to recent climate changes. With respect to species that are already prone to extinction at their equatorial or lower range boundaries, some have either failed to expand poleward or could not to expand due to geographic barriers [10]. These species suffered absolute reductions in range size, thereby increasing the risk of extinction in the future. In the Great Basin of the Western United States, 7 out of 25 re-censused populations of the pika were extinct since they were recorded in the 1930s [57]. Observations and experiments indicated that the adult pika stopped foraging in the midday heat in August and were killed within a half hour at more than 31 • C [58].
Additionally, climate change adaptation and mitigation strategies also include land use management policies. The decision makers can obtain useful references from successful land use management strategies of extant studies on drier and hotter climates [49]. Climate change adaptation strategies include the following: (1) improved water use through sinking shallow wells, water harvesting, increased efficient water use and surface cover; (2) use of adapted crop varieties and practices, such as drought-resistant crops, planting locally adapted species; (3) adopting soil and water conservation practices including soil fertility improvement, crop rotation, farming strategies (mulching, early maturing varieties), and organic manure to enhance soil water-holding capacity. Climate change mitigation strategies also involve reducing atmospheric CO 2 concentrations through methods such as tillage practices. For example, a combination of inorganic fertilizer and organic manure and crop rotation systems are effective in increasing SOC [8]. Finally, public awareness/information dissemination on climate change has also been proven effective in adapting and mitigating climate change.

Conclusions
In this study, we established a methodology and provided an empirical land use optimization method that incorporates climate change adaptation and mitigation. We obtained the following findings: (1) The results indicated that rational adjustments of land use structure and pattern can be an effective strategy to adapt and mitigate climate change at a regional scale in the face of the changing climate. In addition, the present findings revealed that trade-offs exist between climate change adaptation and mitigation. Ultimately, land use planning and policy should jointly consider adaptation and mitigation to cope with climate change and to achieve sustainability. Furthermore, it is necessary to develop methodology to minimize trade-offs between adaptation and mitigation and then achieve maximum efficiency results. (2) The study area corresponded to a typical mountain-basin system wherein the optimized land use pattern characteristics including the following: cultivated land, construction land and most of garden land located in the plain of the study area, a small part of garden land and most of shrub/grassland located in the hills, and the woodland primarily located in the mountainous area. In addition, the changed land use region with different climate change scenarios was primarily located in the transition region of three topography regions (plain, hills, and mountain). In the region, the major changed land use types included shrub/grassland, garden land, and woodland. Considering the topographical characteristics of the study area, it is important to extend the methodology to other regions.
Furthermore, the model limitations should be considered in further research. Future modeling efforts can incorporate a broader scope of adaptation and mitigation into the model framework. In addition, future research can integrate the distribution model and SD model to potentially improve the scientific value of the optimization results.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the privacy agreement among co-authors.