Next Article in Journal
Long-Term Performance Evaluation and Fouling Characterization of a Full-Scale Brackish Water Reverse Osmosis Desalination Plant
Previous Article in Journal
Influencing Factors for the Growth of Cladophora and Its Cell Damage and Destruction Mechanism: Implication for Prevention and Treatment
Previous Article in Special Issue
The WRF-Driven Grid-Xin’anjiang Model and Its Application in Small and Medium Catchments of China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Enhanced Spatio-Temporal Modeling for Rainfall Forecasting: A High-Resolution Grid Analysis

1
ICAR-Central Research Institute for Jute and Allied Fibres, Kolkata 700121, India
2
ICAR-Central Inland Fishery Research Institute, Barrackpore 700120, India
3
ICAR-Indian Agricultural Statistics Research Institute, New Delhi 110012, India
4
Department of PS & IT, Tamil Nadu Agricultural University, Coimbatore 641003, India
*
Authors to whom correspondence should be addressed.
Water 2024, 16(13), 1891; https://doi.org/10.3390/w16131891
Submission received: 19 May 2024 / Revised: 24 June 2024 / Accepted: 28 June 2024 / Published: 1 July 2024

Abstract

:
Rainfall serves as a lifeline for crop cultivation in many agriculture-dependent countries including India. Being spatio-temporal data, the forecasting of rainfall becomes a more complex and tedious process. Application of conventional time series models and machine learning techniques will not be a suitable choice as they may not adequately account for the complex spatial and temporal dependencies integrated within the data. This demands some data-driven techniques that can handle the intrinsic patterns such as non-linearity, non-stationarity, and non-normality. Space–Time Autoregressive Moving Average (STARMA) models were highly known for its ability to capture both spatial and temporal dependencies, offering a comprehensive framework for analyzing complex datasets. Spatial Weight Matrix (SWM) developed by the STARMA model helps in integrating the spatial effects of the neighboring sites. The study employed a novel dataset consisting of annual rainfall measurements spanning over 50 (1970–2019) years from 119 different locations (grid of 0.25 × 0.25 degree resolution) of West Bengal, a state of India. These extensive datasets were split into testing and training groups that enable the better understanding of the rainfall patterns at a granular level. The study findings demonstrated a notable improvement in forecasting accuracy by the STARMA model that can exhibit promising implications for agricultural management and planning, particularly in regions vulnerable to climate variability.

1. Introduction

Rainfall pattern in India is largely influenced by the monsoon system, which serves as the primary source of precipitation across the country. The agricultural sector, which directly sustains a significant portion of the population, relies heavily on rainfall for irrigation and crop cultivation [1,2]. The timely distribution and adequate quantity of rainfall are crucial factors that directly influence crop yields and the livelihoods of millions of farmers [3]. Inconsistent rainfall can lead to droughts, crop failures, and economic distress in rural areas [4,5]. Fluctuations in rainfall are becoming increasingly common due to current environmental changes influenced by factors such as climate change, land use alterations, urbanization, and natural variability [6,7]. These fluctuations complicate the forecasting process, leading to economic uncertainties and social challenges. Consequently, accurate rainfall forecasting has emerged as a pivotal scientific challenge, crucial for effective resource management and disaster preparedness across the country [8,9,10]. Accurate rainfall forecasts play a critical role in agricultural management, providing valuable insights into optimal planting, irrigation, and harvesting schedules for farmers. By aligning farming practices with forecasted rainfall patterns, farmers can optimize crop yields, conserve water resources, and mitigate risks associated with extreme weather conditions like droughts or excessive rainfall [11]. Therefore, reliable rainfall forecasting is essential for planning and economic stability across various sectors.
Rainfall data exhibit significant spatial and temporal variability, which can vary greatly over short distances and time periods [12,13]. This inherent variability, shaped by complex climatic interactions, presents formidable challenges for conventional forecasting approaches. This variability requires models that can capture local and regional nuances to provide accurate predictions. Traditional methods often struggle with accuracy and reliability due to these complexities [14].
Autoregressive Integrated Moving Average (ARIMA) models are conventional linear time series models extensively utilized for predicting hydrological and meteorological phenomena over the past few decades [15,16,17]. However, the ARIMA model operates under constraints that limit its effectiveness in complex environmental applications [18]. Primarily, it is a univariate model that does not account for spatial interdependencies and is inadequate for handling the non-linear and spatially interrelated data [19,20]. Machine learning techniques, such as Artificial Neural Networks (ANNs), Support Vector Machines (SVMs), and Random Forest (RF), have gained prominence in handling complex time series data [21,22,23]. These methods offer substantial improvements in modeling capabilities due to their data-driven nature, which does not necessarily require prior assumptions about the data’s distribution or structure [24]. However, these models often do not account for spatial interactions, which are crucial for accurately forecasting rainfall distribution across different regions [25,26,27].
Given the limitations of existing models and the complexity of rainfall data, which embody temporal, spatial, linear, and non-linear dynamics, Space–Time Autoregressive Moving Average (STARMA) models have been introduced as new alternatives. STARMA models incorporate both lagged effects and neighborhood influences of adjacent regions, offering an improved framework over univariate conventional and machine learning models [28,29]. The development of the spatial weight matrix is key to building these models. This approach promises to be a pivotal tool in enhancing agricultural planning, water management, and disaster preparedness strategies across climatically diverse and sensitive regions like India through efficient rainfall forecasting.

2. Background

