Next Article in Journal
Assessing the Impacts of Climate Change and Water Extraction on Thermal Stratification and Water Quality of a Subtropical Lake Using the GLM-AED Model
Previous Article in Journal
A Half-Century of Human Impact on Nan River Runoff and Sediment Load Supplied to the Chao Phraya River
Previous Article in Special Issue
Estimating a Reliable Water Budget at a Basin Scale: A Comparison between the Geostatistical and Traditional Methods (Foro River Basin, Central Italy)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Towards Groundwater-Level Prediction Using Prophet Forecasting Method by Exploiting a High-Resolution Hydrogeological Monitoring System

1
Dipartimento di Scienze e Ingegneria della Materia, dell’Ambiente ed Urbanistica (SIMAU), Università Politecnica delle Marche, 60131 Ancona, Italy
2
Dipartimento di Ingegneria dell’Informazione, Università Politecnica delle Marche, 60131 Ancona, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Water 2024, 16(1), 152; https://doi.org/10.3390/w16010152
Submission received: 29 November 2023 / Revised: 23 December 2023 / Accepted: 27 December 2023 / Published: 30 December 2023

Abstract

:
Forecasting of water availability has become of increasing interest in recent decades, especially due to growing human pressure and climate change, affecting groundwater resources towards a perceivable depletion. Numerous research papers developed at various spatial scales successfully investigated daily or seasonal groundwater level prediction starting from measured meteorological data (i.e., precipitation and temperature) and observed groundwater levels, by exploiting data-driven approaches. Barely a few research combine the meteorological variables and groundwater level data with unsaturated zone monitored variables (i.e., soil water content, soil temperature, and bulk electric conductivity), and—in most of these—the vadose zone is monitored only at a single depth. Our approach exploits a high spatial-temporal resolution hydrogeological monitoring system developed in the Conero Mt. Regional Park (central Italy) to predict groundwater level trends of a shallow aquifer exploited for drinking purposes. The field equipment consists of a thermo-pluviometric station, three volumetric water content, electric conductivity, and soil temperature probes in the vadose zone at 0.6 m, 0.9 m, and 1.7 m, respectively, and a piezometer instrumented with a permanent water-level probe. The monitored period started in January 2022, and the variables were recorded every fifteen minutes for more than one hydrologic year, except the groundwater level which was recorded on a daily scale. The developed model consists of three “virtual boxes” (i.e., atmosphere, unsaturated zone, and saturated zone) for which the hydrological variables characterizing each box were integrated into a time series forecasting model based on Prophet developed in the Python environment. Each measured parameter was tested for its influence on groundwater level prediction. The model was fine-tuned to an acceptable prediction (roughly 20% ahead of the monitored period). The quantitative analysis reveals that optimal results are achieved by expoiting the hydrological variables collected in the vadose zone at a depth of 1.7 m below ground level, with a Mean Absolute Error (MAE) of 0.189, a Mean Absolute Percentage Error (MAPE) of 0.062, a Root Mean Square Error (RMSE) of 0.244, and a Correlation coefficient of 0.923. This study stresses the importance of calibrating groundwater level prediction methods by exploring the hydrologic variables of the vadose zone in conjunction with those of the saturated zone and meteorological data, thus emphasizing the role of hydrologic time series forecasting as a challenging but vital aspect of optimizing groundwater management.

1. Introduction

Water is indisputably one of the most important natural resources since its availability sustains life, agriculture, industry, and its dependent ecosystems [1,2,3]. A continuously increasing strain on freshwater supplies is observed due to world population growth and climate change [4,5], hence water conservation practice and water availability prediction emerge as a vital objective [6,7]. Indeed, while the causes of climate change are still debated and often divergent, the produced effects are globally recognized. The most common include a rise in dry periods and an increase in intense rainfall events. Those phenomena have a strong impact on both groundwater and surface water resources. In fact, while large amounts of rainfall occurring in a few hours can cause floods, a gradual depletion of groundwater available for human consumption is observed. In this context, some authors suggest that a shift in monitoring focus from the groundwater to the vadose zone is necessary to understand the processes that regulate the decreasing of aquifer recharge [8,9,10,11]. During the last decades, special attention has been paid to monitor the unsaturated zone layers to characterize infiltration processes and thus estimate groundwater recharge vs. run-off generation [12,13]. The unsaturated zone represents the portion that encompasses the topographic surface until the top of the groundwater body, where most of the biogeochemical processes take place and where rainwater infiltrates [14]. At the beginning of the 21st century, the unsaturated zone was included, according to the US National Research Council, in the more elaborated concept of “Critical Zone”, because of the extreme complexity of the phenomena that regulate it [15]. The various heterogeneities in the unsaturated zone (presence of roots, lenses of fine or coarse materials, different gas phases, and redox conditions) are reflected in a heterogeneous response in terms of recharge, transport of pollutants, and their attenuation. Hence, understanding unsaturated-zone processes is crucial to determine the amount and quality of groundwater that is available for human and ecosystem use [16,17]. In this context, the growing interest in groundwater level (GWL) prediction to support water management operations was accompanied by a proliferation of advanced sensors and data collection technologies [18,19]. Usually, Time Domain Reflectometry (TDR) moisture probes are used to monitor the soil hydrologic variables such as water content, soil temperature, and bulk electric conductivity, providing useful information to understand aquifer recharge mechanisms. These sensors can gather vast amounts of hydrologic data at different scales. However, few hydrogeological studies still analyze the aquifer response to meteoric inflow by investigating the vadose zone. Most of the studies deal with precipitation and air temperature data to perform water balance by exploiting the most known empirical methods [20,21,22] or more complex physically based methods [23,24]. A recent work by Berthelin et al., 2023 [25] exploits soil moisture content only measured at a single depth (20 cm) at different locations, to determine the aquifer recharge of the investigated catchment.
The increasing improvement of monitoring systems and the large amount of hydrologic data collected by the scientific community fostered the development of data-driven methods to characterize and understand complex natural phenomena [26]. Several studies are based on the application of machine learning approaches in geosciences and related subjects [27,28,29,30,31]. Utilizing data-driven approaches toward the prediction of groundwater level is not a new phenomenon, and traditionally, numerical methods have been used for groundwater level modeling [32,33,34]. However, recent studies have extensively employed Artificial Intelligence (AI) based techniques [35,36,37]. Due to the inherent non-linear and non-stationary nature of groundwater level time series, intelligent data-driven methodologies have showcased promising results. Across the literature, various popular forecasting approaches have been tested on specific applications of groundwater level forecasting, including Auto-Regressive Integrated Moving Average (ARIMA), Artificial Neural Networks (ANN), Long Short-Term Memory (LSTM) as well as hybrid approaches such as ARIMA-LSTM [38,39,40,41]. Comparative studies have consistently shown that machine learning-based methods outperform traditional numerical approaches [42] with superior prediction performance and capturing complex and non-linear relationships between input and output variables [43]. In a bibliometric study of machine learning and mathematical modeling techniques of forecasting using piezometric data [44], authors find that machine learning techniques such as Random Forest (RF), Support Vector Machine (SVM), and deep learning techniques like ANN achieve higher accuracies if compared to mathematical model techniques. Ren et al. [45] investigated that though deep learning performs well in filling high dynamic gaps, it struggles with reconstructing trends and seasonality-based gaps. On the other hand, ARIMA, a traditional machine learning model, excels in capturing trends and seasonality.
This suggests that traditional models, designed to handle time series with seasonal effects and trends, may be more suitable for groundwater level forecasting tasks. While RFs and SVMs lack built-in mechanisms to handle the temporal nature of groundwater data and explicitly capture seasonality and holiday effects in groundwater-level time series, ANNs require careful architecture design, training, and tuning for time series forecasting. A pressing need exists for a fast, accurate, and tunable forecasting procedure that works best with time series with strong seasonal effects due to the nature of groundwater level oscillation during the hydrologic years. Further, it should also have an easy mechanism of incorporation of exogenous variables while maintaining interpretability. Prophet is an open-source machine learning model specifically designed for time series forecasting, with a particular emphasis on capturing seasonality, trends, and holiday effects [46]. Prophet is robust in automatically handling missing data and outliers, is flexible in incorporating domain knowledge, and can capture seasonality and holiday effects with simplicity. Moreover, Prophet incorporates uncertainty estimation, which is crucial for decision-making in groundwater management. Researchers have tested it for groundwater level estimation, and it has been compared with other models like ARIMA, Multivariate Adaptive Regression (MARS), and Error Trend and Seasonality (ETS) using satellite data to analyze and forecast groundwater level in the Urmia Lake basin (Northwestern Iran) [47]. Prophet consistently outperformed the other models, achieving higher coefficients of determination ( R 2 ) ranging from 0.81 to 0.85, as well as lower Mean Absolute Error (MAE) and Root Mean Squared Error (RMSE). Aguilera [48] conducted a similar study at the Ramsar wetland area of Doñana (Spain) and found that Prophet exhibited strong prediction capabilities for groundwater level time series. The study compared Prophet with various statistical and intelligent models commonly used in time series forecasting, suggesting that methods with additive schemes, such as ARIMA and Prophet, are suitable for modeling groundwater time series. Though it has demonstrated superior performance, the existing literature only focuses on a simple implementation based on regressors such as air temperature and precipitation, which are not the most robust representations of the actual groundwater recharging mechanism.
As far as our best knowledge, no previous research aims to estimate future groundwater availability by combining a vadose zone monitoring system collecting data at different depths into the soil and machine learning approaches. Only a few studies in the literature exploit the Prophet forecasting method to estimate future groundwater levels by exploiting just previously observed groundwater level data [47,48] or precipitation data [42]. To fill this gap in the literature, this research explores the additive machine learning model Prophet (not specifically developed for hydrological purposes) to predict the groundwater level of an alluvial aquifer exploited for drinking purposes by using the vadose zone monitored data. The study site is located in the Conero Mt. Regional Park (central Italy). The monitoring strategy consists of a thermo-pluviometric station, three soil water content, electric conductivity, soil temperature probes, and a piezometer instrumented by a water depth probe. All the variables collected from the monitoring system are tested for their influence on groundwater level prediction. Towards this approach, the proposed research tries to address the following open-ended questions:
  • How do unsaturated zone variables, like soil volumetric water content, bulk electric conductivity, and soil temperature improve groundwater level predictions compared to precipitation and air temperature data and previously observed groundwater levels?
  • Can our proposed Prophet model accurately forecast groundwater levels in a shallow aquifer by integrating high-resolution hydrological monitoring, both from meteorological data and both from the vadose zone data?
  • What is the impact of monitoring the vadose zone at different depths on groundwater level predictions in the study area?

