Revista Geográfica ISSN: 1011-484X Número 62 • Enero-junio 2019 Doi: http://dx.doi.org/10.15359/rgac.62-1.3 Páginas de la 81 a la 103 del documento impreso Recibido: 12/01/2018 • Aceptado: 8/6/2018 URL: www.revistas.una.ac.cr/index.php/geografica/

---

Sensibilidad y escenarios de disponibilidad de agua para consumo humano en la microcuenca del río Porrosatí, Heredia, Costa Rica

Sensitivity and scenarios of availability of water for human consumption in the micro watershed of Porrosati River, Heredia, Costa Rica

Esteban Montero-Sánchez1

Universidad de Wageningen, Holanda

Jorge Herrera-Murillo2

Pablo Ramírez-Granados3

Universidad Nacional, Costa Rica

Resumen La variabilidad climática y el cambio del uso de la tierra representan una de las mayores condicionantes para la disponibilidad de agua para consumo humano en el futuro. Con el objetivo de analizar la disponibilidad de agua para uso y consumo humano en la microcuenca del río Porrosatí, se calcularon balances hídricos con base en la recarga de agua subterránea. Para las estimaciones se utilizaron registros meteorológicos de 15 años, en el periodo comprendido entre el año 2000 a 2014, y análisis de las características hidrológicas. Con el fin de realizar escenarios de recarga hídrica en el mediano (2026-2030) y largo plazo (2051-2055), tomando en consideración escenarios y por el cambio climático se utilizaron las predicciones de variación de temperatura y precipitación elaborados por el Instituto Meteorológico Nacional, bajo un escenario de emisiones A2. Los resultados obtenidos muestran comportamientos atípicos con cambios significativos en el comportamiento estacional de la recarga. Los volúmenes proyectados a futuro dan cuenta de importantes disminuciones, que en algunos casos superan el 15% como media llegando hasta diferencias del 50%. En el escenario de largo plazo se notó un importante cambio en comportamiento estacional, dando como resultado la aparición de dos picos de recarga. Las proyecciones climáticas presentan un escenario complejo con limitaciones a la recarga hídrica de agua subterránea y, por ende, a la disponibilidad de agua para consumo humano en el mediano y largo plazo. Palabras clave: disponibilidad de agua; variabilidad climática; cambio climático; balance hídrico; modelaje. Abstract Climate variability and land-use change represent two of the major constraints to the availability of water for human consumption in the future. With the aim of analyzing the availability of water for human use and consumption in the Porrosati Microbasin, water balances were calculated based on the groundwater recharge. For estimates, meteorological records of 15 years (from 2000 to 2014) and analysis of the hydrological characteristics were used. To set up scenarios of water recharge in the medium (2026-2030) and long-term (2051-2055), taking into account scenarios of climate change, were used predictions of temperature and precipitation variation, developed by the National Meteorological Institute, under a scenario of emissions A2. The results obtained show atypical behaviors with significant changes in the seasonal behavior of recharge. Prospective volumes for the future show significant decreases that, in some cases, exceed the 15% on average, reaching up to 50% differences. In the long-term scenario, a major change is noticed in seasonal behavior, resulting in the emergence of two peaks of recharge. The climate projections present a complex scenario with limitations to the water recharge of groundwater and, hence, the availability of water for human consumption in the medium and long-term. Keywords: water availability; climate variability; climate change; water balance; modeling.

Introducción

La variabilidad climática, potenciada como consecuencia del calentamiento global, está generando importantes impactos sociales, ecológicos y económicos a nivel mundial (van der Velde et al., 2012 y Smith y Katz, 2013). El aumento de la temperatura, la variabilidad de las precipitaciones, el incremento en la intensidad y la frecuencia de fenómenos hidrometeorológicos extremos (inundaciones, sequías), afectan a los sistemas de agua dulce y su gestión (Bates et al., 2008) debido a la intrínseca relación existente entre el sistema climático y el ciclo hidrológico (Zhou et al. 2010). Así, el incremento de la temperatura promedio modifica los rangos de evaporación del agua en el suelo y las superficies acuáticas (Kusangaya et al 2014). En el caso de las superficies terrestres, la temperatura tiene efectos directos sobre los valores de evapotranspiración, siendo este factor a su vez, una limitante de la cantidad de agua en el subsuelo y, por ende, de la recarga de agua subterránea (Schosinsky & Losilla 2000). A su vez, las variaciones en la distribución e intensidad de los eventos de lluvia pueden generar inestabilidad respecto a la planificación de uso del recurso (Ureña, 2005). Los eventos de mucha intensidad, pese al gran volumen de agua caída que representan, no benefician necesariamente la disponibilidad de agua para consumo humano. Una mayor intensidad hace que el suelo no tenga la capacidad suficiente para percolar el líquido a profundidades mayores por lo que su saturación incrementa el porcentaje de agua por escorrentía (Schosinsky 2006).

