![Page 1: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/1.jpg)
www.nr.no
Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland
Orskaug E.a, Scheel I.b, Frigessi A.c,a, Guttorp P.d,a, Haugen J. E.e, Tveito O. E.e, Haug O.a
a Norwegian Computing Center, Oslo, Norway
b Department of Mathematics, University of Oslo, Oslo, Norway
c Department of Biostatistics, University of Oslo, Oslo, Norway
d University of Washington, Seattle, USA
e The Norwegian Meteorological Institute, Oslo, Norway
![Page 2: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/2.jpg)
Motivation
► Climate research produces an increasing number of data sets combining different GCMs, CO2 emission scenarios and downscaling techniques.
► For impact studies, but also as an issue of separate interest, the quality of these data need to be verified.
![Page 3: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/3.jpg)
Goal
► We want to compare downscaled ERA-40 reanalysis data (RCM) against observations of Norwegian precipitation.▪ How good are the RCM data?▪ Where (in the distribution) does the RCM differ from
the observations?▪ Where (geographically) does the RCM perform
best/worst?
![Page 4: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/4.jpg)
Why is this work important?
► It assesses the quality of a dynamic downscaled data and highlights which areas these data capture reality and where there are deviations from the truth.
► Another aim is to show how standard methods of statistical testing may be used to assess dynamic downscaling.
![Page 5: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/5.jpg)
Data
Model data► RCM model, dynamically
downscaled HIRHAM model, forced by ERA-40 reanalysis data from the ENSEMBLES project.
► Spatial resolution of 25 x 25 km2.
► Reliant on the downscaling, still supposed to possess properties similar to real weather locally over longer time periods.
Observations
► Interpolations (1 x 1 km2) from a triangulation of the official measurement stations operated by The Norwegian Meteorological Institute.
► Aggregated to 25 x 25 km2 scale by collecting 1 x 1 km2 grid cells with centre points within the RCM cell, the mean is representing the precipitation within that grid cell.
![Page 6: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/6.jpg)
Data – The RCM
► The RCM from the ENSEMBLES project
![Page 7: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/7.jpg)
Data – properties for both data sets
► Climate variable: precipitation
► Time period: 1961 – 2000
► Time scale: Daily, seasonal
► Resolution: 25 x 25 km2
► Number of grid cells: 777 grid cells covering Norway
![Page 8: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/8.jpg)
Methods for comparison
► Evaluate the distributions1. Global measure:
Kolmogorov Smirnov test
2. Local measures:
![Page 9: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/9.jpg)
Comments
► Drizzle effect avoided: conditioned on wet days; i.e. days with precipitation below a small, positive threshold (0.5 mm/day) are discarded.
► Day-to-day correlation in the RCM is partly lost due to downscaling, hence the distributions have to be compared instead of comparing day by day.
► Separate tests for each grid cell and each season.
![Page 10: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/10.jpg)
Kolmogorov-Smirnov test
► K-S two sample test is used to check whether the empirical distributions from the RCM and the observations are equal.
► To avoid the problem of tied data, a small, random normally distributed number, N(0, σ2), is added to each data point.
σ = 1e-7
![Page 11: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/11.jpg)
Kolmogorov-Smirnov test – Results
► The null hypothesis of equality of the distributions are rejected for almost all grid cells for all the four seasons.
► Global picture: the RCM does not have the same distribution as the observations.
► Next: want to find out where the distributions differ; local measures.
![Page 12: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/12.jpg)
Methods for comparison
► Evaluate the distributions1. Global measure:
Kolmogorov Smirnov test
2. Local measures:
![Page 13: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/13.jpg)
Test equality of quantiles
Construction of the 2 x 2 contigency table
![Page 14: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/14.jpg)
0.05-quantile – Results
► Hardly any rejections of null hypothesis of equality.
► For low quantiles: the RCM reproduces the observations well both season- and nationwide.
![Page 15: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/15.jpg)
0.95-quantile – Results
► Mainly rejections of the null hypothesis of equality.
► Overall picture: the RCM underestimates high precipitation.
![Page 16: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/16.jpg)
Generalized Pareto Distribution (GPD)
►
![Page 17: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/17.jpg)
GPD – Results
► One-year return levels from GPD are more similar than expressed through the Kolmogorov-Smirnov test.
► But still: tendency that the RCM underestimates high precipitation.
![Page 18: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/18.jpg)
Wet day frequency
► Wet day frequency = Proportion of wet days (among all days in the data)
► A wet day is defined to be above 0.5 mm/day for both data sets.
► The equality of the wet day frequency is tested by permutation testing.
![Page 19: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/19.jpg)
Wet day frequency – Results
► Mainly rejections of the null hypothesis of equality.
► Total picture: Wet day frequency of the RCM is greater than for the observations.
![Page 20: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/20.jpg)
Summary
► Small amounts of rainfall: the RCM shows good agreement with the observations.
► When rainfall amounts is beyond the first quartile, the agreement disappear.
► The RCM has too many and too small rain events for all seasons.
► This work is accepted for publication in Tellus A.
An improvement/correction of the RCM is needed.
![Page 21: Www.nr.no Evaluation of a dynamic downscaling of precipitation over the Norwegian mainland Orskaug E. a, Scheel I. b, Frigessi A. c,a, Guttorp P. d,a,](https://reader035.vdocument.in/reader035/viewer/2022062722/56649f285503460f94c406e6/html5/thumbnails/21.jpg)
What to do next?
► We want to add a statistical correction method to the output of the RCM, especially improve the right tail.
► Simple linear regression was tried out, but did not improve the results.
► We are currently working on a more complex transformation with spatial corrections.