2. Materials and Methods

This section presents the materials and methods employed in the investigation of advanced groundwater level prediction by applying the Prophet forecasting method, leveraging the capabilities of a high-resolution hydrological monitoring system. The subsections encompass a comprehensive overview of our research, starting with a description of the study area, the monitoring strategy, and the collected dataset. Subsequently, we delve into the forecasting method and the choice of predictive regressors. Through these carefully integrated methodologies, we aim to contribute to developing an accurate and originally applied forecasting system for groundwater availability and management.

2.1. Hydrogeological Features of the Study Site

The study site is located in the central Mediterranean basin. The research is focused on a sub-catchment (about 5 km2) of the Aspio watershed, situated within the Conero Mt. Regional Park (central Italy) (Figure 1a). The Conero Mt. hydrostructure has been the objective of several studies during past years [49]. Fronzi et al., 2022 [50] identified two main aquifers in the area. The first one, named Scaglia Calcarea aquifer is hosted in the Scaglia Rossa and Scaglia Bianca geological Formations (Fms.), while the second one is hosted in the terraced and alluvial deposits of the Quaternary age (Figure 1b). The Scaglia Calcarea represents the main aquifer and is constituted by stratified mycritic limestones, marly, and flinty limestones [49]. Its hydraulic conductivity is regulated by the presence of fissures and microkarst features, typical of such carbonate geological Fms. [51]. The Scaglia Calcarea aquifer is inferiorly confined by the marly units of the Marne a Fucoidi Fm., (not outcropping in the analyzed basin). At the top, the Scaglia Calcarea aquifer is semi-confined by the presence of a low permeability complex of the Scaglia Variegata Fm., Scaglia Cinerea Fm., the Bisciaro Fm. and the Schlier Fm., which constitute a single aquiclude. This hydrogeological complex, together with Pliocene and Pleistocene clay layers, hydraulically separates the Scaglia Calcarea aquifer from the shallow alluvial aquifer mainly comped by silty-sands, sandy-silts and gravels with sand layers. The alluvial aquifer is responsible for the presence of perennial surface water in the tributaries of the Aspio River and feeds the watercourses through a well-explored groundwater-surface water interaction [50]. As a result, the groundwater body hosted in the alluvial aquifer is responsible for sustaining the aquatic ecosystem of the Regional natural park. The aquifers’ recharge is entirely due to meteoric precipitation which accounts for an average precipitation of about 900 mm/year [49]. Both aquifers play an important role in the drinking water supply. Indeed, they are tapped by the local water management company, providing good quality water in the small towns nearby. However, the meteoric regime of the years 2019, 2020, and 2021 has impacted the groundwater resources in the area with severe groundwater depletion effects and a concomitant drying up of the water course for a prolonged period [50]. This aspect provided the impetus for an in-depth assessment of future groundwater availability by characterizing the infiltration processes and the water movement in the vadose zone towards the alluvial aquifer to support local management companies and authorities. For this reason, starting in January 2022, a high-resolution hydrogeological monitoring system has been set up in the basin. The following chapter describes the implemented monitoring system and the data collection strategy.

2.2. Monitoring System and Data Collection

Starting from November 2021, field operations have been made to implement the monitoring site. Three holes of diameter 0.1 m were manually drilled at a distance of 2 m from each other and inclined 45 degrees to the horizontal. The holes were drilled at 0.6, 0.9, and 1.7 m depth, respectively, into the soil deposits of the alluvial aquifer in the proximity of the wells field employed for drinking water supply (Figure 1b). In the holes, the advanced soil moisture, temperature, and electric conductivity (EC) sensors TEROS 12 (Meter Group Inc., Pullman, WA, USA) were installed at the bottom, taking care to refill and compact the soil previously extracted from the hole over the sensors. This procedure, together with the drilling orientation (45°), was made to minimize the impact of the hole on preferential infiltration flow paths into the soil toward the sensors. The TEROS 12 are advanced sensors that exploit the TDR principle to monitor the vadose zone hydro-physical properties. The soil moisture is expressed as volumetric water content (VWC), measured in m 3 /m 3 , with a resolution of 0.001 m 3 /m 3 and accuracy ±0.03 m 3 /m 3 . The soil temperature (T s o i l ) operational range is −10 to +60 °C, with resolution of 0.1 °C, and accuracy ±0.5 °C from −10 to 0 °C and ±0.3 °C from 0 to 60 °C. The Bulk electrical conductivity (EC) range is 0–20 mS/cm with 0.001 mS/cm resolution and ±5 percent accuracy. The precipitation (P) regime and the air temperature (T a i r ) in the area are monitored through a thermo-pluviometric station, while a standpipe piezometer located in the alluvial aquifer (about 200 m downstream from the vadose zone monitoring point) was equipped with a hydrometric pressure transducer (TD-Diver Eijkelkamp, accuracy ±0.5 cmH2O and resolution 0.2 cmH2O) compensated by atmospheric pressure, for continuously monitoring groundwater level fluctuation (Figure 1b). The monitored period effectively started on 1 January 2022, and all the hydrological parameters were recorded every fifteen minutes for more than one hydrologic year until 30 April 2023, except for the GWL, which was recorded every day at noon. We utilize IoT technology to automatically collect data on the cloud robustly and efficiently, as demonstrated by Galdelli et al., 2019; Galdelli et al., 2021 and Tassetti et al., 2022 [52,53,54]. The obtained high spatial and temporal resolution dataset was used to develop the GWL prediction model, as shown in the next chapter.

2.3. Groundwater Level Forecasting

