On the Smoothing of the Generalized Extreme Value Distribution Parameters Using Penalized Maximum Likelihood: A Case Study on UVB Radiation Maxima in the Mexico City Metropolitan Area

This paper concerns the use and implementation of penalized maximum likelihood procedures to fitting smoothing functions of the generalized extreme value distribution parameters to analyze spatial extreme values of ultraviolet B radiation (UVB) radiation across the Mexico City metropolitan area in the period 2000&#8722;2018. The model was fitted using a flexible semi-parametric approach and the parameters were estimated by the penalized maximum likelihood (PML) method. In order to investigate the performance of the model as well as the estimation method in the analysis of complex nonlinear trends for UVB radiation maxima, a simulation study was conducted. The results of the simulation study showed that penalized maximum likelihood yields better regularization to the model than the maximum likelihood estimates. We estimated return levels of extreme UVB radiation events through a nonstationary extreme value model using measurements of ozone (<inline-formula> <math display="inline"> <semantics> <msub> <mi mathvariant="normal">O</mi> <mn>3</mn> </msub> </semantics> </math> </inline-formula>), nitrogen oxides (<inline-formula> <math display="inline"> <semantics> <msub> <mi>NO</mi> <mi>x</mi> </msub> </semantics> </math> </inline-formula>), particles of 10 <inline-formula> <math display="inline"> <semantics> <mi mathvariant="sans-serif">&#956;</mi> </semantics> </math> </inline-formula>m or less in diameter (<inline-formula> <math display="inline"> <semantics> <msub> <mi>PM</mi> <mn>10</mn> </msub> </semantics> </math> </inline-formula>), carbon monoxide (CO), relative humidity (RH) and sulfur dioxide (<inline-formula> <math display="inline"> <semantics> <msub> <mi>SO</mi> <mn>2</mn> </msub> </semantics> </math> </inline-formula>). The deviance statistics indicated that the nonstationary generalized extreme value (GEV) model adjusted was statistically better compared to the stationary model. The estimated smoothing functions of the location parameter of the GEV distribution on the spatial plane for different periods of time reveal the existence of well-defined trends in the maxima. In the temporal plane, a presence of temporal cyclic components oscillating over a weak linear component with a negative slope is noticed, while in the spatial plane, a weak nonlinear local trend is present on a plane with a positive slope towards the west, covering the entire study area. An explicit spatial estimate of the 25-year return period revealed that the more extreme risk levels are located in the western region of the study area.

Abstract licence: CC BY

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.