Title: | Dataset for Climate Analysis with Data from the Nordic Region |
---|---|
Description: | The Nordklim dataset 1.0 is a unique and useful achievement for climate analysis. It includes observations of twelve different climate elements from more than 100 stations in the Nordic region, in time span over 100 years. The project contractors were NORDKLIM/NORDMET on behalf of the National meteorological services in Denmark (DMI), Finland (FMI), Iceland (VI), Norway (DNMI) and Sweden (SMHI). |
Authors: | Heikki Tuomenvirta [aut], Achim Drebs [aut], Eirik Førland [aut], Ole Einar Tveito [aut], Hans Alexandersson [aut], Ellen Vaarby Laursen [aut], Trausti Jónsson [aut], Jose Gama [cre] |
Maintainer: | Jose Gama <[email protected]> |
License: | GPL (>= 3) |
Version: | 1.2 |
Built: | 2024-11-27 06:47:44 UTC |
Source: | CRAN |
The NORDKLIM data set contains close to 70 000 years of monthly data from 114 stations. The station network covers all five Nordic countries, including data from the Faeroe Islands, Jan Mayen, Bjornoya and Svalbard. There are seven monthly climatic elements describing temperature, two on precipitation and one on air pressure, cloud cover and snow cover.
Project contractors: NORDKLIM/NORDMET on behalf of the National meteorological services in Denmark (DMI), Finland (FMI), Iceland (VI), Norway (DNMI) and Sweden (SMHI)
Package: | nordklimdata1 |
Type: | Package |
Version: | 1.0 |
Date: | 2013-03-10 |
License: | GPL (>= 3) |
Authors@R: c(person(given = "Heikki", family = "Tuomenvirta", role = "aut"), person(given = "Achim", family = "Drebs", role = "aut"), person(given = "Eirik", family = "Forland", role = "aut"), person(given = "Ole", family = "Einar Tveito", role = "aut"), person(given = "Hans", family = "Alexandersson", role = "aut"), person(given = "Ellen", family = "Vaarby Laursen", role = "aut"), person(given = "Trausti", family = "Jonsson", role = "aut"), person(given = "Jose'", family = "Gama", role = c("cre"),email = "[email protected]"))
http://www.smhi.se/hfa_coord/nordklim
Nordklim dataset 1.0 - description and illustrations Norwegian meteorological institute, 08/01 KLIMA, 2001
The NORDKLIM data set - monthly data for 7 climatic elements from 114 stations in 5 Nordic countries.
data(NordklimData)
data(NordklimData)
A data frame with 71329 observations on the following 16 variables.
Nordklim number identifier
Climate element identifier
First year of the dataset
Readings for January
Readings for February
Readings for March
Readings for April
Readings for May
Readings for June
Readings for July
Readings for August
Readings for September
Readings for October
Readings for November
Readings for December
Country code
The NORDKLIM data set has 16 columns, the first three columns are the Nordklim number, climate element number and first year of the dataset, the next 12 columns are twelve months of readings, from January to December and the last column is the country code. Monthly climatic elements in the NORDKLIM data set:
Element number | Climatic element | Unit | Abbreviation |
101 | Mean temperature | 0.1 C | T |
111 | Mean maximum temperature | 0.1 C | Tx |
112 | Highest maximum temperature | 0.1 C | Th |
113 | Day of Th | date | Thd |
121 | Mean minimum temperature | 0.1 C | Tn |
122 | Lowest minimum temperature | 0.1 C | Tl |
123 | Day of Tl | date | Tld |
401 | Mean Pressure | 0.1 hPa | P |
601 | Precipitation Sum | 0.1 mm | R |
602 | Maximum 1-day precipitation | 0.1 mm | Rx |
701 | Number of days with snow cover (> 50% covered) | days | dsc |
801 | Mean cloud cover | % | N |
http://www.smhi.se/hfa_coord/nordklim
Nordklim dataset 1.0 - description and illustrations Norwegian meteorological institute, 08/01 KLIMA, 2001
## Not run: data(NordklimData) str(NordklimData) # get all the country codes countries <- unique(NordklimData$CountryCode) # earliest and latest year of data collection minFirstYear<- min(NordklimData$FirstYear) maxFirstYear<- max(NordklimData$FirstYear) allyears <- min(NordklimData$FirstYear):max(NordklimData$FirstYear) # get the yearly average of all records avgNordk <- cbind(NordklimData[,c('CountryCode','ClimateElement','FirstYear', 'NordklimNumber')], YrAvg=apply(NordklimData[,c('January','February','March','April','May','June', 'July','August','September', 'October','November','December')],1,function(x) {x[x==-9999]<-NA;mean(x,na.rm = TRUE)})) str(avgNordk) # plot the Danish mean temperatures for its 5 stations (for a quick visual # inspection, no need for labels or legends) DanavgNordk <- avgNordk[which(avgNordk$CountryCode=='DK' & avgNordk$ClimateElement==101),c('FirstYear','YrAvg','NordklimNumber')] p <- unique(DanavgNordk$NordklimNumber) for (Dp in p) { plot(DanavgNordk[which(DanavgNordk$NordklimNumber==Dp), c('FirstYear','YrAvg')],type='l',col=( which(Dp==p)), xlim=c(min(DanavgNordk$FirstYear), max(DanavgNordk$FirstYear)), ylim=c(60,120)); if (Dp != p[length(p)]) par(new=T)} # average each country avgNordkCountry=aggregate(YrAvg ~ CountryCode+ClimateElement+FirstYear , data = avgNordk, function(x) {x[x==-9999]<-NA;mean(x,na.rm = TRUE)}) str(avgNordkCountry) # plot the temperatures (mean of all stations) for each country for (country in countries) { plot(avgNordkCountry[ which(avgNordkCountry$CountryCode==country & avgNordkCountry$ClimateElement==101), c('FirstYear','YrAvg')],type='l',col=( which(country==countries)), xlim=c(minFirstYear, maxFirstYear),ylim=c(0,120), main='Mean of yearly means of all stations for each country', xlab='Years',ylab='Mean temperature'); if (country != countries[length(countries)]) par(new=T)} legend('topleft', legend = countries, col=1:5, pch=1, lty=1, merge=TRUE) ## End(Not run)
## Not run: data(NordklimData) str(NordklimData) # get all the country codes countries <- unique(NordklimData$CountryCode) # earliest and latest year of data collection minFirstYear<- min(NordklimData$FirstYear) maxFirstYear<- max(NordklimData$FirstYear) allyears <- min(NordklimData$FirstYear):max(NordklimData$FirstYear) # get the yearly average of all records avgNordk <- cbind(NordklimData[,c('CountryCode','ClimateElement','FirstYear', 'NordklimNumber')], YrAvg=apply(NordklimData[,c('January','February','March','April','May','June', 'July','August','September', 'October','November','December')],1,function(x) {x[x==-9999]<-NA;mean(x,na.rm = TRUE)})) str(avgNordk) # plot the Danish mean temperatures for its 5 stations (for a quick visual # inspection, no need for labels or legends) DanavgNordk <- avgNordk[which(avgNordk$CountryCode=='DK' & avgNordk$ClimateElement==101),c('FirstYear','YrAvg','NordklimNumber')] p <- unique(DanavgNordk$NordklimNumber) for (Dp in p) { plot(DanavgNordk[which(DanavgNordk$NordklimNumber==Dp), c('FirstYear','YrAvg')],type='l',col=( which(Dp==p)), xlim=c(min(DanavgNordk$FirstYear), max(DanavgNordk$FirstYear)), ylim=c(60,120)); if (Dp != p[length(p)]) par(new=T)} # average each country avgNordkCountry=aggregate(YrAvg ~ CountryCode+ClimateElement+FirstYear , data = avgNordk, function(x) {x[x==-9999]<-NA;mean(x,na.rm = TRUE)}) str(avgNordkCountry) # plot the temperatures (mean of all stations) for each country for (country in countries) { plot(avgNordkCountry[ which(avgNordkCountry$CountryCode==country & avgNordkCountry$ClimateElement==101), c('FirstYear','YrAvg')],type='l',col=( which(country==countries)), xlim=c(minFirstYear, maxFirstYear),ylim=c(0,120), main='Mean of yearly means of all stations for each country', xlab='Years',ylab='Mean temperature'); if (country != countries[length(countries)]) par(new=T)} legend('topleft', legend = countries, col=1:5, pch=1, lty=1, merge=TRUE) ## End(Not run)
Information about the Nordklim stations and climate element numbers.
data(NordklimStationCatalogue)
data(NordklimStationCatalogue)
A data frame with 114 observations on the following 31 variables.
Station id
Catalogue id
Station name
Height at sea level
Country
Nordklim id
Lat./Long.
Mean temperature
Mean temperature error
Mean maximum temperature
Mean maximum temperature error
Highest maximum temperature
Highest maximum temperature error
Day of Th
Day of Th error
Mean minimum temperature
Mean minimum temperature error
Lowest minimum temperature
Lowest minimum temperature error
Day of Tl
Day of Tl error
Mean Pressure
Mean Pressure error
Precipitation Sum
Precipitation Sum error
Maximum 1-day precipitation
Maximum 1-day precipitation error
Number of days with snow cover (> 50% covered)
Number of days with snow cover (> 50% covered) error
Mean cloud cover
Mean cloud cover error
The station catalogue has five columns with station information (station name, height at sea level, country code, NORDKLIM number and Lat./Long.) followed by 24 columns, two for each climate element number, the first is the first year of the dataset and the second is the last year.
http://www.smhi.se/hfa_coord/nordklim
Nordklim dataset 1.0 - description and illustrations Norwegian meteorological institute, 08/01 KLIMA, 2001
## Not run: data(NordklimStationCatalogue) str(NordklimStationCatalogue) # 114 stations length(NordklimStationCatalogue$Nordklim.number) # in 5 Nordic countries length(NordklimStationCatalogue$Country) # how many stations per country? table(NordklimStationCatalogue$Country,dnn = list("Number of stations per country")) # how many climate elements recorded per station? climElSta <- rowSums(sign(NordklimStationCatalogue[,c('X101','X111','X112', 'X113','X121','X122','X123','X401','X601','X602','X701','X801')]), na.rm = TRUE) barplot(climElSta,ylab='Climate elements',xlab='Stations', main='Climate elements recorded per station') # how many stations per climate element? staClimEl <- colSums(sign(NordklimStationCatalogue[,c('X101','X111','X112','X113','X121', 'X122','X123','X401','X601','X602','X701','X801')]), na.rm = TRUE) barplot(staClimEl,xlab='Climate element',ylab='Stations', main='Stations per climate element') # how many stations have 1,2,3, ..., 12 climate elements? # (same as Fig. 2 from Nordklim dataset 1.0 - description and illustrations) barplot(table(climElSta),xlab='Climate element',ylab='Stations', main='Number of stations as a function of number of climatic elements') ## End(Not run)
## Not run: data(NordklimStationCatalogue) str(NordklimStationCatalogue) # 114 stations length(NordklimStationCatalogue$Nordklim.number) # in 5 Nordic countries length(NordklimStationCatalogue$Country) # how many stations per country? table(NordklimStationCatalogue$Country,dnn = list("Number of stations per country")) # how many climate elements recorded per station? climElSta <- rowSums(sign(NordklimStationCatalogue[,c('X101','X111','X112', 'X113','X121','X122','X123','X401','X601','X602','X701','X801')]), na.rm = TRUE) barplot(climElSta,ylab='Climate elements',xlab='Stations', main='Climate elements recorded per station') # how many stations per climate element? staClimEl <- colSums(sign(NordklimStationCatalogue[,c('X101','X111','X112','X113','X121', 'X122','X123','X401','X601','X602','X701','X801')]), na.rm = TRUE) barplot(staClimEl,xlab='Climate element',ylab='Stations', main='Stations per climate element') # how many stations have 1,2,3, ..., 12 climate elements? # (same as Fig. 2 from Nordklim dataset 1.0 - description and illustrations) barplot(table(climElSta),xlab='Climate element',ylab='Stations', main='Number of stations as a function of number of climatic elements') ## End(Not run)