Next Article in Journal
Estimation of Spatial and Temporal Groundwater Balance Components in Khadir Canal Sub-Division, Chaj Doab, Pakistan
Next Article in Special Issue
RASPOTION—A New Global PET Dataset by Means of Remote Monthly Temperature Data and Parametric Modelling
Previous Article in Journal
Riparian Vegetation Density Mapping of an Extremely Densely Vegetated Confined Floodplain
Previous Article in Special Issue
Evapotranspiration Trends and Interactions in Light of the Anthropogenic Footprint and the Climate Crisis: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stochastic Analysis of Hourly to Monthly Potential Evapotranspiration with a Focus on the Long-Range Dependence and Application with Reanalysis and Ground-Station Data

by
Panayiotis Dimitriadis
,
Aristoteles Tegos
* and
Demetris Koutsoyiannis
Department of Water Resources and Environmental Engineering, School of Civil Engineering, National Technical University of Athens, Heroon Polytechneiou 5, 15880 Zographou, Greece
*
Author to whom correspondence should be addressed.
Hydrology 2021, 8(4), 177; https://doi.org/10.3390/hydrology8040177
Submission received: 1 November 2021 / Revised: 24 November 2021 / Accepted: 29 November 2021 / Published: 1 December 2021
(This article belongs to the Special Issue Advances in Evaporation and Evaporative Demand)

Abstract

:
The stochastic structures of potential evaporation and evapotranspiration (PEV and PET or ETo) are analyzed using the ERA5 hourly reanalysis data and the Penman–Monteith model applied to the well-known CIMIS network. The latter includes high-quality ground meteorological samples with long lengths and simultaneous measurements of monthly incoming shortwave radiation, temperature, relative humidity, and wind speed. It is found that both the PEV and PET processes exhibit a moderate long-range dependence structure with a Hurst parameter of 0.64 and 0.69, respectively. Additionally, it is noted that their marginal structures are found to be light-tailed when estimated through the Pareto–Burr–Feller distribution function. Both results are consistent with the global-scale hydrological-cycle path, determined by all the above variables and rainfall, in terms of the marginal and dependence structures. Finally, it is discussed how the existence of, even moderate, long-range dependence can increase the variability and uncertainty of both processes and, thus, limit their predictability.

1. Introduction

Evapotranspiration is a paramount element in hydrology, with relevance in many aspects of the geosciences. From hydrological and agronomic perspectives, the potential evapotranspiration (PET) and (potential) evaporation (PEV) are key for water balance estimation, the assessment of crop water demand, and integrated rainfall-runoff modelling. PET [1] is defined as “the amount of water transpired in a given time by a short green crop, completely shading the ground, of uniform height and with adequate water status in the soil profile”. A particular (reference) case thereof is the reference evapotranspiration (ETo), which refers to “the rate of evapotranspiration from a hypothetical reference crop with an assumed crop height of 0.12 m, a fixed surface resistance of 70 s/m, and an albedo of 0.23, closely resembling the evapotranspiration from an extensive surface of green (cool season) grass of uniform height, actively growing, well-watered, and completely shading the ground” [2]. Evaporation is the physical process by which liquid water enters the atmosphere as water vapor. In what follows, when we refer to all of the above processes, we use the acronym PE. We also note that PE may be different from the actual evapo(transpi)ration (in cases where there is not adequate water availability).
For the PE assessment, historically, many models have been developed highlighting the Penman–Monteith model as the most suitable [3]. One of the main shortcomings of estimating PE with the Penman–Monteith model is the requirement of a significant number of meteorological inputs such as, without distinction, temperature, radiation, relative humidity, and wind speed [4,5,6]. In the case that we require a synthetic PE timeseries for risk management (e.g., in a Monte Carlo simulation framework), when the above meteorological inputs are not available for the requested period, then one may use a stochastic model that preserves the important statistical attributes of PE. Additionally, due to the physical complexity of assessing the PE, stochastic modelling provides a solid scientific ground for further consideration in several fields of PE assessment, and the stochastic analysis can contribute to the PE physical interpretation along with other hydrometeorological processes because stochastics is proven as a collection of mathematical tools able to give physical explanations [7]. As highlighted in the aforementioned work [7], the role of stochastics is crucial: (a) to infer dynamics (laws) from past data; (b) to formulate the complex natural system equations; (c) to estimate the involved parameters; and (d) to test any hypothesis regarding the dynamics. There are only limited works providing a thorough stochastic analysis in PE timeseries, even though the necessity of stochastic modelling is of paramount importance. Based on the published literature, a seasonal ARIMA model and Winters’ exponential smoothing model [8] have been investigated for their applicability for forecasting weekly reference crop ETo [9]. Both models demonstrated satisfactory results compared to a simple PE model. Pandey et al. [10] provided a stochastic analysis in assessing black gram evapotranspiration regimes using a long-term pan-evaporation dataset of 23 years in Udaipur, India. Black gram is an important crop of the Udaipur region, and the lack of long-term crop demand assessment led to the need for stochastic analysis using pan-evaporation gauges to predict daily black gram evapotranspiration. As noted by the authors, the new stochastic model for black gram evapotranspiration was found to predict daily black gram evapotranspiration with high accuracy (R2 = 0.94). Dynamic stochastic modelling, with a focus on the marginal probability distribution function (known as cumulative distribution function), has been also used for quantifying the PE uncertainty associated with irrigation scheduling [11,12,13]. Recently, an application of vine copulas with a focus on the short-term structure of the daily evaporation process has been presented [14]. Rainfall-runoff approaches have been presented using stochastic inputs of precipitation and PE to overcome the lack of Penman–Monteith estimates and long-term gauge inputs [15,16].
A substantial amount of previous works have focused on the trend PE assessment [17,18] in conjunction with the well-known term, evaporation paradox [19,20]. The later has been defined as the assumption that, under warming climate and higher temperatures, increased PE rates are expected; however, gauge data show the opposite because observations across the U.S. and the globe show a decreasing trend in pan evaporation. Recent studies recommend the revision of common trend tests through re-evaluation of the statistical significance of an observed trend in a timeseries by assuming a model exhibiting the scaling hypothesis [21], which is shown to be apparent in most key hydrological-cycle processes [22] and provides a more accurate modelling framework than a trend-based approach [23].
The stochastic structure of the PE process, ranging from hourly to climatic scales, is studied here in terms of Hurst–Kolmogorov (HK) dynamics, which describes all processes exhibiting the Hurst phenomenon (i.e., with a power-law autocorrelation function at large scales). Additionally, we focus on the marginal structure of the PE process as fitted through the Pareto–Burr–Feller (PBF) distribution function [24], which includes a large variety of tail-behaviors [25]. Both marginal and second-order dependence structures of the HK dynamics are estimated and compared to the ones identified from global-scale analyses in other key hydrometeorological processes that form the hydrological-cycle path driven by atmospheric turbulence [26], such as temperature, wind, solar radiation, and relative humidity [22,27,28,29].
Because observations for the PE process are usually found on monthly or daily resolutions, here we use two datasets. The first dataset comprises PET timeseries with monthly resolution extracted from the California Irrigation Management Information System (CIMIS) network in California, comprising 41 ground stations. For the second dataset, we extracted gridded reanalysis PEV data of hourly resolution. In particular, we retrieved the reanalysis data for the grid points in the same area of the network of the ground stations, so that we could compare its stochastic structures to the PET records, with a focus on the long-range dependence (LRD) behavior.
In Section 2, we introduce the methodology on the estimation of the marginal and second-order dependence structures, while in Section 3, we present the statistical characteristics of the selected stations as well as the results obtained from the analysis, with a focus on the marginal and the dependence structures. Finally, in Section 4 and Section 5, we summarize our findings, and we discuss how the results may be consistent with the ones obtained from the hydrological-cycle path under HK dynamics, expanding from Gaussian to Pareto-type tail behavior, and from fractal and intermittent behavior at small scales to LRD behavior at large scales.

