A non-stationary extreme value approach for climate projection ensembles: application to snow loads in the French Alps
- 1Univ. Grenoble Alpes, INRAE, UR ETNA, Grenoble, France
- 2Univ. Grenoble Alpes, Grenoble INP, CNRS, IRD, IGE, Grenoble, France
- 3Univ. Grenoble Alpes, Univ. Toulouse, Météo France, CNRS, CNRM, CEN, Grenoble, France
- 1Univ. Grenoble Alpes, INRAE, UR ETNA, Grenoble, France
- 2Univ. Grenoble Alpes, Grenoble INP, CNRS, IRD, IGE, Grenoble, France
- 3Univ. Grenoble Alpes, Univ. Toulouse, Météo France, CNRS, CNRM, CEN, Grenoble, France
Abstract. Anticipating risks related to climate extremes often relies on the quantification of large return levels (values exceeded with small probability) from climate projection ensembles. Current approaches based on multi-model ensembles (MMEs) usually estimate return levels separately for each chain of the MME. By contrast, using MME obtained with different combinations of general circulation model (GCM) and regional climate model (RCM), our approach estimates return levels together from the past observations and all GCM-RCM pairs, considering both historical and future periods. The proposed methodology seeks to provide estimates of projected return levels accounting for the variability of individual GCM-RCM trajectories, with a robust quantification of uncertainties. To this aim, we introduce a flexible non-stationary generalized extreme value (GEV) distribution that includes i) piecewise linear functions to model the changes in the three GEV parameters ii) adjustment coefficients for the location and scale parameters to adjust the GEV distributions of the GCM-RCM pairs with respect to the GEV distribution of the past observations. Our application focuses on snow load at 1500 m elevation for the 23 massifs of the French Alps, which is of major interest for the structural design of roofs. Annual maxima are available for 20 adjusted GCM-RCM pairs from the EURO-CORDEX experiment, under the scenario RCP8.5. Our results show with a model-as-truth experiment that at least two linear pieces should be considered for the piecewise linear functions. We also show, with a split-sample experiment, that eight massifs should consider adjustment coefficients. These two experiments help us select the GEV parameterizations for each massif. Finally, using these selected parameterizations, we find that the 50-year return level of snow load is projected to decrease in all massifs, by −2.9 kN m−2 (−50 %) on average between 1986–2005 and 2080–2099 at 1500 m elevation and RCP8.5. This paper extends to climate extremes the recent idea to constrain climate projection ensembles using past observations.
- Preprint
(2303 KB) -
Supplement
(796 KB) - BibTeX
- EndNote
Erwan Le Roux et al.
Status: final response (author comments only)
-
RC1: 'Comment on esd-2021-79', Anonymous Referee #1, 23 Nov 2021
The comment was uploaded in the form of a supplement: https://esd.copernicus.org/preprints/esd-2021-79/esd-2021-79-RC1-supplement.pdf
-
AC1: 'Reply on RC1', Guillaume Evin, 16 Mar 2022
The comment was uploaded in the form of a supplement: https://esd.copernicus.org/preprints/esd-2021-79/esd-2021-79-AC1-supplement.pdf
-
AC1: 'Reply on RC1', Guillaume Evin, 16 Mar 2022
-
RC2: 'Comment on esd-2021-79', Antonio Speranza, 27 Nov 2021
General comments
In this paper the authors make use of a non-stationary extreme values probability estimation approach based on using both past observations and future projections (from explicit numerical models); the proposed analysis is addressed to evaluating 50 years return levels (for the time period 2019-2100) of annual maxima of snow load at 1500 m elevation for the 23 massifs of the French Alps: past observations and model (GCM-RCM pairs) projections are simultaneously submitted to extreme value probability estimation by means of GEV distributions with piecewise linear trend (in the GEV parameters) and adjustment coefficients suited to making GEV of past observation and future projections compatible.
The paper is, in my opinion, interesting as it addresses problems that are relevant from a both conceptual and applicative point of view:
- estimating extreme distributions of statistical processes embedded in background dynamics evolving on time scales comparable with those of recurrence of extremes themselves is a well known statistical problem and distinguishing “natural variability” of extremes from variability “forced” by background dynamics is a classical scientific-technical challenge;
- projecting in future time variables (temperature, precipitation, snow load, etc.) that are relevant for the management of social-economic activities is a clearly important operational task.
The methodological approach proposed in the paper, based on simultaneous use of observations and model projections, is stimulating in particular when facing problems (like the one addressed in the paper) in which numerical simulation models are characterized by heavy tuning-parameterization (fudge factors) changing in time. However, in order to make the estimation process tractable in the specific application considered in the paper, many ad hoc assumptions have to be introduced and the results are admittedly (Section 5.2) problematic, raising doubts concerning the applicability of the adopted working hypotheses. This situation often occurs in operational statistical estimation: thorough a posteriori analysis is almost invariably required; the authors should critically re-examine their assumptions.
The paper is neat and clean, but here and there not easily legible as it is very concise (“dense”): since the paper proposes issues of potential interest for a wide audience in which “non experienced” readers could find elements of interest I suggest a more “friendly” communication approach; but I leave to the authors deciding whether being concise is more important than being readily accessible for a wider audience.
Specific comments
Many acronyms and “technical slang” words appear in the paper: a glossary may help.
Line 3 “chain of MME”: define in text or in glossary.
Line 7 “with a robust quantiï¬cation of uncertainties.”: this claim appears repeatedly in the paper; I found mathematical definition in Appendix A: Uncertainty estimation a technical quantification uncertainties, but not an analysis-discussion of the “robustness” of the estimation itself.
Line 11 “is of major interest for the structural design of roofs”: not only (skying, avalanches, mobility, etc.); a few more words about applications could help.
Line 24 “EVT makes it possible to robustly estimate return levels”: see Line 7 comment.
Line 29 “estimated separately on each chain of the MME”: see Line 3 comment.
Line 32 “30-year time slices”: perhaps it is worth mentioning that 30 years is the traditional (WMO) “time scale” of “climatological” analysis.
Line 52 “robustly quantify uncertainties”: see above lines 7 and 24 comments.
Line 63 “adjustment coefficients”: a few more words could help.
Line 80 “snow load” see Line 11 comment.
Line 92 “Quantile mapping method ADAMONT”: a few words about it?
Line 98 “Crocus”: ?
Line 204 “For a detailed analysis of the mean logarithmic scores of each parameterization for each massif, see Supplement, Part C.”: what is Supplement, Part C? Where is it?
Fig.4 This figure plays a central role in the paper: some graphical features are too faint.
Line 220 “adjustment adjustment”.
Line 255 “Figure 2.3 of IPCC (2019)”: wouldn’t it be possible to insert this figure or its direct internet link in the text?
Line 273 “because it sometimes leads to prediction failure, i.e. the predictive distribution gives a null probability to some future annual maxima.”: this is not clear to me!
Line 287 “The 90% uncertainty intervals of return levels (Fig. 4) account both for the sampling uncertainty (Appendix A) and the climate model uncertainty (distributions are ï¬tted together from the past observations and all GCM-RCM pairs).”: not easy to distinguish in the figure (see comment to Fig.4 above).
-
AC2: 'Reply on RC2', Guillaume Evin, 16 Mar 2022
The comment was uploaded in the form of a supplement: https://esd.copernicus.org/preprints/esd-2021-79/esd-2021-79-AC2-supplement.pdf
-
RC3: 'Comment on esd-2021-79', Tamas Bodai, 30 Nov 2021
The paper under review "A non-stationary extreme value approach for climate projection ensembles: application to snow loads in the French Alps" estimates 50-year return levels of snow load via fitting annual maxima historical and projection data by a non-stationary GEV distribution, with the global mean surface temperature as the co-variate.
My feeling is that this paper is not providing a good solution to a real problem. They consider several GCM-RCM model pairs ('model' in the following). In the Introduction the authors point out that previous studies evaluated extreme value statistics (EVS) for individual models and in some cases they took then the ensemble mean of return levels. In this regard the authors are concerned that the estimates for separate models are not so reliable because of data scarcity. However, they seem to do this themselves. They introduce the concept of "adjustment coefficients", which is really just a difference of an estimate of the GEV parameter for a particular model (or subset of the data) from the estimate upon lumping all the data. I think we really don't need a name for this, beside the problem that they do what they criticised. On the other hand, lumping all the data together, in order to have seemingly more robust estimates, is also problematic. As i pointed out in some recent publications of mine, a model ensemble (or multi-model ensemble, MME, as it's most often called) does not represent an objective probability distribution. As such, fitting a GEV to MME data is flawed methodology. It has no meaningful probabilistic interpretation. On the contrary, doing this for (converged) initial condition ensembles is fine.
It is actually a profound scientific challenge how to use MME data in a meaningful way. I don't mean to discourage anyone from trying, though, and hope that real progress can be made.
The authors promised a constraining of the estimates/projections using observational data. Emergent constraints is now a popular concept, but it appears to me that the authors did nothing like that. They simply threw the observational data in the mix. However, the information from it is diluted by the large amount of model data.
Obs data is rather used for bias correction. I'm not sure this was done. Or, if it was done, then it seems to have even less use to throw the obs data into the mix for doing EVS.
I share my detailed comments on the manuscript with the authors in an annotated pdf. Hopefully it is useful one way or another. I'm sorry that i cannot be more positive this time. If i misunderstand something, i'm happy to learn from the author's response.
Please note that i always do peer-review non-anonymously, and i never make recommendation for or against publication. I leave this wholly with the editor. If i submit a recommendation, it is only to bypass the rigidity of the electronic submission system of the journal and therefore please consider it void.
Tamas Bodai
-
AC3: 'Reply on RC3', Guillaume Evin, 16 Mar 2022
The comment was uploaded in the form of a supplement: https://esd.copernicus.org/preprints/esd-2021-79/esd-2021-79-AC3-supplement.pdf
-
AC3: 'Reply on RC3', Guillaume Evin, 16 Mar 2022
Erwan Le Roux et al.
Erwan Le Roux et al.
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
632 | 143 | 21 | 796 | 43 | 4 | 5 |
- HTML: 632
- PDF: 143
- XML: 21
- Total: 796
- Supplement: 43
- BibTeX: 4
- EndNote: 5
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1