Research Article  Open Access
Prediction of the Reference Evapotranspiration Using a Chaotic Approach
Abstract
Evapotranspiration is one of the most important hydrological variables in the context of water resources management. An attempt was made to understand and predict the dynamics of reference evapotranspiration from a nonlinear dynamical perspective in this study. The reference evapotranspiration data was calculated using the FAO PenmanMonteith equation with the observed daily meteorological data for the period 1966–2005 at four meteorological stations (i.e., Baotou, Zhangbei, Kaifeng, and Shaoguan) representing a wide range of climatic conditions of China. The correlation dimension method was employed to investigate the chaotic behavior of the reference evapotranspiration series. The existence of chaos in the reference evapotranspiration series at the four different locations was proved by the finite and low correlation dimension. A local approximation approach was employed to forecast the daily reference evapotranspiration series. Low root mean square error (RSME) and mean absolute error (MAE) (for all locations lower than 0.31 and 0.24, resp.), high correlation coefficient (CC), and modified coefficient of efficiency (for all locations larger than 0.97 and 0.8, resp.) indicate that the predicted reference evapotranspiration agrees well with the observed one. The encouraging results indicate the suitableness of chaotic approach for understanding and predicting the dynamics of the reference evapotranspiration.
1. Introduction
Evapotranspiration is one of the most important hydrological variables to consider when estimating hydrologic water balance, performing appropriate water resources allocation, establishing efficient irrigation scheduling and assessing hydrological impact of changing climate conditions [1–4]. Evapotranspiration is affected by some factors including weather parameters, crop characteristics, management, and environmental aspects [5]. Specific devices and accurate measurements of various physical parameters or the soil water balance in lysimeters are required to determine evapotranspiration, and thus direct measurement of evapotranspiration is very difficult for long time and large areas because it is timeconsuming and expensive [6]. Reference evapotranspiration () represents the effect of the climate on the evapotranspiration process. provides an upper limit for actual evapotranspiration and indicates the total available energy in humid and arid climate [7]. Some indirect methods (e.g., [5, 8]) have been proposed for estimating evapotranspiration, and they generally have been developed by multiplying reference evapotranspiration by some correction parameters, such as crop coefficient or soil coefficient [9]. Actual evapotranspiration is a key factor of hydrological cycle and agricultural water management. Therefore, in a sense, the accurate forecast for reference evapotranspiration plays an important role in understanding the hydrological cycle, especially in planning and managing irrigation practices.
Reference evapotranspiration is defined as “the rate of evapotranspiration from a hypothetical reference crop with an assumed crop height of 0.12 m, a fixed crop surface resistance of 70 s m^{−1}, and an albedo of 0.23, closely resembling the evapotranspiration from an extensive surface of green grass cover of uniform height, actively growing, and completely shading the ground with adequate water.” Among several existing methods, the FAO56 PenmanMonteith (PM) method [5] is currently widely used and can be considered as a sort of standard [10]. The PM method overcomes shortcomings of the previous Penman method. The PM method provides reference values of potential evapotranspiration for a uniform grass reference surface worldwide and there is no need for local calibration [5]. However, a major drawback to apply the PM method is its relatively high demand for data. Apart from site location, the method requires air temperature, wind speed, relative humidity, and sunshine duration (or shortwave radiation data). Obtaining large meteorological data sets is difficult, especially in developing countries where the number of the meteorological stations is small and reliable data is scarce [11]. This would lead to the discontinuity for reference evapotranspiration series.
Some new methods for predicting and modeling evapotranspiration are developed or researched considering the limitations and disadvantages of the conventionally applied modeling techniques, such as, simple yearly differencing (YD) model and monthly average (MAV) model [12], seasonal autoregressive integrate moving average (SARIMA) model [13], and transfer function noise (TEN) model [14]. Moreover, in recent years, a large number of researches have been devoted to the application of artificial neural network in the area of modeling and predicting evapotranspiration and reference evapotranspiration (e.g., [15–23]). Most of the methods mentioned above, in fact, still need some basic meteorological values measured in the past as input for the forecast. Few studies have been carried out to forecast reference evapotranspiration using the reference evapotranspiration series alone, in large measure because that reference evapotranspiration can be seen as a complex system impacted by numerous climate factors. The full nonlinear dynamic information for these factors is not easily gained from the reference evapotranspiration series, which impacts on the effects of prediction for reference evapotranspiration.
In recent years, the advent and development of chaos theory provide a new idea for complex systems time series forecast. The study of chaos theory considers that, for a complex system with the chaotic features, a longterm time evolution of any single variable time series from the system contains the evolution information of all the variables in the system. Meanwhile, the selfsimilarity of reference evapotranspiration has been revealed by fractal analysis by Liu et al. [24] and Xie et al. [25]. Therefore, it is now believed that the dynamic evolution of the seemingly complex hydrological processes such as reference evapotranspiration can be best understood and forecasted using nonlinear deterministic chaotic models; that is, the chaotic method is capable of describing the dynamic changes occurring in the underlying evapotranspiration system using only a single variable time series.
In the past decade, chaos theory was applied in a variety of hydrological processes, such as rainfall (e.g., [26–32]) runoff and discharge (e.g., [28, 33–40]), sediment transport [23, 41], and lake volume [42, 43]. Outcomes of these studies are encouraging, as the chaos method has been found to be very useful in providing important information regarding dynamical characteristics of the various hydrological phenomena and their predictabilities.
It is surprising to note that though chaos theory has widespread applications in hydrology, to the knowledge of the author, no study has been attempted to applying the chaos analysis method to reference evapotranspiration thus far. Recognizing the above concerns, an attempt is made in the present study to understand and forecast the dynamics of daily reference evapotranspiration with chaotic theory from several locations in China.
This paper is organized as follows: the next section of this paper provides a brief description about how to reconstruct a multidimensional phasespace with a scalar time series, followed by a description of the local approximation forecasting approach. Section 2.3 gives details of the weather stations and data. The results obtained are analyzed and the performance of the chaotic model in forecasting reference evapotranspiration is investigated and discussed in Sections 3 and 4. The summary and concluding remarks are presented in the final section.
2. Methods and Materials
2.1. PhaseSpace Reconstruction Forecasting Approach
Due to joint influence of many factors, a nonlinear system always exhibits intensive complexity. The nonlinear system is so complex that the systematic information from a series may be lost or cannot effectively be extracted, which is usually very useful for systematic evolution. In fact, it is possible to grasp the information about the dynamics of the entire multivariable system by using rational method [41]. The theory describing the evapotranspiration evolution is unknown and hence the dynamics of evapotranspiration cannot be directly determined. A multidimensional phase space from the onedimensional time series representing the dynamics of the evapotranspiration is reconstructed to map the series to a phasespace trajectory. The basic principles of phasespace reconstruction process are described below. Among a variety of methods available for phasespace reconstruction, the method of delays (e.g., [44]) is the most widely used one.
For a scalar time series , where , the phase space can be reconstructed using the method of delays. The basic idea in this method is that the evolution of any single variable of a system is determined by other variables with which it interacts. Information about the relevant variables is thus implicitly contained in the history of any single variable. On the basis of this an “equivalent” phase space can be reconstructed by assigning an element of the time series and its successive delays as coordinates of a new vector time series where , is referred to as the delay time and for a digitized time series is a multiple of the sampling interval used and is the dimension of the vector , called embedding dimension. If the embedding phasespace dimension is large enough, the chaotic attractor, a set of physical properties toward which a system tends to evolve, regardless of the starting conditions of the system, can be depicted.
An appropriate delay time, , is necessary for the phasespace reconstruction because an optimum selection of gives best separation of neighboring trajectories within the minimum embedding phase space (e.g., [45, 46]). A too large or too small would influence the ability of characterizing dynamical system for reconstructed phase space (e.g., [43, 47]).
Many methods have been developed selecting an appropriate delay time. Well known among these are the autocorrelation function method (ACF) (e.g., [48–50]) and the average mutual information method (AMI) (e.g., [51]). The delay time value is selected as the time when the first zerocrossing of the signal autocorrelation function or has dropped to 1−1/ (, approximately 2.7183 the Euler’s constant, which is the base of natural logarithms) of its initial value or where the MI first attains a minimum. No evidence indicates that one of those two methods must be better than the other one [46]; therefore, both of the two methods were adopted in the present study to determinate the delay time of the reference evapotranspiration series.
It is generally believed that there exists a dimension for which the geometric object formed by is equivalent to the original trajectory . A small value of may not adequately reconstruct the original phase space. However, a large value of would reconstruct the original phase but would lead to a high computational cost [52]. Many algorithms have been formulated for the computation of the correlation dimension of a time series, among which GrassbergerProcaccia algorithm [53] is the most popular one and has been successfully applied in a number of studies; therefore, it is employed in the present study too. According to the embedding theorem [44], a dimensional phase space is required to characterize a dynamic system with an attractor dimension . The attractor dimension can be estimated in the following manner. For a dimensional phase space the correlation function is given by where is the Heaviside step function, with for and for , where , is the radius of sphere centered on or , and is the number of data points. If the time series is characterized by an attractor, then for positive values of the correlation function is related to the radius by the following relation: where is a constant and is the slope or the correlation exponent of versus plot given by:
The slope is generally estimated by a least squares fit of a straight line over a certain range of , called the scaling region. For a finite data set, such as the reference evapotranspiration data series, cannot get too close to zero. To handle this situation, a better way of identifying the scaling region is to plot versus and select the apparently linear portion of the graph. The slope of this portion will approximate , to observe whether chaos exists, and the correlation exponent values are plotted against the corresponding embedding dimension values. If the value of the correlation exponent is finite, low and noninteger, then the system is considered to exhibit lowdimensional chaos, and the saturation value of the correlation exponent is defined as the correlation dimension of the attractor. On the contrary, if the correlation exponent increases without bond with increase in the embedding dimension , then the system under investigation is considered as stochastic.
Many prediction methods based on phase space reconstruction have been developed, and the local approximation method successfully applied to hydrological time series forecast [23, 41, 54] is used here and is explained as follows.
The phase space of the time series can be reconstructed to vector (described in this section), and a functional relationship between the vectors can be assumed, such as where and are vectors of dimension , describing the state of the system at times (e.g., current state) and (e.g., future state), respectively. The problem then is to find an appropriate expression for (e.g., the predictor ).
There are many approaches to obtain , according to the way for approximation, which can be broadly divided into global approximation approach and local approximation approach. In global approximation a function valid over the entire state space is considered. On the other hand, the local approximation approach subdivides the domain of the attractor into many subsets, each of which is represented by an approximation valid only in that same subset. Since the local approximation approach can lead to a considerable reduction in the complexity of the representation without degrading the quality of the forecast, it has been considered that, for a very short term, the local approximation can generally provide better results than those obtained using the global approximation approach [41]. Therefore, the local approximation approach is employed in this study. Local polynomial technique, widely used local approximator, is described in detail as follows.
The identification of the sets in which to subdivide the domain is done by fixing a metric and then, given the starting from which the forecast is initiated, identifying neighbors , which constitute the set corresponding to is done. Based on this, the local functions can be constituted. The next neighborhood is . A local polynomial model of this evolution can be expressed by [41] where , , and are constants, which are learned from the training sets. After the nearest neighbors are found, each of these states can be projected to their respective future states and a local predictor can be constructed by a least squares fit minimizing where , , are the nearest neighbors of and are their corresponding onestep ahead states. In the present study, the local maps are learned using local polynomials (e.g., [42]). Based on those local maps, the forecasts are made from a new point and the nearest neighbor in training set is found. Consequently, the onestep evolution of is found that is denoted as and given by And then the nearest neighbor to is found and the procedure is repeated until the desired prediction horizon is reached.
2.2. The FAO PenmanMonteith Method
The PM method for calculating daily reference evapotranspiration can be summarized as follows [5] where is the net radiation at the crop surface (MJ m^{−2} day^{−1}), is the soil heat flux density (MJ m^{−2} day^{−1}), is the mean daily air temperature at 2 m height (°C), is the wind speed at 2 m height (m s^{−1}), is the saturation vapor pressure (kPa), is the actual vapor pressure (kPa), is the saturation vapor pressure deficit (kPa), is the slope of the vapor pressure (kPa°C^{−1}), and is the psychrometric constant (kPa°C^{−1}).
A complete set of equations to compute the parameters of the above equation and all the data required for the calculation of the reference evapotranspiration are given in [5, chapter 3 in the FAO paper 56].
2.3. Study Area and Data
Data used in this study from four National Meteorological observatory (NMO) stations (i.e., Baotou, Zhangbei, Kaifeng, and Shaoguan) is provided by the National climatic Centre (NCC) of CMA (the China Meteorological Administration), including daily observations of maximum, minimum and mean air temperature, wind speed, relative humidity, and sunshine duration for the period of 40 years (from 1966 to 2005). In terms of different latitudes and longitudes, the spatial distribution of stations selected covers a wide range of climatic region in China. The respective locations are shown in Figure 1 and described in Table 1. The FAO56 PM method was the sole method recommended by FAO to calculate reference evapotranspiration wherever the required input data are available. In the present study, therefore, computed using the FAO56 PM method from observed weather variables is taken as the standard.