Conventional time series models and machine learning models are widely used for univariate time series data, with extensive applications in environmental and weather modeling. Conventional models, such as ARIMA, typically address linear patterns, while machine learning models are adept at handling non-linearities within datasets. Basha et al. combined the ARIMA model with the Generalized Autoregressive Conditional Heteroskedasticity (GARCH) model to predict the summer monsoon rainfall of India [30]. Dotse et al. conducted a comprehensive review on hybrid machine learning models used to enhance rainfall prediction accuracy [9]. Garai et al. explored various machine learning techniques, including Neural Networks (NNs), Support Vector Regression (SVR), and kernel ridge regression, for rainfall forecasting [31]. Sekhar et al. employed both ARIMA and ANN models to forecast rainfall in the coastal Andhra region [32]. Bommala et al. provided an overview of ensemble-based models integrating machine learning and deep learning techniques for rainfall prediction [33]. The integration of Autoregressive (AR) models with Neural Networks (NNs) was explored to enhance monthly rainfall forecasts in the watershed regions of Cuba [34]. This study involved the application of Multi-Layer Perceptron (MLP), Convolutional Neural Network (CNN), Long Short-Term Memory (LSTM) neural networks, ARIMA models, and a hybrid model (ANN + ARIMA) for rainfall prediction. Kaur et al. offered a comprehensive review on the use of AR models in forecasting environmental variables [35]. Zhang et al. explored the forecasting capabilities of ARIMA, Wavelet Neural Network (WNN), and Support Vector Machine (SVM) for drought forecasting [36]. A comparison of various time series models, such as LSTM, Stacked-LSTM, Bidirectional-LSTM Networks, XGBoost, and an ensemble consisting of Gradient Boosting Regressor, Linear Support Vector Regression, and Extra-Trees Regressor, was conducted to forecast hourly rainfall volumes in five major cities of the United Kingdom [37]. Ridwan et al. demonstrated the superiority of Boosted Decision Tree Regression (BDTR) over other models in predicting rainfall in the Terengganu region of Malaysia [38]. Praveen et al. utilized non-parametric models and ANN-MLP to analyze trends and forecast rainfall changes in India [39]. However, these models often fall short when dealing with spatio-temporal data. Handling spatio-temporal data poses significant challenges for conventional and machine learning models due to the inherent complexity of patterns within the data, necessitating the need for alternative models.
The Space–Time Autoregressive Moving Average (STARMA) model stands as a well-established methodology specifically designed for handling such data. Various studies have showcased the efficacy of the STARMA model in effectively managing climatological data, which inherently exhibits spatio-temporal characteristics. Pfeifer and Deutrch proposed a three-stage iterative procedure known as STARIMA models, which integrates ARIMA modeling in the absence of spatial correlation [40]. Subsequently, they extended their work to include Seasonal Space–Time ARIMA modeling [41], along with providing techniques for the identification and interpretation of first-order space–time ARMA models [42]. Pfeifer and Bodily applied space–time ARMA modeling to obtain data from hotels in eight different locations in a major US city [43]. The application of STARMA models extends to comparing monthly mean temperatures from nine meteorological stations across the United Kingdom [44] and studying timber prices [45]. Saha et al. utilized a hybrid spatio-temporal modeling approach by integrating STARMA with ANN and SVM to predict annual precipitation data across six districts in the northern part of West Bengal, India [46]. Rathod et al. applied the Two-Stage Spatio-temporal Time Series model to forecast rice yield [29]. Kumar et al. employed a Hybrid Space–Time modeling approach to forecast the monthly temperatures [47]. Rathod et al. introduced an improved STARMA model for modeling and forecasting the spatio-temporal time series data [48]. Saha et al. employed fuzzy rule-based weighted Space–Time Autoregressive Moving Average models for temperature forecasting [28]. Additionally, the performance of a STARIMA model was compared with a Seasonal ARIMA (SARIMA) model for predicting infectious diseases [49].

3. Materials and Methods

3.1. Data Description

Daily rainfall data of West Bengal, India, were obtained from India Meteorological Department (IMD) web portal (https://www.imdpune.gov.in/lrfindex.php, accessed on 5 December 2022). The dataset consisted of yearly gridded rainfall data of a high-resolution grid (0.25° × 0.25°), encompassing 119 grid points across 23 districts of West Bengal over a period of 50 years (1970–2019) [50]. For each grid, daily data have been converted to annual data using standard procedure. Hence, at each grid point, the spatio-temporal data series included 50 observations. Out of these, data from 45 years were utilized to develop the model, while the remaining observations were reserved for validation purposes. Figure 1 presents the location map of the study region along with the grid points. After calculating the observed and predicted values for each point, these values were interpolated using the Inverse Distance Weighted (IDW) technique to generate comprehensive rainfall maps for the entire state of West Bengal.

3.2. Spatio-Temporal Data

Spatio-temporal time series data are the observations that are recorded with systematic dependencies across different spatial locations over various temporal intervals. Unlike traditional data, the analysis of spatio-temporal data provides comprehensive insights into the subject matter due to its rich and multi-dimensional nature. However, handling such data presents several formidable challenges. These datasets often comprise a massive volume of observations collected across numerous spatial locations and time intervals, necessitating efficient preprocessing and analysis techniques to manage the sheer scale of the data. Furthermore, spatio-temporal datasets exhibit intricate relationships and dependencies across space and time, requiring sophisticated modeling approaches capable of capturing and interpreting these complex interactions accurately. Addressing these challenges is essential for unlocking valuable insights and deriving meaningful information across the domains.

3.3. ARIMA Model

Autoregressive Integrated Moving Average (ARIMA) is a widely used conventional linear time series model [51,52]. There are three parameters of the ARIMA model: autoregressive (p) refers to the information from the past values of the series, order of differencing (d) refers to the number of times the data are differenced to convert them into stationarity, and moving average (q) refers to the information from the past forecast errors. The general equation and working procedure of the ARIMA model are given below.
y t = 1   y t 1   + 2   y t 2   + + p   y t p + ε t θ 1   ε t 1   θ 2   ε t 2   θ q   ε t q
where y t is the time series, i and θ j are model parameters, ε t is a random error, p is the number of lags of autoregressive terms, and q is the number of lags of moving terms.
The stationarity of the data can be checked either by ADF test or graphically observed by using the plots of Autocorrelation Function (ACF) and Partial Autocorrelation Function (PACF). If the time series is not stationary, it can be transformed by appropriately differencing the series. The number of lags of autoregression (p) and moving average (q) can be determined by examining the plots of PACF and ACF [53,54]. If the ACF plot shows a significant autocorrelation at lag k1, and the PACF plot shows a significant partial autocorrelation at lag k2; then, the order of the AR component will be p = k2, and the order of the MA component will be q = k1, respectively. After finding the values of p, d, and q, the ARIMA model is fitted by estimating the parameters of each component using Maximum Likelihood Estimation (MLE). The ARIMA model with the lowest Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) values is usually selected. After fitting a tentative ARIMA model, the adequacy of the model is checked. Confirmation of white noise is performed by examining the residual of the fitted model using Ljung–Box test, which examine the presence of autocorrelation in the residuals. If the model is not white noise, then the procedure is repeated until an adequate model is fitted.

3.4. STARMA Model

