Mathematical Problems in Engineering
Volume 2015 (2015), Article ID 218289, 7 pages
http://dx.doi.org/10.1155/2015/218289

Research Article

Comparison of Stabilization Ability of Models for Hydrological Time Series with a Deterministic Trend

Huantian Xie,¹ Min Xu,² and Dingfang Li²

¹School of Science, Linyi University, Linyi 276005, China
²School of Mathematics and Statistics, Wuhan University, Wuhan 430072, China

Received 1 April 2015; Revised 15 July 2015; Accepted 15 September 2015

Academic Editor: Michael Small

Copyright © 2015 Huantian Xie et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Under influence of climate change and human activities, deterministic trend has been detected and reported in various hydrometeorological observation records. In order to correctly model the stochastic properties, the time series has to be stabilized by removing the trend. Both detrending and differencing have been proposed to fulfill such a task. But the influence of the two stabilizing approaches on the residual series is distinguishing. In this study, ARMA models are constructed based on the above two stabilization approaches for an annual minimum daily discharge series with a deterministic trend. Comparisons are made with respect to stabilization ability, model simulation, and forecasting. Results indicate that the model based on detrending is superior to the one based on differencing in almost all the selected comparison criteria. So detrending is suggested to remove the deterministic trend before using ARMA model to fit the observed data.

1. Introduction

During the last decades, time series analysis and modeling have received considerable attention for simulation, forecasting, and control. In the literature, the linear time series model, particularly the autoregressive moving average (ARMA) model, has been widely applied to different hydrologic and climatic variables such as rainfall [1], runoff [2], droughts [3], and water quality variables [4]. Theoretically, trend of any nonlinear system has its approximation with classical regression models. Numerous studies have been undertaken in different regions of the world with respect to trend estimation by various regression techniques [5–10]. To be consistent with the practical situation in China where most annual discharge series follow Pearson type III (P-III) probability distribution [11, 12], it is wiser to fit the observed data using more flexible polynomials based on P-III distribution than to choose a traditional linear model based on normal distribution blindly.

ARMA model is constructed to fit observed data based on the assumption that the time series is stationary. However, with the development of human society, the climate all over the world is changing gradually, which may result in nonstationarity of the observed data series in the form of deterministic trend. This kind of trend has been detected and reported in various hydrometeorological observation records. In order to correctly model the stochastic properties of the time series, the trend has to be removed to fulfill the stationarity assumption. Traditionally, differencing is utilized to stabilize a nonstationary series, which leads to the autoregressive integrated moving average (ARIMA) model. However, this approach may seriously damage the structure of the residual series [13]. As an alternative, detrending has also been proposed to remove a deterministic trend from a nonstationary time series which does not distort the residual series. Li et al. [14] have applied the approach of detrending to remove the serial correlation of their streamflow data. No result has been documented about the performance of ARMA models based on the two stabilization approaches which is the main objective of our study. As a case study, we take the annual minimum daily discharge series of the Yangtze River at the Hankou hydrological station for the period of 1952–2005 to investigate the stabilization ability.

The paper is organized as follows. First, the ARMA modeling procedure and comparison criteria used in this research are described. Then, ARMA models based on differencing and detrending are constructed to fit the observed discharge series with a deterministic trend. Comparisons are made with different criteria in model calibration and prediction. Finally, suggestions and conclusions are presented at the end of the paper.

2. Methods

2.1. ARMA Model

According to time series theory, difference equations are established based on the observed data series. Among all the equations, ARMA model is the most widely used. The form of an ARMA model is usually written as

where

and is the observed data series, , are model parameters, and are orders of the model. is the backward shift operator (). Random errors are assumed to be white noise. The significance of ARMA model is that the current data can be defined as certain combination of the historical observation data and white noise.

The procedure of building an ARMA model to fit the observed data series often includes model identification, parameter estimation, and model diagnosis. The detailed analysis stages are shown as follows.

(1) Data preparation: transformations of the data (such as square roots or logarithms) can help stabilize the variance in a series where the variation changes with the level. Then, the data are stabilized in the mean level by differencing or detrending. The stabilized data are often easier to model than the original data.

(2) Stationarity and linear correlation checking: stationarity and serial linear correlation checking are necessary to establish ARMA model.

