Next Article in Journal
Energy Dissipation Structures: Influence of Aeration in Supercritical Flows
Previous Article in Journal
Evaluation of Extreme Dry and Wet Conditions Using Climate and Hydrological Indices in the Upper Part of the Gallikos River Basin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Copula Bias Correction for Extreme Precipitation in Reanalysis Data over a Greek Catchment †

by
Georgia Lazoglou
1,*,
Christina Anagnostopoulou
1,
Charalampos Skoulikaris
2 and
Konstantia Tolika
1
1
Department of Meteorology Climatology School of Geology, Aristotle University of Thessaloniki, 541 24 Thessaloniki, Greece
2
Department of Civil Engineering, Aristotle University of Thessaloniki, 541 24 Thessaloniki, Greece
*
Author to whom correspondence should be addressed.
Presented at the 3rd International Electronic Conference on Water Sciences, 15–30 November 2018; Available online: https://ecws-3.sciforum.net.
Proceedings 2019, 7(1), 4; https://doi.org/10.3390/ECWS-3-05817
Published: 15 November 2018
(This article belongs to the Proceedings of ECWS-3)

Abstract

:
The projection of extreme precipitation events with higher accuracy and reliability that engender severe socioeconomic impacts more frequently is considered a priority research topic in the scientific community. Although large-scale initiatives for monitoring meteorological and hydrological variables exist, the lack of data is still evident particularly in regions with complex topographic characteristics. The latter results in the use of reanalysis data or data derived from regional climate models, however both datasets are biased to the observations resulting in nonaccurate results in hydrological studies. The current research presents a newly developed statistical method for the bias correction of the maximum rainfall amount at watershed scale. In particular, the proposed approach necessitates the coupling of a spatial distribution method, namely Thiessen polygons, with a multivariate probabilistic distribution method, namely copulas, for the bias correction of the maximum precipitation. The case study area is the Nestos River basin where the several extreme episodes that have been recorded have direct impacts to the regional agricultural economy. Thus, using daily data by three monitoring stations and daily reanalysis precipitation values from the grids closest to these stations, the results demonstrated that the bias corrected maximum precipitation totals (greater than 90%) is much closer to the real max precipitation totals, while the respective reanalysis value underestimates the real precipitation totals. The overall improvement of the output shows that the proposed Thiessen-copula method could constitute a significant asset to hydrologic simulations.

Published: 15 November 2018

1. Introduction

The extreme precipitation episodes of recent decades have been observed more frequently and their intense is greater compared with the past [1]. The impacts of these extremes are obvious in many different fields such as economy, society, agriculture and hydrology, resulting in the need of reliable projections. According to Mao et al. [2], in order to achieve higher accuracy in hydrological results, it is mandatory to correct the bias between a model’s used real values of several climate parameters and especially precipitations. In cases where a lack of observed data has occurred, but reliable techniques such as the use of reanalysis data, is proposed as scientifically proven solution. Bastola and Misra [3] in particular demonstrate that reanalysis data is useful in simulating realistic hydrological response at watershed scales. At the same time, although precipitation estimates from global reanalysis are dynamically consistent with large-scale circulation, when precipitation is compared with rain gauge estimates the output is poor, since reanalysis products are forecast by the reanalysis system and precipitation is not assimilated [4].
Numerous studies in the fields of insurance and finance have attempted using the copula method [5,6]. The utility of the copula method based on its ability to analyze the dependence of two or more random variables that have not necessarily the same distribution [7]. Additionally, copulas have the advantage to capture the features of the dependence [8] and to examine this dependence not only linearly, as do other indices [2]. As such, copulas have been used widely in hydrology of late. Shiau [9] suggest copulas for drought analysis in order to overcome the problem that drought characteristics have different distributions. Several scientists use copulas for analyzing drought characteristics (e.g., severity and intense) [10,11] or for correlating drought with other climate parameters such as precipitation [12]. Furthermore Golian [13] used a bivariate copula function for studying the rainfall–runoff simulations of a watershed in the region of Iran, while Perera [14] used this method for studying the interdependence between the Kelani River and Kotte Canal in Sri Lanka.
The present study investigates the combination of the copula’s probabilistic distribution method with the Thiessen polygon spatial distribution method to tackle bias correction of extreme precipitations reanalysis data in the important hydrological region of Nestos river basin in Greece. Thiessen polygons have been widely used in hydrology for spatial interpolation since the work of Lee and Schachter in 1980 [15]. The evaluation of the results was achieved using statistical tools such as Taylor diagrams as well as the relative operating characteristics curves (ROC). The latter are popular in clinical epidemiology as they test the accuracy of a diagnosis [16] and in this study is checked the accuracy of the bias corrected values.

