Next Article in Journal
Experimental Analysis of Rock Boulder Impacts on Brick Walls to Support Numerical Modelling of Building Damage
Previous Article in Journal
The Origins of the Hydrogen Sulphide (H2S) Gas in the Triassic Montney Formation, British Columbia, Canada
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Variability Analysis of Rainfall and Water Quality: Insights from Trend Analysis and Wavelet Coherence Approach

by
Syeda Zehan Farzana
1,2,*,
Dev Raj Paudyal
1,*,
Sreeni Chadalavada
2 and
Md Jahangir Alam
2,3
1
School of Surveying and Built Environment, University of Southern Queensland (UniSQ), Toowoomba, QLD 4350, Australia
2
School of Engineering, University of Southern Queensland (UniSQ), Springfield Lakes, QLD 4300, Australia
3
Murray-Darling Basin Authority (MDBA), Canberra, ACT 2601, Australia
*
Authors to whom correspondence should be addressed.
Geosciences 2024, 14(8), 225; https://doi.org/10.3390/geosciences14080225
Submission received: 9 July 2024 / Revised: 9 August 2024 / Accepted: 17 August 2024 / Published: 21 August 2024
(This article belongs to the Section Climate)

Abstract

:
An understanding of the trend and relationship between rainfall patterns and water quality dynamics can provide valuable guidelines for the effective management of water resources. The aim of this study was to reveal the synchronous trends in rainfall and water quality and to explore the potential connection between seasonal variation in rainfall volume and the water quality index. This study scrutinised the seasonal temporal trends of rainfall and water quality parameters of three water supply reservoirs in the Toowoomba region of Australia by applying the modified Mann–Kendall (MMK) test and innovative trend analysis (ITA) methods from data collected over 22 years (2002–2022). The models showed a significant increasing trend of rainfall in two rainfall stations during autumn season. The water quality parameters, such as PO43−, exhibited a significant decreasing trend in all seasons in three reservoirs. On the other hand, the water quality index (WQI) showed a decreasing trend in the Cooby and Cressbrook reservoirs, excepting the Perseverance reservoir, which exhibited an increasing trend. In addition to the detection of trends, this study investigated the potential correlation between seasonal variation of rainfall volume and the water quality index using the wavelet transform coherence (WTC) method. The data of twelve rainfall stations were brought into this analysis. The WTC analysis displayed an apparent correlation between the water quality index and rainfall pattern for 70% of the rainfall stations across 8–16 periods. The highest coherency was noticed in 8–16 periods from 2002–2022, as observed at both the Cooby Creek rainfall station and in the WQI of the Cooby reservoir. This evaluation revealed the intertwined dynamics of rainfall patterns and water quality, providing a deeper understanding of their interdependence and implications, which might be useful for environmental and hydrological management practices.

1. Introduction

Rainfall plays a crucial role in hydrological processes, and its spatial and temporal distribution, along with evolving characteristics, can significantly affect agriculture, ecosystems, and the management of water resources [1]. The understanding of monthly, seasonal, and annual variability of rainfall and as well as the severity of extreme rainfall and how it affects agriculture, ecosystem and water resources is crucial to ensure proper water resource management practices [2]. Rainfall patterns have changed significantly over recent years due to climate change leading to frequent flooding and drought. These events, in turn, have disrupted the ecosystem of water bodies and impacted the water quality [3]. With growing awareness about the escalating impact of climate change, it is important to grasp the characteristics and variations of rainfall, as well as how these patterns shift across various spatial and temporal scales. Previous studies [4,5] have reported that there is global as well as regional scale change in rainfall patterns caused by global climate change.
Inland waters, such as rivers, lakes, reservoirs, and streams, provide diverse ecosystem services, from food sources to transportation channels and recreational sites. The deterioration of water quality due to the impact of extreme rainfall and rapid socio-economic development poses a significant environmental concern for inland water bodies [6]. Increases in extreme hydrological events and circulating (air) temperature are the dominant factors affecting water quality. Prolonged and heavy precipitation after drought periods then causes runoff and erosion that together consequently increase the pollutant release into surface water bodies [7]. Efforts to explore the efficacy of previous initiatives and the primary factors influencing the quality of inland water must continue to be thoroughly examined across different temporal and spatial dimensions. Among various nutrients present in water bodies, Nitrogen (N) and phosphorus (P) are two crucial elements that may disrupt the balance of aquatic ecosystems, thereby influencing diverse services such as navigation, water sports, fisheries and, most importantly, drinking water supply [8].
Time series data analysis is crucial for understanding the dynamics of events, which can play a key role in strategic infrastructure planning and early warning systems [9]. Hydroclimatic time series can display trends, whether upward or downward, and can locate change points that may relate to climate change influences, like increasing greenhouse gas levels, and alterations in land use and management practices [10,11,12]. To consider adaptation and mitigation strategies to counteract the effects of climate change and making informed decisions, it is essential to gain better understanding of trends and variability of hydroclimatic variables [13]. Precipitation concentration index (PCI), seasonality index (SI) and coefficient variation (CV) were extensively applied to examine the variability of rainfall, temperature and streamflow and to analyse the spatial and temporal distribution [14,15,16].
Over the past few years, there has been a growing need to periodically measure various water quality variables to monitor the status of surface water [17]. The steady collection of water quality records can provide a reliable synopsis of water quality status, identifying existing or potential water quality issues, understanding their causes, and evaluating the effectiveness of any implemented measures [17,18]. To date, numerous studies [19,20] have been conducted to explore the influences impacting changes in water quality such as pollution from agriculture, sewage discharge, socio-economic advancements, water consumption, and forest and climate extremes. Among all other variables, climate extremes are expected to significantly influence future hydrological patterns, affecting various sectors that depend on surface water. Delpla et al. [21] showed that the amount of organic matter was likely to increase in streams that experienced an increased volume of rainfall and alterations in drought-rewetting cycles. Consequently, Ponting et al. [22] discovered that inorganic and organic matter were carried into rivers by flood water, whereas Ockenden et al. [23] showed a positive trend and correlation between annual rainfall and total phosphorus load. Moreover, Mortazavi-Naeini et al. [24] have demonstrated a positive relationship between extreme winter rainfall and suspended sediment concentration and an elevated level of phosphorus during the low flow period. Additionally, Watts and Anderson [25] observed that a reduction in summer flow is expected to decrease the dissolved oxygen level and conversely, high summer floods to enhance the concentration of nutrients and pollutants in the river system. Nevertheless, despite the above studies offering basic theoretical frameworks and empirical illustrations of the source of water quality change, there are still notable shortcomings. Because of inadequate data sources and diverse research objectives, the systematic formation of the correlation between water quality parameters and key driving factors has not yet been undertaken. This limitation hinders support for the prevention of water quality deterioration in areas facing significant water quality issues [26,27].
In recent decades, various methods have been developed and used for analysing trends in time series such as the Mann–Kendall (MK) test, seasonal Kendall methods, the modified Mann–Kendall (MMK) test, linear regression analysis, Spearman’s rho (SR) test, the Theil–Sen approach (TSA), Sen’s slope, innovative trend analysis (ITA) and the moving trend method (MTM) [28,29,30,31,32,33,34,35,36]. The climate variability and its trends can be assessed through parametric (t-test, F-test, linear regression) and non-parametric (Mann–Kendall test, Sen’s slope estimator) methods. The application of parametric methods is constrained to time series that follow a normal distribution. The climate time series, specifically precipitation, do not meet the normality requirement and thus the non-parametric method is applied for trend analysis [37], which has been followed in this study. In the case of data with outliers, non-parametric methods are deemed more robust than parametric methods [38]. Although a traditional method, the Mann–Kendall (MK) test, is the most widely utilised approach, the use of this method is constrained by assumptions such as the serial independence of time series data [39] and the inability to analyse trend behaviours across different categories (low, medium and high) of time series [40,41]. The MMK test is an enhancement of the traditional Mann–Kendall test, which offers a statistically rigorous approach to trend the detection of data with autocorrelation [42]. ITA provides visual inspection and allows for the categorisation of trends and is suitable for comprehensive trend analysis [40].
In Australia, the frequency and volume of rainfall are the primary constraining factors in water security levels [43]. There are limitations in the monitoring and reporting of water quality in regional towns. Wyrwoll et al. [44] have reported that there is limited water quality parameter testing and reporting in 24 local councils. The data and reports are inconsistent and do not comply with the Australian Water Quality Guidelines. Furthermore, Queensland government regulations do not require water utilities to furnish thorough quantitative data analysis and reporting for each parameter [7]. Thus, there is a necessity to gain a deep understanding of both the quantitative and qualitative aspects of rainfall and water quality variations, offering a detailed understanding for water quality management. Moreover, the study covering small scale catchments can minimise natural and anthropogenic variability that may evade the impacts of averaging across larger areas [45]. There is a lack of research on the observation of combined trends and correlation of rainfall and water quality. Therefore, a multi method approach of trend and correlation analysis can provide a holistic understanding of the shifts of hydrological and environmental variables.
The aim of this study is to examine the trends of rainfall and water quality parameters and to explore the correlation between the signals of precipitation and water quality parameters. This analysis was carried out for three water supply reservoirs (Cooby, Cressbrook and Perseverance) in the Toowoomba region of Australia. In this study, two trend analysis methods, namely the modified Mann–Kendall (MMK) test and innovative trend analysis (ITA), were applied for rainfall and water quality trend analyses, respectively. Previously, the water quality index (WQI) was calculated based on five water quality parameters—pH, turbidity, phosphate (PO43−), ammonia nitrogen (NH3-N) and total dissolved solids (TDS) [7]—which have been correlated in this study with rainfall patterns through investigation using the wavelet transform coherence (WTC) method. The explanation and reasons behind the selection of five parameters are elaborated on in our previous study [7]. In summary, the specific objectives are as follows:
  • To identify the seasonal trend of rainfall and water quality parameters using a modified Mann–Kendall test (MMK) and innovative trend analysis (ITA), respectively, and to show the spatial distribution of rainfall trend parameters.
  • To detect the temporal patterns and correlation between rainfall and water quality index (WQI) using wavelet transform coherence (WTC) spanning a period of 22 years.

2. Materials and Methods

2.1. Study Area