The STARMA (Space–Time Autoregressive Moving Average) model is a statistical model used for analyzing and forecasting spatio-temporal data. It is an extension of the traditional ARMA (Autoregressive Moving Average) model, adapted to handle univariate time series data across different locations considering both spatial and temporal dependencies [22,39,42]. It considers how observations at different spatial locations influence each other over time and how observations at the same location evolve over time. The model works by estimating parameters that describe the spatial and temporal dependencies within the data. These parameters include coefficients for autoregressive and moving average terms, as well as parameters that describe the spatial structure of the data. Spatial weight matrix with spatial weights helps in considering the neighboring effect. In addition, first, second, third, or higher order may be chosen depending on the neighboring areas. The STARMA model [35,37] is expressed as follows:
z t = k = 1 p l = 0 λ k Φ k l W l z t k + ϵ t k = 1 q l = 0 m k Θ k l W l ε t k
where zi(t) is the observation at site i (i = 1, 2,…, N) and time t (t ϵ Z);
  • z(t) is a N × 1 vector of observations at time t (t = 1, …, T);
  • p is the autoregressive order (AR), and q is the MA order;
  • λk represents the spatial order of the kth AR term, and mk represents the spatial order of the kth MA term;
  • ϕkl is the AR parameter at temporal lag k and spatial lag l (scalar), and θkl is the MA parameter at temporal lag k and spatial lag l (scalar);
  • W(l) is the N × N matrix of weights for spatial order l;
  • ϵ(t) is normally distributed random error at time t with
E   [ ϵ ( t ) ] = 0   and   E t t + s = G ,     s = 0 0 ,     s 0
The STARMA (Space–Time Autoregressive Moving Average) model can be simplified into two distinct forms based on the absence of certain components. Firstly, in the absence of the moving average component, the STARMA model transforms into the Space–Time Autoregressive (STAR) model. This configuration captures the autoregressive relationships between observations at different spatial locations over various time intervals. It delineates how past observations at different locations influence the current observation at a given location, considering both spatial and temporal dependencies. Conversely, when the autoregressive term is absent, the STARMA model becomes the Space–Time Moving Average (STMA) model [40]. In this rendition, the model solely focuses on capturing the moving average relationships between observations at different spatial locations over time. It elucidates how past error terms at different locations contribute to the current observation at a specific location, accounting for both spatial and temporal dependencies. The goodness of the fitted model is examined by checking the white noise in the residual. Multivariate Box–Pierce test is used for checking the autocorrelation of residuals, i.e., the residuals are normally distributed with mean zero and constant variance [23]. This evaluates the model’s ability in capturing the temporal dependencies in the data.

3.5. Model Evaluation Criterion

A comparative modeling and forecasting performance of the proposed hybrid model with the other existing model was assessed by estimating the Root Mean Square Error (MSE), Mean Absolute Percentage Error (MAPE), and Mean Absolute Error (MAE). RMSE, MAPE, and MAE were calculated using the following formula:
R M S E = t = 1 n Y t Y t ^ Y t 2 n
M A P E = t = 1 n | Y t Y t ^ Y t | 2 n 100
M A E = 1 n t = 1 n Y t Y ^ t
where n represents the total number of modeled or forecasted values, Y t represents the actual value at time t, and Y t ^ represents the corresponding modeled or forecasted value.

4. Results and Discussion

The intensity of the rainfall across the 119 locations of the West Bengal region is depicted in Figure 2, which reveals that the northern region experienced abundant rainfall ranging from 3250 to 3750 mm. Conversely, the southern region exhibited lower rainfall levels, with precipitation ranging between 1500 and 1750 mm. The collected rainfall data were analyzed using two different time series models: ARIMA and STARMA. Both the time series models were compared, and the best model was selected using the error metrics. All the prior statistical assumptions are satisfied before proceeding with the statistical techniques.

4.1. ARIMA Fitting

The study begins with fitting the ARIMA model to the rainfall data. Initially, the Augmented Dickey–Fuller (ADF) test was conducted to check the stationarity of the series. The series exhibiting non-stationarity were differenced to achieve stationarity, while some series were stationary in their base form. The ARIMA models were fitted using parameters estimated through Maximum Likelihood Estimation (MLE). Diagnostic checks using the Ljung–Box test confirmed that all fitted models exhibited white noise, indicating a good fit. The ARCH-LM test also showed the absence of volatility in all models, ensuring that the selected models effectively capture the underlying patterns in the rainfall data. Table 1 lists the ARIMA models fitted for the 119 locations.

4.2. STARMA Fitting