2. Materials and Methods

2.1. Data

The present study uses daily precipitation data from four meteorological stations located in the Nestos catchment (Figure 1a). Apart from the observed precipitation records, ERA-Interim reanalysis data with spatial resolution of 12.5 × 12.5 km, were retrieved by the European Centre for Medium-Range Weather Forecasts (ECMWF) for the selected case study area. Thus, for every station the closest continental grid point that presented similar topographic characteristics was selected. Both reanalysis data and observed records cover a time period of nine sequential years, i.e., from 1987 to 1995.

2.2. Methodology

According to Nelsen [17], copulas are multivariate cumulative distribution functions with the ability to mathematically model the dependence between two or more variables using their marginal distributions. Assuming that X and Y are two random variables with F and G marginal distributions respectively, the joint distribution function of the initial variables is H and is equal to the copula function of their marginal (H(x,y) = C(F(x), G(y)). The centered theorem of copulas is the Sklar’s theorem [18]. According to that, if the marginals of the studied variables are continuous, then the copula function can uniquely be defined. Otherwise, C is unique on RanF × RanG (Ran is the range). The importance of this theorem is that every joint distribution function can be decomposed into the marginal of the variables and into a copula function, which completely describes their dependence.
In the present study, the copula method was combined with thiessen triangles, which is an alteration of the thiessen polygons method, to achieve a bias correction of total extreme precipitation between real observations and reanalysis data. Firstly, three (stations 1, 2 and 4) of the four available stations were used for the analysis while the other one (station 3) for evaluation. Specifically, the three selected stations formed a triangle that includes the tested station (Figure 1b).
For the stations located at the triangle vertices, the absolute maximum and the monthly precipitations were used to model the dependence among them. Twelve copula families (Gaussian, Student t, Clayton, Gumbel, Frank, Joe, BB1, BB6, BB7, BB8, Tawn type 1 and 2) coming from both Archimedean and elliptical categories were tested in order to select the one that can best describe the dependence. The final selection was based on Akaike’s information criterion [19] and Bayesian information criterion. Thereafter, using the copula results by the stations located at the triangle vertices and taking also into account the distance between them and the evaluated stations (x-point), a new copula family was defined (Figure 1c). The newly defined copula family can mathematically describe the dependence between the mean and maximum precipitation at the x-point. Consequently, using the new copula family and the reanalysis mean precipitation data, the maximum extreme bias corrected precipitations have been calculated.

3. Results

3.1. Total and Extreme Precipitation Analysis

The dependence between absolute maximum extreme and monthly precipitation was assessed using copulas to three of the four studied stations. After testing several copula families, it was found that the Survival Clayton family reflects the dependence for both station 1 and station 4, while Frank copula reflects the dependence for station 2. The dependence is more power for the Prasinada station as the Kendall’s tau correlation index is almost 0.87, while it is almost 0.78 for the other two stations. Figure 2 visualizes that relationship, confirming the fact that Survival Clayton copula presents upper-tail dependence (stations 1 and 4) while Frank presents no upper or lower tail dependence (station 2).

3.2. Bias Correction Results

Reanalysis data presents important biases from real observations especially for precipitation parameters and extreme events. As a consequence, the main purpose of this study is to reduce extreme rainfall biases between real and reanalysis data at the scale of a hydrological basin. Table 1 presents the observed reanalysis and the bias corrected extreme precipitation indices for the station of Sidironero. Extreme precipitation indices are defined as the 90th, 95th and 99th percentile of the monthly precipitation [20].
As it can be seen from Table 1, for all indices the bias corrected values are closer to the observed ones compared with reanalysis data. Additionally, as the percentile became higher, reanalysis presents higher differences compared to the observations, while the bias corrected were closer to the observed ones. This is also proved by the Taylor diagram (Figure 3a) of the studied data sets. In particular, the Taylor diagram shows that the correlation between observed and reanalysis data was almost zero while after the bias correction the correlation increased to 0.5. Additionally, the root mean square error has been reduced while there is an increase of the variation.
An additional evaluation of the results was conducted with the use of the ROC curves. According to Figure 3b the area under curve—which is an effective measure of accuracy [21]—is bigger for the bias corrected values proving that the Thiessen-copula method can be used for the bias correction of extreme precipitation.

4. Discussion—Conclusions

The bias correction of extreme rainfall with the coupling of the copula method and an alteration of the Thiessen polygon method is presented. The proposed method adjusts the extreme reanalysis precipitation data to observed data in a selected river basin in Greece. The method’s evaluation was achieved after the comparison of the bias corrected values with observed datasets using different statistical and optional methods.
Lafon et al. [22], Teutschbein and Seibert [23] have studied and compared different bias correction methods in hydrology. Methods such as the linear scaling approach, the delta change method [24] or local intensity scaling [25] use simple mathematical equations for bias correction. However, they mainly focus on estimating mean values without expanding to the whole distribution [2].
Additionally, as Yang et al. [26] mention, the accuracy for extremes is much lower even with more dynamic methods. With the same concept, Berg, et al. [27], indicate the importance of bias correction to reanalysis products, regarding soil moisture, runoff, and snow water equivalence at simulations covering the geographic area of North America.
The results of the present study show that the copula method combining with Thiessen triangles can be an accurate tool for rainfall extremes bias correction. The copula method also presented satisfying results in Mao’s [2] study, wherein it was used for the bias correction of the model’s precipitation values in Germany. Additionally, Piani [28] propose the use of copula for a two-dimensional bias correction for the parameters of temperature and precipitation in climate models. In accordance with this study, the success of the method derived from the ability of the copula method to satisfactorily represent the dependence structure of the studied variables. Consequently, the dependence structure is not the same in every station in a specific region, as it is also observed in the three studied stations of the presented case study. Furthermore the bias corrected extreme precipitations not only were much closer to the real ones but they also have higher correlation with the observed extremes, as well as the root mean square error is lower compared with the reanalysis data. The need of bias correction is more obvious in the case of climate change, where the models outputs need to be bias-corrected before it can be used for climate change impact studies. In that case, various methods such as the quantile mapping, the cumulative distribution function transform (CDF-t), and equidistant quantile matching are presented in the literature [29]. In conclusion, this investigation proposes the copula method in combination with the Thiessen triangles technique as a useful tool for the bias correction of extreme precipitations. It is also believed that this method would be a fruitful area for further application in large river basins.

Author Contributions

All authors conceived and initiated the study. G.L. analyzed the data. All authors contributed to the discussion and interpretation of the results and the writing of the manuscript.

Acknowledgments

This work was supported by a grant by the State Scholarship Foundation (IKY) in Greece.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ROCRelative Operating Characteristics curves

References

  1. Pachauri, R.K.; Allen, M.R.; Barros, V.R.; Broome, J.; Cramer, W.; Christ, R.; Dubash, N.K. Climate Change 2014: Synthesis Report; Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC), Ed.; IPCC: Geneva, Switzerland, 2014. [Google Scholar]
  2. Mao, G.; Vogl, S.; Laux, P.; Wagner, S.; Kunstmann, H. Stochastic bias correction of dynamically downscaled precipitation fields for Germany through Copula-based integration of gridded observation data. Hydrol. Earth Syst. Sci. 2015, 19, 1787–1806. [Google Scholar] [CrossRef]
  3. Bastola, S.; Misra, V. Evaluation of dynamically downscaled reanalysis precipitation data for hydrological application. Hydrol. Process. 2014, 28, 1989–2002. [Google Scholar] [CrossRef]
  4. Gehne, M.; Hamill, T.M.; Kiladis, G.N.; Trenberth, K.E. Comparison of Global Precipitation Estimates across a Range of Temporal and Spatial Scales. J. Clim. 2016, 29, 7773–7795. [Google Scholar] [CrossRef]
  5. Giacomini, E. Risk Management with Copula. Master’s Thesis, Hulmboldt University, Berlin, Germany, 2005. [Google Scholar]
  6. Malevergne, Y.; Sornette, D. Testing the Gaussian copula hypothesis for financial assets dependences. Quant. Financ. 2003, 3, 231–250. [Google Scholar] [CrossRef]
  7. Zhang, L.; Singh, V.P. Bivariate rainfall frequency distributions using Archimedean copulas. J. Hydrol. 2007, 332, 93–109. [Google Scholar] [CrossRef]
  8. Canela, M.A.; Collazo, E.P. Modelling dependence in Latin American markets using copula functions. J. Emerg. Mark. Financ. 2012, 11, 231–270. [Google Scholar] [CrossRef]
  9. Shiau, J.T. Fitting drought duration and severity with two-dimensional copulas. Water Resour. Manag. 2006, 20, 795–815. [Google Scholar] [CrossRef]
  10. Yusof, F.; Hui-Mean, F.; Suhaila, J.; Yusof, Z. Characterisation of drought properties with bivariate copula analysis. Water Resour. Manag. 2013, 27, 4183–4207. [Google Scholar] [CrossRef]
  11. Mirabbasi, R.; Fakheri-Fard, A.; Dinpashoh, Y. Bivariate drought frequency analysis using the copula method. Theor. Appl. Climatol. 2012, 108, 191–206. [Google Scholar] [CrossRef]
  12. Fan, L.; Wang, H.; Wang, C.; Lai, W.; Zhao, Y. Exploration of Use of Copulas in Analysing the Relationship between Precipitation and Meteorological Drought in Beijing, China. Adv. Meteorol. 2017. [Google Scholar] [CrossRef]
  13. Golian, S.; Saghafian, B.; Farokhnia, A. Copula-based interpretation of continuous rainfall–runoff simulations of a watershed in northern Iran. Can. J. Earth Sci. 2012, 49, 681–691. [Google Scholar] [CrossRef]
  14. Perera, K.; Sumathipala, P.L.N.S.; Wikramanayake, P.N. Interdependence between the Water Levels of Kotte Canal and Kelani River using Copulas. Sri Lankan J. Appl. Stat. 2015, 16, 146. [Google Scholar] [CrossRef]
  15. Lee, D.T.; Schachter, B.J. Two algorithms for constructing a Delaunay triangulation. Int. J. Comput. Inf. Sci. 1980, 9, 219–242. [Google Scholar] [CrossRef]
  16. Metz, C.E. ROC methodology in radiologic imaging. Investig. Radiol. 1986, 21, 720–733. [Google Scholar] [CrossRef]
  17. Nelsen, R.B. An Introduction to Copulas; Springer: New York, NY, USA, 1999. [Google Scholar]
  18. Sklar, A. Fonctions de Répartition À N Dimensions Et Leurs Marge; Publications de l’Institut de Statistique de L’Universit de Paris: Paris, France, 1959; Volume 8, pp. 229–231. [Google Scholar]
  19. Akaike, H. Information theory and an extension of the maximum likelihood principle. In Second International Symposium on Information Theory; Springer: New York, NY, USA, 1973; pp. 267–281. [Google Scholar]
  20. WMO. Available online: https://www.wmo.int/pages/prog/wcp/ccl/opace/opace2/documents/DraftversionoftheGuidelineso ntheDefinitionandMonitoringofExtremeWeatherandClimateEvents.pdf (accessed on 4 April 2016).
  21. Hanley, J.A.; McNeil, B.J. The meaning and use of the area under a receiver operating characteristic (ROC) curve. Radiology 1982, 143, 29–36. [Google Scholar] [CrossRef]
  22. Lafon, T.; Dadson, S.; Buys, G.; Prudhomme, C. Bias correction of daily precipi-tation simulated by a regional climate model: A comparison of methods. Int. J. Climatol. 2013, 33, 1367–1381. [Google Scholar] [CrossRef]
  23. Teutschbein, C.; Seibert, J. Bias correction of regional climate model simulations for hydrological climate-change impact studies: Review and evaluation of different meth-ods. J. Hydrol. 2012, 456–457, 12–29. [Google Scholar] [CrossRef]
  24. Hay, L.E.; Wilby, R.L.; Leavesley, G.H. A Comparison of Delta Change and Downscaled GCM Scenarios for Three Mountainous Basins in the United States. J. Am. Water Resour. Assoc. 2000, 36, 387–397. [Google Scholar] [CrossRef]
  25. Schmidli, J.; Frei, C.; Vidale, P.L. Downscaling from GCM precipitation: A bench- mark for dynamical and statistical downscaling methods. Int. J. Climatol. 2006, 26, 679–689. [Google Scholar] [CrossRef]
  26. Yang, W.; Andréasson, J.; Phil Graham, L.; Olsson, J.; Rosberg, J.; Wetterhall, F. Distribution-based scaling to improve usability of regional climate model projections for hydrological climate change impacts studies. Hydrol. Res. 2010, 41, 211. [Google Scholar] [CrossRef]
  27. Berg, A.A.; Famiglietti, J.S.; Walker, J.P.; Houser, P.R. Impact of bias correction to reanalysis products on simulations of North American soil moisture and hydrological fluxes. J. Geophys. Res. Atmos. 2003, 108, 4490. [Google Scholar] [CrossRef]
  28. Piani, C.; Haerter, J.O. Two dimensional bias correction of temperature and precipitation copulas in climate models. Geophys. Res. Lett. 2012, 39. [Google Scholar] [CrossRef]
  29. Pierce, D.W.; Cayan, D.R.; Maurer, E.P.; Abatzoglou, J.T.; Hegewisch, K.C. Improved Bias Correction Techniques for Hydrological Simulations of Climate Change. J. Hydrometeorol. 2015, 16, 2421–2442. [Google Scholar] [CrossRef]
