- uncertainty and spatio-temporal trends Juha Aalto, Pentti Pirinen, - - PowerPoint PPT Presentation

uncertainty and spatio temporal trends
SMART_READER_LITE
LIVE PREVIEW

- uncertainty and spatio-temporal trends Juha Aalto, Pentti Pirinen, - - PowerPoint PPT Presentation

Producing a long-term gridded data set in Finland - uncertainty and spatio-temporal trends Juha Aalto, Pentti Pirinen, Kirsti Jylh 10th EUMETNET Data Management Workshop, St. Gallen, Switzerland, 30.10.2015 1. Introduction and aims PLUMES


slide-1
SLIDE 1

Producing a long-term gridded data set in Finland

  • uncertainty and spatio-temporal trends

Juha Aalto, Pentti Pirinen, Kirsti Jylhä 10th EUMETNET Data Management Workshop,

  • St. Gallen, Switzerland, 30.10.2015
slide-2
SLIDE 2
  • 1. Introduction and aims
  • PLUMES consortium
  • Task: create a high-quality daily gridded climate data set of

the key variables across 1961-2010 (”FMI_ClimGrid_1.0”)

  • Focus on interpolation uncertainty
  • Use gridded data to investigate temporal trends in climate
  • Compare the results with existing data (E-OBS)

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 2

slide-3
SLIDE 3
  • 1. Introduction – gridded data
  • Spatially continous data based on a set of observations
  • Most often based on a statistical model
  • Important applications: climate change studies, forest

management, agriculture, biosphere modelling, permafrost ….

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 3

Observations Spatial grid Gridded data

slide-4
SLIDE 4
  • 2. Data – observations
  • Seven climate variables:
  • mean temperature (Tday)
  • maximum temperature (Tmax)
  • minimum temperature (Tmin)
  • precipitation sum (Prec)
  • mean relative humidity (RH)
  • air pressure (P)
  • snow depth (Sn)

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 4

  • Data sources:
  • FMI database
  • ECA&D pan-European database
  • Sweden, Norway, Russian and

Estonia

slide-5
SLIDE 5
  • 2. Data – observations

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 5

slide-6
SLIDE 6
  • 2. Data – quality control
  • National operational QC
  • ”Non-blended” ECA&D series
  • Misscodings, duplicates
  • Local outlier detection protocol:
  • 1. compare each value to local average

and stdev (station in turn excluded)

  • 2. Compare the local stdev to long-term

monthly stdev (1961-2010)

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 6

150km

slide-7
SLIDE 7
  • 2. Data – grid specifications
  • Spatial resolution = 10 km x 10 km
  • Euref-FIN TM35 (epsg: 3067)
  • 5224 points (3364 inside, 1860
  • utside Finland)

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 7

slide-8
SLIDE 8
  • 2. Methods – kriging interpolation
  • Kriging interpolates the value at given point using a weighted

average of the know values inside a neighborhood

  • Weights are assigned by (decreasing) function of the distance,

based on the spatial covariance structure

  • Variogram is used to quantify the spatial dependency in the data

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 8

slide-9
SLIDE 9
  • 2. Data – background data

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 9

  • Used as covariates in the

interpolation model (i.e. trend model)

  • Latitude and longitude
slide-10
SLIDE 10
  • 2. Methods – details
  • Separate trend model was estimated for each day
  • ”Semi” climatological variogram models:
  • Range = monthly means of daily ranges (1961-2010)
  • Separate sill for each day
  • Nugget = 30 % of the measurement precision (e.g. 0.03 for Temp)
  • Exponential variogrammodels
  • Global kriging

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 10

slide-11
SLIDE 11
  • 2. Methods – interpolating

precipitation and snow

  • High and potentially discrete variation vs. sparse
  • bservation network
  • Satellite and radar data might improve
  • Solution: interpolate the probability of precipitation / snow

depth and combine with interpolated amounts

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 11

slide-12
SLIDE 12
  • 2. Methods - evaluation
  • 20 independent evaluation

stations

  • Compare the observed and

interpolated values

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 12

slide-13
SLIDE 13
  • 3. Results – interpolation accuracy

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 13

slide-14
SLIDE 14
  • 3. Results – seasonal variation in

accuracy

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 14

slide-15
SLIDE 15
  • 3. Results – interpolation accuracy

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 15

  • Some meteorological conditions are more challenging

to interpolate than others…

Temperature inversion?

slide-16
SLIDE 16
  • 3. Results - uncertainty

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 16

  • 50 random

permutation / day

  • > 50 different

interpolations

  • Daily uncertainty

estimate for each variable Sources of uncertainty:

  • Observations

(measurements, network, inhomogeneities …)

  • Background variables

(georeferencing, averaging…)

  • Interpolation method
slide-17
SLIDE 17
  • 3. Results – a comparison with E-OBS

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 17

slide-18
SLIDE 18
  • 4. Trends in past climate

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 18

  • Daily grid averages
  • > seasonal / annual aggregates
  • Most uncertain areas excluded

from the trend analysis

slide-19
SLIDE 19
  • 4. Trends in past climate

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 19

slide-20
SLIDE 20
  • 5. Conclusions
  • Long-term gridded dataset were succesfully produced
  • Daily permutation-based uncertainty estimates
  • Clear, but locally varying signal of past climate change
  • Wind and solar radiation in the future
  • The dataset will be made freely available with regular

updates

  • Manuscript in progress…

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 20

slide-21
SLIDE 21

Computing environment

  • R in linux server (FMI supercomputer ”Voima”)
  • Required R-packages: gstat, sp, rgdal, raster, maptools,

PresenceAbsence, Roracle

  • Total time of calculations ~ 6 days / per variable

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 21

slide-22
SLIDE 22

www.fmi.fi

10/30/2015 Juha Aalto | juha.aalto@fmi.fi 22

More information:

juha.aalto@fmi.fi pentti.pirinen@fmi.fi