The primary step in fitting the STARMA model involves constructing weight matrices for the data. This study utilized spatial weight matrices up to the first order. The zero-order spatial weight matrix (W(0)) was treated as an identity matrix to incorporate the temporal pattern within the STARMA model. The first-order spatial weight matrix (W(1)) was constructed by considering the percentage sharing of borders between specific districts and their neighbors. This integration of spatial data from neighboring districts into each district is achieved through a row-normalized first-order weight matrix. The zero-order (W(0)) and first-order (W(1)) spatial weight matrices developed for the study location were represented in the Supplementary Tables S1 and S2. The Multivariate Box–Pierce non-correlation test confirmed significant spatio-temporal correlations within the rainfall data, with a chi-squared value of 69,668.64 and a p-value of <0.01 (Table 2), indicating significant spatio-temporal correlation within the rainfall data. Following the assessment of spatio-temporal correlations, the suitable STARMA model was identified, and its parameters were estimated using Maximum Likelihood Estimation methods. The results of the fitted STARMA model are summarized in Table 2, presenting the estimated parameters alongside their standard errors and probability values. The φ of AR and θ of MA components were derived for both spatial lags 0 and 1. All estimated parameters were statistically significant at the 5% level. The equational form of the STARMA model is given by the following:
Z(t) = 0.87 (W(0)) z(t−1) + 0.13 (W(1)) z(t−1) − 0.48 (W(0)) ε(t−1) − 0.23(W(1)) ε(t−1)
A comparison was conducted between the performance of ARIMA and STARMA models in fitting rainfall data across 119 locations in West Bengal. The best-performing model was determined based on error metrics such as RMSE, MAPE, and MAE. Table 3 displays the error metrics calculated from the observed and predicted values for both models. The results revealed that the STARMA model exhibited lower error rates in predicting rainfall across 119 locations compared to the ARIMA model. Five-year (2015–2019) mean observed and predicted rainfall obtained from the STARMA model was mapped, and it can be seen that the STARMA model predicted rainfall is almost at par with the average observed rainfall for the period 2015–2019 (Figure 2). Thus, the STARMA model consistently outperformed the ARIMA model across all error metrics and locations.
Rainfall data always present a multifaceted challenge for modeling due to their intricate and dynamic nature. One primary complexity lies in its temporal dependency, where rainfall patterns exhibit seasonality, trends, and periodic fluctuations over time. These temporal variations necessitate the models to capture and accommodate such dynamic changes effectively [28,48]. Additionally, rainfall data also possess significant spatial variability, with precipitation levels varying widely across different geographical locations. This variability is influenced by diverse factors that can account for these spatial intricacies. To address these challenges, the study used a sophisticated modeling approach that was capable of capturing the complex interplay of temporal and spatial dynamics inherent to rainfall data. The study initially applied the ARIMA model to all 119 locations, and the candidate ARIMA model was chosen for each dataset. Next, the STARMA model that is highly known for its efficacy in handling spatio-temporal datasets was applied. Prior to fitting the STARMA model, the presence of spatial relationships between the selected locations was confirmed, ensuring the suitability of the model for the study. The construction of spatial weight matrices is a critical step that can influence model accuracy. This study used the first-order weight matrices. STARMA (φ10, θ10, φ11, θ11) model was identified as the best-fitted model for the data. Error metrics such RMSE, MAPE, and MAE were used to compare the performance of the fitted mode. The results show that the univariate ARIMA models fall short in capturing the spatial dependencies. In contrast, the STARMA model that is specifically designed to handle spatio-temporal datasets was found to be more suitable for this study. The STARMA model provides a robust framework for modeling and predicting rainfall in West Bengal. Its superior performance was compared to the ARIMA model, which highlights the importance of considering spatio-temporal dependencies in rainfall data [28,29]. Thus, the STARMA model outperformed all ARIMA models for all study locations, i.e., the grid points. Future research could further explore the application of STARMA models in different regions and climates to validate and extend these findings. Integrating variables like temperature and humidity could enhance predictive power, providing a more comprehensive understanding of rainfall patterns. This has significant implications for agriculture, water management, and disaster preparedness.

5. Conclusions

Every time series model has its own limitations and strong assumptions. Some models excel in handling linear datasets, while others are better suited for non-linear datasets. Therefore, it is crucial to carefully consider the nature of the data under study before selecting a time series model. This study serves as a prime example where the primary objective is to model the rainfall patterns across 119 different locations in West Bengal, India. Rainfall data are highly known for their spatio-temporal variability, which requires some data-driven time series models rather than conventional time series models. This study employed ARIMA and STARMA models where the performance of the fitted model was compared using error metrics such as RMSE, MAPE, and MAE. The results showed that the fitted STARMA (φ10, θ10, φ11, θ11) model consistently yielded lower errors compared to ARIMA across all 119 locations. This provides strong evidence that STARMA models perform well in handling spatio-temporal data by capturing the complex temporal and spatial dynamics inherent in the rainfall data across diverse geographical locations. This helps in efficiently forecasting the rainfall, which ensures the efficient allocation of water resources and the minimization of the potential risks to agriculture. Future studies could enhance forecasting precision by developing hybrid space–time hybrid models integrated with machine learning and deep learning techniques.

Supplementary Materials

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

Author Contributions