2. Metrics of Marginal and Dependence Structures

The estimators and models applied for both the marginal and the second-order dependence structures are part of the stochastic framework of the HK dynamics, with a focus on the LRD behavior [30,31,32,33,34,35], and they have been applied to turbulent and key hydrological-cycle processes of global networks with resolutions spanning from small scales (relevant to the fractal behavior) to climatic scales (for a review, see [26]).
It has been shown that a flexible probability distribution function, which seems to fit well a great variety of key hydrological-cycle processes [25,26], with tail-behaviors ranging from Gaussian to Pareto, is the PBF distribution function [24,36,37,38], i.e.:
F ( x ) = P { x _ x } = 1 ( 1 + ζ ξ ( x d λ ) ζ ) 1 ξ ζ
where x > d , d is a location parameter (in units of x), ζ and ξ are dimensionless shape parameters, and λ is a scale parameter (in units of x). It is noted that here the Dutch convention is adopted, where underlined symbols denote random variables and stochastic processes.
The estimation of the parameters of the PBF distribution function for the identification of the marginal structure of the PE process is based on the first four statistical moments, and particularly on the central moments and coefficients (i.e., mean, variance, skewness, and kurtosis). It is stressed that, although the estimation from the classical moments of high order are unknowable, especially in the presence of heavy tails and LRD [25], the hourly PEV and the monthly PET processes are expected to be close to a light-tail behavior and, therefore, the estimation of skewness and kurtosis coefficients could be, in approximation, reliably estimated from data.
For the dependence structure of the PE processes, we select the climacogram metric, which is defined as the variance of the averaged process at the scale domain [7]. i.e.:
γ ( k ) : = Var [ 0 k x _ ( y ) d y ] / k 2
where k is the scale (in units of x _ ). (See discussion on the origins of the name, mathematical definitions, etc., in [26,39])
It has been shown that the climacogram estimator at the scale domain is a more powerful estimator than the autocovariance function at the lag domain or the power-spectrum at the frequency domain [34], while its classical estimator adjusted for bias is defined as [40]:
γ ^ _ ( κ Δ ) = 1 n / κ i = 1 n / κ ( x _ i ( κ ) μ _ ^ ) 2 + γ ( n / κ κ Δ )
where κ = k / Δ is the dimensionless scale, Δ is the time resolution of the process, μ _ ^ is the mean of the process, n / κ is the integer part of n/κ, and x _ i ( κ ) is the i-th element of the averaged sample of the process at scale κ, i.e.:
x _ i ( κ ) = 1 κ j = ( i 1 ) κ + 1 i κ x _ j
For the climacogram model, contained in the above estimator, we select a generalization of the HK model (for details and more sophisticated models, see [25,26]), which has been shown to well simulate processes from sub-hourly to over-annual resolutions, and from short- to long-term scales associated with fractal and LRD behaviors that exclude the drop of variance at the intermediate scales:
γ ( k ) = a ( 1 + ( k / q ) 2 M   ) ( 1 H ) / M
where a is the variance of the process, q is a scale parameter (in units of the scale k), M is the fractal parameter, and H is the Hurst parameter indicative of the LRD of the process, i.e., for 0.5 < H < 1 the process exhibits LRD behavior, while for 0 < H < 0.5 it exhibits an anti-persistent behavior, and for H = 0.5 a white-noise behavior. Here, the standardized climacogram is used, i.e., γ ^ _ ( k ) / γ ^ _ ( 1 ) , because the effect of the sample variance is already accounted for through the marginal fitting. We also note that a Gaussian process with q 0 and M = 0.5 coincides with the well-known fractional Gaussian noise model (e.g., [41]).