For each of the four locations, the reference evapotranspiration data computed by FAO56 PM method consisted of 40 years’ (1966–2005) daily data records. Table 2 gives some of the important statistics of the daily reference evapotranspiration series computed by FAO56 PM method from the observed weather variables between 1 January, 1966, and 31 December, 2005. As an example, the variation of reference evapotranspiration series from Kaifeng station is shown in Figure 2. With presenting an evident annual and/or seasonal cycle, Figure 2 clearly indicates that the reference evapotranspiration series obviously present significant variation at the same time. Whether the stochastic or chaotic behaviour cannot be directly identified by the seemingly irregular behavior of the reference evapotranspiration series. Therefore, other methods are required to identify whether or not chaos exists. Similar results are also made for the reference evapotranspiration series of Baotou station, Shaoguan station, and Zhangbei station. The data samples consist of 40 years (1966–2005) of daily reference evapotranspiration data computed by FAO56 PM method. Each data set is divided into two parts: the first part of the series is used to investigate the existence of chaos, to find the appropriate embedding parameters in the phase space and to reconstruct the phase space, where the delay time and the correlation dimension are calculated; the forecast is performed on the other part. Therefore, in the present study, the entire data set of 40 years (January 1966–December 2005) is split into two parts: the first 39 years (January 1966–December 2004) of data are used in training and learning (i.e., the phasespace reconstruction) and predictions are made for the subsequent 1 year (January 2005–December 2005) of data.

