Openair 오존 8시간 평균 농도 계산

오존의 대기환경기준은 8시간 평균 60 ppb, 1시간 평균 100 ppb 이하입니다. 

에어코리아에서 제공하는 자료는 1시간 평균(오후 2시 농도는 오후 1시부터 2시까지 5분 간격으로 측정한 값을 평균한 값이며, 오후 2시 당시 농도가 아님)이므로 이 자료에서 100 ppb를 넘는지 바로 확인할 수 있습니다. 

출처: 대기환경연보 부록 14

 

8시간 평균 농도는 연구자가 직접 계산해야 합니다. Openair에서 아래 코드를 이용하면 8시간 평균 농도가 새 열(컬럼)에 생성됩니다. 

The openair book - 26  Utility functions (bookdown.org)

 

mydata <- rollingMean(mydata, pollutant = "O3", width = 8, new.name = "rollingo3", data.thresh = 75)

 

울산 모 측정소 2019년 자료로 오존 그래프를 그렸습니다. 

1시간 기준은 일부 시간대에 초과했지만 8시간 기준은 봄과 여름에 많이 초과했습니다.

 

위 그래프 코드는 다음과 같습니다.

ggplot(data = mydata, aes(x = date, y = O3)) + geom_line(color = "blue") + labs(title = "Ozone (1-hour mean)", x = "Date", y = "Concentration (ppb)") +  geom_hline(yintercept = 100, color = "red", linetype = "dashed")