The implementation of the forecasting algorithm was preceded by a schematic simplification of the natural conditions, according to which the system was divided into three virtual boxes: Box 1 = Atmosphere (and related variables), Box 2 = Vadose zone (and related variables collected at different depths), and Box 3 = Saturated zone (and related variables), following the scheme proposed in Figure 2. The implementation is split into four distinct scenarios that comprehensively assess the model’s predictive capabilities. Initially, variables from Box 1, encapsulating atmospheric conditions, were employed to predict the GWL in Box 3. Subsequently, three additional tests were conducted using variables from Box 2 at depths of 0.6 m, 0.9 m, and 1.7 m, independently to forecast the GWL in Box 3. Each of these scenarios provides valuable insights into the relationships between specific environmental variables and groundwater levels, contributing to an understanding and optimization of the forecasting model for diverse conditions within the system.
Considering the state of the art across the literature we designed a Prophet-based system for forecasting purposes [46]. Prophet is a time series forecasting model developed by the core data science team at Facebook and is an open-source project for analyzing and forecasting time series. Prophet offers numerous strengths, including the robustness of the model in the case of missing data, frequent trend changes, and promising performance even in the presence of outliers. The convenient methods to add additional regressors and robust cross-validation process make it a promising forecasting mechanism in our case study. Prophet represents the time series as the sum of three components: (i) trend, (ii) seasonality, and (iii) holidays, as seen below:
y ( t ) = g ( t ) + s ( t ) + h ( t ) + ϵ t
where:
  • g(t): trend function, models non-periodic changes, it can be logarithmic;
  • s(t): seasonality function, relying on the Fourier series, provides a flexible model of periodic effects to model changes that are repeated at regular time intervals (e.g., weekly and yearly seasonality); it is also possible to have more than one seasonality in the same series;
  • h(t): holidays, models irregular events that temporarily alter the time series;
  • ϵ t : error term, represents changes in the time series that the model does not capture, ϵ t is regarded as a normal distribution.
Prophet adopts a unique approach to forecasting, treating it as a curve-fitting problem in contrast to other methods like ARIMA [55]. Unlike ARIMA or SARIMAX, which relies on autocorrelation and partial autocorrelation to capture temporal dependencies, Prophet decomposes the input time series into additive components. It models trends using a piecewise linear logistic growth curve and incorporates seasonality through Fourier series expansion. The data needs to be converted to a proper time series format to facilitate effective utilization, as the model relies on a structured temporal sequence for accurate trend and seasonality identification. Due to the disparity in temporal resolution between the hydrological variables recorded at fifteen-minute intervals (i.e., P, T a i r , VWC, T s o i l and EC) and the GWL observations, collected at daily scale, a data treatment has deemed necessary. Specifically, the hydrological variables of Box 1 and Box 2 were resampled at the daily scale as follows. The daily mean T a i r and daily cumulative P were computed. Similarly, the daily mean VWC, T s o i l and EC were computed for all different depths, (i.e., 0.6 m, 0.9 m, and 1.7 m). Consequently, the method ensured that all the collected hydrological variables were converted into a consistent daily format, aligning them temporally with the observed GWL.
Training the model involves a combination strategy utilizing grid search for hyperparameter tuning and cross-validation for predictive error assessment on validation data, as depicted in Figure 3. The grid search algorithm systematically explores a range of hyperparameter values, as detailed in Table 1, exhaustively testing all possible combinations for optimization. Hyperparameters such as prior_scale_temperature, prior_scale_rain for the atmosphere, prior_scale_soil_temperature, prior_scale_water_content, and prior_scale_ electric_conductivity and were optimized for this investigation. Moreover, depending on the specific forecasting scenario, weights for exogenous variables acting as regressors were introduced. Such exogenous variables are not directly influenced by the groundwater level, but are important factors affecting it. Introducing these exogenous variables as regressors and assigning weights through hyperparameter tuning, our model can account for external influences, providing a more comprehensive and accurate prediction of GWL. Therefore, to incorporate these variables, depending on the forecasting scenario, hyperparameters such as prior_scale_temperature, prior_scale_rain for the atmosphere (Box 1), prior_scale_soil_temperature, prior_scale_water_content, and prior_scale_electric_conductivity for the vadose zone (Box 2) were defined and custom integrated into the model. These weights allow the model to assign varying degrees of importance to different exogenous variables, enhancing its adaptability to the specific dynamics of the saturated zone (Box 3) namely the GWL prediction task. The systematic selection of optimal hyperparameters, combined with the introduction of these weights, ensures the model’s robust performance and facilitates accurate forecasts for different scenarios.
Cross-validation evaluates the model’s performance across multiple training and validation sets based on forecasting requirements. The determination of the training data, which is used for cross-validation, is guided by key parameters: initial, period, and horizon. These parameters are built into the Prophet and play a crucial role in defining the size of the initial training period and the duration of the forecast window. The period parameter specifies the length of a seasonal cycle, and the horizon parameter defines the duration for which future predictions are made. The initial parameter determines the size of the initial training period. During each cross-validation fold, the model is trained on this initial period, allowing it to learn from historical data. The process iterates across the entire time series, providing a robust assessment of the model’s performance across different data segments. Specifically, our cross-validation approach involves dividing the time series into segments, each equivalent to roughly 10% of the total time series into validation and testing. In the validation process, a timeframe of 2 months was established, mirroring the identified seasonal patterns and effectively capturing the recurring nature of the observed data. Additionally, the test dataset spans 2 months, ensuring a comprehensive evaluation of the model’s performance.
A period of 7 days was deemed appropriate in line with the identified seasonal patterns, and this duration captures the recurring nature of the observed data and aligns with the underlying cyclicality. Additionally, a horizon of 14 days was chosen to represent the forecast window. This window length meets the local groundwater management company’s request to operate in time to ensure the availability of drinking water for the population living nearby. Furthermore, the selection of 14 days ensures that each cross-validation fold covers a significant portion of the time series while allowing for a robust assessment of the model’s performance. This choice allows for a two-week projection, providing a meaningful time frame for anticipating near-future trends and variations without significant loss in forecast accuracy. Therefore, in case of a predicted significant groundwater level decrease the management company can locate water pumps at different depths within the wells, or integrate the amount of extracted groundwater from emergency well fields or different water plants.

2.4. RAPS Method

The Rescaled Adjusted Partial Sums (RAPS) method is a fundamental analytical technique within the realm of hydrogeological time series since it can discern nuanced fluctuations within collected data, making it an essential tool for groundwater studies [56,57]. RAPS works by aggregating the deviations from a specified mean level of hydrological variables and adjusting them against the data’s standard deviation, as reported in Equation (2).
RAPS N = i = 1 N Y i Y ¯ S Y
where:
  • N: the amount of data in the time series;
  • Y i : the value of an individual sample i = 1, 2, …, N;
  • Y ¯ : average value of the observed sample;
  • S Y : the value of the standard deviation of the time series.
The RAPS method was used to analyze observed GWL and the predicted ones at a daily scale by using two time-series lengths. The first one considers the duration of the whole monitoring period, while the second one considers only the testing dataset. This approach was employed to compare trends, periodicities, and fluctuations [58] on observed vs. predicted GWL.

3. Results

3.1. Hydrological Characterization and Collected Dataset

The collected hydrological variables during the monitoring period (1 January 2022–30 April 2023) are reported in the graphs of Figure 4. When examining the graphs, it is possible to observe that the precipitation regime in the area is characterized by wet periods during winter and spring and less rainer periods during summer and early autumn. The summer is marked by sporadic but heavy rainfall events that frequently exceed 10 mm/15 min, while the wet periods are characterized by more continuous rainfall events persisting for several hours to days. The most intense rainfall event can be depicted on 15 September 2022 (12 mm/15 min), recognized also by other authors to be one of the most intense precipitation events affecting the northern Marche Region during the last century [59,60,61].
Concerning the air temperature, it ranges between 0 and 35 °C during the hydrologic year. The higher values are observed during the summer (in July), and the lower are recorded during the winter periods between January and February. The soil temperature mimics the trend of the air temperature at all the monitored depths with a time lag between T a i r and T s o i l increasing as the depth increases. Three thermal stationary periods can be observed in April and October 2022 and between March and April 2023. During the other months, the soil temperature is higher at 0.6 m and lower at 1.7 m during summer. On the contrary, during winter, an inverse temperature gradient is observed, as expected. Indeed, the T s o i l at 0.6 m reflects changes in T a i r faster, while deeper down, the thermal inertia of the soil becomes more sizeable, smoothing temperature variations. T s o i l at 0.6 m ranges between 6.7 and 25.5 °C, T s o i l at 0.9 m ranges between 8.3 and 23.4 °C, while T s o i l at 1.7 m ranges from 10 °C to 20.8 °C (Table 2). A peculiar thermal response to the rainfall events is observed at 0.6 and 0.9 m where precipitation falling in winter tends to increase the soil temperature. In contrast, an opposite trend can be observed after the precipitation events occurring in summer, producing a decrease in soil temperature at both depths. The thermal influence of infiltrating water into the soil is lost at 1.7 m both in the summer and winter periods, as well. As regards the VWC, an increase in its values is observed after rainfall events between January and February 2022, when the vadose zone is near the saturated conditions. Starting from March, the soil progressively dries until September in all the monitored depths. The rainfall events occurring between March and September 2022 do not produce any increase in VWC, and starting from September 2022 the VWC starts to increase at the shallow soils’ layers (0.6 m) and only in October 2022 the infiltration process involve the soil at 0.9 m. Starting from December, the soil at 0.6 and 0.9 m behaves similarly under the effect of precipitation events with a fast increase of the VWC. Only in January 2023, the three monitored depths reach the complete saturated conditions. VWC at 0.6 m ranges between 0.22 and 0.4 m 3 /m 3 , VWC at 0.9 m ranges between 0.24 and 0.4 m 3 /m 3 , while VWC at 1.7 m ranges from 0.24 to 0.36 m 3 /m 3 (Table 2). Regarding the bulk EC, it mimics the VWC behavior at all the monitored depths, with an increase in its values as the VWC increases, except for the depth 1.7 m in which sporadic infiltration events are marked by a decrease in EC values. EC at 0.6 m ranges between 0.22 and 0.64 mS/cm, EC at 0.9 m ranges between 0.31 and 0.79 mS/cm, while EC at 1.7 m ranges from 0.24 to 1.05 mS/cm (Table 2). Eventually, the GWL reach its maximum values during winter (January, February) and early spring (March), and then a depletion phase can be observed starting from April 2022, with sparse recharge effects on May 2022.
The GWL reaches its minimum values during the summer and stands at about −9.7 m below ground level until the beginning of November. Then an increase in GWL is observed until February 2023, under the influence of the rainfall events at different stages, when the vadose zone is completely saturated. The GWL ranges between −2.5 and −9.7 m below ground level during the monitored period. Basic statistics for the monitored hydrological variables are reported in Table 2.