3. Results
3.1. Estimation of the Reconstruction Parameters
Figure 3 shows the feature of the twodimensional phase of the reconstructed attractor for the reference evapotranspiration series from Baotou, Shaoguan, Kaifeng, and Zhangbei, respectively, with delay time , that is, the projection of the attractor on the plane .
(a)
(b)
(c)
(d)
How the series of the scalar reference evapotranspiration is reconstructed in a higher dimensional phase space was illustrated in Figure 3. It can be clearly seen from Figure 3 that a welldefined attractor can be reconstructed for all the stations, except a very few outliers corresponding to some very high values of the reference evapotranspiration series. This implied that the feature of the dynamic changes of the seemingly highly irregular reference evapotranspiration can be seen as a relatively simple evolution using the trajectories of the phasespace diagram.
In order to reconstruct the original phase space, the appropriate delay time and embedding dimension should be estimated firstly. ACF and MI for the reference evapotranspiration were calculated, which were shown in Figure 4. Because of strong seasonality, ACF first attains zeros at the lag time of about 1/4 period, namely, 91, 89, 91, and 91 days for Baotou, Zhangbei, Kaifeng, and Shaoguan, respectively. The MI method gives similar estimation for delay time to the ACF method, about approximately 1/4 annual period. The detailed delay time values are 91, 89, 90, and 90 days for Baotou, Zhangbei, Kaifeng, and Shaoguan, respectively. From the delay time results calculated by ACF method and MI method, we can conclude that for the reference evapotranspiration series it really does not matter whether the autocorrelation function or the mutual information is used.
(a)
(b)
(c)
(d)
In the present study, the correlation function and hence the exponents are computed using the delay times determined by the MI method in the previous section, as explained in Section 2.1, for the reference evapotranspiration series. The delay time values of 91, 90, 91, and 89 are obtained for Baotou, Zhangbei, Kaifeng, and Shaoguan, respectively. Correlation dimension method is a common algorithm to determine the optimum embedding dimension and the correlation dimension, meanwhile, correlation dimension method is a method to investigate the chaotic behavior of the time series by identifying whether the correlation exponent value can “become saturation” with the increased embedding dimension. For the reference evapotranspiration time series from the four stations studied, the relationship between the correlation integral, , and the radius, (i.e., versus ), for embedding dimensions, , from 1 to 18, was shown in Figure 5. Figure 5 indicates explicit scaling regions that allow fairly accurate estimates of the correlation exponents. Figure 6 shows the relationship between the correlation exponent values and the embedding dimension values for the four stations. It can be seen that, for the reference evapotranspiration time series from Baotou, Zhangbei, Kaifeng, and Shaoguan, the correlation exponent value increases with the embedding dimension up to 12, 12, 13, and 12, respectively, and then becomes saturation beyond these values. The saturation of the correlation exponent beyond a certain embedding dimension value is an indication of the existence of deterministic dynamics [32, 52]. The finite and low correlation dimensions are observed in all the cases, indicating that the reference evapotranspiration time series from the four stations exhibit lowdimensional chaotic behaviors. The correlation dimension results achieved for the reference evapotranspiration time series from different stations are listed in Table 3.

