- Research
- Open access
- Published:
Decoding the impact of environmental shifts on snail density dynamics in the Yangtze River basin: a 26-year study
Parasites & Vectors volume 18, Article number: 156 (2025)
Abstract
Background
With the intensification of climate change and human engineering activities, environmental changes have affected schistosome-transmitting snails. This study explored the influence of environmental changes on the evolution of snail populations.
Methods
Data from annual snail surveys and related factors such as hydrology, temperature, vegetation, etc., on nine bottomlands from 1997 to 2022 were collected retrospectively from multiple sources. Interpretable machine learning and the Bayesian spatial-temporal model assessed the relationship between environmental change and snail density.
Results
Between 1997 and 2003, mean snail density was in a high-level fluctuation stage. From 2003 to 2012, it declined significantly from 0.773/0.1 m2 to 0.093/0.1 m2. However, it increased by 27.6% between 2013 (0.098/0.1 m2) and 2022 (0.125/0.1 m2). Since operation of the Three Gorges Dam (TGD) began in 2003, the duration of bottomland flooding decreased from 122 days (1997–2003) to 57 days (2003–2012) and then rebounded in 2012–2022, which was noticeable in the Anhui Section. The ground surface temperature and night light index of the bottomlands increased from 1997 to 2022. After adjusting for confounding factors (e.g. rainfall, temperature, and vegetation), the relative risk (RR) of increased snail density rose with flooding duration of between 20 and 100 days but decreased with flooding duration of > 100 days. Snail density showed an “L”-shaped relationship with the night light index, and the RR of increased snail density was lower at a higher night light index. Compared with bottomlands in the first quartile cluster of ground surface temperature, bottomlands in the second, third, and fourth quartile clusters of ground surface temperature had higher snail density RR values of 1.271 (95% CI 1.082–1.493), 1.302 (95% CI 1.146–1.480), and 1.278 (1.048, 1.559), respectively.
Conclusions
The TGD lowered the water level and flooding duration, which were not conducive to snail population growth. However, over time, the inhibitory effect of the TGD on snails may have been weakening, especially in areas far from the TGD. In recent years, the rebound of snail density may have been related to the rise in water levels and the change in the microenvironment. Establishing an efficient monitoring and response system is crucial for precisely controlling snails.
Graphical Abstract