3.2. Forecasting Scenarios

The implementation, which is split into four distinct scenarios as described earlier, is considered for prediction purposes, ensuring a robust comparison against atmospheric and vadose zone variables at different depths. Each scenario utilizes the optimized Prophet-based model, which includes components for trend, seasonality, and additional regressors chosen from the specific to three virtual boxes to capture the particular influences of the chosen input variables at a given depth. The testing period is set from the start of March 2023 to the end of April 2023, and therefore, the model makes 61 prediction intervals (days) representing the GWL into the future. The forecasted groundwater level is then compared with actual values to calculate performance metrics. The forecasting performance was evaluated through the following reported metrics:
  • Mean Absolute Error (MAE), which is the absolute value of the difference between paired accurate and predicted data;
  • Mean Absolute Percentage Error (MAPE), which is the percentage expression of MAE obtained through the normalization by the real data;
  • Root Mean Squared Error (RMSE) is the average difference between predicted and real data;
  • Correlation is the strength and direction of the linear relationship between predicted and actual values.
The quantitative results for all the scenarios are presented in Table 3.
These independent scenarios offer a comprehensive approach to understanding and predicting groundwater levels under different conditions. We ran a total of four tests, and the results of each scenario will be presented as the virtual boxes defined earlier.

3.2.1. From Atmosphere (Box 1) to Saurated Zone (Box 3)

The regressors used in this scenario are T a i r and P. The resulting plot is observed in Figure 5, and the observed values of MAE, MAPE, RMSE, and correlation are 0.299, 0.101, 0.356, and 0.857. As determined by cross-validation, the best hyperparameters were a changepoint prior scale of 0.05, a changepoint range of 0.9, daily seasonality enabled, a prior scale for precipitation of 0.1, and a prior scale for temperature (air) of 0.01. The seasonality mode was set to ‘multiplicative’ with a seasonality prior scale of 12.0, and both weekly and yearly seasonality enabled. There is a positive correlation between observed and predicted values, and considering the errors, the results are promising.

3.2.2. From Vadose Zone (Box 2) to Saurated Zone (Box 3)

  • Depth 0.6 m: The regressors used in this scenario are T s o i l 0.6 m, VWC 0.6 m and EC 0.6 m. The resulting plot is observed in Figure 6, and the observed values of MAE, MAPE, RMSE, and correlation are 0.255, 0.086, 0.270, and 0.850, where the errors are lesser than the atmospheric variables. The optimal hyperparameters marked a changepoint prior scale of 0.9, a changepoint range of 0.9, a prior scale for electric conductivity of 0.5, a prior scale for soil temperature of 0.000001, a prior scale for water content of 1.0, a seasonality mode set to ‘multiplicative’ a seasonality prior scale of 15.0. The daily, weekly, and yearly seasonality are set to ‘True’. On visually observing, the forecasted values follow a very similar trend to the observed data points.
  • Depth 0.9 m: Further increasing the depth till 0.9 m, we use the regressors as T s o i l 0.9 m, VWC 0.9 m and EC 0.9 m, respectively. The resulting plot is observed in Figure 7, and the observed values of MAE, MAPE, RMSE, and correlation are 0.274, 0.090, 0.374, and 0.833. The overall performance has been depreciated compared to 0.6 m depth; however, the performance is slightly better than the atmospheric variables. The optimal hyperparameters marked a changepoint prior scale of 0.9, a changepoint range of 0.95, a prior scale for electric conductivity of 0.7, a prior scale for soil temperature of 0.000001, a prior scale for water content of 1.0, a seasonality mode set to ‘additive’ a seasonality prior scale of 0.01. The daily, weekly, and yearly seasonality are set to ‘True’.
  • Depth 1.7 m: The final scenario is implemented at 1.7 m, and here we utilize the regressors T s o i l 1.7 m, VWC 1.7 m, and EC 1.7 m. The resulting plot is observed in Figure 8. The observed values of MAPE, MAE, and correlation are 0.189, 0.062, 0.244 and 0.923. The calculated errors are considerably lower and offer the most efficient performance when compared between different depths of the vadose zone. The hyperparameters used to produce the forecast included a changepoint prior scale of 0.7, a changepoint range of 0.95, a prior scale for electric conductivity of 0.000001, a prior scale for soil temperature of 0.5, a prior scale for water content of 0.5, a seasonality mode set to ‘additive’ a seasonality prior scale of 0.01. The daily, weekly, and yearly seasonality are set to ‘True’.

3.2.3. RAPS Method and Trend Analysis

The analysis of the time series of RAPS data for observed GWL vs. predicted ones for the entire dataset (Figure 9) indicates the existence of three sub-periods. The first sub-period is characterized by an upward trend until May 2022, the second sub-period is marked by a downward trend until December 2022, while the last sub-period is characterized again by an upward trend.
The RAPS data for all the analyzed time series follow the same trend, periodicity, and oscillation along the entire timespan. When analyzing the time series of RAPS data of the observed GWL vs. forecasted ones for the testing period (Figure 10) it is possible to observe a general congruence between the obtained graphs with an upward trend until March 2023 and a decreasing trend in all the analyzed time series starting from the second half of March 2023. However, even if marked by the same trends, the periodicity and the oscillation of the RAPS data obtained from the analysis display different results between the analyzed variables. The periodicity and the oscillation of RAPS for the predicted GWL obtained by exploiting Box 2 at 1.7 m (dashed green line) seem to be the most correlated with the observed one (black line). On the contrary, the RAPS for the predicted GWL, obtained by exploiting Box 1 and Box 2 at 0.9 m depth (dashed red line and dashed blue line respectively), display a shift of about 10 days in their peak if compared to the RAPS of the observed GWL. Eventually, the RAPS data for the predicted GWL, obtained by exploiting Box 2 at 0.6 m depth (dashed pink line), deviates more from the RAPS data of observed GWL especially in the first two weeks (i.e., until 15 March 2023).

4. Discussion

