The Association between Meteorological Drought and the State of the Groundwater Level in Bursa, Turkey

: This study addressed the intricate interplay between meteorological droughts and ground-water level ﬂuctuations in the vicinity of Mount Uludag in Bursa, Turkey. To achieve this, an exhaustive analysis encompassing monthly precipitation records and groundwater level data sourced from three meteorological stations and eight groundwater observation points spanning the period from 2007 to 2018 was performed. Subsequently, this study employed the Standard Precipitation Index (SPI) and Standard Groundwater Level (SGL) metrics, meticulously calculating the temporal extents of drought events for each respective time series. Following this, a judicious application of both the Thiessen and Support Vector Machine (SVM) methodologies was undertaken to ascertain the optimal groundwater observation wells and their corresponding SGL durations, aligning them with SPI durations tied to the selected meteorological stations. The SVM technique, in particular, excelled in the identiﬁcation of the most pertinent observation wells. Additionally, the Elman Neural Network (ENN) and its optimized version through the Fireﬂy Algorithm (ENN-FA), demonstrated their prowess in accurately predicting SPI durations based on SGL durations. The results were favorable, as evidenced by the commendable performance metrics of the Normalized Root Mean Square Error (NRMSE), the Nash–Sutcliffe Efﬁciency (NSE), the product of the coefﬁcient of determination and the slope of the regression line (bR 2 ), and the Kling–Gupta Efﬁciency (KGE). Consequently, the favorable simulation results were construed as evidence supporting the presence of a discernible association between SGL and the duration of the SPI. As we substantiate the concordance between the temporal extent of meteorological droughts and the perturbations in groundwater levels, this unmistakably underscores the fact that the historical ﬂuctuations in groundwater levels within the region were predominantly attributable to climatic inﬂuences, rather than being instigated by anthropogenic activities. Nevertheless, it is imperative to underscore that this revelation should not be misconstrued as an endorsement of future heedless exploitation of groundwater resources.


Introduction
Extreme weather events are well-known threats around the globe [1].In particular, drought represents a multifaceted and intricate phenomenon, predominantly contingent upon the temporal and spatial availability of water resources [2,3].It has the potential to strain agro-food and socio-economic stability to the point of collapse [4,5].Its genesis lies in the diminishing or outright absence of precipitation, constituting meteorological drought, which subsequently evolves into more intricate hydro-environmental manifestations [6].These include the depletion of soil moisture, culminating in agricultural drought, as well as the reduction in streamflow rates, representing hydrological drought.Moreover, the dearth of accessible water resources for essential human activities constitutes socio-economic drought [7].Throughout history, the climate has undergone significant transformations, primarily driven by natural forces.However, the contemporary climate crisis is inexorably linked to rapid industrialization and the pervasive emission of greenhouse gases.This climate metamorphosis engenders fluctuations in water resource availability, instigating extreme climatic events such as droughts and floods.Nonetheless, the classification of drought hinges on the nature of the available water resources within a specific locale, dictating the norms of climatic conditions in terms of both temporal and spatial dimensions and what can be accepted as a drought for a specific region.
Groundwater, a vital component of the Earth's freshwater resources, constitutes approximately 30 percent of the available freshwater worldwide [8].In regions where surface water resources are either scarce or unreliable, groundwater emerges as the optimal alternative to meet various water demands.The controlled release of aquifer-stored water through natural or engineered means, including wells, drainage systems, aqueducts, springs, and pumping, underscores the importance of groundwater as a resource.Paradoxically, the abundant availability of groundwater in aquifers can also be a catalyst for droughts, often intertwined with meteorological factors such as precipitation and snow cover, hydrological variables like streamflow, and human-induced activities.It is worth noting that groundwater levels are intricately linked to the recharge processes facilitated by the infiltration of precipitation and snowmelt, particularly in high-altitude areas [9].Consequently, the nexus between meteorological drought and groundwater drought (excluding fossil or deep reservoirs), exemplified by fluctuations in groundwater levels, plays a pivotal role in shaping local groundwater conditions within a given region [10].A profound comprehension of the intricate interactions between meteorological phenomena and groundwater dynamics (aquifers) is of paramount significance, particularly in regions characterized by hydrological variability [11][12][13].For instance, a comprehensive investigation conducted by Bloomfield and Marchant [14] addressed the relationship between the Standard Precipitation Index (SPI) and the Groundwater Level Index, using a dataset spanning 103 years.Their findings underscored the dependence of Groundwater Level Index autocorrelation on the dominating aquifer flow and storage characteristics.Similarly, Liu et al. [15] conducted a study on the state of groundwater drought and its response to meteorological drought in China.By employing the standardized groundwater level index in conjunction with SPI, they conducted a cross-correlation analysis and identified trends in the respective indices.Their conclusion revealed that the duration and magnitude of groundwater droughts were notably protracted and substantial.
In recent years, substantial strides have been made in predictive modeling techniques harnessing the capabilities of machine learning methodologies [16,17].In particular, the Elman Neural Network (ENN), enhanced through nature-inspired optimization algorithms, has garnered considerable attention due to its adeptness in deciphering intricate relationships within input and output datasets [18,19].The ENN, a recurrent neural network architecture, excels in modeling sequential data by accounting for temporal dependencies.Augmenting its capabilities, the Firefly Algorithm (FA) offers an optimization approach inspired by natural phenomena, resulting in refined model architectures and heightened predictive accuracy [20,21].Similarly, the Support Vector Machine (SVM), a versatile machine learning algorithm, is renowned for its prowess in feature selection and classification tasks [22,23].Within the context of our study, SVM adeptly identifies the most influential meteorological stations, streamlining our analysis and yielding invaluable insights.
The interplay between groundwater abundance and its correlation with meteorological phenomena in the Bursa region of Turkey presents a distinct case that has not received the depth of investigation it warrants.While prior research has delved into aspects such as thermal waters and local aquifer modeling [24][25][26], comprehensive analysis of this relationship remains underexplored.Existing studies have primarily focused on the biological and chemical attributes of water content in the region, leaving a gap in our understanding of the intricate dynamics governing the connection between meteorological drought and groundwater level fluctuations.In general, the records pertaining to regional aquifer groundwater and the hydrogeology of the area have not undergone sufficient investigation.Historically, groundwater has not been the primary water source in this region.However, in light of recent droughts, there is significant potential for groundwater to assume a leading role as the primary water source.However, the traditional approaches toward probing the association between meteorological drought and groundwater level changes have predominantly relied on rudimentary correlations, overlooking the nuanced temporal and spatial intricacies that underlie these associations.
In response to this knowledge gap, we present an exhaustive study employing advanced methodologies to unravel the intricate connections between the SPI and the so-called Standard Groundwater Level (SGL) durations within the vicinity of Mount Uludag.This research also endeavors to elucidate the link between disturbances in the local groundwater table and precipitation deficits within the broader Bursa district.This is primarily imperative in comprehending and strategizing for the potential utilization of subterranean aquifers as a dynamic water resource and fortification against climate perturbations within the locale.To achieve this, we apply the following steps.
(i) The SGL and SPI for groundwater level and precipitation records are determined.(ii) The duration curves of SGL and SPI are scrutinized to ascertain the relevance of their respective durations.(iii) The SVM and Thiessen methods are used alternately in the selection of the most relevant stations.(iv) The presence of an association between SGL and SPI in the selected stations previously determined either by means of SVM or the Thiessen was tested using cross-correlation analysis and the estimation of SPI durations with the help of SGL durations using ENN empowered by FA.
As a result, by confirming the presence of an association between groundwater disturbances and meteorological droughts, this study is expected to provide grounds for future studies and avoiding encroachments on groundwater sources.

Study Area
The study area is situated in the northwestern region of Turkey and boasts a notably extensive coastline along the Marmara Sea.Geographically, it spans between the latitudes of 39 • 34 and 40 • 38 N and the longitudes of 28 • 05 and 29 • 60 E, while its topographical variation ranges from sea level up to an elevation of 2543 m above sea level (mASL) at the summit of Mount Uludag.In accordance with the Koppen-Geiger classification system, as elucidated by Kottek et al. [27], the climate within this region falls into two distinct categories: either a warm temperate climate characterized by dry summers and cold-arid temperatures (Csa), or a warm temperate climate marked by dry summers and warm summer temperatures (Csb).The primary river systems coursing through this region include the Mustafakemalpasa, Nilufer, Goksu, Kocadere, Karadere, and Aksudere, all of which play a pivotal role in shaping the riverine hydrology and water cycle within the area.For a comprehensive overview, Figure 1 illustrates the study area, climate regions based on the Koppen-Geiger classification map (attainable in: https://koeppen-geiger.vuwien.ac.at/present.htm(accessed on 29 October 2023)), and the locations of the selected meteorological stations and observation wells.It also showcases the Thiessen polygons utilized in the meticulous selection process of the optimal meteorological station, a matter we will delve into later.
Consequently, the analysis employed the monthly time series of precipitation data gathered from three meteorological sites (Osmangazi, Uludag, and Keles) in conjunction with the groundwater level time series data recorded at eight observation wells (Adakoy, Aticilar, Cayirkoy, Delicay, Kursunlu, Narlidere, Yeniceabat, and Yenice) throughout the water year, spanning from 2007 to 2018.For further details regarding the precise locations of the meteorological stations and wells, as well as the fundamental statistical properties of the monthly time series utilized in this analysis, refer to Table 1.Notably, all of the meteorological stations are strategically positioned within the lofty Uludag mountain range at high altitudes, while the selected observation wells are situated on hillsides at considerably lower elevations.This positioning was implemented to naturally establish a hydraulic gradient from the meteorological stations towards the groundwater wells on the hillsides.Bearing this key fact in mind, our analysis endeavored to examine the degree of association between precipitation and groundwater time series by investigating both the simultaneous and lagged behaviors of these two crucial variables.Consequently, the analysis employed the monthly time series of precipitation data gathered from three meteorological sites (Osmangazi, Uludag, and Keles) in conjunction with the groundwater level time series data recorded at eight observation wells (Adakoy, Aticilar, Cayirkoy, Delicay, Kursunlu, Narlidere, Yeniceabat, and Yenice) throughout the water year, spanning from 2007 to 2018.For further details regarding the precise locations of the meteorological stations and wells, as well as the fundamental statistical properties of the monthly time series utilized in this analysis, refer to Table 1.As depicted in Figure 1 and through the meticulous application of Thiessen polygons, we were able to establish associations between the data records obtained from Uludag, Os-mangazi, and Keles meteorological stations and from specific locations, at Adakoy, Delicay, Narli Dere, and Yenice; Aticilar, Cayirkoy, and Yeniceabat; and Kursunlu, respectively.This crucial criterion subsequently served as a pivotal determinant in selecting the most suitable input stations for the development of the models, a subject that will be expounded upon in detail later in this analysis.
Table 1 also provides a comprehensive breakdown of the monthly averaged precipitation time series, coupled with the mean groundwater level time series derived from observation wells (from unconfined aquifers) dispersed across our study region.Notably, the station registering the highest mean precipitation rate was Uludag, while the lowest values were associated with the Osmangazi station.Intriguingly, despite Osmangazi's lower average, it recorded maximum precipitation levels surpassing those observed at the Keles station.Turning our attention to the groundwater data, the highest mean groundwater level elevations were documented at the Kursunlu and Yenice stations, while the lowest elevation was recorded at Yeniceabat.
In light of this information, the time series data spanning from January 2007 to November 2018, encompassing a total of 143 months, served as a foundational framework for evaluating the temporal congruence between groundwater and precipitation records.

Methods
The primary objective of this study was to ascertain the relationship between precipitation patterns at higher elevations and the fluctuations in groundwater levels along hillside regions.This investigation centered on the concurrent dynamics of meteorological drought and perturbations in groundwater levels.The methodology outlined in Figure 2 served as the cornerstone of this research endeavor.

Standard Precipitation Index (SPI)
The Standard Precipitation Index (SPI), a renowned drought index, establishes the water abundance within a particular location by scrutinizing the precipitation deficit or surplus.This index is formulated based on the probability of outcomes derived from prefitted probability distributions such as Gamma, Gumbel, among others [28].These distributions are subsequently employed for standardization, facilitating the quantification of water deficit and surplus intensities.
Within the context of this study, we focused on monthly SPI rates, also referred to as SPI-1 values, in our analysis.However, it is worth noting that the incorporation of longterm deficit/surplus responses can also be achieved using the moving average operator.The preference for SPI-1 over long-term SPI values in this study is attributed to two key factors: the limited availability of data records pertaining to groundwater levels and the inherent complexity associated with short-term drought prediction, which took precedence.
Given the well-established nature of the methodology and equations pertaining to SPI determination, we abstain from reiterating the same details in this section.Interested In the initial phase, we computed two pivotal indices, the SPI and the SGL.Subsequently, we arranged the time series data for SGL and SPI in ascending order to construct duration curves, thereby determining the corresponding percentiles for both SGL and SPI.Employing a reverse analogical approach, we assessed whether the duration of the SGL index correlates with that of the meteorological drought (i.e., SPI), aiming to confirm or refute any link between groundwater table disturbances and meteorological droughts.To minimize the influence of potential confounding factors, such as neighboring aquifers and surface water, which may impact groundwater conditions on the hillside, we employed SGL as an input parameter and SPI as the output in our simulations.
The identification of the most representative observation wells alternated between employing Thiessen polygons (as depicted in Figure 1) and the SVM as input variables for our simulation models.Subsequently, the scenarios derived from either Thiessen polygons or SVM served as the basis for forecasting monthly SPI values at meteorological sites.Initially, SPI and SGL functioned as the output and input variables, respectively, within our models.The selection of training and testing data was carried out through a randomized process, aligning with the concurrent duration of SPI and SGL.The ENN was subsequently employed, with the results fine-tuned using the FA, ultimately leading to the introduction of the ENN-FA approach.This approach, along with Thiessen (referred to as Thiessen-ENN-FA) and SVM (known as SVM-ENN-FA) methodologies, will be elaborated upon in greater detail later in this study.

Standard Precipitation Index (SPI)
The Standard Precipitation Index (SPI), a renowned drought index, establishes the water abundance within a particular location by scrutinizing the precipitation deficit or surplus.This index is formulated based on the probability of outcomes derived from pre-fitted probability distributions such as Gamma, Gumbel, among others [28].These distributions are subsequently employed for standardization, facilitating the quantification of water deficit and surplus intensities.
Within the context of this study, we focused on monthly SPI rates, also referred to as SPI-1 values, in our analysis.However, it is worth noting that the incorporation of longterm deficit/surplus responses can also be achieved using the moving average operator.The preference for SPI-1 over long-term SPI values in this study is attributed to two key factors: the limited availability of data records pertaining to groundwater levels and the inherent complexity associated with short-term drought prediction, which took precedence.
Given the well-established nature of the methodology and equations pertaining to SPI determination, we abstain from reiterating the same details in this section.Interested readers seeking a comprehensive understanding of SPI calculation may refer to the studies conducted by McKee et al. [28], Aksoy [29], and Vaheddoost and Safari [30] for in-depth insights and further applications of the SPI.

Standard Groundwater Level (SGL)
Diverging from prior approaches presented in earlier studies such as those by Bloomfield and Marchant [14] and Liu et al. [15], which introduced the groundwater drought index by employing probability distributions like SPI, this study proposes a simpler standardization method similar to the Z-Score Index (ZSI) for the determination of SGL, as elucidated in Vaheddoost and Safari [30].The ZSI primarily centers on identifying authentic anomalies within the dataset.
In this context, classical standardization, using the formula z = (x -µ) σ , was implemented and served either as a drought indicator or as a means to highlight anomalies in the temporal correspondence of data.In this equation, the normalized values (z) are derived by subtracting the long-term mean groundwater level (µ) from the monthly groundwater level records (x), thus revealing groundwater level anomalies.These anomalies are then divided by the sample standard deviation of the groundwater level time series (σ).
The acquired SGL time series was subsequently employed in conjunction with the SPI time series in a multitude of analyses aimed at probing the depth of the relationship between SPI and SGL.

Association between SPI and SGL
The monthly time series, coupled with the corresponding duration curves for both SGL and SPI, were instrumental in pursuing this objective.It is postulated that the temporal and spatial congruence of disruptions in the SGL and SPI time series can be regarded as indicative of the propensity for an association between SGL and SPI events.In essence, this implies that if the duration of SPI can be elucidated through SGL duration, and the selection of the study period is conducted on a randomized basis, one can infer that meteorological factors constitute the primary drivers behind groundwater disturbances at the chosen sites.
Hence, in order to validate this hypothesis concerning the relationship between SGL and SPI, the durations of randomly selected SPI events were predicted using the same ranking criteria associated with the corresponding SGL events.As delineated earlier, the Thiessen and SVM methods came into play in the determination of the most optimal scenario-stations and the identification of the most suitable inputs for the models.The SVM stands as a formidable machine learning algorithm renowned for its effectiveness across a spectrum of tasks, encompassing classification, regression, and feature selection, as evidenced by Gunn [31].In our study, SVM assumes the role of a feature selection algorithm, its purpose being the identification of the most pertinent meteorological stations for dissecting the interrelation between SPI and SGL durations.This approach enables a meticulous focus on the meteorological stations wielding the greatest influence over predictive accuracy.At the heart of SVM lies the pivotal concept of locating a hyperplane that optimally segregates data points belonging to distinct classes, a notion substantiated by the works of Sugumaran et al. [32] and Cervantes et al. [33].In the realm of feature selection, SVM undertakes the task of pinpointing the subset of features-in this instance, meteorological stations-that maximizes the demarcation between classes.Within our study, these classes align with SPI and SGL durations.
The SVM process encompasses several crucial phases: (i) Data preparation: each meteorological station's SPI and SGL durations are represented as feature vectors, serving as the foundation for constructing a training dataset for SVM.(ii) Kernel functions: SVM leverages kernel functions to map input data into a higher-dimensional space.This mapping enhances data separability, simplifying the quest for a hyperplane that efficiently distinguishes classes.In our investigation, we explored six distinct kernel functions, including dot, radial, polynomial, neural, ANOVA (analysis of variance), and Epachnenikov [34,35].(iii) Margin maximization: SVM's primary aim is the discovery of a hyperplane that maximizes the margin, defined as the distance between the hyperplane and the nearest data points from each class.This strategy fosters superior generalization to unseen data.(iv) Support vectors: the data points situated closest to the hyperplane assume the designation of support vectors, playing a pivotal role in determining the hyperplane's position and, consequently, the decision boundaries.(v) Soft margin and regularization: SVM accommodates scenarios where data lacks linear separability by introducing a soft margin, permitting a degree of misclassification.A regularization parameter governs the equilibrium between margin maximization and misclassification minimization.
In our case, SVM undertook the task of ranking the significance of meteorological stations in forecasting SPI and SGL durations.The various kernel functions were systematically applied to the training dataset, resulting in feature rankings that facilitated the selection of the most influential stations.Through the selection of meteorological stations bearing the highest rankings, we aspired to elevate the precision of our predictive model and unearth insights into the primary drivers underpinning the connection between meteorological drought and groundwater level fluctuations.

Application of Elman Neural Network (ENN) Coupled by Firefly Algorithm (FA)
The ENN represents a recurrent neural network architecture with a track record of applications in modeling soil, water, and meteorological time series, as demonstrated in the works of Zhang et al. [36] and Li et al. [37].ENN's distinct advantage lies in its capacity to capture temporal dependencies within data, rendering it particularly adept for tasks involving sequential patterns and relationships, as underscored by Lu et al. [38].In this study, we harnessed the ENN as a predictive tool to establish the nexus between SPI and SGL durations, thereby delving into the correlation between meteorological drought and fluctuations in groundwater levels.
The ENN is structured around interconnected layers of processing units, commonly known as neurons, as expounded upon by Li et al. [39].These neurons are organized into input, hidden, and output layers.Notably, the ENN incorporates a context layer, introducing a feedback loop that allows the network to retain information from previous time steps [40,41].This contextual information plays a pivotal role in modeling temporal relationships within the data.The ENN comprises four core layers: the input layer, hidden layer, undertake layer, and output layer [42].The input layer receives SGL data, which are utilized to predict SPI durations.The hidden layer, distinguished by its recurrent connections, captures temporal dependencies and patterns within the data, enabling the ENN to adeptly model sequential information.undertake layer acts as a memory for the hidden layer output.Ultimately, the output layer generates predictions for SPI durations, predicated upon the acquired knowledge of relationships between input variables.The ENN's recurrent connections enable it to factor in the influence of past observations, rendering it well-suited for capturing dynamic associations between meteorological and groundwater events.
Further refinement of this architecture was achieved through the FA optimization process, enhancing the accuracy of SPI predictions and determining the most suitable network architecture and weightings for the ENN.The FA draws inspiration from nature, specifically the flashing behavior of fireflies, a concept initially proposed by Yang and He [43].It simulates the movement of fireflies in their quest for optimal solutions within a multidimensional parameter space.In this study, the parameters subjected to optimization encompassed the network architecture (e.g., the number of hidden layers and neurons per layer) and the weights associated with network connections.This tuning process unfolds in several key steps: (i) Initialization: the FA begins by initializing a population of fireflies, where each firefly corresponds to a potential solution, defined by a set of network parameters.(ii) Brightness calculation: each firefly's brightness is determined based on its fitness value, reflecting how well its associated ENN configuration performs in predicting the association between SPI and SGL durations.(iii) Firefly movement: fireflies gravitate towards brighter ones within the search space, guided by a mathematical expression that considers their relative positions and brightness values.This movement aims to converge towards optimal configurations.(iv) Parameter updating: as fireflies move, their parameters are updated to mirror their new positions.These updates translate into adjustments to the ENN's architecture and weights.(v) Iterative process: the movement and parameter updating steps iterate through multiple cycles, enabling fireflies to refine their positions and approach superior solutions.(vi) Convergence: the FA continues to refine firefly positions until a predefined stopping criterion is met, signaling the convergence of the optimization process.
By harnessing the ENN tuned via the Firefly Algorithm, our objective was to discern the optimal network architecture and weightings that facilitate precise predictions of SPI durations predicated on associated SGL durations.For an in-depth exploration of the approach involving the tuning of ANN models via FA, we recommend consulting the work by Mohammadi [21].

Performance Indicators
In order to comprehensively assess the accuracy and reliability of our predictive models, we employed a battery of four distinct performance indicators: the Normalized Root Mean Square Error (NRMSE), the Nash-Sutcliffe Efficiency (NSE), the coefficient of determination multiplied by the slope of the regression line (bR 2 ), and the Kling-Gupta Efficiency (KGE) [44].These metrics collectively offer insights into various facets of model performance, encompassing aspects such as error magnitude, precision, linearity, and the overarching concordance between simulated and observed values.In essence, these performance indicators collaboratively furnish a comprehensive evaluation of the models' predictive aptitude, the extent of modeling error, and the overall consistency between the predicted and observed durations of SPI and SGL.

Results and Discussion
As elucidated previously, the initial assessment of the selected stations' SGL and SPI values hinged upon the careful examination of time series plots and the geographical congruence between meteorological stations and observation wells.As depicted in Figure 3, the spatial dimension, defined by Thiessen polygons (as previously illustrated in Figure 1), and the temporal dimension, characterized by the resemblance between time series, do not readily reveal a pronounced similarity between monthly SGL and SPI data.However, the broader-scale association between these selected time series does underscore the pertinent connection between precipitation and groundwater abundance.In accordance with the findings of Han et al. [45], who probed the effect of meteorological drought propagation on groundwater level perturbations, it may take several months to discern the groundwater level disturbances linked to an ongoing meteorological drought.Nevertheless, the augmentation of soil moisture expedites the response of groundwater to precipitation during surplus periods due to the stored-full runoff mechanism.Conversely, intensified evapotranspiration and heightened heat wave rates amplify drought propagation during deficit periods.In a recent study by Gong et al. [46], they explored the effects of lag, pooling, lengthening, and attenuation in monitoring the propagation of meteorological drought on groundwater table abundance.Their conclusions highlighted that the impact of lag time in discerning the connection between meteorological and groundwater drought is influenced by factors such as aquifer depth, soil permeability, anthropogenic activities, temporal land use shifts, and prevailing climate conditions.
Therefore, bearing in mind the inherent challenges of scrutinizing short-term periods (i.e., monthly data) in establishing the link between SGL and SPI, the insights gleaned from Figure 3 do provide partial indications of the correlation between precipitation and groundwater levels in the region.Furthermore, considering the local geomorphology, the hydraulic gradient between high-altitude regions and hillsides serves as a compelling rationale for the existence of a relationship between these two variables, namely groundwater and precipitation.This intriguing proposition is subject to further exploration in the ensuing cross-correlation diagram presented in Figure 4.
In the majority of cases, as illustrated in Figure 4, the lagging and leading effects appear to be inconsequential, with the linear correlation scarcely surpassing a meager 0.3.This trend is particularly evident for stations like Delicay, Cayirkoy, and Kursunlu; however, these stations exhibit a somewhat stronger association.Consequently, the limited linear correlation between SPI and SGL time series fails to furnish compelling evidence for a robust linear relationship between precipitation and groundwater levels.
As previously delineated, to enrich our analysis and substantiate the connection between groundwater abundance and meteorological occurrences, we delve into the comparison and prediction of SPI event durations in tandem with SGL durations.Given that neither groundwater nor surface water are likely to exert a significant influence on meteorological precipitation, the prediction of SPI based on SGL offers an alternative perspective on ongoing events within the region.In pursuit of this, we initially derive duration curves for both SGL and SPI, as illustrated in Figure 5.The results reveal a noteworthy resemblance and symmetry between the durations of SGL and SPI.Nonetheless, SPI exhibits some outliers, primarily for the most frequent and extreme events.Evidently, unexpected occurrences, such as extreme or short-term precipitation events, exhibit a weaker association with disturbances in groundwater levels.However, when we consider the broader context, it becomes evident that a substantial correlation and association prevail between hydro-meteorological events, which may occasionally be perturbed by anthropogenic or climate-induced factors.In the subsequent phase, aimed at substantiating the link between groundwater and precipitation, SPI durations are prognosticated utilizing SGL durations randomly drawn from the duration curves (as depicted in Figure 5).During this process, the Thiessen polygons and SVM techniques are employed interchangeably to determine the ensemble of groundwater observation wells (Scenarios) utilized as predictors for SPI duration at the meteorological stations.
The Thiessen polygons approach suggests (Figure 1) the inclusion of the Narlidere, Delicay, Adakoy, and Yenice observation wells in estimating SPI duration at the Uludag meteorological station.Similarly, according to the Thiessen method, the Aticilar, Yeniceabat, and Cayirkoy observation wells should be incorporated to ascertain SPI durations in Osmangazi, while the Kursunlu observation well is identified as being pertinent for determining SPI duration at the Keles station.In the majority of cases, as illustrated in Figure 4, the lagging and leading effects appear to be inconsequential, with the linear correlation scarcely surpassing a meager 0.3.This trend is particularly evident for stations like Delicay, Cayirkoy, and Kursunlu; however, these stations exhibit a somewhat stronger association.Consequently, the limited linear correlation between SPI and SGL time series fails to furnish compelling evidence for a robust linear relationship between precipitation and groundwater levels.
As previously delineated, to enrich our analysis and substantiate the connection between groundwater abundance and meteorological occurrences, we delve into the comparison and prediction of SPI event durations in tandem with SGL durations.Given that neither groundwater nor surface water are likely to exert a significant influence on meteorological precipitation, the prediction of SPI based on SGL offers an alternative perspective on ongoing events within the region.In pursuit of this, we initially derive duration curves for both SGL and SPI, as illustrated in Figure 5.The results reveal a noteworthy resemblance and symmetry between the durations of SGL and SPI.Nonetheless, SPI exhibits some outliers, primarily for the most frequent and extreme events.Evidently, unexpected occurrences, such as extreme or short-term precipitation events, exhibit a weaker association with disturbances in groundwater levels.However, when we consider the broader context, it becomes evident that a substantial correlation and association prevail between hydro-meteorological events, which may occasionally be perturbed by anthropogenic or climate-induced factors.In the subsequent phase, aimed at substantiating the link between groundwater and precipitation, SPI durations are prognosticated utilizing SGL durations randomly drawn from the duration curves (as depicted in Figure 5).During this process, the Thiessen polygons and SVM techniques are employed interchangeably to determine the ensemble of groundwater observation wells (Scenarios) utilized as predictors for SPI duration at the Conversely, the SVM approach recommends the utilization of specific observation wells (as detailed in Table 2) for forecasting SPI durations at the Keles, Uludag, and Osmangazi stations, based on various kernel functions, including dot, radial, polynomial, neural, ANOVA, and Epachnenikov kernels.According to the results outlined in Table 2, the SPI duration at the Uludag station exhibits a significant association with most of the groundwater level records, except those from Yenice and Yeniceabat.Despite the geographic distance between Yeniceabat and the Uludag station, the findings align with those obtained via the Thiessen polygons method, which suggests the inclusion of SGL records from the Yenice observation well in estimating SPI durations in Uludag.Similarly, the SVM method posits that the Narlidere observation-well records should be included in determining the SPI in Osmangazi, while Narlidere, Yeniceabat, and Kursunlu observationwell records should be excluded from the estimation of SPI durations at the Keles station.Notably, despite Narlidere's geographical distance from the Osmangazi station, the results pertaining to the Keles station exhibit substantial disparities, notably in rejecting the role of SPI duration in the SGL duration at the Kursunlu observation well-a unique groundwater observation well introduced in the prior Thiessen-based stage.Hence, in order to conduct further scrutiny and validate the outcomes derived from the Thiessen and SVM methodologies, we proceeded to estimate SPI durations using the corresponding SGL durations at the recommended observation wells through the ENN-FA approach.Initially, the ENN was employed to gauge SPI durations, and subsequently, the FA algorithm was deployed to bolster the performance metrics and assess the potential for scenario-models to yield superior results.In this context, Table 3 furnishes a comprehensive evaluation of the ENN-FA approach's performance in delineating the relationship between SPI and SGL during the testing phase.The resultant findings are also graphically illustrated in the box and whisker plots featured in Figure 6.As delineated in Table 3, it is apparent that SVM generally yields more favorable outcomes in comparison to the Thiessen method.This essentially underscores the fact that geographic proximity and mathematical approaches may fall short in encapsulating the intricate hydrogeological dynamics.Nevertheless, in scenarios where no alternatives are As delineated in Table 3, it is apparent that SVM generally yields more favorable outcomes in comparison to the Thiessen method.This essentially underscores the fact that geographic proximity and mathematical approaches may fall short in encapsulating the intricate hydrogeological dynamics.Nevertheless, in scenarios where no alternatives are available or simpler methodologies are preferred, the Thiessen polygon method retains its reliability.
When we focus on SVM as the most effective scenario-model predictor, it becomes evident that Radial-ENN-FA for Uludag, ANOVA-ENN-FA for Osmangazi, and Dot-ENN-FA for Keles can be seen as the best-performing combinations, displaying minimal bias and high efficacy.However, when addressing outliers and probability distributions (including mean, median, quartiles, etc.), Polynomial-ENN-FA for Uludag, Epachnenikov-ENN-FA for Osmangazi, and Dot-ENN-FA for Keles emerge as the superior choices.Given the close correlation among the performance metrics for the scenario-models presented in Table 3, one may find it prudent to employ a scatter plot approach to ensure compliance with both outcome probability and outlier criteria.Consequently, in cases where information is lacking or geological data are absent, this approach may serve as a robust means for achieving successful predictions.
The degree of association and predictability of SPI with the assistance of SGL also attests to the natural connection between groundwater and precipitation in the region.Consequently, it is more plausible that groundwater disturbances in the region are a consequence of climatic factors rather than anthropogenic influences.This holds particular relevance in regions like Bursa, where groundwater is not heavily relied upon as the primary water source.Nevertheless, given the increasing prevalence of drought and population growth in the area, it is foreseeable that more comprehensive investigations will be necessitated in the near future to elucidate the behavior of groundwater levels in the region.

Conclusions
In this research, we conducted an in-depth investigation into the nexus between meteorological drought and fluctuations in groundwater levels in the vicinity of Mount Uludag in Bursa, Turkey.Employing a rigorous analytical approach encompassing meteorological and groundwater data from multiple observation sites, our objective was to unveil the underlying intricacies that connect SPI events and SGL disturbances within the region.Our empirical findings unveil a nuanced and intricate relationship between hydro-meteorological phenomena within the study area.
Spatial and temporal analyses have indicated that although an overt direct spatial and temporal correlation between SGL and SPI time series at the monthly scale may not be readily discernible, a more extensive-scale relationship indeed exists, underscoring the significance of the interplay between precipitation patterns and groundwater abundance.We conducted a meticulous comparison of the duration curves for SGL and SPI events, introducing the concept of employing SGL duration as a predictive measure for SPI duration.
Through the application of both the Thiessen polygon and SVM techniques, we identified the optimal observation wells that demonstrate the highest predictive capacity for SPI duration at meteorological sites.Our results showcased the consistency and complementarity of insights provided by both methodologies in elucidating the connection between meteorological and groundwater occurrences.Furthermore, we harnessed the power of the ENN fine-tuned via the FA to forecast SPI durations utilizing SGL durations.This novel approach enabled us to explore the feasibility of accurately forecasting meteorological events based on groundwater behavior.Our analysis unveiled that the ENN-FA models yielded promising results, furnishing enhanced insights into the intricate relationship between SPI and SGL durations.
The present study significantly contributes to an improved comprehension of the complex dynamics governing the interplay between meteorological drought and groundwater level fluctuations.While the detection of short-term associations may pose challenges, our findings robustly support the existence of a broader nexus between precipitation patterns and groundwater abundance.The amalgamation of innovative methodologies, encompassing SVM, ENN-FA, and duration curve analysis, represents a comprehensive approach for unraveling these intricate relationships.Future research endeavors could delve deeper into the influence of various hydrogeological factors on these observed associations, thereby facilitating the development of more robust predictive models.

Figure 1 .
Figure 1.Global positioning of the study area (up left), climate classification [25] (down left), and the selected stations together with the Thiessen polygons (right) used in the analysis.

Figure 1 .
Figure 1.Global positioning of the study area (up left), climate classification [25] (down left), and the selected stations together with the Thiessen polygons (right) used in the analysis.

Figure 2 .
Figure 2. Flowchart of the methods used in the analysis.

Figure 2 .
Figure 2. Flowchart of the methods used in the analysis.

Figure 3 .
Figure 3.Time series plot of the SPI and SGL defined for the (a) Uludag, (b) Osmangazi, and (c) Keles regions determined by means of Thiessen polygon classification detailed in Figure 1.

Figure 3 .
Figure 3.Time series plot of the SPI and SGL defined for the (a) Uludag, (b) Osmangazi, and (c) Keles regions determined by means of Thiessen polygon classification detailed in Figure 1.

Figure 5 .
Figure 5. Drought duration curve for SGL and SPI.

Figure 5 .
Figure 5. Drought duration curve for SGL and SPI.

Figure 6 .
Figure 6.Box and whisker plots of the obtained results in the test stage for the (a) Uludag, (b) Osmangazi, and (c) Keles stations.

Figure 6 .
Figure 6.Box and whisker plots of the obtained results in the test stage for the (a) Uludag, (b) Osmangazi, and (c) Keles stations.

Table 1 .
Properties of the stations and associated data.

Table 1 .
Properties of the stations and associated data.

Table 2 .
Selection of the best groundwater observation wells for the prediction of SPI based on the SVM method (grey highlights indicate the significant kernels larger than 0.75 that used in selection of the stations).

Table 3 .
Results obtained for the ENN-FA models based on the groundwater selection scenarios at the test stage.