| AirQual {SwissAir} | R Documentation |
Air Quality Data of Switzerland
Description
Dataset contains Ozone (= O3), Sum of nitrogen monoxide and nitrogen dioxide (= NOx), nitrogen monoxid (= NO), ambient temperature (T), dew point (Td), wind speed (WS) and wind direction (WD) at 3 sites around lake of Lucerne in Central Switzerland in 30 min time resolution for the year 2004.
Usage
data(AirQual)Format
A data frame with 17568 observations on the following 22 variables.
startstart time of observation (GMT+1h, no day saving time)
ad.O3Mean concentration [ppb] of O3 in ad
ad.NOxMean concentration [ppb] of NOx in ad
ad.NOMean concentration [ppb] of NO in ad
ad.WSWind speed [m/s] in ad
ad.WDWind direction in ad
ad.TMean ambient temperature [deg C] in ad
ad.TdMean dew point [deg C] in ad
lu.O3Mean concentration [ppb] of O3 in sz
lu.NOxMean concentration [ppb] of NOx in sz
lu.NOMean concentration [ppb] of NO in sz
lu.WSWind speed [m/s] in sz
lu.WDWind direction in sz
lu.TMean ambient temperature [deg C] in sz
lu.TdMean dew point [deg C] in sz
sz.O3Mean concentration [ppb] of O3 in sz
sz.NOxMean concentration [ppb] of NOx in sz
sz.NOMean concentration [ppb] of NO in sz
sz.WSWind speed [m/s] in sz
sz.WDWind direction in sz
sz.TMean ambient temperature [deg C] in sz
sz.TdMean dew point [deg C] in sz
Details
The 3 sites are
- ad
Site in Altdorf is located 100 m east of motorway A2 from Bale to Chiasso, on an open field at the beginning of a more than 2000 m deep valley, at 438 m altitude.
- lu
Site is located in Sedel next to town of Lucerne 35m above and 250m south of motorway A14 on a small hill with free 360 degree panorama at 484 m altitude.
- sz
Site is located in Schwyz in an aerea of medium density of buildings next to a shopping center.
NO and O3 react in the atmosphere within seconds to NO2. The production and destruction of the sum of O3 and NO2 (= Ox) takes place on a much lower time scale, so that the spatial (cf. ipairs output) and temporal (cf. ilagplot) correlation of Ox is much more pronounced than the correlation of O3 or NO.
Note
Type of variable start has changed in Version 1.08 from factor to character to save memory.
Source
The data are collected by inNet corporation https://www.innetag.ch/ on behalf of in-Luft, an association for air quality controle, of the Swiss Cantons Aargau, Luzern, Nidwalden, Obwalden, Schwyz, Uri und Zug https://in-luft.ch/.
References
Rene Locher, Andreas Ruckstuhl; Plausibilisierung von Ozon, Stickoxiden und PM10: Statistische Methoden zur Effizienz- und Qualitaetssteigerung der Messdatenplausibilisierung; inLuft 2003
See Also
For viewing large datasets see library(IDPmisc)
Examples
str(AirQual)
sapply(AirQual, function(x) sum(is.na(x)))
cbind(min = sapply(AirQual[,-1], min, na.rm = TRUE),
median = sapply(AirQual[,-1], median, na.rm = TRUE),
max = sapply(AirQual[,-1], max, na.rm = TRUE))
if (require(IDPmisc) && require(grid)) {
## low correlation, density on logarithmic scale
ipairs(AirQual[,c("ad.O3","lu.O3","sz.O3")],
ztrans = function(x){x[x<1] <- 1; log2(x)*10})
## Not run:
ipairs(AirQual[,c("ad.NOx","lu.NOx","sz.NOx")],
ztrans = function(x){x[x<1] <- 1; log2(x)*10})
## End(Not run)
lags <- c(1, 2, 4, 8)
ilagplot(AirQual[,c("ad.O3")], set.lags = lags,
ztrans = function(x){x[x<1] <- 1; log2(x)*10})
Ox <- AirQual[,c("lu.O3","sz.O3")]+
AirQual[,c("lu.NOx","sz.NOx")]-
AirQual[,c("lu.NO","sz.NO")]
names(Ox) <- c("lu","sz")
## high correlation, density on logarithmic scale
ipairs(Ox,
ztrans = function(x){x[x<1] <- 1; log2(x)*10})
ilagplot(Ox$lu, set.lags = lags,
ztrans = function(x){x[x<1] <- 1; log2(x)*10})
dat <-
data.frame(month =as.numeric(substr(AirQual$start,4,5)),
hour = as.numeric(substr(AirQual$start,12,13)),
WD = AirQual$ad.WD,
NOx = AirQual$ad.NOx,
O3 = AirQual$ad.O3,
Ox = AirQual$ad.O3+AirQual$ad.NOx-AirQual$ad.NO)
med.dayrose <- rose(dat[,c("NOx","O3","Ox")],
subset = dat$month > 4 & dat$month < 10,
cyclVar =dat$hour, n.cyclVar = 24, circle = 24,
FUN = median, na.rm = TRUE)
## NOx (= NO+NO2) and ozone (O3) have a distinct diurnal variation of
## concentration, whereas Ox (= NO2+O3) varies only very slightly
grid.newpage()
plot(med.dayrose,
general = general.control(lwd=2),
grid =
grid.control(ray.n = 12,
circ.n = 2,
circ.sub.n = 2,
cyclVar.lab = seq(0,by = 2,to = 22)),
title =
title.control(text = "Day Rose of Medians\nduring summer time"))
grid.newpage()
plot(med.dayrose,
general = general.control(lwd = 3),
grid =
grid.control(ray.n = 12,
circ.n = 2,
circ.sub.n = 2,
cyclVar.lab = seq(0, by = 2, to = 22)),
title =
title.control(text = "Day Rose of Medians\nduring summer time"))
## exploration of upslope (North) downslope (South) wind system
## during summer time in the valley north of Gotthard
ncol <- 3
grid.newpage()
pushViewport(viewport(layout = grid.layout(nrow = 3, ncol = ncol),
width = 0.98, height = 0.98))
for (hour in seq(0, 21, 3)) {
windrose <-
rose(dat$WD, cyclVar = dat$WD, circle = 360, n.cyclVar = 32,
subset = dat$hour >= hour & dat$hour < hour+2 &
dat$month > 4 & dat$month < 10,
FUN = function(x) sum(!is.na(x)),
warn = FALSE)
pushViewport(viewport(layout.pos.col = (hour/3)%%ncol+1,
layout.pos.row = (hour/3)%/%ncol+1))
pushViewport(viewport(width = 0.9, height = 0.9))
plot(windrose,
general =
general.control(lwd = 3),
grid =
grid.control(circ.r = seq(0, 150, 50),
circ.sub.r = seq(25, 150, 25),
circ.between = -0.2,
circ.cex = 0.5,
cyclVar.cex = 0.8,
ray.lim = c(0, 150)),
title = title.control(
text = paste(hour, "-", hour+3),
between = 0.3, cex = 1.2))
popViewport(n = 2)
} ## end for
} else print("Package IDPmisc is not available")