(a)
(b)
(c)
(d)
A data set with high variability in the values provides a high dimension and indicates the higher complexity of the dynamics of the phenomenon. The dimension results obtained in the present study indicate that the reference evapotranspiration from Shaoguan station has higher variability and irregularity than those from other stations. The nearest integer above the saturation value provides the minimum number of phase spaces or variables necessary to model the dynamics of the attractor. Taking the reference evapotranspiration time series from Shaoguan station as an example, the correlation dimension value about 4.13 indicates that the minimum number of variables necessary to characterize the dynamic feature of the reference evapotranspiration is 5. For other stations in the present study, the minimum number of variables necessary to characterize the dynamic feature of reference evapotranspiration shown in Table 3 all is 4.
3.2. Prediction Results with Chaotic Method
In view of the above analysis, it should be clear that the dynamics of reference evapotranspiration from all the stations under different climates in China can be well described by a chaotic approach. Namely, there are the possibilities of achieving good forecasts for the reference evapotranspiration series, especially with the use of local approximation methods.
The local approximation method explained in Section 2.1 is now employed to forecast the reference evapotranspiration series from the four selected stations in the present study. The delay time and optimum embedding dimensions calculated in the previous section for the reference evapotranspiration series from each station are used for reconstructing the phase space, and predictions are made for one time step ahead (i.e., lead time = 1). The performances of the local approximation prediction method are evaluated using some statistical indices including mean absolute error (MAE), root mean square error (RMSE), correlation coefficient (CC), and modified coefficient of efficiency. The detailed definitions of these error indicators (absolute and relative) are described as follows:
(a) Mean absolute error (MAE) is expressed by
(b) The root mean square error (RMSE) is expressed by where is the number of observations and and are, the th observed (calculated with the FAO56 PM method) and predicted data (using the local approximation prediction method).
(c) Correlation coefficient (CC) is expressed by where and are the mean of the observed and predicted data in the evaluation set, respectively.
(d) Modified coefficient of efficiency is The correlation coefficient ranges from 0 to 1, with higher values indicating better agreement between the observed and the predicted data, whereas the modified coefficient of efficiency ranges from to 1, with higher positive values indicating better agreement.
Figure 7 presents a comparison of the reference evapotranspiration time series for all the stations between the modeled values and the observed values. The scatter plot of the modeled and the observed values is presented in Figure 8, with the solid 1 : 1 (diagonal) line for reference.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
Table 4 presents a summary of forecasted results, in terms of MAE, RMSE, CC, and modified coefficient of efficiency (), for each of the four stations in this study using the chaotic model.

In respect to each station, these results are the best results among the ones achieved with different number of neighbors for the optimal embedding dimension. For each station, the number of neighbors that yields the best results, which can be also called the optimum number of neighbors, is also presented in Table 4. The statistics in Table 4 indicate that reasonably good forecasts are achieved for all the weather stations selected in the present study. Results for CC and values are larger than 0.97 and 0.8 for all the stations, respectively. All MAE and RMSE values are smaller than 0.35 mm day^{−1} and 0.25 mm day^{−1}. A closer look at the statistics shows the best and the relative bad predictions are achieved (in terms of the four evaluation indices) in Baotou station (with , , , ) and Shaoguan station (with , , , ). A comparison of the observed values with the predicted ones in the form of time series plots is presented in Figure 7. As can be seen, the predicted values are, in general, in good agreement with the observed one for all the four stations. A closer look at the fig reveals that the local approximation method can capture very well the major trends and minor fluctuations in the reference evapotranspiration series and even most of extreme values are fairly well predicted except a few very high values (reference evapotranspiration in summer) at Baotou station and Kaifeng station. In addition, the performance mentioned above of the chaotic model can be revealed by the scatter plots of forecasted and observed reference evapotranspiration series shown in Figure 8, where the solid 1 : 1 line is plotted for reference. The regression relationship between the forecasted and observed reference evapotranspiration series at all the stations was shown in Figure 9. Also, the reference evapotranspiration series values forecasted by chaotic model obviously agree well with those computed with fullobserved meteorological data sets. The variance of residuals of the regression is about the same for high, middle, and low values for Baotou and Kaifeng stations, but it is smaller for those two stations (i.e., arid and semihumid location) than the other two stations, Shaoguan and Zhangbei stations (i.e., semiarid and humid location). Moreover, for comparison, autoregressive (AR) model and back propagation (BP) neural network model are also employed in this study to forecast the reference evapotranspiration. The three network layers are adopted in BP neural network model. The input variables include daily mean temperature, sunshine duration, relative humidity, and wind speed in 2 m. The output is daily . The summary of forecasted results by AR model and BP neural network model is also presented in Table 4. It can be apparently seen that the results by chaotic model are similar to that calculated by BP neural network model based on weather data, but it is better than that simulated by AR model. This encouraging result indicates that the chaos theory provides a new method for prediction without weather data.
(a)
(b)
(c)
(d)
4. Discussion
Reference evapotranspiration is a complex variable influenced by a lot of factors, with a significant system dynamic feature, and the suitability of the phasespace reconstruction and local approximation prediction approaches for understanding those dynamics is indicated by the accurate prediction results measured via the statistic indices, the time series plots, the scatter diagrams, and the regression relationship for the forecasted and observed values. The minimum number of variable necessary to describe the dynamics of the reference evapotranspiration can be determined by the correlation dimension values computed by the correlation dimension method. In the present study, the correlation dimension results achieved for the four stations are 3 or 4 (Table 3), which means that the reference evapotranspiration dynamics is dominantly influenced by 4 or 5 variables. The results are similar to the research for other hydrological series with the chaos theory; for example, only three variables dominantly influenced the suspended sediment concentration dynamics, even though a large number of variables may have influence on its dynamics [41]. What are these 4 or 5 variables? The wellknown facts are that the several important impact factors of the reference evapotranspiration are temperature, air humidity, wind speed, and sunshine hours (or solar radiation). With regard to the present results, a possible implication is that temperature, air humidity, wind speed, and sunshine hours are exactly the 4 or 5 variables achieved by the correlation dimension approach. Nevertheless, a further confirmation and verification still should be made.
In any case, a multivariable time series phasespace reconstruction for a single variable time series can recover the dynamics of reference evapotranspiration, which have been echoed by the results that reference evapotranspiration series can be well forecasted by chaotic model.
Different from the reported reference evapotranspiration forecast methods, which are established based on the statistics of the time series variation, such as yearly differencing (YD) or monthly average model (MAV) [12], the seasonal autoregressive integrated moving average model (SARIMA) [16], transfer function noise (TFN) [14], and ANN method [16], chaotic forecast model is based on the representation of the reference evapotranspiration dynamics in the reconstructed phasespace step by step in local neighborhoods.
5. Summary
The analytical and forecasted method, based on the chaos theory, used to understand and predict the dynamics of reference evapotranspiration series, was tested for four weather stations (namely, Baotou, Zhangbei, Kaifeng, and Shaoguan) located at different latitudes and longitudes with different climates. The studied reference evapotranspiration time series of the four weather stations is over a period of 40 years (January 1966–December 2005), which is calculated by FAO56 PenmanMonteith model with the daily weather data. The delay time of the reference evapotranspiration series was first estimated with both the zerocrossing of the autocorrelation function and the first minimum of the mutual information. For representing the dynamics of the reference evapotranspiration, a multidimensional phase space was reconstructed for the scalar reference evapotranspiration series. The correlation dimension method was used to investigate the presence of chaos in the reference evapotranspiration series at the four locations selected in this study, and then a local approximation method was employed for making predictions of the reference evapotranspiration.
The results have shown that chaotic characteristics obviously exist in the reference evapotranspiration series at the four different locations (i.e., Baotou, Zhangbei, Kaifeng, and Shaoguan), duo to the finite and low correlation dimension (i.e., 3.41, 3.76, 3.83, and 4.13, resp.). The good performance for the chaotic approach to forecast the reference evapotranspiration series was presented by the accurate predictions for the four locations and measured via statistics indices (i.e., low RSME and MAE; high CC and modified coefficient of efficiency values for all the four stations), suggesting the effectiveness of the chaotic approach for reference evapotranspiration forecast. These results indicated that, seemingly, the variation of the reference evapotranspiration is irregular, in fact, which can be seen as the result of a simple deterministic system influenced by several dominant variables, the number of which determined by the correlation dimension value of the series. The correlation dimension value are 4, 4, 4, and 5, respectively, for Baotou, Zhangbei, Kaifeng, and Shaoguan. In regard to this, those dominant variables may be the temperature, relative humidity, wind speed, and solar radiation, which are considered as the important factors based on physical mechanism.
To the knowledge of the author, this study was the first ever study conducted on the analysis of reference evapotranspiration series based on the chaos theory. However, the present investigation is only the first step to lead to more complete understanding of the dynamics of the hydrological processes for the reference evapotranspiration. Our attempt should continue by applying the chaos concept to time series of the physical factors of the reference evapotranspiration, such as temperature, relative humidity, solar radiation, elevation, and wind speed, which may allow us to provide general and strong conclusion regarding the use of such a concept, for example, the physical explanation to the chaos phenomenon. Moreover, comparing the chaotic approach with other existing forecast approaches to determine a better approach for the evapotranspiration system is the future work. In a word, the authors believe that the present study is a beneficial attempt in these directions, and the further studies in this field are under way and will be addressed in forthcoming papers.
Conflict of Interests
The authors of the paper do not have a direct financial relation with the commercial identity mentioned in this papert that might lead to a conflict of interests.
Acknowledgments
This study was sponsored by the Development Program of China during the 12th Fiveyear Plan Period (2013BAC10B01), the National Natural Science Foundation of China (no. 51379057, no. 51009046, and no. 51109056), Qing Lan Project, the Fund of Advanced Science and Technology Innovation in Colleges and Universities in Jiangsu Province, the Program for Distinguished Talents in Hohai University, and a Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD). The authors thank the National Meteorological Information Centre of China for providing the meteorological data. Finally, cordial thanks should be extended to the Editor and anonymous referees for their valuable comments which greatly improved the quality of this paper.
References
 L. Tong, S. Kang, and L. Zhang, “Temporal and spatial variations of evapotranspiration for spring wheat in the Shiyang river basin in northwest China,” Agricultural Water Management, vol. 87, no. 3, pp. 241–250, 2007. View at: Publisher Site  Google Scholar
 D. Chen, M. Webber, J. Chen, and Z. Luo, “Emergy evaluation perspectives of an irrigation improvement project proposal in China,” Ecological Economics, vol. 70, no. 11, pp. 2154–2162, 2011. View at: Publisher Site  Google Scholar
 W. Wang, S. Peng, T. Yang, Q. Shao, J. Xu, and W. Xing, “Spatial and temporal characteristics of reference evapotranspiration trends in the Haihe River basin, China,” Journal of Hydrologic Engineering, vol. 16, no. 3, pp. 239–252, 2011. View at: Publisher Site  Google Scholar
 W. Wang, Q. Shao, S. Peng et al., “Reference evapotranspiration change and the causes across the Yellow River Basin during 1957–2008 and their spatial and seasonal differences,” Water Resources Research, vol. 48, no. 5, Article ID W05530, 2012. View at: Publisher Site  Google Scholar
 R. G. Allen, L. S. Pereira, D. Raes, and M. Smith, “Crop evapotranspiration: guidelines for computing crop water requirements,” FAO Irrigation and Drainge Paper 56, Rome, Italy, 1998. View at: Google Scholar
 Y. Dinpashoh, “Study of reference crop evapotranspiration in I.R. of Iran,” Agricultural Water Management, vol. 84, no. 12, pp. 123–129, 2006. View at: Publisher Site  Google Scholar
 L. Gong, C.Y. Xu, D. Chen, S. Halldin, and Y. D. Chen, “Sensitivity of the PenmanMonteith reference evapotranspiration to key climatic variables in the Changjiang (Yangtze River) basin,” Journal of Hydrology, vol. 329, no. 34, pp. 620–629, 2006. View at: Publisher Site  Google Scholar
 M. Beyazgül, Y. Kayam, and F. Engelsman, “Estimation methods for crop water requirements in the Gediz Basin of western Turkey,” Journal of Hydrology, vol. 229, no. 12, pp. 19–26, 2000. View at: Publisher Site  Google Scholar
 G. W. Kite and P. Droogers, “Comparing evapotranspiration estimates from satellites, hydrological models and field data,” Journal of Hydrology, vol. 229, no. 12, pp. 3–18, 2000. View at: Publisher Site  Google Scholar
 C.Y. Xu, L. Gong, T. Jiang, D. Chen, and V. P. Singh, “Analysis of spatial distribution and temporal trend of reference evapotranspiration and pan evaporation in Changjiang (Yangtze River) catchment,” Journal of Hydrology, vol. 327, no. 12, pp. 81–93, 2006. View at: Publisher Site  Google Scholar
 P. Droogers and R. G. Allen, “Estimating reference evapotranspiration under inaccurate data conditions,” Irrigation and Drainage Systems, vol. 16, no. 1, pp. 33–45, 2002. View at: Publisher Site  Google Scholar
 J. C. Tracy, M. A. Marino, and S. A. Taghavi, “Predicting water demand in agricultural regions using time series forecasts of reference crop evapotranspiration,” in Water Resources Planning and Management: Saving a Threatened Resourcein Search of Solutions, M. Karamouz, Ed., pp. 50–55, American Society of Civil Engineer, New York, NY, USA, 1992. View at: Google Scholar
 M. A. Marino, J. C. Tracy, and S. A. Taghavi, “Forecasting of reference crop evapotranspiration,” Agricultural Water Management, vol. 24, no. 3, pp. 163–187, 1993. View at: Publisher Site  Google Scholar
 T. Hameed, M. A. Marino, and R. H. Shumway, “Evapotranspiration transferfunctionnoise modeling,” Journal of Irrigation and Drainage Engineering, vol. 121, no. 2, pp. 159–169, 1995. View at: Publisher Site  Google Scholar
 K. P. Sudheer, A. K. Gosain, and K. S. Ramasastri, “Estimating actual evapotranspiration from limited climatic data using neural computing technique,” Journal of Irrigation and Drainage Engineering, vol. 129, no. 3, pp. 214–218, 2003. View at: Publisher Site  Google Scholar
 S. Trajkovic, B. Todorovic, and M. Stankovic, “Forecasting of reference evapotranspiration by artificial neural networks,” Journal of Irrigation and Drainage Engineering, vol. 129, no. 6, pp. 454–457, 2003. View at: Publisher Site  Google Scholar
 Ö. Kisi, G. Yildirim, K. P. Sudheer, A. K. Gosain, and K. S. Ramasastri, “Discussion of “Estimating actual evapotranspiration from limited climatic data using neural computing technique” by K. P. Sudheer, A. K. Gosain, and K. S. Ramasastri,” Journal of Irrigation and Drainage Engineering, vol. 131, no. 2, pp. 219–220, 2005. View at: Publisher Site  Google Scholar
 O. Kisi and G. Yildirim, “Discussion of “forecasting of reference evapotranspiration by artificial neural networks” by S. Trajkovic, B. Todorovic, and M. Stankovic,” Journal of Irrigation and Drainage Engineering, vol. 131, pp. 390–391, 2005. View at: Google Scholar
 S. Trajkovic, “Temperaturebased approaches for estimating reference evapotranspiration,” Journal of Irrigation and Drainage Engineering, vol. 131, no. 4, pp. 316–323, 2005. View at: Publisher Site  Google Scholar
 O. Kisi, “Evapotranspiration modelling from climatic data using a neural computing technique,” Hydrological Processes, vol. 21, no. 14, pp. 1925–1934, 2007. View at: Publisher Site  Google Scholar
 J. B. Cai, Y. Liu, T. W. Lei, and L. S. Pereira, “Estimating reference evapotranspiration with the FAO PenmanMonteith equation using daily weather forecast messages,” Agricultural and Forest Meteorology, vol. 145, no. 12, pp. 22–35, 2007. View at: Publisher Site  Google Scholar
 G. Landeras, A. OrtizBarredo, and J. J. López, “Comparison of artificial neural network models and empirical and semiempirical equations for daily reference evapotranspiration estimation in the Basque Country (Northern Spain),” Agricultural Water Management, vol. 95, no. 5, pp. 553–565, 2008. View at: Publisher Site  Google Scholar
 B. Sivakumar and A. W. Jayawardena, “An investigation of the presence of lowdimensional chaotic behaviour in the sediment transport phenomenon,” Hydrological Sciences Journal, vol. 47, no. 3, pp. 405–416, 2002. View at: Publisher Site  Google Scholar
 B. Liu, D. Shao, and X. Shen, “Research on temporal fractal features of reference evapotranspiration,” Journal of Irrigation Drainage, vol. 25, no. 5, pp. 9–13, 2006 (Chinese). View at: Google Scholar
 X. Xie, Y. Cui, and Y. Zhou, “Longterm correlation and multifractality of reference crop evapotranspiration time series,” Journal of Hydraulic Engineering, vol. 39, no. 12, pp. 1327–1333, 2008 (Chinese). View at: Google Scholar
 I. RodriguezIturbe, F. B. De Power, M. B. Sharifi, and K. P. Georgakakos, “Chaos in rainfall,” Water Resources Research, vol. 25, no. 7, pp. 1667–1675, 1989. View at: Publisher Site  Google Scholar
 R. Berndtsson, K. Jinno, A. Kawamura, J. Olsson, and S. Xu, “Dynamical systems theory applied to longterm temperature and precipitation time series,” Trends Hydrology, vol. 1, pp. 291–297, 1994. View at: Google Scholar
 A. W. Jayawardena and F. Lai, “Analysis and prediction of chaos in rainfall and stream flow time series,” Journal of Hydrology, vol. 153, no. 1–4, pp. 23–52, 1994. View at: Publisher Site  Google Scholar
 D. Koutsoyiannis and D. Pachakis, “Deterministic chaos versus stochasticity in analysis and modeling of point rainfall series,” Journal of Geophysical Research D, vol. 101, no. 21, pp. 26441–26451, 1996. View at: Publisher Site  Google Scholar
 C. E. Puente and N. Obregón, “A deterministic geometric representation of temporal rainfall: results for a storm in Boston,” Water Resources Research, vol. 32, no. 9, pp. 2825–2839, 1996. View at: Publisher Site  Google Scholar
 B. Sivakumar, S. Liong, C. Liaw, and K. Phoon, “Singapore rainfall behavior: chaotic?” Journal of Hydrologic Engineering, vol. 4, no. 1, pp. 38–48, 1999. View at: Publisher Site  Google Scholar
 B. Sivakumar, R. Berndttson, J. Olsson, and K. Jinno, “Evidence of chaos in the rainfallrunoff process,” Hydrological Science Journal, vol. 5, pp. 645–651, 2001. View at: Google Scholar
 A. Porporato and L. Ridolfi, “Nonlinear analysis of river flow time sequences,” Water Resources Research, vol. 33, no. 6, pp. 1353–1367, 1997. View at: Publisher Site  Google Scholar
 Q. Liu, S. Islam, I. RodriguezIturbe, and Y. Le, “Phasespace analysis of daily streamflow: characterization and prediction,” Advances in Water Resources, vol. 21, no. 6, pp. 463–475, 1998. View at: Publisher Site  Google Scholar
 I. Krasovskaia, L. Gottschalk, and Z. W. Kundzewicz, “Dimensionality of Scandinavian river flow regimes,” Hydrological Sciences Journal, vol. 44, no. 5, pp. 705–723, 1999. View at: Publisher Site  Google Scholar
 N. Lambrakis, A. S. Andreou, P. Polydoropoulos, E. Georgopoulos, and T. Bountis, “Nonlinear analysis and forecasting of a brackish karstic spring,” Water Resources Research, vol. 36, no. 4, pp. 875–884, 2000. View at: Publisher Site  Google Scholar
 A. Elshorbagy, U. S. Panu, and S. P. Simonovic, “Analysis of crosscorrelated chaotic streamflows,” Hydrological Sciences Journal, vol. 46, no. 5, pp. 781–794, 2001. View at: Publisher Site  Google Scholar
 A. Elshorbagy, S. P. Simonovic, and U. S. Panu, “Estimation of missing streamflowdata using principles of chaos theory,” Journal of Hydrology, vol. 255, no. 1–4, pp. 123–133, 2002. View at: Publisher Site  Google Scholar
 Y. Zhou, Z. Ma, and L. Wang, “Chaotic dynamics of the flood series in the Huaihe River Basin for the last 500 years,” Journal of Hydrology, vol. 258, no. 1–4, pp. 100–110, 2002. View at: Publisher Site  Google Scholar
 W. W. Ng, U. S. Panu, and W. C. Lennox, “Chaos based Analytical techniques for daily extreme hydrological observations,” Journal of Hydrology, vol. 342, no. 12, pp. 17–41, 2007. View at: Publisher Site  Google Scholar
 B. Sivakumar, “A phasespace reconstruction approach to prediction of suspended sediment concentration in rivers,” Journal of Hydrology, vol. 258, no. 1–4, pp. 149–162, 2002. View at: Publisher Site  Google Scholar
 H. D. I. Abarbanel and U. Lall, “Nonlinear dynamics of the Great Salt Lake: system identification and prediction,” Climate Dynamics, vol. 12, no. 4, pp. 287–297, 1996. View at: Google Scholar
 T. B. Sangoyomi, U. Lall, and H. D. I. Abarbanel, “Nonlinear dynamics of the Great Salt Lake: dimension estimation,” Water Resources Research, vol. 32, no. 1, pp. 149–159, 1996. View at: Publisher Site  Google Scholar
 F. Takens, “Detecting strange attractors in turbulence,” in Dynamical systems and Turbulence, D. A. Rand and L. S. Young, Eds., vol. 898 of Lecture Note in Mathematics, pp. 366–381, Springer, Berlin, Germany, 1981. View at: Publisher Site  Google Scholar  MathSciNet
 T. Frison, “Nonlinear data analysis techniques,” in Tranding on the Edge: Neural, Genetic, and Fuzzy Systems for Chaotic Financial Markets, G. J. Deboeck, Ed., pp. 280–296, Wiley, New York, NY, USA, 1994. View at: Google Scholar
 B. Sivakumar, “Chaos theory in hydrology: important issues and interpretations,” Journal of Hydrology, vol. 227, no. 1–4, pp. 1–20, 2000. View at: Publisher Site  Google Scholar
 J. W. Havstad and C. L. Ehlers, “Attractor dimension of nonstationary dynamical systems from small data sets,” Physical Review A, vol. 39, no. 2, pp. 845–853, 1989. View at: Publisher Site  Google Scholar  MathSciNet
 J. Holzfuss and G. MayerKress, “An approach to errorestimation in the application of dimension algorithms,” in Dimensions and Entropies in Chaotic System, G. MayerKress, Ed., pp. 114–122, Springer, New York, NY, USA, 1986. View at: Publisher Site  Google Scholar  MathSciNet
 H. G. Schuster, Deterministic Chaos, VCH, Weinheim, Germany, 1998.
 A. A. Tsonis and J. B. Elsner, “The weather attractor over very short timescales,” Nature, vol. 333, no. 6173, pp. 545–547, 1988. View at: Publisher Site  Google Scholar
 A. M. Fraser and H. L. Swinney, “Independent coordinates for strange attractors from mutual information,” Physical Review A, vol. 33, no. 2, pp. 1134–1140, 1986. View at: Publisher Site  Google Scholar  MathSciNet
 P. Shang, X. Na, and S. Kamae, “Chaotic analysis of time series in the sediment transport phenomenon,” Chaos, Solitons and Fractals, vol. 41, no. 1, pp. 368–379, 2009. View at: Publisher Site  Google Scholar
 P. Grassberger and I. Procaccia, “Measuring the strangeness of strange attractors,” Physica D: Nonlinear Phenomena, vol. 9, no. 12, pp. 189–208, 1983. View at: Publisher Site  Google Scholar  MathSciNet
 B. Sivakumar, S. Sorooshian, H. Vijai Gupta, and X. Gao, “A chaotic approach to rainfall disaggregation,” Water Resources Research, vol. 37, no. 1, pp. 61–72, 2001. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Weiguang Wang 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.