Abstract
Introduction:
Pollinators are vital for rural ecosystems and agriculture, yet their diversity is declining, threatening biodiversity and ecosystem stability.
Method:
This study evaluates pollinator habitats in Pukou District, Nanjing, using an optimized Maximum Entropy (Maxent) model and landscape pattern analysis to map spatial distribution and assess village-level habitat suitability. By combining field survey data with multi-source database data, we optimized Maxent parameters and applied entropy-based weighting with K-means clustering to evaluate 66 villages.
Results:
Results reveal a three-tiered gradient pattern (mountainforest core, agricultural transition, urban fringe), with 19 high-, 28 moderate-, and 19 low-suitability villages. Natural–artificial edges form synergistic hotspots, while roads support 69%–71% of habitats, and settlements suppress 33%–47%. High-suitability villages benefit from low patch density and high heterogeneity, unlike low-suitability villages with intensive agriculture. We propose core protection, ecological corridors, and community-based conservation, offering a framework for rural biodiversity management in urbanizing regions.
1 Introduction
Pollinator diversity is significant for rural biodiversity and is an important resource for sustainable agriculture. Currently, more than 87.5% of the pollination of angiosperms, more than 70% of the pollination of crops, and more than 35% of the pollination of food crops rely on pollinating insects (). Rural areas are hotspots for biodiversity and home to a rich variety of pollinators. However, with the advance of modernization, changes in land use (; ), intensification of agricultural activities (), and shifts in the lifestyles of rural residents () have destroyed the habitats of pollinators. Studies in rural-urban gradients show that urbanization reduces functional diversity of different species of bees (). Therefore, it is urgent to conduct research on the precise identification and graded protection of pollinator habitats at the village level in the background of rapid urbanization to coordinate ecological security and urban development.
The conservation and enhancement of pollinator habitats are increasingly becoming a focal point in international research. Regional-scale research emphasizes the identification of potential habitats () and assessment of existing habitat types () in order to protect and restore pollinator habitats with higher landscape connectivity. Landscape-scale research concentrates on the spatial distribution patterns of habitats (), enhancing landscape heterogeneity and connectivity (), with an emphasis on the rational allocation of land use types, land cover (), and plant species (). Patch-scale research aims to enhance the quality of the monotypic habitat () and constructing microhabitat networks ().
Although considerable research has addressed habitat conservation for pollinating insects, several limitations persist. The identification of potential pollinator habitats currently relies primarily on remote sensing imagery and field investigations (; ), as the main methods of habitat assessment. While these methods provide data with reasonable accuracy, they are unable to effectively integrate the diverse habitat requirements of pollinators or support continuous assessment across multiple spatial scales, owing to constraints in image data resolution and the extensive expertise required for field surveys. Moreover, while existing studies suggest that landscape heterogeneity and connectivity determine the movement, genetic exchange, and resource accessibility of pollinator insects, and that patch-scale measures such as establishing wildflower strips () and hedgerows () can effectively enhance local pollinator abundance and diversity, most research focuses solely on the proportion of natural habitat without accurately quantifying the gradual impacts of human activities on pollinator habitats. Additionally, these studies fail to optimize the spatial allocation of conservation measures to maximize their efficacy.
To address these challenges, this study examines Hymenoptera (HYM), Lepidoptera (LEP), and Diptera (DIP) as the primary pollinator orders, distinguished by their high flower-visitation rates (). Their scales or body hair facilitate pollen transfer during flower visits, delivering effective pollination services and enhancing fruit set rates (; ). While Coleoptera (beetles) can contribute to pollination in some systems, they were excluded due to lower flower-visitation frequencies and pollination efficiency () compared to HYM, LEP, and DIP in rural Chinese landscapes.
This study aims to: (1) map the spatial distribution of pollinator habitats in rural areas using an optimized MaxEnt model; (2) assess village-level habitat suitability through landscape pattern analysis and entropy-based weighting; and (3) propose targeted conservation strategies for biodiversity management in urbanizing regions. By integrating field surveys with multi-source data, we provide a framework for rural pollinator conservation, offering theoretical support for ecological planning in rapidly urbanizing areas like China.
2 Materials and methods
2.1 Selection of typical areas
The rural areas of Pukou District, Nanjing, are situated in the northwestern and western regions of the district, encompassing 178 natural village sites, covering a total area of 589.27 km2 (Figure 1). Within the study area, substantial areas of farmland, forest, and other habitat types remain, which are essential for pollinator habitat conservation. Moreover, the study area includes key ecological nodes, such as Laoshan, Guanyinshan, and Lvshuiwan, which are vital for conserving pollinator diversity in Pukou District’s rural areas. The Laoshan Nature Reserve hosts 2,081 species, representing 47.66% of Jiangsu Province’s total species, including 718 insect species. Among them, Lvshuiwan has complex terrain and topography, with a variety of wetland types, providing favorable conditions for the growth of numerous plants and animals.
FIGURE 1
Pukou District, a core component of Nanjing Jiangbei New Area, faces a conflict between rapid urbanization and ecological conservation. According to the Nanjing Pukou District Territorial Spatial Plan (2021–2035) and related policies, major industrial projects in Jiangbei New Area have driven the expansion of construction land, with 61.9% occupying arable land, encroaching on pollinator habitats. Concurrently, modern urban development has intensified regional pollution and disturbances, exacerbating pollinator habitat fragmentation, thereby heightening the vulnerability of ecosystems.
Thus, selecting Pukou District as the study area not only captures the characteristic features of rural pollinator habitats but also serves as a critical case study for biodiversity conservation under intense urbanization.
Pollinator data were sourced from field surveys conducted in spring 2025 in Pukou District’s rural areas using the cross-sectional counting method (), supplemented by the Global Biodiversity Information Facility (GBIF, https://www.gbif.org/, accessed February 2025), the Integrated Insect Specimen Database (https://insectaintegration.com/, accessed February 2025), and literature reviews, resulting in 300 records with precise geographic coordinates.
Habitat information was compiled from the Chinese Animal Theme Database (http://www.zoology. csdb.cn, accessed February 2025) and Fauna Sinica Insecta, summarizing the habitat characteristics of pollinators in Pukou District’s rural areas. Ultimately, habitat features for 37 families and 313 pollinator species were compiled, with key attributes extracted, including host plants, suitable temperature ranges, and altitude (Table 1).
TABLE 1
| Order | Family | Habitat keywords |
|---|---|---|
| HYM | Apidae, Megachilidae, Scoliidae, etc. | Interview plants, temperature, forest canopy closure, terrain |
| LEP | Pieridae, Papilionidae, Lycaenidae, etc. | Temperature, host plants, terrain, vegetation structure, forest canopy closure, water source |
| DIP | Syrphidae, Bombyliidae, etc. | Larval parasitism, temperature, forest canopy closure, water source |
Habitat characteristics of pollinator species in Pukou district.
Most crops serve as targets for pollinators. Syrphidae larvae prey on aphids, thrips, and other agricultural pests. Most LEP insects congregate in rice paddies and wetlands near small streams to access water. LEP host plants exhibit varied environmental requirements, necessitating consideration of factors such as forest canopy closure and terrain. Given the varied flight abilities of pollinators, habitat proximity is a key factor influencing their distribution. Concurrently, human activities significantly impact pollinators. Although all three insect orders are influenced by temperature, climate change was deemed insignificant due to the study area’s limited spatial extent, and thus excluded as an environmental factor. In summary, 13 environmental factors were selected (Table 2).
TABLE 2
| Category | No. | Environmental factors | Source | HYM PC/PI (%) | LEP PC/PI (%) | DIP PC/PI% |
|---|---|---|---|---|---|---|
| Natural environmental factors | E1 | DEM/m | ASTER GDEM (Version 1) and GDEM (Version 3) datasets of Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences | 2.6/8.5 | 1.6/7.1 | 1.6/2.3 |
| E2 | NDVI | National Ecosystem Science Data Center, National Science & Technology Infrastructure of China | 7.6/10.6 | 4.3/3.4 | 0/0.1 | |
| E3 | slope | Calculated using DEM data | 1.5/1 | 3.2/3.1 | 0.6/2.7 | |
| E4 | slope direction | Calculated using DEM data | 1.7/2.8 | 1.6/0.4 | 3.7/11.4 | |
| Habitat factors | E5 | sparse forest/m | Calculated the Euclidean distance between sparse forests and each grid in LUCC | 1.6/0.9 | 6.4/11.5 | 4.2/2.9 |
| E6 | dry land/m | Calculated the Euclidean distance between dry land and each grid in LUCC | 5.4/6.2 | 2.4/5.6 | 31.6/3.6 | |
| E7 | dense forest/m | Calculated the Euclidean distance between dense forests and each grid in LUCC | 12.1/9.5 | 21.1/21 | 20.7/12.5 | |
| E8 | river/m | Calculated the Euclidean distance between the OSM river vector data and each grid | 10.5/4.6 | 15.2/11.1 | 6.9/12.8 | |
| E9 | ponds, lakes/m | Calculated the Euclidean distance between the OSM ponds, lakes vector data and each grid | 2.9/6.6 | 10.8/11.5 | 3.1/11.6 | |
| E10 | paddy field/m | Calculated the Euclidean distance between paddy field and each grid in LUCC | 15.5 | /1.5 | 8.2/6.5 | |
| Artificial interference factors | E11 | towns/m | Calculated the Euclidean distance between towns and each grid in LUCC | 25.4 | /41.8 | 14.9/10.6 |
| E12 | rural settlements/m | Calculated the Euclidean distance between rural settlement and each grid in LUCC | 10.7/2.5 | 4.9/4.4 | 1.4/2.9 | |
| E13 | road/m | Calculated the Euclidean distance between the OSM road vector data and each grid | 2.5/3.5 | 5.6/3.8 | 4.7/6 |
Environmental factors and their contributions to Maxent modeling in Pukou district.
E2 and E3 were excluded due to high correlation () and low contribution rate.
PC (Percent contribution) indicates the direct contribution percentage of a factor to the model prediction; PI (Permutation importance) indicates the percentage impact on model accuracy after randomly replacing the factor values.
Terrain data were sourced from the Geospatial Data Cloud, Computer Network Information Center, Chinese Academy of Sciences (http://www.gscloud.cn). The Normalized Difference Vegetation Index (NDVI) was obtained from the National Ecosystem Science Data Center, National Science & Technology Infrastructure of China (http://www.nesdc.org.cn). Land cover data (LUCC) were sourced from the Resource and Environment Science Data Platform (https://www.resdc.cn/Default.aspx). Road network and river data were sourced from OpenStreetMap (https://www.openstreetmap.org/). All data were collected for 2020, with distance factors calculated using Euclidean distances to land use types. Data processing and analysis were conducted using ArcGIS 10.8, with all raster files standardized to a 30 m 30 m pixel resolution.
2.2 Research methods
2.2.1 Data and model Parameter optimization
1. Sample distribution and data cleaning: Approximately 45% of records from online databases were concentrated in the Laoshan forest area, with fewer records from farmland and urban fringe areas. This is likely due to Laoshan’s high ecological integrity and the convenience of field surveys, resulting in sampling bias. This uneven distribution could cause the Maxent model to overfit the Laoshan habitat, thereby reducing prediction accuracy in other areas. To address this issue, the “Trim duplicate occurrences” function in ENMTools was applied to retain only one record per 30 m 30 m grid, minimizing spatial autocorrelation. After cleaning, 82 valid records for HYM, 151 for LEP, and 46 for DIP were obtained.
2. Environmental factor screening: Percent Contribution (PC) and Permutation Importance (PI) were used in the Maxent model to evaluate the influence of ecological factors on species distribution (; ). Initially, PC and PI values for 13 environmental factors were derived from a preliminary experiment. Subsequently, the ENMTools software package was employed to conduct correlation analysis on the 13 environmental factors. If the correlation coefficient between any two environmental factors was, factors with lower PC and PI values were excluded.
3. Model parameter optimization: Biological species exhibit diverse physiological mechanisms and ecological niche adaptations shaped by evolution. Using generic algorithms with preset parameters for modeling often limits applicability. Adjusting Feature Combination (FC) and Regularization Multiplier (RM) significantly enhances the Maxent model’s accuracy (). This study optimized these parameters using the ENMeval package in R. The RM was set within the range [0.5, 4]. Literature and sample complexity suggest that a lower RM (0.5–1) is suitable for simple models to prevent overfitting, whereas a higher RM (3–4) is appropriate for complex landscapes to enhance generalization (). Feature combinations comprised six types (L, LQ, LQP, H, LQH, LQHP), spanning simple to complex features, and were tailored to the ecological niche differences among insect orders ().The ENMevaluate function was used to conduct a full permutation analysis of 48 parameter combinations, with optimal parameters selected based on the minimum delta AICc criterion.
2.2.2 Species habitat extraction
Filtered species data and environmental factors were imported into the model, optimal parameters were applied, and model simulation training was conducted. Pollinator distribution data were divided into a 3:1 ratio, with 75% allocated to a modeling subset and 25% to an independent validation set (). Optimized RM and FC parameters were used to perform 10 independent simulations. The Area Under the Curve (AUC) and True Skill Statistic (TSS) were used to assess the accuracy of the model’s predictions (). An AUC value 0.7 and a TSS value 0.5 indicate that the model’s predictions are reliable. The Maxent model formula is presented below:
is an environmental variable; is the model coefficient; is the probability of species occurrence; is the expected distribution; is the normalization constant.
2.2.3 Suitable area classification and village suitability evaluation
In the model’s predictions, the distribution probability (p) of pollinators ranged from 0 to 1. Using ArcGIS 10.8, distribution maps of suitable habitats for the three insect orders were generated. The SDMtools toolbox was used to convert suitable habitats for the three pollinator orders into binary files, with areas having a distribution probability >0.3 designated as suitable habitats.
The shape index (; ), landscape fragmentation index (; ), and habitat advantage (; ) served as habitat patch evaluation indicators for each village. Furthermore, Maskell () noted that changes in land use and landscape patterns can impact pollinator habitats, with landscape fragmentation, homogenization (; ), and single land use types (; ) adversely affecting pollinator habitats. Consequently, patch density, CONTAG index, and Shannon diversity index were selected as evaluation indicators for village landscape patterns based on land use data (Table 3).
TABLE 3
| Scale | Primary indicator | Secondary indicator | Meaning | Normalization |
|---|---|---|---|---|
| Edge complexity | FRAC MN | 1 indicates a square; the larger the value, the more complex the shape of the patch | positive | |
| ED | Unit area boundary length; the larger the value, the more complex the shape of the patch | positive | ||
| Habitat patch | Fragmentation degree | PD | The number of patches per unit area; the higher the value, the finer the patch segmentation | negative |
| 1/AREA MN | The higher the value, the greater the degree of fragmentation | negative | ||
| Dominance | LPI | The ratio of the area of the largest patch to the total area of the landscape; the higher the value, the more obvious the patch dominance | positive | |
| PLAND | Percentage of habitat patch area; the higher the value, the more obvious the patch dominance | positive | ||
| Village landscape pattern | PD | Reflecting the degree of landscape fragmentation, high PD values indicate small and scattered patches | negative | |
| CONTAG | Reflecting landscape aggregation, low CONTAG values indicate patch dispersion | positive | ||
| SHDI | Reflecting landscape heterogeneity, high SHDI values correspond to diverse land use types | positive |
Evaluation indicator system.
1/AREA MN is obtained by calculating AREA MN from fragstats and then applying excel to calculate the inverse.
Indicator weights were calculated using the Entropy Weight Method, which assigns weights based on information entropy, minimizing human bias and enhancing the scientific rigor of the evaluation (). Subsequently, both indicator types were normalized in Excel, and suitability values for the three pollinator orders were calculated using weighted summation. The arithmetic mean of the suitability values for the three pollinator orders was calculated to derive the comprehensive suitability index for the 66 villages. Finally, ArcGIS 10.8 was used to visualize the 66 administrative villages/communities, classified into three suitability categories—low, medium, and high—using the Jenks natural breaks method. The normalization formula is presented below:
In Formulas 2, 3, X1 is the normalized value of the pattern indicator, X2 is the inverse value of the pattern indicator, Xt is the actual value of the pattern indicator, Xmax is the maximum value of the pattern indicator, and Xmin is the minimum value of the pattern indicator.
The research methodology centered on optimizing the Maxent model, improving habitat distribution prediction accuracy through parameter optimization and multi-source data integration Landscape pattern analysis and entropy-based methods were integrated to thoroughly assess village suitability, ensuring result objectivity and comparability (Figure 2).
FIGURE 2
3 Results
3.1 Data and model optimization results
Using ENMTools, Pearson correlation coefficients were calculated for 13 environmental factors. The results revealed that certain factors were highly correlated, with correlation coefficients of r = 0.834 for E3 and E1, and r = 0.804 for E2 and E7. Based on the Percent Contribution (PC) and Permutation Importance (PI) values from the preliminary experiment (Table2), factors E3 and E2 were excluded. Ultimately, 11 environmental factors were retained for modeling.
Using the ENMeval package, optimal parameters were determined for the three pollinator orders. The results indicated that the optimal parameters were: RM = 3 and FC = LQ for HYM; RM = 3 and FC = LQHPT for LEP; and RM = 2.5 and FC = LQHP for DIP (Figure 3). Optimized parameters were applied to the Maxent model, yielding improved accuracy compared to default settings. The average AUC values were 0.793 0.023 (TSS = 0.53) for HYM, 0.803 0.016 (TSS = 0.56) for LEP, and 0.853 0.028 (TSS = 0.55) for DIP, validating the reliability of the predictions.
FIGURE 3
3.2 Spatial distribution of species habitats
Maxent model results (Figure 4) indicate that pollinator habitats are primarily distributed within the Laoshan Mountain forest area, extending gradually to surrounding farmland and river habitats. Habitat area varies by pollinator order, with DIP exhibiting the smallest and LEP the largest. (1) HYM habitats are centered on the Laoshan forest area, with dispersed patches across rural areas in the region. The northern region, dominated by high-standard farmland, exhibits reduced ecological carrying capacity for HYM due to habitat homogenization. Conversely, the southwestern economic forest area supports HYM habitats due to its diverse vertical vegetation structure. (2) LEP habitats are primarily concentrated in the Laoshan forest area, with secondary aggregation centers in the northern Yongning-Tangquan ecological zone and southwestern village clusters. This distribution reflects ecological niche differentiation among LEP families: Pieridae predominantly occupy open agricultural areas, Papilionidae favor forest ecosystems with complex terrain and dense vegetation, and Lycaenidae are concentrated in grasslands with robust herbaceous layers and marginal habitats. The study area’s multidimensional habitat gradients effectively sustain LEP community diversity. (3) DIP habitats are primarily located in Laoshan Mountain forests, with scattered patches in rural areas and communities to the north and west. This distribution is attributed to the saprophytic nature of DIP larvae, as forests provide optimal saprophytic substrates and microhabitats for their development.
FIGURE 4
3.3 Habitat patches and village landscape patterns
Landscape pattern characteristics of the 66 villages were quantitatively analyzed (Figure 5). The southwestern region displays high fragmentation, with the Shannon Diversity Index (SHDI) exhibiting a concentric pattern of low values in the central core and moderate values in the outer transition zone. The CONTAG index is inversely correlated with SHDI. In the northern region, fragmentation increases along the north–south axis, while SHDI decreases, with low CONTAG values indicating landscape dispersion in the paddy field–fish pond system. The central area shows low fragmentation and high CONTAG, but single land cover results in low SHDI. The southeastern region, however, exhibits both low fragmentation and high diversity.
FIGURE 5
Six landscape indices for habitat patches of the three pollinator orders were calculated across the study area and normalized (Figure 6). LEP habitat patches exhibit greater habitat advantage than those of the other two orders, corresponding to their larger suitable habitat area. LEP habitats also demonstrate higher integrity, with elevated normalized PD and 1/AREA_MN values. HYM habitats show reduced edge complexity and greater patch fragmentation. DIP habitat patches exhibit the lowest dominance among the three orders, yet display low landscape fragmentation, with patches concentrated in the central mountainous areas and surrounding regions.
FIGURE 6
PLAND, PD, LPI, ED, 1/AREA_MN, and FRAC_MN were assigned weights of 0.111, 0.147, 0.150, 0.405, 0.100, and 0.087, respectively. Higher normalized PD and 1/AREA_MN values indicate greater habitat integrity. Edge complexity, habitat integrity, and habitat advantage for the three insect orders were derived through weighted calculations. The results are presented below:
1. HYM core habitats are concentrated in the central mountainous forest area, exhibiting significant ecological advantages and structural integrity (Figure 7). These habitats are characterized by high vegetation cover, well-developed vertical structure, and low edge complexity, forming a stable ecological unit. The southwestern sub-core area displays heterogeneous habitat quality with variable dominance and integrity. The northern area exhibits a fragmented habitat mosaic. Although small continuous habitats exist, the overall pattern limits stable population development.
2. LEP central core habitats demonstrate high ecological advantages, robust structural integrity, and minimal edge effects, serving as primary breeding and dispersal sources (Figure 8). Secondary core areas display a heterogeneous habitat quality gradient, with local patches exhibiting high dominance and integrity. However, at a broader scale, agricultural interface expansion results in habitat patches with fragmented to semi-continuous characteristics, with moderate dominance and integrity.
3. DIP habitats form a significant advantage center in the central region, reflecting continuous high-quality habitat units (Figure 9). Northeastern patches exhibited high edge complexity, whereas southwestern patches demonstrated relatively high integrity. Overall, DIP habitat patches display a “core aggregation–edge fragmentation” pattern.
FIGURE 7
FIGURE 8
FIGURE 9
3.4 Suitability evaluation
The weights for CONTAG, SHDI, and PD were 0.170, 0.214, and 0.616, respectively. The arithmetic mean of the suitability values for the three insect orders was calculated to generate the final suitability map (
Figure 10).
1. High-suitability villages (0.612–0.713): 19 administrative units, including Chalukou and Huaguang. These villages are primarily located around mountains and forests or along the Lvshuiwan wetlands of the Yangtze River. Forest cover exceeds 20%, with undulating terrain and dispersed settlements. The water system comprises primarily scattered ponds and ditches. Landscape pattern analysis reveals low patch density and high SHDI, indicating high habitat heterogeneity and ecological integrity, which support abundant habitats and breeding resources for pollinators. Southeastern villages benefit from robust shoreline protection policies, further enhancing habitat stability.
2. Moderate-suitability villages (0.518–0.612): 28 administrative units, including Daqiao and Yongning. Using the factoextra package in R, K-means clustering was applied to moderate-suitability villages/communities, yielding optimal results with four categories and a silhouette coefficient of 0.42 (Figure 11). The first category, comprising eight administrative units, exhibits a complex landscape pattern with high heterogeneity and connectivity, but structural deficiencies in pollinator habitats limit insect dispersal. The second category, comprising seven administrative units, displays “dual degradation”, with both the landscape pattern and pollinator habitats in a low-functioning state. The third category, comprising six administrative units, demonstrates “landscape–habitat synergy optimization”. Its landscape pattern creates a synergistic effect with the pollinator habitat system, offering strong potential for improvement. The fourth category, comprising seven administrative units, has a moderately functional landscape pattern but exhibits robust habitat patch patterns, indicating functional enhancements through optimized local ecological processes.
3. Low-suitability villages (0.353–0.518): 19 administrative units, including Liugong and Wujiang. These villages are characterized by intensive agriculture and plantations, with high patch density, low SHDI, severe habitat homogenization, and low suitability for the three insect orders. Intensive agriculture has resulted in insufficient habitats, constraining pollinator survival and reproduction.
FIGURE 10
FIGURE 11
The findings revealed a “three-tiered gradient pattern” of pollinator habitats in Pukou District, synergistic hotspots at natural–artificial interfaces, and differential impacts of roads and settlements on habitat distribution. These findings demonstrate pollinators’ ecological adaptability to habitat heterogeneity and underscore the critical role of human activities in shaping habitat patterns. Future research will investigate the ecological mechanisms driving habitat distribution differentiation and synergy, examine the influence of landscape patterns on insect diversity, and develop targeted conservation and restoration strategies based on village suitability evaluations to enhance rural biodiversity conservation in Pukou District.
4 Discussion
4.1 Optimized modeling enhances habitat prediction accuracy
By optimizing the Maxent model’s parameters, this study markedly enhanced habitat prediction accuracy compared to studies using default settings. This study revealed that diverse RM and FC settings are required for different pollinator orders, reflecting their distinct ecological niches and environmental adaptations. The high AUC and TSS values obtained demonstrate that the optimized model is highly applicable to the study area and serves as a reliable tool for habitat assessment in regions with comparable ecological conditions. However, the limited DIP sample size may constrain prediction reliability, and future studies should include larger sample sizes to enhance model generalization.
4.2 Ecological adaptability reflected in habitat distribution patterns
The findings demonstrate that pollinator habitat distribution follows a “three-tiered gradient pattern” (mountain–forest core, agricultural transition zone, urban fringe), closely tied to the ecological adaptability of different pollinator orders. For instance, LEP’s extensive suitable habitat area reflects its robust adaptability to diverse habitats, particularly those with high vegetation diversity and complex terrain. This aligns with , who reported that LEP species exploit diverse ecological niches due to their varied host plant requirements. Conversely, DIP habitats are primarily concentrated in forest areas, reflecting the saprophytic larval requirements for specific microhabitats. The synergistic pattern at natural–artificial interfaces highlights pollinators’ ecological flexibility in exploiting resources at habitat edges. This finding aligns with , , who noted that habitat edges frequently serve as ecological hotspots for pollinators due to resource concentration. However, the urban fringe exerts a pronounced inhibitory effect on habitat distribution due to fragmentation and pollution, particularly for HYM. Thus, future conservation efforts should prioritize enhancing habitat connectivity at the urban–rural interface to mitigate urbanization’s adverse impacts.
4.3 Habitat distribution synergy
Kernel density analysis (Figure 12) reveals that habitats of the three pollinator orders form synergistic hotspots at natural–artificial interfaces, including Xiangtang Village, the Laoshan–university campus transition zone, Hou Chong Wetland Park, and the Fenghuang Mountain area. The high SHDI index in the hot area suggests that habitat heterogeneity fosters interspecific coexistence. Moreover, Synergistic hotspots arise from resource sharing and interspecific interactions. The diverse landscapes of wetlands and transition zones offer varied food sources, supporting the foraging behavior of HYM and LEP. Water sources boost pollinator activity by sustaining humidity gradients and vegetation diversity. Conversely, DIP rely on saprophytic environments, restricting overlap with HYM and LEP, yet they exhibit Local synergies in ecological transition zones at forest edges.
FIGURE 12
Compared with European agricultural landscape studies (), Pukou District’s synergistic hotspots are predominantly located in areas of high ecological integrity, emphasizing the pivotal role of forests and wetlands over single agricultural areas. This synergy bolsters ecosystem stability and sustains pollination services for agriculture. However, urban expansion threatens connectivity in these hotspots, underscoring the need to prioritize protection of core habitats and ecological corridors.
4.4 Landscape patterns reveal habitat suitability heterogeneity
Landscape pattern analysis reveals pronounced spatial heterogeneity in habitat suitability across the 66 villages in Pukou District. High-suitability villages, primarily located in mountainous and forested areas, exhibit low patch density and high SHDI, indicating that complex landscape patterns support pollinator diversity. This aligns with , who reported that highly heterogeneous landscapes sustain greater species richness. Conversely, low-suitability villages, characterized by intensive agriculture, exhibit high patch density and low SHDI, resulting in habitat homogenization and diminished pollinator suitability. Moderate-suitability villages display varied landscape patterns, with some demonstrating landscape–habitat synergy, while others are constrained by structural deficiencies. K-means clustering results underscore the potential for targeted improvements in these villages, particularly those with high landscape heterogeneity but limited habitat functionality. These findings highlight the critical role of landscape pattern optimization in enhancing habitat suitability and offer a scientific foundation for village-level conservation planning.
4.5 Effects of human disturbance gradients on habitat distribution
This study finds roads and settlements are primary human disturbance factors, given their pronounced spatial gradient effects in Pukou District’s rural landscape. Within a 0–300 m road buffer, habitat proportions for the three pollinator orders were 69%, 65%, and 71%, respectively, suggesting that roadside vegetation strips serve as foraging and migration corridors for HYM and LEP. Conversely, within a 0–300 m settlement buffer, habitats accounted for 47%, 40%, and 33%, reflecting the suppressive impact of intense human activity (Figure 13). Gradient analysis reveals that anthropogenic activities shape pollinator distribution by modifying habitat integrity and resource availability, and that road buffers can function as ecological corridors. These findings offer a scientific foundation for village suitability assessments and conservation strategies.
FIGURE 13
4.6 Protection and restoration recommendations
To realize mutual benefits for pollinator habitat conservation and rural development in Pukou District, we suggest a four-tier strategy of core protection, ecological corridor connection, agricultural optimization, and community participation.
1. Core Protection Areas. For high-suitability villages (e.g., Chalukou, Huaguang), conservation should: (a) designate forests and wetlands as permanent ecological zones, prohibiting road expansion or mining; (b) establish 50–100 m wide vegetation buffers at forest–town edges; and (c) adopt sustainable land-use policies aligned with Nanjing’s 2021–2035 Spatial Plan.
2. Ecological Corridor and Stepping-Stone Habitat Construction. For moderate-suitability villages (e.g., Daqiao, Yongning), restoration should: (a) plant native forbs in 10–20 m wide strips; (b) reduce patch density by merging fragments; (c) convert roads into green corridors with 10–20 m vegetation belts; (d) establish riverine corridors; and (e) spacing wildflower strips 0.five to two km to match pollinator dispersal.
3. Agricultural Optimization. For low-suitability villages (e.g., Liugong, Wujiang), restoration should: (a) reduce pesticide use; (b) plant diverse native species; (c) adopt crop rotation and intercropping; (d) create 20–50 m buffer zones; and (e) monitor non-point source pollution.
4. Community Participation and Policy Support. Community programs should: (a) promote conservation awareness; (b) develop eco-tourism routes, funding farmer-led wildflower maintenance; (c) establish monitoring networks with universities; and (d) integrate corridors into Nanjing’s green infrastructure plans. Landscape-scale planning should optimize road networks and encourage dispersed settlements to reduce fragmentation.
In summary, high-, moderate- and low-suitability areas were proposed as core protection, target areas for quality improvement and pilot areas for restoration, respectively. Based on policy incentives, technological support and public participation, this study offers a scientific basis for the effective protection of pollinators in Pukou District and a methodological guideline for the biodiversity conservation in the rapidly urbanized areas through optimized Maxent modeling and landscape-scale approaches.
5 Conclusion
This study employed an optimized Maxent model to map the habitat distribution of three pollinator orders (HYM, LEP, and DIP) in Pukou District’s rural areas, identifying a “three-tiered gradient pattern” and synergistic effects at natural–artificial interfaces. The findings demonstrate that landscape patterns profoundly influence habitat suitability, with high-suitability villages thriving due to complex landscape patterns, while low-suitability villages experience habitat homogenization from intensive agriculture. These findings offer a scientific foundation for developing targeted conservation strategies, including habitat preservation in high-suitability villages, restoration in moderate-suitability villages, and extensive habitat restoration in low-suitability villages. The findings theoretically validated the applicability of the optimized Maxent model for small-scale habitat prediction, methodologically integrated entropy-based methods and landscape pattern analysis to provide a versatile framework for multi-scale ecological assessments, and practically informed ecological conservation planning in Pukou District and biodiversity management in rapidly urbanizing regions.
However, the limited DIP sample size and the restricted study area scope may constrain the generalizability of the findings. Future research should integrate recent remote sensing data to monitor dynamic habitat changes, validate model predictions with small-scale field surveys, and investigate the integration of Maxent with other ecological models. Furthermore, examining microhabitat characteristics within patches will enhance understanding of pollinator ecological adaptations. This study offers methodological frameworks for rural ecological conservation, applicable to other rapidly urbanizing regions to support global biodiversity conservation goals.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
YY: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Software, Writing – original draft. JQ: Conceptualization, Methodology, Writing – review and editing. XX: Data curation, Writing – review and editing.
Funding
The author(s) declare that no financial support was received for the research and/or publication of this article.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AyersA. C.RehanS. M. (2025). Determining bee community response to urbanization through multi-year monitoring. Urban Ecosyst.28, 177. 10.1007/s11252-025-01797-1
2
BahnV.McGillB. J. (2013). Testing the predictive performance of distribution models. Oikos122, 321–331. 10.1111/j.1600-0706.2012.00299.x
3
Barahona-SegoviaR. M.Gatica-BarriosP.Durán-SanzanaV.Smith-RamírezC. (2023). No wild bees? don’t worry! non-bee flower visitors are still hard at work: the edge effect, landscape, and local characteristics determine taxonomic and functional diversity in apple orchards. Agric. Ecosyst. and Environ.354, 108554. 10.1016/j.agee.2023.108554
4
BartholoméeO.LavorelS. (2019). Disentangling the diversity of definitions for the pollination ecosystem service and associated estimation methods. Ecol. Indic.107, 105576. 10.1016/j.ecolind.2019.105576
5
BegoshA.SmithL. M.McMurryS. T. (2022). Major land use and vegetation influences on potential pollinator communities in the high plains of Texas. J. Insect Conservation26, 231–241. 10.1007/s10841-022-00380-9
6
BergholzK.SittelL.-P.RistowM.JeltschF.WeissL. (2022). Pollinator guilds respond contrastingly at different scales to landscape parameters of land-use intensity. Ecol. Evol.12, e8708. 10.1002/ece3.8708
7
BirkenbachM.StraubF.KieselA.AyasseM.WilfertL.KupplerJ. (2024). Land-use affects pollinator-specific resource availability and pollinator foraging behaviour. Ecol. Evol.14, e11061. 10.1002/ece3.11061
8
BrownK. S.FreitasA. V. L. (2000). Atlantic forest butterflies: indicators for landscape conservation 1. Biotropica32, 934–956. 10.1111/j.1744-7429.2000.tb00631.x
9
BuchoriD.RizaliA.LarasatiA.HidayatP.NgoH.Gemmil-HerrenB. (2019). Natural habitat fragments obscured the distance effect on maintaining the diversity of insect pollinators and crop productivity in tropical agricultural landscapes. Heliyon5, e01425. 10.1016/j.heliyon.2019.e01425
10
CaitanoB.ChavesT. P.DodonovP.DelabieJ. H. C. (2020). Edge effects on insects depend on life history traits: a global meta-analysis. J. Insect Conservation24, 233–240. 10.1007/s10841-020-00227-1
11
CavigliassoP.PhiferC. C.KnowltonJ. L.LicataJ. A.FlaspohlerD. J.WebsterC. R.et al (2022). Influence of landscape composition on wild bee communities: effects of functional landscape heterogeneity. Agric. Ecosyst. and Environ.340, 108150. 10.1016/j.agee.2022.108150
12
ElithJ.PhillipsS. J.HastieT.DudíkM.CheeY. E.YatesC. J. (2011). A statistical explanation of maxent for ecologists. Divers. distributions17, 43–57. 10.1111/j.1472-4642.2010.00725.x
13
GeppertC.HassA.FöldesiR.DonkóB.AkterA.TscharntkeT.et al (2020). Agri-environment schemes enhance pollinator richness and abundance but bumblebee reproduction depends on field size. J. Appl. Ecol.57, 1818–1828. 10.1111/1365-2664.13682
14
GillespieM. A.BaudeM.BiesmeijerJ.BoatmanN.BudgeG. E.CroweA.et al (2024). Flowering plant communities mediate the effects of habitat composition and configuration on wild pollinator communities. Funct. Ecol.38, 2576–2594. 10.1111/1365-2435.14666
15
GomezJ. M.Munoz-PajaresA. J.AbdelazizM.LoriteJ.PerfecttiF. (2014). Evolution of pollination niches and floral divergence in the generalist plant erysimum mediohispanicum. Ann. Bot.113, 237–249. 10.1093/aob/mct186
16
GoodwinE. K.RaderR.Encinas-VisoF.SaundersM. E. (2021). Weather conditions affect the visitation frequency, richness and detectability of insect flower visitors in the australian alpine zone. Environ. Entomol.50, 348–358. 10.1093/ee/nvaa180
17
HarisA.JózanZ.RollerL.ŠimaP.TóthS. (2024). Changes in population densities and species richness of pollinators in the carpathian basin during the last 50 years (hymenoptera, diptera, lepidoptera). Diversity16, 328. 10.3390/d16060328
18
Huang ChunWeiH. C.McDonaldR.SetoK. (2018). The importance of land governance for biodiversity conservation in an era of global urban expansion. Landsc. Urban Plan.173, 44–50. 10.1016/j.landurbplan.2018.01.011
19
KhalafS. M.AlqahtaniM. S.AliM. R.AbdelalimI. T.HodhodM. S. (2024). Using maxent modeling to analyze climate change impacts on pseudomonas syringae van hall, 1904 distribution on the global scale. Heliyon10, e41017. 10.1016/j.heliyon.2024.e41017
20
Kovács-HostyánszkiA.FöldesiR.BáldiA.EndrédiA.JordánF. (2019). The vulnerability of plant-pollinator communities to honeybee decline: a comparative network analysis in different habitat types. Ecol. Indic.97, 35–50. 10.1016/j.ecolind.2018.09.047
21
KremenC.WilliamsN. M.AizenM. A.Gemmill-HerrenB.LeBuhnG.MinckleyR.et al (2007). Pollination and other ecosystem services produced by mobile organisms: a conceptual framework for the effects of land-use change. Ecol. Lett.10, 299–314. 10.1111/j.1461-0248.2007.01018.x
22
LajosK.SamuF.BihalyÁ. D.FülöpD.SárospatakiM. (2021). Landscape structure affects the sunflower visiting frequency of insect pollinators. Sci. Rep.11, 8147. 10.1038/s41598-021-87650-9
23
LarssonM. (2005). Higher pollinator effectiveness by specialist than generalist flower-visitors of unspecialized knautia arvensis (dipsacaceae). Oecologia146, 394–403. 10.1007/s00442-005-0217-y
24
LevéM.BaudryE.Bessa-GomesC. (2019). Domestic gardens as favorable pollinator habitats in impervious landscapes. Sci. Total Environ.647, 420–430. 10.1016/j.scitotenv.2018.07.310
25
LichtenbergE. M.HeiserJ.BaumK. A.NeffJ. L.JhaS. (2024). Pollinators differentially respond to local and landscape grassland features. Insect Conservation Divers.18, 417–428. 10.1111/icad.12816
26
LiuM.WangX.YangC. (2025). Exploring the potential distribution areas of changnienia amoena and its pollinators in China based on maxent and gtwr models. J. Nat. Conservation86, 126946. 10.1016/j.jnc.2025.126946
27
Louis-LucasT.ClavelJ.BortolamiolS.BlancN.GrésillonE.ClauzelC. (2025). Exploring the relationships between landscape connectivity and urban biodiversity: insights from citizen science on pollinators and birds in paris, France. Biol. Conserv.309 (1), 111271. 10.1016/j.biocon.2025.111271
28
MariniL.ÖckingerE.BergmanK.-O.JaukerB.KraussJ.KuussaariM.et al (2014). Contrasting effects of habitat area and connectivity on evenness of pollinator communities. Ecography37, 544–551. 10.1111/j.1600-0587.2013.00369.x
29
MarshallE. J.MoonenA. C. (2002). Field margins in northern europe: their functions and interactions with agriculture. Agric. Ecosyst. and Environ.89, 5–21. 10.1016/s0167-8809(01)00315-2
30
MaskellL. C.CroweA.DunbarM. J.EmmettB.HenrysP.KeithA. M.et al (2013). Exploring the ecological constraints to multiple ecosystem service delivery and biodiversity. J. Appl. Ecol.50, 561–571. 10.1111/1365-2664.12085
31
Mateos-FierroZ.GarrattM. P.FountainM. T.AshbrookK.WestburyD. B. (2023). The potential of wildflower strips to enhance pollination services in sweet cherry orchards grown under polytunnels. J. Appl. Ecol.60, 1044–1055. 10.1111/1365-2664.14394
32
MerowC.SmithM. J.Silander JrJ. A. (2013). A practical guide to maxent for modeling species’ distributions: what it does, and why inputs and settings matter. Ecography36, 1058–1069. 10.1111/j.1600-0587.2013.07872.x
33
MurakoshiN.ItagakiT.OguroM.SakaiS. (2024). Effects of floral display size, local open raceme density, patch size, and distance between patches on pollinator behaviour in salvia nipponica. Sci. Rep.14, 967. 10.1038/s41598-024-51327-w
34
NatuharaY.ImaiC.TakahashiM. (1998). Evaluation of community indices in seasonal assemblages of butterflies (lepidoptera) at different frequency of transect count. Biodivers. and Conservation7, 631–639. 10.1023/a:1008804418188
35
ParmentierL.Van KerckvoordeA.CouckuytJ.van CalsterH.SmaggheG.HaesaertG. (2025). Sinus management: meandering mowing as a novel method to improve pollinator biodiversity and habitat heterogeneity in mesic grasslands. Agric. Ecosyst. and Environ.382, 109478. 10.1016/j.agee.2025.109478
36
PhillipsS. J.DudíkM. (2008). Modeling of species distributions with maxent: new extensions and a comprehensive evaluation. Ecography31, 161–175. 10.1111/j.0906-7590.2008.5203.x
37
PhillipsS. J.AndersonR. P.SchapireR. E. (2006). Maximum entropy modeling of species geographic distributions. Ecol. Model.190, 231–259. 10.1016/j.ecolmodel.2005.03.026
38
PramanikM.DiwakarA. K.DashP.SzaboS.PalI. (2021). Conservation planning of cash crops species (garcinia gummi-gutta) under current and future climate in the western ghats, India. Environ. Dev. Sustain.23, 5345–5370. 10.1007/s10668-020-00819-6
39
RaderschallC. A.LundinO.LindströmS. A.BommarcoR. (2022). Annual flower strips and honeybee hive supplementation differently affect arthropod guilds and ecosystem services in a mass-flowering crop. Agric. Ecosyst. and Environ.326, 107754. 10.1016/j.agee.2021.107754
40
RanalliR.GalimbertiA.LabraM.BiellaP. (2025). Forest spatial configuration and local management influence bee pollinator biodiversity in urban and rural landscapes. J. Environ. Manag.377, 124672. 10.1016/j.jenvman.2025.124672
41
RenP.DidhamR. K.MurphyM. V.ZengD.SiX.DingP. (2023). Forest edges increase pollinator network robustness to extinction with declining area. Nat. Ecol. and Evol.7, 393–404. 10.1038/s41559-022-01973-y
42
SchochK.TschumiM.LutterS.RamseierH.ZinggS. (2022). Competition and facilitation effects of semi-natural habitats drive total insect and pollinator abundance in flower strips. Front. Ecol. Evol.10, 854058. 10.3389/fevo.2022.854058
43
ShiX.MaC.de KrakerJ.GongS.HodgsonJ. A.LuoS.et al (2024). Influence of agricultural intensification on pollinator pesticide exposure, food acquisition and diversity. J. Appl. Ecol.61, 1905–1917. 10.1111/1365-2664.14701
44
ShwartzA.MuratetA.SimonL.JulliardR. (2013). Local and management variables outweigh landscape effects in enhancing the diversity of different taxa in a big metropolis. Biol. Conserv.157, 285–292. 10.1016/j.biocon.2012.09.009
45
SimkinR. D.SetoK. C.McDonaldR. I.JetzW. (2022). Biodiversity impacts and conservation implications of urban land expansion projected to 2050, Proc. Natl. Acad. Sci. U. S. A.119, e2117297119, 10.1073/pnas.2117297119
46
TommasiN.BiellaP.MaggioniD.FallatiL.AgostinettoG.LabraM.et al (2023). Dna metabarcoding unveils the effects of habitat fragmentation on pollinator diversity, plant-pollinator interactions, and pollination efficiency in maldive islands. Mol. Ecol.32, 6394–6404. 10.1111/mec.16537
Appendix A
Appendix A is a supplemental explanation of Figure 11, as adding the names of all the villages in Figure 11 would have made the picture difficult to see, the choice was made to assign a number to each village.the information relating to the villages to which each number corresponds can be viewed in Table A1.
TABLE A1
| Number | name | Composite index of habitat patch | Composite index of village landscape pattern |
|---|---|---|---|
| 1 | Da qiao | 0.695274167 | 0.599019167 |
| 3 | Ban qiao | 0.573669 | 0.511438 |
| 4 | Gao hua | 0.552449167 | 0.4777045 |
| 8 | Xi shan | 0.552449167 | 0.702424 |
| 10 | Tang quan | 0.586693667 | 0.778125 |
| 14 | Luo qiao | 0.554407667 | 0.710823 |
| 21 | Wan long | 0.519828667 | 0.700424 |
| 23 | Gao li | 0.592540833 | 0.734857 |
| 27 | Jiu long | 0.530468333 | 0.654346 |
| 30 | Yong ning | 0.5982705 | 0.734036 |
| 31 | Gao wang | 0.581279 | 0.638604 |
| 32 | Xi ge | 0.519916333 | 0.748575 |
| 33 | Long shan | 0.522916333 | 0.712577 |
| 34 | Da yan | 0.612298167 | 0.6546 |
| 39 | Shan xi | 0.523938333 | 0.646698 |
| 43 | Di shui zhu | 0.556272333 | 0.616311 |
| 45 | Hou wei | 0.5716145 | 0.666227 |
| 50 | Shuang miao | 0.537579167 | 0.572241 |
| 53 | Quan dong | 0.525931667 | 0.702539 |
| 56 | Wang cun | 0.518004833 | 0.713882 |
| 57 | Gao miao | 0.521771167 | 0.616174 |
| 59 | Qin shan | 0.559383833 | 0.762872 |
| 60 | Lian he | 0.592244667 | 0.736619 |
| 61 | Da huang | 0.572153833 | 0.650409 |
| 63 | Wa dian | 0.566870833 | 0.584317 |
| 64 | Nan yi | 0.524158667 | 0.610723 |
Information of villages with moderate suitability.
Summary
Keywords
pollinators, Maxent model, landscape pattern, rural habitat, biodiversity
Citation
Yin Y, Qian J and Xu X (2025) Study on mapping and assessment of pollinator habitats in rural areas, using an optimized Maxent model. Front. Environ. Sci. 13:1668003. doi: 10.3389/fenvs.2025.1668003
Received
17 July 2025
Accepted
02 September 2025
Published
11 September 2025
Volume
13 - 2025
Edited by
Juergen Pilz, University of Klagenfurt, Austria
Reviewed by
Claudia Bita-Nicolae, Institute of Biology Bucharest of the Romanian Academy, Romania
Thomas Tscheulin, University of the Aegean, Greece
Updates
Copyright
© 2025 Yin, Qian and Xu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jing Qian, Jinn0220@njtech.edu.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.