Assessing the Indicators of Water Stress in Usangu Agricultural Wetland, Tanzania ()
1. Introduction
The relationship between Land Surface Temperature (LST), Normalized Difference Vegetation Index (NDVI), Standardized Precipitation Index (SPI), and precipitation is critical for understanding climate change (Gaikwad et al., 2022; Khan et al., 2018). Landsat and Sentinel imagery have been extensively used for climate analysis and hydrological studies over several decades, providing valuable insights into water resource management (Cheng et al., 2024; Hoai et al., 2025). The Usangu wetland, a significant water catchment within the Ruaha River hydrological system (Mwita, 2016), plays a vital role in regional water supply. However, research indicates a notable annual decline in the Ruaha River water flow, which adversely affects downstream socio-ecological systems (MOW, 2015). The Usangu Wetland is not a Ramsar site under the Convention on Wetlands (Kashaigili et al., 2006). Its conservation relies entirely on Tanzania’s national policies and stewardship commitment, lacking international safeguards.
Remote Sensing (RS) techniques are widely recognized for their ability to evaluate water stress resulting from climatic variations, offering a robust framework for monitoring environmental changes. NDVI trend analysis is frequently employed to assess vegetation responses to climate change (Muhury et al., 2024), with studies demonstrating a strong linear correlation between NDVI and climate change impacts (Zhang et al., 2024). LST serves as a key indicator of the energy balance between land surfaces and water bodies (Song et al., 2025), and its variations are closely tied to topographic factors, land cover, and vegetation indices. Literature suggests that NDVI fluctuations are influenced by both climate change and anthropogenic activities (Tuoku et al., 2024). LST significantly impacts NDVI variations due to energy imbalances (Karnieli et al., 2010; Xu et al., 2011). Latitude and altitude variations further affect vegetation growth, with a negative LST-NDVI relationship associated with altitude often indicating soil water stress and increased evapotranspiration rates (Nemani & Running, 1997; Stisen et al., 2007).
Monitoring water presence in wetlands like the Usangu wetland often involves indices such as Normalised Difference Water Index (NDWI), Modified Normalised Difference Water Index (MNDWI), Modified Normalised Difference Water Index 2 (MNDWI2), and the Water in Wetland (WIW) method (Acharya et al., 2018). For instance, a study on the Renuka wetland, a Ramsar site, utilized the Normalized Difference Water Index (NDWI) through Google Earth Engine (GEE) to map water bodies using Landsat imagery. NDWI effectively identifies water distribution, even in areas with minimal vegetation. The WIW method (Lefebvre et al., 2019) provides a comprehensive RS approach for detecting water in vegetated regions, leveraging Near Infrared (NIR) and Short-Wave Infrared 2 (SWIR2) bands, making it adaptable across diverse vegetation types.
2. Materials and Methods
Study Area
The Usangu wetland, located in Mbarali District, Mbeya Region, Tanzania, lies between latitudes 7˚45'S to 9˚00'S and longitudes 33˚40'E to 35˚00'E and consist of five subsystems (Franks, 2003). The Usangu wetland is the study area within Mbarali District, though its hydrological influence extends slightly beyond, including the riparian reach through the Ruaha National Park and the Mtera/Kidatu hydroelectric system. The Usangu sub-catchment drainage area spans approximately 20,800 km2, with 4840 km2 (23% of the plains) situated below 1100 m above mean sea level (msl), while the remaining 77%, representing the higher catchment, ranges from 1100 m to 2000 m above msl (Hyandye & Martz, 2017).
The Usangu Plains as shown in Figure 1 and Figure 2 consist of western and eastern wetlands connected by the Great Ruaha River, with other major perennial rivers draining the plains, including the Mbarali, Kimani, Chimala, and Ndembera (Kashaigili et al., 2006). Smaller rivers, such as the Umrobo, Mkoji, Lunwa, Mromboji, Ipatagwa, Mambi, Kioga, Mjenje, Kimbi, Itumbo, and Mswiswi, also contribute to the drainage network. The plains feature extensive alluvial deposits that support irrigation, livestock keeping, and dryland farming, as illustrated in Figure 3. In the Usangu wetland, mean rainfall ranges from 500mm to 700mm, and mean annual temperature ranges from 18˚C to 28˚C (Thomas et al., 2024). The wetland’s outflow sustains the Ruaha National Park, providing essential ecosystem services for wildlife conservation, serves as a water source for the Kidatu Dam, and contributes approximately 15% to the 4 million capacity Stigler Gorge Hydro Power Dam (Franks et al., 2004).
3. Methodology
Wetland areas condition was assessed using 30m SRTM DEM, SPI (PERSIANN-CDR/CHIRPS), LST (Landsat/Sentinel-2), WIW, NDVI (MODIS), NDWI/MNDWI. Validation of the CHIRPS data done using historical CMIP6. The datasets used in analysis and their processing are listed in Table 1.
3.1. SPI Data
The precipitation dataset sourced from the Climate Data Record (CDR) of NOAA PERSIANN, which is accessible through the cloud-based Google Earth Engine (GEE) platform using JavaScript API code. This dataset has a spatial resolution of 0.25 degrees. The dataset was clipped and filtered to match the wetland area. The Standardized Precipitation Index (SPI) was calculated to assess temporal and spatial patterns of drought and wet conditions, with SPI values ranging from -3 (extreme drought) to 2 (extreme wet). For visualization, a timeseries plot was generated to examine the temporal consistency of SPI across the study period (as shown in ). Spatial consistency was validated by comparing SPI with precipitation clustering, derived using the K-Means algorithm, and presented in .
Figure 1. Usangu catchment map.
Figure 2. The wetland area.
Figure 3. Cropland in the wetland area.
3.2. LST Data
Land Surface Temperature (LST) modelling for the Usangu wetland was conducted using thermal infrared bands ST_B10 from Landsat 8 and 9. The Landsat derived LST was then statistically compared to Sentinel-2 spectral bands (B2-B12) and topographic features elevation, slope and NASA SRTM Digital Elevation Model (DEM) data to predict high-resolution (10m) LST across the wetland, processed within Google Earth Engine (GEE). To ensure data quality, cloud filtering techniques were applied to mask clouds and generate cloud-free images, minimizing noise in the dataset for training and testing. The images were then processed using a median function to create composite images and cropped to fit the study area.
A Random Forest model was employed for LST prediction, trained using machine learning techniques with input features including spectral bands, vegetation indices, and elevation data derived from the SRTM DEM. The model was configured to use 70% of the dataset for training and 30% for testing. Model performance was evaluated by comparing the predicted LST to reference LST data, with errors quantified using the Mean Absolute Error (MAE). The relationship between reference and predicted LST was assessed through a linear correlation, visualized with a trend line and quantified using the Pearson correlation coefficient.
3.3. Precipitation Data
Precipitation data for the study area were sourced from the Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS) dataset (Zhang et al., 2025), accessed as a spatial dataset via the cloud based Google Earth Engine (GEE) platform. The CHIRPS dataset provides a high spatial resolution of 0.05 degrees and consists of daily precipitation records, covering the analysis period since 1981 (Agha Kouchak et al., 2011). To examine trends, seasonal patterns, and residual components in the precipitation data, statistical time series decomposition was performed using the Seasonal-Trend Decomposition using LOESS (STL) method (He et al., 2022; RB, 1990). The STL approach was selected for its ability to elucidate natural variability in precipitation data, aiding in the interpretation of trends, seasonal cycles, and time series characteristics, as well as identifying anomalies (He et al., 2022).
Coupled Model Intercomparison Project Phase 6 (CMIP6) historical precipitation data was bias corrected using CHIRPS observations via linear scaling (Lafon et al., 2013; Teutschbein & Seibert, 2012). The resulting time series Figure 4 reveals severe dry spells (≤200 mm/month) in the Usangu wetland region during 2010-2015.
Figure 4. Bias corrected CMIP6 precipitation.
K-Means clustering of precipitation data was conducted using machine learning techniques within the GEE Python environment, as illustrated in Figure 10 (Reddy et al., 2023) and equation 1. This process utilized geographical data from 299 satellite images obtained from the Global Precipitation Measurement (GPM) dataset, specifically NASA/GPM_L3/IMERG_MONTHLY_V07, which has a spatial resolution of 0.1 degrees. The GPM data, available monthly, spanned the period from 2000 to 2024. The K-Means algorithm was applied to categorize the precipitation data into four classes: Heavy, Moderate, High, and Low. A timeseries plot was generated to compare the classified precipitation data with the Standardized Precipitation Index (SPI) for validation.
(1)
Ci: Data points in cluster, μi: Mean (centroid) of cluster, x: Precipitation value.
3.4. NDVI Data
The Normalized Difference Vegetation Index (NDVI) data were sourced from the MODIS MOD13A1.061 product (Terra Vegetation Indices 16-Day Global 500 m), covering the period from 2019 to 2023 using equation 2. The NDVI data were clipped to the wetland area. To improve data quality, the NDVI infrared bands were filtered to remove cloud cover and then multiplied by the appropriate scale factor provided in the MODIS dataset to enhance visualization. Monthly averages of NDVI were calculated for each year from 2019 to 2023 to capture temporal variations in vegetation health across the study period Figure 13.
, (2)
3.5. Visualization of MDWI, ANDWI, NDVI and WIW
Visualization of water conditions in the Usangu wetland was conducted using spectral indices, including the Normalized Difference Vegetation Index (NDVI), Modified Normalized Difference Water Index (MDWI), Normalized Difference Water Index (NDWI), Water in Wetland (WIW), and Automated Normalized Difference Water Index (ANDWI). The data were sourced from Landsat 8 Level 2 (surface reflectance) images, covering the period from 2020 to 2024. These images were filtered to include only those with less than 10% cloud cover to ensure data quality for the wetland area. A composite median image was derived from the surface reflectance bands of the filtered satellite images to reduce temporal variability and enhance consistency.
The NDVI map was generated in Figure 15(a) using the near-infrared (NIR) and red bands to assess vegetation health, calculated as NDVI = (NIR − Red)/(NIR + Red). To detect water amidst the wetland complex vegetation cover, multiple water indices were employed, leveraging various spectral bands. NDWI, as described in prior studies (Ashok et al., 2021; Laonamsai et al., 2023; Xu, 2006), was calculated using the NIR and green bands to identify water bodies. The MDWI (Laonamsai et al., 2023) utilized the short-wave infrared (SWIR) band alongside the green band to enhance water detection in vegetated areas. The WIW index method (Lefebvre et al., 2019) was applied using surface reflectance bands in the NIR and SWIR ranges, incorporating a logical thresholding approach to improve water detection reliability, particularly when other indices were less effective. The ANDWI (Rad et al., 2021) was computed using red, green, blue (RGB), NIR, and both SWIR1 and SWIR2 bands to improve water demarcation and visibility in vegetated wetlands.
(3)
(4)
(5)
(6)
(6a)
Table 1. Description of dataset.
Dataset |
Description |
Year |
Source |
DEM |
Digital Elevation Model data used to analyse elevation impacts on LST and other variables, with a spatial resolution of 30 m. |
2024 |
NASA SRTM (Shuttle Radar Topography Mission) DEM, accessed via Google Earth Engine (GEE). |
SPI |
Standardized Precipitation Index calculated to assess drought and wet conditions, derived from daily precipitation data with a spatial resolution of 0.25 degrees. |
2000-2024 |
NOAA PERSIANN-CDR (Climate Data Record), accessed via Google Earth Engine (GEE) (Chang et al., 2022). |
LST |
Land Surface Temperature data used to model temperature variations and their impacts on water and vegetation, derived from thermal infrared bands with a resolution of 30 m. |
2020-2024 |
Landsat 8 TIRS and predicted with indices from Sentinel-2 Level 2A, processed via Google Earth Engine (GEE). |
WIW |
Water in Wetland index used to detect water pixels in vegetated areas, leveraging NIR and SWIR bands with logical thresholds, applied to map water extent. |
2020-2024 |
Landsat 8 Level 2 and Sentinel-2, processed via Google Earth Engine (GEE) (Lefebvre et al., 2019). |
NDVI |
Normalized Difference Vegetation Index used to assess vegetation health, derived from NIR and red bands, with a spatial resolution of 500 m. |
2019-2023 |
MODIS MOD13A1.061 (Terra Vegetation Indices 16-Day Global 500m), accessed via Google Earth Engine (GEE). |
NDWI |
Normalized Difference Water Index used to detect water bodies, calculated using green and NIR bands, applied to map water distribution in the wetland. |
2020-2024 |
Landsat 8 Level 2, processed via Google Earth Engine (GEE). |
MDWI (MNDWI) |
Modified Normalized Difference Water Index, a variant of NDWI using SWIR bands (e.g., SWIR2), applied to enhance water detection in vegetated areas. |
2020-2024 |
Landsat 8 Level 2, processed via Google Earth Engine (GEE) (Bourletsikas et al., 2023; Laonamsai et al., 2023). |
Precipitation (CHIRPS) |
Daily precipitation data used for trend analysis and SPI calculation, with a spatial resolution of 0.05 degrees. |
1981-2024 |
Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS), accessed via Google Earth Engine (GEE) (Zhang et al., 2025). |
Precipitation (GPM) |
Monthly precipitation data used for K-Means clustering to classify precipitation patterns, with a spatial resolution of 0.1 degrees. |
2000-2024 |
NASA/GPM_L3/IMERG_MONTHLY_V07 (Global Precipitation Measurement), accessed via Google Earth Engine (GEE) (Reddy et al., 2023). |
4. Results
4.1. Land Surface Temperature (LST)
Studies show that rivers drain approximately 75% of the surface water on land, underscoring the critical role of river systems in wetland hydrology (Bourletsikas et al., 2023). In agricultural regions, land surface emissivity has been identified as a crucial factor for energy balance, influencing thermal dynamics in such areas (Das et al., 2021). Also, research has demonstrated that Land Surface Temperature (LST) significantly affects surface water temperature and river flow, with higher LST values often leading to increased water temperatures and altered flow patterns (Lefebvre et al., 2019). These findings from other regions provide a comparative framework for interpreting the hydrological and thermal dynamics observed in the Usangu wetland Figure 5.
Figure 5. LST spatial variation in the study area.
The linear regression r value of 0.789 in the correlation of LST reference against prediction in Figure 6 suggests that the temperature values are significantly captured. The wetland ecosystem is positively affected by the LST due to an imbalance of the energy (Xu et al., 2011). The regression line data points exhibit significant clustering between 20˚C and 30˚C. The alteration in Land Surface Temperature (LST) trends has affected land cover changes, desertification, and water management. Figure 7 illustrates the area influenced by the increase in LST within the basin.
The r correlation in Figure 7 illustrates the MDWI (Modified Normalized Difference Water Index) spectral index, which indicates that water detection is lower than that of MNDWI, which indicates the availability of water in the wetland but is obscured by vegetation. According to the Histogram, the correlation results in the wetland exhibit identical greenness and dryness values. The negative correlation between the LST and elevation in Figure 8 indicates that there is opposite relationship between the two variables. Consequently, the wetland is characterized by a higher temperature than the rest of the areas as shown in Figure 5.
Figure 6. LST baseline (Reference) against prediction.
Figure 7. Correlation coefficient with LST.
Figure 8. LST against elevation.
4.2. Precipitation
4.2.1. STL
The precipitation data was analysed using the ML of the stats model in the statistics package. The analysis of the trend indicates a long-term increase in precipitation from the 1990s to the 2020s, characterized by below-average levels in the 1990s and subsequent growth in the 2000s. The timeseries decomposition Figure 9 indicates years below the average from 2000 to 2003 and from 2010 to 2012. This timeframe aligns with the SPI Figure 11 depicting the dry spell in the wetland area.
Figure 9. Timeseries decomposition of precipitation.
4.2.2. Standard Precipitation Index (SPI) in the Wetland Area
In K Means, the precipitation data categorised into four classifications: heavy (red), high (orange), low (blue), and moderate (green) as depicted in Figure 10. The years characterised by moderate precipitation, correlated with the SPI, are predominantly represented by green bars (moderate) and exhibit fewer red, orange, and blue bars: 2003, 2008, 2013, 2017, and 2023. Mean precipitation during the years with moderate rainfall in the study area is 731 mm. We selected 4 clusters to capture the meteorological significance of precipitation intensity distribution, achieving a silhouette score of 0.695 despite it being a lower score for the two clusters and three clusters.
The index indicates a long-term deficiency in rainfall. Figure 11 illustrates extended droughts occurring in the years 2003, 2008, 2013, 2017, and 2023. The SPI drought trend of the data is compared with Figure 9, which represents the precipitation dataset using K-Means.
Figure 10. K-means precipitation temporal clustering.
Figure 11. Standardized precipitation index (SPI) with drought years.
4.3. NDVI of the Vegetation Cover
The Usangu wetland, an agricultural area characterized by extensive irrigation development, experiences a drop in NDVI from April to May, coinciding with the paddy harvesting season. This seasonal decline explains the negligible correlation between drought anomaly years in the SPI index (shown in Figure 11) and NDVI values in the wetland in Figure 13. For instance, the drought year 2023 had a mean NDVI of 0.50, reflecting minimal impact. NDVI exhibits a decreasing annual trend from April to October, coinciding with rising LST in the study area shown in Figure 14, which peaks during this period. The NDVI map Figure 12 further illustrates the spatial distribution of vegetation across the region.
NDVI and NDWI of the Wetland
The NDVI Figure 15(a) indicate that the wetland is covered with healthy vegetation, with the maximum positive NDVI reaching 0.5. Normalised Difference Water Index (NDWI) (Xu, 2006) Figure 15(b) indicates a maximum negative value of −0.25, suggesting that the wetland is vegetated but deprived of water.
Figure 12. NDVI map showing spatial variation in the wetland area.
Figure 13. NDVI values show seasonal changes in vegetation health across the wetland throughout the year.
Figure 14. LST show seasonal changes in temperature across the wetland area throughout the year NDWI.
Figure 15. (a) Usangu Wetland NDVI; (b) Usangu Wetland NDWI.
4.4. Analysis of Water Availability in the Wetland Using ANDWI
The Modified Normalized Difference Water Index (MDWI) (Acharya et al., 2018), a variant of the Normalized Difference Water Index (NDWI) that incorporates the SWIR2 band, was applied to assess water and vegetation dynamics in the Usangu Wetland. Analysis of Figure 17(a) (ANDWI) and Figure 17(b) (NDWI) indicates water index values between −0.6 and −0.3, suggesting mostly non-water features within the wetland, as shown in Figure 16. This range of values, particularly the prevalence of negative NDWI values, indicates that the wetland areas are characterized by healthy vegetation, as negative values typically reflect higher vegetation presence due to the influence of the SWIR2 band in reducing water detection sensitivity in vegetated regions.
4.5. Water in Wetland Method (WIW)
The Water in Wetland (WIW) approach (Lefebvre et al., 2019), using single spectral bands from Landsat and Sentinel data alongside logical thresholds, was used to identify water pixels in the wetland area. Figure 18(a) and Figure 18(b) indicate that this methodology identified water pixels, hence demonstrating the presence of water
Figure 16. Histogram of wetland NDWI.
Figure 17. (a) ANDWI values in wetland indicate water coverage in the study area; (b) NDWI values in the wetland.
Figure 18. (a) WIW values using a relaxed classification logic, show the pixel-based water extent in the wetland study areas; (b) WIW values using a restricted classification logic, show the pixel-based water extent in the wetland study areas.
bodies within the wetland. Subsequent analysis in Figure 18(a) reveals an increased pixel density, signifying that the WIW method with moderate logical criteria detected a greater concentration of water pixels in contrast to a more restrictive approach in Figure 18(b).
The Water in Wetland (WIW) method employs spectral thresholds applied to Near-Infrared (NIR) and Short-Wave Infrared 2 (SWIR2). Two distinct logical conditions govern water detection. The stricter condition equation 4a requires both thresholds to be satisfied simultaneously and is best suited for mapping permanent water bodies where precision is critical. The relaxed condition equation 4b requires at least one threshold to be satisfied and is best suited for wetland delineation where omission errors must be minimised.
5. Conclusion
The presented analysis of the wetlands condition, including WIW, NDVI, LST, and SPI, can be applied for various purposes such as water resource management, biodiversity conservation and agricultural planning.
A statistically significant warming trend in the study area, as indicated by rising LST values, poses a risk to wetland management. Higher temperatures may increase evaporation rates, reduce water availability, and stress the wetland ecosystem.
A decrease in greenness, reflected by declining NDVI values across multiple spectral bands, suggests a reduction in vegetation health within the wetland. This decline negatively impacts biodiversity by reducing habitat quality and food availability for wetland species.
The negative correlation coefficient (r) between LST and Digital Elevation Model (DEM) values directly affects the wetland. This negative correlation, as noted in Figure 8, indicates that lower elevations experience higher temperatures. Since the wetland is situated at lower elevations where water accumulates, as supported by the WIW analysis, this warming trend exacerbates heat stress in the wetland, potentially leading to increased evaporation and reduced water levels, which further threatens the ecosystem.
Indicators derived from wetland analysis, including Land Surface Temperature (LST), Standardized Precipitation Index (SPI), vegetation indices (NDVI), and water indices (MNDWI, NDWI, WIW), provide critical data for future environmental planning in the region. However, comprehensive Land Use/Land Cover (LULC) analysis is required to identify specific threats to wetland sustainability under prevailing trend of the Climate Change.
Acknowledgements
We acknowledge NASA/USGS for providing the Landsat satellite images used in this study. We also thank the Department of Earth and Climate Science of the University of Nairobi for the technical editing of this manuscript.