La escasez de agua y los efectos de fenómenos extremos plantean un reto para todos los países, ya que exacerban la vulnerabilidad de los recursos hídricos socavando las posibilidades reales de desarrollo económico (Kusangaya et al., 2014). Si a los efectos derivados del cambio climático se le suman desafíos históricos que ha enfrentado la gestión integral del recurso hídrico, derivados entre otros, de la impermeabilización de zonas de recarga, la ausencia de planificación territorial y el crecimiento demográfico, hace que se potencie la incertidumbre y complejidad respecto a los escenarios de disponibilidad del agua en el corto y mediano plazo (Ureña, 2005).

En la región centroamericana se pronostican variaciones climáticas de severa magnitud como cambios y anomalías estacionales y aumento en la magnitud y frecuencia de eventos extremos como sequías e inundaciones (Christensen et al., 2007). Tan sólo para el Valle Central de Costa Rica se pronostica una reducción de entre -15 a -35 % en el nivel de precipitación con respecto a los datos actuales (Alvarado et al., 2012). Este comportamiento puede ser similar al presentado en épocas del fenómeno El Niño/ Oscilación del Sur (ENOS), en donde los eventos de sequía o precipitación suceden de manera extrema.

La generación de escenarios es una herramienta vital para la toma de decisiones y acciones sobre el manejo, uso y preservación del agua. Pese a la incertidumbre que implica realizar escenarios de variables como el clima o el uso del suelo en el futuro, los esfuerzos en esta dirección son valiosos en cuanto se van perfeccionando las técnicas y manejo de datos. En este sentido, el análisis de la disponibilidad de agua, utilizando la cuenca hidrográfica como unidad de análisis permite una visión integral de todos los datos que interaccionan con relación a la recarga de agua subterránea y la disponibilidad de agua para la población inmersa dentro de esta área (Zhou et al., 2010).

En el presente artículo se muestran los resultados de los escenarios de disponibilidad hídrica obtenidos con datos climáticos y de uso del suelo, basados en el análisis de las dinámicas de cambio hidrológico en los años 2000-2014 y la realización de proyecciones para el mediano y largo plazo en la microcuenca en estudio. Los resultados derivados del análisis realizado están dirigidos a los diferentes actores involucrados en la gestión del agua para uso y consumo humano en la microcuenca y pretende ofrecer un respaldo científico para apoyar el establecimiento de un plan de acción que garantice la protección del recurso hídrico a partir del desarrollo de medidas de adaptación al cambio climático en la microcuenca del río Porrosatí.

Metodología

Área de estudio

La microcuenca del río Porrosatí se encuentra ubicada en el Valle Central y cuenta con una gran cantidad de afloramientos naturales que abastecen de forma directa a cerca de 50 mil personas y probablemente a una cifra mayor de manera indirecta (Sibaja 2013). El punto más alto se encuentra en las faldas del volcán Barva, siendo los 2870 m sobre el nivel del mar el punto de mayor altitud. En las partes medias y bajas, la altitud desciende hasta cerca de los 500 m sobre el nivel del mar. De acuerdo con la figura 1, las partes altas presentan un relieve más accidentado, dando paso a terrenos más planos en donde se asientan importantes centros de población. La mayor área de la microcuenca, cerca del 80 %, se encuentra en un rango de pendientes del 8 a 30 %, lo que da cuenta de un terreno escarpado, principalmente, en las partes altas y medias junto a las zonas de cañón del cauce principal. Un 15 % se encuentra dentro del rango de 3 a 8 %, lo que da cuenta de terrenos planos ondulados hacía las partes medias y bajas.

La edad del material rocoso data del Cuaternario y se caracteriza por facies proximales de rocas volcánicas recientes, compuestas por coladas de lava, aglomerados, lahar y ceniza volcánica (USGS, 1997). Dentro de la microcuenca se encuentran formaciones hidrogeológicas del miembro Porrosatí-Carbonal, los cuales están formados por arenas volcánicas gruesas y tobas arcillosas meteorizadas, formando acuitardos de gran extensión, donde subyacen los acuíferos locales Barva Superior (SENARA, 1985). Estas pertenecen a la Formación Barva, constituida por coladas de lavas andesiticas y andesito basalticas, los cuales, en general, por fracturación de la roca favorecen los procesos de infiltración y percolación de agua hacia los acuíferos (Denyer & Kussmaul, 2000).

Figura. 1. Modelos de elevación digital de la microcuenca del río Porrosatí.

Fuente: Elaboración propia.

La mayor parte de la microcuenca presenta geomorfologías asociadas al volcán Barva y se le clasifica como un escudo andesíticas o estratovolcán (Denyer & Kussmaul, 2000). En cuanto a las partes medias y bajas de la microcuenca, se caracterizan por un piemonte de relieve ondulado a plano-ondulado, conformado por depósitos de lavas andesíticas del Cuaternario con capas de cenizas y piroclástos de origen lahárico (Mata & Ramírez, 1999; Bergoeing, 2007). La totalidad de la microcuenca está compuesta por suelos de tipo andisoles.

