Universidad de Costa Rica Facultad de Ciencias Sociales Escuela de Geografía
Universidad de Costa Rica Facultad de Ciencias Sociales Escuela de Geografía
Escuela de Geografía
Setiembre, 2023
1
Tribunal examinador
7
~ Dr. Acfolfo Quesada•Román
Lector
t
2
Agradecimientos
A mi mamá Norma, por sembrar en mí la curiosidad por la naturaleza y los ríos. A Ricardo
por apoyarme en todo momento y motivarme a estudiar, a Andrey, mi primo por ayudarme
A Cristian Aguilar, Omar Mejía, Tomás Fernández, Jesús y Manuel Céspedes, por las
Sharon Camacho, Ana María Araya, María José Carpio y Sara Blanco, por la consejería,
A mis amigos de Aserrí: Alejandro Fallas, Danny Mora y Luis Carvajal por ser amigos
A Jasson Mora y Helber Guillén y Manuel Peralta por el apoyo durante los muestreos y
trabajo campo. A Osvaldo Ureña, por al apoyo y la motivación para seguir estudiando.
Quesada por los consejos y las oportunidades de aprendizaje. Ana María Durán del
3
Dedicatoria
A Rex, que me acompaño la mitad de mi vida, en las noches, tardes y días de escritura y
estudio.
A mi mamá Norma.
A mis amigos, Omar, Cristian, Tomás, Jesús, Manuel, Danny, Luis, Alejandro, Andrey.
4
Tabla de contenidos.
5
3.3 Balance hidrológico simulado: .............................................................................. 57
Capítulo 4. Evaluación de influencia de las variables hidrogeomorfológicas e isotópicas
en el síndrome de la corriente urbana en la microcuenca media-alta de Quebrada
Negritos............................................................................................................................ 62
Cambios en el paisaje: geomorfología y usos de suelo: .............................................. 62
Cambios hidrológicos: ................................................................................................. 65
Capítulo 5. Discusión de resultados, conclusiones y recomendaciones.......................... 70
5.1 Discusión de resultados: ........................................................................................ 70
5.2 Conclusiones. ......................................................................................................... 75
5.3 Recomendaciones: ................................................................................................. 78
Referencias................................................................................................................... 80
Lista de tablas.
Lista de figuras.
Figura 6. Población de los distritos del cantón San Pedro con base en datos del INEC 1973,
6
Figura 8. Terrazas Aluviales (TA). .................................................................................. 44
Figura 15. Balance hidrológico diario simulado para el periodo de estudio 2019, Quebrada
Figura 16. Aportes al caudal total simulado de diferentes fuentes de agua, Quebrada
Negritos, 2019.................................................................................................................. 60
Figura 17. Modificaciones del paisaje por cambios en el uso de suelo. .......................... 64
Figura 18 Figura 18. Fracciones de agua nueva respecto al caudal muestreado. ........... 69
Lista de abreviaturas.
7
Gwcontr Contribución de agua subterránea
IMN Instituto Meteorológico Nacional.
INEC Instituto Nacional de Estadística y Censos.
INEC Insitutito Nacional de Estadísticas y Censo
LBP Ladera de pendiente baja
LEL Línea Evaporativa Local
LPM Ladera de pendiente moderada
LSC Ladera escarpadas
MDS Modelo Digital de Superficie
MIVAH Ministerio de Vivienda y Asentamiento Humanos.
MJO Oscilación Madden-Julien
ODS Objetivos de Desarrollo Sostenible
ODS Objetivos de Desarrollo Sotenible
ONU Organización de las Naciones Unidas.
P Precipitación
QN Quebrada Negritos.
Qs Caudal simulado
Qup Escorrentía Superficial
SIG Sistema de Información geográfica
TA Terrazas Aluviales
TI Área total impermeabilizada
UCR Universidad de Costa Rica
ZCIT Zona de Convergencia Intertropical
8
Resumen
producido cambios en el uso del suelo que han alterado las funciones de regulación y el
afectando los flujos de base en los ríos urbanos. Además, la expansión de la red vial y el
están condicionados por las formas del relieve, promueven la escorrentía superficial hacia
impermeable y solo el 17% es permeable. Esto implica una conversión rápida de la lluvia
total del caudal, lo que conlleva problemas como la erosión lateral de las laderas del río e
QN.
9
Capítulo 1. Marco introductorio
1.1 Introducción
En las cuencas urbanas, los cambios en el uso del suelo, el incremento de la red pluvial y
desfogan directamente al canal. El aumento del caudal y los flujos pico durante las
concepto del síndrome de la corriente urbana (SCU) en adelante SCU. Este concepto,
principales factores de afectación y cambio del ciclo hidrológico en ríos urbanos (Walsh
et al., 2005).
El propósito de esta investigación fue conocer cuáles factores del SCU se encuentran en
los miembros finales o las fuentes geográficas de agua, y el nexo del régimen hidrológico
10
con el uso del suelo y la geomorfología. Con ello, se pretende detallar los mecanismos del
alcantarillado pluvial, los cuales han determinado los patrones de escorrentía que inciden
deslizamientos en la QN.
Negritos?
1.3 Problema
2016). La escorrentía rápida y directa a los canales y afluentes principales puede provocar
almacenamiento de agua en el suelo y los acuíferos disminuyen, lo que modifica los flujos
base en los ríos. El desfogue de aguas residuales y pluviales también puede afectar
2009).
Estas alteraciones del ciclo hidrológico abordadas desde el SCU permiten comprender las
11
elementos como la geomorfología, hidrología y la ecohidrología. Debido a la amplitud de
et al., 2016, 2018; Francis, 2012; Glenn et al., 2015; Grimmond & Oke, 1999; Poff &
Ward, 1989).
cauce. Además, los grandes volúmenes de escorrentía superficial conducidas por la red
inundación relacionados con la expansión urbana en las zonas cercanas a la cuenca del río
QN (CNE, 1990). Así Además, se han evidenciado problemas de erosión que afectan la
12
impactado en los componentes locales del ciclo hidrológico en la zona. Este tipo de
potencial para utilizarse como indicador de eficacia de infraestructura verde para el control
de inundaciones en zonas urbanas, así como cuencas de mayor tamaño o incluso a nivel
cantonal.
1.4 Antecedentes
factores afectan la forma del canal simplificándolo y estrechándolo (Vietz et al., 2015).
superficial (Wenger et al., 2009). Además, el efecto del confinamiento, los cambios en su
Para estimar los efectos de la urbanización, y determinar los tiempos de estancia y edades
del agua en una cuenca urbana, relacionado a aspectos como la influencia de los usos de
13
de escorrentía superficial. En función de los aportes de diversas fuentes de agua, distintas
Oxígeno 18 (18O) y el Deuterio (2H) para conocer la dinámica hidrológica y sus procesos
desarrollado modelos con los que se pueden discriminar las fuentes geográficas del agua
(en adelante FGA) mediante análisis hidroquímicos, hidrométricos e isotópicos (Birkel &
Soulsby, 2015; Hooper, 2003; McDonnell et al., 2010; Ogle et al., 2014).
A través de los modelos de mezcla de fuentes geográficas del agua (EMMA, por sus siglas
en inglés), se plantea que una corriente de agua es una combinación de varias FGA. Por
Los modelos de mezcla se han usado ampliamente para conocer la concentración isotópica
e hidroquímica del agua y así evaluar el estado ecohidrológico de una cuenca al emplear
la concentración de solutos como indicador para la escogencia de las FGA (Correa et al.,
El incremento del área impermeable tiende a saturar suelos poco profundos, lo que
ser tan eficaces. Estas medidas no sustituyen la función de retención de las áreas naturales
14
o permeables, pues solo atienden una pequeña porción del área impermeabilizada y
convierten la escorrentía en una descarga directa hacia las aguas superficiales, lo que
infraestructura verde, se han empleado isótopos estables del agua para evaluar las medidas
de control de aguas pluviales y sus efectos en los aportes de agua en una cuenca con un
40% de su área urbana en Carolina del Norte (Jefferson, 2021). Se determinó que la firma
principal. Esto se debe a mayores niveles evaporación y una mejor mezcla de agua dentro
de las piscinas. Además, se estimó que las piscinas aportan al caudal de hasta un 10% del
flujo total durante la crecida del río en el momento de mayor caudal. Adicionalmente, las
piscinas contribuyeron en un promedio de entre 32% y hasta un 54% del flujo en el arroyo,
posterior a la crecida del río. Esta investigación refleja que las piscinas y biorretenedores
tienen un comportamiento adecuado para aumentar los tiempos de residencia del agua en
Otros enfoques de investigación señalan que se deben replantear las relaciones entre la
estructura y función ecológica del paisaje, en relación con los usos del suelo en zonas
por lo que usualmente se utiliza un porcentaje del área total impermeabilizada (TI) como
15
cambios en la morfología del canal, disminución de la riqueza ecológica, la simplificación
y homogeneización del ecosistema, como los más importantes (Walsh et al., 2005).
Al seguir las líneas de análisis del SCU, se toman en cuenta los cambios en la dinámica
clave en la manifestación del síndrome (Booth et al., 2016; Groffman et al., 2003, 2005).
Precipitación
t + Evapotranspiración
t -Evapotranspiración
i
Figura 1. Alteraciones del ciclo hidrológico en cuencas urbanas. A la izquierda se muestran los
procesos hidrológicos en equilibrio en cuencas naturales o forestadas. A la derecha se detalla como
el proceso de urbanización afecta el equilibrio de los procesos hidrológicos.
16
El uso de la técnica de separación de hidrogramas ha sido ampliamente utilizado. En el
realizadas en cuencas con poca alteración en climas templados y fríos. Estos estudios se
nieve. Sin embargo, al revisar los estudios realizados en cuencas urbanas, se ha encontrado
que las fracciones de agua nueva son superiores al 60%. Esto se debe a la impermeabilidad
del suelo, lo que produce una mayor escorrentía superficial en lugar de infiltración en el
suelo. Estas fracciones de agua pueden variar según las condiciones específicas de la
tener en cuenta las limitaciones asociadas con el uso de isótopos estables en la separación
espacial y temporal de los isótopos, la combinación de fuentes de agua y los altos costos
de muestreo y análisis.
urbana, sin dejar de lado las variables hidrológicas sobre la percepción de usos del agua y
participación comunitaria en la gestión del recurso (Capps et al., 2016; Pataki et al., 2011).
oferta y demanda de agua, basados en mapeo de usos del suelo, combinado con entrevistas,
buscan conocer la cantidad de agua usada para regar zonas verdes o campos agrícolas
Estudios relacionados con el SCU en Costa Rica, Gervasi (2012) analizó el caso de la
17
densidad de viviendas y caminos durante el período 1989-2010, determinó que se
acrecentó el pico de caudal debido a los aportes inmediatos de escorrentía superficial. Este
del canal (bosque ripario), cuando esta pasó de 16 metros de ancho a 100 metros de ancho.
Para la QN, se encuentran las investigaciones de Chaves (2018) donde cuantifica y modela
hidráulico sin calibración y la subsecuente zonificación de estas amenazas. Por otra parte,
Mora et al. (2015) hacen una caracterización geotécnica de las laderas de la QN que
atraviesan la Universidad de Costa Rica, para esos fines determinaron las zonas con mayor
1.5 Justificación
al caudal. Los aspectos morfológicos y los cambios en el uso del suelo se manifiestan en
una escorrentía superficial de alta intensidad y frecuencia, lo que lleva a un aumento súbito
de los caudales de los ríos urbanos. Esto puede convertirse en una amenaza para la
18
seguridad y la salud de las personas, así como para los ecosistemas y las obras civiles
urbanas en el 2019 fue de 56%, el mismo ente estimó para Costa Rica en ese mismo año
Nación (2019) señala que en el Gran Área Metropolitana (en adelante GAM) en 1986
habían 66.3 habitantes por hectárea (ha) urbana incrementando a un 75.9 para el 2019. La
y aspectos morfológicos (Camacho Lizano et al., 2004; Chaves, 2018; Mora et al., 2015,
Diversas investigaciones sobre el SCU toman algunos de los efectos del cambio del
hidrología isotópica con el marco interpretativo del SCU, abre una perspectiva para el
estudio de los cambios hidrológicos en ríos urbanos (Fletcher et al., 2011a; Vietz et al.,
institucionales y legales locales, como la Constitución Política remite al Artículo 50: “El
Estado procurará el mayor bienestar a todos los habitantes del país, organizando y
19
derecho a un ambiente sano y ecológicamente equilibrado. Por ello, está legitimada para
denunciar los actos que infrinjan ese derecho y para reclamar la reparación del daño
Desarrollo Sostenible (en adelante, ODS) que conciernen al ambiente el agua y las
ciudades. A saber los ODS 6 ,11, 13 y 15 que buscan garantizar la disponibilidad de agua
y su gestión sostenible y el saneamiento. Por otro lado, pretenden lograr que las ciudades
adopción de medidas para combatir el cambio climático y sus efectos, así como velar por
1.6 Objetivos
20
2). Evaluar el régimen hídrico en términos de cantidad de agua de diferentes
o de agua en relación con la pendiente, orientación del terreno y la rugosidad del terreno
(Davie, 2008). Además, estudia los materiales, formas y procesos de formación del
relieve, siendo un área importante para la planificación y para el análisis del paisaje (Otto
Geográfica (SIG) que impulsó el desarrollo de modelos y otros productos que sirven para
el análisis multivariado del terreno (Otto et al., 2018). Para ello se emplean los archivos
21
tipo ráster, estos son representaciones numéricas del espacio que contienen información
compuestas por píxeles de igual tamaño, ordenados en una matriz en columnas y filas que
contienen la información de las entidades y sus coordenadas (ESRI, s/f-a; McDonnell &
Kemp, 1995). Dentro de los archivos ráster se encuentran los Modelos Digitales de
Superficie (MDS) como el AW3D Estándar, con una resolución espacial de 1x1 metros.
contemplan los objetos que sobresalen del terreno como infraestructura, vegetación o
cualquier otro objeto (ESRI, s/f-b). El potencial del software SIG y la información de alta
pequeña y a la vez extraer información derivada del procesamiento del (MDS), como
como insumos para la caracterización geomorfológica. Esto permite entender los procesos
áreas como medio ambiente, ecología, geografía e hidrología (Otto & Smith, 2016).
Como procesos hidrológicos, entendemos una serie de fenómenos físicos donde los flujos
e intercambios energéticos del ciclo hidrológico, que se encuentran supeditados a los usos
22
escorrentía y almacenamiento en el suelo (Naghettini, 2016). Para el cálculo de caudal se
acude al método área-velocidad, para el cual se tiene una sección transversal del canal
Q= (V x A) Ec.1
temperatura, precipitación y velocidad del viento. Las características del suelo como
textura, estructura, color, y densidad, así como las características de la vegetación, altura
ρCp δe
Rn ∆ +
r
Ea = λ⌊∆+γ(1+r a/r )⌋ Ec.2
s a
específico del aire a presión constante (J kg-1 °C-1), de es el déficit de presión de vapor
23
aerodinámica en s m-1, que es función de 208/U2 donde U2 es la velocidad del viento a
del área foliar (Dunn & Mackay, 1995). Para ra, debido a la ausencia de datos, se utilizó
meteorológica utilizada en este estudio (CIGEFI, UCR) está situada en una zona urbana,
a 245 metros del área de estudio y se considera la simplificación válida de Kc ~1 para que
Ea ≅ ET0.
Otra de las formas para conocer la dinámica hidrológica, es mediante los trazadores
isotópicos, estos muestran los procesos de formación de agua en los diversos estadios del
ciclo hidrológico. Según el Instituto Mexicano de Tecnología del Agua (2009) (en
adelante IMTA) la hidrología isotópica es una técnica que utiliza isótopos estables,
radiactivos naturales y artificiales que sirven para identificar los distintos cambios
fisicoquímicos del agua dentro del ciclo hidrológico. Los isótopos son átomos cuyos
autor muestra que los isótopos tienen tres tipos fundamentales: radioactivos, radiogénicos
isótopos estables (18O, 13C, 34S) la abundancia de estos en el agua responde a variaciones
hidrología isotópica, pues las concentraciones en el agua nos indican bajo cuales procesos
24
y estadios pasaron las partículas de agua dentro del ciclo hidrológico (Instituto Geológico
Los isótopos fungen como un indicador sobre la dinámica de una cuenca en términos del
tiempo de estancia del agua y los procesos de mezcla de diversas fuentes geográficas de
Deuterio (δ2H), estos son utilizados en estudios hidrológicos, debido a su origen natural y
a sus propiedades conservativas (Birkel et al., 2021; Correa et al., 2019; James & Roulet,
2006). Además permiten observar la variación isotópica de las fuentes geográficas del
Respecto a la separación del hidrograma, este se basa en una separación de dos o más
balance de masa, donde se asume que, la mezcla de agua del caudal, está constituida por
diversas fuentes geográficas de agua (lluvia, agua de suelo, agua subterránea, aguas
pluviales). Esta separación (Ec. 3 – 5), se aplica a través de un balance de masa con al
menos 2 FGA identificadas para extraer el aporte de cada fuente al total del caudal
Qt = Qn + Qv Ec. 3
Ct = cn + cv. Ec. 4
25
Donde Qt es la mezcla de aguas en el caudal total, Qn es aporte de agua nueva (lluvia),
además Qv es el aporte de aguas viejas (flujo base). El caudal total Qt y sus componentes
menos dos FGA y su contenido isotópico (δn y δv) podemos inferir el aporte de estas a la
mezcla (δt). Por lo tanto, podemos conocerlos mediante la unión de las ecuaciones 3 y 4
𝑸𝒏 (𝜹𝒕)𝜹𝒗)
= (𝜹𝒏)𝜹𝒗)
Ec. 5
𝑸𝒕
Donde Qt es el caudal total, (δt) trazador isotópico del agua del agua de la corriente, (δv)
el trazador isotópico del agua vieja (flujo base), (δn) trazador isotópico del agua nueva
Este balance de masa nos permite conocer la fracción de agua nueva en el caudal en
almacenada. El SCU asume que en cuencas urbanas (Figura 1), la modificación de las
descargas directas de aguas pluviales por el aumento del área impermeabilizada de suelo
que incide en un aumento en la frecuencia y magnitud de los flujos pico en el caudal. Con
crecidas más repentinas, al tiempo disminuyen los flujos base por la disminución de agua
almacenada en el suelo (Booth et al., 2016; Hession, 2004; Walsh et al., 2005).
26
incidiendo en una mayor cantidad de nutrientes en el río, propiciando el crecimiento de
causando eutrofización en cuerpos superficiales de agua (Booth et al., 2016; Walsh et al.,
2005) y su relación con la salud ambiental y equilibrio ecológico del río (UNESCO, 2018;
(en adelante GIRH) que incluye actores como la academia, sociedad civil y tomadores
modelos lluvia-escorrentía permiten conocer entre otros aspectos cómo responde una
los cuales se estiman los flujos de agua hacia un río principal formando su caudal (Figura
2) para un esquema conceptual del modelo aplicado en este estudio). Dichos flujos se
27
En este estudio se utiliza un modelo de dos reservorios que conceptualiza el suelo y el
acuífero como fuentes de agua basado en los coeficientes de almacenamiento lineal que
almacenamiento lineal del flujo lateral en los humedales ribereños. y se calculan los
. FIUJOS
superficiales
(Qs)
Flujos de agua
de suelo
laterales (Qup)
Flu¡os de agua
subterránea
(Gwflow)
Figura 22.Modelo
Modelo
Modeloconceptual
conceptual
conceptuallluvia-escorrentía
lluvia-escorrentía.
lluvia-escorrentía.Modelo conceptual de dos reservorios que detalla las
entradas (P) y salidas (ET y Qsim) para el cálculo del balance hidrológico.
Mediante los balances de masa, se calculan los flujos de agua, la entrada como la
definido con base en el percentil 99 de los registros de lluvia horaria de la estación CIGEFI
metodología que contempla trabajo de campo con el uso de sensores, análisis de series de
de Quebrada Negritos, San Pedro de Montes de Oca serán procesados con Sistemas de
modelo digital de superficie (MDS) “AW3D Estándar” con una resolución espacial de
1x1 metros, ortofotos a escala 1:1000 de los años 2015-2018 suministradas por Sistema
29
2019. Además, información vectorial como curvas nivel cada 5 metros del SNIT y un
procesamiento del (MDS) AW3D, del se extraerán los ráster de acumulación y dirección
2016; Otto et al., 2018; Quesada-Román, 2018) que consisten en 3 fases. La primera de
un mapa preliminar con las principales unidades identificadas, y finalmente el post mapeo,
setiembre del 2019, donde se realizaron aforos manuales en una sección definida de 2,6
30
(GLOBAL WATER FP111), la velocidad y altura de la columna de agua mediante marcas
estos eventos aforados se analizó por isotopos estables, mediante espectrometría láser de
un instrumento de tipo PICARRO L2120-i δD/δ18O Isotopic Water Analyzer con una
Con base en los datos recolectados y suministrados, se elaboró una matriz de datos en el
los aforos por medio de dos criterios estadísticos (RMSE y MAE). El promedio de la raíz
cuadrada del error cuadrado (RMSE, siglas en inglés) y el error absoluto (MAE, siglas en
algorithm, NSGA2, siglas en inglés) para minimizar los errores hacia una solución óptima.
31
El algoritmo usó un total de 10,000 iteraciones para encontrar una solución óptima en
mapa de cobertura con una clasificación basada en dos zonas (permeable e impermeable),
respectivo al área total de la QN. Este mapa se elaboró con base en fotografías aéreas
premapeo para poder conocer las condiciones del área y posteriormente hacer
la visita campo.
fenómeno observado en esas fuentes y con los resultados analizados en el área de estudio.
Estos análisis buscarán evaluar cuál es el aporte de las FGA, en el tiempo de respuesta de
32
Capítulo 2. Caracterización geomorfológica de la sección
media-alta de la microcuenca de la Quebrada Negritos, San
Pedro de Montes de Oca.
2.1 Área de estudio:
distritos San Pedro y Sabanilla del cantón Montes Oca, provincia San José. Su afluente
principal tiene una longitud de 1495 metros y es tributario del Río Torres. La QN discurre
entre valles fluviales de baja disección con pendientes entre los 8° y 16° (Camacho Lizano
et al., 2004). La alta permeabilidad de la roca y los suelos propician una alta probabilidad
zonas de pendiente media y baja. (Camacho Lizano et al., 2004; Chaves, 2018; Mora et al.,
33
496020E
N
• CIGEFI
O 100 200 m
496020E
Contexto
Nacional
Simbología
Puntos de muestro D Distritos
e CIGEFI
Altimetría
Punto de aforo (m snm)
1.330
- - Drenaje Principal
1.120
O cuenca
Sistema de coordenadas Proyectadas: CRTMOS - Datum: WGS 84
Figura 3. Mapa de ubicación de sitios de muestreo de la Quebrada Negritos, dentro de los distritos
San Pedro y Sabanilla de Montes de Oca.
34
2.2 Clima:
diciembre a marzo, con abril como mes de transición a lluvioso, y la temporada de lluvias,
región es de 2322 mm anuales con 148 días con lluvias, siendo los meses de setiembre y
octubre los de mayor precipitación (> a 200mm). Las temperaturas promedio oscilan entre
los 26.7 °C y 16.8 °C (IMN, 2020; (IMN), 2021)(IMN, 2020; (IMN), 2021).
Por otro lado, el sitio se ve influenciado por las fases del ENOS. El fenómeno El Niño
representa la fase cálida, mientras que La Niña representa la fase fría, con la fase neutral
caracterizada por oscilaciones de temperatura entre -0,5 °C y 0,5 °C. En el caso de Costa
sinóptica como la Oscilación Madden-Julian (MJO) , que opera en ciclos que abarcan de
región. La influencia de los eventos de El Niño y La Niña se reflejan en los años con
entre abril y junio y los más fríos en noviembre, diciembre, enero y febrero
Los años deficitarios, guardan relación con el fenómeno El Niño, en 2009 (-17%), 2011
(-5%), 2012 (-43%), 2013 (-16%), 2014 (-26%), 2015 (-32%), 2018 (-17%) y 2019 (-7%),
siendo el año 2012 y 2015 los de mayor déficit. A pesar de que el índice El Niño 3-4
indicara que la fase neutra fuera la de mayor incidencia, a nivel nacional durante 2012,
ambas vertientes sufrieron una sequía importante que ocasionó una disminución
entre los 20,1 °C y 25,7 °C, lo que ocasionó una disminución en la precipitación.
La influencia de La Niña es notoria en los años con superávit, 2006, 2007, 2008, 2010,
2016 y 2017, este superávit de precipitación fue de (32,5%), (60%), (28,8%), (3,3%) y
(34,8%) respectivamente. El superávit más bajo, durante el 2016, tuvo una fuerte
influencia de El Niño durante los primeros cuatro meses del año. Por otra parte, el año
2008 el acumulado anual fue de 2802 mm, un 60% más respecto al promedio anual de la
promedios de precipitación de enero a abril de 11,3 mm, 6,5 mm, 8,5 mm y 44,9 mm,
respectivamente, lo que muestra una transición gradual hacia la temporada lluviosa. Este
36
estableciéndose la primera parte de la temporada de lluviosa. Durante julio, por efecto de
temporada lluviosa, siendo setiembre y octubre los meses más lluviosos con promedios
del sitio tienen poca variación diaria donde el percentil 99 de los promedios diarios es de
22,62 °C. La temperatura promedio diaria de toda la serie oscila entre la mínima 14,60
°C y máxima de 26,08 °C. El promedio de los meses más cálidos se presenta en abril
20,7 °C, mayo 20,9 °C y junio con 20,7 °C. Así mismo las temperaturas más bajas se dan
50 500
40 400
30 300
20 200
~ 10 100 ~
C)
;s·
~ o
:, o ~
-;¡¡
- -
C)
ai
c.
E 50
2013 20 14
- 20 15 2016 20 17 2018 2019
5003
a:
:,
~ 2-
40 400
30 300
20
-- ... ,,,..... - 200
10
o
EFMAMJJASOND EFMAMJJASOND
1
1
J
EFMAMJJASOND
1
EFMAMJJASOND
Mes
EFMAMJJASOND EFMAMJJASOND EFMAMJJASOND
1
100
37
en los meses de noviembre 19,6 °C, diciembre 19,2 °C, enero 18,9 °C y febrero 19,2°C.
El plan regulador del cantón de Montes de Oca, estimó que en los años 1989 y 1998 hubo
una disminución en las áreas destinadas a cultivo y un aumento del uso forestal y urbano.
hectáreas urbanas en los años 1989 y 1998 respectivamente. Hacia 1998, el uso urbano es
de un 84% en el distrito San Pedro y un 66% en Sabanilla (PRODUS, 2001, pp. 2–4).
un 75%.
la red vial y el 17% restante, áreas permeables como bosques y pastos. La distribución del
red de carreteras.
cuenca y en las orillas del río destacando las áreas de protección de la Universidad de
Costa Rica y los parques que se ubican a lo largo de la ribera. En el resto de la cuenca, los
parches más pequeños de áreas verdes, como parques, patios o zonas verdes, se encuentran
38
4960Z0E
Usos de la tierra
Urbano/Impermeable -
Bosques y pastos/Penneable -
Carreteras -
Simbología
Cuenca
Red fluvial
Quebrada Negrttos -
~óecoordeNdas~. CflTM05•
01tu111; wc;s114
F\.lel'(ti:$dll!fflilleOOl'IIIOlll'tOffl'11(:i1611
T~1(5NIT}
Contexto Nacional
(t.boralopor;Gcóg, "'-1hcdSolltSC.,!1rv
""
2.4 Población:
Montes de Oca en 1973 fue de 27.536 habitantes y para el 2011, 49.132, lo que representa
los censos 1973, 1984, 2000 y 2011 han colocado a San Pedro como un distrito urbano,
En el caso de Sabanilla para el Censo de 1973 se consideraba como un distrito rural con
una población de 3475 personas y para los años 1984, 2000 y 2011 cambió a distrito
20000·
(/)
2e
2 Distrito
:.ol'IJ
.i:::
■ Mercedes
Q) ■ Montes de Oca
"O ■ Sabanilla
"O
l'IJ ■ San Pedro
"O ■ San Rafael
'E 10000·
l'IJ
u
O·
Figura 6. Condición geográfica de la población ubicada en los distritos del cantón San Pedro con
base en datos del INEC 1973, 1984, 2000 y 2011.
drenaje artificial. Estos cambios aumentan la conectividad con el canal principal, estos
procesos inciden en picos de caudal más abruptos, pero de menor duración que generan
40
tienden a disminuir los flujos de base debido a la incapacidad del sistema de infiltrar agua
de la Cordillera Volcánica Central. La capa más superficial está constituida por mantos de
cenizas y suelos Andisoles con texturas limo-arenosas que son fácilmente erodables
(Camacho Lizano et al., 2004). Estos depósitos conocidos como lavinas se caracterizan
por presentar un espesor variable. Estas avalanchas y flujos de detritos dan lugar a colinas
alargadas que se extienden en dirección este-oeste, como resultado del colapso del
Estas unidades sobreyacen los depósitos laháricos que se pueden visualizar en el fondo de
pendientes entre 0,4°-20° y una media de 4,41°, lo que permite distinguir diversas
los nombres de las algunas unidades cercanas al área de estudio propuestas por Camacho
Lizano et al. (2004): los fondos de valle (FV), fondos de valle confinados (FVC), terrazas
aluviales (TA), laderas de pendiente moderada (LPM), laderas de baja pendiente (LBP) y
41
2.5.1 Unidades
forma
la parte alta de la cuenca. Esta unidad posee una pendiente media de 3,88°, donde las
menor incidencia en la profundidad y ancho del canal. El punto de menor altitud (aforo)
con un relieve de baja energía, donde convergen los flujos de escorrentía superficial
producto de la precipitación, además hay diversos puntos donde desfogan las aguas del
Debido al aporte de agua por flujos de escorrentía rápidos, la acción erosiva en las laderas
del río causa que el material coluvial y aluvial alojado en sus laderas, se precipite al cauce
de desecho como ladrillos y bloques de cemento que modifican la rugosidad del canal.
Esta disección de la quebrada, se acentúa aguas abajo, desde la parte media en las
42
inmediaciones de la Escuela Nueva Laboratorio - Emma Gamboa, en este punto se
con profundidades promedio 15 m, una mínima de 2 m y una máxima de 18m (Figura 7).
Figura 7. Fondos de Valle (FV) de la Quebrada Negritos (QN). (A) Vista noreste (aguas arriba) de las
laderas del FV. (B) Vista hacia el oeste (aguas abajo) donde se aprecia material de coluvial (MC) y
aluvial (MA) del cauce principal de QN, la línea con triángulos naranja, detalla la ladera sur del río
que está sometido a procesos erosivos. Fuente: Trabajo de campo.
estudio, tienen una pendiente media de 5,52°. Estas unidades son de baja energía y
colindan con la base de los frentes de la lavina en las zonas ribereñas. En las partes más
43
encañonadas de las laderas de las TA, se pueden encontrar materiales de origen aluvial y
Dichas laderas están sometidas a la acción erosiva del río en las partes medias y bajas,
estas áreas se han utilizado para el desarrollo de infraestructura, parques y zonas verdes al
margen del río, en ellas se divisan los distintos niveles (4) de terrazas en la zona más baja
Figura 8. Terrazas Aluviales (TA). (A) Vista hacia el noreste de las TA de nivel 3 y 4 de la QN, sector
cuenca baja, Parque del Duende. (B) Vista aérea sobre las TA, en la cuenca baja de la QN, donde se
aprecia cobertura forestal y unidades de infraestructura sobre las TA de nivel T1 y T2. Fuente:
Trabajo de campo.
44
Laderas de Baja Pendiente (LBP) & Laderas de Pendiente Moderada (LPM):
Estas unidades conformadas por los lomos de los frentes de la lavina, son zonas casi planas
y anchas, poseen una pendiente media de 3,54°, de baja o poca disección ligeramente
ondulada y corresponden a un 42% del área de estudio, estas unidades con orientación E-
flujo tiende a ser de media, lo que favorece flujos superficiales rápidos en las vías y en el
alcantarillado pluvial. Hacia el sur del afluente principal, estas unidades colindan de
manera más abrupta con la QN, debido que el sector sur de la cuenca tiende a ser más alto
que el sector norte. (Ver Figura 12. Perfiles A-B, C-D, E-F y G-H).
Figura 9. Laderas de Baja Pendiente y Pendiente Moderada (LBP - LPM). (A) Vista hacia el norte
desde parte media alta de la cuenca de la QN, donde se aprecia el relieve de baja energía con ligeras
ondulaciones. (B) Vista al noreste desde la cuenca media, se aprecian LPM, elongadas y con dirección
45
N-S. Nótese una cobertura urbana tota, además de la orientación de la pendiente que promueve el
movimiento de la escorrentía superficial aguas abajo. Fuente: Trabajo de campo.
Las (LPM) Figura 9 B ,comprenden el 46% del área de estudio y tienen una pendiente
media de 4,53° y son las laderas paralelas al cauce principal de la QN, al igual que las
del flujo en esta unidad es media-alta, en ellas se captan y conducen las aguas provenientes
de las (LBP), estas áreas se encuentran en la parte alta de la cuenca y tienen una mayor
área de captación, aportando una mayor cantidad de agua al canal principal (Figura 9, A).
Estas unidades ocupan el 6% del área de estudio y responden a las zonas de mayor
pendiente (9° en promedio y máximas de 20°) de los frentes de la lavina. Dos sectores
son paralelos a los fondos de valle en la parte baja de la cuenca, en las inmediaciones del
elevación de este sector. Las otras tres unidades se encuentran en la parte media y alta de
la cuenca donde son visibles escarpes que corresponden a los límites de los flujos de la
lavina y cubiertos por capas de ceniza, en ellas se encuentran parques y áreas verdes de
46
Figura 10. Laderas de Pendiente Moderada y Escarpadas. (LPM-LSC). (A) Pequeño escarpe de unos
50 metros de longitud. (B) Vista hacia el oeste desde la parte alta de la QN, se muestran diferencias
de pendiente entre la ladera sur (izquierda) y la ladera norte (derecha). (C). Ladera de pendiente
moderada, urbanización San Marino, la diferencia de alturas entre las viviendas muestra que el
terreno tuvo que someterse al terraceo perpendicular a los frentes elongados de la lavina, la flecha
celeste indica la dirección de la ladera.
Red Hidrográfica:
sugieren la existencia de una red de drenaje más amplia con una extensión de 1872 m.
información extraída del MDS, hacia el extremo noreste de la cuenca, se infieren más
47
afluentes, sin embargo, en campo estos no se pueden comprobar debido al desarrollo
urbano. No obstante, el trazado vial se considera como afluentes artificiales que se activan
con la precipitación, además sistema de alcantarillado pluvial que transfiere esas aguas
cuadras y vías son en general, perpendiculares a los frentes de la lavina. Este patrón urbano
contribuye a que las vías, caños y red de drenaje tengan desfogues de agua perpendiculares
al cauce principal, al ser aguas pluviales estos flujos son rápidos e inciden en un aumento
Figura 11. Modificaciones de la red hidrográfica. (A) Caja de registro al pie de una (LPM),
urbanización San Marino. (B) Confinamiento del cauce principal de la QN, en las inmediaciones de
la Escuela Emma Gamboa. (C) Desfogue de alcantarillado pluvial so sobre la QN en el Colegio
Monterrey.
48
496020E
z .....
o o
~ ~
.....
~ ~
....
o
z
123B 1270
1236
1234
1260
1232
1230
122B 1250
1226
100 200 30(
1296
1294
1280
1292
1290
l270 [D 1288
1286
1284
Simbología
Puntos de muestreo Unidades Geomorfológicas. ■ cuenca
Punto de aforo ■ Fondo de valle
• Estación CIGER ~ Valle confinado
• desfogues Terrazas
Laderas de baja pendiente
Canales princi piles
- - Abierto ■ Laderas de pendiente mocerada
-confinado Laderas escarpadas
Figura 12. Distribución y localización de las unidades geomorfológicas, red drenaje y desfogues pluviales de la
QN.
49
Capítulo 3. Evaluación del régimen hídrico en términos de
cantidad de agua de diferentes fuentes geográficas de agua en
la sección media-alta de la microcuenca de Quebrada
Negritos, San Pedro de Montes de Oca, mediante técnicas de
hidrología isotópica.
dinámica del agua en la zona y su relación con factores climáticos y geográficos. Para
El año 2019 tuvo un comportamiento normal dentro del contexto climático del Valle
Central y la Vertiente del Pacífico. Las anomalías de precipitación respecto a los otros
años de la serie se vieron afectadas por eventos a escala sinóptica como la Zona de
Oscilación del Sur (ENOS) en su fase cálida (El Niño), fría (La Niña) y neutra. La
precipitación anual fue de 1625,1 mm, con un déficit del 7%. La distribución interanual
50
mm, para febrero y marzo de 0,1 mm. Durante abril se establece la transición a la
temporada lluviosa con un acumulado de 17,6 mm. En ese primer cuatrimestre del 2019,
la influencia de El Niño fue notoria, habiendo déficit de precipitación del 85%, 98%, 99%
vertiente del Caribe y en un alto déficit en la vertiente del Pacífico (IMN, 2019a, 2019b,
2019c).
Para abril, el patrón seco se mantuvo gobernado por El Niño y la Oscilación Madden-
En contraste, mayo fue el mes más lluvioso del 2019, con una precipitación acumulada de
(MJO) estuvo en fase divergente y frente a las costas del Pacífico de Costa Rica y
Nicaragua, se ubicó un sistema de baja presión, lo que permitió el ingreso de humedad del
Durante el mes de junio, se registró una precipitación acumulada de 208,1 mm, lo que
representó un déficit del 20% en comparación con la climatología de la serie. Este patrón
un total de 128,7 mm y un déficit del 32%. Este comportamiento fue similar al observado
51
en el Valle Central, donde los déficits de precipitación oscilaron entre el 11% y el 74% en
ambos meses, y se atribuyó a la fase neutra del ENOS (IMN, 2019f, 2019g). En agosto, la
precipitación fue de 194,2mm con un déficit del 6%. Al igual que en los meses anteriores,
Valle Central, presentándose déficit entre 16% y 28% y superávit de un 18% en distintas
acumuladas más altas del año, con 284,1 mm y 274,6 mm respectivamente, y un déficit
del 10% y 3% cada uno. Este patrón se repitió en otras estaciones del Valle Central, donde
el déficit en la precipitación fue de hasta el 60% (IMN, 2019i, 2019j). Durante noviembre
transición a la temporada seca (IMN, 2019k). Para diciembre se registró una precipitación
Central con déficits entre el 30% y 70% (IMN 2019l), ver figura 6.
52
40 400
30 300
-o
~ al
()
-5·
~
:::;¡
-¡¡; 20 200 ~
5:
ai
o.
:::;¡
E 3
~
-ª-
10 100
o o
Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Die
Mes
99 de la precipitación horaria, para cada mes de la serie, tomando en cuenta los meses de
estudio marzo-setiembre del 2019. Además, tomando el máximo de los eventos superiores
al percentil 99, demás, al analizar el máximo de los eventos que superaron dicho percentil
En mayo, se registró una precipitación máxima horaria de 33,8 mm, la cual superó
ampliamente el valor del percentil 99 para mayo, que era de 17,5 mm. Este hallazgo
característica de ser el mes más lluvioso del año (ver Figura 13).
En junio, se registró una precipitación máxima horaria de 31,5 mm, la cual también superó
significativamente el valor del percentil 99 para junio, que estaba establecido en 4,1 mm.
53
En julio, la precipitación máxima horaria alcanzó los 16,7 mm, superando con creces el
percentil 99 para julio, que era de 6,6 mm. Estos valores también clasifican como eventos
En agosto, se registró una precipitación máxima de 25,8 mm en una hora, la cual superó
horaria alcanzó los 40,5 mm, superando el percentil 99 de 8,0 mm para este mes. Este
La caracterización isotópica permite conocer los estadios del agua dentro del ciclo
la Línea Evaporativa Local (en adelante LEL) se tomaron muestras de FGA con 14
14,61‰. Esto indica que hay un empobrecimiento del trazador a medida que avanza la
época lluviosa. Si tomamos como referencia los primeros 6 muestreos entre el 20/05/2019
y los últimos 2 el 30/8/2019, vemos que los primeros 6 tienen valores entre -13,60‰ y -
10,64‰, mientras que los últimos se encuentran entre -7,04‰ y -3,18‰. Esto se relaciona
con la estacionalidad de la lluvia, donde las primeras muestras tienden a perder oxígeno
54
precipitación durante el mes de mayo favoreció procesos más intensos de evaporación y
mayor influencia de aguas oceánicas incorporadas por tormentas u ondas tropicales. Por
18
otra parte, el valor promedio de O‰ encontrado en el caudal es de -9,85‰, con un
y= 0.12 X+ -1.53 •
-30 R"2 = 0.9
y= 0.12 X+ -1.24 ♦
R"2 = 0.96 ♦
----
~
0
-60 ♦
I
N
IO
•
-90
•
•
-15 -12 -6 -3
Figura 14 . Línea evaporativa y meteórica local, Quebrada Negritos, 2019, muestra una correlación
casi idéntica entre los isótopos en la lluvia (puntos anaranjados) R2 de 0.96 y en el caudal R2 0.90
(puntos morados).
escorrentía, con muy poca mezcla con el agua previamente almacenada en la cuenca. Los
55
datos isotópicos indican que la conversión de la lluvia en escorrentía es veloz debido a las
períodos: antes del evento de lluvia (6), durante el evento (13) y posterior al evento (10).
a setiembre de 2019.
Los caudales aforados oscilaron entre los 0,07 m3/s y los 41,6 m3/s, con un promedio de
4,50 m3/s. En los muestreos antes de la precipitación, el caudal promedio es de 0,18 m3/s,
fluctuando entre los 0,07 m3/s a los 0,61 m3/s, con precipitaciones entre los 0,0 mm y los
1,30 mm.
El caudal promedio durante los eventos de precipitación fue de 9,03 m3/s, con un mínimo
de 0.07 m3/s y un máximo de 41,6 m3/s. Las precipitaciones promedio fueron de 11,39
mm, oscilando entre un mínimo de 1,60 mm y un máximo de 33,80 mm. En los muestreos
durante la precipitación, se identificaron los caudales más altos, 32,5 m3/s y 41,6 m3/s.
fracciones de agua nueva en promedio fueron del 66%, lo que indica que más de la mitad
56
3.3 Balance hidrológico simulado:
diciembre del 2019 donde se pueden distinguir diferentes FGA y sus porcentajes en
relación a la descarga total. Además, se realizó una estimación del balance hidrológico
para todo el año 2019, sin embargo, es importante señalar que no se consideraron los
meses de enero, febrero, marzo y los primeros 8 días de abril que corresponde a la fase
para inicializar el modelo. En general, el modelo llegó a una precisión alta durante la
calibración y el RMSE indica un ajuste de 5,53 E+07 y el MAE 2,87 E+08 (Figura 15).
40
~
E 30
.s
e
:g 20
.'!!
·a.
-~ 10
a::
o
01/01/19 01/04/19 01/07/19 01/10/19 01/01/20
~
E
E
-; 0.4
•O
·¡;
~
·a.
(/)
~ 02
oc.
>
w
"' o.o
01/01/19 01/04/19 01/07/19 01/10/19 01/01/20
8 ••
6
~
~ 4 • •
•
ª 2
o
01/01/19 01/04/19 01/07/19 01/10/19 01/01/20
Figura 15. Balance hidrológico diario simulado para el periodo de estudio 2019, Quebrada
Negritos, las barras azules indican los valores horarios de precipitación, las verde claro la
evapotranspiración horaria y verde oscuro el caudal simulado, los puntos corresponden a los
caudales observados .
57
Cabe mencionar que para poder cerrar el balance hidrológico del modelo se tuvieron que
de m3/s a mm/h. Esto debido a que los aforos realizados superaban tres veces el volumen
de lluvia caída.
QN estimadas en mm responden a los meses con más lluvias tienen descargas mayores y
contexto de la evaluación del modelo con respecto a los eventos extremos, se observa que
los valores de RMSE y MAE representan el error promedio entre las predicciones del
modelo y los valores observados. Estos valores numéricos indican que la precisión de las
correlación positiva entre las simulaciones del modelo y los datos observados, lo que
sugiere que el modelo sigue las tendencias generales de los eventos extremos. Esta
correlación positiva respalda la capacidad del modelo para capturar patrones similares a
los valores reales, así como la dinámica del hidrograma en relación a los eventos de
El caudal mensual presenta un promedio de 174,9 mm, siendo octubre, mayo, junio,
setiembre y agosto los meses con descargas superiores al promedio, con valores de 306,7
mm, 300,1 mm, 287,1 mm, 228,8 mm y 177,4 mm, respectivamente. En cambio, abril,
16 mm, 24,9 mm, 107,9 mm y 125,6 mm, respectivamente. En relación a las fracciones
aporta en promedio un 2% del total del caudal, con un máximo de 4% en mayo y valores
58
de 3% en los meses de agosto y setiembre, 2% en abril, junio, julio, octubre y noviembre
y un 0% en diciembre.
Por su parte, los aportes laterales de agua del suelo al caudal representan en promedio el
15% del total del caudal. Los aportes de agua al caudal, tanto subterráneos como
superficiales, presentan una distribución desigual a lo largo del año. Los aportes de agua
por escorrentía lateral tienen mayor proporción en mayo, con un 21%, seguidos de abril y
agosto con 17% cada uno, julio con un 16%, de octubre a noviembre con un 14%, junio
En cuanto a los aportes subterráneos, estos representan el mayor porcentaje del caudal,
con un promedio del 83%. Los mayores aportes se presentan en junio con un 86%, julio
con un 82%, octubre y noviembre con un 84% y diciembre con un 91%. Por otro lado, los
meses con menores aportes de agua subterránea son abril con un 81%, mayo con un 76%,
agosto con un 80% y setiembre con un 77%. Los datos muestran que los aportes de las
distintas FGA se correlacionan con los meses de mayor o menor precipitación, indicando
precipitación.
En el caso de los flujos de agua subterránea, se observa una ligera inversión en algunos
meses, como en abril, mayo, agosto y setiembre, estos son los meses con mayor
precipitación pero menos aportes de agua subterránea. Esto sugiere que los procesos de
aporte de agua subterránea al caudal total durante los meses de menor precipitación. Ver
figura 16.
59
01 /04/19 01/07/19 01/10/19 01 /01 /20
Mes
Figura 16. Separación del hidrograma según los aportes de la escorrentía superficial simulada, flujos
laterales simulados y flujos subterráneas al caudal simulado, Quebrada Negritos, 2019.
60
Tabla 2: Datos de muestro por evento. Elaboración propia con base en trabajo de campo y datos del (CIGEFI).
Fecha y hora Muestra P Q P_δO18 P_δ2H Q_δO18 Q_δ2H δv MQ_δO1 δn MP_δO18 FAN
(‰)
20/5/2019 10:00 1 14,70 0,148 -13,36 -93,98 -11,72 -82,00 -9,38 4 -13,36 1 0,59
20/5/2019 14:00 2 4,20 0,207 -13,61 -98,64 -13,41 -94,40 -11,72 1 -13,61 2 0,89
20/5/2019 15:00 3 0,20 0,177 -13,00 -94,58 -11,98 -84,68 -11,72 1 -13,00 3 0,20
21/5/2019 15:00 5 1,50 0,266 -10,64 -77,40 -9,48 -63,71 -9,38 4 -10,64 5 0,08
21/5/2019 17:00 6 1,60 0,354 -12,06 -83,65 -10,24 -72,16 -9,48 5 -12,06 6 0,29
24/5/2019 16:00 11 0,20 0,148 -8,15 -66,39 -9,50 -66,57 -9,50 11 -8,15 11 0,00
27/5/2019 15:00 13 12,10 5,664 -11,50 -78,29 -10,77 -72,51 -9,50 11 -11,50 13 0,64
28/5/2019 17:00 15 3,60 0,071 -9,31 -63,79 -9,13 -61,58 -11,18 14 -9,31 15 1,10*
12/6/2019 16:00 23 0,00 0,000 -14,61 -109,89 -13,98 -108,56 -10,79 22 -14,61 23 0,83
3/7/2019 12:00 24 4,10 1,062 -12,21 -104,04 -6,15 -50,15 -13,98 23 -12,21 24 NA
8/7/2019 16:00 25 16,70 2,124 -8,88 -61,95 -8,25 -58,58 -6,15 24 -8,88 25 0,77
21/8/2019 17:00 27 9,70 9,912 -5,59 -35,95 -7,10 -40,87 -10,22 26 -5,59 27 0,67
P: precipitación en milímetros por hora, Q: caudal en metros cúbicos por segundo, P_δO18(‰) concentración isotópica por mil de Oxígeno
18 en muestras de aguas de precipitación, P_δ2H (‰) concentración isotópica por mil de Deuterio en muestras de aguas de precipitación.
Q_δO18 (‰) concentración isotópica por mil de Oxígeno 18 en muestras de aguas de escorrentía y Q_δ2H(‰) concentración isotópica por mil
de Deuterio en muestras de aguas de escorrentía, δv: Isótopo de agua vieja, M Q_δO18(‰): muestra de referencia para el agua vieja, δn: Isótopo
+, .-)./
de agua nueva, MP_δO18(‰): Muestra de referencia para agua vieja, FAN: Fracción de Agua Nueva: +- = . * La muestra 15 fisicamente
.,)./
no es posible, solo matemáticamente.
61
Capítulo 4. Evaluación de influencia de las variables
hidrogeomorfológicas e isotópicas en el síndrome de la
corriente urbana en la microcuenca media-alta de Quebrada
Negritos.
Las observaciones realizadas durante 2019 permiten realizar una síntesis para evaluar la
modificación del paisaje que desencadena las manifestaciones del SCU en la QN.
En la Figura 17 (A) y B), correspondiente a la ortofoto del IGN de 1974 y 1989, se puede
Sociales, anteriormente era una zona permeable con pastos y árboles dispersos. Sin
al norte, hacia las terrazas de la QN, lo que propicia que la escorrentía superficial sea
usos de suelo son aún más drásticos, entre las fotografías aéreas de 1974 y 1989, se puede
observar cómo las áreas circundantes a la QN dejaron de ser zonas permeables, como
62
cafetales, pastos y bosques dispersos, para convertirse en una urbanización de alta
densidad. Esto incrementó la conectividad con el canal de la QN, debido al aumento del
Además, esta área se caracteriza por laderas de pendiente moderada con orientación sur,
y de baja pendiente con orientación norte, este efecto entre las laderas y su
un área permeable de pastos con un cauce tributario de la QN. Sin embargo, hacia 1989,
extiende de norte a sur sobre laderas de pendiente moderada que conducen y facilitan el
flujo de la escorrentía superficial y las aguas provenientes de la red pluvial, las cuales son
conducidas aguas abajo para ser descargadas en la QN. En este punto, el afluente fue
parques recreativos.
pastizales entre 1974 y 1989, consolidándose la condición urbana hasta el año 2019. En
hacia la QN por la red pluvial, y en el flanco sur no hay conexión directa con la QN, sino
que estas aguas son conducidas por la red de alcantarillado hacia aguas debajo de la QN.
63
En el punto 5, Urbanización Carmiol, la cobertura urbana estaba presente durante la
orientación sur de las laderas de baja y moderada pendiente colindan con las terrazas del
río. En el punto 6, en los alrededores del bar La Bamba, se presentan cambios de uso,
pasando de bosques y pastos a una zona urbana, en este sector las laderas de baja
pendiente, tiene un cambio abrupto hacia laderas escarpadas que paulatinamente vuelven
Simbología
- Cuenca Quebrada Negritos
Puntos de referencia
A: Ortofoto IGN 1974
B: Ortofoto IGN 1989
C: Imagen Planet 2019
1: Facultad de Ciencias Sociales
2: Escuela Nueva Laboratorio
3: Urbanización San Marino
4: Cementerio de Sabanilla
5: Urbanización Carmiol
6: La Bamba
o 500 1000 m
64
Como se evidencia, los cambios en el uso de suelo han modificado el paisaje,
incrementa los flujos rápidos y puede causar cambios como erosión lateral y vertical,
del canal del río, resulta en la aceleración del flujo debido a las variaciones en el área del
cauce, causando erosión diferencial entre la infraestructura y el material de las laderas, así
Cambios hidrológicos:
Al ser una metodología experimental, uno de los retos del muestreo y la toma de datos es
12 para el cálculo de la FAN: las muestras 1, 2, 3, 4, 5, 6, 11, 13, 15, 23, 28 y 29.
Los cambios hidrológicos generados por las transformaciones en el paisaje urbano tienen
65
tormentas. La frecuencia y magnitud de los eventos de precipitación afectan directamente
modelo de mezcla.
En las muestreas 1, 2 y 3 con fecha del 20/05/2019 entre las 12:00 hrs y las 15:00 hrs, el
evento de precipitación previa a la muestra 1 11:00 hrs es de 0,1 mm. Al momento del
nueva fue del 59%, en un caudal de 0,148m3/s. A las 13:00 hrs la precipitación fue de
0mm y para la muestra 2 a las 14:00 hrs, la precipitación fue de 4,20mm, con un aporte
de agua nueva del 89% y un caudal de 0.207m3/s. Además, durante la muestra 3 15:00 hrs
restableciéndose el flujo base y que el aporte de agua nueva en el evento 2, está regido
Para las muestras 4, 5 y 6 del 21/05/2019 entre las 10:00 hrs y las 17:00 hrs. Donde antes
9:00 hrs y durante el muestreo 4, a las 10:00 hrs, la precipitación fue de 0mm.
Consecuentemente, la fracción de agua nueva fue de 0%, con un caudal de 0,078m3/ que
precipitación. Para el muestreo 5, 15:00 hrs la precipitación fue de 0,3mm, con un caudal
de 0,266m3/s y el aporte de agua nueva al caudal fue del 8%. Con la muestra 6, la
66
precipitación fue de 1,60mm, lo que incrementó el caudal a 0,354m3/s, así mismo el aporte
Durante la muestra 11 del día 24/05/2019 a las 16:00 hrs, la precipitación fue de 0,20 mm
con un caudal de 0,14 m3/s. Este volumen pertenece al flujo base ya que no hay aportes
de aguas nuevas, el 0% de agua nueva se explica debido a que en las tres horas previas al
muestreo 11 la precipitación no superó los 0,80 mm. Esto responde a que la intensidad y
Para el 27/05/2019 a las 15:00 hrs, muestra 13, correspondió a una precipitación de 12,10
mm, con un caudal de 5,66 m3/s. El 64% del pertenecen a aportes de agua nueva, en este
17:00 hrs, en las 5 horas previas al muestreo, la precipitación acumulada fue de 1,6mm,
sin superar los 0,5 mm/h. Al momento del muestreo 15, la precipitación fue de 3,60 mm,
con un caudal de 0,71m3/s, donde el 100%, del agua en el caudal pertenece a agua nueva,
sin embargo, este dato puede pertenecer a un error en la medición o en las características
67
Durante las 4 horas previas a la muestra 23 del 12/06/2019 a las 16:00 hrs, la precipitación
acumulada fue de 22 mm con un evento principal de 14,2 mm a las 13:00 hrs. Durante el
muestreo fue de 0mm, pero el aporte de agua nueva al caudal fue de 83%. Es decir, la
16:00 hrs, tomada durante el evento de lluvia, cuya precipitación fue de 16,70 mm, con
un caudal de 2,12m3/s. El aporte de agua nueva fue de un 77%, ya que en las tres horas
previas al muestreo la precipitación no superó los 1,9 mm, esto muestra que en la mezcla
La muestra 27 del 21/8/2019 a las 17:00 hrs, la precipitación fue de 9,70 mm, con un
El día 30/8/2019 se tomaron las muestras 28 13:00 hrs y 29 15:00 hrs, en las 3 horas
previas al muestreo 28, la precipitación no fue mayor a los 0,1 mm, con un evento principal
de 6,90 mm y los aportes de agua nueva fueron de 53% y 25%. La cantidad de agua se
relaciona al evento de lluvia principal, donde los aportes de agua nueva dominan en la
mezcla de la corriente, en la muestra 28 el caudal fue de 8,14m3/s. Sin embargo, dos horas
retención y regulación del paisaje por los cambios de uso y cobertura del suelo. El
68
volumen de agua nueva presente en la mezcla del caudal se asocia directamente a la
10.0
7.5
5.0
-¡¡¡
"O
:::,
ü"'
2.5
o.o
o
---
FAN
--Eventos muestreados
o ■ 25 ■ 50 ■ 75
10
■ 100
- 15
Figura 18. Fracciones de agua nueva (FAN) respecto al caudal muestreado las barras muestran en
color azul los eventos muestreados donde la FAN es superior al 50% hasta el100%, en amarillo los
eventos con valores menores al 50% de FAN en el caudal..
69
Capítulo 5. Discusión de resultados, conclusiones y
recomendaciones.
5.1 Discusión de resultados:
en función a las variables estudiadas, y cuales fueron identificadas en la QN. Los síntomas
Tabla 3. Síntomas del SCU identificados en la Quebrada Negritos según Walsh et al.
(2005)
superficial. superficial.
70
Geomorfología + Anchura del canal + Anchura del canal
urbanos como la QN, aunque con elevados porcentajes escorrentía superficial proveniente
de la lluvia (hasta un 89%) y conducida por la red pluvial hacia el canal, aún hay
abajo.
71
Un efecto relacionado al aumento de los volúmenes del caudal que requiere atención es la
podrían verse potencialmente afectadas por la socavación lateral del río en los fondos de
alta no hay mecanismos suficientes para la gestión de aguas pluviales, combinada con el
Esto ha propiciado aumentos en los niveles del caudal de hasta 32,45 m3/s y 41,6 m3/s. La
cauce principal del río, ha generado una disminución en el tiempo de retardo del flujo
máximo, así mismo el río vuelve a su nivel base rápidamente, lo que muestra alteraciones
zonas vulnerables.
Según (Fletcher et al., 2011b) en áreas urbanas la escorrentía se activa con precipitaciones
72
superiores a 1 mm, siempre y cuando estén conectadas directamente con el canal principal,
condición urbana imperante desde 1989. Esto afecta directamente la estabilidad del canal
el 20% del total de la cuenca, la escorrentía se duplica. Si el área impermeable oscila entre
el 100%, la escorrentía superficial se multiplica por más de cinco en comparación con las
el aporte de la escorrentía superficial a la mezcla del caudal. Ya que las FAN indican que
lluvia-escorrentía. Además, los tiempos de residencia del agua son relativamente cortos,
esto es evidente en los cambios de fracción entre eventos muestreados el mismo día, como
el caso de las muestras 2 y 3 donde en 1 hora la fracción de agua nueva osciló entre 89%
nueva pasó de 8% a 29%. Esto indica lo dinámico del sistema en cuento a la incorporación
73
de agua de lluvia en el caudal y la conducción y evacuación del agua de la escorrentía
de lluvia diaria en la serie 2006-2019 los eventos que superan el percentil 99, con valores
aumentos en la magnitud del flujo pico, atenuado por las aguas provenientes de las aguas
dentro de un plan regulador, estas áreas tienen una importancia fundamental para la
ser una prioridad para las autoridades encargadas de la gestión territorial y la planificación
urbana.
En relación a las Zonas de Protección y el régimen de propiedad, debe ser entendido como
cuanto a la zona de protección del río y las zonificaciones propuestas dentro de un plan
74
regulador puede ser un factor estructural que contribuye a la manifestación del SCU. Es
necesario fortalecer las figuras institucionales que regulan los bienes de uso público y la
5.2 Conclusiones.
QN, así como el desconocimiento sobre distribución, localización y tamaño de red pluvial
del cantón de Montes de Oca, impide una mejor gestión de estas, además del continuo
75
redujo las áreas permeables y de almacenamiento de agua en el suelo, lo que limita su
la mezcla del caudal. Los tiempos de residencia del agua son relativamente cortos, lo que
alcantarillado pluvial.
aguas pluviales, para mejorar la gestión territorial y la planificación urbana del cantón de
Montes de Oca.
embargo presenta limitaciones a nivel metodológico, por lo que se deben integrar diversos
76
implementado para el fraccionamiento de agua, probó su utilidad para estimar cómo se
hidrológico, esto sugiere el aporte de agua de otras fuentes por medio trasvases o red
pluvial de otros sectores que desfogan en la QN. Así mismo, la incorporación del análisis
infraestructura condicionado por las características del terreno, por lo que provee un
del río limita la capacidad de conservar las áreas verdes cercanas a la QN. Esto es resultado
Montes de Oca y las comunidades que podrían verse afectadas. Un enfoque en la gestión
de las aguas pluviales mediante obras de retención y la recuperación de áreas verdes, junto
77
con la planificación urbana, podría constituir esfuerzos dirigidos a mitigar el riesgo y
la QN.
5.3 Recomendaciones:
1. Establecer una red de monitoreo y otros puntos de muestreo tanto de caudal como
2. El monitoreo del SCU, podría ser un insumo para diversas instituciones o actores
en zonas urbanas y tomar acciones para disminuir o revertir los efectos negativos
zonas con mayor afectación, donde localizar infraestructura verde y sitios que se
urbanas.
3. Generar un modelo hidrológico conceptual calibrado con los datos isotópicos para
estimar las fracciones de agua nueva en el río y simular los aportes de la lluvia a
la descarga total.
78
4. Ampliar las áreas de estudio, para comprobar si la metodología es replicable en
heterogéneas.
QN, con ello desarrollar y colocar infraestructura verde que compense los efectos
escorrentía rápida.
que puede ser eficaz para reducir los efectos del cambio hidrológico. Este modelo
para reducir los riesgos de inundaciones y mejorar la calidad del agua en la QN.
agua.
79
Referencias
Bhaskar, A. S., Beesley, L., Burns, M. J., Fletcher, T. D., Hamel, P., Oldham, C. E., &
Roy, A. H. (2016). Will it rise or will it fall? Managing the complex effects of
urbanization on base flow. Freshwater Science, 35(1), 293–310.
https://doi.org/10.1086/685084
Bhaskar, A. S., Hogan, D. M., Nimmo, J. R., & Perkins, K. S. (2018). Groundwater
recharge amidst focused stormwater infiltration. Hydrological Processes, 32(13),
2058–2068. https://doi.org/10.1002/hyp.13137
Birkel, C., Correa Barahona, A., Duvert, C., Granados Bolaños, S., Chavarría Palma, A.,
Durán Quesada, A. M., Sánchez Murillo, R., & Biester, H. (2021). End member
and Bayesian mixing models consistently indicate near-surface flowpath
dominance in a pristine humid tropical rainforest. Hydrological Processes, 35(4).
https://doi.org/10.1002/hyp.14153
Booth, D. B., Roy, A. H., Smith, B., & Capps, K. A. (2016). Global perspectives on the
urban stream syndrome. Freshwater Science, 35(1), 412–420.
https://doi.org/10.1086/684940
Camacho Lizano, D., Chaves Gamboa, J., & Murillo Montoya, D. (2004). Diagnóstico
de la susceptibilidad al deslizamiento e inundación, análisis del contexto sísmico y
reconocimiento de peligros volcánicos para el cantón de Montes de Oca: Una
contribución al plan regulador del municipio. Geología, 1(1), 55.
Capps, K. A., Bentsen, C. N., & Ramírez, A. (2016). Poverty, urbanization, and
environmental degradation: Urban streams in the developing world. Freshwater
Science, 35(1), 429–435. https://doi.org/10.1086/684945
80
Christophersen, N., Neal, C., Hooper, R. P., D., V. R., & Sjur Andersen. (1990).
MODELLING STREAMWATER CHEMISTRY AS A MIXTURE OF SOILWATER
END-MEMBERS -- A STEP TOWARDS SECOND- GENERATION
ACIDIFICATION MODELS. 116, 307–320.
CNE. (1990). Inundaciones en las zonas adyacentes a la Quebrada Los Negritos, 1989.
https://www.cne.go.cr/CEDO-CRID/CEDO-
CRID%20V4/pdf/spa/dec88/dec88.htm
Correa, A., Breuer, L., Crespo, P., Célleri, R., Feyen, J., Birkel, C., Silva, C., &
Windhorst, D. (2019). Spatially distributed hydro-chemical data with temporally
high-resolution is needed to adequately assess the hydrological functioning of
headwater catchments. Science of the Total Environment, 651, 1613–1626.
https://doi.org/10.1016/j.scitotenv.2018.09.189
Daniel Poleo Brito, Eladio Solano León, & Werner Stolz España. (2014). La Oscilación
atmosférica Madden-Julian (MJO) y las lluvias en Costa Rica. 5–19.
http://cglobal.imn.ac.cr/documentos/revista/topicosmet20142/index.html
Dunn, S. M., & Mackay, R. (1995). Spatial variation in evapotranspiration and the
influence of land use on catchment hydrology. En ELSEVIER Journal of Hydrology
(Vol. 171).
Fletcher, T. D., Walsh, C. J., Bos, D., Nemes, V., RossRakesh, S., Prosser, T., Hatt, B.,
& Birch, R. (2011a). Restoration of stormwater retention capacity at the allotment-
scale through a novel economic instrument. Water Science and Technology, 64(2),
494–502. https://doi.org/10.2166/wst.2011.184
81
Fletcher, T. D., Walsh, C. J., Bos, D., Nemes, V., RossRakesh, S., Prosser, T., Hatt, B.,
& Birch, R. (2011b). Restoration of stormwater retention capacity at the allotment-
scale through a novel economic instrument. Water Science and Technology, 64(2),
494–502. https://doi.org/10.2166/wst.2011.184
Francis, R. A. (2012). Positioning urban rivers within urban ecology. Urban Ecosystems,
15(2), 285–291. https://doi.org/10.1007/s11252-012-0227-6
Gervasi, A. (2012). Land cover , urbanization , and slope effects on tropical stream
hydrology : the need for “ buff - er ” buffer zones. July 2012, 0–13.
https://doi.org/10.13140/RG.2.2.10154.31687
Glenn, D. T., Endter-Wada, J., Kjelgren, R., & Neale, C. M. U. (2015). Tools for
evaluating and monitoring effectiveness of urban landscape water conservation
interventions and programs. Landscape and Urban Planning, 139, 82–93.
https://doi.org/10.1016/j.landurbplan.2015.03.002
Groffman, P. M., Bain, D. J., Band, L. E., Belt, K. T., Brush, G. S., Grove, J. M.,
Pouyat, R. V., Yesilonis, I. C., & Zipperer, W. C. (2003). Down by the Riverside:
Urban Riparian Ecology. Frontiers in Ecology and the Environment, 1(6), 315.
https://doi.org/10.2307/3868092
Gurnell, A., Lee, M., & Souch, C. (2007). Urban Rivers: Hydrology, Geomorphology,
Ecology and Opportunities for Change. Geography Compass, 1(5), 1118–1137.
https://doi.org/10.1111/j.1749-8198.2007.00058.x
Hale, R. L., Scoggins, M., Smucker, N. J., & Suchy, A. (2016). Effects of climate on the
expression of the urban stream syndrome. Freshwater Science, 35(1), 421–428.
https://doi.org/10.1086/684594
82
Hidalgo, P. J., Alvarado, G. E., & Linkimer, L. (2011). La Lavina del Valle Central
(Costa Rica): ¿lahar o Debris Avalanche? Revista Geológica de América Central,
30. https://doi.org/10.15517/rgac.v0i30.7261
Hooper, R. P. (2003). Diagnostic tools for mixing models of stream water chemistry.
Water Resources Research, 39(3), 1–13. https://doi.org/10.1029/2002WR001528
Instituto Geológico y Mienro de España. (2002). Metodología isotópica 2.1 (Vol. 1, pp.
181–185).
83
Instituto Meteorológico Nacional. (2019g). Boletín Meteorológico Mensual, Julio 2019.
https://www.imn.ac.cr/documents/10179/474355/JULIO
Instituto Meteorológico Nacional (IMN). (2021). Región Central: Dos Valles (p. 4).
https://www.imn.ac.cr/documents/10179/31165/Clima+Valle+Central/9b5f5941-
3708-4d18-991f-ada9cf795864
Jefferson, A. (2021). Calculate the amount of new event water in the stream at each time
point using isotope and discharge data.
https://serc.carleton.edu/hydromodules/steps/153843.html
Kendall, C., McDonnell, J. J., & Gu, W. (2001). A look inside “black box” hydrograph
separation models: A study at the hydrohill catchment. Hydrological Processes,
15(10), 1877–1902. https://doi.org/10.1002/hyp.245
Klaus, J., & McDonnell, J. J. (2013). Hydrograph separation using stable isotopes:
Review and evaluation. En Journal of Hydrology (Vol. 505, pp. 47–64).
https://doi.org/10.1016/j.jhydrol.2013.09.006
84
Konrad, C. P., & Booth, D. B. (2005). Hydrologic changes in urban streams and their
ecological significance. American Fisheries Society Symposium, 2005(47), 157–
177.
McDonnell, J. J., McGuire, K., Aggarwal, P., Beven, K. J., Biondi, D., Destouni, G.,
Dunn, S., James, A., Kirchner, J., Kraft, P., Lyon, S., Maloszewski, P., Newman,
B., Pfister, L., Rinaldo, A., Rodhe, A., Sayama, T., Seibert, J., Solomon, K., …
Wrede, S. (2010). How old is streamwater? Open questions in catchment transit
time conceptualization, modelling and analysis. Hydrological Processes, 24(12),
1745–1754. https://doi.org/10.1002/hyp.7796
Mora, R., González, L. A., Fabián, G., Montero, C., Rolando, M. S., Chinchilla, M.,
Fabián, B., & Montero, C. (2015). DESLIZAMIENTOS EN CONDICONES
ESTÁTICAS Y DINÁMICAS EN LOS ALREDEDORES DE LA QUEBRADA LOS
NEGRITOS EN LA CUIDAD UNIVERSITARIA RODRIGO FACIO , MONTES DE
OCA , SAN JOSÉ. February 2020, 0–33.
Ogle, K., Tucker, C., & Cable, J. M. (2014). Beyond simple linear mixing models:
Process-based isotope partitioning of ecological processes. Ecological Applications,
24(1), 181–195. https://doi.org/10.1890/12-1970.1
Otto, J. C., Prasicek, G., Blöthe, J., & Schrott, L. (2018). GIS Applications in
Geomorphology. Comprehensive Geographic Information Systems, 3(2000), 81–
111. https://doi.org/10.1016/B978-0-12-409548-9.10029-6
85
Pataki, D. E., Boone, C. G., Hogue, T. S., Jenerette, G. D., Mcfadden, J. P., & Pincetl, S.
(2011). Socio-ecohydrology and the urban water challenge. Ecohydrology, 4(2),
341–347. https://doi.org/10.1002/eco.209
Soulsby, C., Birkel, C., & Doerthe, T. (2014). Assessing urbanization impacts on
catchment transit times. Geophysical Prospecting, April, 6298–6305.
https://doi.org/10.1002/2013GL058716.
Vietz, G. J., Rutherfurd, I. D., Fletcher, T. D., & Walsh, C. J. (2016). Thinking outside
the channel: Challenges and opportunities for protection and restoration of stream
morphology in urbanizing catchments. Landscape and Urban Planning, 145, 34–
44. https://doi.org/10.1016/j.landurbplan.2015.09.004
Vietz, G. J., Walsh, C. J., & Fletcher, T. D. (2015). Urban hydrogeomorphology and the
urban stream syndrome: Treating the symptoms and causes of geomorphic change.
Progress in Physical Geography, 40(3), 480–492.
https://doi.org/10.1177/0309133315605048
Walsh, C. J., Roy, A. H., Feminella, J. W., Cottingham, P. D., Groffman, P. M., &
Morgan, R. P. (2005). The urban stream syndrome: Current knowledge and the
search for a cure. Journal of the North American Benthological Society, 24(3), 706–
723. https://doi.org/10.1899/04-028.1
Wenger, S. J., Roy, A. H., Jackson, C. R., Bernhardt, E. S., Carter, T. L., Filoso, S.,
Gibson, C. A., Hession, W. C., Kaushal, S. S., Martí, E., Meyer, J. L., Palmer, M.
A., Paul, M. J., Purcell, A. H., Ramírez, A., Rosemond, A. D., Kate, A., Sudduth,
E. B., Walsh, C. J., … Walsh, C. J. (2009). Twenty-six key research questions in
urban stream ecology : an assessment of the state of the science Twenty-six key
research questions in urban stream ecology : an assessment of the state of the
science. https://doi.org/10.1899/08-186.1
Zalewski, M., & Wagner, I. (2005). Ecohydrology - The use of water and ecosystem
processes for healthy urban environments. Ecohydrology and Hydrobiology, 5(4),
263–268.
87