Skip to main content

Advertisement

Log in

Assessing the risk of disruption of wind turbine operations in Saudi Arabia using Bayesian spatial extremes

  • Published:
Extremes Aims and scope Submit manuscript

Abstract

Saudi Arabia has been seeking to reduce its dependence on oil by diversifying its energy portfolio, including the largely underused energy potential from wind. However, extreme winds can possibly disrupt the wind turbine operations, thus preventing the stable and continuous production of wind energy. In this study, we assess the risk of disruptions of wind turbine operations, based on return levels with a hierarchical spatial extreme modeling approach for wind speeds in Saudi Arabia. Using a unique Weather Research and Forecasting dataset, we provide the first high-resolution risk assessment of wind extremes under spatial non-stationarity over the country. We account for the spatial dependence with a multivariate intrinsic autoregressive prior at the latent Gaussian process level. The computational efficiency is greatly improved by parallel computing on subregions from spatial clustering, and the maps are smoothed by fitting the model to cluster neighbors. Under the Bayesian hierarchical framework, we measure the uncertainty of return levels from the posterior Markov chain Monto Carlo samples, and produce probability maps of return levels exceeding the cut-out wind speed of wind turbines within their lifetime. The probability maps show that locations in the South of Saudi Arabia and near the Red Sea and the Persian Gulf are at very high risk of disruption of wind turbine operations.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Similar content being viewed by others