The present work shows an AI-based application to GWL prediction focused on the Prophet forecasting method. Interestingly, the best prediction is obtained by using the sensor placed in the vadose zone (Box 2) at 1.7 m. Indeed, as reported in Table 3 this scenario is characterized by maximum correlation values (0.923) and minimum associated errors. On the other hand, the prediction obtained by exploiting the regressors collected from Box 1 (P and T a i r ) and Box 2 (VWC, EC, and T s o i l at 0.6 and 0.9 m) are comparable if considering only the correlation values, and slightly better for the Box 2 regressors’ if considering the associated errors. The obtained results match the observed hydrological processes. In fact, if using Box 1 regressors (P and T a i r ) to forecast the GWL, the effective infiltration processes through the unsaturated layers towards the aquifer are not taken into account. By exploiting the surface data the prediction (even if accurate) is not as accurate as the ones obtained by exploiting Box 2 collected hydrological variables, supporting the fact that not all the precipitation events are correlated by effective recharge processes, due to the evapotranspiration occurring between April and October in our study site. On the other hand the prediction increases by observing the relationship between Box 2 and Box 3, as the probes installed in the vadose zone can catch the effective infiltration patterns occurring in the saturated zone. However, considering Box 2 at 0.6 and 0.9 m two different processes can be observed: (1) the shallower portion of the soil is still influenced by the evapotranspiration effects, with a decreasing influence as the depth increases; (2) T s o i l acts as a nonlinear regressor during the hydrologic year (i.e., increasing after rainfall events in winter and decreasing with the same magnitude under the effect of similar rainfall events occurring in summer). This aspect supports the prior scale for T s o i l automatically obtained by the forecasting method equal to 0.000001 proving that the T s o i l at 0.6 and 0.9 m impact on GWL prediction is lower, or in other terms, it generates noise. The best prediction is obtained by using all the regressors (VWC, EC, and T s o i l ) of Box 2 at 1.7 m. At this depth, the T s o i l variations under the effect of single rainfall events are not caught by the TDR due to the fact that the temperature is recognized as a non-conservative tracer and the T s o i l variations are smoothed during the infiltration processes [62,63]. At this depth, only the seasonal thermal effect can be observed and those are much more correlated to the seasonal GWL oscillation. In fact, the last scenario (Box 2 to Box 3 at 1.7 m) is characterized by a prior scale for T s o i l automatically obtained by the forecasting method equal to 0.5, higher than the ones obtained for the other depths, strengthening the fact that the seasonal GWL oscillation is linked to the seasonal temperature of the soil recorded at 1.7 m. On the contrary, the VWC at 1.7 m impact on the prediction is lower (0.5) with respect to the other depths because the VWC variations at lower depths are slow, and once the saturated conditions are reached the VWC still remains high and constant for a prolonged period. The opposite behavior is observed during summer when the VWC still remains low for a prolonged period and no sudden variations of VWC are recorded under the effect of rainfall events. The EC affects the prediction at all depths in a non-conservative way. This outcome can be related to the heterogeneous dissolution processes occurring in the soil at different depths when the percolation occurs. In fact, when the water infiltrates into the soil it can produce an effective solute transport whose magnitude is related to the rainfall intensity and duration, the soil’s granulometric size, and the bio-geochemical processes occurring during the monitored period, resulting in a non-linear behavior of the EC under the effect of each single precipitation event at different depth.
Generally, in our study site, the water infiltrating into the soil produces an increase in EC at all the monitored depths except at 1.7 m where the EC sometimes decreases after the precipitation. Eventually, we exploited two “virtual boxes” (Box 1 and Box 2) with their collected hydrological variables fetched as a time series and used to predict GWL. The findings underscore the significance of hydrological variables from the vadose zone, revealing their capacity to enhance prediction accuracy when compared to the thermo-pluviometric variables. This observation highlights the added value of vadose zone variables in refining GWL forecasting. Our analysis indicates that the optimum performance within the vadose zone is achieved at a depth of 1.7 m, emphasizing the critical influence of this specific depth on the predictive accuracy of GWL models in our study site. The RAPS method applied to the entire dataset (Figure 9) demonstrates the efficiency of the exploited forecasted method highlighting a general correlation both in trends, periodicities, and oscillations along the training, validation, and testing periods. If focusing on the testing dataset (Figure 10) the best correlation of RAPS data for observed GWL is obtained with the TDR placed at 1.7 m depth into the vadose zone, highlighting and strengthening how this specific depth affects the predictive accuracy of the GWL in our case study. While this research has provided valuable insights by concentrating on a singular model and employing grid search optimization, it is essential to acknowledge its limitations. Future remarks may lay in a more comprehensive understanding of the forecasting that can be achieved by exploring different approaches to hyperparameter tuning. This implies that the study could be extended to investigate alternative methods beyond grid search, such as Bayesian optimization [64], to identify the most effective configuration for the Prophet model.
Future developments may lie in comparing the obtained results with other state-of-the-art time series models like ARIMA and SARIMAX, and obtain benchmarks for the selected site [42]. Additionally, incorporating emerging deep learning techniques like NeuralProphet or transformer-based models would contribute to a more comprehensive evaluation [65,66,67]. These models have garnered global attention for their ability to capture complex temporal dependencies and patterns, potentially offering alternative solutions to time series forecasting challenges. The suggested benchmarking across various models may enhance a broader perspective on the strengths and weaknesses of different approaches. Another limitation of the study lies in the duration of the monitoring period. Indeed, since the exploited time series is 16 months long, they cannot catch year seasonality on the observed hydrological processes. This aspect may infer the forecasting process and limit the forecasting time period. Further development of the study can be done only if two or more complete hydrologic years of data are collected. Eventually, a further approach may involve the use of effective precipitation as a regressor of Box 1 by excluding the contribution of evapotranspiration, which can be calculated at a daily scale exploiting the Hargreaves and Samani method [68], starting from the collected P and T a i r data. However, it has to be pointed out that at this stage, the characterization of the aquifer recharge processes occurring in the study area is out of the scope of the present research, which aims to improve and explore the Prophet forecasting method by using hydrological variables collected in the vadose zone.

5. Conclusions

The implemented high-resolution monitoring system coupled with the Prophet forecasting method, allows for the early detection of changes and trends in groundwater level, offering a proactive approach to managing this valuable resource. This capability is vital in promptly identifying and addressing declining water levels and groundwater availability, facilitating more effective and sustainable groundwater resource management practices. Indeed, water management companies can proactively respond by strategically placing water pumps at increased depths within wells or by supplementing the aqueduct system through integrated or emergency water plants. Besides, in this study, a high-resolution hydrogeological monitoring system has been demonstrated imperative because it provided precise and quantitative data on crucial hydrological processes. Such monitoring offers unparalleled accuracy in predicting groundwater levels. The ultimate goal of this study is to highlight the importance of combining continuous field monitoring of the unsaturated zone to understand the recharge mechanisms occurring in heterogeneous media in a continuously changing climatic context characterized by an increase in hydrogeological extremes.

Author Contributions

Conceptualization, D.F.; methodology, D.F., G.N. and A.G.; software, G.N. and A.G.; validation, D.F., A.G., A.P., A.M. and A.T.; formal analysis, D.F., G.N. and A.G.; investigation, D.F. and A.P.; resources, D.F. and A.T.; data curation, D.F.; writing—original draft preparation, D.F., G.N. and A.G.; writing—review and editing, D.F., A.G., A.P., A.M. and A.T.; visualization, D.F. and G.N.; supervision, A.M. and A.T.; project administration, A.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available within the article.

Acknowledgments

