# Almacenamientos al 1 de octubre de las Presas Abelardo L. Rodríguez
# y El Molinito, información recabada en CONAGUA Oficinas de Hermosillo
# el día 28 de Agosto del 2014.
##################################################################
# Limpio mi espacio de trabajo
rm(list=ls())
# Cargo este programa
library(ggplot2)
##################################################################
# Almacenamiento reportados al mes de Octubre, en Hm3, del 2000 al 2017,
molinito<-c(147.889223, 26.031212, 41.323024, 79.051130,29.638671, 10.086477, 8.495519,
20.35, 12.17, 7.984, 13.387, 5.213, 8.454, 41.224, 36.262,
37.254, 8.628, 36.278, 13.130, 13.357, 17.367, 90,26.3, 25.194, 19.759
)
alr<-c(144.978338, 139.098043, 249.053681, 136.863843, 59.370044, 3.458480, 3.398603,
1.63, 13.85, 0, 7.687, 2.284, 0, 32.923, 2.85, 0.006, 0, 2.379, 0.2,
1.8, 3.2, 24.52, 22.79, 1.01,0.003)
length(alr)
# los años que me interesa graficar
years<-c(1993:2017)
length(years)
##### Figura 2
df3 <-data.frame (year=years, ALR= alr, RFV=molinito)
df3
p <- ggplot(df3, aes(year), height=600, width=1000) +
geom_area(aes(y = RFV), fill="green", alpha = .9)+
geom_area(aes(y = ALR), fill="black", alpha = .8)+
geom_line(aes(y = RFV), size=1.5, linetype="dashed", colour = "white")+
geom_line(aes(y = ALR), size=0.1)+
scale_x_continuous(breaks = c(1993:2017)) + xlab("A\u{F1}os") +
ylab(parse(text='Hm^3'))
p= p + theme(axis.text = element_text(colour = "black", size=16, angle = 90),
axis.title.x = element_text(colour="black",size=16,angle=0,hjust=.5,vjust=0,face="plain"),
axis.title.y = element_text(colour="black",size=18,angle=90,hjust=.5,vjust=.5,face="plain")
)
# elaboramos manualmente la leyenda
rect2 <- data.frame(xmin=2009, xmax=2010, ymin=220, ymax=230)
rect3 <- data.frame(xmin=2009, xmax=2010, ymin=200, ymax=210)
####Tandeos
tandeo1 <- data.frame(xmin=1998, xmax=1999, ymin=135, ymax=150)
tandeo2 <- data.frame(xmin=2005, xmax=2006, ymin=135, ymax=150)
tandeo3 <- data.frame(xmin=2010, xmax=2012, ymin=135, ymax=150)
p + geom_rect(data=rect2, aes(xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax),
fill="green", size=1.2, alpha = .9, linetype="dashed", color="white",
inherit.aes = FALSE) + annotate("text", x = 2011.2, y = 225, label = "El Molinito")+
geom_rect(data=rect3, aes(xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax),
fill="black", alpha = .8,inherit.aes = FALSE) +
annotate("text", x = 2012.5, y = 205, label = "Abelardo L. Rodríguez.") +
geom_rect(data=tandeo1, aes(xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax),
fill="red", alpha = .9,inherit.aes = FALSE) +
geom_rect(data=tandeo2, aes(xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax),
fill="red", alpha = .9,inherit.aes = FALSE) +
geom_rect(data=tandeo3, aes(xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax),
fill="red", alpha = .9,inherit.aes = FALSE)+
annotate("text", x = 2015, y = 140, label = "Tandeos (racionamientos)")
##################################################################
##################################################################
Gráfica de líneas y áreas sólidas (2)
Dr. Luis Alan Navarro Navarro
Centro de Estudios en Gobierno y Asuntos Públicos
El Colegio de Sonora
La siguiente gráfica es una versión actualizada de la que se encuentra en la Figura 2 de mi artículo titulado Cambios en el paisaje arbolado en Hermosillo: escasez de agua y plantas nativas (2016) que hace referencia al almacenamiento de las presas Abelardo L. Rodríguez (ALR) y Rodolfo Félix Valdez ("El Molinito") localizadas sobre el río Sonora, en el Noroeste de México, "El Molinito" ubicada aguas arriba de la presa ALR.
Adicionalmente, la gráfica marca con rectángulos de color rojo, los años en los que la ciudad ha sufrido "tandeos" o racionamientos de agua
(Haro-Velarde y colaboradores 2016).
Los tandeos han sido la medida más drástica en la que la población ha sentido la escasez de agua en cuenca del río Sonora, reflejada en los bajos volúmenes de almacenamiento de las presas.
Abajo se muestra el código de R con la información necesaria para replicar esta gráfica.
R code