La microcuenca del río Porrosatí destaca por una gran cantidad de afloramientos de agua subterránea y un cauce superficial permanente en todas las épocas del año. Pese a su relativa poca extensión permite el abastecimiento de importantes de actividades humanas de la provincia de Heredia, en donde se encuentran al menos 161 concesiones entre manantiales, pozos y captaciones superficiales según datos de la Dirección de Aguas del MINAE (2015). Dentro de los entes que se abastecen de fuentes de agua dentro de la microcuenca para suministro humano se encuentran ASADAS, Acueductos Municipales, la Empresa de Servicios Públicos de Heredia y múltiples usuarios de carácter privado. Los caudales asignados van desde los 6 hasta los 30 L/s, siendo común encontrar múltiples concesiones de bajo caudal para un mismo ente operador.

Análisis de datos meteorológicos

Se identificaron un total de cuatro estaciones meteorológicas en las inmediaciones de la microcuenca, ubicadas en: Santa Bárbara, Aeropuerto Internacional Juan Santamaría, Santa Lucía, pertenecientes al Instituto Meteorológico Nacional (IMN) y la estación Monte de la Cruz perteneciente al Servicio Nacional de Aguas Subterráneas Riego y Avenamiento (SENARA). A partir de esta selección, se consolido una base de datos con 15 años de información meteorológica a escala diaria, la cual fue ordenada para identificar cualquier error o ausencia en los datos. Tomando en cuenta que las bases de datos presentaban datos faltantes que varían de un día a un mes completo, se procedió a realizar la identificación de vacíos en los datos para su rellenado y corrección.

Para la estimación y corrección de datos meteorológicos se utilizó el método de razón de valores normales (Guevara 2003), el cual emplea el promedio de al menos 3 estaciones con condiciones fisiográficas y climáticas (ecuación 1) que se consideren representativas de la estación por estimar. Cada valor es corregido por un factor basado en la relación de comportamiento entre la estación por estimar y la estación de referencia (Alfaro & Pacheco, 2000). (1)

Donde

p0 = valor de precipitación mensual de la estación 0 (mm).

G = número de estaciones utilizadas en la estimación de p0 (mm)

P0 = precipitación anual de la estación 0 (mm).

Pg = precipitación anual de la estación g

Los criterios de elegibilidad para la aplicación del método de las razones normales para la estimación de registros inexistentes o incompletos se describen en la Tabla 2.

Tabla 2. Criterios para la escogencia de estaciones para interpolaciones de datos meteorológicos.

Criterio	Justificación	Método utilizado
Distancia	La menor cercanía entre estaciones representa la condición deseable, pues se asume que las variables que afectan el comportamiento climático tendrán mayor similitud y, por ende, un comportamiento más homogéneo entre las estaciones.	Despliegue de las estaciones en un SIG para visualizar las distancias entre cada una de ellas.
Altitud	La altitud es una variable determinante, principalmente asociada a la conformación orográfica.	Mediante el SIG, se despliegan los datos de altitud de las estaciones. Se realizó un modelo de elevación digital, en donde se obtienen rangos de distribución de la información climática sobre el área de estudio.
Patrones de comportamiento	Cuando no se cumplieron los criterios anteriores, se procedió a tomar en cuenta variables basadas en la observación y análisis del comportamiento entre estaciones.	Además del cálculo de promedios de las sumatorias, se efectuaron análisis por mínimos cuadrados para conocer la relación de comportamiento entre las estaciones.

Fuente: Elaboración propia.

Determinación del Balance Hídrico para el período 2010-2014

Para el cálculo de la recarga potencial de acuíferos, se utilizó el modelo propuesto por Schosinsky (2006), el cual fue declarado oficial para la estimación de cálculos hidrogeológicos en Costa Rica, según el Acuerdo MINAET 60- 2012. Este modelo es una combinación entre los métodos de precipitación que infiltra y balance de humedad de suelos.

Tabla 3. Características hidrológicas en las zonas de balance para la microcuenca del río Porrosatí

Zona de balance	Características
Zona 1	Zona alta de la cuenca. Uso: bosque. Pendiente: mayor a 7 %. Estación meteorológica: Monte de la Cruz.
Zona 2	Zona alta de la cuenca. Uso: plantaciones de ciprés y pasto. Pendiente: mayor a 7 %. Estación meteorológica: Monte de la Cruz.
Zona 3	Zona media de la cuenca. Uso: pasto y bosque. Pendiente: menor a 7 %. Estación meteorológica: Santa Lucía.
Zona 4	Zona media de la cuenca. Uso: cultivo. Pendiente: menor a 7 %. Estación meteorológica: Santa Lucía.
Zona 5	Zona media y baja de la cuenca. Uso: cultivo. Pendiente: menor a 7 %. Estación meteorológica: Santa Bárbara.
Zona 6	Zona media y baja de la cuenca. Uso: cultivo y bosque. Pendiente: mayor a 7 %. Estación meteorológica: Santa Bárbara.