Conceptualization, N.M.A.; analysis, N.M.A., M.R. and C.J., methodology N.M.A., M.R. and C.J.; data curation, C.J. and S.V.S., supervision, S.M., S.K.P. and G.K., resources, S.V.S.; writing—review and editing, S.G., S.P.M. and R.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data will be made available on reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Pandey, D. Agricultural Sustainability and Climate Change Nexus. In Contemporary Environmental Issues and Challenges in Era of Climate Change; Springer: Singapore, 2020; pp. 77–97. [Google Scholar]
  2. Rao, C.S.; Lal, R.; Prasad, J.V.; Gopinath, K.A.; Singh, R.; Jakkula, V.S.; Sahrawat, K.L.; Venkateswarlu, B.; Sikka, A.K.; Virmani, S.M. Potential and challenges of rainfed farming in India. Adv. Agron. 2015, 133, 113–181. [Google Scholar]
  3. Maharana, P.; Dimri, A.P. The Indian Monsoon: Past, present and future. Proc. Indian Natl. Sci. Acad. 2019, 85, 403–420. [Google Scholar]
  4. Singh, R.B. Environmental consequences of agricultural development: A case study from the Green Revolution state of Haryana, India. Agric. Ecosyst. Environ. 2000, 82, 97–103. [Google Scholar] [CrossRef]
  5. Haile, M. Weather patterns, food security and humanitarian response in sub-Saharan Africa. Philos. Trans. R. Soc. B Biol. Sci. 2005, 360, 2169–2182. [Google Scholar] [CrossRef] [PubMed]
  6. Portela, M.M.; Espinosa, L.A.; Zelenakova, M. Long-term rainfall trends and their variability in mainland Portugal in the last 106 years. Climate 2020, 8, 146. [Google Scholar] [CrossRef]
  7. Turner, A.G.; Annamalai, H. Climate change and the South Asian summer monsoon. Nat. Clim. Change 2012, 2, 587–595. [Google Scholar] [CrossRef]
  8. Mondini, A.C.; Guzzetti, F.; Melillo, M. Deep learning forecast of rainfall-induced shallow landslides. Nat. Commun. 2023, 14, 2466. [Google Scholar] [CrossRef]
  9. Dotse, S.Q.; Larbi, I.; Limantol, A.M.; DeSilva, L.C. A review of the application of hybrid machine learning models to improve rainfall prediction. Model. Earth Syst. Environ. 2024, 10, 19–44. [Google Scholar] [CrossRef]
  10. Allawi, M.F.; Abdulhameed, U.H.; Adham, A.; Sayl, K.N.; Sulaiman, S.O.; Ramal, M.M.; Sherif, M.; El-Shafie, A. Monthly rainfall forecasting modelling based on advanced machine learning methods: Tropical region as case study. Eng. Appl. Comput. Fluid Mech. 2023, 17, 2243090. [Google Scholar] [CrossRef]
  11. Palanisami, K.; Meinzen-Dick, R.; Giordano, M. Climate change and water supplies: Options for sustaining tank irrigation potential in India. Econ. Political Wkly. 2010, 45, 183–190. [Google Scholar]
  12. Sahu, R.T.; Verma, M.K.; Ahmad, I. Impact of long-distance interaction indicator (monsoon indices) on spatio-temporal variability of precipitation over the Mahanadi River basin. Water Resour. Res. 2023, 59, e2022WR033805. [Google Scholar] [CrossRef]
  13. Tariq, A.; Qin, S. Spatio-temporal variation in surface water in Punjab, Pakistan from 1985 to 2020 using machine-learning methods with time-series remote sensing data and driving factors. Agric. Water Manag. 2023, 30, 108228. [Google Scholar] [CrossRef]
  14. De Gooijer, J.G.; Hyndman, R.J. 25 years of time series forecasting. Int. J. Forecast. 2006, 22, 443–473. [Google Scholar] [CrossRef]
  15. Chodakowska, E.; Nazarko, J.; Nazarko, Ł.; Rabayah, H.S.; Abendeh, R.M.; Alawneh, R. Arima models in solar radiation forecasting in different geographic locations. Energies 2023, 16, 5029. [Google Scholar] [CrossRef]
  16. Al-Selwi, H.F.; Abas, F.S.; Ahmad, K.A.; Noor, N.M.; Abdul, R.S.F. Attention Based Spatial-Temporal GCN with Kalman filter for Traffic Flow Prediction. Int. J. Technol. 2023, 14, 1299. [Google Scholar] [CrossRef]
  17. Kontopoulou, V.I.; Panagopoulos, A.D.; Kakkos, I.; Matsopoulos, G.K. A review of ARIMA vs. Machine learning approaches for time series forecasting in data driven networks. Future Internet 2023, 15, 255. [Google Scholar] [CrossRef]
  18. Zhang, G.P. Time series forecasting using a hybrid ARIMA and neural network model. Neurocomputing 2003, 50, 159–175. [Google Scholar] [CrossRef]
  19. Kim, J.; Kim, T.; Ryu, J.G.; Kim, J. Spatiotemporal graph neural network for multivariate multi-step ahead time-series forecasting of sea temperature. Eng. Appl. Artif. Intell. 2023, 126, 106854. [Google Scholar] [CrossRef]
  20. Iung, A.M.; Cyrino, O.F.L.; Marcato, A.L. A review on modeling variable renewable energy: Complementarity and spatial–temporal dependence. Energies 2023, 16, 1013. [Google Scholar] [CrossRef]
  21. Daviran, M.; Shamekhi, M.; Ghezelbash, R.; Maghsoudi, A. Landslide susceptibility prediction using artificial neural networks, svms and random forest: Hyperparameters tuning by genetic optimization algorithm. Int. J. Environ. Sci. Technol. 2023, 20, 259–276. [Google Scholar] [CrossRef]
  22. Masini, R.P.; Medeiros, M.C.; Mendes, E.F. Machine learning advances for time series forecasting. J. Econ. Surv. 2023, 37, 76–111. [Google Scholar] [CrossRef]
  23. Ifaei, P.; Nazari-Heris, M.; Charmchi, A.S.; Asadi, S.; Yoo, C. Sustainable energies and machine learning: An organized review of recent applications and challenges. Energy 2023, 1, 126432. [Google Scholar] [CrossRef]
  24. Zheng, Z.; Wang, F.; Gong, G.; Yang, H.; Han, D. Intelligent technologies for construction machinery using data-driven methods. Autom. Constr. 2023, 1, 104711. [Google Scholar] [CrossRef]
  25. Baig, F.; Ali, L.; Faiz, M.A.; Chen, H.; Sherif, M. How accurate are the machine learning models in improving monthly rainfall prediction in hyper arid environment? J. Hydrol. 2024, 1, 131040. [Google Scholar] [CrossRef]
  26. Pirone, D.; Cimorelli, L.; Del, G.G.; Pianese, D. Short-term rainfall forecasting using cumulative precipitation fields from station data: A probabilistic machine learning approach. J. Hydrol. 2023, 1, 128949. [Google Scholar] [CrossRef]
  27. Nocentini, N.; Rosi, A.; Segoni, S.; Fanti, R. Towards landslide space-time forecasting through machine learning: The influence of rainfall parameters and model setting. Front. Earth Sci. 2023, 13, 1152130. [Google Scholar] [CrossRef]
  28. Saha, A.; Singh, K.N.; Ray, M.; Rathod, S.; Dhyani, M. Fuzzy rule–based weighted space–time autoregressive moving average models for temperature forecasting. Theor. Appl. Climatol. 2022, 150, 1321–1335. [Google Scholar] [CrossRef]
  29. Rathod, S.; Saha, A.; Patil, R.; Ondrasek, G.; Gireesh, C.; Anantha, M.S.; Rao, D.V.; Bandumula, N.; Senguttuvel, P.; Swarnaraj, A.K.; et al. Two-stage spatiotemporal time series modelling approach for rice yield prediction & advanced agroecosystem management. Agronomy 2021, 11, 2502. [Google Scholar] [CrossRef]
  30. Basha, S.J.; Prasad, G.L.; Vivek, K.; Kumar, E.S.; Ammannamma, T. Leveraging Ensemble Time-series Forecasting Model to Predict the amount of Rainfall in Andhra Pradesh. In Proceedings of the 2022 2nd International Conference on Artificial Intelligence and Signal Processing (AISP), Vijayawada, India, 12–14 February 2022; pp. 1–7. [Google Scholar]
  31. Garai, S.; Paul, R.K.; Yeasin, M.; Roy, H.S.; Paul, A.K. Machine learning algorithms for predicting rainfall in India. Curr. Sci. 2024, 126, 360. [Google Scholar]
  32. Sekhar, P.H.; Kesavulu, P.K.; Bhupathi, M. Modelling and prediction of coastal Andhra rainfall using ARIMA and ANN models. Int. J. Stat. Appl. Math. 2020, 5, 104–110. [Google Scholar]
  33. Bommala, H.; Kumar, G.B.; Babu, B.S.; Sivajyothi, M.; Kumar, M.V.; Basha, S.J. A Review on Ensemble Model-Based Rainfall Prediction Techniques. In Proceedings of the 2nd International Conference on Intelligent Data Communication Technologies and Internet of Things (IDCIoT), Bengaluru, India, 4–6 January 2024; pp. 1474–1478. [Google Scholar]
  34. Pérez-Alarcón, A.; Garcia-Cortes, D.; Fernández-Alvarez, J.C.; Martínez-González, Y. Improving monthly rainfall forecast in a watershed by combining neural networks and autoregressive models. Environ. Process. 2022, 9, 53. [Google Scholar] [CrossRef]
  35. Kaur, J.; Parmar, K.S.; Singh, S. Autoregressive models in environmental forecasting time series: A theoretical and application review. Environ. Sci. Pollut. Res. 2023, 30, 19617–19641. [Google Scholar] [CrossRef] [PubMed]
  36. Zhang, Y.; Yang, H.; Cui, H.; Chen, Q. Comparison of the ability of ARIMA, WNN and SVM models for drought forecasting in the Sanjiang Plain, China. Nat. Resour. Res. 2020, 29, 1447–1464. [Google Scholar] [CrossRef]
  37. Barrera-Animas, A.Y.; Oyedele, L.O.; Bilal, M.; Akinosho, T.D.; Delgado, J.M.; Akanbi, L.A. Rainfall prediction: A comparative analysis of modern machine learning algorithms for time-series forecasting. Mach. Learn. Appl. 2022, 15, 100204. [Google Scholar] [CrossRef]
  38. Ridwan, W.M.; Sapitang, M.; Aziz, A.; Kushiar, K.F.; Ahmed, A.N.; El-Shafie, A. Rainfall forecasting model using machine learning methods: Case study Terengganu, Malaysia. Ain Shams Eng. J. 2021, 12, 1651–1663. [Google Scholar] [CrossRef]
  39. Praveen, B.; Talukdar, S.; Shahfahad, M.S.; Mondal, J.; Sharma, P.; Islam, A.R.; Rahman, A. Analyzing trend and forecasting of rainfall changes in India using non-parametrical and machine learning approaches. Sci. Rep. 2020, 10, 10342. [Google Scholar] [CrossRef]
  40. Pfeifer, P.E.; Deutrch, S.J. A three-stage iterative procedure for space-time modeling phillip. Technometrics 1980, 22, 35–47. [Google Scholar] [CrossRef]
  41. Pfeifer, P.E.; Deutsch, S.J. Seasonal space-time ARIMA modeling. Geogr. Anal. 1981, 13, 117–133. [Google Scholar] [CrossRef]
  42. Pfeifer, P.E.; Deutrch, S.J. Identification and interpretation of first order space-time ARMA models. Technometrics 1980, 22, 397–408. [Google Scholar] [CrossRef]
  43. Pfeifer, P.E.; Bodily, S.E. A test of space-time arma modelling and forecasting of hotel data. J. Forecast. 1990, 3, 255–272. [Google Scholar] [CrossRef]
  44. Subba, R.T.; Costa, A.A.M. Spatio-temporal modelling of temperature time series: A comparative study. In Time Series Analysis and Applications to Geophysical Systems: Part I; Springer New York: New York, NY, USA, 2004; pp. 123–150. [Google Scholar]
  45. Zhou, M.; Buongiorno, J. Space-time modeling of timber prices. J. Agric. Resour. Econ. 2006, 1, 40–56. [Google Scholar]
  46. Saha, A.; Singh, K.N.; Ray, M.; Rathod, S. A hybrid spatio-temporal modelling: An application to space-time rainfall forecasting. Theor. Appl. Climatol. 2020, 142, 1271–1282. [Google Scholar] [CrossRef]
  47. Kumar, R.R.; Sarkar, K.A.; Dhakre, D.S.; Bhattacharya, D. A Hybrid Space–Time Modelling Approach for Forecasting Monthly Temperature. Environ. Model. Assess. 2023, 28, 317–330. [Google Scholar] [CrossRef]
  48. Rathod, S.; Gurung, B.; Singh, K.N.; Ray, M. An improved space-time autoregressive moving average (STARMA) model for modelling and forecasting of spatio-temporal time-series data. J. Indian Soc. Agric. Stat. 2018, 72, 239–253. [Google Scholar]
  49. Piantari, E.; Rabbani, I.M.; Megasari, R. Seasonal space-time based model for infectious disease prediction. In Proceedings of the AIP Conference Proceedings, Copenhagen, Denmark, 30 September 2023; AIP Publishing: Long Island, NY, USA, 2023. [Google Scholar]
  50. Pai, D.S.; Rajeevan, M.; Sreejith, O.P.; Mukhopadhyay, B.; Satbha, N.S. Development of a new high spatial resolution (0.25 × 0.25) long period (1901–2010) daily gridded rainfall data set over India and its comparison with existing data sets over the region. Mausam 2014, 65, 1–8. [Google Scholar] [CrossRef]
  51. Shankar, S.V.; Chandel, A.; Gupta, R.K.; Sharma, S.; Chand, H.; Kumar, R.; Mishra, N.; Ananthakrishnan, S.; Aravinthkumar, A.; Kumaraperumal, R.; et al. Exploring the dynamics of arrivals and prices volatility in onion (Allium cepa) using advanced time series techniques. Front. Sustain. Food Syst. 2023, 31, 1208898. [Google Scholar]
  52. Alam, N.M.; Sarkar, S.K.; Jana, C.; Raizada, A.; Mandal, D.; Kaushal, R.; Sharma, N.K.; Mishra, P.K.; Sharma, G.C. Forecasting meteorological drought for a typical drought affected area in India using stochastic models. J. Indian Soc. Agric. Stat. 2016, 70, 71–81. [Google Scholar]
  53. Shankar, S.V.; Ajaykumar, R.; Ananthakrishnan, S.; Aravinthkumar, A.; Harishankar, K.; Sakthiselvi, T.; Navinkumar, C. Modeling and forecasting of milk production in the western zone of Tamil Nadu. Asian J. Dairy Food Res. 2023, 42, 427–432. [Google Scholar] [CrossRef]
  54. Alam, N.M.; Mishra, P.K.; Jana, C.; Adhikary, P.P. Stochastic model for drought forecasting for Bundelkhand region in Central India. Indian J. Agric. Sci. 2014, 82, 79–84. [Google Scholar] [CrossRef]