(3) Model structure identification: after data preparation and checking, the order number of the model can be got from the autocorrelation function (ACF) and partial autocorrelation function (PACF) plots.

(4) Parameter estimation: this stage means finding the values of the model coefficients that provide the best fit to the data. There are sophisticated computational algorithms designed to do this. Here, the maximum likelihood method is adopted.

(5) Model checking: this step involves testing the assumptions of the model to identify any areas where the model is inadequate. If the model is found to be inadequate, it is necessary to go back to Step (3) and try to identify a better model.

(6) Model prediction: forecasting is what the whole procedure is designed to accomplish. Once the model has been selected, estimated, and checked, it is usually a straightforward task to compute forecasts. Of course, computer does this.

2.2. Comparison Approach

The ARMA models built in this study are compared through a multicriteria comparison by applying a set of evaluation metrics [15]. The evaluation metrics can be classified into three groups: (1) metrics that calculate the absolute error, (2) metrics that calculate the relative error, and (3) the dimensionless metrics [16].

The used evaluation metrics of the above categories in this study are listed as follows:

(1) Metrics for calculating absolute errors are as follows:

 Absolute maximum error (AME)

 Peak difference (PD)

 Mean absolute error (MAE)

 Root mean squared error (RMSE)

(2) Metrics for calculating relative errors are as follows:

 Relative absolute error (RAE)

 Mean relative error (MRE)

(3) Dimensionless metrics are as follows:

 Coefficient of determination (-squared)

 Coefficient of efficiency (CE)

 Index of agreement (IA)

In the above equations, is the observed time series, is the predicted time series, and and are the mean of the observed and predicted time series, respectively.

3. Application and Comparison

In this section, we illustrate the application of ARMA model for an annual minimum daily discharge series with a statistically significant deterministic trend. ARMA models are constructed to fit the residual time series after the trend is removed from the original series.

In the part of detrending, we adopt three types of polynomials to analyze the trend of annual minimum discharge series of Yangtze River at the Hankou hydrological station: linear trend based on the nonparametric approach of TSA [13], linear regression based on P-III distribution, and quadratic polynomial regression based on P-III distribution [17], which are denoted as Model 1, Model 2, and Model 3, respectively, hereinafter. Besides, ARMA model built on differencing is denoted as Model 4. Estimates of the coefficients are found by maximum likelihood estimate with Newton-Raphson algorithm [18].

Then, comparisons are made with regard to stabilization ability, simulation, and forecasting accuracy between the above four models.

3.1. Model Construction

The annual minimum daily discharge series from Hankou hydrological station in Yangtze River basin of China is chosen as the case study. The discharge series for the period of 1952–2005 is illustrated in Figure 1. The time series plot indicates a deterministic trend present in the mean level of the series. In fact, the value in the linear regression test [19], Mann-Kendall (MK) trend test, and Spearman product-moment correlation test for the study series indicates that an increasing trend is significant at the significance level of 0.05 (Table 1). Meanwhile, the ACF and PACF of the original series (Figure 2) indicate no seasonality in the used series.

Table 1: Trend test result of observed data.

Method
Statistic
value
Trend

Mann-Kendall
3.48
0.01
Increase

Spearman test
3.67
0.01
Increase

Linear regression
34.92
0.01
Increase

Figure 1: Annual minimum flow series of Hankou station.

Figure 2: ACF and PACF for annual minimum flow series for Hankou station.

In order to establish ARMA model, the original data series has to be stabilized by removing the trend. After detrending and differencing, there is no visual trend in the series which is confirmed by the MK trend test. The stabilized series are depicted in Figure 3. From all the above, it can be concluded that both detrending and differencing can remove the trend from a nonstationary time series successfully. However, Yue and Pilon [13] have documented that detrending can remove the trend without distorting the residual series while differencing may seriously damage the existing AR (1) process. To investigate the influence of the two stabilization approaches, ARMA models based on the above stabilized series will be built and compared in multicriteria. The ACF and PACF of the residual time series after detrending (Figure 4) suggest that the stabilized series are stochastic variables, which is verified through model diagnosis. So Model 1, Model 2, and Model 3 are composed of the stochastic process and the deterministic trend component, which leads to the following fitted series [14]:

where is the combined series generated by model of detrending, comes from the residual series model ARMA, and is the deterministic trend estimated by the nonparametric approach of TSA [13], linear and quadratic polynomial regression based on P-III distribution, where form a sequence of Pearson type III distributed random variables.

Figure 3: Stabilized time series by removing trend and differencing.

Figure 4: ACF and PACF for the residual series using detrending.

After first-order difference, the ACF and PACF of the corresponding residual series are presented in Figure 5, and ARMA is selected as the adequate model after model checking. So ARIMA become the time series model for the original data series, which has been denoted as Model 4.

Figure 5: ACF and PACF for the residual series after first-order difference.

3.2. Model Comparison

This section provides the multicriteria for model performance evaluation and comparison of Model 1, Model 2, Model 3, and Model 4 fitted to the discharge data at Hankou station.

In Table 2, the multicriteria introduced in Section 2.2 are considered. The comparison of results listed in it suggests that the detrending-based models perform much better than the differencing-based model according to the evaluation metrics. Meanwhile, Model 3 seems to relatively outperform Model 1 and Model 2 according to evaluation metrics. The evaluation metrics also indicate slightly better performance of Model 2 than Model 1. These findings also confirm that regression model based on P-III distribution is more effective and quadratic polynomial is more flexible to analyze the trend of discharge series in Hankou station.

Table 2: Model criteria for the annual minimum daily discharge series.

Model
AME
PD
MAE
RMSE
RAE
MRE

CE
IA

Model 1
2291
−2291
754
968
0.85
0.02
0.9
0.24
0.90

Model 2
2256
−2256
754
968
0.85
0.02
0.9
0.24
0.90

Model 3
2047
−1781
719
887
0.81
0.02
0.9
0.36
0.91

Model 4
2576
−2576
787
994
0.88
0.01
0.9
0.20
0.89

To further investigate the accuracy of the studied models, the time series plots are depicted in Figure 6. Although the performance of the models looks almost the same for most of the annual minimum daily discharge, it can be seen that, for peak values, Model 4 based on differencing usually reaches the biggest value while Model 3 based on detrending has the largest error for valley data series, which indicates that the accuracy of the two models is poor at extreme points of the series.

Figure 6: Comparison of fitted values for Models 1–4.

The prediction results are presented in Table 3. In the forecast, the first 52 pieces of data are used to determine model parameters and the data of the last 5 years are used for comparison. The forecasting error of Table 3 is the relative error between the prediction value and the real value. The comparisons suggest that model based on detrending is a better model for one- and two-year prediction of the annual minimum daily discharge, while model based on differencing performs better for three- to five-year forecast for the studied series.

Table 3: Predictions for the annual minimum daily discharge series.

Year
Observed
Predicted (relative error: %)

Model 1
Model 2
Model 3
Model 4

2006
7222
7675 (6.3)
7718 (6.9)
8633 (6.9)
8179 (13.2)

2007
7483
7709 (3.0)
7753 (3.6)
8766 (3.6)
8124 (8.6)

2008
8191
7742 (−5.5)
7788 (−4.9)
8903 (−4.9)
8121 (−0.9)

2009
9010
7776 (−13.7)
7823 (−13.2)
9043 (−13.2)
8121 (−0.1)

2010
8681
7809 (−10)
7858 (−9.5)
9187 (−9.5)
8121 (−6.5)

4. Discussion and Conclusion

This study investigates the influence of stabilizing approaches, detrending and differencing, on the performance of the subsequent time series models. ARMA models are constructed based on the above two stabilization approaches for an annual minimum daily discharge series with a deterministic trend. As a study case, the annual minimum daily discharge series of the Yangtze River at the Hankou hydrological station for the period 1954–2005 is analyzed using different stabilizing approaches. Comparisons are made with respect to stabilization ability, model simulation, and forecasting. Results indicate the following:

(1) The model based on detrending is superior to the one based on differencing in almost all the selected comparison criteria. So detrending is suggested to remove the deterministic trend before using ARMA model to represent the observed data.

(2) For the stabilization ability and model simulation, detrending model of the quadratic polynomial regression based on P-III distribution is more effective than the linear regression based on P-III distribution, and the linear regression based on P-III distribution outperforms the linear regression based on TSA method.