Fuente: Elaboración propia

Para la construcción de las zonas de balance hídrico, se analizó la información contenida en los mapas de pendientes reclasificadas, según los valores determinados por el modelo; el uso del suelo de los años 1998, 2005 y 2015 y la distribución de la precipitación en el área de estudio. Se delimitaron por la similitud de las propiedades analizadas. Las características de cada polígono que formaría cada zona de balance, se describen en la tabla 3.

Figura 2. Zonas de recarga hídrica en la microcuenca del río Porrosatí.

Fuente: Elaboración propia.

Muestreo de suelo

Se realizó a lo largo de toda la microcuenca, para cada uno de los diferentes usos encontrados, utilizando cilindros para obtener las respectivas muestras. En total se seleccionaron 6 sitios de muestreo debidamente georreferenciados. A las muestras colectadas se les ejecutaron los siguientes análisis: conductividad hidráulica, capacidad de campo, punto de marchitez y granulometría, siguiendo los procedimientos descritos por Henríquez & Cabalceta (1999). Las pruebas se realizaron en el Laboratorio de Hidrogeología y Recursos Hídricos de la Escuela de Ciencias Ambientales en la Universidad Nacional de Costa Rica (Tabla 4).

Tabla 4. Parámetros y métodos empleados para la determinación de condiciones hidráulicas del suelo.

Parámetro	Método	Ecuación
Conductividad hidráulica	Determinación en laboratorio mediante la construcción de permeámetro de flujo constante con cilindro de muestra.	(2)
Capacidad de campo	Determinación en laboratorio. Saturación de cilindros para su posterior extracción de agua mediante aplicación de presión a 0.33 atm durante 72 horas.	(3)
Punto de marchitez	Determinación en laboratorio. Saturación de cilindros para su posterior extracción de agua mediante aplicación de presión a 15 atm durante 72 horas.	(4)

Fuente: Elaboración propia.

Escenarios de recarga y demanda hídrica en los tractos temporales 2025-2030 y 2050-2055

Para determinar los escenarios de recarga hídrica, se utilizó el modelo propuesto por Schosinsky (2006), variando los parámetros de uso del suelo y empleando los escenarios de temperatura y precipitación elaborados por Alvarado et al. (2012). Las principales propiedades del suelo como la textura, capacidad de campo, punto de marchitez y conductividad hidráulica son constituidas a través de complejos procesos pedogénicos durante prolongados periodos de tiempo, por lo que se tomaron como constantes para su utilización dentro de las modelaciones (Pritchett 1986; Núñez 1981).

Con los datos de climatología se incorporó la incidencia del cambio climático, tomando como base las proyecciones de temperatura y precipitación modeladas por el IMN en el escenario de emisiones A2 del IPCC (Alvarado et al. 2012), en el cual se proyecta una población mundial en continuo crecimiento con un panorama de desarrollo económico y tecnológico fragmentado sin mayores cambios respecto a acciones de mitigación. La climatología base 1950 al 2000 utilizada por Alvarado et al (2012) se obtuvo en cuadrículas de 1 kilómetro por 1 kilómetro. Una vez seleccionado el valor correspondiente a la ubicación geográfica de cada una de las tres estaciones, se utilizó interpolación por polígonos de Thiessen. A los valores mensuales de la climatología base se les sumó el valor de anomalía generado por Alvarado et al (2012) con el modelo PRECIS, los cuales se presentaron en escala de 50 x 50 km. Una vez obtenidos los escenarios de temperatura, se calcularon los escenarios de evapotranspiración por el método Thornthwaite(1948).

Uso de la tierra

Para obtener los porcentajes de cobertura de uso de la tierra, se realizó una clasificación del uso de la tierra mediante la fotointerpretación en sistemas de información geográfica. Con el fin de obtener datos de alta precisión y poder evidenciar en los resultados el efecto de los cambios de uso, principalmente la impermeabilización por efecto de la urbanización, se utilizaron imágenes áreas de los años Terra 1998, Carta 2005 e imágenes georreferenciadas de alta resolución de 2015 obtenidas con de Bing Maps, con el software SAS Planet. La metodología empleada se basó en la realización de mapas de polígonos vectoriales, en donde se utilizó una generalización entre los usos Bosque, Plantación de ciprés, Pastos, Cultivos y Urbano.

Para la determinación del avance en el cambio del uso del suelo, tanto para el periodo comprendido entre 2000 a 2014 como los escenarios en el mediano y largo plazo, se calculó un coeficiente de cambio de área de uso. Dicho coeficiente se determinó a partir del área de uso urbano, en donde se establecieron tres años base para el cálculo, los cuales fueron 1998, 2005 y 2015. La escogencia de los años base estuvo limitada a la disponibilidad de imágenes áreas. El valor es multiplicado por la cantidad de años que se desea conocer y es expresado matemáticamente de la siguiente manera:

Cimp: coeficiente de impermeabilización (m2)

Ca: área de recarga en el año base inicial (m2)