Figure 1. Location map depicting the 119 grids of data points.
Figure 1. Location map depicting the 119 grids of data points.
Water 16 01891 g001
Figure 2. Average (A) observed (B) predicted annual rainfall for the period (2015–2019) in West Bengal.
Figure 2. Average (A) observed (B) predicted annual rainfall for the period (2015–2019) in West Bengal.
Water 16 01891 g002
Table 1. Candidate ARIMA models fitted for the data.
Table 1. Candidate ARIMA models fitted for the data.
LocationsModelLocationsModelLocationsModelLocationsModelLocationsModel
WB1ARIMA(0,1,1)WB26ARIMA(0,1,2)WB51ARIMA(0,1,1)WB76ARIMA(1,2,0)WB101ARIMA(1,1,0)
WB2ARIMA(1,1,0)WB27ARIMA(0,1,2)WB52ARIMA(0,1,1)WB77ARIMA(1,3,2)WB102ARIMA(0,1,1)
WB3ARIMA(1,1,1)WB28ARIMA(0,1,2)WB53ARIMA(0,1,1)WB78ARIMA(1,2,0)WB103ARIMA(2,1,0)
WB4ARIMA(3,2,1)WB29ARIMA(1,1,2)WB54ARIMA(1,1,0)WB79ARIMA(1,0,0)WB104ARIMA(1,1,0)
WB5ARIMA(0,1,1)WB30ARIMA(1,1,0)WB55ARIMA(1,3,2)WB80ARIMA(0,1,1)WB105ARIMA(0,1,1)
WB6ARIMA(0,1,1)WB31ARIMA(2,1,0)WB56ARIMA(2,2,2)WB81ARIMA(2,2,0)WB106ARIMA(1,3,2)
WB7ARIMA(2,1,0)WB32ARIMA(0,1,1)WB57ARIMA(2,2,0)WB82ARIMA(1,2,0)WB107ARIMA(4,2,0)
WB8ARIMA(1,1,1)WB33ARIMA(0,1,1)WB58ARIMA(2,2,2)WB83ARIMA(1,1,0)WB108ARIMA(2,2,0)
WB9ARIMA(1,1,1)WB34ARIMA(0,1,1)WB59ARIMA(2,2,0)WB84ARIMA(0,1,1)WB109ARIMA(0,1,1)
WB10ARIMA(3,2,0)WB35ARIMA(3,2,1)WB60ARIMA(2,2,0)WB85ARIMA(0,1,1)WB110ARIMA(0,1,1)
WB11ARIMA(0,1,2)WB36ARIMA(3,2,2)WB61ARIMA(0,1,2)WB86ARIMA(0,1,1)WB111ARIMA(1,1,0)
WB12ARIMA(1,1,0)WB37ARIMA(0,1,1)WB62ARIMA(0,1,2)WB87ARIMA(2,2,0)WB112ARIMA(2,2,0)
WB13ARIMA(1,1,0)WB38ARIMA(0,1,1)WB63ARIMA(1,1,0)WB88ARIMA(2,2,0)WB113ARIMA(1,3,2)
WB14ARIMA(1,1,1)WB39ARIMA(1,1,0)WB64ARIMA(1,3,2)WB89ARIMA(2,2,2)WB114ARIMA(1,3,2)
WB15ARIMA(2,1,0)WB40ARIMA(1,1,1)WB65ARIMA(0,1,1)WB90ARIMA(3,2,0)WB115ARIMA(1,3,2)
WB16ARIMA(1,1,1)WB41ARIMA(1,1,0)WB66ARIMA(2,2,0)WB91ARIMA(1,2,0)WB116ARIMA(0,1,1)
WB17ARIMA(3,2,0)WB42ARIMA(1,1,0)WB67ARIMA(2,2,0)WB92ARIMA(1,3,2)WB117ARIMA(2,2,0)
WB18ARIMA(3,2,1)WB43ARIMA(1,2,0)WB68ARIMA(2,2,0)WB93ARIMA(1,2,0)WB118ARIMA(0,1,1)
WB19ARIMA(3,2,1)WB44ARIMA(0,1,1)WB69ARIMA(2,2,0)WB94ARIMA(1,2,0)WB119ARIMA(2,3,0)
WB20ARIMA(3,2,1)WB45ARIMA(0,1,1)WB70ARIMA(4,2,0)WB95ARIMA(1,2,0)
WB21ARIMA(3,1,0)WB46ARIMA(2,2,0)WB71ARIMA(2,2,2)WB96ARIMA(1,2,0)
WB22ARIMA(1,1,0)WB47ARIMA(2,2,0)WB72ARIMA(1,1,0)WB97ARIMA(1,2,0)
WB23ARIMA(1,1,0)WB48ARIMA(2,2,0)WB73ARIMA(1,0,0)WB98ARIMA(2,2,0)
WB24ARIMA(1,1,0)WB49ARIMA(2,2,0)WB74ARIMA(3,2,0)WB99ARIMA(1,2,0)
WB25ARIMA(0,1,1)WB50ARIMA(0,1,1)WB75ARIMA(2,2,0)WB100ARIMA(0,1,1)
Table 2. Parameters of STARMA model.
Table 2. Parameters of STARMA model.
Spatial LagSlag 0Slag 1
ModelARMAARMA
φ10θ10φ11θ11
Estimate0.87−0.480.13−0.23
Standard error0.040.060.050.06
t-value19.47−8.332.83−3.71
p-value<0.01<0.01<0.01<0.01
Multivariate Box–Pierce Non-Correlation Test69,668.64(<0.01)
Table 3. Comparison of models based on error metrics.
Table 3. Comparison of models based on error metrics.
ModelARIMASTARMA
LocationsRMSEMAPEMAERMSEMAPEMAE
WB1110.735.2485.210.140.010.12
WB2452.0425.92446.370.130.010.10
WB3156.628.56147.480.200.010.17
WB4245.1811.71212.980.260.010.25
WB5182.0011.42171.000.180.010.15
WB689.083.9769.340.110.000.09
WB7201.989.38169.360.240.010.16
WB8204.1010.68147.500.240.010.17
WB9273.2412.83242.490.300.010.25
WB10603.5831.72524.650.330.020.28
WB11342.3017.26303.930.260.010.24
WB12262.9512.80218.260.290.010.24
WB13233.2013.96216.390.280.020.25
WB14257.9814.20186.960.290.020.26
WB15200.3310.52151.210.210.010.19
WB16121.495.8894.700.220.010.18
WB17592.7127.07484.050.340.020.31
WB18431.6221.05402.870.330.020.31
WB19463.2024.51421.660.290.020.27
WB20392.5817.62342.820.280.280.26
WB21166.699.69143.330.230.010.21
WB22224.0310.06170.910.250.010.20
WB23205.8411.09133.200.260.020.24
WB24329.5421.44254.300.330.020.28
WB25230.5211.81211.360.270.010.25
WB26607.1028.76483.530.660.040.58
WB27419.0021.52358.860.290.020.26
WB28445.1822.80382.420.290.020.27
WB29274.3416.77205.860.190.010.15
WB30260.1314.53186.700.150.010.11
WB31148.109.66142.890.150.010.13
WB32202.4111.09172.930.260.010.23
WB33199.1912.74167.200.260.020.21
WB34127.608.10117.350.180.010.16
WB35446.5925.55417.080.210.010.17
WB36639.3133.70591.390.370.020.29
WB37274.6216.49259.480.300.020.28
WB38281.1516.31260.490.290.020.27
WB39236.8218.96198.950.270.020.21
WB40301.8822.74252.600.220.020.20
WB41425.8022.53351.540.480.020.34
WB42201.4511.72174.050.190.010.14
WB43838.0663.28808.370.210.010.19
WB44326.4720.50282.520.410.030.33
WB45199.5611.23152.790.230.010.16
WB46988.2569.62851.290.300.020.25
WB47998.1468.86860.950.320.020.25
WB48809.5749.31665.160.340.020.30
WB49995.1761.59860.780.310.020.29
WB50247.3014.34198.490.210.010.19
WB51231.1313.49163.220.250.020.23
WB52295.4017.60232.510.370.020.31
WB53685.5123.57486.690.800.030.59
WB54282.1412.45201.750.290.010.20
WB55314.4319.94268.140.270.020.25
WB56379.0519.22307.620.420.030.37
WB57621.4653.50529.940.450.030.32
WB58295.8118.20250.270.410.020.30
WB59861.0155.38726.710.310.020.30
WB60836.4752.33695.290.320.020.28
WB61277.3520.11226.560.240.020.20
WB62286.1115.59235.350.310.020.25
WB63343.6916.12248.140.430.020.33
WB64431.2718.49304.670.480.040.42
WB65374.5926.64334.100.460.030.38
WB661120.06116.75986.040.470.040.37
WB671011.0385.23894.890.350.020.28
WB68499.4035.72455.680.260.010.16
WB69789.8953.21639.640.310.020.25
WB70258.8215.18222.450.250.020.23
WB71334.0315.48245.820.230.020.27
WB72434.9937.31390.530.460.030.36
WB73445.3342.24395.770.480.040.39
WB74890.2593.81828.400.490.040.37
WB75982.2882.60867.870.350.020.26
WB761192.5594.771045.350.320.020.22
WB77374.3823.03338.140.280.020.23
WB78412.4130.12364.460.310.020.23
WB79489.5748.92469.720.490.040.43
WB80330.7122.40240.090.410.030.36
WB811206.94114.071043.410.400.030.33
WB821358.38120.811198.150.350.020.28
WB83273.1318.18209.610.360.020.27
WB84320.8120.10257.390.380.020.27
WB85312.0623.19256.500.390.030.32
WB86326.0722.03244.670.400.030.32
WB871135.30108.46984.010.400.030.32
WB88640.1758.30559.740.390.030.31
WB89226.7917.38197.190.240.020.20
WB90455.3833.49393.610.180.010.16
WB911425.33162.051323.240.330.030.26
WB92211.7114.69155.160.230.020.18
WB931477.81170.731366.800.360.040.30
WB941591.61232.691501.940.450.060.38
WB95746.7784.91665.010.400.040.35
WB96792.6779.47694.530.390.040.33
WB97402.5640.17348.080.450.040.37
WB98395.7318.60310.420.270.010.23
WB991282.4392.151184.960.390.030.35
WB100495.0228.93478.990.310.020.27
WB101132.293.3799.330.180.000.13
WB102196.475.25149.010.300.010.22
WB103409.139.19274.920.480.010.40
WB104436.6710.94385.790.500.010.38
WB105314.999.11264.530.280.010.27
WB1061641.2052.221578.990.470.010.38
WB107516.4413.49449.970.180.000.16
WB1081104.7023.48981.120.660.010.58
WB109412.808.69354.590.390.010.33
WB110463.2910.46413.250.530.010.52
WB111493.468.04324.200.460.010.35
WB1121491.9539.281445.750.560.010.52
WB113487.4612.53443.590.390.010.33
WB114346.728.46323.850.250.010.22
WB1151067.6126.131045.750.330.010.29
WB116363.956.87286.080.570.010.49
WB1172285.0560.292100.020.730.020.64
WB118521.0011.49473.180.660.010.59
WB1191887.6738.061385.060.470.010.40
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