References

  • Banerjee, S, Carlin, B P, Gelfand, A E: Hierarchical Modeling and Analysis for Spatial Data, 2nd ed. Monographs on Statistics and Applied Probability. Chapman and Hall/CRC, Boca Raton (2014)

    Google Scholar 

  • Besag, J, Kooperberg, C: On conditional and intrinsic autoregressions. Biometrika 82(4), 733–746 (1995)

    MathSciNet  MATH  Google Scholar 

  • Besag, J, Green, P, Higdon, D, Mengersen, K: Bayesian computation and stochastic systems. Stat. Sci. 10(1), 3–41 (1995)

    MathSciNet  MATH  Google Scholar 

  • Chen, W, Castruccio, S, Genton, M G, Crippa, P: Current and future estimates of wind energy potential over Saudi Arabia. J. Geophys. Res.: Atmos. 123(12), 6443–6459 (2018)

    Article  Google Scholar 

  • Coles, S: An Introduction to Statistical Modeling of Extreme Values, vol. 208. Springer, Berlin (2001)

    Book  MATH  Google Scholar 

  • Cooley, D, Sain, S R: Spatial hierarchical modeling of precipitation extremes from a regional climate model. J. Agric. Biol. Environ. Stat. 15(3), 381–402 (2010)

    Article  MathSciNet  MATH  Google Scholar 

  • Cooley, D, Nychka, D, Naveau, P: Bayesian spatial modeling of extreme precipitation return levels. J. Am. Stat. Assoc. 102(479), 824–840 (2007)

    Article  MathSciNet  MATH  Google Scholar 

  • Davison, A C, Smith, R L: Models for exceedances over high thresholds. J. R. Stat. Soc. Ser. B (Methodol.) 52(3), 393–442 (1990)

    MathSciNet  MATH  Google Scholar 

  • Davison, A C, Padoan, S A, Ribatet, M: Statistical modeling of spatial extremes. Stat. Sci. 27(2), 161–186 (2012)

    MathSciNet  MATH  Google Scholar 

  • de Haan, L, Ferreira, A: Extreme Value Theory: an Introduction. Springer, New York (2006)

    Book  MATH  Google Scholar 

  • de Winter, R, Sterl, A, Ruessink, B: Wind extremes in the North Sea Basin under climate change: An ensemble study of 12 CMIP5 GCMs. J. Geophys. Res.: Atmos. 118(4), 1601–1612 (2013)

    Article  Google Scholar 

  • Dumat Al Jandal wind farm: https://dumataljandal.com7, Accessed: 20 May 2019 (2017)

  • DuMouchel, W H: Estimating the stable index α in order to measure tail thickness: a critique. Ann. Stat. 11(4), 1019–1031 (1983)

    Article  MathSciNet  MATH  Google Scholar 

  • Emeis, S: How well does a power law fit to a diabatic boundary-layer wind profile? DEWI Mag. 26, 59–62 (2005)

    Google Scholar 

  • Fisher, R A, Tippett, L H C: Limiting forms of the frequency distribution of the largest or smallest member of a sample. In: Mathematical Proceedings of the Cambridge Philosophical Society, vol. 24, pp 180–190. Cambridge University Press (1928)

  • Gaetan, C, Grigoletto, M: A hierarchical model for the analysis of spatial rainfall extremes. J. Agric. Biol. Environ. Stat. 12(4), 434–449 (2007)

    Article  MathSciNet  MATH  Google Scholar 

  • Gelfand, A E, Vounatsou, P: Proper multivariate conditional autoregressive models for spatial data analysis. Biostatistics 4(1), 11–15 (2003)

    Article  MATH  Google Scholar 

  • Giani, P, Tagle, F, Genton, M G, Castruccio, S, Crippa, P: Closing the gap between wind energy targets and implementation for emerging countries. Appl. Energy 269, 115085 (2020)

    Article  Google Scholar 

  • Gnedenko, B: Sur la distribution limite du terme maximum d’une série aléatoire. Ann. Math. 44, 423–453 (1943)

    Article  MathSciNet  MATH  Google Scholar 

  • Heffernan, J E, Tawn, J A: A conditional approach for multivariate extreme values (with discussion). J. R. Stat. Soc.: Ser. B (Stat. Methodol.) 66 (3), 497–546 (2004)

    Article  MathSciNet  MATH  Google Scholar 

  • Hosking, J, Wallis, J R: Regional Frequency Analysis: an Approach Based on L-Moments. Cambridge University Press, Cambridge (1997)

    Book  Google Scholar 

  • Huser, R, Davison, A: Space–time modelling of extreme events. J. R. Stat. Soc.: Ser. B (Stat. Methodol.) 76(2), 439–461 (2014)

    Article  MathSciNet  Google Scholar 

  • Jalbert, J, Favre, A C, Bélisle, C, Angers, J F: A spatiotemporal model for extreme precipitation simulated by a climate model, with an application to assessing changes in return levels over North America. J. R. Stat. Soc.: Ser. C (Appl. Stat.) 66(5), 941–962 (2017)

    Article  MathSciNet  Google Scholar 

  • KA-CARE: King Abdullah City for Atomic and Renewable Energy: Building the Renewable Energy Sector in Saudi Arabia. Report (2012)

  • Kavanagh, L, Lee, D, Pryce, G: Is poverty decentralizing? Quantifying uncertainty in the decentralization of urban poverty. Ann. Am. Assoc. Geogr. 106(6), 1286–1298 (2016)

    Google Scholar 

  • Kharin, V V, Zwiers, F W: Changes in the extremes in an ensemble of transient climate simulations with a coupled atmosphere–ocean GCM. J. Clim. 13(21), 3760–3788 (2000)

    Article  Google Scholar 

  • Kriegel, HP, Pfeifle, M: Density-based clustering of uncertain data. In: Proceedings of the Eleventh ACM SIGKDD International Conference on Knowledge Discovery in Data Mining, pp 672–677. ACM (2005a)

  • Kriegel, HP, Pfeifle, M: Hierarchical density-based clustering of uncertain data. IEEE (2005b)

  • Kumar, D, Mishra, V, Ganguly, A R: Evaluating wind extremes in CMIP5 climate models. Clim. Dyn. 45(1-2), 441–453 (2015)

    Article  Google Scholar 

  • Ledford, A W, Tawn, J A: Statistics for near independence in multivariate extreme values. Biometrika 83(1), 169–187 (1996)

    Article  MathSciNet  MATH  Google Scholar 

  • Ledford, A W, Tawn, J A: Modelling dependence within joint tail regions. J. R. Stat. Soc.: Ser. B (Stat. Methodol.) 59(2), 475–499 (1997)

    Article  MathSciNet  MATH  Google Scholar 

  • Lee, D: CARBayes: an R package for Bayesian spatial modeling with conditional autoregressive priors, vol. 55 (2013)

  • Leroux, B.G., Lei, X., Breslow, N.: Estimation of disease rates in small areas: a new mixed model for spatial dependence. In: Statistical Models in Epidemiology, the Environment, and Clinical Trials, pp 179–191. Springer, Berlin (2000)

  • Mann, E M, Steinman, B A, Miller, S K: Absence of internal multidecadal and interdecadal oscillations in climate model simulations. Nat. Commun. 11(49), 1–9 (2020)

    Google Scholar 

  • Martins, E S, Stedinger, J R: Generalized maximum-likelihood generalized extreme-value quantile estimators for hydrologic data. Water Resour. Res. 36(3), 737–744 (2000)

    Article  Google Scholar 

  • NEOM Project: https://www.neom.com/, Accessed: 28 Aug 2018 (2017)

  • Nikulin, G, Kjellströ, M E, Hansson, U, Strandberg, G, Ullerstig, A: Evaluation and future projections of temperature, precipitation and wind extremes over Europe in an ensemble of regional climate simulations. Tellus A: Dyn. Meteorol. Oceanogr. 63(1), 41–55 (2011)

    Article  Google Scholar 

  • Pickands, J: Statistical inference using extreme order statistics. Ann. Stat. 3(1), 119–131 (1975)

    MathSciNet  MATH  Google Scholar 

  • Pryor, S, Barthelmie, R: Assessing climate change impacts on the near-term stability of the wind energy resource over the United States. Proc. Natl. Acad. Sci. 108(20), 8167–8171 (2011)

    Article  Google Scholar 

  • Ramos, A, Ledford, A: A new class of models for bivariate joint tails. J. R. Stat. Soc.: Ser. B (Stat. Methodol.) 71(1), 219–241 (2009)

    Article  MathSciNet  MATH  Google Scholar 

  • Rehman, S, Ahmad, A: Assessment of wind energy potential for coastal locations of the Kingdom of Saudi Arabia. Energy 29(8), 1105–1115 (2004)

    Article  Google Scholar 

  • Rehman, S, El-Amin, I, Ahmad, F, Shaahid, S, Al-Shehri, A, Bakhashwain, J: Wind power resource assessment for Rafha, Saudi Arabia. Renew. Sustain. Energy Rev. 11(5), 937–950 (2007)

    Article  Google Scholar 

  • Reich, B J, Shaby, B A: A spatial Markov model for climate extremes. J. Comput. Graphic. Stat. 28(1), 117–126 (2019)

    Article  MathSciNet  Google Scholar 

  • Resnick, S I: Extreme Values, Regular Variation and Point Processes. Springer, New York (1987)

    Book  MATH  Google Scholar 

  • Rue, H, Held, L: Gaussian Markov Random Fields: Theory and Applications. Chapman and Hall/CRC, Boca Raton (2005)

    Book  MATH  Google Scholar 

  • Sang, H, Gelfand, A E: Hierarchical modeling for extreme values observed over space and time. Environ. Ecol. Stat. 16(3), 407–426 (2009)

    Article  MathSciNet  Google Scholar 

  • Sang, H, Gelfand, A E: Continuous spatial process models for spatial extreme values. J. Agric. Biol. Environ. Stat. 15(1), 49–65 (2010)

    Article  MathSciNet  MATH  Google Scholar 

  • Scarrott, C, MacDonald, A: A review of extreme value threshold estimation and uncertainty quantification. REVSTAT–Stat. J. 10(1), 33–60 (2012)

    MathSciNet  MATH  Google Scholar 

  • Shaahid, S, Al-Hadhrami, L M, Rahman, M: Potential of establishment of wind farms in western province of Saudi Arabia. Energy Procedia 52, 497–505 (2014)

    Article  Google Scholar 

  • Smith, R L: Maximum likelihood estimation in a class of nonregular cases. Biometrika 72(1), 67–90 (1985)

    Article  MathSciNet  MATH  Google Scholar 

  • Smith, R L: Extreme value analysis of environmental time series: an application to trend detection in ground-level ozone. Stat. Sci. 4(4), 367–377 (1989)

    MathSciNet  MATH  Google Scholar 

  • Tagle, F, Castruccio, S, Crippa, P, Genton, M G: A non-Gaussian spatio-temporal model for daily wind speeds based on a multivariate skew-t distribution. J. Time Ser. Anal. 40, 312–326 (2019)

    Article  MATH  Google Scholar 

  • Telesca, L, Guignard, F, Laib, M, Kanevski, M: Analysis of temporal properties of wind extremes. arXiv:180808847 (2018)

  • Tobler, W. R.: A computer movie simulating urban growth in the Detroit region. Econ. Geogr. 46(sup1), 234–240 (1970)

    Article  Google Scholar 

  • Towe, R P, Tawn, J A, Lamb, R, Sherlock, C G: Model-based inference of conditional extreme value distributions with hydrological applications. Environmetrics 30, e2575 (2019)

    Article  MathSciNet  Google Scholar 

  • Turkman, K F, Turkman, M A, Pereira, J: Asymptotic models and inference for extremes of spatio-temporal data. Extremes 13(4), 375–397 (2010)

    Article  MathSciNet  MATH  Google Scholar 

  • Vision 2030: http://vision2030.gov.sa/en/node/87, Accessed: 28 Aug 2018 (2016)

  • Weiss, J, Bernardara, P, Benoit, M: Modeling intersite dependence for regional frequency analysis of extreme marine events. Water Resour. Res. 50(7), 5926–5940 (2014)

    Article  Google Scholar 

  • Winter, H C, Tawn, J A, Brown, S J, et al.: Modelling the effect of the El Nino-Southern Oscillation on extreme spatial temperature events over Australia. Ann. Appl. Stat. 10(4), 2075–2101 (2016)

    Article  MathSciNet  MATH  Google Scholar 

  • Xie, L, Gu, Y, Zhu, X, Genton, M G: Short-term spatio-temporal wind power forecast in robust look-ahead power system dispatch. IEEE Trans. Smart Grid 5(1), 511–520 (2014)

    Article  Google Scholar 

  • Yip, C M A: Statistical Characteristics and Mapping of Near-Surface and Elevated Wind Resources in the Middle East. PhD thesis, King Abdullah University of Science and Technology (2018)

  • Yip, C M A, Gunturu, U B, Stenchikov, G L: Wind resource characterization in the Arabian Peninsula. Appl. Energy 164, 826–836 (2016)

    Article  Google Scholar 

  • Zhu, X, Genton, M G, Gu, Y, Xie, L: Space-time wind speed forecasting for improved power system dispatch. TEST 23(1), 1–25 (2014)

    Article  MathSciNet  MATH  Google Scholar 

