Integrating Modelling and Expert Knowledge for Evaluating Current and Future Scenario of Large Cardamom Crop in Eastern Nepal

Large Cardamom (Amomum subulatum Roxb.) is one of the most valuable cash crop of the Himalayan mountain region including Nepal, India, and Bhutan. Nepal is the world’s largest producer of the crop while the Taplejung district contributes a 30%–40% share in Nepal’s total production. Large cardamom is an herbaceous perennial crop usually grown under the shade of the Uttis tree in very specialized bioclimatic conditions. In recent years, a decline in cardamom production has been observed which is being attributed to climate-related indicators. To understand the current dynamics of this under-canopy herbaceous crop distribution and its future potential under climate change, a combination of modelling, remote sensing, and expert knowledge is applied for the assessment. The results suggest that currently, Uttis tree cover is 10,735 ha in the district, while 50% (5198 ha) of this cover has a large cardamom crop underneath. When existing cultivation is compared with modelled suitable areas, it is observed that the cultivatable area has not yet reached its full potential. In a future climate scenario, the current habitat will be negatively affected, where mid elevations will remain stable while lower and higher elevation will become infeasible for the crop. Future changes are closely related to temperature and precipitation which are steadily changing in Nepal over time.


Introduction
Large cardamom (Amomum subulatum Roxb.) is mostly being cultivated in the Himalayan mountain region of Nepal, India, and Bhutan [1]. There is an increasing demand of the spices from local to global markets which has fascinated farmers in its commercial cultivation [2]. The high-value crop cultivation has substantially improved the livelihood of farmers residing in these mountain regions. However, several studies have highlighted threats of cardamom farming which include problems in disease management, change in climatic conditions, and human activities such as infrastructure development in cardamom growing areas [2,3]. Alongside this, comprehensive information on its current distribution and potentially suitable areas for cultivation is unavailable, which is essential for crop planning and management. Further, for its distribution and current suitability, information on the impacts of climate change is also vital for future planning as climate has a significant impact on the geographical distribution of plant species and also alters habitat conditions [4,5].
Large cardamom is an herbaceous perennial crop usually grown under shade. Uttis (Alnus nepalensis) trees provide excellent shade, supply a good amount of litter from twigs and leaves, and nitrogen from the root nodules to understory cardamom when they are young. The crop is generally grown at an altitude of 700 to 2000 m above sea level in humid conditions, with temperature ranges between 4-20 • C and annual precipitation around 2000-2500 mm [6].
Although satellite remote sensing data can be used for accurate, timely, and consistent information on the agricultural productivity at local and regional scales [7,8], detection of understory plants is inhibited due to canopy cover, canopy gap shadowing, and terrain variability [9]. Differentiating signals of understory vegetation from overstory canopy is still challenging due to complex interactions between overstory and understory vegetation [10]. Several studies have been conducted to distinguish understory vegetation using remote sensing approaches. In this regard, evergreen understory vegetation were identified from Landsat images of leaf off-season in deciduous forest [11,12]. Phenological difference between overstory and understory vegetation was also used to detect understory vegetation using Landsat images [13,14]. However, such studies require multi-temporal data which are generally available in coarser spatial resolutions. Leduc et al. [15] have mapped wild leek which grows on forest floors using low flying Unmanned Aerial Vehicle (UAV). Other studies have used data from active sensor, mainly LIDAR data, to map understory plants in boreal forests [16,17]. Combinations of LIDAR data and hyperspectral images were used to map understory invasive species in tropical forests [9], and the use of LIDAR data and high-resolution IKONOS imagery to identify understory plant invasion in urban forests [18]. Nevertheless, there are limitations on the use of these data due to its high cost and low availability, mainly in developing countries [19].
Understory vegetation are usually hard to identify only from remote sensing, thus a substantial efforts have been made to identify these vegetation through the integration of several approaches. For instance, Wang et al. [20] have mapped understory bamboo by integrating neural network and Geographic Information System (GIS) expert system, and Tuanmu et al. [10] have detected understory vegetation using phenology metrics derived from time series of Moderate Resolution Imaging Spectroradiometer (MODIS) data together with a suitability model like the maximum entropy (Maxent) model. The Maxent model is a species distribution model (SDM) that provides an approach for making predictions from existing distribution and a set of predictors [21]. Such models predicting the potential distribution of species are useful for several applications in conservation biology [22,23]. Other approaches, such as object-based image analysis (OBIA) on multispectral and hyperspectral data, were found to be effective in identifying understory plant species [24] due to its characteristics of using contextual relationships together with spectral information. In addition, expert knowledge [25] and other ancillary data such as elevation, slope, and aspect can also be used to improve classification accuracy. Moreover, the combination of participatory mapping with remote sensing technique can further improve the accuracy as the two methods complement and validate each other [26]. Participatory mapping is an effective tool to obtain accurate baseline data of the field [27], which can be integrated in remote sensing analyses to produce more accurate maps [28].
Several species distribution models such as the genetic algorithm for rule-set production (GARP), ecological niche factor analysis (ENFA), bioclimatic modeling (BIOCLIM), CLIMEX (climate change experiment), domain environmental envelope (DOMAIN), and Maxent (maximum entropy) have been used for species distribution prediction. Among these, Maxent is widely used as it can perform better even with small sample sizes compared to other modelling methods [29][30][31]. In addition, it has other merits, such as: it requires only presence data of the species and environmental variables for the study area, it can use both categorical and continuous data and can incorporate interactions among different variables, it can produce spatially explicit habitat suitability maps, and the importance of each environmental variable on the model can be evaluated using the built-in jackknife test.
This study aims to evaluate existing spatial distribution, potential suitable areas for cultivation under various scenarios, and core distributional shift with future climate condition in the Taplejung district of Nepal. The study uses a combination of expert knowledge and species modelling approach to capture all aspects of habitat conditions in the complex environmental conditions in the mountain region.

