Three Gorges Dam: Impact of Water Level Changes on the Density of Schistosome-Transmitting Snail Oncomelania hupensis in Dongting Lake Area, China.

BACKGROUND: Schistosomiasis remains an important public health issue in China and worldwide. Oncomelania hupensis is the unique intermediate host of schistosoma japonicum, and its change influences the distribution of S. japonica. The Three Gorges Dam (TGD) has substantially changed the ecology and environment in the Dongting Lake region. This study investigated the impact of water level and elevation on the survival and habitat of the snails.
METHODS: Data were collected for 16 bottomlands around 4 hydrological stations, which included water, density of living snails (form the Anxiang Station for Schistosomiasis Control) and elevation (from Google Earth). Based on the elevation, sixteen bottomlands were divided into 3 groups. ARIMA models were built to predict the density of living snails in different elevation areas.
RESULTS: Before closure of TGD, 7 out of 9 years had a water level beyond the warning level at least once at Anxiang hydrological station, compared with only 3 out of 10 years after closure of TGD. There were two severe droughts that happened in 2006 and 2011, with much fewer number of flooding per year compared with other study years. Overall, there was a correlation between water level changing and density of living snails variation in all the elevations areas. The density of living snails in all elevations areas was decreasing after the TGD was built. The relationship between number of flooding per year and the density of living snails was more pronounced in the medium and high elevation areas; the density of living snails kept decreasing from 2003 to 2014. In low elevation area however, the density of living snails decreased after 2003 first and turned to increase after 2011. Our ARIMA prediction models indicated that the snails would not disappear in the Dongting Lake region in the next 7 years. In the low elevation area, the density of living snails would increase slightly, and then stabilize after the year 2017. In the medium elevation region, the change of the density of living snails would be more obvious and would increase till the year 2020. In the high elevation area, the density of living snails would remain stable after the year 2015.
CONCLUSION: The TGD influenced water levels and reduced the risk of flooding and the density of living snails in the study region. Based on our prediction models, the density of living snails in all elevations tends to be stabilized. Control of S. japonica would continue to be an important task in the study area in the coming decade.

Introduction

Schistosomiasis remains a serious public health problem worldwide, affecting more than 200 million people in approximately 76 countries with a loss of 1.53 million disability-adjusted life years (DALYs) [1]. Schistosomiasis japonica (S. japonica) is distributed in 12 provinces in China and 11.6 million people have been infected since 1949 [2-4]. S. japonica causes the most hazardous schistosomiasis, and is difficult to prevent and treat [5, 6]. After continued implementations of comprehensive control measures from the mid-1950s to 1980s, endemic regions were circumscribed in certain core areas in China, especially in the Dongting Lake region at the middle reaches of the Yangtze River [7-11]. Most of these core areas are in the downstream of Three Gorges Dam (TGD). The TGD is a world-class water conservancy project. It began to impound water and sediment discharge in 2003, and is one of several tremendous engineering projects transforming China’s ecology and natural environment. The construction and the operation of TGD have obviously affected the ecological environment. The water level of TGD in 2003 was 135m and it reached 160m in 2011. It is believed that the TGD project reduces the frequency of major flooding in the downstream areas from once every ten years to once every 100 years; however it threatens the living of aquatic animals in the Yangtze River including the river dolphin, or baiji, and finless porpoise, or jiangzhu [12]. Similarly the TGD project can also influence the survival of O. hupensis, the unique intermediate host of Schistosoma japonicum [13, 14].

The TGD has been completed for over 10 years, and its impact on the transmission of S. japonica has to be evaluated [15, 16]. This study aimed to determine the impact of changes in water level and number of flooding per year on the snail (O. hupensis) density, and to predict the changes in the density of living snails in the Dongting Lake region using the autoregressive integrated moving average (ARIMA) model, which combines the advantages of autoregressive (AR) model and moving average (MA) model.

Methods

1. Study area and sampling

This study was conducted in the 16 bottomlands near four hydrological stations in Anxiang County of the Dongting Lake region (Fig 1). These bottomlands are outside the embankment. The Dongting Lake is located at 28°30′–30°20′N and 111°40′–113°40′E in the northeastern part of Hunan Province and covers a water surface area of 2,681 km2, and plays an important role in regulating the amount of water in the Yangtze River [17]. It is a typical S. japonica endemic area with an at-risk population of about 429,000 people in marshland and lake regions [18]. Anxiang County is located in the Dongting Lake area and at the downstream of the TGD. There are 4 hydrological stations located in Anxiang County and for each hydrological station we selected 4 bottomlands as our study fields (Fig 1): 1) Guanyuan hydrological station—Qinglong, Dongbao, Nanyang and Yangshutan; 2) Zizhiju hydrological station—Yandoukou, Dongdi, Yongheyuan, Yifenju; 3) Shiguishan hydrological station—Qiangkou, Huangjiatai, Fuxing, Wuyang; and 4) Anxiang hydrological station—Yucheng, Liujiao, Zhulin, Anxiang. A total of sixteen bottomlands were divided into 3 groups according to three elevation: Low (<33m), medium (33 to 35m) and high (>35 m).

Fig 1

Map of sampling area and hydrological station in Anxiang county, Hunan province, People’s Republic of China.

2. Data sources

The routine data of snails which could be obtained and used was snail density and snail area. The change in snail areas was not large in this study field and period. The density of living snails is an important index reflecting the survival and reproduction of snails. Thus the density of living snails was used to assess the impact of TGD on snails. The Anxiang Station for Schistosomiasis Control provided data of density of living snails in every bottomland of Anxiang County for the period from 2004 to 2014. Snail surveys were conducted in each spring and were implemented using a traditional method of random quadrant sampling (0.11 m2-sized frames, 20 m apart between frames) [19].

We downloaded data of water levels from the Hunan Flood Prevention Information System, including daily water levels at 8:00 am at the 4 hydrological stations from 1995 to 2013. The elevation for the 4 hydrological stations was 33m for the Guanyuan station, 31m for the Zizhiju and Anxiang stations, and 30m for the Shiguishan station. Data of elevation were collected through the Google Earth (Google Ltd, USA, http://www.google.com/earth/) by sampling 50 points (5*10, 200 m apart between points) for each bottomland (Fig 2). A weighted average elevation was calculated for each bottomland.

Fig 2

Sketch of collecting points in a bottomland in Anxiang county, Hunan province, People’s Republic of China.

3. Analysis

We calculated number of flooding per year, which is calculated in the basis of the difference between the daily water level at each hydrological station and the weighted average elevation of the bottomland near the station. If the difference was greater than 0, the bottomland around the hydrological station was considered to be flooded.

We used the density of living snails time-series data to fit an optimum ARIMA model and predict the density of living snails in the future. ARIMA model is a component of the Box-Jenkins approach to time-series modeling [20-22]. An ARIMA model is derived by combining the three techniques: autoregressive modeling (AR), moving average modeling (MA) and differencing. They are presented as ARIMA(p,d,q), in which p and q represent the orders of AR and MA models, respectively, and d denotes the order of differencing. In ARIMA models, we assume a stationary time series, which means that the data vary around a constant mean and variance over time [Chin J Parasit Dis Con. 2004 ">23-25]. Nonstationary time series variables can be converted into stationary ones. An ARIMA(p,d,q) model can be written as [26, 27]: where c is a constant, y’t = yt-yt-1 represents the differenced series, y’t-p are lagged value and zt is a white noise process.

The ARIMA modeling procedure consists of three iterative steps: selecting a candidate model, estimating the model and performing diagnostic tests and forecasting.

Before fitting an ARIMA model, selecting a candidate model is the process of identifying randomness, stationarity and seasonality using the autocorrelation functions (ACF) and partial autocorrelation functions (PACF). If data do not meet these requirements, a transformation of data should be implemented [28]. ACF is a statistical tool that measures whether an earlier value in the series has some relation to a later value. PACF captures the amount of correlation between a variable and a lag of this variable that is not explained by correlation at all low-order lags. Parameters in the ARIMA model(s) are estimated with the conditional least squares (CLS) method [29]. Since there are three parameters in an ARIMA model, different parameter combinations will lead to various results. The second step is to fit an optimum ARIMA model based on the Bayesian information criterion (BIC), for which the less BIC is, the better the model fits the data. Finally, the fitted model can be used to forecast the density of living snails and its confidential interval [30]. IBM SPSS 20.0 (IBM Corporation USA, http://www-01.ibm.com/software/analytics/spss/) was used for all the analysis.

Results

1. Univariate analyses

Fig 3 shows the water levels at the Anxiang hydrological station from 1995 to 2013, with a warning water level for flooding of 37.0 m. June 1st 2003 is the day that TGD began to impound. Before the year of 2003, 7 out of 9 years had a water level beyond the warning level at least once at this station compared with only 3 out of 10 years after the year of 2003 (2003 not included. Fig 3 also shows similar decreasing trends for the occurrence frequency of flooding and water level after 2003.

Fig 3

Plot of water level of Anxiang hydrological station.

Fig 4 shows the number of flooding per year and density of living snails in low, medium and high elevation areas. The data of 2009 were not included due to its too many outliers. The density of living snails was correlated with the water level a year before, which is caused by a lag effect of the water level. There were two severe droughts happened in 2006 and 2011.

Fig 4

Number of flooding per year and density of living snails in different elevation areas.

The density of living snails in the high and medium elevation areas had a continuously decreasing trend from 2004 to 2014. In the low elevation area, the density of living snails had a similar trend between 2004 and 2012 but turned to increase after 2012. In the medium elevation areas, the density of living snails decreased sharply from 2004 to 2008 and stabilized after 2008.

Table 1 presents number of flooding per year before and after 2003 in 3 different elevations. After the construction of TGD, the number of flooding significantly decreased in the areas of both medium elevation (t = 4.519, p<0.000) and high elevation areas (t = 3.475, p = 0.001), but increased slightly in the low elevation area (t = -1.428, p = 0.158).

Table 1

Number of flooding per year before and after 2003 in different elevation areas.

	Low elevation		Medium elevation		High elevation
	Before 2003	After 2003	Before 2003	After 2003	Before 2003	After 2003
Minimum	88	43	1	2	3	3
Maximum	166	199	118	130	91	54
Mean	135.9	144.4	67.56	47.58	50.83	31.88
Standard deviation	19.01	32.62	28.9589	27.96	24.84	16.35

2. Time-series analysis

ARIMA model was fitted to predict the density of living snails in the future in the 3 different elevation areas in the basis of historical time-series data. Autocorrelation graphs of the density of living snails indicated that the data were randomized and smooth. The parameter group (1,0,0) was preferred according to the BIC values of the density of living snails in the low, medium and high elevation areas (-4.12, -0.61 and -2.71).

Fig 5 shows the results of ARIMA (1,0,0) model in the low elevation area, which generates predicted values based on historical time-series data. The predicted value of the density of living snails in the low elevation area would increase from 0.37 /0.11 m2 in 2015 to 0.40/0.11 m2 in 2019 and after.

Fig 5

Observed and predicted density of living snails in the low elevation area.

Similarly, Fig 6 shows the results of ARIMA (1,0,0) model for the medium elevation area, and the predicted value of the density of living snails would increase from 0.42/0.11m2 in 2015 to 0.71/0.11m2 in 2020. Fig 7 presents the findings of ARIMA (1,0,0) model for the high elevation area, and the predicted values would not change markedly between 2015 and 2020.

Fig 6

Observed and predicted density of living snails in the medium elevation area.

Fig 7

Observed and predicted density of living snails in the high elevation area.

Discussion

After the impoundment of TGD in 2003, the dry season in Dongting Lake region was reported to arrive earlier and to be longer than before, and the water level was close to the lowest level in history for several times [31-34]. Anxiang County of Hunan Province is located in Dongting Lake region and at the downstream of TGD. Our results indicated that the water levels in all the elevation areas and number of flooding decreased after 2003 and that there was a severe drought in 2006.

Overall, there was a correlation between water level changing and density of living snails variation in all the elevations areas. The density of living snails in all elevations areas was decreasing after the TGD was built. The relationship between number of flooding per year and the density of living snails was more complicated. In the medium and high elevation areas, the density of living snails kept decreasing from 2003 to 2014. In low elevation area, however, the density of living snails decreased after 2003 first and turned to increase after 2011. The data of 2014 show that, the value of the density of living snails high in the low elevation area, low in the high elevation area and in between in the medium elevation area. Number of flooding per year decreased after 2003, so was the density of snails in the medium and high elevation areas. The association did not reach statistical significance in the low elevation area. Oncomelania is an amphibious snail, and its larva needs to live in water. When it grows into adult stage, it tends to inhabit in a humid region, like grass. Water is one of the necessary conditions for growth and reproduction of the snails that have to live in water or a wet place, and it is difficult to survive in the dry environment [35-37]. In the medium and high elevation areas, the reduced number of flooding per year might result in droughts in some months, which led to the decreased density of living snails. It might not be the case in the low elevation area. It is reported that the density of living snails in middle reaches of the Yangtze River including Dongting Lake Region is sharply reduced from 2003 to 2014 [38]. Another research finds that the density of living snails in Poyang Lake Region is declined after the impoundment of TGD [39]. The results in this study are consistent with the previous findings.

Droughts in this area are associated with the density of living snails. The severe drought in 2006 caused a decline of the density of living snails in 2007. The frequency of droughts in the medium and high elevation areas after 2003 was more than that before 2003, and as a consequence the density of living snails had a marked change. In the low elevation area however, the frequency of droughts after 2003 was less than that before 2003 and the decline of the density of living snails in this region was less obvious as compared with other areas.

Based on our prediction models, the snails would not disappear in the Dongting Lake region in near future. In the low elevation area, the density of living snails would increase slightly and then stabilize after the year of 2017. In the medium elevation region, the change of the density of living snails would be more obvious and would increase till the year of 2020. In the high elevation area, the density of living snails would remain stable after the year of 2015. Control of the snails would continue to be a challenge in the study area in the coming decade, and some measures, such as snail surveillance and fencing bovines on the marshland, need to continue to be implemented.

There are some limitations in this study. Firstly, data of the density of living snails were not available before 2003 and we were not able to make a comparison of the density of living snails before and after TGD was built. Secondly, only water level and elevation were studied while other factors such as climate conditions and plant types were not considered. However, we did not identify any important factors that had changed dramatically during the study period.

In conclusion, hydrology is an important determinant for the density of living snails. The density of living snails changed in different elevation areas, which was correlated with the variation of hydrology. TGD influenced the hydrology in the study areas and further changed the density of living snails. Based on the results from ARIMA models, we predicted that controlling snails would continue to be a challenge in the study area in the coming decade although TGD might lead to a reduction in the density of living snails in the region.