Cb: área de recarga en el año base final (m2)

Cc: número de años entre Ca y Cb (años)

|En la tabla 1 se sintetizan los resultados de la clasificación del uso del suelo en la microcuenca del río Porrosatí para el año 2015. La mayor área de uso está dedicado a cultivos como café, tomate, cebolla, ornamentales y algunas otras hortalizas. La clasificación de bosques es la segunda en importancia. Los pastos representan el 12,7 % del área de la microcuenca en donde aún persisten terrenos pequeños para la producción de leche y queso. El área urbana cubre un 10,4 %, con focos dispersos en las zonas media y baja. La cuenca está situada en una dinámica rural con tendencias a la urbanización, similar a lo encontrado por Ureña (2005) para la microcuenca del río Ciruelas.

Tabla 1. Proporción en usos del suelo en la microcuenca del río Porrosatí, 2015

Uso de la tierra	Área (km2)	Porcentaje (%)
Bosque	6,7	24,0
Plantación	3,8	13,7
Pasto	3,6	12,7
Cultivo	10,1	36,1
Urbano	3,8	13,5
Total	28	100

Fuente: Elaboración propia con base en Imagen satelital Bing Maps (2015).

Resultados

Balance y comportamiento hidrológico en el periodo 2000-2014 en la microcuenca del río Porrosatí

En la figura 3 se muestra el comportamiento mensual de la recarga potencial, obtenido a partir de los registros climáticos de precipitación y temperatura para el período 2000-2014. Entre los meses de diciembre y abril se muestran valores de recarga potencial que no exceden 1 m3 en ninguno de los años del periodo 2000–2014, mientras que de mayo a noviembre aumentan considerablemente conforme avanza la época lluviosa, siendo los meses de setiembre a octubre los que presentan los valores más altos.

Figura 3. Valores de recarga potencial mensual en el periodo 2000–2014 en la microcuenca del río Porrosatí.

Fuente: Elaboración propia.

Por su parte, el análisis de la dinámica de cambio de uso de la tierra en el periodo en estudio muestra la transformación progresiva de terrenos dedicados a cultivos hacia áreas urbanas. El área urbanizada tiende a aumentar de manera más acelerada en las zonas de la 1 a la 5 en el periodo 2000–2005, mientras que en el polígono 6, el cambio en el intervalo de 2005 a 2014 fue bastante más elevado en relación con los demás. El incremento del área urbana repercute en la impermeabilización inmediata de la tierra, viéndose reducida el área efectiva de recarga como se muestra en la tabla 5.

Tabla 5. Área efectiva de recarga para cada zona de balance hídrico en la microcuenca del río Porrosatí (km2)

	Zonas de recarga
Año	1	2	3	4	5	6
2000	3,63	6,60	3,04	2,58	6,69	2,87
2001	3,63	6,59	3,03	2,56	6,68	2,87
2002	3,63	6,58	3,02	2,55	6,67	2,86
2003	3,63	6,57	3,02	2,53	6,66	2,86
2004	3,63	6,56	3,01	2,51	6,65	2,85
2005	3,63	6,54	3,00	2,50	6,64	2,85
2006	3,63	6,54	3,00	2,49	6,60	2,83
2007	3,63	6,54	3,00	2,48	6,57	2,82
2008	3,63	6,54	3,00	2,46	6,53	2,80
2009	3,63	6,54	3,00	2,45	6,50	2,79
2010	3,63	6,54	2,99	2,44	6,46	2,77
2011	3,63	6,54	2,99	2,43	6,43	2,76
2012	3,63	6,54	2,99	2,42	6,39	2,74
2013	3,63	6,54	2,99	2,41	6,35	2,73
2014	3,63	6,54	2,99	2,40	6,32	2,71

Fuente: Elaboración propia

El decrecimiento de las áreas de recarga se dio principalmente en las zonas medias y bajas de la microcuenca (zonas 2 a 6). El volumen de recarga anual promedio del agua subterránea fue de 3,04 km3 con una desviación estándar de 0,58 km3. La localización geográfica y las condiciones hidrogeológicas de la microcuenca, la ubican en un área de alta recarga, subyaciendo en uno de los reservorios de agua subterráneos de más relevancia en Centroamérica como lo es el sistema acuífero Barva-Colima (Reynolds 2002; Ramírez 2007). Sobresalen los años 2007, 2008 y 2010 como los de mayor recarga anual, mientras que los años 2000, 2001 y 2009 fueron los de menor recarga. El año de máxima recarga se dio en el año 2008 sobrepasando notablemente los valores de los demás años, mientras que el año 2009 representó el año de menor recarga. No se nota un efecto en los valores de recarga con respecto al aumento al área impermeabilizada hacia los últimos años en análisis como podría ser esperable considerando las áreas. Por tanto, en la sensibilidad de la determinación de los volúmenes de recarga resultan de mayor peso las condiciones climáticas que las áreas efectivas de recarga.