Alam, N.M.; Mitra, S.; Pandey, S.K.; Jana, C.; Ray, M.; Ghosh, S.; Paul Mazumdar, S.; Shankar, S.V.; Saha, R.; Kar, G. Enhanced Spatio-Temporal Modeling for Rainfall Forecasting: A High-Resolution Grid Analysis. Water 2024, 16, 1891. https://doi.org/10.3390/w16131891

AMA Style

Alam NM, Mitra S, Pandey SK, Jana C, Ray M, Ghosh S, Paul Mazumdar S, Shankar SV, Saha R, Kar G. Enhanced Spatio-Temporal Modeling for Rainfall Forecasting: A High-Resolution Grid Analysis. Water. 2024; 16(13):1891. https://doi.org/10.3390/w16131891

Chicago/Turabian Style

Alam, Nurnabi Meherul, Sabyasachi Mitra, Surendra Kumar Pandey, Chayna Jana, Mrinmoy Ray, Sourav Ghosh, Sonali Paul Mazumdar, S. Vishnu Shankar, Ritesh Saha, and Gouranga Kar. 2024. "Enhanced Spatio-Temporal Modeling for Rainfall Forecasting: A High-Resolution Grid Analysis" Water 16, no. 13: 1891. https://doi.org/10.3390/w16131891

APA Style

Alam, N. M., Mitra, S., Pandey, S. K., Jana, C., Ray, M., Ghosh, S., Paul Mazumdar, S., Shankar, S. V., Saha, R., & Kar, G. (2024). Enhanced Spatio-Temporal Modeling for Rainfall Forecasting: A High-Resolution Grid Analysis. Water, 16(13), 1891. https://doi.org/10.3390/w16131891

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