The authors kindly acknowledge Stefano Palpacelli for the effort made to support the implementation of the field monitoring site, Consorzio Gorgovivo Azienda Speciale, and Viva Servizi S.p.A. for permission to occupy part of their well-field.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ayers, R.S.; Westcot, D.W. Water Quality for Agriculture; Food and Agriculture Organization of the United Nations: Rome, Italy, 1985; Volume 29. [Google Scholar]
  2. Kløve, B.; Ala-Aho, P.; Bertrand, G.; Boukalova, Z.; Ertürk, A.; Goldscheider, N.; Ilmonen, J.; Karakaya, N.; Kupfersberger, H.; Kvœrner, J.; et al. Groundwater dependent ecosystems. Part I: Hydroecological status and trends. Environ. Sci. Policy 2011, 14, 770–781. [Google Scholar] [CrossRef]
  3. Klemeš, J.J. Industrial water recycle/reuse. Curr. Opin. Chem. Eng. 2012, 1, 238–245. [Google Scholar] [CrossRef]
  4. Vorosmarty, C.J.; Green, P.; Salisbury, J.; Lammers, R.B. Global water resources: Vulnerability from climate change and population growth. Science 2000, 289, 284–288. [Google Scholar] [CrossRef] [PubMed]
  5. Arnell, N.W. Climate change and global water resources. Glob. Environ. Chang. 1999, 9, S31–S49. [Google Scholar] [CrossRef]
  6. Delgado, J.A.; Groffman, P.M.; Nearing, M.A.; Goddard, T.; Reicosky, D.; Lal, R.; Kitchen, N.R.; Rice, C.W.; Towery, D.; Salon, P. Conservation practices to mitigate and adapt to climate change. J. Soil Water Conserv. 2011, 66, 118A–129A. [Google Scholar] [CrossRef]
  7. Maier, H.R.; Dandy, G.C. Neural networks for the prediction and forecasting of water resources variables: A review of modelling issues and applications. Environ. Model. Softw. 2000, 15, 101–124. [Google Scholar] [CrossRef]
  8. Dahan, O. Vadose zone monitoring as a key to groundwater protection. Front. Water 2020, 2, 61. [Google Scholar] [CrossRef]
  9. Poulain, A.; Watlet, A.; Kaufmann, O.; Van Camp, M.; Jourde, H.; Mazzilli, N.; Rochez, G.; Deleu, R.; Quinif, Y.; Hallet, V. Assessment of groundwater recharge processes through karst vadose zone by cave percolation monitoring. Hydrol. Process. 2018, 32, 2069–2083. [Google Scholar] [CrossRef]
  10. Singh, G.; Kaur, G.; Williard, K.; Schoonover, J.; Kang, J. Monitoring of water and solute transport in the vadose zone: A review. Vadose Zone J. 2018, 17, 1–23. [Google Scholar] [CrossRef]
  11. Harter, T.; Hopmans, J.W.; Feddes, R. Role of Vadose Zone Flow Processes in Regional Scale Hydrology: Review, Opportunities and Challenges; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2004; Volume 6, pp. 179–208. [Google Scholar]
  12. Seiler, K.P.; Gat, J.R. Groundwater Recharge from Run-Off, Infiltration and Percolation; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2007; Volume 55. [Google Scholar]
  13. Nogueira, G.E.; Gonçalves, R.D. Groundwater recharge in phreatic aquifers, a case study: Modeling unsaturated zone and recharge rates of the Rio Claro Aquifer using Hydrus-1D. Holos Environ. 2021, 21, 402–422. [Google Scholar] [CrossRef]
  14. Cassiani, G.; Binley, A.; Ferré, T.P. Unsaturated zone processes. In Applied Hydrogeophysics; Springer: Berlin/Heidelberg, Germany, 2006; pp. 75–116. [Google Scholar]
  15. Lin, H. Earth’s Critical Zone and hydropedology: Concepts, characteristics, and advances. Hydrol. Earth Syst. Sci. 2010, 14, 25–45. [Google Scholar] [CrossRef]
  16. Goss, M.J.; Ehlers, W.; Unc, A. The role of lysimeters in the development of our understanding of processes in the vadose zone relevant to contamination of groundwater aquifers. Phys. Chem. Earth Parts A/B/C 2010, 35, 913–926. [Google Scholar] [CrossRef]
  17. Tanner, J.L. Understanding and Modelling of Surface and Groundwater Interactions. Ph.D. Thesis, Rhodes University, Rhodes, Greece, 2013. [Google Scholar]
  18. Calderwood, A.J.; Pauloo, R.A.; Yoder, A.M.; Fogg, G.E. Low-cost, open source wireless sensor network for real-time, scalable groundwater monitoring. Water 2020, 12, 1066. [Google Scholar] [CrossRef]
  19. Barzegar, M.; Blanks, S.; Gharehdash, S.; Timms, W. Development of IOT-based low-cost MEMS pressure sensor for groundwater level monitoring. Meas. Sci. Technol. 2023, 34, 115103. [Google Scholar] [CrossRef]
  20. Mammoliti, E.; Fronzi, D.; Mancini, A.; Valigi, D.; Tazioli, A. WaterbalANce, a WebApp for Thornthwaite–Mather Water Balance Computation: Comparison of applications in two European watersheds. Hydrology 2021, 8, 34. [Google Scholar] [CrossRef]
  21. Donker, N. WTRBLN: A computer program to calculate water balance. Comput. Geosci. 1987, 13, 95–122. [Google Scholar] [CrossRef]
  22. Mammoliti, E.; Fronzi, D.; Palpacelli, S.; Biagiola, N.; Tazioli, A. Assessment of urban landslide groundwater characteristics and origin using artificial tracers, hydro-chemical and stable isotope approaches. Environ. Earth Sci. 2023, 82, 211. [Google Scholar] [CrossRef]
  23. Mao, W.; Yang, J.; Zhu, Y.; Ye, M.; Liu, Z.; Wu, J. An efficient soil water balance model based on hybrid numerical and statistical methods. J. Hydrol. 2018, 559, 721–735. [Google Scholar] [CrossRef]
  24. Loliyana, V.D.; Patel, P.L. A physics based distributed integrated hydrological model in prediction of water balance of a semi-arid catchment in India. Environ. Model. Softw. 2020, 127, 104677. [Google Scholar] [CrossRef]
  25. Berthelin, R.; Olarinoye, T.; Rinderer, M.; Mudarra, M.; Demand, D.; Scheller, M.; Hartmann, A. Estimating karst groundwater recharge from soil moisture observations–a new method tested at the Swabian Alb, southwest Germany. Hydrol. Earth Syst. Sci. 2023, 27, 385–400. [Google Scholar] [CrossRef]
  26. Fronzi, D.; Di Curzio, D.; Rusi, S.; Valigi, D.; Tazioli, A. Comparison between periodic tracer tests and time-series analysis to assess mid-and long-term recharge model changes due to multiple strong seismic events in carbonate aquifers. Water 2020, 12, 3073. [Google Scholar] [CrossRef]
  27. Lary, D.J.; Alavi, A.H.; Gandomi, A.H.; Walker, A.L. Machine learning in geosciences and remote sensing. Geosci. Front. 2016, 7, 3–10. [Google Scholar] [CrossRef]
  28. Mammoliti, E.; Di Stefano, F.; Fronzi, D.; Mancini, A.; Malinverni, E.S.; Tazioli, A. A machine learning approach to extract rock mass discontinuity orientation and spacing, from laser scanner point clouds. Remote Sens. 2022, 14, 2365. [Google Scholar] [CrossRef]
  29. Maniar, H.; Ryali, S.; Kulkarni, M.S.; Abubakar, A. Machine-learning methods in geoscience. In Proceedings of the SEG International Exposition and Annual Meeting, SEG, Tokyo, Japan, 12–14 November 2018; p. SEG-2018. [Google Scholar]
  30. Karpatne, A.; Ebert-Uphoff, I.; Ravela, S.; Babaie, H.A.; Kumar, V. Machine learning for the geosciences: Challenges and opportunities. IEEE Trans. Knowl. Data Eng. 2018, 31, 1544–1554. [Google Scholar] [CrossRef]
  31. Dramsch, J.S. 70 years of machine learning in geoscience in review. Adv. Geophys. 2020, 61, 1–55. [Google Scholar]
  32. Shirmohammadi, B.; Vafakhah, M.; Moosavi, V.; Moghaddamnia, A. Application of several data-driven techniques for predicting groundwater level. Water Resour. Manag. 2013, 27, 419–432. [Google Scholar] [CrossRef]
  33. Sarma, R.; Singh, S. A comparative study of data-driven models for groundwater level forecasting. Water Resour. Manag. 2022, 36, 2741–2756. [Google Scholar] [CrossRef]
  34. Amiri, S.; Rajabi, A.; Shabanlou, S.; Yosefvand, F.; Izadbakhsh, M.A. Prediction of groundwater level variations using deep learning methods and GMS numerical model. Earth Sci. Inform. 2023, 1–15. [Google Scholar] [CrossRef]
  35. Rajaee, T.; Ebrahimi, H.; Nourani, V. A review of the artificial intelligence methods in groundwater level modeling. J. Hydrol. 2019, 572, 336–351. [Google Scholar] [CrossRef]
  36. Vadiati, M.; Rajabi Yami, Z.; Eskandari, E.; Nakhaei, M.; Kisi, O. Application of artificial intelligence models for prediction of groundwater level fluctuations: Case study (Tehran-Karaj alluvial aquifer). Environ. Monit. Assess. 2022, 194, 619. [Google Scholar] [CrossRef]
  37. Sharafati, A.; Asadollah, S.B.H.S.; Neshat, A. A new artificial intelligence strategy for predicting the groundwater level over the Rafsanjan aquifer in Iran. J. Hydrol. 2020, 591, 125468. [Google Scholar] [CrossRef]
  38. Najafabadipour, A.; Kamali, G.; Nezamabadi-pour, H. The Innovative Combination of Time Series Analysis Methods for the Forecasting of Groundwater Fluctuations. Water Resour. 2022, 49, 283–291. [Google Scholar] [CrossRef]
  39. Yan, Z.; Lu, X.; Wu, L. Exploring the Effect of Meteorological Factors on Predicting Hourly Water Levels Based on CEEMDAN and LSTM. Water 2023, 15, 3190. [Google Scholar] [CrossRef]
  40. Dadhich, A.P.; Goyal, R.; Dadhich, P.N. Assessment and prediction of groundwater using geospatial and ANN modeling. Water Resour. Manag. 2021, 35, 2879–2893. [Google Scholar] [CrossRef]
  41. Khozani, Z.S.; Banadkooki, F.B.; Ehteram, M.; Ahmed, A.N.; El-Shafie, A. Combining autoregressive integrated moving average with Long Short-Term Memory neural network and optimisation algorithms for predicting ground water level. J. Clean. Prod. 2022, 348, 131224. [Google Scholar] [CrossRef]
  42. Galdelli, A.; Narang, G.; Migliorelli, L.; Izzo, A.D.; Mancini, A.; Zingaretti, P. An AI-Driven Prototype for Groundwater Level Prediction: Exploring the Gorgovivo Spring Case Study. In Proceedings of the International Conference on Image Analysis and Processing, Udine, Italy, 11–15 September 2023; Springer: Cham, Switzerland, 2023; pp. 418–429. [Google Scholar]
  43. Khan, J.; Lee, E.; Balobaid, A.S.; Kim, K. A Comprehensive Review of Conventional, Machine Leaning, and Deep Learning Models for Groundwater Level (GWL) Forecasting. Appl. Sci. 2023, 13, 2743. [Google Scholar] [CrossRef]
  44. Afrifa, S.; Zhang, T.; Appiahene, P.; Varadarajan, V. Mathematical and Machine Learning Models for Groundwater Level Changes: A Systematic Review and Bibliographic Analysis. Future Internet 2022, 14, 259. [Google Scholar] [CrossRef]
  45. Ren, H.; Cromwell, E.; Kravitz, B.; Chen, X. Using long short-term memory models to fill data gaps in hydrological monitoring networks. Hydrol. Earth Syst. Sci. 2022, 26, 1727–1743. [Google Scholar] [CrossRef]
  46. Taylor, S.J.; Letham, B. Forecasting at scale. PeerJ 2018, 72, 37–45. [Google Scholar] [CrossRef]
  47. Zarinmehr, H.; Tizro, A.T.; Fryar, A.E.; Pour, M.K.; Fasihi, R. Prediction of groundwater level variations based on gravity recovery and climate experiment (GRACE) satellite data and a time-series analysis: A case study in the Lake Urmia basin, Iran. Environ. Earth Sci. 2022, 81, 180. [Google Scholar] [CrossRef]
  48. Aguilera, H.; Guardiola-Albert, C.; Naranjo-Fernández, N.; Kohfahl, C. Towards flexible groundwater-level prediction for adaptive water management: Using Facebook’s Prophet forecasting approach. Hydrol. Sci. J. 2019, 64, 1504–1518. [Google Scholar] [CrossRef]
  49. Mussi, M.; Nanni, T.; Tazioli, A.; Vivalda, P.M. The Mt Conero limestone ridge: The contribution of stable isotopes to the identification of the recharge area of aquifers. Ital. J. Geosci. 2017, 136, 186–197. [Google Scholar] [CrossRef]
  50. Fronzi, D.; Gaiolini, M.; Mammoliti, E.; Colombani, N.; Palpacelli, S.; Marcellini, M.; Tazioli, A. Groundwater-surface water interaction revealed by meteorological trends and groundwater fluctuations on stream water level. Acque Sotter.-Ital. J. Groundw. 2022, 11, 19–28. [Google Scholar] [CrossRef]
  51. Aquilanti, L.; Clementi, F.; Nanni, T.; Palpacelli, S.; Tazioli, A.; Vivalda, P.M. DNA and fluorescein tracer tests to study the recharge, groundwater flowpath and hydraulic contact of aquifers in the Umbria-Marche limestone ridge (central Apennines, Italy). Environ. Earth Sci. 2016, 75, 1–17. [Google Scholar] [CrossRef]
  52. Tassetti, A.N.; Galdelli, A.; Pulcinella, J.; Mancini, A.; Bolognini, L. Addressing Gaps in Small-Scale Fisheries: A Low-Cost Tracking System. Sensors 2022, 22, 839. [Google Scholar] [CrossRef] [PubMed]
  53. Galdelli, A.; Mancini, A.; Tassetti, A.N.; Ferrà Vega, C.; Armelloni, E.; Scarcella, G.; Fabi, G.; Zingaretti, P. A Cloud Computing Architecture to Map Trawling Activities Using Positioning Data. In Proceedings of the 15th IEEE/ASME International Conference on Mechatronic and Embedded Systems and Applications, Beijing, China, 12–15 October 2019; Volume 9, p. V009T12A035. [Google Scholar]
  54. Galdelli, A.; Mancini, A.; Frontoni, E.; Tassetti, A.N. A Feature Encoding Approach and a Cloud Computing Architecture to Map Fishing Activities. In Proceedings of the 17th IEEE/ASME International Conference on Mechatronic and Embedded Systems and Applications (MESA), Virtual, 17–19 August 2021; Volume 7, p. V007T07A003. [Google Scholar]
  55. Ho, S.; Xie, M. The use of ARIMA models for reliability forecasting and analysis. Comput. Ind. Eng. 1998, 35, 213–216. [Google Scholar] [CrossRef]
  56. Šrajbek, M.; Đurin, B.; Sušilović, P.; Singh, S.K. Application of the RAPS Method for Determining the Dependence of Nitrate Concentration in Groundwater on the Amount of Precipitation. Earth 2023, 4, 266–277. [Google Scholar] [CrossRef]
  57. Fiorillo, F.; Petitta, M.; Preziosi, E.; Rusi, S.; Esposito, L.; Tallini, M. Long-term trend and fluctuations of karst spring discharge in a Mediterranean area (central-southern Italy). Environ. Earth Sci. 2015, 74, 153–172. [Google Scholar] [CrossRef]
  58. Garbrecht, J.; Fernandez, G.P. Visualization of Trends and Fluctuations in Climatic Records 1. J. Am. Water Resour. Assoc. 1994, 30, 297–306. [Google Scholar] [CrossRef]
  59. Santangelo, M.; Althuwaynee, O.; Alvioli, M.; Ardizzone, F.; Bianchi, C.; Bornaetxea, T.; Brunetti, M.; Bucci, F.; Cardinali, M.; Donnini, M.; et al. Inventory of landslides triggered by an extreme rainfall event in Marche-Umbria, Italy, on 15 September 2022. Sci. Data 2023, 10, 427. [Google Scholar] [CrossRef]
  60. Torcasio, R.C.; Papa, M.; Del Frate, F.; Dietrich, S.; Toffah, F.E.; Federico, S. Study of the Intense Meteorological Event Occurred in September 2022 over the Marche Region with WRF Model: Impact of Lightning Data Assimilation on Rainfall and Lightning Prediction. Atmosphere 2023, 14, 1152. [Google Scholar] [CrossRef]
  61. Morelli, S.; Bonì, R.; Guidi, E.; De Donatis, M.; Pappafico, G.; Francioni, M. L’alluvione delle Marche del 15 settembre 2022, cause e conseguenze. Cult. Territ. Linguaggi 2023, 24, 136–147. [Google Scholar]
  62. Luhmann, A.J.; Covington, M.D.; Alexander, S.C.; Chai, S.Y.; Schwartz, B.F.; Groten, J.T.; Alexander, E.C., Jr. Comparing conservative and nonconservative tracers in karst and using them to estimate flow path geometry. J. Hydrol. 2012, 448, 201–211. [Google Scholar] [CrossRef]
  63. Westhoff, M.; Bogaard, T.; Savenije, H. Quantifying the effect of in-stream rock clasts on the retardation of heat along a stream. Adv. Water Resour. 2010, 33, 1417–1425. [Google Scholar] [CrossRef]
  64. Owolabi, O.O.; Sunter, D.A. Bayesian Optimization and Hierarchical Forecasting of Non-Weather-Related Electric Power Outages. Energies 2022, 15, 1958. [Google Scholar] [CrossRef]
  65. Triebe, O.; Hewamalage, H.; Pilyugina, P.; Laptev, N.; Bergmeir, C.; Rajagopal, R. NeuralProphet: Explainable Forecasting at Scale, 2021.
  66. Mancini, A.; Cosoli, G.; Galdelli, A.; Violini, L.; Pandarese, G.; Mobili, A.; Blasi, E.; Tittarelli, F.; Revel, G.M. A monitoring platform for the built environment: Towards the development of an early warning system in a seismic context. In Proceedings of the 2023 IEEE International Workshop on Metrology for Living Environment (MetroLivEnv), Milano, Italy, 29–31 May 2023; pp. 102–106. [Google Scholar]
  67. Wang, J.; Li, C.; Li, L.; Huang, Z.; Wang, C.; Zhang, H.; Zhang, Z. InSAR time-series deformation forecasting surrounding Salt Lake using deep transformer models. Sci. Total Environ. 2023, 858, 159744. [Google Scholar] [CrossRef]
  68. Hargreaves, G.H.; Samani, Z.A. Reference crop evapotranspiration from temperature. Appl. Eng. Agric. 1985, 1, 96–99. [Google Scholar] [CrossRef]
Figure 1. Location of the study area with a (a) simplified geological map of the analyzed basin and (b) schematic hydrogeologic cross-section. Modified from [50].
Figure 1. Location of the study area with a (a) simplified geological map of the analyzed basin and (b) schematic hydrogeologic cross-section. Modified from [50].
Water 16 00152 g001
Figure 2. Supporting scheme of natural conditions with related collected variables, input for the Prophet method.
Figure 2. Supporting scheme of natural conditions with related collected variables, input for the Prophet method.
Water 16 00152 g002
Figure 3. Rolling cross-validation strategy applied towards hyperparameter tuning.
Figure 3. Rolling cross-validation strategy applied towards hyperparameter tuning.
Water 16 00152 g003
Figure 4. Collected hydrological variables for the monitoring period.
Figure 4. Collected hydrological variables for the monitoring period.
Water 16 00152 g004
Figure 5. Prediction scenario of Box 1 to Box 3 using surface level variables.
Figure 5. Prediction scenario of Box 1 to Box 3 using surface level variables.
Water 16 00152 g005
Figure 6. Prediction scenario of Box 2 to Box 3 at a depth of 0.6 m.
Figure 6. Prediction scenario of Box 2 to Box 3 at a depth of 0.6 m.
Water 16 00152 g006
Figure 7. Prediction scenario of Box 2 to Box 3 at a depth of 0.9 m.
Figure 7. Prediction scenario of Box 2 to Box 3 at a depth of 0.9 m.
Water 16 00152 g007
Figure 8. Prediction scenario of Box 2 to Box 3 at a depth of 1.7 m.
Figure 8. Prediction scenario of Box 2 to Box 3 at a depth of 1.7 m.
Water 16 00152 g008
Figure 9. Time series of RAPS data of observed GWL (black line) and predicted GWL data span of the entire dataset across each forecast scenario, encompassing Box 1 to Box 3 and Box 2 (at 0.6 m, 0.9 m, and 1.7 m depths) to Box 3.
Figure 9. Time series of RAPS data of observed GWL (black line) and predicted GWL data span of the entire dataset across each forecast scenario, encompassing Box 1 to Box 3 and Box 2 (at 0.6 m, 0.9 m, and 1.7 m depths) to Box 3.
Water 16 00152 g009
Figure 10. Time series of RAPS data of observed GWL (black line) and predicted GWL data span of the testing dataset for each predicted scenario, including Box 1 to Box 3 and Box 2 (at 0.6 m, 0.9 m, and 1.7 m depths) to Box 3.
Figure 10. Time series of RAPS data of observed GWL (black line) and predicted GWL data span of the testing dataset for each predicted scenario, including Box 1 to Box 3 and Box 2 (at 0.6 m, 0.9 m, and 1.7 m depths) to Box 3.
Water 16 00152 g010
Table 1. Hyperparameters grid optimized using rolling cross-validation strategy for different scenarios: Atmosphere to Saturated zone (Box 1 to Box 3) and Vadose zone to Saturated zone (Box 2 to Box 3).
Table 1. Hyperparameters grid optimized using rolling cross-validation strategy for different scenarios: Atmosphere to Saturated zone (Box 1 to Box 3) and Vadose zone to Saturated zone (Box 2 to Box 3).
ScenarioHyperparametersValues
Box 1 to Box 3changepoint_prior_scale[0.01, 0.05, 0.1, 0.7, 0.9]
seasonality_mode[’additive’, ’multiplicative’]
seasonality_prior_scale[0.01, 2.0, 12.0, 15.0, 75.0]
prior_scale_temperature[0.000001, 0.01, 0.5, 0.7, 1.0]
prior_scale_rain[0.000001, 0.01, 0.5, 0.7, 1.0]
changepoint_range[0.70, 0.85, 0.90, 0.95]
daily_seasonality[’True’, ’False’]
weekly_seasonality[’True’, ’False’]
yearly_seasonality[’True’, ’False’]
Box 2 to Box 3changepoint_prior_scale[0.01, 0.05, 0.1, 0.7, 0.9]
seasonality_mode[’additive’, ’multiplicative’]
seasonality_prior_scale[0.01, 2.0, 12.0, 15.0, 75.0]
prior_scale_water_content[0.000001, 0.01, 0.5, 0.7, 1.0]
prior_scale_soil_temperature[0.000001, 0.01, 0.5, 0.7, 1.0]
prior_scale_electric_conductivity[0.000001, 0.01, 0.5, 0.7, 1.0]
changepoint_range[0.70, 0.85, 0.90, 0.95]
daily_seasonality[’True’, ’False’]
weekly_seasonality[’True’, ’False’]
yearly_seasonality[’True’, ’False’]
Table 2. Basic statistics of the considered time-series.
Table 2. Basic statistics of the considered time-series.
Box n°Hydrological VariableMeasurement UnitMeanMin25thMedian75thMax
1Pmm/15min0.02000012
T a i r °C14.107.1212.5319.6135
2T s o i l 0.6 m°C14.66.79.712.719.225.5
VWC 0.6 mm3/m30.310.220.270.330.350.4
EC 0.6 mmS/cm0.40.220.290.460.490.64
T s o i l 0.9 m°C14.58.310.112.718.823.4
VWC 0.9 mm3/m30.320.240.280.330.350.4
EC 0.9 mmS/cm0.510.310.420.550.580.79
T s o i l 1.7 m°C14.51011.413.318.120.8
VWC 1.7 mm3/m30.290.240.250.30.320.36
EC 1.7 mmS/cm0.510.240.260.460.651.05
3GWLm−6.39−9.7−9.2−6.42−4.19−2.5
Table 3. Performance metrics of the model in different scenarios.
Table 3. Performance metrics of the model in different scenarios.
ScenarioDepthAdditive RegressorsMAEMAPERMSECorrelation
Box 1 to Box 3SurfaceP, T a i r 0.2990.1010.3560.857
Box 2 to Box 3Depth 0.6 mT s o i l 0.6 m, VWC 0.6 m, EC 0.6 m0.2550.0860.2700.850
Box 2 to Box 3Depth 0.9 mT s o i l 0.9 m, VWC 0.9 m, EC 0.9 m0.2740.0900.3740.833
Box 2 to Box 3Depth 1.7 mT s o i l 1.7 m, VWC 1.7 m, EC 1.7 m0.1890.0620.2440.923
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