El comportamiento del volumen de recarga en contraposición con el Índice Oceánico del Niño (ONI por sus siglas en inglés) en el periodo 2000-2015 (figura 4) muestra la sensibilidad del proceso de recarga ante las variaciones climáticas inducidas por los fenómenos ENOS. El fenómeno del Niño se presentó en los años 2004-2005 y 2006-2007 en forma débil y en 2002-2003 y 2009-2010 en forma moderada. A su vez, los periodos comprendidos entre 2000-2001 y 2011-2012 fueron clasificados como eventos de la Niña, débiles, mientras que los años 2007-2008 y 2010-2011 fueron moderados. En el Valle Central de Costa Rica, el fenómeno del Niño está relacionado con períodos secos más prolongados de lo habitual, mientras que los años La Niña presentan estaciones lluviosas de mayor intensidad. Esta relación muestra que los picos y depresiones en los valores de recarga estuvieron relacionados con años de ocurrencia de la Niña y el Niño respectivamente.

Figura 4. Relación de volumen de recarga en la microcuenca del río Porrosatí y el Índice Oceánico del Niño en el período 2000-2015.

Fuente: Elaboración propia y datos NOAA (2016).

Escenarios de recarga hídrica en los periodos 2026-2030 y 2051-2055

Los escenarios de recarga hídrica se elaboraron en periodos de 5 años, para el mediano (2025-2030) y largo plazo (2050-2055). Se obtuvieron perspectivas del comportamiento de la recarga ante el cambio de uso del suelo y las afectaciones derivadas de las variaciones en los patrones climáticos. La tendencia en cuanto al uso del suelo indica aumentos significativos en el área urbanizada. El mayor crecimiento del área impermeabilizada se da en las zonas 5 y 6, ubicadas en la parte baja de la microcuenca. El principal cambio de uso ocurre en las áreas de uso de pasto y agrícolas hacia uso urbano, mientras que en las zonas 1 y 2 es bajo (tabla 6).

Tabla 6. Escenarios de área de recarga por zonas para el mediano plazo, periodo 2026-2030 (km2).

	2026	2027	2028	2029	2030
Zona 1	3,62	3,62	3,62	3,62	3,62
Zona 2	6,74	6,73	6,73	6,72	6,72
Zona 3	3,22	3,21	3,21	3,21	3,20
Zona 4	2,54	2,52	2,51	2,50	2,49
Zona 5	6,73	6,68	6,63	6,59	6,54
Zona 6	2,91	2,90	2,89	2,88	2,87

Fuente: Elaboración propia

El volumen de recarga en el mediano plazo (2026-2030) presenta una tendencia a la homogenización de los valores de recarga anuales (figura 5). La poca fluctuación responde a la inexistencia de años sobresalientemente lluviosos o secos. El efecto del aumento del área impermeabilizada provocó grandes pérdidas de recarga, siendo una limitante para la disponibilidad de recurso subterráneo en la microcuenca.

Figura 5. Escenario de volumen de recarga (hm3) para los periodos 2025-2030 y 2051-2055 en la microcuenca del río Porrosatí. Año 1; 2026, 2051: año 2; 2027, 2052; año 3: 2028, 2053; año 4: 2029, 2054; año 5: 2030, 2055

Fuente: Elaboración propia.

En el periodo 2051-2055 el área de recarga disminuye drásticamente en todas las zonas a excepción de la zona 1 (Tabla 7). El supuesto utilizado, caracterizado por la ausencia de regulaciones para las áreas medias y bajas de la microcuenca, también es un factor que incide en este resultado.

Tabla 7. Escenarios de área de recarga por zonas para el largo plazo, periodo 2051-2055 (km2)

	2051	2052	2053	2054	2055
Zona 1	3,62	3,62	3,62	3,62	3,62
Zona 2	6,62	6,62	6,61	6,61	6,60
Zona 3	3,11	3,10	3,10	3,09	3,09
Zona 4	2,23	2,21	2,20	2,19	2,18
Zona 5	5,55	5,50	5,45	5,40	5,36
Zona 6	2,66	2,65	2,64	2,63	2,62

Fuente: Elaboración propia

El comportamiento del volumen de recarga en el largo plazo (2051-2055), es influenciado principalmente por el decrecimiento de las áreas de recarga y las variaciones climáticas proyectadas (figura 5). La comparación permite visualizar una variabilidad más drástica en este periodo con respecto al mediano plazo (2026-2030). Así, se observan mínimos de hasta 42 hm3.

Los volúmenes de recarga sufren reducciones importantes por el notorio efecto del incremento del área impermeabilizada y la variabilidad climática.

Al comparar el promedio de los 15 años de registros con los escenarios de recarga propuestos, se muestra un comportamiento similar en los tres periodos con respecto a la estacionalidad de la recarga (figura 6). Sin embargo, es notorio que el pico iniciado en mayo con la entrada de la época lluviosa aumenta su intensidad en el mediano plazo y de forma más notoria en el escenario a largo plazo, llegando a ser similar al pico de setiembre octubre. La existencia de dos picos de lluvia de similar intensidad durante el año interrumpidos por la canícula de julio, indica la existencia de dos épocas de mayor recarga en el año. Esto sin duda modificará la respuesta hidrológica de los manantiales dentro de la microcuenca y el nivel de la lámina de agua con ello, la alteración de la disponibilidad y suministro de agua en la microcuenca.

