Calidad del aire en Hermosillo y movilidad
Introducción
La población en general está expuesta a los óxidos de nitrógeno principalmente al respirarlos en el aire.
La exposición a niveles más altos puede producirse en las proximidades de las cocinas de gas, al quemar madera o combustibles o si se fuma.
Los óxidos de nitrógeno son liberados al aire por los tubos de escape de los vehículos motorizados, en los procesos de combustión del carbón, petróleo, o gas natural, y en tratamientos industriales (soldadura al arco, galvanoplastía, grabado de metales y detonación de dinamita).
Efectos
Los efectos del NO2 en la salud humana se centran sobre todo en el aparato respiratorio
Existen numerosos estudios realizados con personas que sufren asma, enfermedades crónicas obstructoras del pulmón y bronquitis crónicas que han demostrado efectos a bajas concentraciones de NO2. Estos efectos son reducción del volumen forzado de expiración o incrementos en la resistencia de las vías respiratorias.
Bibliotecas
setwd("~/R Scripts") # Directorio de trabajo.
library("pacman") # Importa biblioteca "pacman". Se utiliza para hacer una mejor gestión de paquetes.
p_load("base64enc", "htmltools", "mime", "xfun", "prettydoc","readr", "knitr","DT","dplyr", "ggplot2","plotly", "gganimate","gifski","scales", "readxl", "tidyverse","cluster", "factoextra","NbClust","tidyr") # Paquetes necesarios para la elaboración.Descarga este documento
Para fines de reproducibilidad se incluye el código para su descarga.
xfun::embed_file("U1A15.Rmd")Datos movilidad en Sonora obtenidos de Google
- Descargar datos
xfun::embed_file("MovilidadSonora.xlsx")- Leer datos
mov <- read_excel("MovilidadSonora.xlsx")Gráficas de movilidad
Venta de productos y lugares de recreación
movilidad <- ggplot(mov)+
geom_line(aes(x=date,y=retail_and_recreation_percent_change_from_baseline,colour="Venta de productos y lugares de recreación"))+
labs(colour = "") +
labs(title="Reporte de movilidad en Sonora",x="Fecha",y="Cambio sobre la línea base")
ggplotly(movilidad)En esta gráfica se puede apreciar que el cambio sobre la línea base (la cual es cero). Podemos observar que a partir del 13 de marzo la movilidad empezó a bajar hasta llegar a un valor de -71, pero a partir del 10 de abril empieza a fluctuar. Se puede apreciar que en el mes de diciembre empieza a subir y esto podría ser debido a las compras navideñas y el año nuevo, dado que el 25 de diciembre y el 1 de enero tiene una movilidad muy baja.
Tiendas y farmacias
movilidad <- ggplot(mov)+
geom_line(aes(x=date,y=grocery_and_pharmacy_percent_change_from_baseline,colour="Tiendas y farmacias"))+
labs(colour = "") +
labs(title="Reporte de movilidad en Sonora",x="Fecha",y="Cambio sobre la línea base")
ggplotly(movilidad)En esta gráfica los valores que más se destacan son el 24 de diciembre y el 31 de enero por ser muy altos. Se puede deducir que estos valores se deben a las preparaciones para las cenas de navidad y año nuevo respectivamente. Los días con movilidad más baja son el 25 de diciembre y el 1 de enero.
Parques
movilidad <- ggplot(mov)+
geom_line(aes(x=date,y=parks_percent_change_from_baseline,colour="Parques"))+
labs(colour = "") +
labs(title="Reporte de movilidad en Sonora",x="Fecha",y="Cambio sobre la línea base")
ggplotly(movilidad)La movilidad en los parques de Sonora se ha mantenido por debajo de la línea base durante la mayor parte de nuestra gráfica.
Tránsito
movilidad <- ggplot(mov)+
geom_line(aes(x=date,y=transit_stations_percent_change_from_baseline,colour="Tránsito"))+
labs(colour = "") +
labs(title="Reporte de movilidad en Sonora",x="Fecha",y="Cambio sobre la línea base")
ggplotly(movilidad)Podemos apreciar que hubo una gran disminución de la movilidad a partir del 17 de marzo y empezo una lenta subida, alcanzando sus valores máximos justamente el 23 de diciembre, justamente 33 puntos sobre la línea base. Justamente un año desde que inicio la cuarentena la movilidad se pisicionó en -31 puntos sobre la línea base.
Lugares de trabajo
movilidad <- ggplot(mov)+
geom_line(aes(x=date,y=workplaces_percent_change_from_baseline,colour="Lugares de trabajo"))+
labs(colour = "") +
labs(title="Reporte de movilidad en Sonora",x="Fecha",y="Cambio sobre la línea base")
ggplotly(movilidad)Podemos aperciar que la movilidad ha fluctuado demasiado desde el inicio de nuestra gráfica el 7 de marzo de 2020 hasta un año después, pero generalmente se mantiene por debajo de la línea base y esta bajó más en los días festivos.
Lugares residenciales
movilidad <- ggplot(mov)+
geom_line(aes(x=date,y=residential_percent_change_from_baseline,colour="Lugares residenciales"))+
labs(colour = "") +
labs(title="Reporte de movilidad en Sonora",x="Fecha",y="Cambio sobre la línea base")
ggplotly(movilidad)En esta gráfica las personas que permanecen en case ha subido más de la línea base inicial, dado a que muchas actividades se han suspendido o reducido por la pandemia. Podemos apreciar mayores valores de en los primeros meses, pero se ha mantenido fluctuando. Podemos apreciar que para el 7 de marzo de 2021 no falta tanto para llegar a la línea inicial.
Datos sobre calidad del aire en Hermosillo obtenidos de la UNAM
- Descargar datos
xfun::embed_file("Calidad_aire_hmo.xlsx")Download Calidad_aire_hmo.xlsx
- Leer datos
cal <- read_excel("Calidad_aire_hmo.xlsx")Gráfica calidad del aire
gcalidad_aire <- ggplot(cal)+
geom_line(aes(x=Fecha,y=O3,colour="O3"))+
geom_line(aes(x=Fecha,y=NO2,colour="NO2"))+
labs(colour = "Elementos") +
xlab("Fecha") +
ylab("Calidad (ppb)") +
labs(title="Calidad de aire en Hermosillo")
ggplotly(gcalidad_aire)¿Por que mejoro la calidad del aire?
La reducción en la actividad económica por inicios de la pandemia fue la razón por la que mejoró la calidad del aire, al estar los ciudadanos en sus casas, la actividades de varias industrias se han frenado y con ellas su emisión de contaminantes, las mejoras de la calidad del aire se pueden comprobar por la reducción de NO2, lo cierto que es que el aumento de la calidad del aire solo se ha dado por el paro de la actividad económica, con la reanudación de estas de manera normal la contaminación retoma los mismos niveles.
Promediode 03 por Hora
El por hora de concenctracón de O3 durante el mes de marzo del 2020
p <- read.csv("hora.csv")
horas <- ggplot(p)+
geom_line(aes(x=hora, y=PromedioHora, colour="PromedioHora"))+
labs(colour = "---") +
xlab("Hora") +
ylab("O3 (ppb)") +
labs(title="Promedio por hora en el mes de marzo de 2020")
ggplotly(horas)Promedio de O3 diario
El promedio diario de concenctracón de O3 durante el mes de marzo del 2020
p <- read.csv("promedioDia.csv")
horas <- ggplot(p)+
geom_line(aes(Fecha, Promedio_Dia) , color="red")+
xlab("Dia del mes") +
ylab("O3 (ppb)")
ggplotly(horas)Conclusión
La mejora del aire fue una buena noticia, pero lo cierto es que se dio solo por el drastico paro de la actividad economica, esto fue una oportunidad para no solo la industria, tambien para los ciudadanos para intentar emitir menos contaminantes al aire y asi siquiera no permitir que los niveles se volvieran a elevar tanto, sin embargo se ha visto como con el pasar el tiempo, despues de los primeros dias de cuarentena las actividades retomaron lentamente su ritmo, contaminando al final de igual manera.
Fuentes
Noticias ONU. (2020, 22 abril). La cuarentena por el coronavirus mejora la calidad del aire, pero no. https://news.un.org/es/story/2020/03/1471562#:%7E:text=Las%20observaciones%20han%20demostrado%20que,cuarentenas%20en%20China%20e%20Italia.&text=El%20NO2%20permanece%20en%20la,otros%20gases%20en%20la%20atm%C3%B3sfera.
Consejería de Salud de la Región de Murcia. (2010). Dióxido de nitrógeno. marzo 12, 2020, de Consejería de Salud de la Región de Murcia Sitio web: https://www.murciasalud.es/pagina.php?id=180252&idsec=1573#