Download references

Acknowledgements

This publication is based on research supported by the King Abdullah University of Science and Technology (KAUST) Office of Sponsored Research (OSR) under Award No: OSR-2018-CRG7-3742 and in part by the Center of Excellence for NEOM Research at KAUST. We are grateful to Professor Georgiy Stenchikov’s group, the Atmospheric and Climate Modeling group at KAUST, for producing and providing the high-resolution WRF dataset. Many thanks also to Professor Daniel Cooley for providing the codes used in the article by Cooley and Sain (2010).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Wanfang Chen.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Appendices

Appendix A: Supplementary figures, tables and discussions

Fig. 8
figure 8

Cluster means for GMLEs of the shape parameter ξ with different values of k and different assigned weights for longitude, latitude and GMLEs for ξ, respectively

Table 2 Range of cluster sizes with different values of k and different weights. The two “best” configurations are in bold
Fig. 9
figure 9

Example for a cluster neighbor. Each color represents one single cluster, and all the seven clusters represent the cluster neighbor of the central cluster in blue

Fig. 10
figure 10

QQ-plots of marginal GPD fitting for 9 randomly selected locations, using 95% quantile as the threshold

1.1 Discussion on CAR and IAR specifications for spatial random effects

There have been discussions on whether the proper or improper CAR (referred to as IAR) specification should be used in practice (e.g., Besag et al. 1995, Besag and Kooperberg 1995, Gelfand and Vounatsou 2003 and Banerjee et al. 2014). As Besag et al. (1995) showed, the marginal maximal bivariate correlation that can be captured with a proper Gaussian field is around 0.6. Besag and Kooperberg (1995) pointed out a common disadvantage of a proper CAR that appreciable correlations between the spatial random effects at neighboring sites require parameter values extremely close to a particular boundary of the parameter space. Gelfand and Vounatsou (2003) demonstrated that IAR is analogous to the nonstationary or random walk case in familiar autoregressive time series models and can be advantageous in accommodating more irregular spatial behaviors. Banerjee et al. (2014) also claimed that the breadth of spatial patterns may be too limited if the proper CAR is used, and the improper IAR choice may actually enable a wider scope for posterior spatial patterns.