Study Area
The study area, the Taplejung district, lies in eastern Nepal. The district ( Forest is the major land cover of the district followed by agriculture, bare land, grass, snow, shrub, water body, and built-up, respectively. The district is the major producer of large cardamom in the country, which is the major source of income for farmers in the district. In Nepal, an estimated 12,000 ha of land in over 40 mid-hills districts are under cardamom cultivation with estimated annual production of 6000 metric tons. While the Taplejung district contribute around 4500 ha of land and yearly production of around 2600 metric tons [32]. The study broadly followed two broad approaches for Cardamom habitat mapping including species modelling-and expert knowledge-based assessment ( Figure 1).

Study Area
The study area, the Taplejung district, lies in eastern Nepal. The district (27°57′10′′-27°16′5′′ N and 87°26′40′′-88°12′6′′ E) with an area of 3646 km 2 physiographically lies in the mid-hill to high-himal region. It lies at an elevation of 498 m to 8464 m. Forest is the major land cover of the district followed by agriculture, bare land, grass, snow, shrub, water body, and built-up, respectively. The district is the major producer of large cardamom in the country, which is the major source of income for farmers in the district. In Nepal, an estimated 12,000 ha of land in over 40 mid-hills districts are under cardamom cultivation with estimated annual production of 6000 metric tons. While the Taplejung district contribute around 4500 ha of land and yearly production of around 2600 metric tons [32].
The study broadly followed two broad approaches for Cardamom habitat mapping including species modelling-and expert knowledge-based assessment ( Figure 1).

Participatory Mapping
Participatory mapping with the community in developing countries is found to be effective for understanding local level phenomena [33]. A two-day exercise was conducted with a community representative to map existing large cardamom farming areas on the high-resolution data of Google Earth images available in digital format as well as on the printed maps. During the exercise, large clusters of cardamom cultivated fields were identified. However, detailed delineation of cardamom fields including small patches was not possible. The data was gathered to get existing farming areas on the district which was also cross-checked with outputs obtained from remote sensing data. A part of the pockets identified by the farmers in the participatory exercise was further verified through ground validation exercise. This validation also provided GPS samples of cardamom field points along with other vegetation in the area.

Uttis (Alnus nepalensis) Mapping Using High-Resolution Satellite Data
During the participatory mapping and the field activity, it was observed that the Uttis is the major tree species grown for shade to cultivate large cardamom. Since the understory crop could not

Participatory Mapping
Participatory mapping with the community in developing countries is found to be effective for understanding local level phenomena [33]. A two-day exercise was conducted with a community representative to map existing large cardamom farming areas on the high-resolution data of Google Earth images available in digital format as well as on the printed maps. During the exercise, large clusters of cardamom cultivated fields were identified. However, detailed delineation of cardamom fields including small patches was not possible. The data was gathered to get existing farming areas on the district which was also cross-checked with outputs obtained from remote sensing data. A part of the pockets identified by the farmers in the participatory exercise was further verified through ground validation exercise. This validation also provided GPS samples of cardamom field points along with other vegetation in the area.

Uttis (Alnus nepalensis) Mapping Using High-Resolution Satellite Data
During the participatory mapping and the field activity, it was observed that the Uttis is the major tree species grown for shade to cultivate large cardamom. Since the understory crop could not be directly mapped using remote sensing data, delineation of Uttis cover was consider as proxy to map the large cardamom crop. Our focus was to map Uttis in the entire district, which can provide information on prospective cultivation of large cardamom. Uttis tree cover was mapped by using Sentinel-2A Level 1C product image acquired in February 2016. The image has very high-resolution spectral coverage that includes 12 bands (coastal aerosol, blue green, red, 3 vegetation red edge, NIR, vegetation red edge water vapor, SWIR Cirrus, and 2 SWIR respectively). The spatial resolution of blue green, red and NIR is 10 m. The resolution of the vegetation red edge band is 20 m and for the rest of the bands it is 60 m. These features of the sensor are suited for agricultural monitoring systems [34]. Level 1C product is a Top of Atmosphere (TOA) product for which atmospheric correction had to be done to get reflectance values of the image so that the image can be used for mapping. SNAP (Sentinel toolbox) software (version 5.0.0) was used for atmospheric correction of the image.
Spectral separability of forest types such as coniferous forest, broadleaf forest, Uttis, and shrub was studied before the classification. The mean pixel values of the abovementioned forest types and shrub were plotted against eight bands of Sentinel-2A image to evaluate the potential of image spectral separability before the classification. The object-based image analysis (OBIA) classification approach was adopted for Sentinel-2A image classification. The technique uses spectral and contextual information in an integrative way [25]. The fundamental technique of OBIA is the segmentation of satellite images which overcome the salt and pepper effect [35]. In this study, the chessboard and multi-resolution segmentation algorithm in eCognition software (version 8.7,) was used to develop image objects [36] using scale parameter = 80, shape = 0.1, and compactness = 0.8.
After segmentation, Assign Class algorithm was used to classify general classes (agriculture, conifer forest, broadleaf forest, shrub, water, and snow). This was done to filter the land cover which were not the focus of the study. Several features such as NDVI, brightness, slope, elevation, and field information were used in this step. For the rest of the unclassified image objects, Nearest Neighbor Classification was applied, which is a powerful approach [37] to map Uttis trees. Uttis trees grown between the elevation range of 800-2200 m and slope up to 45 degrees were mapped. This is also the elevation and the slope range where large cardamom is cultivated [38]. Seventy percent (70%) of field data were used to train the samples while the remaining 30% were used as test data for accuracy assessment. Further, in order to reduce error and improve classification accuracy, interactive visual analysis was done on a classified map using Google Earth images [39].
The remote sensing-derived Uttis cover and participatory mapping-based identification of cardamom crop areas were overlaid to get the existing large cardamom farming area. The obtained large cardamom maps in the Taplejung district were further analyzed based on elevation, slope, aspect, and Village Development Committees (VDCs). The analysis was done to understand the pattern of cultivation and environment suitability conditions in the district and the information could be useful in further planning and management of large cardamom farming in the district. Elevation range of the study area was categorized into 9 ranges (below 800 m, 800-1000 m, 1000-1200 m, 1200-1400 m, 1400-1600 m, 1600-1800 m, 1800-2000 m, 2000-2200 m, and above 2200 m). In order to comprehend the farming area slope-wise, the slope was categorized into four gradient levels (0

Environmental Variables and Species Occurrence Records
The habitat suitability model was primarily developed based on the variables related to climate, soil, terrain, and vegetation type. Initially, 24 variables were selected to model the current distribution pattern. These comprised of 19 bioclimatic variables with 30 arc seconds (~1 km) spatial resolution from WorldClim dataset (http://www.worldclim.org/), 12.5 m resolution digital elevation model (DEM) generated by The Japan Aerospace Exploration Agency (JAXA) using ALOS PALSAR RTC products [40], slope and aspect layers generated from the DEM using spatial analyst tools in ArcGIS 10.6.1 soil pH with the resolution of 250 m was downloaded from https://soilgrids.org/#!/?layer=TAXNWRB_250m& vector=1 and Uttis cover thematic layer. All bioclimatic variables, topographic layers, and pH were resampled into 20 m spatial resolution to make them compatible with the resolution of Uttis cover. This was done in ArcGIS 10.6.1 with the nearest neighbor resampling technique. For future scenarios (year 2050), we initially selected DEM, slope, aspect, and 19 bioclimatic variables for RCP2.6 (the minimum greenhouse gas emission scenario), and RCP8.5 (the maximum greenhouse gas emission scenario) as adopted by the Intergovernmental Panel on Climate Change (IPCC) in its fifth assessment report (AR5). The climatic variables were also resampled into 20 m to make the variables uniform at one resolution.
Multicollinearity between predictor variables was tested as it can lead to inaccurate prediction by excluding significant explanatory variables [41]. The test was conducted calculating Pearson's Correlation Coefficient (r) to assess the cross-correlation and the one variable from any pair of variables with a cross-correlation coefficient value of > ±0.8 was excluded [42]. Variables were chosen based on biological relevancy to the species. For example, pH was highly correlated with temperature seasonality (BIO4) (r = 0.88), from this pair BIO4 (temperature seasonality) was removed as pH plays a crucial role in cardamom plantation [2]. In addition, the variation inflation factor (VIF) was also used to check collinearity among the variables in R software (version 3.6). Variables with VIF values greater than 10 were excluded for modelling [43]. Out of 24 variables, 8 environmental variables (Uttis cover, pH, slope, aspect, isothermality, maximum temperature of warmest month, minimum temperature of coldest month, and precipitation of wettest month) were selected for the current scenario. For the future scenario of RCP2.6-2050 and for RCP8.5-2050, 6 variables (slope, aspect, isothermality, maximum temperature of warmest month, minimum temperature of coldest month, and precipitation of wettest month) were selected for modelling. As part of the species presence data, a total of 102 occurrence records (GPS coordinates of points) of large cardamom were collected randomly during the field visit conducted in May 2016.

Spatial Modelling and Statistical Analysis
The maximum entropy modelling approach using Maxent software (version 3.3.3k) was applied in this study for predicting habitat suitability of the large cardamom. Maxent is a machine-learning method with a simple and precise mathematical formulation. It uses a maximum entropy algorithm to produce a model that shows the probability of presence of the species that varies from 0 to 1, i.e., from the lowest to the highest probability [21]. We selected 70% of the data for training and the remaining 30% for testing. Area under the ROC (receiver operating characteristic) curve (AUC) was used to evaluate the model performance which ranges from 0 to 1. The curve is plotted with True Positive Rate (sensitivity) at the vertical axis and False Positive Rate (1-specificity) at the horizontal axis. The jackknife was used to evaluate the importance of the variables on the model. The model used logistic format. The final distribution maps have values ranges from 0 to 1 which were grouped into four classes of suitable habitat viz., unsuitable (<0.2), marginally suitable (0.2-0.4), moderately suitable (0.4-0.6), and highly suitable (>0.6). Further, analysis of existing farming area in comparison to current habitat suitability map was performed to understand the gaps and opportunities. In order to assess the change between current and future suitable areas, we quantified the areas (ha) of classes of habitat suitability under different scenarios across different elevation ranges.

Participatory Mapping of Large Cardamom
The spatial location of large cardamom fields in the Taplejung District was recorded on the basis of local people knowledge. Out of 50 VDCs in the district, 49 VDCs were found cultivating ( Figure 2) large cardamom. However, the farming area varies from VDC to VDC. Participatory mapping has given a general overview of farming area in the district.

Uttis (Alnus nepalensis) Cover Mapping and Delineation of Accurate Large Cardamom Map
Spectral separability of major vegetation of the study area using Sentinel-2A image shows that the vegetation are largely separable from each other in NIR and Red Edge bands (Figure 3a). In this study, integration of ancillary data with Sentinel-2A image was applied to map Uttis using OBIA (Figures 3b and 3c). The classification approach in eCognition has helped to use expert knowledge in differentiating trees within agricultural land and forests, resulting in better classification accuracy. During post classification improvement, small tree patches of Uttis within agricultural land were separated. Shadow in high mountain areas was affecting the classification which was improved during post classification using Google Earth images. There are several studies done on forest type classification [44,45] and land cover classification [39], however, these studies do not include separation of tree species within the forest cover classes.

Uttis (Alnus nepalensis) Cover Mapping and Delineation of Accurate Large Cardamom Map
Spectral separability of major vegetation of the study area using Sentinel-2A image shows that the vegetation are largely separable from each other in NIR and Red Edge bands (Figure 3a). In this study, integration of ancillary data with Sentinel-2A image was applied to map Uttis using OBIA (Figure 3b,c). The classification approach in eCognition has helped to use expert knowledge in differentiating trees within agricultural land and forests, resulting in better classification accuracy. During post classification improvement, small tree patches of Uttis within agricultural land were separated. Shadow in high mountain areas was affecting the classification which was improved during post classification using Google Earth images. There are several studies done on forest type classification [44,45] and land cover classification [39], however, these studies do not include separation of tree species within the forest cover classes.

Uttis (Alnus nepalensis) Cover Mapping and Delineation of Accurate Large Cardamom Map
Spectral separability of major vegetation of the study area using Sentinel-2A image shows that the vegetation are largely separable from each other in NIR and Red Edge bands (Figure 3a). In this study, integration of ancillary data with Sentinel-2A image was applied to map Uttis using OBIA (Figures 3b and 3c). The classification approach in eCognition has helped to use expert knowledge in differentiating trees within agricultural land and forests, resulting in better classification accuracy. During post classification improvement, small tree patches of Uttis within agricultural land were separated. Shadow in high mountain areas was affecting the classification which was improved during post classification using Google Earth images. There are several studies done on forest type classification [44,45] and land cover classification [39], however, these studies do not include separation of tree species within the forest cover classes.  Uttis cover in the Taplejung district is about 10,735 ha ( Figure 4a). Overall, accuracy of the classification was 80% with producer's and user's accuracies at 89% and 84%, respectively ( Table 1). The obtained Uttis cover and participatory maps were overplayed to get fine resolution large cardamom farming area (Figure 4b). The overlay produced 5198 ha of existing farming area in the district. Uttis cover in the Taplejung district is about 10,735 ha ( Figure 4a). Overall, accuracy of the classification was 80% with producer's and user's accuracies at 89% and 84%, respectively ( Table 1). The obtained Uttis cover and participatory maps were overplayed to get fine resolution large cardamom farming area (Figure 4b). The overlay produced 5198 ha of existing farming area in the district.  The highest cultivation area is found at the elevation range of 1600-1800 m followed by 1800-2000 m, 1400-1600 m, 1200-1400 m, 2000-2200 m, 1000-1200 m, and 800-1000 m, respectively ( Figure  5a). Although appropriate elevation for large cardamom is 800-2200 m [38], the farming is also observed below 800 m and above 2200 m. However, there are several species of large cardamom which can be farmed based on altitude [2]. The optimal productivity is possible if large cardamom species is chosen based on elevation.
The largest cultivation area is found in the range of 15°-30° slope while the least is found in the range of 45° and above (Figure 5b). The cultivation area is nearly doubled on the slope of 30°-45° compared to the 0°-15° slope. Several literatures [46,47] stated that the best aspect for cardamom cultivation is North and North-East facing slope. The field data collected for the study showed that cardamom is being cultivated in all aspects. There is no such strict constraint in the selection of aspect for farming. The result shows that the farming area is primarily found in South-West followed by West, North, East, South, North-East, South-East, and North-West (Figure 5c). This indicates that large cardamom can be grown in all kinds of aspects. An evaluation on productivity of large cardamom cultivated at various aspects is essential to identify the best aspect for the crop farming.
Out of 50 VDCs of the Taplejung district, Dhungesanghu, Hangdewa, Hangpang, Phurumbu, Sikecha, and Thukinba are the top six VDCs which have the highest cultivation area (Figure 5d). Among those VDCs, Phurumbu and Hangdewa consist of more than 400 ha of farming area whereas Dhungesanghu and Hangpang contain approximately 285 ha to 305 ha of cultivation area. Sikecha and Thukimba hold around 200 ha of cardamom field.   Overall accuracy: 80% The highest cultivation area is found at the elevation range of 1600-1800 m followed by 1800-2000 m, 1400-1600 m, 1200-1400 m, 2000-2200 m, 1000-1200 m, and 800-1000 m, respectively (Figure 5a). Although appropriate elevation for large cardamom is 800-2200 m [38], the farming is also observed below 800 m and above 2200 m. However, there are several species of large cardamom which can be farmed based on altitude [2]. The optimal productivity is possible if large cardamom species is chosen based on elevation.
The largest cultivation area is found in the range of 15 • -30 • slope while the least is found in the range of 45 • and above (Figure 5b). The cultivation area is nearly doubled on the slope of 30 • -45 • compared to the 0 • -15 • slope. Several literatures [46,47] stated that the best aspect for cardamom cultivation is North and North-East facing slope. The field data collected for the study showed that cardamom is being cultivated in all aspects. There is no such strict constraint in the selection of aspect for farming. The result shows that the farming area is primarily found in South-West followed by West, North, East, South, North-East, South-East, and North-West (Figure 5c). This indicates that large cardamom can be grown in all kinds of aspects. An evaluation on productivity of large cardamom cultivated at various aspects is essential to identify the best aspect for the crop farming. Out

Current Suitable Habitat
The Maxent model provided a comprehensive understanding of the distribution of large cardamom (Figure 6a). Currently, the highly suitable area in the district is about 13,679 ha and the moderately suitable habitat is about 27,778 ha. The most suitable habitat for large cardamom was predicted in the southern part of Taplejung and its distribution is almost continuous. Suitability decreases with an increase in altitude. The maxent-predicted model had high accuracy with an AUC value of 0.941 for training data and 0.934 for test data. Jackknife results showed that Bio6 (minimum temperature of coldest month) among the eight variables considered for the model had the highest predictive power. The Bio5 (maximum temperature of warmest month) is the second most important variable followed by Bio13 (precipitation of wettest month), pH, Uttis cover, and Bio3 (Isothermality) (Figure 6b).

Current Suitable Habitat
The Maxent model provided a comprehensive understanding of the distribution of large cardamom (Figure 6a). Currently, the highly suitable area in the district is about 13,679 ha and the moderately suitable habitat is about 27,778 ha. The most suitable habitat for large cardamom was predicted in the southern part of Taplejung and its distribution is almost continuous. Suitability decreases with an increase in altitude. The maxent-predicted model had high accuracy with an AUC value of 0.941 for training data and 0.934 for test data. Jackknife results showed that Bio6 (minimum temperature of coldest month) among the eight variables considered for the model had the highest predictive power. The Bio5 (maximum temperature of warmest month) is the second most important variable followed by Bio13 (precipitation of wettest month), pH, Uttis cover, and Bio3 (Isothermality) (Figure 6b). SDM are influenced by several factors, such as data quality [48,49] and decisions taken during the model fitting [50], sample size [51], multicollinearity [52], and selection of independent variables [53]. Despite these, SDM are increasingly used for various purposes, such as species conservation planning [42] and risk analysis [50]. In this study, we have dealt with some of these issues, such as multicollinearity by removing highly correlated variables, selection of important variables for the species, and considered default settings in Maxent as it provided the best model. Maxent performs best among other modelling methods and even performs better with small sample sizes compared to other modelling methods [30,54,55].

Current Cardamom Cultivation and Habitat Suitability Analysis
Almost 33% of the current farming area is found in the highly suitable class, 44% of the area is found in the moderately suitable class, and nearly 23% of the area is in the marginally suitable class. This indicates that the cultivation has not yet reached its full potential range. Therefore, the farmers need sensitization on potential areas for the farming which would eventually improve productivity (Figure 7). SDM are influenced by several factors, such as data quality [48,49] and decisions taken during the model fitting [50], sample size [51], multicollinearity [52], and selection of independent variables [53]. Despite these, SDM are increasingly used for various purposes, such as species conservation planning [42] and risk analysis [50]. In this study, we have dealt with some of these issues, such as multicollinearity by removing highly correlated variables, selection of important variables for the species, and considered default settings in Maxent as it provided the best model. Maxent performs best among other modelling methods and even performs better with small sample sizes compared to other modelling methods [30,54,55].

Current Cardamom Cultivation and Habitat Suitability Analysis
Almost 33% of the current farming area is found in the highly suitable class, 44% of the area is found in the moderately suitable class, and nearly 23% of the area is in the marginally suitable class. This indicates that the cultivation has not yet reached its full potential range. Therefore, the farmers need sensitization on potential areas for the farming which would eventually improve productivity (Figure 7).
The Uttis map, existing cultivated area, and current habitat suitability of large cardamom provided the means to look at the gaps and opportunities for cardamom cultivation at the VDC level. Figure 7 demonstrates that the current cultivation area of large cardamom is much less in most of the VDCs, though there are more Uttis areas and highly suitable areas. The cardamom cultivated area and Uttis area in Ambegudin, Hangdewa, Phurumbu, and Thukimba VDCs are much closer, whereas current highly suitable areas in Hangdewa, Ikhabu, Nalbu, Phurumbu, and Tellok VDCs are less than the existing cultivated area. The statistics describing gaps and opportunities for cardamom crop cultivation will support management plans across the administrative area.
Suitability maps help in selecting the areas for better success in cultivating any particular crop. However, productivity of a crop is not solely dependent on suitable site but also heavily depend on

Suitable Habitat under Climate Change Scenerio
The predicted future habitat suitability for large cardamom in the district under RCP2.6-2050 and RCP8.5-2050 is shown in Figure 8  The Uttis map, existing cultivated area, and current habitat suitability of large cardamom provided the means to look at the gaps and opportunities for cardamom cultivation at the VDC level. Figure 7 demonstrates that the current cultivation area of large cardamom is much less in most of the VDCs, though there are more Uttis areas and highly suitable areas. The cardamom cultivated area and Uttis area in Ambegudin, Hangdewa, Phurumbu, and Thukimba VDCs are much closer, whereas current highly suitable areas in Hangdewa, Ikhabu, Nalbu, Phurumbu, and Tellok VDCs are less than the existing cultivated area. The statistics describing gaps and opportunities for cardamom crop cultivation will support management plans across the administrative area.
Suitability maps help in selecting the areas for better success in cultivating any particular crop. However, productivity of a crop is not solely dependent on suitable site but also heavily depend on many other factors, such as infrastructure and investments (like irrigation, fertilization, and farm management), biological (insect, pests, and disease), and governance (farmers' training on growth management). A map on Large Cardamom suitability classes and crop yield data, of year 2011, across the VDCs of Taplejung District is given in the Supplementary Materials S1 and S2.

Suitable Habitat under Climate Change Scenerio
The predicted future habitat suitability for large cardamom in the district under RCP2.6-2050 and RCP8.5-2050 is shown in Figure 8. the VDCs of Taplejung District is given in the supplementary material S1-S2.

Suitable Habitat under Climate Change Scenerio
The predicted future habitat suitability for large cardamom in the district under RCP2.6-2050 and RCP8.5-2050 is shown in Figure 8.

Projected Changes in the Suitable Habitat Area
Compared to the currently suitable areas, the total areas tend to decrease in both the future scenarios. However, there is no major difference in total areas of the highly suitable class (2%) in RCP 2.6 and a relatively high change in RCP 8.5 (12%). Elevation-wise, we can see there is loss of highly suitable area from the current scenario to both future scenarios, which for the lower is 80% and in higher elevation it is 94% ( Table 2). The area tends to maintain a nearly steady-state at elevation range 1000-2000 m, which is ± 0.2%. In conclusion, total potentially suitable areas maintain a nearly steadystate under medium elevation range (1000-2000 m) for all scenarios, whereas the areas are decreasing under future scenarios for the low (500-1000 m) and high elevation (2000-3000 m) ( Table 2).

Projected Changes in the Suitable Habitat Area
Compared to the currently suitable areas, the total areas tend to decrease in both the future scenarios. However, there is no major difference in total areas of the highly suitable class (2%) in RCP 2.6 and a relatively high change in RCP 8.5 (12%). Elevation-wise, we can see there is loss of highly suitable area from the current scenario to both future scenarios, which for the lower is 80% and in higher elevation it is 94% ( Table 2). The area tends to maintain a nearly steady-state at elevation range 1000-2000 m, which is ± 0.2%. In conclusion, total potentially suitable areas maintain a nearly steady-state under medium elevation range (1000-2000 m) for all scenarios, whereas the areas are decreasing under future scenarios for the low (500-1000 m) and high elevation (2000-3000 m) ( Table 2). This is the first study to explore current crop cultivation as well as the climate change impacts on the potential distribution of large cardamom crop. Models with AUC values more than 0.75 are considered robust, acceptable, and potentially useful for ecological niche model interpretation [56,57]. Our models obtained AUC values greater than 0.9 which shows our models are satisfactory. The spatial distribution trend under climate change may vary with species causing shifts, contraction, or expansions [58][59][60][61]. Our prediction showed suitable habitat for the crop is diminishing under future conditions compared to the current scenario. Future changes in suitable habitat area are closely related to temperature and precipitation, which are steadily changing in Nepal over time. It is really important to understand the dynamics of these environmental variables. Stations data showed that the maximum temperature in Nepal increased at the rate of 0.056 • C/year for the period of 1975-2014 [62]. The rate was recorded higher at stations placed at higher altitude, whereas very wet and extremely wet days are diminishing significantly in the northern districts of Nepal [63]. Climate projections for temperature developed for entire Nepal by the Ministry of Forest and Environment (MoFE) [63] showed that the temperature trend is higher in high mountains than other regions in both medium-(2016-2045 for RCP 4.5 and RCP 8.5) and long-term (2036-2065 for RCP 4.5 and RCP 8.5) scenarios. The study found central and western parts of the country will be wetter than the eastern region despite the fact that average annual precipitation change tends to increase by 2.1% in the medium-term period and increase by 7.9% in the long-term period. As large cardamom needs cool and humid climatic conditions [39], most of the current habitat will be negatively affected in future considering the climate projection and data recorded in the past. Particularly, lower altitude would not be favorable for the crop as the temperature tends to reach higher degrees in these elevations. Although annual precipitation change increases in the country, warmer temperature in High Mountain would lead to it becoming drier as it would stimulate more evapotranspiration [58]. These projected drier conditions will negatively influence the habitat suitability of cardamom in the study area. Further, Gudade et al. [64] observed that declining soil nutrition and soil moisture are key drivers in reduction in crop productivity both in India and Nepal.

Conclusions
Mapping understory crops and respective habitat conditions is challenging particularly in a complex mountain terrain. The study demonstrated using a combination of expert knowledge-based mapping and species modelling based on bio-climatic factors to understand crop dynamics in a complex mountain terrain. The work applied object-based image analysis on freely available high-resolution satellite data of Sentinel-2A to identify Uttis tree cover in the study area. Use of textures and related indices made it possible to identify a particular tree species within the land cover map. The methodology developed in this study can be used for Uttis cover mapping in other areas of Nepal related to cardamom mapping. This is the first study to demonstrate potential distribution of large cardamom in the region and the impact of climate change on the distribution. The study showed that crop habitat distribution patterns of understory crop could be modelled using Maxent and coupled with expert knowledge-based mapping to understand the environmental constraints and dynamics. The model under the current scenario showed potentially suitable areas for large cardamom, which has created scope to fill the gaps between farmers' practice in selection of farming land, and planners to develop sustainable land use planning and management. Current cultivation area is shown the highest in the medium suitability class. Farmers can change the farming area on high suitability area. The planning agencies can use the cardamom map and other information produced in this work for land use planning and crop management to optimize the crop sown area. Models under future scenarios showed that the overall suitable habitat is shrinking to mid-hill elevations (1000-2000 m). Therefore, in response to the projected loss of potential distribution, appropriate spatial planning for crop management and livelihood strategies needs to be developed for large cardamom farmers.