Figure 1. (a) The map of Greece—the red circle includes the studied region (b) zoom of the Nestos region and the four studied stations with the studied triangle (1 = Achladia, 2 = Prasinada, 3 = Sidironero, 4 = Toxotes), (c) St1 to St3 indicate the three stations at the triangle vertices. The x-point is the unknown station and dist1 to dist3 are the distances between x-point and the three stations.
Figure 1. (a) The map of Greece—the red circle includes the studied region (b) zoom of the Nestos region and the four studied stations with the studied triangle (1 = Achladia, 2 = Prasinada, 3 = Sidironero, 4 = Toxotes), (c) St1 to St3 indicate the three stations at the triangle vertices. The x-point is the unknown station and dist1 to dist3 are the distances between x-point and the three stations.
Proceedings 07 00004 g001
Figure 2. Structure of values upon which dependence is modelled by Survival Clayton or Frank Copulas.
Figure 2. Structure of values upon which dependence is modelled by Survival Clayton or Frank Copulas.
Proceedings 07 00004 g002
Figure 3. (a) Taylor diagram of extreme precipitation at Sidironero. The blue circle presents the reanalysis data and red circle the bias corrected data. (b) ROC - Relative Operating Characteristics curves of extreme precipitation at Sidironero. The black line concerns the observed extreme values and red line the bias corrected values.
Figure 3. (a) Taylor diagram of extreme precipitation at Sidironero. The blue circle presents the reanalysis data and red circle the bias corrected data. (b) ROC - Relative Operating Characteristics curves of extreme precipitation at Sidironero. The black line concerns the observed extreme values and red line the bias corrected values.
Proceedings 07 00004 g003
Table 1. Extreme precipitation indices (mm) from observed, reanalysis and bias-corrected data.
Table 1. Extreme precipitation indices (mm) from observed, reanalysis and bias-corrected data.
90%95%99%
Observations32.54867
Reanalysis26.331.537.4
Bias Corrected36.342.549