Fig. 11
figure 11

Values of threshold and KL divergence for marginal extremes using 90% (top) and 95% (bottom) quantile, respectively, as a rule of thumb in the GPD model

Here we first implement our hierarchical model with the proper CAR prior used for the spatial random effects ϕ. The multivariate CAR model (Kavanagh et al. 2016) we used is given by:

$$ \boldsymbol{\phi} \sim \mathcal{N}_{3N}\left( \boldsymbol{0},[\boldsymbol{Q}(\boldsymbol{W},\rho)\otimes \boldsymbol{\Sigma}^{-1}]^{-1} \right), $$

where Q(W,ρ) = ρ[diag(W1) −W] + (1 − ρ)I (1 is the N × 1 vector of ones, I is the N × N identity matrix) is the N × N precision matrix for the joint distribution corresponding to the CAR prior proposed by Leroux et al. (2000), while Σ is a 3 × 3 cross-variables covariance matrix. The matrix Q(W,ρ) controls the spatial autocorrelation structure of the random effects, and is based on a non-negative symmetric N × N neighborhood (or adjacency) matrix W, and a spatial dependence parameter ρ. We use the common binary specification for W, where its entry wij = 1 if the grid cells i and j are adjacent, and is zero otherwise. The parameter ρ is a spatial autoregressive parameter, with ρ close to one corresponding to strong spatial dependence and ρ = 0 corresponding to independence in space. When ρ = 1, we obtain the multivariate IAR model. The correlation structure is specified via the full conditionals:

$$ (\boldsymbol{\phi}_{i}|\boldsymbol{\phi}_{-i},\boldsymbol{W},\boldsymbol{\Sigma},\rho) \sim \mathcal{N}_{3}\left( \frac{\rho{\sum}_{k=1}^{N}w_{ki}\boldsymbol{\phi}_{i}}{\rho{\sum}_{k=1}^{N}w_{ki}+1-\rho}, \frac{\boldsymbol{\Sigma}}{\rho{\sum}_{k=1}^{N}w_{ki}+1-\rho}\right), $$

where ϕi denotes the set of spatial random effects except those at the i th location. With the choice of the matrix W, the conditional expectation of spatial random effect at one location is a weighted average of the random effects in its adjacent locations, and the covariance is weighted by the number of adjacent locations.

In our Bayesian hierarchical model, a Uniform[0, 1] prior is assigned to the spatial autoregressive parameter ρ, as the negative spatial autocorrelation is rarely seen in practice in spatial areal unit data (Tobler 1970), and ρ ∈ [0, 1) is a sufficient condition for the covariance matrix of the joint distribution to be nonsingular (Banerjee et al. 2014). The parameter ρ is updated with the MH algorithm, where the candidate for ρ is drawn from a truncated normal distribution in the unit interval so as to bound ρ in [0, 1). Other settings for priors and computational details are the same as in the main text. The posterior density from MCMC samples for ρ is peaked near ρ = 1 for all subregions (see Figs. 12 and 13), suggesting that there is more spatial dependence in the data than the model can capture. Therefore, we replace the multivariate CAR with the multivariate IAR in order to capture the irregular and strong spatial dependence in our high-resolution data. Although the IAR is improper, we are only using it as a prior; the posterior will typically still emerge as proper, so Bayesian inference can still proceed. On the other hand, we can impose a sum-to-zero constraint on ϕ as a remedy of impropriety, which is numerically convenient in the MCMC sampling procedure.

Fig. 12
figure 12

Trace plot of 10,000 MCMC iterations for the parameter ρ in the MCAR model in a random spatial cluster as selected in Fig. 1 in the main text

Fig. 13
figure 13

Posterior mean of the spatial autoregressive parameter ρ in the MCAR model from MCMC samples with Bayesian hierarchical spatial extremes modeling

Fig. 14
figure 14

Potential best locations for siting wind farms over Saudi Arabia where the wind speeds exceed 9 m/s for at least half of the time in Summer, and the risk of disruption of wind turbine operations is lower than 1% for the setting of a k = 200, Threshold = 95% quantile and b k = 250, Threshold = 95% quantile. “prob” stands for “probability of the 30-year return levels that exceed 25 m/s”

Fig. 15
figure 15

Cluster means of shape GMLEs with a all data and b locations with highly uncertain marginal shape estimates (i.e., where the length of temporal cluster maxima, \(N_{u_{i}}\), is less than 50) removed

Appendix B: Supplementary R codes

The datasets generated and/or analyzed during the current study were used under license, and so are not publicly available. Data are available however from the corresponding author upon reasonable request. The R codes related to this article can be found online at the github repository: https://github.com/wanruofenfang123/Bayesian-Hierarchical-Spatial-Extremes.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, W., Castruccio, S. & Genton, M.G. Assessing the risk of disruption of wind turbine operations in Saudi Arabia using Bayesian spatial extremes. Extremes 24, 267–292 (2021). https://doi.org/10.1007/s10687-020-00384-1

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s10687-020-00384-1

Keywords

AMS 2000 Subject Classifications

Navigation