La tendencia de la homogeneidad de los datos en el mediano y largo plazo puede deberse a distintas razonas. Por un lado, las posibilidades de los modelos climáticos para predecir la ocurrencia y magnitud de eventos atmosféricos como los de El Niño o La Niña se ven comprometidas al tratarse de anomalías oceánicas y atmosféricas de difícil predicción (Bindoff et al. 2013).

Figura 6. Comportamiento mensual del volumen de recarga en periodos de 5 años.

Fuente: Elaboración propia.

Los volúmenes de recarga proyectados a futuro dan cuenta de importantes decrecimientos, que en algunos casos superaron el 15 % como media y llegando hasta diferencias del 50 % con respecto al comportamiento promedio en los años 2000 a 2014. En general, las canículas en el mediano y largo plazo no parecen aumentar su intensidad en cuanto a lo observado en 2000-2014.

En primera instancia, se podría prever un comportamiento con anomalías de poca a mediana intensidad durante el año. La relativa homogeneidad de los escenarios dificulta el pronóstico de la respuesta de la producción de agua de los manantiales y la tabla de agua subterránea. El promedio baja de 25,4 hm3 en el promedio 2000-2014 a 23,6 y 21,9 hm3 en el mediano y largo plazo respectivamente. Así también la desviación estándar es menor (tabla 18).

Tabla 1. Promedio del volumen de recarga entre los períodos 2000-2014 y los escenarios de mediano y largo plazo

	Periodos de recarga (hm3)
	2000-2004	2005-2009	2010-2014	2000-2014	2025-2029	2050-2054
Promedio	23,3	26,6	26,2	25,4	23,6	21,9
Desv. Std.	21,4	23,7	25,5	21,3	18,3	17,7

Fuente: Elaboración propia

Conclusiones

Los sistemas de agua subterránea del Valle Central sobresalen por ser uno de los sistemas acuíferos más relevantes en Centroamérica. En el Gran Área Metropolitana (GAM), más del 65 % de la población es abastecida con agua proveniente de fuentes subterráneas (Reynolds 2002; Reynolds-Vargas & Fraile 2009). En este sentido, lo encontrado en los resultados coincide con Ramírez (2007) y Castro (2011), en resaltar las partes altas del macizo del Barva como las zonas en donde ocurren los mayores volúmenes de recarga al agua subterránea.

La respuesta obtenida en los periodos extremos secos o lluviosos de los años 2000 a 2014, muestra la sensibilidad del modelo y con ello del comportamiento del sistema hidrológico estudiado. La incidencia de los fenómenos ENOS sobre la variabilidad de la recarga quedó demostrada pese a que en el periodo 2000-2014 solo se dieron eventos de intensidad moderada y débil.

En cuanto a las proyecciones de la recarga en el mediano y largo plazo, la existencia de dos picos de lluvia de similar intensidad durante el año interrumpidos por la canícula de julio, indica la potencial alteración de la estacionalidad del ciclo hidrológico y con ello de la recarga de agua subterránea. Esto sin duda modificará la respuesta hidrológica y el nivel de la lámina de agua. A la vez, ésta conllevara potencialmente a la alteración de la dinámica del abastecimiento y uso del agua en la microcuenca.

Los volúmenes de recarga proyectados a futuro dan cuenta de importantes reducciones que superan el 15 % como media y alcanzan diferencias de hasta 50 % con respecto al comportamiento promedio en los años 2000 a 2014. En general, las canículas en el mediano y largo plazo no parecen aumentar su intensidad en cuanto a lo observado en 2000-2014.

En los escenarios de los periodos 2025-2030 y 2050-2055, la incidencia de condiciones extremas no se visualiza con claridad, posiblemente por la dificultad que representa el modelado de fenómenos climáticos complejos como ENOS. Sin embargo, la comunidad científica internacional ha estimado como probable el aumento en intensidad y magnitud de eventos climáticos extremos (Bindoff et al. 2013).

Como se mencionó en las secciones respectivas, los datos regionalizados están proyectados bajo un escenario de emisiones A2, en el cual, continuaran la emisión de gases con efecto invernadero a un ritmo similar al actual, siendo el escenario menos optimista. La evaluación del escenario A2 permite prever las condiciones más drásticas y realizar ejercicios de planificación bajo este umbral.

Las principales limitantes metodológicas están relacionadas con la incertidumbre asociada a los escenarios climáticos, el uso del suelo y variables hidrogeológicas, para que permitan alimentar el modelo hidrológico de la manera más precisa posible.

Agradecimientos