Fronzi, D.; Narang, G.; Galdelli, A.; Pepi, A.; Mancini, A.; Tazioli, A. Towards Groundwater-Level Prediction Using Prophet Forecasting Method by Exploiting a High-Resolution Hydrogeological Monitoring System. Water 2024, 16, 152. https://doi.org/10.3390/w16010152

AMA Style

Fronzi D, Narang G, Galdelli A, Pepi A, Mancini A, Tazioli A. Towards Groundwater-Level Prediction Using Prophet Forecasting Method by Exploiting a High-Resolution Hydrogeological Monitoring System. Water. 2024; 16(1):152. https://doi.org/10.3390/w16010152

Chicago/Turabian Style

Fronzi, Davide, Gagan Narang, Alessandro Galdelli, Alessandro Pepi, Adriano Mancini, and Alberto Tazioli. 2024. "Towards Groundwater-Level Prediction Using Prophet Forecasting Method by Exploiting a High-Resolution Hydrogeological Monitoring System" Water 16, no. 1: 152. https://doi.org/10.3390/w16010152

APA Style

Fronzi, D., Narang, G., Galdelli, A., Pepi, A., Mancini, A., & Tazioli, A. (2024). Towards Groundwater-Level Prediction Using Prophet Forecasting Method by Exploiting a High-Resolution Hydrogeological Monitoring System. Water, 16(1), 152. https://doi.org/10.3390/w16010152

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