The Toowoomba region is situated in the Darling Downs, Queensland, Australia [46], at an elevation ranging from 125 m to 740 m with an average of approximately 450 m above mean sea level (MSL) [47]. The study area is located in part of the mountain chain that forms the Great Dividing Range of eastern Australia [48]. The average annual rainfall of this region is 735 mm, which peaks in the warm season [49]. The average annual temperature is 19°, while the average summer temperature is 25 °C, and the average winter temperature is 12 °C. Average temperatures across the region are currently 1 °C higher than they were 100 years ago [50]. Three water supply reservoirs, namely Cooby (27.3825° S, 151.9244° E), Cressbrook (27.2641° S, 152.1959° E) and Perseverance (27.2883° S, 152.1239° E), serve as the main source of water supply in this region. The water supply capacity of the three reservoirs is 19.7 megalitres (ML), 78.8 megalitres (ML) and 26.9 megalitres (ML), respectively, and the storage area is 306 ha, 517 ha and 250 ha, respectively. The total catchment area of Cressbrook, including Perseverance reservoir, is 320 km2 and that of the Cooby catchment is 159 km2 [51]. The catchments are situated in a warm/humid climatic zone of subtropical Australia, featuring distinctly cool, dry winters and warmer, wetter summers [7]. The topography of the catchments features gentle slopes at lower elevations and hills at higher elevations [48]. Twelve rainfall stations surrounding the three dam reservoir catchments were considered in order to observe the trends of rainfall within this region. The selected rainfall stations and water supply reservoirs are illustrated in Figure 1.

2.2. Data Collection and Compilation