Share and Cite

MDPI and ACS Style

Lazoglou, G.; Anagnostopoulou, C.; Skoulikaris, C.; Tolika, K. Copula Bias Correction for Extreme Precipitation in Reanalysis Data over a Greek Catchment. Proceedings 2019, 7, 4. https://doi.org/10.3390/ECWS-3-05817

AMA Style

Lazoglou G, Anagnostopoulou C, Skoulikaris C, Tolika K. Copula Bias Correction for Extreme Precipitation in Reanalysis Data over a Greek Catchment. Proceedings. 2019; 7(1):4. https://doi.org/10.3390/ECWS-3-05817

Chicago/Turabian Style

Lazoglou, Georgia, Christina Anagnostopoulou, Charalampos Skoulikaris, and Konstantia Tolika. 2019. "Copula Bias Correction for Extreme Precipitation in Reanalysis Data over a Greek Catchment" Proceedings 7, no. 1: 4. https://doi.org/10.3390/ECWS-3-05817

APA Style

Lazoglou, G., Anagnostopoulou, C., Skoulikaris, C., & Tolika, K. (2019). Copula Bias Correction for Extreme Precipitation in Reanalysis Data over a Greek Catchment. Proceedings, 7(1), 4. https://doi.org/10.3390/ECWS-3-05817

Article Metrics

Back to TopTop