A las instituciones Instituto Meteorológico Nacional (IMN), Servicio Nacional de Aguas Subterráneas Riego y Avenamiento (SENARA) y Dirección de Aguas de Ministerio de Ambiente y Energía (MINAE), por la atenta colaboración y disponibilidad en el uso de datos.

Referencias

Alvarado, L., Contreras, C., Alfaro, J. y Jiménez, S. (2012). Escenarios de cambio climático regionalizados para Costa Rica. San José, Costa Rica, IMN. 1060 p.

Bates, B.C., Kundzewicz, Z.W., Wu, S. y Palutikof, J.P. (2008). Climate Change and Water. Technical Paper of the Intergovernmental Panel on Climate Change, IPCC Secretariat, Geneva (2008) 210 pp.

Bergoeing, W. (2007). Geomorfología de Costa Rica. 2da edición. San José, Costa Rica, Librería Francesa. 328 p.

Christensen, J.H., Hewitson, B., Busuioc, A., Chen, A., Gao, X., Held, I., Jones, R., Kolli, R. K., Kwon, W.-T., Laprise, R., Magaña Rueda, V., Mearns, L., Menéndez, C.G., Räisänen, J., Rinke, A., Sarr, A. and Whetton, P. (2007). Regional Climate Projections. In: Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change [Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K.B., Tignor, M. and H.L. Miller (eds.)]. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA.

Denyer, P. y Kussmaul, D. (2000). Geología de Costa Rica. Cartago, Costa Rica, Editorial Tecnológica. 515 p.

Díaz, J. M. G. (2003).Métodos de estimación y ajuste de datos climáticos. (80). CDCH UCV.

Henríquez, J. y Cabalceta, E. (1999). Guía práctica para el estudio introductorio de los suelos con un enfoque agrícola. Asociación Costarricense de la Ciencia del Suelo. 111 p.

Hijmans, R., Cameron, S., Parra, J., Jones, P. & Jarvis, A. (2005). Very High Resolution Interpolated Climate Surfaces for Global Land Areas. International Journal of Climatology 25, 1965-1978.

Kusangaya, S., Warburton, M.L., van Garderen, E.A. & Jewitt, G.P. (2014). Impacts of climate change on water resources in southern Africa: a review. Phys. Chem. Earth 67, 47-54.

Mata, R. y Ramírez, J. (1999). Estudio de caracterización de suelos y su relación con el manejo del cultivo de Café en la provincia de Heredia. San José, Costa Rica, Instituto del Café de Costa Rica. 92 p.

Núñez, J. (1981). Fundamentos de edafología. San José, Costa Rica, EUNED. 220 p.

Pritchett, W. (1986). Suelos forestales: Propiedades, conservación y mejoramiento. México, Editorial Limusa.

Ramírez, C. (2007). Recarga potencial del acuífero Colima y Barva, Valle Central. Costa Rica, SENARA. 41 p.

Reynolds, J. (2002). Manejo integrado de aguas subterráneas: un reto para el futuro. San José, Costa Rica, EUNED. 325 p.

Schosinsky, G. y Losilla, M. (2000). Modelo analítico para determinar la infiltración con base en la lluvia mensual. Revista Geológica de América Central (23): 43-55.

Schosinsky, G. (2006). Cálculo de la recarga potencial de acuíferos mediante un balance hídrico de suelos. Revista Geológica de América Central, (34-35): 13-30.

SENARA, B. (1985). Mapa hidrogeológico del Valle Central de Costa Rica. Escala, 1(50), 000.

Sibaja, P. (2013). El crecimiento urbano y la producción hídrica en la cuenca del río Porrosatí, provincia de Heredia, Costa Rica. (Tesis de Licenciatura de Ciencias Geográficas). San José, Costa Rica, Universidad Nacional de Costa Rica. 150 p.

Smith, A.B. y Katz, R.W. (2013). US billion-dollar weather and climate disasters: data sources, trends, accuracy and biases. Nat. Hazards, 67, 387–410.

Thornthwaite, C. W. (1948). An approach toward a rational classification of climate. Geographical review, 38(1), 55-94.

U.S. Geological Survey (USGS). (1997). Geology and resource assessment of Costa Rica at 1:500,000 scale; a digital representation of maps of the U.S. Geological Survey's 1987 folio I-1865. Dirección General de Geología y Minas e Hidrocarburos; Universidad de Costa Rica

Ureña, N. (2005). Efectos del aumento poblacional y del cambio del uso del suelo sobre los recursos hídricos en la microcuenca del río Ciruelas, Costa Rica. (Tesis de Maestría). Turrialba, Cartago, Costa Rica, CATIE. 170 p.

Van der Velde, M., Tubiello, F.N., Vrieling, A. y Bouraoui, F. (2012). Impacts of extreme weather on wheat and maize in France: evaluating regional crop simulations against observed data. Clim. Change 113, 751–765.

Zhou, Y., Zwahlen, F., Wang, Y., Li, Y. (2010). Impact of climate change on irrigation requirements in terms of groundwater resources. Hydrogeology Journal 18: 1571-1582.

---