The water quality data of three dam reservoirs were collected from the Toowoomba Regional Council (TRC), which is responsible for the management of the dam reservoirs and potable water supply in the Toowoomba region. Twenty-two years (2001–2022) of weekly water quality data were collected from TRC. Five water quality parameters affected by rainfall runoff, such as pH, turbidity, total dissolved solids (TDS), ammonia nitrogen (NH3-N) and phosphate (PO43−), were considered for computing the Water Quality Index (WQI) [7]. In conjunction with the water quality data collection, daily rainfall data were collected over the same time span from the Bureau of Meteorology (BOM), Australia (http://www.bom.gov.au/, accessed on 7 March 2024).
There are four seasons in Australia (summer, autumn, winter, and spring) [52]. The summer season includes the three hottest months December, January and February and the winter season includes the three coldest months June, July and August. The trends of water quality parameters and rainfall were scrutinised across four seasons throughout the timeframe (2001–2022) using the modified Mann–Kendall (MMK) test and the innovative trend analysis (ITA) method. The combined use of these two methods provided a comprehensive understanding of seasonal trends (both monotonic and nonmonotonic) of rainfall and water quality. The dataset was checked, and missing values were handled appropriately. Additionally, wavelet transform coherence (WTC) was utilised to examine the correlation between rainfall and WQI. The trend analysis and wavelet transform were conducted using R 4.3.2 software, and the spatial distribution of rainfall was analysed in ArcGIS Pro software v3.03. The applied techniques are described in the following sections, with the flow diagram showing the research method (Figure 2).

2.3. Precipitation Concentration Index (PCI)

The term precipitation concentration index (PCI) was proposed by Oliver [53] to enumerate the temporal variability of rainfall, also used in soil erosion studies [54]. This index, generally applied on annual scales, represents the monthly precipitation and precipitation heterogeneity within a year. A value between 8.3 and 10 indicates uniform distribution, with a moderate level of uniformity in the distribution of precipitation falling in the range from 10 to 15, and values from 15 to 20 indicating an irregular distribution of precipitation. A value greater than 20 characterises a strong concentration of irregular precipitation over a limited period [53]. In this study, the PCI values were calculated using Equation (1) in Table 1 to examine the spatial and temporal variability of precipitation applied to monthly precipitation [1] data from 2001 to 2022.

2.4. Coefficient of Variance (CV)

The coefficient of variance (CV) or normalised root mean square deviation is a relative measurement of deviation of frequency distribution. It is calculated by dividing the standard deviation of a series by the mean of the series [55]. The formula to calculate CV (Equation (2)) is provided in Table 1. The degree of variability of rainfall is classified by the coefficient of variation (CV) [56]. The rainfall events’ variability is classified as high (CV > 30), moderate (20 < CV < 30) or low (CV < 20) [56]. CV was computed in this study to assess the relative annual variability of rainfall in the region.

2.5. Mann–Kendall (MK) and Modified Mann–Kendall (MMK) Tests

The Mann–Kendall (MK) test is a non-parametric statistical test applied to determine the monotonic upward and downward trend in time series data [57]. This test exhibits upward and downward trends with statistical significance [58]. The MK statistics increase by one if the data point of a year is greater than the preceding year and, if the converse occurs, it decreases by one. This test performs well if there is no autocorrelation in the data. This test provides a test variable (S) and a variance (Equations (3) and (4)) in Table 1.
The monotonic trend in the time series can be determined by the Z value. The positive value indicates an increasing trend, whereas the negative value suggests a decreasing trend. The statistical significance of the trend is determined by the specific p-value, with a value less than 0.05 indicating statistical significance [59].
The use of the MK method is constrained by autocorrelation effects in data. If there is autocorrelation in the time series, the MK method may indicate a trend if no actual trend exists [42]. A sequential Mann–Kendall test is used to define the starting time of a significant trend in time series [60]. Whereas, in the MMK test, the modified VAR(S) is applied to detect the trend, as shown in Equation (5) (Table 1). In this study, the MMK test was applied to assess the seasonal trends of both rainfall and water quality parameters. The significance of a trend was determined based on a p-value. The trend with a p-value less than 0.05, which corresponds 5 percent uncertainties or less, was considered to be a significant trend.

2.6. Sen’s Slope Estimator

The magnitude of the trend slope is estimated by a non-parametric procedure known as Sen’s slope estimator [61]. The equation of Sen’s slope (Equation (6)) is provided in Table 1. In this study, Sen’s slope estimator was used to quantify the amount of change in a season per year.

2.7. Innovative Trend Analysis (ITA)

The innovative trend analysis (ITA) is a non-parametric graphical method first introduced by Şen [62]. This method is capable of handling autocorrelation and outliers within time series data, detecting both monotonic trends and sub-trends, and identifying combinations of trends across different periods [63]. In the ITA method, firstly the times series is divided into two sub-series each containing an equal number of observations. Subsequently, both series are sorted in ascending order and plotted against each other in a Cartesian coordinate system. In the plot, the first half is plotted on the X axis while the second half is plotted on the Y axis. In the later phase, a straight line is fitted in the scatter plot, which displays either a monotonic trend or no trend. If most scattered points lie above the 1:1 (45°) line, there is an increasing trend in the time series; conversely, if the points cluster below the 1:1 line, this indicates a decreasing trend. Similarly, no trend is signified in the data series, if the points clustered along the trend line [62]. A trend detector was used to examine whether the trend was increasing or decreasing (Equation (7) in Table 1). The ITA method calculations were performed on both rainfall and water quality data series. The innovative trend detection (ITD) was shown as increasing or decreasing using up arrows or down arrows, respectively. Another approach, innovative polygon trend analysis (IPTA) has also been proposed by Şen et al. [64]. This can identify the trend in a particular sequence by providing a trend polygon for the linguistic and numerical interpretation of trends. However, we applied the ITA method in this study to easily identify straightforward trends and shifts rather than using complex visualisation techniques.

2.8. Wavelet Analysis

Wavelet analysis is a dominant analytical tool for analysing processes exhibiting a time–frequency representation of signal in the time domain [65]. Wavelet transform provides localised frequency decomposition, offering insights into frequency components, and, as a consequence, wavelets have a notable dominance over Fourier analysis [66]. In this study, continuous wavelet tools, particularly wavelet coherence, were applied to measure and examine the degree of local correlation between two time series of rainfall and water quality indices in the time–frequency domain, along with wavelet coherence phase differences. The fundamental concepts of wavelet and wavelet transform coherence are explained in the following sections.

2.8.1. Wavelet

A wavelet can be defined as per [66]. It is a true valued square integrable function and can be interpreted as illustrated in Equation (8) (Table 1). A wavelet consists of two control parameters, location, and scale parameter. The location parameter (u) determines the precise position of the wavelet, whereas the scale parameter(s) defines the extent of the expansion of the wavelet. Scale is inversely related to frequency; therefore, a lower scale means a more compressed wavelet, and a higher scale means a less compressed wavelet. A more compressed wavelet can detect higher frequencies of a time series [66].

2.8.2. Wavelet Transform Coherence (WTC)

A bivariate framework known as wavelet coherence explores the interaction between two time series. Wavelet transform coherence (WTC) locates the sections in the time–frequency space where two distinct time series are certainly not powerful. The wavelet coherence of the time series is calculated by using Equation (9) in Table 1. The value of the squared wavelet coherence coefficient (Rn in Equation (9)) represents the correlation [66]. Values close to 1 indicate a strong correlation whereas values close to 0 exhibit a weak correlation. Thus, the squared wavelet coherence estimates the local linear correlation between two time series that is homogeneous to the squared correlation coefficient in linear regression [66]. The equations of the proper smoothing operator (S) for Morlet wavelet coherence are provided in Table 1 (Equations (10)–(12)). For the Morlet wavelet, a factor of 0.6 is experimentally designated as the scale decorrelation length [67,68,69].
In wavelet coherence plots, the thick black contour represents the 5% significance level. The cone of influence (COI) is the region in the plot where the edge affects the wavelet power (lighter shading in the plot) and the phase is indicated by arrows. Zero phase difference means the analysed time series move together. Arrows pointing to the right indicate in-phase, while those pointing to the left indicate anti-phase. Arrows pointing up indicate that the first time series leads the second one by 90° and arrows pointing down indicate that the second one is leading the first time series. The stronger correlation is displayed as red and yellow colours in the plot and the blue colour specifies a weak correlation. In the present study, a combination of arrows pointing in different directions was observed in the plot.
Table 1. List of equations used in the analysis.
Table 1. List of equations used in the analysis.
EquationSN.DescriptionsReferences
PCI annual = i = 1 12 Ri 2 ( i = 1 12 Ri ) 2 × 100 (1)PCI = Precipitation concentration index
Ri = Rainfall in month i.
[54]
  CV = σ μ (2)σ = Standard deviation
μ = Mean of the series
[55]
S = i = 1 n 1 j = i + 1 n sign X j   X i (3)n = Number of data; Xi, Xj = data points of time series i and j[58]
VAR S = 1 18   n n 1 2 n + 5 k = 1 m t k t k   1 2 t k + 5 (4)N = Number of data points
m = Number of tied groups
tk = Number of ties in kth group
[58]
VAR S = n n 1 2 n + 5 18 . n n e * (5) n n e * = Correction factor to adjust the auto correlated data[42]
Q i   = x j x k j k (6)xj, xk = Values of data at time j and k (j > k)
i = 1, 2, 3…N
[61]
B = 1 n   i = 1 n 10   x j x k x ¯ (7)B = Trend detector
n = Size of individual sub series
xj, xk = The values of sub series x ¯ = Average of first sub series xk.
[62]
ψ u , s   t = 1 s   φ   t u s (8)1/√s = A normalisation factor ensuring the unit variance of the wavelet, with || ψ u , s   ||2 = 1
u = Location parameter
s = Scale parameter
[66]
R n 2 u , s = | S   ( S 1 W n XY   u , s | 2 S   S 1   W n x u , s 2   S   S 1 W n y u , s 2 (9)Rn = Squared wavelet coherence coefficient ranging from 0 to 1
n = Time index
S = Smoothing operator
[65,66]
S W =   S scale   S time   W n S (10)Sscale = Smoothing along wavelet scale
Stime = Time parameter
[65,66]
S time   W   S = W n   S   ·   S 1 1 2 2 s 2 S (11)Stime = Time parameter[65,66]
S time   W   S = W n   S   ·   C 2 Π   0.65 (12)C2 = Normalisation factor, and Π = rectangle function[65,66]
In this study, the MMK and ITA tests were applied to rigorously detect trends in rainfall and water quality. The MMK test was chosen for its robustness in identifying statistically significant trends while accounting for serial correlation and seasonality, making it ideal for environmental data with autocorrelation. Meanwhile, the ITA was utilised for its visual and intuitive approach. Additionally, WTC was employed to examine the correlation between rainfall and water quality, capturing the time and frequency dimensions of their relationship.

3. Results and Discussion

3.1. Rainfall Distribution

The basic annual rainfall statistics, such as mean, standard deviation (SD), coefficient of variation (CV) and precipitation concentration index (PCI), of 12 rainfall stations around the three dam catchments were analysed from 2001 to 2022. The mean of the long-term annual rainfall varied from 575 mm to 890 mm. The range of SD and CV was 171 mm to 350 mm and 28% to 39%, respectively. The highest mean rainfall was found at Mount Kynoch station, where SD and CV values were also at the highest levels. As discussed in Section 2.4, CV (%) greater than 30 was classified as high variability and between 20 to 30 was classified as moderate variability of rainfall. We found that the value of CV (%) was greater than 30 in 8 rainfall stations, which indicated high variability of rainfall, and the data from the rest four stations could be classified as moderate to highly variable. The PCI value ranged from 15–20, except for Pechey Forestry, which was 14.2 indicating an irregular distribution of precipitation in the region. The findings are presented in Table 2.
The inverse distance weighting (IDW) method was applied to show the spatial distribution of rainfall in this study. The spatial distribution of mean, CV and PCI is illustrated in Figure 3. The mean rainfall varied from 813 to 890 mm in the stations near the dam catchment area. The value of CV was high (39%) in Mount Kynoch and low (28–30%) in Oakey Aero, Perseverance Dam and Cooby Creek. The PCI value was in the range of 14–15 (moderate) in the northeast side of the study area and from 15.1–18 (irregular) on the other side.

3.2. Auto Correlation Function (ACF)

Table 3 presents the values and plot of lag 1 autocorrelation in seasonal rainfall data for 22 years (2000–2022). In detecting the trend in time series data, the calculation of autocorrelation is a primary step, the magnitude of sequential correlation increases with the variance of the MK test statistic. If there is a positive sequential correlation in the time series data, this amplifies the likelihood of type I errors (false positives) and identifies a noteworthy trend even if there is no genuine trend [70]. In this study, a lag 1 autocorrelation coefficient was calculated at a 95% significance level which was 0.4. It is found that the autumn season rainfall data of Doctors Creek shows autocorrelation at a 95% significance level (highlighted in Table 3). Because autocorrelation exists in the data series in one station, the modified Mann–Kendall (MMK) test was applied to detect the trend of seasonal rainfall.

3.3. Seasonal Rainfall Trends

3.3.1. Autumn Trend

In scrutinising the trend of rainfall in the autumn season for the entire period of 2001–2022, distinct patterns and fluctuations in precipitation are revealed (Table 4). The Z values showed increasing rainfall trends in all stations, varying between 0.75 to 2.26 and the maximum is at Cooby Creek station. The trend was significant in Cooby Creek and Moyola stations (p < 0.05).
Similarly, the ITA method showed an increasing trend. The indicator values of ITA varied from 1.55–7.38, with the highest at Perseverance Dam (7.38). The detection of the trend is shown by using upward or downward arrows (ITD column) for simplicity.
All other coefficients, such as MK tau Sen’s slope and ITA indicator, showed increasing trend (Table 4). According to the value of Sen’s slope, the highest change of rainfall in this season was 7.86 mm/yr in Mount Kynoch station whereas according to ITA, Perseverance Dam showed the highest significant trend.
Figure 4 shows the spatial distribution of autumn rainfall trend parameters. Similar and consistent patterns were observed across the parameters, with higher values concentrated in the lower region and lower values in the upper region, which suggests that lower areas were more affected in the region.
Table 4. Rainfall trend (autumn).
Table 4. Rainfall trend (autumn).
StationZ StatisticsMK Taup ValueSen’s Slope
(mm/yr)
ITDB (ITA)
Cooby Creek2.260.300.024.89Geosciences 14 00225 i0013.39
Cressbrook Dam1.000.210.724.77Geosciences 14 00225 i0011.55
Doctors Creek1.400.300.165.57Geosciences 14 00225 i0013.81
Glenaven1.500.270.134.29Geosciences 14 00225 i0012.43
Goombungee PO1.880.310.075.84Geosciences 14 00225 i0014.51
Haden PO0.750.150.453.56Geosciences 14 00225 i0013.63
Moyola2.200.300.035.61Geosciences 14 00225 i0014.90
Mount Kynoch1.790.280.077.86Geosciences 14 00225 i0015.24
Oakey Aero1.110.210.274.10Geosciences 14 00225 i0013.69
Pechey Forestry1.310.260.195.32Geosciences 14 00225 i0014.12
Perseverance Dam1.360.250.176.34Geosciences 14 00225 i0017.38
Tamba1.690.300.097.03Geosciences 14 00225 i0016.36
N.B: ITD = Innovative trend detection. The red arrow indicates an increasing trend. A significant trend is highlighted by the yellow colour, (p value < 0.05).

3.3.2. Winter Trend

The MMK and ITA methods revealed the same scenarios regarding the winter rainfall trend in the region. Of 12 stations, a non-significant increasing trend was observed in 8 whereas a decreasing trend was found in the remaining 4 stations. The value of Z statistics and ITA indicator were higher at Pechey Forestry, at 0.75 and 0.77 respectively. In Haden PO, the value of Z statistics was minimum (−1.10) indicating a decreasing trend. The average change of rainfall was 0.16 mm/yr. Table 5 illustrates rainfall trends in winter and Figure 5 shows the spatial distribution of winter rainfall trend parameters

3.3.3. Spring Trend

The rainfall trend in the spring season was characterised by low rainfall trends. Cooby Creek and Glenaven exhibited negligible trends, with Z statistics close to zero, suggesting stable conditions. Cressbrook Dam, Doctors Creek, Goombungee PO, Moyola, Mount Kynoch, Oakey Aero, and Pechey Forestry implied increasing trends according to the value of Z statistics and ITA indicators. Based on the value of Sen’s slope, the average change of rainfall was 0.57 mm/yr. Conversely, Haden PO and Tamba demonstrated negative trends. Additionally, ITA indicator values indicated low-to-moderate variability in precipitation trends. Table 6 illustrates the rainfall trend in the spring season and Figure 6 shows the spatial distribution of spring rainfall trend parameters.

3.3.4. Summer Trend

In examining the summer season’s rainfall trend, MMK analysis shows a non-significant increasing trend in seven stations and a decreasing trend in the remaining five stations, as illustrated in Table 7 and Figure 7. Notably, Haden PO stood out, with a negative Sen’s slope and Z statistics, indicating a substantial decrease in precipitation. The average change of rainfall in this was 0.13 mm/yr. The value of the ITA indicator varied from −2.11 to 2.25.

3.4. Trend of Seasonal Water Quality

The aim of this analysis was to elucidate temporal variations in water quality across four distinct seasons spanning 22 years (2001–2022). This comprehensive study examined trends in five water quality parameters—ammonia nitrogen (NH3-N), pH, phosphate (PO43−), total dissolved solids (TDS), and turbidity along—with the water quality index (WQI), which was determined based on these five parameters. The selection of parameters, computation of WQI, and its analysis are discussed in detail in a previous study [7]. The WQI trends detected by the MMK and ITA method. The findings are illustrated in the following sections.

3.4.1. Cooby Reservoir

The trend results of both MMK and ITA provided a combination of both increasing and decreasing trend of water quality parameters. Out of five water quality parameters, only PO43− exhibited a significant decreasing trend in all four seasons. There was no change in the concentration based on the value of Sen’s slope. The trend of NH3-N was non-significantly increasing, pH showed an increasing trend in winter and spring and a decreasing trend in autumn and summer. The TDS showed almost no trend in winter but an increasing trend in the other seasons. The turbidity of water exhibited a non-significant increasing trend in all four seasons.
As shown in Table 8 and the WQI of this reservoir exhibited a decreasing trend in all four seasons. The average change was −1.74 according to the value of the ITA indicator. In WQI calculation [7], the maximum weightage (5 on a scale of 1–5) was assigned to NH3-N and PO43−, and the minimum weightage 3 was assigned to turbidity. Both TDS and pH had a weightage of 4. The relative effects were observed with WQI trend, with an increasing trend of NH3-N, TDS and pH in autumn and summer, resulting in an overall decreasing trend of WQI. Hence, the effects of the decreasing trend of PO43− might not have much impact on the WQI trend. The ITA plots of Cooby reservoir is shown in Figure 8.

3.4.2. Cressbrook Reservoir

For Cressbrook reservoir, it was observed from MMK and ITA outcomes that NH3-N had an increasing trend in all seasons (Table 9). pH and turbidity demonstrated a non-significant increasing trend in autumn, winter and summer and a decreasing trend in spring. However, PO43− exhibited a significant declining trend entirely in all four seasons, where the minimum Z statistics value was noted in autumn (−4.88). In the case of Cooby, turbidity showed an increasing trend in all seasons, but in Cressbrook, the trend was decreasing in spring season. Like Cooby, TDS showed no trend in winter. Furthermore, WQI showed a non-significant declining trend across all four seasons. There was a slight increase in spring in the second half of the data (2012–2022) (Figure 9). In relation to the WQI of Cooby reservoir, the increasing trend of NH3-N, turbidity and pH also affected the overall trend of WQI of this reservoir. Based on the value of Sen’s slope, the WQI is decreasing by −0.01 on average per year (Table 9).

3.4.3. Perseverance Reservoir

A noticeable difference in trend was observed in Perseverance reservoir compared with the other two reservoirs. Most significantly, WQI showed a significant increasing trend in three seasons, except for spring, while the other two exhibited a declining trend. This matched with the pH trend as non-significantly decreasing in summer while increasing in the other three seasons. Similarly, PO43− revealed a significant decreasing trend, while TDS showed a significant increasing trend in winter and spring and turbidity was significantly increasing in spring. In comparison with the WQI of the other two reservoirs, the WQI was very poor (<10) in this reservoir [7]. The slight increasing trend was positive for this reservoir water quality and the average increase was 0.06, based on the value of Sen’s slope (Table 10). The WQI trends detected by the ITA method is shown in Figure 10 and Figure 11 is the Z statistic plot of the water quality parameters of three reservoirs.

3.5. Comparative Analysis of MMK and ITA Trend Methods

In this study, the observation of the rainfall and water quality data revealed distinct and significant insights into trends through the application of both traditional and statistical method MMK, and the more recent and graphical ITA method.
In terms of identifying the trend of rainfall, the trend detection from both MMK and ITA was mostly consistent. Both methods had identified notable trends where these were present, complementing each other’s findings and providing a comprehensive understanding of seasonal rainfall variations. For example, in the autumn season, Cooby Creek and Moyola rainfall stations showed a significant increasing trend in rainfall, with ITA indicators of 3.39 and 4.90 values respectively.
Similarly, in the analysis of water quality trends across the Cooby, Cressbrook and Perseverance reservoirs, both methods demonstrated the same alignment. In all reservoirs, PO43− showed a significant decreasing trend in all seasons and the average ITA indicator values were −4.40 in Cooby −1.27 in Cressbrook and −0.70 in Perseverance reservoir. The WQI showed a significant increasing trend in the Perseverance reservoir, which was also matched with the ITA indicator.
Compared with the MMK method, the ITA method offers several distinct advantages. It is more effective in identifying trends, particularly those subtle variations in data that conventional methods such as MMK may overlook. ITA’s visual and graphical representation of trends make it easier to identify extreme events such as significant rainfall leading to flooding or minimal rainfall resulting in drought [39]. In this study, the ITA graphs were generated for WQI only, in order to examine how WQI has changed over various seasons during the period of 22 years. These plots clearly show that WQI was on a decreasing trend in the first half (2001–2011) of the time series in Cooby and Cressbrook reservoirs. In this period, Toowoomba region experienced drought from 2000–2006 and extreme flooding from December 2010 to January 2011 [71]. Although the WQI showed a slight increasing trend in Perseverance reservoir, this was noticed in the second half of the data (2012–2022). ITA graphs allowed easy recognition of recurring patterns or periodic fluctuations of WQI. This visual approach can be particularly useful when identifying changes that might be masked by the MMK test.
Therefore, the MMK test’s statistical findings of decreasing or increasing trends in specific seasons were mathematically and visually confirmed by the ITA, which displayed corresponding shifts in the data. This dual approach not only validated the results, but also enhanced the interpretability of the data, making the observations accessible and understandable for broader audiences.

3.6. Wavelet Transform Coherence (WTC) Analysis

Following the trend analysis of rainfall and water quality parameters, the WTC method was applied in this study to identify the significant temporal relationship and co-movement of rainfall and WQI data series and to determine the prevailing frequencies or periods. The WTC plots of rainfall of twelve rainfall stations and WQI are presented in Figure 12, Figure 13 and Figure 14.

3.6.1. WTC of Rainfall and WQI (Cooby Reservoir)

The temporal relationship between rainfall and the WQI of Cooby reservoir is represented in Figure 12. The individual WTC plot of each station can be found in the link of Supplementary Files. The very high coherency (dark red) was visible in 8–16 periods from 2002–2022 between rainfall and WQI in the Cooby Creek rainfall stations and the Haden Post Office rainfall station (Figure 12a,f). This revealed notably strong coherence at extended frequencies, corresponding to cycles with periodicity in 8 to 16 and exceeding 16 within the time frame spanning from 2005 to 2018 in most of the stations. High correlation was significant in all stations during 4 to 8 periods from 2014 to 2020. On the other hand, low coherence is depicted in the Pechey Forestry Rainfall and WQI plot (Figure 12k), here only 2–4 periods from 2014–2018 exhibited high coherence. In the Doctor’s Creek WTC plot, a significant region from 2001 to 2008 with 16 periods was apparent, which showed the rainfall anomaly leading to WQI (straight down arrows in 5% significance level in Figure 12c,d). Based on Figure 12i, a significant region from 2004 to 2014 with an 8 to 16 period is evident. Most of the regions with a 5% significance level were scattered in all plots and are relatively small in Figure 12j,k.
Figure 12. WTC plot of Rainfall and WQI (Cooby). The coherency is shown by colour code, ranging from red (high coherency, value close to 1) to blue (low coherency, value close to 0). Statistically significant periodicity is displayed with black contours indicating a 5% significance level. The time scale ‘Year’ is shown in x axis and vertical axis (y axis) presents ‘Period’ across seasons in a year. The twelve rainfall stations are marked with lettering ranging (al). Individual plot of each station is available in the link of Supplementary Files.
Figure 12. WTC plot of Rainfall and WQI (Cooby). The coherency is shown by colour code, ranging from red (high coherency, value close to 1) to blue (low coherency, value close to 0). Statistically significant periodicity is displayed with black contours indicating a 5% significance level. The time scale ‘Year’ is shown in x axis and vertical axis (y axis) presents ‘Period’ across seasons in a year. The twelve rainfall stations are marked with lettering ranging (al). Individual plot of each station is available in the link of Supplementary Files.
Geosciences 14 00225 g012

3.6.2. WTC of Rainfall and WQI (Cressbrook Reservoir)

WTC plots of rainfall and WQI of Cressbrook reservoir revealed a mid-term overall behaviour of two variables exhibiting a mixed combination of coherency. In all plots of Figure 13, high coherency was displayed from 2015 to 2022 with the period 4 to 16. The arrows are pointing upward and in phase (right). The largest significant region was observed at the Oakey Aero rainfall station (Figure 13i) during the period from 4 to 16, exhibiting a phase status from 2018 to 2022 with a 5% significance level. However, from 2007 to 2012, a significant 8 to 16 period showed an anti-phase status and from 2019 to 2022, with the 8 to 16 period exhibiting a phase status (Figure 13a). Regions with a 5% significant level were scattered across all plots. The individual WTC plot of each station can be found in the link of Supplementary Files.
Figure 13. WTC plot of Rainfall and WQI (Cressbrook). The coherency is shown by colour code ranging from red (high coherency, value close to 1) to blue (low coherency, value close to 0). Statistically significant periodicity is displayed with black contours indicating a 5% significance level. The time scale ‘Year’ is shown in x axis and vertical axis (y axis) presents ‘Period’ across seasons in a year. The twelve rainfall stations are marked with lettering ranging from (al). Individual plot of each station is available in the link of Supplementary Files.
Figure 13. WTC plot of Rainfall and WQI (Cressbrook). The coherency is shown by colour code ranging from red (high coherency, value close to 1) to blue (low coherency, value close to 0). Statistically significant periodicity is displayed with black contours indicating a 5% significance level. The time scale ‘Year’ is shown in x axis and vertical axis (y axis) presents ‘Period’ across seasons in a year. The twelve rainfall stations are marked with lettering ranging from (al). Individual plot of each station is available in the link of Supplementary Files.
Geosciences 14 00225 g013

3.6.3. WTC of Rainfall and WQI (Perseverance Reservoir)

In Figure 14, the WTC is specified for rainfall of twelve stations (a to l) and the WQI for Perseverance reservoir. The individual WTC plot of each station can be found in the link of Supplementary Files. There was a large and significant area in the Haden Post Office plot (Figure 14f) from 2001 to 2016, with 8 to 16 periods that had combinations of upward and downward arrows and other regions with 4 to 8 periods from 2017 to 2020. The 5% significant zones are mostly visible in Oakey Aero (i), Glenaven (d), Tamba (l) and Cooby Creek (a). In the Oakey aero plot, significant coherence was visible in the periods from 2 to 16 between 2018 to 2022 and from 8 to 16 between 2001 to 2012. When comparing the overall coherency of rainfall and WQI across all three reservoirs, relatively low coherence was observed in the case of Perseverance reservoir.
Figure 14. WTC plot of Rainfall and WQI (Perseverance). The coherency is shown by colour code ranging from red (high coherency, value close to 1) to blue (low coherency, value close to 0). Statistically significant periodicity is displayed with black contours indicating a 5% significance level. The time scale ‘Year’ is shown in x axis and vertical axis (y axis) presents ‘Period’ across seasons in a year. The twelve rainfall stations are marked with lettering ranging from (al). Individual plot of each station is available in the link of Supplementary Files.
Figure 14. WTC plot of Rainfall and WQI (Perseverance). The coherency is shown by colour code ranging from red (high coherency, value close to 1) to blue (low coherency, value close to 0). Statistically significant periodicity is displayed with black contours indicating a 5% significance level. The time scale ‘Year’ is shown in x axis and vertical axis (y axis) presents ‘Period’ across seasons in a year. The twelve rainfall stations are marked with lettering ranging from (al). Individual plot of each station is available in the link of Supplementary Files.
Geosciences 14 00225 g014

4. Conclusions and Future Works

This study investigated the trend and correlations of rainfall and water quality parameters using statistical, graphical and mathematical methods. Previously, we estimated WQI based on five water quality parameters and this study used statistical trend analysis and applied wavelet transform techniques to demonstrate the pronounced seasonal variation of those parameters and their associated WQI along with rainfall patterns, with the intention of illuminating the complex interactions between climatic factors and water quality dynamics. A thorough exploration of spatial and temporal trends of rainfall and the potential time–frequency relationship between precipitation and water quality variations has been presented. The trends of rainfall across four seasons for the period of 22 years (2001–2022) were investigated using MMK and ITA methods. The MMK test was applied to check if the trend was significant, though the magnitude of the trend was calculated with Sen’s slope estimator. The ITA indicator was calculated to observe whether the trend was increasing or decreasing. Both methods (MMK and ITA) were also applied to examine the trend of water quality parameters. In the span of data of 22 years, the MMK test and ITA provided almost the same results in trend analysis. Moreover, wavelet transformation was used to see the correlation between seasonal rainfall and seasonal WQI. The major findings of the study can be summarised as follows:
  • A significant increasing trend of rainfall was observed in two rainfall stations (Cooby Creek and Moyola) in autumn and the combination of non-significant increasing and decreasing trends in the other three seasons. The rainfall during autumn increased by 5.43 mm/yr. The highest value of Sen’s slope was also observed in the autumn, (7.86 mm/yr) at Mount Kynoch.
  • In water quality trends, NH3-N showed increasing trends across four seasons in all three reservoirs. PO43− showed a significant decreasing trend in all three reservoirs in four seasons and was not affecting the trend of WQI. The declining trend of PO43− levels was a positive indicator of improved water quality in the reservoir. It indicated improved agricultural practice reducing phosphate runoff and better environmental regulations aimed at controlling phosphate discharge. The increasing trend of NH3-N, pH and turbidity affected the overall trend of WQI in Cooby and Cressbrook reservoirs. A slight increasing trend of WQI was observed in the Perseverance reservoir.
  • The MMK test provided statistical robustness and identified significant trends in rainfall and water quality data and ITA offered a detailed graphical representation highlighting subtle and nonlinear trends that MMK might have overlooked. In the case of water quality, WQI showed a decreasing trend in the first half of the data (2001–2012) in Cooby and Cressbrook reservoirs and slight increasing trend in the second half of the data in Perseverance reservoir. ITA’s visual output was particularly valuable in illustrating gradual changes and potential breakpoints, complementing the statistical validation provided by MMK.
  • In wavelet transform coherence (WTC) results, a notable correlation was found between seasonal rainfall and WQI. At Cooby Creek station, a very high coherence was noticed between rainfall and WQI in 8 to 16 periods ranging from 2002 to 2022. The regions with a 5% significance level were scattered in all plots. Among the twelve stations, the higher coherency was visible in Cooby Creek and at the Haden Post Office rainfall station with WQI of all reservoirs. However, comparatively low coherency was observed amidst rainfall and WQI of Perseverance reservoir. In terms of water quality trends and correlation analysis, Perseverance reservoir showed a slightly different scenario. The causes of this difference may include land use practices, catchment size, different outflow and inflow patterns and the variation in distribution and intensity of rainfall. This can be thoroughly investigated in future studies.
A few limitations of this study are worthwhile to mention, MMK relies on more than 30 years of data for optimal outcomes; however, this study utilised 22 years of temporal data, which may affect the overall results. Yang et al. [72] used 20 years of data to observe the changes in water quality and to find the effectiveness of water treatment. The previous studies [32,44] have suggested that both MMK and ITA methods provided reliable results; however, we only considered the rainfall effects due to the availability of data. Future research should explore the impact of runoff on water quality fluctuations.
Overall, this study provided a deeper understanding of the impacts of rainfall variability on water quality fluctuations, capturing both statistically significant changes and visually apparent patterns that provide a clear picture of environmental dynamics, thus informing decision making in water resource management and environmental conservation. Proactive conservation measures and community engagement can be collectively committed towards long-term sustainability and can ensure equitable access to water resources for present and future generations.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/geosciences14080225/s1.

Author Contributions

Conceptualisation: S.Z.F.; data collection and compilation: S.Z.F.; methods and analysis: S.Z.F.; writing—original draft preparation: S.Z.F.; writing—review and editing: D.R.P., M.J.A. and S.C.; supervision: D.R.P., M.J.A. and S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data should be available from first author after reasonable request.

Acknowledgments

The authors would like to thank the TRC for providing the water quality data and BOM for rainfall data. This research has been supported by the Graduate Research School, University of Southern Queensland and this is a part of the first author’s Ph.D. project entitled “An Integrated Spatial Decision Support Framework for Monitoring and Management of Surface Water Quality Influenced by Climate Change.”

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mathew, M.M.; Sreelash, K.; Mathew, M.; Arulbalaji, P.; Padmalal, D. Spatiotemporal variability of rainfall and its effect on hydrological regime in a tropical monsoon-dominated domain of Western Ghats, India. J. Hydrol. Reg. Stud. 2021, 36, 100861. [Google Scholar] [CrossRef]
  2. Ashcroft, L.; Karoly, D.J.; Dowdy, A.J. Historical extreme rainfall events in Southeastern Australia. Weather Clim. Extrem. 2019, 25, 100210. [Google Scholar] [CrossRef]
  3. Mohan, P.R.; Srinivas, C.; Yesubabu, V.; Baskaran, R.; Venkatraman, B. Simulation of a heavy rainfall event over Chennai in Southeast India using WRF: Sensitivity to microphysics parameterization. Atmos. Res. 2018, 210, 83–99. [Google Scholar] [CrossRef]
  4. Sun, Q.; Miao, C.; Duan, Q.; Wang, Y. Temperature and precipitation changes over the Loess Plateau between 1961 and 2011, based on high-density gauge observations. Glob. Planet. Chang. 2015, 132, 1–10. [Google Scholar] [CrossRef]
  5. Onyutha, C. Variability of seasonal and annual rainfall in the River Nile riparian countries and possible linkages to ocean–atmosphere interactions. Hydrol. Res. 2016, 47, 171–184. [Google Scholar] [CrossRef]
  6. Liu, S.; Qiu, Y.; Fu, R.; Liu, Y.; Suo, C. Identifying the water quality variation characteristics and their main driving factors from 2008 to 2020 in the Yellow River Basin, China. Environ. Sci. Pollut. Res. 2023, 30, 66753–66766. [Google Scholar] [CrossRef] [PubMed]
  7. Farzana, S.Z.; Paudyal, D.R.; Chadalavada, S.; Alam, M.J. Prediction of Water Quality in Reservoirs: A Comparative Assessment of Machine Learning and Deep Learning Approaches in the Case of Toowoomba, Queensland, Australia. Geosciences 2023, 13, 293. [Google Scholar] [CrossRef]
  8. Alam, M.J.; Dutta, D. Modelling of nutrient pollution dynamics in river basins: A review with a perspective of a distributed modelling approach. Geosciences 2021, 11, 369. [Google Scholar] [CrossRef]
  9. Shumway, R.H.; Stoffer, D.S. Time series regression and exploratory data analysis. In Time Series Analysis and Its Applications: With R Examples; Springer: Cham, Switzerland, 2017; pp. 45–74. [Google Scholar] [CrossRef]
  10. Almazroui, M.; Şen, Z. Trend analyses methodologies in hydro-meteorological records. Earth Syst. Environ. 2020, 4, 713–738. [Google Scholar] [CrossRef]
  11. Şen, Z. Innovative Trend Methodologies in Science and Engineering; Springer: Cham, Switzerland, 2017. [Google Scholar]
  12. Şen, Z. Innovative trend significance test and applications. Theor. Appl. Climatol. 2017, 127, 939–947. [Google Scholar] [CrossRef]
  13. Phuong, D.N.D.; Hai, L.M.; Dung, H.M.; Loi, N.K. Temporal trend possibilities of annual rainfall and standardized precipitation index in the Central Highlands, Vietnam. Earth Syst. Environ. 2022, 6, 69–85. [Google Scholar] [CrossRef]
  14. Imteaz, M.A.; Hossain, I. Climate change impacts on ‘seasonality index’and its potential implications on rainwater savings. Water Resour. Manag. 2023, 37, 2593–2606. [Google Scholar] [CrossRef]
  15. Berihun, M.L.; Tsunekawa, A.; Haregeweyn, N.; Tsubo, M.; Yasuda, H.; Fenta, A.A.; Dile, Y.T.; Bayabil, H.K.; Tilahun, S.A. Examining the past 120 years’ climate dynamics of Ethiopia. Theor. Appl. Climatol. 2023, 154, 535–566. [Google Scholar] [CrossRef]
  16. Fenta, A.A.; Yasuda, H.; Shimizu, K.; Haregeweyn, N.; Kawai, T.; Sultan, D.; Ebabu, K.; Belay, A.S. Spatial distribution and temporal trends of rainfall and erosivity in the Eastern Africa region. Hydrol. Process. 2017, 31, 4555–4567. [Google Scholar] [CrossRef]
  17. Naddafi, K.; Honari, H.; Ahmadi, M. Water quality trend analysis for the Karoon River in Iran. Environ. Monit. Assess. 2007, 134, 305–312. [Google Scholar] [CrossRef]
  18. Antonopoulos, V.Z.; Papamichail, D.M.; Mitsiou, K.A. Statistical and trend analysis of water quality and quantity data for the Strymon River in Greece. Hydrol. Earth Syst. Sci. 2001, 5, 679–692. [Google Scholar] [CrossRef]
  19. Xu, H.; Gao, Q.; Yuan, B. Analysis and identification of pollution sources of comprehensive river water quality: Evidence from two river basins in China. Ecol. Indic. 2022, 135, 108561. [Google Scholar] [CrossRef]
  20. Xiao, J.; Lv, G.; Chai, N.; Hu, J.; Jin, Z. Hydrochemistry and source apportionment of boron, sulfate, and nitrate in the Fen River, a typical loess covered area in the eastern Chinese Loess Plateau. Environ. Res. 2022, 206, 112570. [Google Scholar] [CrossRef]
  21. Delpla, I.; Jung, A.-V.; Baures, E.; Clement, M.; Thomas, O. Impacts of climate change on surface water quality in relation to drinking water production. Environ. Int. 2009, 35, 1225–1233. [Google Scholar] [CrossRef]
  22. Ponting, J.; Kelly, T.J.; Verhoef, A.; Watts, M.J.; Sizmur, T. The impact of increased flooding occurrence on the mobility of potentially toxic elements in floodplain soil—A review. Sci. Total Environ. 2021, 754, 142040. [Google Scholar] [CrossRef]
  23. Ockenden, M.C.; Hollaway, M.J.; Beven, K.J.; Collins, A.; Evans, R.; Falloon, P.; Forber, K.J.; Hiscock, K.; Kahana, R.; Macleod, C. Major agricultural changes required to mitigate phosphorus losses under climate change. Nat. Commun. 2017, 8, 161. [Google Scholar] [CrossRef]
  24. Mortazavi-Naeini, M.; Bussi, G.; Elliott, J.A.; Hall, J.W.; Whitehead, P.G. Assessment of risks to public water supply from low flows and harmful water quality in a changing climate. Water Resour. Res. 2019, 55, 10386–10404. [Google Scholar] [CrossRef]
  25. Watts, G.; Anderson, M. Water Climate Change Impacts Report Card 2016 Edition. Living with Environmental Change. 2016, pp. 9–10. Available online: https://www.ukri.org/wp-content/uploads/2021/12/091221-NERC-LWEC-WaterClimateChangeImpacts-ReportCard2016.pdf (accessed on 8 August 2024).
  26. Chen, J.; Gao, Y.; Qian, H.; Jia, H.; Zhang, Q. Insights into water sustainability from a grey water footprint perspective in an irrigated region of the Yellow River Basin. J. Clean. Prod. 2021, 316, 128329. [Google Scholar] [CrossRef]
  27. Quan, J.; Xu, Y.; Ma, T.; Wilson, J.P.; Zhao, N.; Ni, Y. Improving surface water quality of the Yellow River Basin due to anthropogenic changes. Sci. Total Environ. 2022, 836, 155607. [Google Scholar] [CrossRef] [PubMed]
  28. Becker, S.; Gemmer, M.; Jiang, T. Spatiotemporal analysis of precipitation trends in the Yangtze River catchment. Stoch. Environ. Res. Risk Assess. 2006, 20, 435–444. [Google Scholar] [CrossRef]
  29. Gemmer, M.; Fischer, T.; Jiang, T.; Su, B.; Liu, L.L. Trends in precipitation extremes in the Zhujiang River basin, South China. J. Clim. 2011, 24, 750–761. [Google Scholar] [CrossRef]
  30. Şen, Z. Moving trend analysis methodology for hydro-meteorology time series dynamic assessment. Water Resour. Manag. 2024, 1–15. [Google Scholar] [CrossRef]
  31. Sonali, P.; Kumar, D.N. Review of trend detection methods and their application to detect temperature changes in India. J. Hydrol. 2013, 476, 212–227. [Google Scholar] [CrossRef]
  32. Sun, W.; Mu, X.; Song, X.; Wu, D.; Cheng, A.; Qiu, B. Changes in extreme temperature and precipitation events in the Loess Plateau (China) during 1960–2013 under global warming. Atmos. Res. 2016, 168, 33–48. [Google Scholar] [CrossRef]
  33. Tabari, H.; Marofi, S.; Ahmadi, M. Long-term variations of water quality parameters in the Maroon River, Iran. Environ. Monit. Assess. 2011, 177, 273–287. [Google Scholar] [CrossRef]
  34. Wang, Y.; Xu, Y.; Lei, C.; Li, G.; Han, L.; Song, S.; Yang, L.; Deng, X. Spatio-temporal characteristics of precipitation and dryness/wetness in Yangtze River Delta, eastern China, during 1960–2012. Atmos. Res. 2016, 172, 196–205. [Google Scholar] [CrossRef]
  35. Westra, S.; Alexander, L.V.; Zwiers, F.W. Global increasing trends in annual maximum daily precipitation. J. Clim. 2013, 26, 3904–3918. [Google Scholar] [CrossRef]
  36. Zang, C.; Liu, J. Trend analysis for the flows of green and blue water in the Heihe River basin, northwestern China. J. Hydrol. 2013, 502, 27–36. [Google Scholar] [CrossRef]
  37. Mahmood, R.; Jia, S.; Zhu, W. Analysis of climate variability, trends, and prediction in the most active parts of the Lake Chad basin, Africa. Sci. Rep. 2019, 9, 6317. [Google Scholar] [CrossRef] [PubMed]
  38. Zhang, Q.; Liu, C.; Xu, C.-y.; Xu, Y.; Jiang, T. Observed trends of annual maximum water level and streamflow during past 130 years in the Yangtze River Basin, China. J. Hydrol. 2006, 324, 255–265. [Google Scholar] [CrossRef]
  39. Wang, Y.; Xu, Y.; Tabari, H.; Wang, J.; Wang, Q.; Song, S.; Hu, Z. Innovative trend analysis of annual and seasonal rainfall in the Yangtze River Delta, eastern China. Atmos. Res. 2020, 231, 104673. [Google Scholar] [CrossRef]
  40. Güçlü, Y.S. Multiple Şen-innovative trend analyses and partial Mann-Kendall test. J. Hydrol. 2018, 566, 685–704. [Google Scholar] [CrossRef]
  41. Dabanlı, İ.; Şen, Z.; Yeleğen, M.Ö.; Şişman, E.; Selek, B.; Güçlü, Y.S. Trend assessment by the innovative-Şen method. Water Resour. Manag. 2016, 30, 5193–5203. [Google Scholar] [CrossRef]
  42. Hamed, K.H.; Rao, A.R. A modified Mann-Kendall trend test for autocorrelated data. J. Hydrol. 1998, 204, 182–196. [Google Scholar] [CrossRef]
  43. Miles, K.; Joeland, B.; Bannon, K. Review of Regional Water Quality and Security 663, Volume 1—Review and Reform Strategy; Infrastucture Australia: Sydney, NSW, Australia, 2010; pp. 1–68. [Google Scholar]
  44. Wyrwoll, P.R.; Manero, A.; Taylor, K.S.; Rose, E.; Quentin Grafton, R. Measuring the gaps in drinking water quality and policy across regional and remote Australia. npj Clean Water 2022, 5, 32. [Google Scholar] [CrossRef]
  45. Willems, P.; Olsson, J. Impacts of Climate Change on Rainfall Extremes and Urban Drainage Systems; IWA PublishingAlliance House: London, UK, 2012. [Google Scholar]
  46. Queensland Government. Queensland Place Names. Available online: https://www.qld.gov.au/environment/land/title/place-names (accessed on 23 June 2024).
  47. Queensland Government. Business Queensland. Available online: https://www.business.qld.gov.au/running-business/support-assistance/mapping-data-imagery/maps/topographic-maps (accessed on 20 June 2024).
  48. Amos, M.; Baxter, G.; Finch, N.; Lisle, A.; Murray, P. I just want to count them! Considerations when choosing a deer population monitoring method. Wildl. Biol. 2014, 20, 362–370. [Google Scholar] [CrossRef]
  49. Bureau of Meteorology. Climate Statistics for Australian Locations. Available online: http://www.bom.gov.au/climate/averages/tables/cw_041103.shtml (accessed on 23 June 2024).
  50. Queensland Government. In Climate Change in the Eastern Downs Regions, Version 1, 2019; Department of Environment and Science, Queensland Government: Brisbane, QLD, Australia, 2019.
  51. Toowoomba Regional Council. Where Our Water Comes from. Available online: https://www.tr.qld.gov.au/environment-water-waste/water-sources-supply-storage-levels/dams-bores/13244-where-our-water-comes-from (accessed on 20 June 2024).
  52. Bureau of Meteorology. Climate Glossary. Available online: http://www.bom.gov.au/climate/glossary/seasons.shtml (accessed on 23 June 2024).
  53. Oliver, J.E. Monthly precipitation distribution: A comparative index. Prof. Geogr. 1980, 32, 300–309. [Google Scholar] [CrossRef]
  54. Siqueira, B.; Teixeira Nery, J. Spatial and temporal variability of precipitation concentration in northeastern Brazil. Investig. Geográficas 2021. [Google Scholar] [CrossRef]
  55. Abdi, H. Coefficient of variation. Encycl. Res. Des. 2010, 1, 169–171. [Google Scholar]
  56. Hare, W. Assessment of knowledge on impacts of climate change–contribution. Arctic 2003, 100, 25–35. [Google Scholar]
  57. Zaiontz, C. Mann Kendall Test. Available online: https://real-statistics.com/time-series-analysis/time-series-miscellaneous/mann-kendall-test/ (accessed on 23 June 2024).
  58. Gedefaw, M.; Yan, D.; Wang, H.; Qin, T.; Girma, A.; Abiyu, A.; Batsuren, D. Innovative trend analysis of annual and seasonal rainfall variability in Amhara regional state, Ethiopia. Atmosphere 2018, 9, 326. [Google Scholar] [CrossRef]
  59. Gao, F.; Wang, Y.; Chen, X.; Yang, W. Trend analysis of rainfall time series in Shanxi Province, Northern China (1957–2019). Water 2020, 12, 2335. [Google Scholar] [CrossRef]
  60. Tosunoglu, F. Trend analysis of daily maximum rainfall series in Çoruh Basin, Turkey. J. Inst. Sci. Technol. 2017, 7, 195–205. [Google Scholar] [CrossRef]
  61. Sen, P.K. Estimates of the regression coefficient based on Kendall’s tau. J. Am. Stat. Assoc. 1968, 63, 1379–1389. [Google Scholar] [CrossRef]
  62. Şen, Z. Innovative trend analysis methodology. J. Hydrol. Eng. 2012, 17, 1042–1046. [Google Scholar] [CrossRef]
  63. Pour, S.H.; Abd Wahab, A.K.; Shahid, S.; Ismail, Z.B. Changes in reference evapotranspiration and its driving factors in peninsular Malaysia. Atmos. Res. 2020, 246, 105096. [Google Scholar] [CrossRef]
  64. Şen, Z.; Şişman, E.; Dabanli, I. Innovative polygon trend analysis (IPTA) and applications. J. Hydrol. 2019, 575, 202–210. [Google Scholar] [CrossRef]
  65. Ebrahimi, A.; Rahimi, D.; Joghataei, M.; Movahedi, S. Correlation wavelet analysis for linkage between winter precipitation and three oceanic sources in Iran. Environ. Process. 2021, 8, 1027–1045. [Google Scholar] [CrossRef]
  66. Vacha, L.; Barunik, J. Co-movement of energy commodities revisited: Evidence from wavelet coherence analysis. Energy Econ. 2012, 34, 241–247. [Google Scholar] [CrossRef]
  67. Grinsted, A.; Moore, J.C.; Jevrejeva, S. Application of the cross wavelet transform and wavelet coherence to geophysical time series. Nonlinear Process. Geophys. 2004, 11, 561–566. [Google Scholar] [CrossRef]
  68. Torrence, C.; Compo, G.P. A practical guide to wavelet analysis. Bull. Am. Meteorol. Soc. 1998, 79, 61–78. [Google Scholar] [CrossRef]
  69. Torrence, C.; Webster, P.J. Interdecadal changes in the ENSO–monsoon system. J. Clim. 1999, 12, 2679–2690. [Google Scholar] [CrossRef]
  70. Das, J.; Mandal, T.; Rahman, A.S.; Saha, P. Spatio-temporal characterization of rainfall in Bangladesh: An innovative trend and discrete wavelet transformation approaches. Theor. Appl. Climatol. 2021, 143, 1557–1579. [Google Scholar] [CrossRef]
  71. Queensland Government. Queensland’s Extended Wet and Dry Periods; Ecosciences Precinct, Queensland Government: Brisbane, QLD, Australia, 2019.
  72. Yang, Q.; Chen, J.; Liu, Y.; Li, Y.; Zhang, H.; Zhang, J.; Sun, X.; Lu, M.; Ma, R.; Wu, Y. Analysis of changes in water quality and treatment effectiveness of seven major river basins in China from 2001 to 2020. Front. Environ. Sci. 2024, 12, 1340994. [Google Scholar] [CrossRef]
Figure 1. Study area map. The left map shows the area of rainfall stations, including the location of Cooby, Cressbrook and Perseverance reservoirs. The right map is the map of Australia, in which the blue highlighted area is the state of Queensland and the red marks indicate the approximate location of reservoirs.
Figure 1. Study area map. The left map shows the area of rainfall stations, including the location of Cooby, Cressbrook and Perseverance reservoirs. The right map is the map of Australia, in which the blue highlighted area is the state of Queensland and the red marks indicate the approximate location of reservoirs.
Geosciences 14 00225 g001
Figure 2. The schematic framework of the research method.
Figure 2. The schematic framework of the research method.
Geosciences 14 00225 g002
Figure 3. Spatial distribution of rainfall: (a) mean (RF), (b) coefficient of variation (CV) (%), and (c) precipitation concentration index (PCI).
Figure 3. Spatial distribution of rainfall: (a) mean (RF), (b) coefficient of variation (CV) (%), and (c) precipitation concentration index (PCI).
Geosciences 14 00225 g003
Figure 4. Spatial distribution of autumn rainfall trend parameters. Map (a) Z statistics, map (b) Sen’s Slope, and map (c) B (ITA).
Figure 4. Spatial distribution of autumn rainfall trend parameters. Map (a) Z statistics, map (b) Sen’s Slope, and map (c) B (ITA).
Geosciences 14 00225 g004
Figure 5. Spatial distribution of Winter rainfall trend parameters. Map (a) is Z statistics, Map (b) is Sen’s Slope, and Map (c) is B (ITA).
Figure 5. Spatial distribution of Winter rainfall trend parameters. Map (a) is Z statistics, Map (b) is Sen’s Slope, and Map (c) is B (ITA).
Geosciences 14 00225 g005
Figure 6. Spatial distribution of spring rainfall trend parameters. Map (a) Z statistics, map (b) Sen’s slope, and map (c) B (ITA).
Figure 6. Spatial distribution of spring rainfall trend parameters. Map (a) Z statistics, map (b) Sen’s slope, and map (c) B (ITA).
Geosciences 14 00225 g006
Figure 7. Spatial distribution of summer rainfall trend parameters. Map (a) Z statistics, map (b) Sen’s slope, and map (c) B (ITA).
Figure 7. Spatial distribution of summer rainfall trend parameters. Map (a) Z statistics, map (b) Sen’s slope, and map (c) B (ITA).
Geosciences 14 00225 g007
Figure 8. Plot of the innovative trend of WQI time series of Cooby reservoir (2001–2022). The blue dots represent the WQI in different seasons.
Figure 8. Plot of the innovative trend of WQI time series of Cooby reservoir (2001–2022). The blue dots represent the WQI in different seasons.
Geosciences 14 00225 g008
Figure 9. Plot of the innovative trend of the WQI time series of Cressbrook reservoir (2001–2022). The blue dots represent the WQI in different seasons.
Figure 9. Plot of the innovative trend of the WQI time series of Cressbrook reservoir (2001–2022). The blue dots represent the WQI in different seasons.
Geosciences 14 00225 g009
Figure 10. Plot of the innovative trend of WQI time series of Perseverance reservoir (2001–2022). The blue dots represent the WQI in different seasons.
Figure 10. Plot of the innovative trend of WQI time series of Perseverance reservoir (2001–2022). The blue dots represent the WQI in different seasons.
Geosciences 14 00225 g010
Figure 11. Modified Mann–Kendall test Z statistic plot, (top) Cooby, (middle) Cressbrook, (bottom) Perseverance.
Figure 11. Modified Mann–Kendall test Z statistic plot, (top) Cooby, (middle) Cressbrook, (bottom) Perseverance.
Geosciences 14 00225 g011
Table 2. Summary statistics of annual rainfall of 12 rainfall stations (2001–2022).
Table 2. Summary statistics of annual rainfall of 12 rainfall stations (2001–2022).
StationsLocationMSL (m)Annual Rainfall
Latitude (°)Longitude (°)Mean (mm)SD (mm)CV (%)PCI
Cooby Creek −27.3825151.924449764819329.815.9
Cressbrook Dam−27.2641152.195929576524732.315.8
Doctors Creek −27.2067151.846761266519629.515.2
Glenaven−27.1882151.963461270922631.915.7
Goombungee PO−27.3072151.850649760019232.016.1
Haden PO−27.2242151.883364059719833.217.7
Moyola−27.5233151.881955964522234.416.0
Mount Kynoch−27.5094151.954773989035039.316.2
Oakey Aero−27.4034151.741340657517129.715.9
Pechey Forestry−27.3042152.054266776326134.214.2
Perseverance Dam−27.2883152.123947077622128.515.5
Tamba−27.4722151.948164287831435.815.1
Table 3. Lag 1 autocorrelation of seasonal rainfall over the years (2001–2022) (yellow highlight indicates that an autocorrelation exists in the data series).
Table 3. Lag 1 autocorrelation of seasonal rainfall over the years (2001–2022) (yellow highlight indicates that an autocorrelation exists in the data series).
StationLag 1 ACF
Autumn
Lag 1 ACF
Winter
Lag 1 ACF
Spring
Lag 1 ACF
Summer
Cooby Creek0.04−0.06−0.21−0.14
Cressbrook Dam0.340.04−0.29−0.18
Doctors Creek0.40−0.02−0.25−0.02
Glenaven0.23−0.04−0.240.14
Goombungee PO0.26−0.18−0.300.09
Haden PO0.210.16−0.060.11
Moyola−0.02−0.01−0.010.05
Mount Kynoch0.180.04−0.010.13
Oakey Aero0.250.08−0.320.28
Pechey Forestry0.33−0.05−0.110.23
Perseverance Dam0.230.08−0.22−0.35
Tamba0.280.01−0.100.11
Table 5. Rainfall trend (winter).
Table 5. Rainfall trend (winter).
StationZ StatisticsMK Taup ValueSen’s Slope
(mm/yr)
ITDB (ITA)
Cooby Creek0.500.620.080.75Geosciences 14 00225 i0010.22
Cressbrook Dam0.360.060.720.53Geosciences 14 00225 i0010.18
Doctors Creek0.430.070.670.75Geosciences 14 00225 i0010.23
Glenaven−0.11−0.020.91−0.24Geosciences 14 00225 i002−0.48
Goombungee PO0.620.080.530.58Geosciences 14 00225 i0010.20
Haden PO−1.10−0.210.27−2.04Geosciences 14 00225 i002−1.82
Moyola−0.16−0.030.87−0.42Geosciences 14 00225 i002−0.76
Mount Kynoch−0.05−0.010.96−0.27Geosciences 14 00225 i002−0.51
Oakey Aero0.170.030.860.37Geosciences 14 00225 i0010.23
Pechey Forestry0.750.110.450.53Geosciences 14 00225 i0010.77
Perseverance Dam0.390.070.700.52Geosciences 14 00225 i0010.31
Tamba0.570.090.570.90Geosciences 14 00225 i0010.53
N.B: ITD = Innovative trend detection. The red arrow indicates an increasing trend, and the blue arrow shows a decreasing trend.
Table 6. Rainfall trend (spring).
Table 6. Rainfall trend (spring).
StationZ StatisticsMK Taup ValueSen’s Slope
(mm/yr)
ITDB (ITA)
Cooby Creek0.020.031.00−0.02Geosciences 14 00225 i002−0.86
Cressbrook Dam0.280.040.780.47Geosciences 14 00225 i0020.59
Doctors Creek0.880.110.381.69Geosciences 14 00225 i0010.41
Glenaven−0.07−0.010.95−0.20Geosciences 14 00225 i002−0.28
Goombungee PO0.720.080.472.33Geosciences 14 00225 i0011.45
Haden PO−0.24−0.040.81−0.76Geosciences 14 00225 i001−0.99
Moyola0.320.050.751.51Geosciences 14 00225 i0010.65
Mount Kynoch0.290.050.771.13Geosciences 14 00225 i0020.56
Oakey Aero0.220.030.830.57Geosciences 14 00225 i0010.60
Pechey Forestry0.180.030.860.74Geosciences 14 00225 i0010.75
Perseverance Dam0.070.010.950.36Geosciences 14 00225 i0010.92
Tamba−0.06−0.010.95−0.91Geosciences 14 00225 i002−0.82
N.B: ITD = Innovative trend detection. The red arrow indicates an increasing trend, and the blue arrow shows a decreasing trend.
Table 7. Rainfall trend (summer).
Table 7. Rainfall trend (summer).
StationZ StatisticsMK Taup ValueSen’s Slope
(mm/yr)
ITDB (ITA)
Cooby Creek−0.24−0.040.81−0.89Geosciences 14 00225 i002−0.64
Cressbrook Dam0.620.080.532.80Geosciences 14 00225 i0010.20
Doctors Creek0.430.070.672.00Geosciences 14 00225 i0011.48
Glenaven0.230.040.820.48Geosciences 14 00225 i0010.77
Goombungee PO−0.10−0.020.92−0.26Geosciences 14 00225 i002−0.53
Haden PO−1.46−0.210.14−5.48Geosciences 14 00225 i002−2.11
Moyola0.100.020.920.14Geosciences 14 00225 i0010.49
Mount Kynoch0.470.080.643.31Geosciences 14 00225 i0012.25
Oakey Aero−0.42−0.050.68−1.40Geosciences 14 00225 i002−0.77
Pechey Forestry−0.30−0.060.77−1.71Geosciences 14 00225 i002−0.75
Perseverance Dam0.530.060.600.87Geosciences 14 00225 i0010.24
Tamba0.520.090.601.73Geosciences 14 00225 i0010.10
N.B: ITD = Innovative trend detection. The red arrow indicates an increasing trend, and the blue arrow shows a decreasing trend.
Table 8. Details of water quality parameters trend (Cooby reservoir).
Table 8. Details of water quality parameters trend (Cooby reservoir).
ParameterSeasonZ StatisticMK Taup ValueSen’s SlopeITDB (ITA)
NH3-NAutumn0.330.120.740.01Geosciences 14 00225 i0010.07
Winter0.350.130.730.01Geosciences 14 00225 i0012.37
Spring0.500.110.620.01Geosciences 14 00225 i0010.27
Summer0.080.030.940.01Geosciences 14 00225 i0010.30
pHAutumn−0.25−0.040.810.00Geosciences 14 00225 i002−0.06
Winter0.840.110.400.01Geosciences 14 00225 i0010.00
Spring0.610.090.540.00Geosciences 14 00225 i0010.08
Summer−1.33−0.190.18−0.01Geosciences 14 00225 i002−0.10
PO43−Autumn−3.19−0.390.000.00Geosciences 14 00225 i003−4.10
Winter−2.09−0.380.040.00Geosciences 14 00225 i003−3.83
Spring−4.55−0.640.000.00Geosciences 14 00225 i003−3.95
Summer−4.22−0.630.000.00Geosciences 14 00225 i003−4.30
TDSAutumn0.060.030.952.64Geosciences 14 00225 i0011.89
Winter0.000.001.000.73Geosciences 14 00225 i0011.93
Spring0.040.020.972.55Geosciences 14 00225 i0011.81
Summer0.130.050.902.49Geosciences 14 00225 i0011.81
TurbidityAutumn1.240.180.210.04Geosciences 14 00225 i0010.09
Winter0.870.150.380.03Geosciences 14 00225 i0010.03
Spring0.420.140.680.02Geosciences 14 00225 i0010.95
Summer0.200.040.840.01Geosciences 14 00225 i0020.09
WQIAutumn−0.06−0.030.95−0.15Geosciences 14 00225 i002−1.83
Winter−0.04−0.020.97−0.05Geosciences 14 00225 i001−1.82
Spring−0.06−0.030.95−0.05Geosciences 14 00225 i002−1.60
Summer−0.04−0.020.97−0.03Geosciences 14 00225 i002−1.73
N.B: ITD = Innovative trend detection. The red arrow indicates the increasing trend, and the blue arrow shows the decreasing trend. A significant trend is highlighted with the yellow colour (p value < 0.05).
Table 9. Details of water quality parameters trend (Cressbrook reservoir).
Table 9. Details of water quality parameters trend (Cressbrook reservoir).
ParameterSeasonZ StatisticMK Taup ValueSen’s SlopeITDB (ITA)
NH3-NAutumn0.950.160.340.01Geosciences 14 00225 i0010.06
Winter1.730.270.080.01Geosciences 14 00225 i0010.05
Spring0.580.950.560.01Geosciences 14 00225 i0010.62
Summer0.400.110.690.01Geosciences 14 00225 i0010.06
pHAutumn0.220.040.820.00Geosciences 14 00225 i0010.01
Winter0.870.190.380.01Geosciences 14 00225 i0010.03
Spring−0.25−0.040.810.00Geosciences 14 00225 i002−0.09
Summer0.190.040.850.00Geosciences 14 00225 i0010.03
PO43−Autumn−4.88−0.510.000.00Geosciences 14 00225 i002−1.06
Winter−4.59−0.520.000.00Geosciences 14 00225 i002−1.52
Spring−3.90−0.580.000.00Geosciences 14 00225 i002−1.22
Summer−3.55−0.610.000.00Geosciences 14 00225 i002−1.31
TDSAutumn−0.24−0.080.81−0.88Geosciences 14 00225 i002−0.98
Winter0.000.001.00−0.02Geosciences 14 00225 i002−0.08
Spring−0.09−0.040.93−0.46Geosciences 14 00225 i002−0.84
Summer−0.16−0.060.88−0.48Geosciences 14 00225 i002−0.59
TurbidityAutumn2.080.330.040.05Geosciences 14 00225 i0010.69
Winter1.290.190.200.03Geosciences 14 00225 i0010.33
Spring−0.30−0.080.76−0.01Geosciences 14 00225 i002−0.08
Summer0.880.200.380.02Geosciences 14 00225 i002−0.68
WQIAutumn−0.05−0.020.96−0.01Geosciences 14 00225 i002−0.01
Winter−0.13−0.050.89−0.01Geosciences 14 00225 i002−0.71
Spring−0.10−0.040.92−0.01Geosciences 14 00225 i002−0.56
Summer−0.12−0.040.91−0.02Geosciences 14 00225 i002−0.72
N.B: ITD = Innovative trend detection. The red arrow indicates the increasing trend and the blue arrow shows the decreasing trend. A significant trend is highlighted with the yellow colour (p value < 0.05).
Table 10. Details of water quality parameters trend (Perseverance reservoir).
Table 10. Details of water quality parameters trend (Perseverance reservoir).
ParameterSeasonZ StatisticMK Taup ValueSen’s SlopeITDB (ITA)
NH3-NAutumn0.870.220.390.01Geosciences 14 00225 i0010.53
Winter0.380.160.710.00Geosciences 14 00225 i0010.99
Spring0.400.140.690.00Geosciences 14 00225 i0010.05
Summer0.170.080.860.01Geosciences 14 00225 i0010.05
pHAutumn0.100.030.920.00Geosciences 14 00225 i0010.01
Winter0.410.110.680.01Geosciences 14 00225 i0010.06
Spring0.270.050.790.01Geosciences 14 00225 i0010.08
Summer−0.06−0.020.950.00Geosciences 14 00225 i002−0.01
PO43−Autumn−3.76−0.540.00−0.01Geosciences 14 00225 i002−1.00
Winter−2.33−0.400.020.00Geosciences 14 00225 i002−1.25
Spring−2.35−0.510.020.00Geosciences 14 00225 i002−0.04
Summer−2.72−0.630.010.00Geosciences 14 00225 i002−0.52
TDSAutumn1.610.380.111.43Geosciences 14 00225 i0011.11
Winter1.970.380.051.34Geosciences 14 00225 i0011.65
Spring2.050.440.041.47Geosciences 14 00225 i0012.01
Summer1.770.490.081.60Geosciences 14 00225 i0011.63
TurbidityAutumn1.820.400.070.16Geosciences 14 00225 i0010.63
Winter2.180.330.290.14Geosciences 14 00225 i0010.89
Spring2.140.420.030.16Geosciences 14 00225 i0010.99
Summer1.640.420.100.17Geosciences 14 00225 i0010.93
WQIAutumn2.050.450.040.06Geosciences 14 00225 i0010.40
Winter1.870.460.060.05Geosciences 14 00225 i0010.41
Spring1.690.420.090.05Geosciences 14 00225 i0010.33
Summer2.130.570.030.07Geosciences 14 00225 i0010.44
N.B: ITD = Innovative trend detection. The red arrow indicates the increasing trend, and the blue arrow shows the decreasing trend. A significant trend is highlighted with the yellow colour (p value < 0.05).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Farzana, S.Z.; Paudyal, D.R.; Chadalavada, S.; Alam, M.J. Spatiotemporal Variability Analysis of Rainfall and Water Quality: Insights from Trend Analysis and Wavelet Coherence Approach. Geosciences 2024, 14, 225. https://doi.org/10.3390/geosciences14080225

AMA Style

Farzana SZ, Paudyal DR, Chadalavada S, Alam MJ. Spatiotemporal Variability Analysis of Rainfall and Water Quality: Insights from Trend Analysis and Wavelet Coherence Approach. Geosciences. 2024; 14(8):225. https://doi.org/10.3390/geosciences14080225

Chicago/Turabian Style

Farzana, Syeda Zehan, Dev Raj Paudyal, Sreeni Chadalavada, and Md Jahangir Alam. 2024. "Spatiotemporal Variability Analysis of Rainfall and Water Quality: Insights from Trend Analysis and Wavelet Coherence Approach" Geosciences 14, no. 8: 225. https://doi.org/10.3390/geosciences14080225

APA Style

Farzana, S. Z., Paudyal, D. R., Chadalavada, S., & Alam, M. J. (2024). Spatiotemporal Variability Analysis of Rainfall and Water Quality: Insights from Trend Analysis and Wavelet Coherence Approach. Geosciences, 14(8), 225. https://doi.org/10.3390/geosciences14080225

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop