Long-Term Changes of Aquatic Invasive Plants and Implications for Future Distribution: A Case Study Using a Tank Cascade System in Sri Lanka

: Climate variability can inﬂuence the dynamics of aquatic invasive alien plants (AIAPs) that exert tremendous pressure on aquatic systems, leading to loss of biodiversity, agricultural wealth, and ecosystem services. However, the magnitude of these impacts remains poorly known. The current study aims to analyse the long-term changes in the spatio-temporal distribution of AIAPs under the inﬂuence of climate variability in a heavily infested tank cascade system (TCS) in Sri Lanka. The changes in coverage of various features in the TCS were analysed using the supervised maximum likelihood classiﬁcation of ten Landsat images over a 27-year period, from 1992 to 2019 using ENVI remote sensing software. The non-parametric Mann–Kendall trend test and Sen’s slope estimate were used to analyse the trend of annual rainfall and temperature. We observed a positive trend of temperature that was statistically signiﬁcant ( p value < 0.05) and a positive trend of rainfall that was not statistically signiﬁcant ( p values > 0.05) over the time period. Our results showed ﬂuctuations in the distribution of AIAPs in the short term; however, the coverage of AIAPs showed an increasing trend in the study area over the longer term. Thus, this study suggests that the AIAPs are likely to increase under climate variability in the study area.


Introduction
Freshwater ecosystems are characterised by high levels of biodiversity, and they also provide vital ecosystem services to sustain life on earth [1]. These important ecosystems account for merely 0.8% of the land surface on earth, but a vast number of species (100,000 species or 6% of total species) that live on the planet depend on this vital ecosystem [2]. Moreover, these aquatic systems contain the highest number of species threatened with extinction along with more imperiled habitats; thus, they are considered as severely endangered ecosystems on earth [3,4]. Negative influences on this ecosystem have been attributed to a variety of causes i.e., land use changes, nutrient and chemical pollution, climate change, and invasive plants [5]. The ecological impacts caused by aquatic invasive alien plants (AIAPs) on freshwater habitats is widely acknowledged. It is one of the key drivers of biodiversity depletion in fresh water systems. Aquatic invasive alien plants cause a significant impact on the richness of native plants and animals across many aquatic habitats [6][7][8]. They change structure and function of aquatic systems and rapidly change habitats [1]. Ecosystem services are disrupted resulting in direct economic losses [9]. Certain AIAPs (i.e., Salvinia molesta and Eichhornia crassipes) proliferate rapidly and produce a mat on the water surface, which eventually impedes light and oxygen supply to the lower level plants and animals, leading to destruction of the entire ecosystem [10,11]. In and its relation to the distribution of AIAPs. The specific aims of this study were to (i) evaluate the long-term spatial and temporal spread dynamics of AIAPs, (ii) analyse the longterm trends of rainfall and temperature, (iii) examine how IAPS distribution may change in response to climate change and investigate potential impacts on agriculture. To the best of our knowledge, no systematic study has been undertaken to understand the influence of climate variations on the range dynamics of AIAPs in Sri Lanka. Identifying changes in surface cover of AIAPs is vital to fully understand how they respond to climate change and further to provide implications and impetus towards evidence-based decision making.

Study Area
The present study was conducted in the Nachchaduwa tank and three cascade systems, located in the Anuradhapura district encompassing five divisional secretariats of Thirappane, Nachchaduwa, Nuwaraam Palatha East, Ipalogama, and Kekirawa on latitude 8 • 6 -8 • 18 N and 80 • 27 -80 • 36 E, Sri Lanka ( Figure 1). A tank cascade system consists of small man-made reservoirs that are connected by canals and spillways and drain into a large water reservoir [39,40]. The Nachchaduwa tank cascade system (here after TCS) comprises of 67 tanks of different sizes with a total area of approximately 26 km 2 (blue colour area in Figure 1) that spreads over an area located along a generally flat, undulating to rolling terrain. This area receives runoff water mainly from the northeast monsoon rains in the months of September to February. This unique agricultural system provides an outstanding service by conserving biodiversity and ecosystem services in the landscape.
AIAPs, (ii) analyse the long-term trends of rainfall and temperature, (iii) examine how IAPS distribution may change in response to climate change and investigate potential impacts on agriculture. To the best of our knowledge, no systematic study has been undertaken to understand the influence of climate variations on the range dynamics of AIAPs in Sri Lanka. Identifying changes in surface cover of AIAPs is vital to fully understand how they respond to climate change and further to provide implications and impetus towards evidence-based decision making.

Study Area
The present study was conducted in the Nachchaduwa tank and three cascade systems, located in the Anuradhapura district encompassing five divisional secretariats of Thirappane, Nachchaduwa, Nuwaraam Palatha East, Ipalogama, and Kekirawa on latitude 8°6′-8°18′ N and 80°27′-80°36′ E, Sri Lanka ( Figure 1). A tank cascade system consists of small man-made reservoirs that are connected by canals and spillways and drain into a large water reservoir [39,40]. The Nachchaduwa tank cascade system (here after TCS) comprises of 67 tanks of different sizes with a total area of approximately 26 km 2 (blue colour area in Figure 1) that spreads over an area located along a generally flat, undulating to rolling terrain. This area receives runoff water mainly from the northeast monsoon rains in the months of September to February. This unique agricultural system provides an outstanding service by conserving biodiversity and ecosystem services in the landscape.
The TCS, which has been sustained for centuries, is presently confronted with challenges from aquatic plant invasion, which is particularly linked to impacts of climate change. At present, a majority of the cascade tanks are covered with some of the world's worst invasive plants: S. molesta and E. crassipes. Salvinia molesta and E. crassipes were introduced to Sri Lanka as early as 1930 and 1905 respectively and soon dispersed to the water bodies across the country [11,41]. By the early 1980s, S. molesta and E. crassipes had infested the tanks system in the dry zone, where the majority of the tanks were entirely covered with a thick growth [11,42].  The TCS, which has been sustained for centuries, is presently confronted with challenges from aquatic plant invasion, which is particularly linked to impacts of climate change. At present, a majority of the cascade tanks are covered with some of the world's worst invasive plants: S. molesta and E. crassipes. Salvinia molesta and E. crassipes were introduced to Sri Lanka as early as 1930 and 1905 respectively and soon dispersed to the water bodies across the country [11,41]. By the early 1980s, S. molesta and E. crassipes had infested the tanks system in the dry zone, where the majority of the tanks were entirely covered with a thick growth [11,42].

Landsat Data
Landsat 4-5, 7, and 8 remote sensing data were employed in this study. Multi-spectral Landsat images were downloaded in GeoTIFF format for the years 1992, 1996, 2000, 2003, 2007, 2010, 2013, 2016, 2017, and 2019 from USGS Earth explorer Website (https: //earthexplorer.usgs.gov/). Thus, the gap between data acquisitions was approximately similar to make a comparative analysis of the changes. The images were of good quality with relatively few clouds over the area of interest. Clouds and cloud shadows influence image classifications, and it was challenging to find cloud-free Landsat images, which is a frequent problem in tropical countries [33,43,44]. Thus, accessibility to good quality Landsat images was an issue in this study region, which restricted the number of images that can be used in the analysis. The selected Landsat images were only good quality images during the first quarter (JFM) that were not affected (or minimally affected) by clouds and cloud shadows during the study period 1992-2019. The details of the imageries used in this study are presented in Table 1. The northeast monsoon brings rain to the dry zone from September to February, and water bodies in this area become full by the first few months of the year. Thus, images acquired during the first quarter of the year (wet season) were selected to avoid possible influences attributed to seasonal changes (i.e., water level fluctuations) of AIAPs distribution over the period of time and to examine long-term changes. We received the updated geospatial boundary demarcation vector data (shapefile) of TCS from the Department of Agriculture Sri Lanka, which was used to mask the study area.

Image Processing
Image processing, which involves the conversion of digital number values of pixels to useful information, is undertaken in three consecutive steps: pre-processing, image enhancement, and classification [45]. We used ENVI version 5.5 (Exelis Visual Information Solutions, Boulder, CO, USA) for satellite image processing and further analyses. ArcMap version 10.4.1 was employed for producing maps. Image pre-processing is vital to improving image quality and performance, which is done by suppressing unnecessary distortions while enhancing some image features [46]. Digital image processing techniques, radiometric calibration, and dark object subtraction were used to enhance Landsat images before classification. Radiometric correction was undertaken as the study compared multiple images across a long period of time [47]. This technique corrects for errors in digital number values of image pixels and improves the interpretability and quality of remotely sensed data [33,48]. Furthermore, we used dark-object subtraction to remove atmospheric scattering effects before analysis for efficient extraction of land cover information [49]. It is a simple but widely used and accepted atmospheric correction approach [34].
Image classification assigns pixels of the image into groups of similar land use classes/clusters on the earth's surface [45,50]. Landsat satellite remote sensing data contain different spectral variations for different feature types [51]. Thus, features in the study area were classified using the supervised classification technique. This technique assigns specific pixel value ranges for land use classes [45]. Therefore, for each of the pre-processed images, local knowledge of the study area and visual interpretation by referring to historical Google Earth images were used to identify the major land cover types that were common in all the time series images. Through this process, the four features identified and delineated as classification classes were water, AIAPs, non-aquatic plants, and open areas ( Table 2). Aquatic plants found on the water surface of the TCS generally represent E. crassipes and S. molesta; however, other native and non-native aquatic plants could be find occasionally. For each of the training class, regions of interest (ROIs) were defined using a specified colour code. Visualisation of the various features during the creation of ROIs was enhanced through false colour composite by variation of the Red-Green-Blue (RGB) bands. ROIs were randomly selected as much as possible for better capturing pixel variability. For image classification, maximum likelihood supervised classification was used. This method assigns pixels to the nearest classification class (highest probability of belonging to a class) while running the classification.

Accuracy Assessment
In maximum likelihood classification, some pixels can be classified incorrectly due to an inconsistent distribution of data [52]. In order to provide reliable results, image classification studies should produce a high degree of overall accuracy, which provides a quantitative measure of how effectively pixels are categorised into respective feature classes [53]. Thus, accuracy assessment of the classified images plays an important role in order to assess the reliability of remote sensing-derived information contained in the maps and to provide reliable input data for subsequent scientific investigations and policy decisions [43,54]. In the study, ground truth ROIs were randomly acquired for accuracy assessment by making reference to historical Google Earth imagery of the same year, as field-based data collection is not possible for accuracy assessment of historical images. The number of ground reference pixels that were developed by drawing polygons around the features of interest varied temporarily for different land cover classes, and this was based on ease of identification and variabilities within the classes. Therefore, the assessment was carried out using 117, 163, 144, 173, 206, 155, 229, 280, 125, and 187 pixels for the 1992,1996,2000,2003,2007,2010,2013,2016,2017, and 2019 images, respectively. Statistical analysis of the classified images and the ground reference pixels for the various years was done using confusion matrix (error matrix) in ENVI. Accuracy of the classification was assessed in terms of overall accuracy, kappa coefficient, producer's accuracy, and user's accuracy. The kappa statistic, that ranges from 0 (no agreement) to 1 (complete agreement), is highly recommended for assessing classification accuracy as it provides interclass discrimination specifically than overall accuracy [55]. Using ArcMap software, the areas covered by each land cover type were extracted from the classified images.

Long-Term Trend Analysis of Climate Variables
To understand and describe the variations in climate parameters in the study area, we analysed temporal climate data of the Anuradhapura district over the period 1990-2019, obtained from the Department of Meteorology, Sri Lanka. These climate data are from one weather station (Anuradhapura) that covers the TCS and surrounding areas. We analysed time-series data of climate parameters, annual average mean temperature, and total annual rainfall ranging from 1990 to 2019 to understand the variability and trends that could influence AIAPs distribution in the long term. We tested the statistical significance of the trend (i.e., α = 0.10, 0.05) using the widely used non-parametric Mann-Kendall trend analysis statistical test [56,57] and Sen's slope estimator [58] using XLSTAT software (Addinsoft 2020; https://www.xlstat.com). The Mann-Kendall test determines whether there is a significant positive or negative trend of the variables over time. The Sen's slope estimator assesses the magnitude and direction of the detected trend in the sample of N pairs of data, which results in a positive value (>0) indicating an increasing trend and a negative value (<0) indicating a decreasing trend [59,60]. The relationship between the climate variability and AIAPs distribution in the study area over the period was examined.

Accuracy Assessment and Land Use
The results of the accuracy assessment of the classified images for the ten-time periods are displayed in Table 3. We obtained good accuracy estimates for all ten images, with >87% for overall accuracy and >0.8 for the Kappa coefficient for each of the images. According to Lea and Curtis [61], these accuracy levels are satisfactory for the calculation of land use/cover areas. Thus, our classification results are robust with good performance.   The distribution of land uses in the study area over the period 1992-2019 is presented in Figure 2, whereas the temporal changes in coverage of various land uses during the study period are presented in Table 4. According to the figures, the coverages of AIAPs vary during the study period. The range dynamics of AIAPs is illustrated graphically in Figure 3, which shows a positive trend during the study period. The trends were statistically significant with p values of 0.01 at 95% confidence level. Additional details relevant to area changes are provided in Table 5. Accordingly, we observed an overall increase of 28.6% of the coverage of AIAPs over the period from 1992 to 2019.

Long-Term Trend Analysis of Climate Variables
The results of climate data analysis showing the magnitude of temporal variability of rainfall and temperature are presented in Figure 4. The annual average temperature and total annual rainfall show upward or increasing trends during the study period (Figure 4a). Figure 4b shows that the highest amount of average monthly rainfall is recorded in the months of October to December, and the lowest is in the months of June to September. The average monthly temperature ranges from 26.1 to 29.6 °C. During the first quarter of the study years, the total rainfall showed an increasing trend, whereas average temperature shows a decreasing trend (Figure 4c).
The summary results on long-term trend analysis based on Mann-Kendall and Sen's slope statistical tests for total annual rainfall and annual average temperature for the period of 1990-2019 in the TCS are presented in Table 6. At a significance level of 0.05, the p-value shows that the null hypothesis is rejected for annual average temperature, suggesting a significant trend in the time series; however, there is no trend for total rainfall at this significance level. We observed a significant trend for rainfall at a lower significance level of 0.1. According to Mahmood et al. [62], trends at α = 0.05 signify a strong signal of trend whereas trends at α = 0.1 signify a weak signal of trend. Thus, the analysis of climate data over the past three decades indicates a strong increasing signal for average temperature and weak increasing signal for total rainfall in the study area. We did not consider the seasonal changes of climate variables, as all images downloaded belonged to the same season of the year over the study period.

Long-Term Trend Analysis of Climate Variables
The results of climate data analysis showing the magnitude of temporal variability of rainfall and temperature are presented in Figure 4. The annual average temperature and total annual rainfall show upward or increasing trends during the study period (Figure 4a). Figure 4b shows that the highest amount of average monthly rainfall is recorded in the months of October to December, and the lowest is in the months of June to September. The average monthly temperature ranges from 26.1 to 29.6 • C. During the first quarter of the study years, the total rainfall showed an increasing trend, whereas average temperature shows a decreasing trend (Figure 4c).
The summary results on long-term trend analysis based on Mann-Kendall and Sen's slope statistical tests for total annual rainfall and annual average temperature for the period of 1990-2019 in the TCS are presented in Table 6. At a significance level of 0.05, the p-value shows that the null hypothesis is rejected for annual average temperature, suggesting a significant trend in the time series; however, there is no trend for total rainfall at this significance level. We observed a significant trend for rainfall at a lower significance level of 0.1. According to Mahmood et al. [62], trends at α = 0.05 signify a strong signal of trend whereas trends at α = 0.1 signify a weak signal of trend. Thus, the analysis of climate data over the past three decades indicates a strong increasing signal for average temperature and weak increasing signal for total rainfall in the study area. We did not consider the seasonal changes of climate variables, as all images downloaded belonged to the same season of the year over the study period.

Discussion
There have been very few published studies in the literature relevant to the range dynamics of AIAPs in freshwater aquatic systems. Although it is challenging, we explored the potential of medium-resolution satellite remote sensing data for assessing the longer-

Discussion
There have been very few published studies in the literature relevant to the range dynamics of AIAPs in freshwater aquatic systems. Although it is challenging, we explored the potential of medium-resolution satellite remote sensing data for assessing the longer-term changes of AIAPs spread in the freshwater aquatic system. In this study, we considered all species of AIAPs as one class, and individual species were not considered due to insufficient resolution of Landsat data. Thus, we did not undertake the accuracy assessment for individual species. As such, the high level of overall accuracy and Kappa coefficient we received for the images are explained. According to the findings, all fourland use classes (i.e., non-aquatic plants, AIAPs, open areas, water) display variations in area coverage. An analysis of the spread of AIAPs in the TCS show variable trend in coverage over the period 1992 to 2019. The changes fluctuate over short-term periods but increase over the 27-year period. We observed declines in the distribution of AIAPs in 1996 (0.18% per year), in 2016 (0.14% per year), and in 2019 (0.06% per year) relevant to the previous assessment year. However, the overall distribution of AIAPs from 1992 to 2019 showed a positive response with the increasing temporal variability of total annual rainfall and annual average temperature. The spread of AIAPS in the TCS has become a challenging issue for ecosystem functions in the study area and thus, the land managers should consider the negative consequences, particularly to agriculture. Distribution of land uses over the past two decades also signifies a general decline in water cover and increased distribution of non-aquatic plants in the study area. The siltation and sedimentation of tanks due to soil erosion because of intensive agriculture in the area is a major impediment for sustenance of the dry zone TCS [38]. Consequently, various native and non-native plants acquire these silted areas in the tank. Illegal cultivations on the tank bunds is another issue when water level fluctuates during the year. The distribution of open land area seems to be quite varying over the time period due to water level fluctuations.
Global warming is manifested by rising temperatures in air, water, and the earth's surface; and increased frequency of temperature extremes (i.e., hot days, heat waves) [63]. The rising temperature has implications on the distribution of non-native species under climate change [64,65]. The projected effects can make changes in aquatic systems by altering their pathways of introduction, establishment, spread, and distribution [19,64]. Under changing climate, ranges of non-native species are likely to expand into higher altitudes and latitudes. Experimental studies have confirmed that the CO 2 enrichment and growth enhancement of AIAPs is strongly influenced by warmer temperatures [66,67]. Elevated temperature intensifies biomass production of aquatic macrophytes and influences the rapid development of shoots [68]. Thus, water temperature can be a crucial driver that influences the establishment, distribution, and impact of AIAPs [64]. The results of climate data analysis using regression analysis for the study period, Mann-Kendall test, and Sen's slope estimator showed strong significant positive trends for annual average temperature and a weaker significant positive trend for total annual rainfall in the study area at a significance level of 0.05. The temperature increase in the study area can be due to various reasons, such as global warming, urbanisation, or atmospheric circulation, and such changes can lead to increased evapotranspiration and unusually dry spells [60].
The increasing spread of AIAPs that was observed in this study area can be attributed to several factors. Nitrogen deposition increases the abundance of exotics in an ecosystem while depleting its species richness [69]. Due to the increased usage of nitrogen fertiliser in agriculture, the excessive runoffs are added to the water bodies, resulting in increased growth of AIAPs [17]. Nitrogen-rich sediment deposits in the TCS might have resulted in the steady increase in non-aquatic plants cover in the study area. In TCS, the excess water of a tank upstream is diverted to the tank downstream through the paddy fields. Although this process enables an efficient use of water resources, it may ultimately cause increased nutrient content deposition in water bodies, resulting in a luxuriant growth of AIAPs [70]. The study revealed that small tanks are heavily infested by AIAPs, and the major tank is relatively free from invasions ( Figure 2). This could be due to the continuous maintenance and cleaning of the major tank. The statistics given in Table 4 for AIAPs is mainly contributed by the infestations occurring in the small tanks. The distribution of AIAPs in the small tank system is significant relevant to the size of the tanks and thus, the effective management of AIAPs should be an important concern in the sustenance of the TCS. AIAPs coverage in the water tanks increased by around 28.6% over the 1992-2019 period.
In its 5th assessment report, the IPCC announced that the projected surface temperature change of the earth is expected to exceed 1.5 • C by 2100 relative to 1850-1900, and more intense and more frequent extreme precipitation events are predicted [71]. Climate plays a vital role in the geographic distribution of plant species [72]. Climate change makes alterations in the freshwater habitats, and such affected habitats have the potential to be highly vulnerable and adversely impacted by biological invaders [7,73,74]. We have observed positive relationships between climate variables and range dynamics of AIAPs in the study area during this period. Experimental studies have illustrated that the spatiotemporal distribution of plants is affected by variations in temperature and rainfall [75,76]. Through the bioclimatic modeling of three invasive plants, Egeria densa, Myriophyllum aquaticum, and Ludwigia spp., Gillard, et al. [77] highlighted that long-term changes occurring in the climate may have a substantial impact on the potential distribution of AIAPs in freshwater systems.
The results of this study imply a long-term decreasing trend of the area of surface water in the TCS. In addition, the findings suggest that the AIAPs have already expanded their coverage over the last few decades. We are not aware of the relative contribution of each variable to the AIAPs spread; however, in view of the relation of climate variability with the area of spread of AIAPs in the TCS, this expansion will continue in the future. Therefore, managing the potential negative consequences of AIAPs would be an increasing challenge to the land managers. The vegetative growth rate of AIAPs, such as Water hyacinth and Salvinia, is substantially high, and they constantly decay and produce a considerable amount of debris (biomass), which can cause the conversion of AIAPs areas to open areas over time [17]. The change in the land use pattern of the surrounding areas (for example, conversion from forest areas to agriculture) may also contribute to a high level of surface runoff, leading to increased soil erosion and silt deposition in the water bodies; the result of this is a reduction in the water-holding capacity of the tank system. AIAPs (such as water hyacinth) cause high levels of evapotranspiration, leading to significant water losses in water bodies [78]. In addition, they block irrigation canals and interfere with water supply to agriculture. Thus, the expansion of AIAPS would have severe negative consequences for agricultural production, as TCS serve as the main source of water for irrigated agriculture in this area.
Therefore, the continuous monitoring and evaluation of potential impacts of climate changes on freshwater aquatic systems is of great importance. Remote sensing image analysis, based on the supervised maximum likelihood classification method, provides important information for land managers, such as spatiotemporal spread patterns, range changes, and potential distributions of AIAPs. Such information is reliable and can be included in the formulation of policies and management actions for better land use management. In the TCS, propagules can disperse from one tank to the other easily, especially during rainy periods, as tanks are interconnected. Thus, this should be considered while formulating necessary management actions. The findings of the present study may be important for conservation planners and water resource managers to be vigilant on AIAPs issue as a potential impact of climate variability and climate change. The results of this study also provide implications for the use of Landsat data for similar ecological studies in the region.

Limitations and Challenges of the Study
The information provided by the moderate resolution Landsat data is not sufficient to fully differentiate AIAPs to the species level. Aquatic invaders may share the same habitat niche with other non-AIAPs; thus, perhaps our results may overestimate the spread of AIAPs. However, the primary occupants in the tank system are AIAPs (personal observations by 1st and 4th authors), and non-AIAPs have always been a minor component. Thus, we presume that the interference by non-AIAPs would be minimal for the analysis. Further, we have taken all possible steps to correct and enhance the quality and improve the appearance of Landsat images. Practically, classification is sometimes challenging for some pixels, as the classification algorithm may not identify land uses correctly due to certain reasons (i.e., when water bodies are turbid or depending on the level of chlorophyll-a contents in water), and such situations may interfere with the classifications [51]. Heavy rainfall can inundate shallow open areas, and thus, such areas can be classified as water. While classifying the images, there was a challenge with the weeds (non-aquatic or nonnative plants) growing on the silted areas because of the similar spectral reflectance with AIAPs. We minimised this error by using the Google Earth image of the same year, as the AIAPs could be identified through their distinct growing pattern and coloration, which was unique and different from other plants. Although Google Earth is a useful image analysis tool, such historical images are also not sharp enough, particularly those before 2011.
In this study, we did not intend to look at seasonal variations in the spread of the AIAPs; our key aim was to look at long-term changes. That is the reason why we chose all the images of the same season. Hence, even if there were seasonal variations, by keeping all the images to one season, we would still be able to pick the long-term trends. Seasonal changes and effects are removed by selecting images of the same time period. This is an appropriate and accepted methodology for looking at long-term changes. Our results of post-classification represent only one land use in one pixel, and it may overlook small sub-pixel level conversions [79]. However, in reality, some pixels represent several land uses. The pixel size we used in these Landsat images (30 m × 30 m) is not at fine enough resolution to overcome these issues. As such, a thorough knowledge of the study area is needed to minimise such issues. Furthermore, freshwater ecosystems are invaded by multiple plant invaders [80]. However, AIAPs observed in this study area (S. molesta, E. crassipes) generally form monocultures, and thus, the possibility to assign erroneous classification is relatively low. In addition to climate variability, several other factors, such as nutrient level, waves and water currents, and human influences, may contribute to the AIAPs distribution at various scales [81,82]. However, in this study, we did not consider other potential causes, as climate variability is considered as the most influential factor for plant invasions [72].

Conclusions
Freshwater aquatic systems are vulnerable to biological invaders under changing climate conditions. We conducted this research to examine how AIAPs' distribution responds to changes in climate variability in a selected TCS in Sri Lanka using ten historical time periods from 1992 to 2019. The findings revealed that the AIAPs infestation is increasing in the water systems over the long term, though there are shorter-term fluctuations. Based on the results of the Mann-Kendall test, Sen's slope estimator, and linear regression, we demonstrated positive trends in the annual average temperature and total rainfall in the study area. The output of the study shows a potential relationship between AIAPs and the changes in temperature, which could be highly important in understanding the future spatial and temporal distribution of AIAPs under climate change.