3. Data Extraction and Processing

For the analysis of the hourly PEV process, we use the reanalysis ensemble data extracted (access date at 29/10/2021; with coordinates S32-N42 and W115-E125) from the ERA5 [42] of the Centre for Medium-Range Weather Forecasts (ECMWF; https://cds.climate.copernicus.eu/ accessed on 1 October 2021) across California (Figure 1) and for the period 1979–today (Table 1).
Observations for the PE processes are usually available in monthly or daily resolutions and usually only for short periods, while a global gridded dataset based on the ERA5 data has been recently released [43]. Here, we use two datasets and compare the marginal and dependence structures of the reanalysis PEV timeseries with the PET timeseries of coarser monthly resolution, extracted from a network of 41 ground stations (see details in Table 1 and Figure 2 and Figure 3). Particularly, the monthly Penman–Monteith dataset of the CIMIS network is used, in which reference evapotranspiration and potential evapotranspiration coincide due to local surface and vegetation conditions. The samples of 41 meteorological stations (https://cimis.water.ca.gov/, accessed on 1 October 2021) are well-distributed across California (Figure 1) for the period 1983–2013 (Table 1), which corresponds to a maximum of 372 monthly values. The meteorological network has been developed in co-operation with Davis University, and the local environment of the meteorological stations allow accurate estimation of the PET.
To account for the impact of the double periodicity (diurnal and seasonal) of the PE processes on the dependence structure, we simulate the transformed process by applying a double standardization on the original timeseries. Particularly, we subtract the hourly and monthly means (Figure 4 and Figure 5) and then we divide with the hourly and monthly standard deviations (Figure 6 and Figure 7). Other transformation methods could be applied that take into consideration higher moments (e.g., [26]) such as skewness (Figure 8) and kurtosis (Figure 9) coefficients, or even more sophisticated ones [44]; however, as can be derived from Table 1 and Figure 6 and Figure 7, the PE processes (especially the aggregated PET process) is close to a light-tail distribution, and therefore we do not expect any significant differences by applying those methods. After the double standardization, we de-standardize each timeseries based on the total mean and standard deviation of the original timeseries (Table 1 and Figure 10). Finally, we fit the marginal and dependence models described in the previous section to each transformed timeseries, and the results are depicted and described in the next section.

4. Results

The PBF marginal distribution function is fitted to each transformed timeseries (e.g., Figure 11), and the parameters for all transformed timeseries can be seen in Table 2. Note that the fit of the PBF to all timeseries is exceptionally good. From Table 2, it can be observed that the transformed PEV and PET processes exhibit a light-tail behavior. The average values of the shape parameters are estimated as ξ ≈ 0.04 and ζ ≈ 5.7 for the CIMIS dataset, and ξ = 0.08 and ζ = 7.6 for the ERA5 transformed timeseries. It has been shown [25] that the tail index, ξ, does not depend on the averaging scale. Therefore, the slight differences in the estimated values are either due to statistical uncertainty or to differences in the nature of the data.
Additionally, the combined climacogram from all the empirical ones for the CIMIS transformed timeseries is depicted in Figure 12 and compared to the one from the ERA5 transformed timeseries depicted in Figure 13. It can be observed that a Hurst–Kolmogorov behavior is detected in both data sources, with a Hurst parameter of approximately 0.65. Specifically, the estimated parameters for the CIMIS dataset are H = 0.64 and q = 1.17 months (M is assumed to be 0.5 because the empirical climacogram is very close to an fGn process), and for the ERA5 timeseries they are H = 0.69, q = 19.7 h, and M = 0.8.

5. Discussion

Here we discuss how the above results can contribute to the existing literature relating to the potential evaporation and evapotranspiration from the point of view of stochastics and, in particular, of the HK dynamics.
The stochastic analysis of the potential evaporation (PEV) and potential evapotranspiration (PET) presented is useful (a) to highlight the stochastic similarities between them, (b) to quantify the variability and uncertainty of these processes, and (c) to develop a stochastic model capable of simulating important stochastic characteristics, for purposes such as forecasting and risk management. The PEV timeseries is extracted in hourly resolution as a reanalysis ensemble over California and through the ERA5 network, while for the PET, the high-quality CIMIS dataset with 41 stations is used over the same area for comparison.
The analysis of the above three tasks is performed based on the stochastic metrics and Hurst–Kolmgorov (HK) dynamics. Moreover, the marginal structures and second-order dependence structures are compared to the structures of each other and of other key hydrological-cycle processes such as temperature, relative humidity, wind speed, streamflow, and precipitation, as analyzed from a global network of stations in [25].
In particular, and similar to the global analysis, it is illustrated how the Pareto–Burr–Feller (PBF) probability distribution function may well describe the marginal structure of both the hourly PEV and monthly PET. Additionally, both processes are shown to exhibit a light-tail behavior. However, it is noted that the shape parameters of the PBF (i.e., ξ and ζ), which characterize the type of the tail, are slightly smaller in the CIMIS data (i.e., overall mean from stations 0.04 and 5.7, respectively) as compared to the reanalysis data (i.e., 0.08 and 7.6, respectively), indicating a heavier tail for the latter.
Additionally, it is found that, similarly to the other key hydrological-cycle processes mentioned above, both PEV and PET processes exhibit long-range dependence, with a Hurst parameter of medium strength. In particularly, H is estimated as 0.65 and 0.68 for the PET and PEV processes, respectively, which is weaker than the ones for temperature, relative humidity, solar radiation, and wind speed (0.80–0.85 [25]) and stronger than the one for precipitation (i.e., 0.61 [25]) for the examined range of scales spanning from the hourly resolution to the climatic scales. This can be interpreted as an indication that the PET and PEV processes have a wider predictability time window than precipitation’s, and narrower than the rest (i.e., entailing a higher degree of long-term unpredictability).
As a final remark, the need to apply a suitable stochastic model to reproduce important characteristics, such as LRD behavior, is stressed. The work shows the robust use of a stochastic framework to simulate the variability and uncertainty of a hydrometeorological process in emerging new practices and challenges:
  • Stochastic modelling of evapotranspiration at a fine time scale (e.g., hourly) is considered to be useful for numerous agronomist applications because it is strongly connected to the forecast of the plant water demands. In recent years of micro-farm techniques, the stochastic modelling of evapotranspiration, with sound physical-interpretation, has tracked the attention of the scientific community in order to simulate more accurately the water-food-energy nexus.
  • A proper stochastic model for the simulation of the evapotranspiration should be based at a wide range of spatio-temporal scales and meteorological conditions; thus, a global-scale analysis is important in order to identify stochastic similarities so as to improve the simulation techniques.
  • Stochastic simulation of the error analysis between the modelled and the measured Penman–Monteith assessment could highly contribute to improving potential evapotranspiration estimates.
  • Stochastic PET modeling could offer a solid probabilistic frame for identifying the long-term trend of hydrometeorological components in horizons greater than the available records and thus is of potential interest for climatological studies.

6. Conclusions

A stochastic model is presented for hourly potential evaporation (PEV) and monthly potential evapotranspiration (PET) based on the ERA5 hourly reanalysis data and the Penman–Monteith model applied to the well-known CIMIS network.
It was found that both the marginal probability distributions of PEV and PET are light-tailed when estimated through the Pareto–Burr–Feller distribution function. Additionally, the long-range dependence of both the PEV and PET is found to be of moderate strength, quantified through a Hurst parameter of 0.64 and 0.69, respectively.
The above results reveal the stochastic similarities between the ground and reanalysis data series. Additionally, the results are shown to be consistent to the hydrological-path of the marginal and dependence structures of Hurst–Kolmogorov dynamics. In particular, both PET and PEV can be placed between the stochastic structures of temperature, relative humidity, solar radiation, and wind speed (i.e., strong LRD and light- to medium-tail) and the precipitation’s structures (i.e., weak LRD and heavy tail). Finally, it is discussed how the existence of, even moderate, long-range dependence and tail distribution increase the variability and uncertainty of both processes, and thus limit their predictability.

Author Contributions

Conceptualization, P.D.; methodology, P.D. and D.K.; formal analysis, P.D.; investigation, P.D. and A.T.; data curation, A.T.; writing—review and editing, D.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Hersbach, H. et al., (2018) was downloaded from the Copernicus Climate Change Service (C3S) Climate Data Store. The results contain modified Copernicus Climate Change Service information 2020. Neither the European Commission nor ECMWF is responsible for any use that may be made of the Copernicus information or data it contains.

Acknowledgments

We thank both the anonymous reviewers for their suggestions and comments, and the editors for handling the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alexandris, S.; Proutsos, N. How significant is the effect of the surface characteristics on the Reference Evapotranspiration estimates? Agric. Water Manag. 2020, 237, 106181. [Google Scholar] [CrossRef]
  2. Doorenboos, J.; Pruitt, W.O. Guidelines for Predicting Crop Water Requirements. FAO Irrigation and Drainage Paper No.24; Food and Agricultural Organization of the United Nations: Rome, Italy, 1977. [Google Scholar]
  3. McMahon, T.A.; Finlayson, B.L.; Peel, M.C. Historical developments of models for estimating evaporation using standard meteorological data. Wiley Interdiscip. Rev. Water 2016, 3, 788–818. [Google Scholar] [CrossRef]
  4. Tegos, A.; Efstratiadis, A.; Koutsoyiannis, D. A Parametric Model for Potential Evapotranspiration Estimation Based on a Simplified Formulation of the Penman-Monteith Equation; InTech: Rijeka, Croatia, 2013; pp. 143–165. [Google Scholar]
  5. Tegos, A.; Malamos, N.; Koutsoyiannis, D. A parsimonious regional parametric evapotranspiration model based on a simplification of the Penman–Monteith formula. J. Hydrol. 2015, 524, 708–717. [Google Scholar] [CrossRef]
  6. Tegos, A.; Malamos, N.; Efstratiadis, A.; Tsoukalas, I.; Karanasios, A.; Koutsoyiannis, D. Parametric Modelling of Potential Evapotranspiration: A Global Survey. Water 2017, 9, 795. [Google Scholar] [CrossRef] [Green Version]
  7. Koutsoyiannis, D. HESS opinions. A random walk on water. Hydrol. Earth Syst. Sci. 2010, 14, 585–601. [Google Scholar] [CrossRef] [Green Version]
  8. Winters, P.R. Forecasting Sales by Exponentially Weighted Moving Averages. Manag. Sci. 1960, 6, 324–342. [Google Scholar] [CrossRef]
  9. Mohan, S.; Arumugam, N. Forecasting weekly reference crop evapotranspiration series. Hydrol. Sci. J. 1995, 40, 689–702. [Google Scholar] [CrossRef]
  10. Pandey, P.K.; Pandey, V.; Singh, R.; Bhakar, S.R. Stochastic Modelling of Actual Black Gram Evapotranspiration. J. Water Resour. Prot. 2009, 1, 448–455. [Google Scholar] [CrossRef] [Green Version]
  11. Rhenals, A.E.; Bras, R.L. The irrigation scheduling problem and evapotranspiration uncertainty. Water Resour. Res. 1981, 17, 1328–1338. [Google Scholar] [CrossRef]
  12. Uliana, E.M.; Da Silva, D.D.; Fraga, M.D.S.; Lisboa, L. Estimate of reference evapotranspiration through continuous probability modelling. Eng. Agríc. 2017, 37, 257–267. [Google Scholar] [CrossRef] [Green Version]
  13. Khanmohammadi, N.; Rezaie, H.; Montaseri, M.; Behmanesh, J. Regional probability distribution of the annual reference evapotranspiration and its effective parameters in Iran. Theor. Appl. Climatol. 2018, 134, 411–422. [Google Scholar] [CrossRef]
  14. Pouliasis, G.; Torres-Alves, G.A.; Morales-Napoles, O. Stochastic Modeling of Hydroclimatic Processes Using Vine Copulas. Water 2021, 13, 2156. [Google Scholar] [CrossRef]
  15. Pham, M.T.; Vernieuwe, H.; De Baets, B.; Willems, P.; Verhoest, N.E.C. Stochastic simulation of precipitation-consistent daily reference evapotranspiration using vine copulas. Stoch. Environ. Res. Risk Assess. 2016, 30, 2197–2214. [Google Scholar] [CrossRef] [Green Version]
  16. Pham, M.T.; Vernieuwe, H.; De Baets, B.; Verhoest, N.E.C. A coupled stochastic rainfall–evapotranspiration model for hydrological impact analysis. Hydrol. Earth Syst. Sci. 2018, 22, 1263–1283. [Google Scholar] [CrossRef] [Green Version]
  17. Zhang, Y.; Liu, C.; Tang, Y.; Yang, Y. Trends in pan evaporation and reference and actual evapotranspiration across the Tibetan Plateau. J. Geophys. Res. Atmos. 2007, 112. [Google Scholar] [CrossRef]
  18. Stefanidis, S.; Alexandridis, V. Precipitation and Potential Evapotranspiration Temporal Variability and Their Relationship in Two Forest Ecosystems in Greece. Hydrology 2021, 8, 160. [Google Scholar] [CrossRef]
  19. Hobbins, M.; Ramírez, J.A.; Brown, T.C. Trends in pan evaporation and actual evapotranspiration across the conterminous U.S.: Paradoxical or complementary? Geophys. Res. Lett. 2004, 31. [Google Scholar] [CrossRef] [Green Version]
  20. Cong, Z.T.; Yang, D.W.; Ni, G.H. Does evaporation paradox exist in China? Hydrol. Earth Syst. Sci. 2009, 13, 357–366. [Google Scholar] [CrossRef] [Green Version]
  21. Tegos, A.H.; Tyralis, D.; Koutsoyiannis, D.; Hamed, K.H. An R function for the estimation of trend signifcance under the scaling hypothesis—Application in PET parametric annual time series. Open Water J. 2017, 4, 66–71. [Google Scholar]
  22. Koutsoyiannis, D. Revisiting the global hydrological cycle: Is it intensifying? Hydrol. Earth Syst. Sci. 2020, 24, 3899–3932. [Google Scholar] [CrossRef]
  23. Iliopoulou, T.; Koutsoyiannis, D. Projecting the future of rainfall extremes: Better classic than trendy. J. Hydrol. 2020, 588, 125005. [Google Scholar] [CrossRef]
  24. Koutsoyiannis, D.; Dimitriadis, P.; Lombardo, F.; Stevens, S.; Tsonis, A.A. From Fractals to Stochastics: Seeking Theoretical Consistency in Analysis of Geophysical Data. In Advances in Nonlinear Geosciences; Springer International Publishing: Berlin/Heidelberg, Germany, 2017; pp. 237–278. [Google Scholar]
  25. Koutsoyiannis, D. Stochastics of Hydroclimatic Extremes—A Cool Look at Risk; National Technical University of Athens: Athens, Greece, 2021; 330p. [Google Scholar]
  26. Dimitriadis, P.; Koutsoyiannis, D.; Iliopoulou, T.; Papanicolaou, P. A Global-Scale Investiga-tion of Stochastic Similarities in Marginal Distribution and Dependence Structure of Key Hydrological-Cycle Processes. Hydrology 2021, 8, 59. [Google Scholar] [CrossRef]
  27. Glynis, K.T.; Iliopoulou, P.; Dimitriadis, D. Koutsoyiannis, Stochastic investigation of daily air temperature extremes from a global ground station network. Stoch. Environ. Res. Risk Assess. 2021, 35, 1585–1603. [Google Scholar] [CrossRef]
  28. Katikas, L.; Dimitriadis, P.; Koutsoyiannis, D.; Kontos, T.; Kyriakidis, P. A stochastic simulation scheme for the long-term persistence, heavy-tailed and double periodic behavior of observational and reanalysis wind time-series. Appl. Energy 2021, 295, 116873. [Google Scholar] [CrossRef]
  29. Koudouris, G.P.; Dimitriadis, T.; Iliopoulou, N.; Mamassis, N.; Koutsoyiannis, D. A stochastic model for the hourly solar radiation process for application in renewable resources management. Adv. Geosci. 2018, 45, 139–145. [Google Scholar] [CrossRef] [Green Version]
  30. Hurst, H.E. Long-Term Storage Capacity of Reservoirs. Trans. Am. Soc. Civ. Eng. 1951, 116, 770–799. [Google Scholar] [CrossRef]
  31. Kolmogorov, A.N. Wiener Spirals and Some Other Interesting Curves in a Hilbert Space. Selected Works of Kolmogorov, A.N.; Tikhomirov, V.M., Ed.; Wolters Kluwer: Dordrecht, The Netherlands, 1991; pp. 303–307. [Google Scholar]
  32. Koutsoyiannis, D. Hurst–Kolmogorov dynamics as a result of extremal entropy production. Phys. A Stat. Mech. Appl. 2011, 390, 1424–1432. [Google Scholar] [CrossRef]
  33. Koutsoyiannis, D. Hurst-Kolmogorov dynamics and uncertainty. J. Am. Water Resour. Assoc. 2011, 47, 481–495. [Google Scholar] [CrossRef]
  34. Dimitriadis, P.; Koutsoyiannis, D. Climacogram versus autocovariance and power spectrum in stochastic modelling for Markovian and Hurst–Kolmogorov processes. Stoch. Environ. Res. Risk Assess. 2015, 29, 1649–1669. [Google Scholar] [CrossRef]
  35. Dimitriadis, P.; Iliopoulou, T.; Sargentis, G.-F.; Koutsoyiannis, D. Spatial Hurst–Kolmogorov Clustering. Encyclopedia 2021, 1, 77. [Google Scholar] [CrossRef]
  36. Singh, S.K.; Maddala, G.S. A Function for Size Distribution of Incomes: Reply. Economic 1978, 46, 461. [Google Scholar] [CrossRef]
  37. Burr, I.W. Cumulative Frequency Functions. Ann. Math. Stat. 1942, 13, 215–232. [Google Scholar] [CrossRef]
  38. Feller, W. Law of large numbers for identically distributed variables. Introd. Probab. Theory Appl. 1971, 2, 231–234. [Google Scholar]
  39. Beran, J. Statistical Methods for Data with Long-Range Dependence. Stat. Sci. 1992, 7, 404–416. [Google Scholar]
  40. Dimitriadis, P.; Koutsoyiannis, D. Stochastic synthesis approximating any process dependence and distribution. Stoch. Environ. Res. Risk Assess. 2018, 32, 1493–1515. [Google Scholar] [CrossRef]
  41. Mandelbrot, B.B.; Van Ness, J.W. Fractional Brownian motions, fractional noises and applications. J. Soc. Ind. Appl. Math. 1968, 10, 422–437. [Google Scholar] [CrossRef]
  42. Hersbach, H.; Bell, B.; Berrisford, P.; Biavati, G.; Horányi, A.; Muñoz-Sabater, J.; Nicolas, J.; Peubey, C.; Radu, R.; Rozum, I.; et al. ERA5 Hourly Data on Single Levels from 1979 to Present; Copernicus Climate Change Service (C3S) and Climate Data Store (CDS): Bonn, Germany, 2018. [Google Scholar] [CrossRef]
  43. Singer, M.B.; Asfaw, D.T.; Rosolem, R.; Cuthbert, M.O.; Miralles, D.G.; MacLeod, D.; Quichimbo, E.A.; Michaelides, K. Hourly potential evapotranspiration at 0.1° resolution for the global land surface from 1981-present. Sci. Data 2021, 8, 224. [Google Scholar] [CrossRef]
  44. Vavoulogiannis, S.; Iliopoulou, T.; Dimitriadis, P.; Koutsoyiannis, D. Multiscale temporal irreversibility of streamflow and its stochastic modelling. Hydrology 2021, 8, 63. [Google Scholar] [CrossRef]
Figure 1. Map of Locations of the selected stations.
Figure 1. Map of Locations of the selected stations.
Hydrology 08 00177 g001
Figure 2. Total recorded and missing values of the PET timeseries for each station.
Figure 2. Total recorded and missing values of the PET timeseries for each station.
Hydrology 08 00177 g002
Figure 3. Minimum and maximum values of the PET timeseries for each station.
Figure 3. Minimum and maximum values of the PET timeseries for each station.
Hydrology 08 00177 g003
Figure 4. Monthly means of the PET timeseries for each station.
Figure 4. Monthly means of the PET timeseries for each station.
Hydrology 08 00177 g004
Figure 5. Hourly means of the PEV timeseries for each month.
Figure 5. Hourly means of the PEV timeseries for each month.
Hydrology 08 00177 g005
Figure 6. Monthly standard deviations of the PET timeseries for each station.
Figure 6. Monthly standard deviations of the PET timeseries for each station.
Hydrology 08 00177 g006
Figure 7. Hourly standard deviations of the PEV timeseries for each month.
Figure 7. Hourly standard deviations of the PEV timeseries for each month.
Hydrology 08 00177 g007
Figure 8. Monthly skewness coefficients of the PET timeseries for each station.
Figure 8. Monthly skewness coefficients of the PET timeseries for each station.
Hydrology 08 00177 g008
Figure 9. Skewness and kurtosis coefficients of the PET timeseries for all stations.
Figure 9. Skewness and kurtosis coefficients of the PET timeseries for all stations.
Hydrology 08 00177 g009
Figure 10. Observed and transformed PET timeseries at the Davis station.
Figure 10. Observed and transformed PET timeseries at the Davis station.
Hydrology 08 00177 g010
Figure 11. Observed and theoretical results of the PBF marginal distribution function of the PET transformed timeseries at the Davis station.
Figure 11. Observed and theoretical results of the PBF marginal distribution function of the PET transformed timeseries at the Davis station.
Hydrology 08 00177 g011
Figure 12. Observed and theoretical climacograms through the HK model for all the available PET transformed timeseries adjusted for bias, with the 25% and 75% quantiles (note that the coefficient of determination for the model is R2 = 0.993).
Figure 12. Observed and theoretical climacograms through the HK model for all the available PET transformed timeseries adjusted for bias, with the 25% and 75% quantiles (note that the coefficient of determination for the model is R2 = 0.993).
Hydrology 08 00177 g012
Figure 13. Observed and theoretical climacograms through the HK model for the ERA5 transformed timeseries adjusted for bias (note that the coefficient of determination for the model is R2 = 0.997).
Figure 13. Observed and theoretical climacograms through the HK model for the ERA5 transformed timeseries adjusted for bias (note that the coefficient of determination for the model is R2 = 0.997).
Hydrology 08 00177 g013
Table 1. Information on the selected stations and the reanalysis data.
Table 1. Information on the selected stations and the reanalysis data.
Sequence NumberNameProcessTemporal ResolutionTime PeriodNumber of Data ValuesMean (mm)Standard Deviation (mm)Skewness Coefficient
1Five PointsPETmonthly1982–2013363131.573.70.0
2DavisPETmonthly1982–2013372120.668.70.0
3Firebaugh TelesPETmonthly1982–2013370118.168.90.1
4GerberPETmonthly1982–2013370117.367.90.1
5DurhamPETmonthly1982–2013369107.861.70.1
6CarminoPETmonthly1982–2013369116.868.80.3
7StratfordPETmonthly1982–2013369128.275.40.0
8CastorvillePETmonthly1982–201336879.932.00.1
9KettlemanPETmonthly1982–2013368130.473.90.0
10BishopPETmonthly1983–2013363125.560.90.0
11ParlierPETmonthly1983–2013362112.566.00.1
12CalipatriaPETmonthly1983–2013360151.265.2−0.1
13Mc_ArthurPETmonthly1983–2013357101.266.20.2
14UC_RiversidePETmonthly1985–2013337121.947.00.1
15BrentwoodPETmonthly1985–2013327115.868.10.1
16San_Luis_ObispoPETmonthly1986–2013327107.539.5−0.1
17Blackwells_cornerPETmonthly1987–2013321128.973.10.2
18Los BanosPETmonthly1988–2013301119.870.40.1
19BuntigvillePETmonthly1986–2013325112.967.80.1
20TemeculaPETmonthly1986–2013320113.539.90.0
21Santa_YnezPETmonthly1986–2013320105.146.30.0
22SeeleyPETmonthly1987–2013314159.769.1−0.1
23MantecaPETmonthly1987–2013308109.764.70.1
24ModestoPETmonthly1987–2013312110.764.90.1
25IrvinePETmonthly1987–2013309105.039.40.1
26OakvillePETmonthly1989–2013292103.855.50.0
27PomonaPETmonthly1989–2013291103.444.70.1
28Frenso_StatePETmonthly1988–2013297117.771.20.1
29Santa_RosaPETmonthly1990–201328293.950.90.0
30Browns_ValleyPETmonthly1989–2013291112.265.40.1
31LindcovePETmonthly1989–2013290110.465.90.1
32MelolandPETmonthly1989–2013283153.366.5−0.1
33AlturasPETmonthly1989–201329197.060.70.3
34CuyamaPETmonthly1989–2013289128.461.40.1
35TulelakePETmonthly1990–201329196.460.60.2
36Goleta_foothills *PETmonthly1990–201319799.134.80.0
37WindsorPETmonthly1990–201326696.453.60.1
38De_LaveagaPETmonthly1990–201327488.639.4−0.1
39WestlandsPETmonthly1992–2013255131.276.00.0
40Sanel_ValleyPETmonthly1990–2013269107.262.80.1
41Santa_MonicaPETmonthly1993–201324699.134.90.0
42CIMIS (overall)PETmonthly1983–201312985114.463.50.2
44ERA5PEVhourly1979–20210.93 × 1060.080.111.5
* There is a large gap in timeseries from 1995–2001.
Table 2. Parameters of the marginal probability distribution function for all transformed timeseries of each station (note that the squared correlation coefficient is R2 > 0.99 for all models). The symbols ξ, ζ, λ, and d correspond to Equation (1).
Table 2. Parameters of the marginal probability distribution function for all transformed timeseries of each station (note that the squared correlation coefficient is R2 > 0.99 for all models). The symbols ξ, ζ, λ, and d correspond to Equation (1).
Sequence NumberNameξζλ (mm)d (mm)
1Five Points0.1004.5240.0–105.0
2Davis0.0949.2355.1–236.0
3Firebaugh Teles0.0718.0353.2–227.5
4Gerber0.0788.2349.6–227.9
5Durham0.0636.3285.2–167.6
6Carmino0.0345.8326.1–191.1
7Stratford0.0548.0427.1–286.5
8Castorville0.0495.4132.7–45.7
9Kettleman0.0184.7308.4–154.0
10Bishop0.06710.4304.9–171.3
11Parlier0.0426.4326.8–199.4
12Calipatria0.0937.7285.2–131.2
13Mc_Arthur0.0386.8340.7–223.3
14UC_Riverside0.0714.2145.7–14.8
15Brentwood0.0727.4344.0–221.1
16San_Luis_Obispo0.0774.6138.0–25.1
17Blackwells_corner0.0015.1352.9–196.8
18Los Banos0.0567.2367.4–235.1
19Buntigville0.0256.0327.6–193.9
20Temecula0.0745.1145.7–26.0
21Santa_Ynez0.0125.1199.9–78.5
22Seeley0.0857.7285.3–116.6
23Manteca0.0464.3233.0–107.6
24Modesto0.0133.7231.3–100.0
25Irvine0.0314.1132.3–15.5
26Oakville0.0023.5188.2–65.3
27Pomona0.0256.0208.3–91.1
28Frenso_State0.0313.3210.4–72.4
29Santa_Rosa0.0263.6169.3–61.1
30Browns_Valley0.0044.4279.5–143.7
31Lindcove0.0456.2315.7–190.4
32Meloland0.0295.2268.0–94.7
33Alturas0.0195.0259.4–142.3
34Cuyama0.0306.8343.4–197.5
35Tulelake0.0225.2262.2–146.8
36Goleta_foothills0.0174.7129.3–18.5
37Windsor0.0013.1166.1–51.4
38De_Laveaga0.0015.5195.3–90.9
39Westlands0.0183.2230.2–76.6
40Sanel_Valley0.0016.6385.2–252.7
41Santa_Monica0.0164.7144.3–33.4
42CIMIS (meanl)0.0405.7260.8–132.4
43ERA5-PEV0.0767.60.63–0.54
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dimitriadis, P.; Tegos, A.; Koutsoyiannis, D. Stochastic Analysis of Hourly to Monthly Potential Evapotranspiration with a Focus on the Long-Range Dependence and Application with Reanalysis and Ground-Station Data. Hydrology 2021, 8, 177. https://doi.org/10.3390/hydrology8040177

AMA Style

Dimitriadis P, Tegos A, Koutsoyiannis D. Stochastic Analysis of Hourly to Monthly Potential Evapotranspiration with a Focus on the Long-Range Dependence and Application with Reanalysis and Ground-Station Data. Hydrology. 2021; 8(4):177. https://doi.org/10.3390/hydrology8040177

Chicago/Turabian Style

Dimitriadis, Panayiotis, Aristoteles Tegos, and Demetris Koutsoyiannis. 2021. "Stochastic Analysis of Hourly to Monthly Potential Evapotranspiration with a Focus on the Long-Range Dependence and Application with Reanalysis and Ground-Station Data" Hydrology 8, no. 4: 177. https://doi.org/10.3390/hydrology8040177

APA Style

Dimitriadis, P., Tegos, A., & Koutsoyiannis, D. (2021). Stochastic Analysis of Hourly to Monthly Potential Evapotranspiration with a Focus on the Long-Range Dependence and Application with Reanalysis and Ground-Station Data. Hydrology, 8(4), 177. https://doi.org/10.3390/hydrology8040177

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