Background
China has significantly reduced schistosomiasis infections by implementing flexible control strategies at different times and is currently promoting the elimination of schistosomiasis [1, 2]. However, as a natural epidemic disease, the risk of schistosomiasis transmission has existed for a long time because the natural and social environment in the epidemic area has not undergone fundamental changes [3]. The Oncomelania hupensis snail is the only intermediate host involved in the transmission of Schistosoma japonicum [1, 3]. Snail control is an effective measure recommended by WHO to eliminate the hidden dangers of schistosomiasis [4]. However, the area of snail habitats has remained at a high level of 3.6 billion m2 since 2000, with emerging habitats continuing to appear [5].
The dynamics of the snail life cycle are influenced by both natural environmental and human activities. Natural factors mainly include climate change and major flood events, while human activity involves water conservancy projects (e.g. dam construction) [6, 7]. Temperature and precipitation are key limiting factors regarding snail distribution. Studies by Zhou et al. [8] suggest that with global warming, snail habitats may expand northward into the Yangtze River basin, reaching northern Jiangsu and Anhui. Following the 1998 Yangtze River flood event, the areas of re-emerging and newly established snail habitats in the basin were 2.6 and 2.7 times larger, respectively, than in years with normal hydrological conditions [9]. Water conservancy projects, such as dams and irrigation systems, often alter water distribution and create favorable habitats for snails [10, 11]. In the Senegal River basin, evidence suggests that after the construction of the Diama Dam, ecological changes promoted the growth and reproduction of freshwater snails [12]. The proportion of schistosome-transmitting snails among collected specimens increased from 22 to 60% [12]. This chain reaction—dam construction, environmental change, and proliferation of freshwater snails—has drawn significant attention [10,11,12].
The Three Gorges Dam (TGD) is a large-scale hydroelectric power project built in the upper reaches of the Yangtze River and has been in operation since 2003 [13]. Due to the artificial regulation of water resources by the TGD, the downstream water regime has undergone significant changes [13, 14]. Li et al. [15] showed that following the TGD, the ecological environment (such as volume of runoff, sediment volume, and sedimentation rate) has changed in Dongting Lake, causing the occurrence rate of frames with living snails and living snail density to decrease by > 80% from 2003 to 2015. Wu et al. [16] showed that the TGD reduced the risk of flooding and living snail density. Fundamental changes in water environment inevitably impact various aspects such as the microenvironment (vegetation, sunshine light, and temperature) and utilization of the bottomland. These ecological elements are also important components of the survival conditions required for snails, but research on the impact of the water environment and related factors on snail distribution are still inadequate.
Artificial intelligence (AI) and new model algorithms play a significant role in monitoring snail habitats and exploring shifts in snail density in their environment [17]. First, AI can process remote sensing data from satellites, drones, and other geospatial sources in real time to track key environmental variables such as water body, temperature, and land cover [18]. Second, new model algorithms can detect spatial and temporal patterns that may not be apparent using traditional analysis. For instance, Qureshi et al. [19] used machine learning and SHapley Additive exPlanations (SHAP) to analyze the flight trajectory characteristics of mosquitoes on mosquito nets and provide direction for mosquito interventions. Third, AI-driven automated monitoring systems can analyze real-time data from environmental sensors, field surveys, and image recognition techniques to detect changes in disease vectors. Explainable machine learning and Bayesian models are increasingly applied in vector-borne disease analysis because of their flexibility and efficiency in handling complex data. For instance, Villela et al. [20] constructed a Bayesian hierarchical model to estimate the abundance and spatial density of Aedes aegypti, demonstrating the potential of related models in assessing vector populations and their environmental dynamics.
This study aimed to explore the long-term evolutionary characteristics of snail density related to the construction of TGD using explainable machine learning and a Bayesian spatial-temporal model by collecting continuous snail survey data from nine bottomlands from 1997 to 2022. This study could provide a reference for understanding the evolution of snails associated with environmental change and ecological snail control.
Methods
Study area
The Yangtze River basin is characterized by a humid subtropical climate with distinct monsoon influences [21]. Precipitation increases from southwest to northeast across the basin, with annual rainfall ranging from 800 to 2000 mm. The region has small temperature variations and high humidity, with an average annual temperature of 15 °C to 19 °C. The topography is primarily flat, consisting of plains and low hills, interspersed with numerous rivers, lakes, and wetlands. The region is influenced by both monsoon and tropical climate systems, leading to significant water level fluctuations and frequent flooding [22]. These conditions support the basin's rich vegetation and biodiversity, including snail populations.
The Yangtze River, spanning a colossal 6300 km and draining an expansive 1.8 million km2 from the Tibetan Plateau to eastern China, is the largest river in Asia. Situated in the upper reaches of the Yangtze River, the TGD began filling with water in June 2003 [13]. In this study, we focused on the downstream region of the TGD, selecting nine snail-infested sites from Dongting Lake, Poyang Lake, and the Anhui section (Fig. 1). In each region, three study sites were chosen, representing high-, medium-, and low-elevation bottomlands (Additional file 1: Table S1).
Data collection and preprocessing
The snail survey was conducted annually in the nine bottomlands each spring from 1997 to 2022 (Fig. 1; Additional file 1: Table S1). This was performed by using the Chinese traditional method of systematic sampling (0.1 m2-sized frames) [23, 24]. All snails within the square frames were collected and brought to the laboratory. These collected snails were then counted and crushed to determine their vitality [23, 24]. The overall mean density of living snails and living snail density at different elevations were obtained by dividing the number of living snails by the number of snail survey frames.
Six hydrological stations in proximity to the study bottomland within the Dongting Lake area, Poyang Lake area, and Anhui section of the Yangtze River were selected. From 1997 to 2022, daily water level data at 8 a.m. were retrospectively collected from these hydrological stations. The average elevation of each bottomland was derived from elevation grid data (https://www.resdc.cn/). A bottomland was classified as submerged if the difference between the daily water level at the hydrological station and the average elevation of the beach was > 0 [25]. The flooding duration (FD; unit: days) was defined as the number of days that a bottomland was submerged in water during a year.
Vegetation is an important factor affecting snails, which can usually regulate the temperature and humidity of snail habitat and serve as their food. The normalized difference vegetation index (NDVI) was an accurate gauge of vegetation greenness, photosynthetic intensity, metabolic activity, and seasonal and interannual variations. This study employed the maximum composites method [26] on monthly NDVI grid data from the PKU GIMMS NDVI [27], enabling the derivation of annual NDVI grid data between 1997 and 2022.
The nighttime light conditions of the bottomland and its surroundings can reflect the development and utilization of these areas, as well as the infrastructure construction, which indirectly affects the breeding and reproduction of snails [28]. The nighttime light index (NL) used in this study is derived from the corrected DMSP-OLS-like data for mainland China (continuously updated since 1992), which were obtained by integrating DMSP-OLS and SNPP-VIIRS data [29]. The data were sourced from the National Earth System Science Data Center (https://www.geodata.cn/data/).
The sunshine hour (SSH/hour), ground surface temperature (GST/℃), precipitation (Pre/mm), air temperature (Tem/℃), and relative humidity (RH/%) were acquired from the National Meteorological Data Service Platform (http://cdc.cma.gov.cn/). Kriging interpolation techniques were applied to generate raster data [30]. Subsequently, the data for the nine bottomlands were extracted from this interpolated dataset. Extreme temperature conditions were gauged using the average minimum temperature in January (Tmin/℃) and the average maximum temperature in July (Tmax/℃), sourced from the National Earth System Science Data Center (http://www.geodata.cn/).
There may be correlations between the above variables, which may cause collinearity and modeling errors, so we conducted a correlation analysis to exclude variables with high correlations. The Spearman correlation coefficients between the other variables were all < 0.5, with only those between GST and Tem, as well as FD and RH, > 0.7 (Additional file 1: Table S2). To avoid collinearity, Tem and RH were excluded from the analysis. The final variables included in the analysis were FD, NDVI, NL, Tmin, GST, Pre, SSH, and Tmax.
Statistical analysis
Previous theoretical studies have suggested potential differences in water and sediment dynamics during various operational stages of the TGD [31]. Therefore, this study divided the research time into pre-TGD (1997–2002) and post-TGD (the 2003–2012 and 2013–2022 periods). In addition, microenvironments such as temperature and sunshine may also change over time. The Shapiro-Wilk normality test was employed to assess the distribution of environmental factors, which revealed a departure from the normal distribution. The P-values from the Shapiro-Wilk test for all variables were < 0.05 (Additional file 1: Table S3), indicating that the data did not follow a normal distribution. Therefore, non-parametric (Kruskal-Wallis test) methods were used to compare the environmental factors among different stages. Bonferroni correction was applied to mitigate the risk of Type I errors associated with multiple comparisons [32].
This study used the light gradient boosting machine (LightGBM) model with the SHAP method to discern the principal drivers influencing vegetation and amphibious freshwater snails, elucidating their relationships with core ecological factors. LightGBM continuously refines the model by fitting residuals from the current learner and employs the forward distribution algorithm for iterative training [33]. The SHAP method, a model explanation tool rooted in Shapley values, functions as a post hoc explainability scheme [34]. It elucidates the significance of each feature and the direction in which each feature influences the decision-making process of the model [35].
Although SHAP values identified the relationship between individual environmental factors and snail density, these relationships might be confounded by multiple factors, and no statistical tests have been conducted. We built a multi-factor Bayesian space-time model to further examine the snail density and related factors. The general form of the model is given as follows [36, 37]:
In this model, \({y}_{it}\) represents the observed number of snails, and \({e}_{it}\) denotes the expected number of snails. The model assumes that \({y}_{it}\) follows a negative binomial distribution with a mean \({u}_{it}={e}_{it}{\theta }_{it}\), where \({\theta }_{it}\) represents the relative risk (RR) at site \(i\) in year \(t\). The parameter \(r\) accounts for dispersion. \({\beta }_{k}\) is a vector of regression parameters, and \({X}_{itk}\) represents the covariate. \(S({X}_{itq})\) denotes the smoothing spline function terms used to approximate the nonlinear association between environmental variables and snail density. The terms \({u}_{i}\) and \({\varphi }_{t}\) represent the spatially unstructured and temporally structured effects, respectively; both are assume to follow a normal distribution. The spatially structured effect \({v}_{i}\) was modeled using a conditional autoregressive (CAR) process, with the CAR adjacency matrix \(\text{w}\) defined by the K nearest neighbors. The snail densities in the same region were similar, so K was set to 2 in this study. The temporally structured component \({\gamma }_{t}\) was modeled using a first-order autoregression. The spatiotemporal interaction effect \({\delta }_{it}\) was also assumed to follow a normal distribution. Hyperpriors for the precisions of the random effects were assigned as follows: \({\tau }_{u}\), \({\tau }_{\text{v}}\), \({\tau }_{\varphi }\), \({\tau }_{\gamma }\), \({\tau }_{\delta }\). Their prior distributions were \({\tau }_{u}\sim (\text{1,0.0005})\), \({\tau }_{\text{v}}\sim (\text{1,0.0005})\), \({\tau }_{\varphi }\sim (1, 0.00005)\), \({\tau }_{\gamma }\sim (\text{1,0.00005})\), \({\tau }_{\delta }\sim (\text{1,0.00005})\).
The models were fitted using the Integrated Nested Laplace Approximation method and evaluated based on the deviance information criterion (DIC). The model with the lowest DIC score was chosen as the final model for this study. The analyses were conducted in R 4.4.1 (R core team, Vienna, Austria) and ArcGIS 10.2 (Esri, Redlands, CA, USA). We considered a two-sided P-value < 0.05 as statistically significant.
Results
The characteristics of snail density
Between 1997 and 2003, the snail density was in a high-level fluctuation stage, and the overall snail density dropped from 1.357 in 1998 to 0.610 in 2000 and then rebounded to 0.773 in 2003 (Fig. 2A). During the 1st decade operation of the TGD (2003–2012), the snail density in high-elevatione bottomlands began to decrease in 2003, and the snail density in medium-elevation bottomlands started to decrease in 2004 (Fig. 2B). The snail density in low-elevation bottomlands increased from 2003 to 2005 and then began to decrease sharply in 2006. However, during the 2nd decade of TGD (2013–2022), snail density experienced a rebound from 0.098/0.11 m2 in 2013 to 0.125/0.11 m2 in 2022. Notably, over half (six of nine) of the bottomlands witnessed an increase in snail density between 2013 and 2022, namely Matian, Xingang, and Junshan Park in Dongting Lake, Ganyu in Poyang Lake, and Xin and Lao in the Anhui section. The snail density in the Anhui section of the Yangtze River, which is far away from the dam, increased significantly from 0.059/0.1 m2 to 0.256/0.1 m2 (Fig. 2C, D).
The trends in snail density from 1997 to 2022: A average snail density downstream from bottomlands of the dam, B trends in snail density at high-, medium-, and low-elevation bottomlands, C trends in snail density in the Dongting Lake area, Poyang Lake area, and Anhui section, D trends in snail density at high-, medium-, and low-elevation bottomlands in the Anhui section
Changes in the water regime, vegetation, and microclimate
The average FD in the nine downstream bottomlands of the TGD was 73.56 ± 44.31 days for the range of years covered in this study, with the bottomlands experiencing a minimum of 5 days of flooding and a maximum of 161 days of flooding. The maximum value of annual highest water level was 41.80 m (occurring before the TGD at Shiguishan in Dongting Lake), and the minimum value of the annual lowest water level was 1.40 m (occurring after the TGD at Xinqiaozhao in Anhui section). The NDVI and NL range was from 0.31 to 0.83 and 69.31 to 256.83 for nin3 bottomlands, with the lowest values of NDVI and NL occurring before the operation of the dam. The range of RH and SSD was 76.25–84.38% and 1730–2058 h (Additional file 1: Table S4).
Following the 1st decade of the TGD between 2003 and 2012, there was a noticeable decline in both water levels and their fluctuations, as depicted in Fig. 3. However, during the 2nd decade of the TGD (2013–2022), there was a slight resurgence in water levels and their fluctuations. During 2013–2022, the annual average water levels in the Dongting Lake region, Poyang Lake region, and Anhui section of the Yangtze River rose by 0.25 m, 0.32 m, and 0.53 m, respectively, with the rise becoming more pronounced far away from the TGD. A more detailed analysis reveals that during the periods 1997–2002, 2003–2012, and 2013–2022, the average water levels in the Dongting Lake area, Poyang Lake area, and Anhui section exhibited variations of 28.65–28.24–28.49, 12.00–11.11–11.43, and 8.26–7.43–7.96, respectively. Similarly, the annual fluctuations in water levels during these periods in the Dongting Lake area, Poyang Lake area, and Anhui section showed changes of 11.41–9.97–10.09, 11.39–10.62–10.92, and 9.18–7.91–8.51, respectively.
The FD on the bottomlands decreased from 122 days during 1997–2002 to 57 days during 2003–2012 (Fig. 4A) and then rose again to 65.5 days in 2013–2022 (P < 0.001). In the Dongting Lake and Poyang Lake regions, flooding duration fluctuated during 2003–2022 (Fig. 4B, C). However, in the Anhui Section far from the TGD, the FD increased steadily between 2008 and 2022, with durations of 38, 53, and 88 days during 2008–2012, 2013–2017, and 2018–2022, showing an upward trend (Fig. 4D).
As shown in Fig. 5, the bottomland environmental factors experienced various changing trends from 1997 to 2022. The NDVI increased from 0.53 between 1997 and 2002 to 0.64 between 2003 and 2012 and kept steady between 2013 and 2022. Similarly, the NL rose from 34.77 between 1997 and 2002 to 61.82 between 2003 and 2012 and reached 92.89 between 2013 and 2022 (χ2 = 30.02, P < 0.05). The Tmax showed a rising trend, with values increasing from 32.80 ℃ between 1997 and 2002 to 33.95 ℃ between 2003 and 2012, slightly decreasing to 33.70 ℃ between 2013 and 2022 (χ2 = 19.13, P < 0.05). The GST exhibited an obvious increasing trend with values of 19.24 ℃, 19.49 ℃, and 19.68 ℃ from 1997 to 2002, 2003 to 2012, and 2013 to 2022 (χ2 = 19.10, P < 0.05), respectively. However, the differences in the Tmin, precipitation, and SSH between 1997–2002, 2003–2012, and 2013–2022 were not statistically significant.
Changes in micro-environmental factors of the snail-breeding bottomland during the periods of 1997–2002, 2003–2012, and 2013–2022, including NDVI (A), NL (B), Tmin (C), Tmax (D), GST (E), Pre (F), and SSH (G). NDVI, normalized difference vegetation index; NL, night light index; Tmin, average minimum temperature in January; Tmax, average maximum temperature in July; GST, ground surface temperature; Pre, precipitation; SSH, sunshine hours
Impact of environmental change on snails
The Tmin, NL (a high night light value reflects the high development and utilization of the bottomland), FD, and NDVI were the four most important factors affecting snail density (Fig. 6). Based on the LightGBM model, a SHAP diagram was conducted to assess the relationship between individual factors and snail density (Fig. 7). At an FD of 50–120 days and NDVI of 0.45 and 0.70, the SHAP values were generally positive, indicating conditions favorable for snail reproduction. Conversely, at a Tmin < 0.3 ℃ and NL > 40, the SHAP values were generally negative, indicating unfavorable conditions for snail breeding. Furthermore, conditions such as higher rainfall and ground surface temperature, and a sunshine duration of 1700 to 1800 h were conducive to snail breeding.
Beeswarm (A) and bar (B) plot showing the contribution of a single factor affecting snail density. Tmin the average minimum temperature in January; NL night light index, FD flooding duration, NDVI normalized difference vegetation index, Tmax average maximum temperature in July, GST ground surface temperature, SSH sunshine hour, Pre precipitation
SHAP plot illustrates the relationships between snail density and the factors FD (A), NDVI (B), NL (C), Tmin (D), Pre (E), GST (F), SSH (G), and Tmax (H). FD flooding duration, NDVI normalized difference vegetation index, NL night light index, Tmin the average minimum temperature in January, Pre precipitation, GST ground surface temperature, SSH sunshine hour, Tmax average maximum temperature in July
Next, we fitted a total of 24 models including different variables and different nonlinear approximated terms for the environmental factors (Table 1). Model 1 to model 4 included the results of the univariate models and model 5 to model 24 included the results of the multivariate models. Model 20 was selected as the final reported model because it had the lowest DIC value.
The fixed effects results of model 20 are listed in Fig. 8 and Table 2. After adjusting for related confounding factors and the space-time effect (Fig. 8), the RR of increased snail density rose with a flooding duration between 20 to 100 days. However, it decreased as FD increased over 100 days. The relationship between night lighting and snail density was “L” shaped. At low night lighting levels, the RR of increased snail density was higher. As the NDVI increased from 0 to 0.6, the RR of increased snail density also rose. However, the RR then decreased as the NDVI continued to increase. When the Tmin ranged from − 1 to 2 °C, the RR of increased snail density grew as it increased. An increasing trend was detected among bottomlands with a greater GST compared with bottomlands in the first quartile cluster; bottomlands in the second, third, fourth quartile clusters of the GST had higher RR values for snail density, which were 1.271 (95%CI 1.082–1.493), 1.302 (95%CI 1.146–1.480), and 1.278 (1.048, 1.559), respectively.
RR of increased snail density associated with FD (A), NL (B), NDVI (C), and Tmin (D) in the posterior distribution after adjusting for related factors. RR relative risk, FD flooding duration, NL night light index, NDVI normalized difference vegetation index, Tmin average minimum temperature in January
Discussion
Compared to previous studies [15, 16] conducted during the first decade of the TGD construction, this research revealed that during the second decade of the dam’s operation, various water level indicators (average water level and annual water level fluctuation range) showed an increasing trend (Fig. 3). These changes can be attributed to two primary factors: the increase in extreme rainfall and the varying impacts of the dam's operational phases on downstream water and sediment dynamics. On the one hand, the intensity and frequency of heavy rainstorms in the Yangtze River basin increased under the influence of global climate change [38, 39], contributing to rising water levels. On the other hand, some studies have shown that in the early stages of the TGD, discharge of clean water caused strong scouring of the downstream river channel and promoted decreasing water levels [40, 41]. However, as time went on, the sediment content in the river increased and sediment siltation in the river may have occurred, resulting in a rising water level and establishing a new water-sediment balance [40, 41]. The new water-sediment balance and the change pattern in precipitation might explain why the water level rebounded in 2013–2022 compared with 2003–2012. The recovery of water levels during this period showed regional differences, with a significant rise observed in the Anhui section of the Yangtze River (Fig. 3). Spatially, areas closer to the dam experienced stronger riverbed scouring, while regions farther downstream transitioned more easily to a sediment deposition phase and a rise in water levels [40, 41]. In addition, Guo et al. [42] showed that the most significant impact of the TGD was concentrated on the river section near the TGD, and the degree of impact was five times that of downstream section such as the Anhui section. Therefore, the water levels in the Anhui section are influenced by local tributaries and precipitation more than Dongting Lake. These reasons might be why water level recovery is more obvious in these regions.
The Bayesian spatial-temporal model showed that after adjusting for confounding factors (temperature, vegetation, rainfall, etc.), the RR of increased snail density rose with flooding duration between of 20 and 100 days, while it decreased when FD increased to > 100 days. The continuous change in water level as a fundamental element affected the abundance and distribution of snails. After the construction of the TGD, the water level gradually decreased, and the duration of flooding on the bottomlands was much shorter, affecting the water environment required for normal hatching of snail eggs and snail development [43, 44]. After 2003, the decline in water levels first impacted high-elevation bottomlands, rendering them unsuitable for snail breeding. This was followed by the medium-elevation bottomlands, aligning with the observed decrease in snail density (which began in 2003 for high-elevation areas and in 2004 for medium-elevation areas). Interestingly, snail density in low-elevation bottomlands increased from 2003 to 2005 but then sharply declined starting in 2006. This pattern may be attributed to the prolonged submersion under water of low-elevation bottomlands (e.g. above 100 days), which exceeded the optimal conditions for snail breeding. Initially, the reduction in water levels and shorter flooding duration promoted an increase in snail density. However, as water levels continued to drop beyond the snails' tolerance, their density began to decrease. In addition, the night light index, which reflects the influence of human activities, infrastructure, and other factors on the environment [45], steadily increased in areas surrounding the bottomlands following the construction of the TGD. The study revealed that snail density showed an “L” shape with NL. The RR of increased snail density was lower at higher NL. Since the TGD filled with water, the lower water levels extended the time in which bottomlands remained unflooded. These changes facilitated activities such as farming, cultivation, and tourism and spurred infrastructure development in surrounding areas [46]. Utilization of bottomland resources also contributed to the decline in snail density.
Hydrological changes are not the sole factor influencing Oncomelania snail density; micro-environmental conditions, including vegetation and surface temperature, also contribute to the dynamic fluctuations in snail populations. This study revealed that snail density showed a slight rebound, and over half (six/nine) of the bottomlands witnessed a fluctuating increment in snail density between 2013 and 2022. Zhou et al. [47] observed the effects of summer flooding on snail reproduction and found that the number of snail offspring in the 3-month flooding group was 2.5 times that of the non-flooding group. The recovery of the water level facilitated the extension of the flooding duration of the snail habitat, improving the condition of snail breeding and promoting juvenile snails to develop into adult snails [44, 46]. In addition, the ground surface temperature of the bottomland exhibited a rising trend from 1997 to 2022, and an increasing trend of snail density was detected among bottomlands with a greater GST. Laboratory research showed that the development rate of snail eggs could be accelerated with increasing temperature [48, 49]. Previous studies [48, 49] also have shown that extremely high temperatures may limit normal activities of snails such as opening their shells and foraging. However, suitable vegetation could counteract the limiting effect of high temperatures on snails [21, 50]. The NDVI of the bottomland has a small increase within the suitable range for snail survival, which can regulate the extremely high temperature of snail habitats and protect snails from strong sunlight. Therefore, improving hydrological conditions of snail habitats in conjunction with the changing microenvironment contributed to the recovery of snail density. There is also concern in other regions (such as the Philippines) that environmental changes related to dam construction may affect the prevalence of schistosomiasis by affecting freshwater snails [51, 52]. Hence, constructing a real-time monitoring and response system (involving key environmental factors and snails) is crucial for controlling potential risks of schistosomiasis transmission.
This study has some limitations. First, our data came from multiple sources such as the meteorological bureau, hydrological stations, and remote sensing inversion, and there may be measurement bias. Second, studies on the changes in microenvironment, such as vegetation, temperature, etc., need to consider the impact of climate change. In the future, advanced monitoring equipment (such as temperature and humidity monitors, soil sensors, artificial video, etc.) can be installed in snail-breeding bottomlands to precisely collect data covering multiple dimensions of environmental changes in real time. Using big data and artificial intelligence technology, the collected data can be deeply analyzed to predict the reproduction and spreading trends of snails. When the monitoring system detects that environmental changes may cause a surge of snails, an emergency plan should be immediately activated and corresponding control measures should be taken, such as environmental modification, chemical drug release, or biological control. This can not only improve early warning capabilities but also provide a scientific basis for formulating prevention and control measures.
Conclusions
After construction of the dam, the reduction of the water level and flooding duration of the bottomlands and the increase in the utilization of the bottomland were conducive to reducing snail density. However, over time, the inhibitory effect of the TGD on the snail may be weakening, which is mainly manifested in the rising and unstable water level trend, especially in areas far from the TGD. Environmental changes in the snail-breeding bottomland may bring uncertainty to the snail population. In recent years, the rebound of the snail density in local bottomlands may be related to the rise in water levels and the change in the microenvironment. Establishing an efficient monitoring and response system is crucial for the precise control of snails and the elimination of schistosomiasis.
Availability of data and materials
Supporting data for the conclusions of this article are included within the article. The raw data supporting the conclusions of this article will be made available upon reasonable request.
Abbreviations
- TGD:
-
The Three Gorges Dam
- FD:
-
Flooding duration
- NDVI:
-
Normalized difference vegetation index
- NL:
-
Nighttime light index
- SSH:
-
Sunshine hour
- GST:
-
Ground surface temperature
- Pre:
-
Precipitation
- Tem:
-
Temperature
- RH:
-
Relative humidity
- Tmin:
-
Average minimum temperature in January
- Tmax:
-
Average maximum temperature in July
- RR:
-
Relative risk
- AI:
-
Artificial intelligence
- SHAP:
-
Shapley additive explanations
- LightGBM:
-
Light gradient boosting machine
- DIC:
-
Deviance information criterion
References
Li Q, Zheng JX, Jia TW, Feng XY, Lv C, Zhang LJ, et al. Optimized strategy for schistosomiasis elimination: results from marginal benefit modeling. Parasit Vectors. 2023;16:419.
Gong YF, Feng JX, Luo ZW, Xue JB, Guo ZY, Zhang LJ, et al. Spatiotemporal heterogeneity of schistosomiasis in mainland China: evidence from a multi-stage continuous downscaling sentinel monitoring. Asian Pac J Trop Med. 2022;15:26–34.
Gong YF, Zhu LQ, Li YL, Zhang LJ, Xue JB, Xia S, et al. Identification of the high-risk area for schistosomiasis transmission in China based on information value and machine learning: a newly data-driven modeling attempt. Infect Dis Poverty. 2021;10:88.
WHO guideline on control and elimination of human schistosomiasis. Geneva: World Health Organization; 2022. https://www.who.int/publications/i/item/9789240041608#:~:text=This%20guideline%20provides%20evidence-based%20recommendations%20in%20the%20following,of%20WASH%20and%20snail%20control%20activities%20to%20su. Accessed 8 Nov 2024.
Gong Y, Tong Y, Jiang H, Xu N, Yin J, Wang J, et al. Three Gorges Dam: the changing trend of snail density in the Yangtze River basin between 1990 and 2019. Infect Dis Poverty. 2023;12:45.
McCreesh N, Booth M. Challenges in predicting the effects of climate change on Schistosoma mansoni and Schistosoma haematobium transmission potential. Trends Parasitol. 2013;29:548–55.
Lv SB, He TT, Hu F, Li YF, Yuan M, Xie JZ, et al. The impact of flooding on snail spread: the case of endemic schistosomiasis areas in Jiangxi Province, China. Trop Med Infect Dis. 2023;8:259.
Zhou XN, Yang GJ, Yang K, Wang XH, Hong QB, Sun LP, et al. Potential impact of climate change on schistosomiasis transmission in China. Am J Trop Med Hyg. 2008;78:188–94.
Guo SY, Li L, Zhang LJ, Li YL, Li SZ, Xu J. From the One Health perspective: Schistosomiasis Japonica and flooding. Pathogens. 2021;10:1538.
Lanza GR, Upatham S, Chen A. A Place-Based Conceptual Model (PBCM) of Neotricula aperta/Schistosoma mekongi habitat before and after dam construction in the Lower Mekong River. PLoS Negl Trop Dis. 2023;17:e0011122.
Boon NAM, Mbow M, Paredis L, Moris P, Sy I, Maes T, et al. No barrier breakdown between human and cattle schistosome species in the Senegal River basin in the face of hybridisation. Int J Parasitol. 2019;49:1039–48.
Sokolow SH, Jones IJ, Jocque M, La D, Cords O, Knight A, et al. Nearly 400 million people are at higher risk of schistosomiasis because dams block the migration of snail-eating river prawns. Philos Trans R Soc Lond B Biol Sci. 2017;372:20160127.
Zhou YB, Liang S, Chen Y, Jiang QW. The Three Gorges Dam: does it accelerate or delay the progress towards eliminating transmission of schistosomiasis in China? Infect Dis Poverty. 2016;5:63.
Zeng L, Panizzo VN, Wang Z, Huang X, Chen X, McGowan S. Hydrological isolation accelerates algal blooms in floodplain lakes: Biomarker evidence from Dongting Lake, China and its satellite lake. Water Res. 2025;279:123430.
Li F, Ma S, Li Y, Tan H, Hou X, Ren G, et al. Impact of the Three Gorges project on ecological environment changes and snail distribution in Dongting Lake area. PLoS Negl Trop Dis. 2017;11:e0005661.
Wu JY, Zhou YB, Chen Y, Liang S, Li LH, Zheng SB, et al. Three Gorges Dam: impact of water level changes on the density of schistosome-transmitting snail Oncomelania hupensis in Dongting Lake Area, China. PLoS Negl Trop Dis. 2015;9:e0003882.
Zhou Y. Application of artificial intelligence in parasitic diseases and parasitology. Zhongguo Xue Xi Chong Bing Fang Zhi Za Zhi. 2024;36:551–4.
Yang L, Driscol J, Sarigai S, Wu Q, Lippitt CD, Morgan M. Towards synoptic water monitoring systems: a review of AI methods for automating water body detection and water quality monitoring using remote sensing. Sensors. 2022;22:2416.
Qureshi YM, Voloshin V, Facchinelli L, McCall PJ, Chervova O, Towers CE, et al. Finding a husband: using explainable AI to define male mosquito flight differences. Biology. 2023;12:496.
Villela DA, Codeço CT, Figueiredo F, Garcia GA, Maciel-de-Freitas R, Struchiner CJ. A Bayesian Hierarchical model for estimation of abundance and spatial density of Aedes aegypti. PLoS ONE. 2015;10:e0123794.
Gong Y, Tong Y, Jiang H, Xu N, Yin J, Wang J, et al. Three Gorges Dam: Potential differential drivers and trend in the spatio-temporal evolution of the change in snail density based on a Bayesian spatial-temporal model and 5-year longitudinal study. Parasit Vectors. 2023;16:232.
Chen S, Li YL, Duan L, Liu JB, Zhou J, Lin DD, et al. Assessment of the influence of levees along Yangtze River on Oncomelania hupensis, the intermediate host of Schistosoma japonicum. Parasit Vectors. 2024;17:291.
Zhou YB, Liang S, Chen GX, Rea C, He ZG, Zhang ZJ, et al. An integrated strategy for transmission control of Schistosoma japonicum in a marshland area of China: findings from a five-year longitudinal survey and mathematical modeling. Am J Trop Med Hyg. 2011;85:83–8.
Li YL, Dang H, Guo SY, Zhang LJ, Feng Y, Ding SJ, et al. Molecular evidence on the presence of Schistosoma japonicum infection in snails along the Yangtze River, 2015–2019. Infect Dis Poverty. 2022;11:70.
Yang Y, Zheng SB, Yang Y, Cheng WT, Pan X, Dai QQ, et al. The Three Gorges Dam: does the flooding time determine the distribution of schistosome-transmitting snails in the Middle and Lower Reaches of the Yangtze River, China? Int J Environ Res Public Health. 2018;15:1304.
Shi Y, Jin N, Ma X, Wu B, He Q, Yue C, et al. Attribution of climate and human activities to vegetation change in China using machine learning techniques. Agric For Meteorol. 2020;294:108146.
Li M, Cao S, Zhu Z, Wang Z, Myneni RB, Piao S. Spatiotemporally consistent global dataset of the GIMMS Normalized Difference Vegetation Index (PKU GIMMS NDVI) from 1982 to 2022. Earth Syst Sci Darta. 2023;5:4181–203.
Zhang X, Zhong Z, Zhang M, Zhao F, Wu Y, Sun Y, et al. Analysis of anthropogenic disturbance and spatial and temporal changes of bird communities in plateau wetlands fusing bird survey and nighttime light remote sensing data. J Environ Manage. 2025;375:124349.
Wu Y, Shi K, Chen Z, Liu S, Chang Z. Developing improved time-series DMSP-OLS-like data (1992–2019) in China by integrating DMSP-OLS and SNPP-VIIRS. IEEE Trans Geosci Remote Sens. 2021;60:4407714.
Kuo PF, Huang TE, Putra IGB. Comparing Kriging estimators using weather station data and local greenhouse sensors. Sensors. 2021;21:1853.
Chen YY, Cai SX, Xiao Y, Shan XW, Zhang J, Liu JB. Impact of implementation of Three Gorges Project on schistosomiasis endemic situation in Hubei Province. Zhongguo Xue Xi Chong Bing Fang Zhi Za Zhi. 2014;26:498–503.
Sedgwick P. Multiple hypothesis testing and Bonferroni’s correction. BMJ. 2014;349:g6284.
Yan J, Xu Y, Cheng Q, Jiang S, Wang Q, Xiao Y, et al. LightGBM: accelerated genomically designed crop breeding through ensemble learning. Genome Biol. 2021;22:271.
Wang K, Tian J, Zheng C, Yang H, Ren J, Liu Y, et al. Interpretable prediction of 3-year all-cause mortality in patients with heart failure caused by coronary heart disease based on machine learning and SHAP. Comput Biol Med. 2021;137:104813.
Valiente Fernández M, de Lesmes González Aledo A, Delgado Moya FP, Badía Martín I. SHAP model explainability in ECMO-PAL mortality prediction: a critical analysis. Intensive Care Med. 2023;49:1559.
Liao J, Qin Z, Zuo Z, Yu S, Zhang J. Spatial-temporal mapping of hand foot and mouth disease and the long-term effects associated with climate and socio-economic variables in Sichuan Province, China from 2009 to 2013. Sci Total Environ. 2016;563–564:152–9.
Blangiardo M, Cameletti M, Baio G, Rue H. Spatial and spatio-temporal models with R-INLA. Spat Spatiotemporal Epidemiol. 2013;7:39–55.
Wu X, Wang L, Niu Z, Jiang W, Cao Q. More extreme precipitation over the Yangtze River Basin, China: insights from historical and projected perspectives. Atmos Res. 2023;292:106883.
Bai L, Rong Y. Reanalysis of the characteristics of extreme rainfall in the Yangtze River basin during recent 50 years. Water Resour Res. 2015;4:88–100.
Zheng J, Gu XG, Xu YL, Ge JH, Yang XX, He CH, et al. Relationship between the transmission of schistosomiasis japonica and the construction of the Three Gorge Reservoir. Acta Trop. 2002;82:147–56.
Zhu HM, Xiang S, Yang K, Wu XH, Zhou XN. Three Gorges Dam and its impact on the potential transmission of schistosomiasis in regions along the Yangtze River. EcoHealth. 2008;5:137–48.
Guo S, Xiong L, Zha X, Zeng X, Cheng L. Impacts of the Three Gorges Dam on the streamflow fluctuations in the downstream region. J Hydrol. 2021;598:126480.
Zheng Y, Shi Y, Li S, Deng W, Zhou J, Hu B, et al. Environmental determinants for snail density in Dongting Lake region: an ecological study incorporating spatial regression. Am J Trop Med Hyg. 2022;107:1178–84.
Yang Y, Gao J, Cheng W, Pan X, Yang Y, Chen Y, et al. Three Gorges Dam: polynomial regression modeling of water level and the density of schistosome-transmitting snails Oncomelania hupensis. Parasit Vectors. 2018;11:183.
Han G, Zhou T, Sun Y, Zhu S. The relationship between night-time light and socioeconomic factors in China and India. PLoS ONE. 2022;17:e0262503.
Zhang GM, Ma SC, Cao ZG, Wang Y, Zhang SQ, Wang TP, et al. Water level changes caused from the Three Gorges Dam impact on the breeding and distribution of Oncomelania snails in the marshland of Chenqiaozhou, Hexian County, Anhui Province. Chin J Dis Control Prevent. 2018;22:943–7.
Zhou XN, Huang JZ, Zong ZM, Wu Zheng, Kan JG, Yang HM, et al. Comparative study on the effects of flooding on snail proliferation during summer flood season. Zhongguo Xue Xi Chong Bing Fang Zhi Za Zhi. 1989; (02):23–25.
Kalinda C, Chimbari M, Mukaratirwa S. Implications of changing temperatures on the growth, fecundity and survival of intermediate host snails of schistosomiasis: a systematic review. Int J Environ Res Public Health. 2017;14:80.
Manyangadze T, Chimbari MJ, Gebreslasie M, Ceccato P, Mukaratirwa S. Modelling the spatial and seasonal distribution of suitable habitats of schistosomiasis intermediate host snails using Maxent in Ndumo area, KwaZulu-Natal Province. South Africa Parasit Vectors. 2016;9:572.
Liu MM, Feng Y, Yang K. Impact of micro-environmental factors on survival, reproduction and distribution of Oncomelania hupensis snails. Infect Dis Poverty. 2021;10:47.
Diakité NR, Winkler MS, Coulibaly JT, Guindo-Coulibaly N, Utzinger J, N’Goran EK. Dynamics of freshwater snails and Schistosoma infection prevalence in schoolchildren during the construction and operation of a multipurpose dam in central Côte d’Ivoire. Infect Dis Poverty. 2017;6:93.
Semenza JC, Rocklöv J, Ebi KL. Climate change and cascading risks from infectious disease. Infect Dis Ther. 2022;11:1371–90.
Acknowledgements
We acknowledge all the teachers, reviewers, and editorial teams who provided assistance in the writing, revision, and editing of this article.
Funding
This work is supported by the National Natural Science Foundation of China (nos. 82273754, 81673236). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Author information
Authors and Affiliations
Contributions
YB Zhou, SZ Li, and QW Jiang conceived and designed the study. YB Zhou, SQ Zhang, DD Lin, and Y Cai collected and managed the data. SQ Zhang, SB Lv, M Zheng, and BJ Hu conducted the snail survey. YF Gong, XL lei, JH Huang, Y Zhou, LY Zhu, JM Wang, and N Xu performed statistical analyses. YF G and YB Z drafted and revised the manuscript. YB Zhou, Y Chen, and SZ Li revised the manuscript and provided important intellectual content. All authors have approved the final manuscript for publication.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Ethical review and approval were not required for the study in accordance with the local legislation and institutional requirements.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
13071_2025_6782_MOESM1_ESM.docx
Additional file 1: Table S1. Information on the nine bottomlands. Table S2. Correlation analysis of variables. Table S3. Results of Shapiro-Wilk test. Table S4. Characteristics of variables in the bottomlands in the middle and low reaches of the Yangtze River basin.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Gong, Y., Zhang, S., Lin, D. et al. Decoding the impact of environmental shifts on snail density dynamics in the Yangtze River basin: a 26-year study. Parasites Vectors 18, 156 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13071-025-06782-3
Received:
Accepted:
Published:
DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13071-025-06782-3