(3) As the stabilization power is sensitive to the method utilized to estimate the trend, before we stabilize the series by detrending, it is very necessary to find out the probability distribution type that fits the studied data and flexibly choose a fitted regression curve that represents the real long-term trend best.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research was financially supported by the National Science Foundation of China (51190094), the Natural Science Foundation of Shandong Province (ZR2013FL006, BS2015DX012), AMEP, and Dr. Scientific Research Foundation of Linyi University.

References

A. Astel, J. Mazerski, Ż. Polkowska, and J. Namieśnik, “Application of PCA and time series analysis in studies of precipitation in Tricity (Poland),” Advances in Environmental Research, vol. 8, no. 3-4, pp. 337–349, 2004.

T. B. M. J. Ouarda, J. W. Labadie, and D. G. Fontane, “Indexed sequential hydrologic modeling for hydropower capacity estimation,” Journal of the American Water Resources Association, vol. 33, no. 6, pp. 1337–1349, 1997.

R. Modarres, “Streamflow drought time series forecasting,” Stochastic Environmental Research and Risk Assessment, vol. 21, no. 3, pp. 223–233, 2007.

F. Worrall and T. P. Burt, “A univariate model of river water nitrate time series,” Journal of Hydrology, vol. 214, no. 1–4, pp. 74–90, 1999.

C. Ambrosino, R. E. Chandler, and M. C. Todd, “Southern African monthly rainfall variability: an analysis based on generalized linear models,” Journal of Climate, vol. 24, no. 17, pp. 4600–4617, 2011.

B. C. Bates, R. E. Chandler, and A. W. Bowman, “Trend estimation and change point detection in individual climatic series using exible regression methods,” Journal of Geophysical Research D: Atmospheres, vol. 117, Article ID D16106, 2012.

A. W. Bowman, A. Pope, and B. Ismail, “Detecting discontinuities in nonparametric regression curves and surfaces,” Statistics and Computing, vol. 16, no. 4, pp. 377–390, 2006.

C. M. Hurvich, J. S. Simonoff, and C.-L. Tsai, “Smoothing parameter selection in nonparametric regression using an improved Akaike information criterion,” Journal of the Royal Statistical Society. Series B: Statistical Methodology, vol. 60, no. 2, pp. 271–293, 1998.

N. R. Draper and H. Smith, Applied Regression Analysis, John Wiley & Sons, Hoboken, NJ, USA, 3rd edition, 1998.

J. Krampe, J.-P. Kreiss, and E. Paparoditis, “Hybrid wild bootstrap for nonparametric trend estimation in locally stationary time series,” Statistics and Probability Letters, vol. 101, pp. 54–63, 2015.

D. J. Zhan and S. Z. Ye, Engineering Hydrology, Chinese Water Power Press, Beijing, China, 2000 (Chinese).

L. H. Xiong and S. L. Guo, “Trend test and change-point detection for the annual discharge series of the Yangtze River at the Yichang hydrological station,” Hydrological Sciences Journal, vol. 49, no. 1, pp. 99–114, 2004.

S. Yue and P. Pilon, “Interaction between deterministic trend and autoregressive process,” Water Resources Research, vol. 39, no. 4, p. 1077, 2003.

D. Li, H. Xie, and L. Xiong, “Temporal change analysis based on data characteristics and nonparametric test,” Water Resources Management, vol. 28, no. 1, pp. 227–240, 2014.

R. Modarres, “Multi-criteria validation of artificial neural network rainfall-runoff modeling,” Hydrology and Earth System Sciences, vol. 13, no. 3, pp. 411–421, 2009.

C. W. Dawson, R. J. Abrahart, and L. M. See, “HydroTest: a web-based toolbox of evaluation metrics for the standardised assessment of hydrological forecasts,” Environmental Modelling and Software, vol. 22, no. 7, pp. 1034–1052, 2007.

C. Jiang and L. H. Xiong, “Trend analysis on the annual minimum monthly runoff series of Yangtze river at Hankou hydrological station based on Pearson type III probability distribution,” Journal of Yangtze River Scientific Research Institute, vol. 30, no. 7, pp. 16–21, 2013.

R. A. Rigby and D. M. Stasinopoulos, “Generalized additive models for location, scale and shape,” Applied Statistics, vol. 54, no. 3, pp. 507–554, 2005.

J. M. Utts and R. F. Heckard, Mind on Statistics, China Machine Press, Beijing, China, 2002.