Efectos del cambio climático sobre la disponibilidad de recursos hídricos en la región litoral-mesopotamia (Argentina)

Ingenier´ıa del Agua, Vol. 17, No 2, Junio 2010 EFECTOS DEL CAMBIO CLIM´ ATICO SOBRE LA DISPONIBILIDAD DE RECURSOS H´ IDRICOS EN LA REGI´ ON LITORALMESOPOTAMIA (ARGENTINA) Ra´ul Pedraza, Ma. del Valle Morresi, Graciela Pusineri y M´onica D’El´ıa Facultad de Ingenier´ıa y Ciencias H´ıdricas, Universidad Nacional del Litoral Ciudad Universitaria, 3000 Santa Fe, Argentina rp[email protected]r Resumen: La regi´on Litoral-Mesopotamia argentina (LM) es la zona m´as rica del pa´ıs en cuanto a sus recursos h´ıdricos. El objetivo de este trabajo fue evaluar los efectos del Cambio Clim´atico sobre la disponibilidad de los recursos h´ıdricos de la regi´on, con un horizonte temporal hasta la d´ecada 2081-2090. Para evaluar tales efectos, se aplic´o un modelo de balance h´ıdrico mensual en 22 estaciones representativas de la regi´on, para el a˜no medio del per´ıodo de referencia 1981/1990, representativo de la situaci´on actual, y para el a˜no medio del per´ıodo 2081/2090, considerando para este ´ultimo dos escenarios de emisi´on de gases de efecto invernadero, A2 y B2, definidos por el Panel Intergubernamental de Expertos sobre el Cambio Clim´atico (IPCC). A tal efecto, se utilizaron datos observados de precipitaci´on y temperatura mensuales en dichas estaciones, as´ı como datos sint´eticos de esas variables calculados con el Modelo Clim´atico Regional del Centro de Investigaci´on del Mar y la Atm´osfera (CIMA). Los resultados obtenidos indican una dr´astica disminuci´on de los excesos y un aumento generalizado del d´eficit h´ıdrico en la regi´on, siendo m´as severos los efectos para el Escenario A2. Estos cambios se producir´ıan, fundamentalmente, por aumento de la evapotranspiraci´on, debido al aumento de la temperatura del aire. Los cambios implicar´ıan una reducci´on tanto de la escorrent´ıa directa como de la recarga subterr´anea. Los cursos aut´octonos, cuya escorrent´ıa se genera por lluvias ca´ıdas dentro de la regi´on en estudio, ver´ıan reducidos sus caudales medios, lo que implicar´ıa un cambio en la tendencia de aumento registrada en las ´ultimas 3 d´ecadas. Se identifican sectores de la regi´on donde el severo aumento del d´eficit podr´ıa generar dificultades para el abastecimiento de agua para consumo humano y la posible extinci´on o modificaci´on de ecosistemas. Se concluye en la necesidad de planificar medidas de adaptaci´on para mitigar los efectos del cambio clim´atico en la regi´on. INTRODUCCI´ ON Este estudio se desarroll´o en el marco del Proyecto “Vulnerabilidad de los Recursos H´ıdricos en el Litoral – Mesopotamia (LM)” (FICH, 2006), desarrollado como parte de las “Actividades Habilitantes para la 2a. Comunicaci´on Nacional del Gobierno de la Rep´ublica Argentina a las Partes de la Convenci´on Marco de las Naciones Unidas sobre Cambio Clim´atico - TF 51287/AR”. El objetivo del proyecto fue evaluar el impacto sobre el clima y los recursos h´ıdricos de la regi´on LM que producir´a el Cambio Clim´atico, con un horizonte temporal hasta la d´ecada 2081-2090. Para ello se consideraron los escenarios de emisi´on de gases de efecto invernadero, A2 y B2, definidos por el Panel Intergubernamental de Expertos sobre el Cambio Clim´atico (IPCC, 2000). Estos escenarios no representan condiciones extremas, sino condiciones medias altas y medias bajas de emisi´on de di´oxido de carbono, respectivamente, y se muestran en la Figura 1. La familia de escenarios A2 supone un mundo heterog´eneo cuyas caracter´ısticas m´as distintivas son: a) la preservaci´on de las identidades locales, b) una alta tasa de crecimiento poblacional, c) un desarrollo econ´omico regional y d) un desarrollo econ´omico menor que en otras l´ıneas narrativas. © Fundaci´ on para el Fomento de la Ingenier´ ıa del Agua ISSN: 1134–2196 Recibido: Diciembre 2009 Aceptado: June 2010 104 R. Pedraza, M. del Valle, G. Pusineri y M. D’El´ıa Figura 1. Escenarios de Emisiones de CO2(fuente: IPCC, 2000) P: precipitación, Etr: evapotranspiración real, Exc: exceso de agua (Exc = Esc + Per), R: reserva de agua utilizable por las plantas, Esc: escorrentía directa, Per: percolación al acuífero P Etr R Esc Per Exc Almacenamiento Subterráneo Nivel Terreno Volumen de control Figura 2. Esquema del volumen de control y variables consideradas en el balance h´ıdrico La familia de escenarios B2 supone un mundo con ´enfasis en las soluciones locales cuyas caracter´ısticas m´as importantes son: a) un aumento menor de la poblaci´on que en A2 y b) niveles intermedios de desarrollo econ´omico. OBJETIVOS El objetivo general de este estudio fue evaluar los efectos del cambio clim´atico sobre la disponibilidad de los recursos h´ıdricos en la regi´on LM, con un horizonte temporal hasta la d´ecada 20812090. Los objetivos particulares fueron evaluar los efectos sobre la distribuci´on areal de excesos y d´eficit h´ıdricos en la regi´on. MODELO DE BALANCE H´ IDRICO APLICADO Para evaluar los efectos mencionados, se aplic´o un modelo de balance h´ıdrico mensual en estaciones representativas de la regi´on, para el a˜no medio del per´ıodo de referencia 1981/1990, representativo de la situaci´on actual y para el a˜no medio del per´ıodo 2081/2090, considerando para este ´ultimo los escenarios A2 y B2. El modelo aplica la ecuaci´on de continuidad a un volumen de control ubicado en la subzona de evapotranspiraci´on de la zona no saturada del suelo, comprendida entre la superficie del terreno y la profundidad radicular media. Figura 2 muestra un esquema del volumen de control y las variables consideradas en el balance h´ıdrico. Efectos del cambio clim´atico sobre la disponibilidad de recursos h´ıdricos en la regi´on litoralmesopotamia (Argentina) 105 La ecuaci´on de continuidad aplicada al volumen de control en un intervalo de tiempo ∆t es: P=Etr +Exc + ∆R(1) donde ∆Res la variaci´on de la reserva de agua utilizable por las plantas. Deben cumplirse las siguientes restricciones: Etr ≤Etp (2) 0≤R≤Rm´ax (3) donde Etp es la evapotranspiraci´on potencial y Rm´ax la reserva de agua m´axima utilizable por las plantas. El valor de Rm´ax se calcula: Rm´ax =Wcc −Wpmp (4) donde Wcc es la capacidad de campo del suelo yWpmp el punto de marchitez permanente. Para este estudio se adopt´o un intervalo de tiempo de c´alculo mensual. Teniendo como datos los valores mensuales de Etp yP, para cada mes se calcularon los valores de Etr,R,Exc y D´ef. De acuerdo a la relaci´on entre PyEtp en cada mes, se consideraron 2 situaciones: a) P < Etp (P−Etp < 0): no se satisface Etp con P, el almacenamiento suministra agua para satisfacerla y disminuye la reserva (∆Rnegativo). Rtiene un l´ımite m´ınimo igual a cero. Cuando la reserva se agota y a´un no se ha satisfecho Etp, se produce un d´eficit (D´ef). b) P≥Etp (P−Etp ≥0): se satisface Etp (Etr =Etp) y el sobrante de agua aumenta la reserva de agua ´util (∆Rpositivo). Si se completa Rm´ax y sigue existiendo sobrante, ´este pasa a ser exceso (escurrimiento directo m´as percolaci´on). Por tratarse el a˜no medio de un ciclo cerrado, se aplic´o un procedimiento de c´alculo iterativo. Para un mes gen´erico i, se calcul´o ∆Rcomo se describe a continuaci´on: a) Si (P−Etp)i<0: Si |P−Etp|i=Ri−1⇒∆Ri= (P−Etp)i(5) Si |P−Etp|i> Ri−1⇒(6) b) Si (P−Etp)i= 0: Si (P−Etp)i= (Rm´ax −Ri−1)⇒∆Ri = (P−Etp)i(7) Si (P−Etp)i>(Rm´ax −Ri−1)⇒∆Ri =Rm´ax −Ri−1(8) En cualquiera de estos casos, se actualiz´o el valor de Ral final del mes icon: Ri=Ri−1+ ∆Ri(9) Finalmente, se calcularon los valores de Exc, Etr yD´ef del mes i: Si (P−Etp)i= 0 ⇒Exci= 0 (10) Si Ri=Rm´axy (P−Etp)i>0⇒Exci = (P−Etp)i−∆Ri(11) Si (P−Etp)i<0⇒Etri=Pi+|∆Ri|(12) Si (P−Etp)i= 0 ⇒Etri=Etpi(13) Defi=Etpi−Etri(14) Los valores de Etp se calcularon en base a los datos mensuales de temperatura del aire, aplicando el m´etodo de Thornthwaite (Thornthwaite y Holzman, 1942). Estos autores correlacionaron valores de temperatura del aire medias mensuales con valores de Etp determinados por balance h´ıdrico en cuencas de la costa Este central de EUA (zona h´umeda), con pasturas bajas. Sus resultados fueron obtenidos bajo condiciones en las que la tasa de evaporaci´on no estaba influenciada por advecci´on de aire h´umedo o seco (Thornthwaite y Mather, 1955). Debido a que estas condiciones se cumplen en la mayor parte de la regi´on LM y dado que no se dispuso de datos simulados a futuro de otras variables meteorol´ogicas, como heliofan´ıa, velocidad del viento y humedad relativa, se adopt´o el m´etodo mencionado. APLICACI´ ON DEL MODELO AL PER´ IODO DE REFERENCIA 1981/1990 A efectos de determinar una distribuci´on areal de excesos y d´eficit h´ıdricos en la regi´on LM, representativa de la situaci´on actual, se aplic´o el modelo de balance h´ıdrico en 22 estaciones de la regi´on (Figura 3), en base a datos de precipitaci´on y temperatura observados para el a˜no medio del per´ıodo de referencia 1981/1990. 106 R. Pedraza, M. del Valle, G. Pusineri y M. D’El´ıa Figura 3. Ubicaci´on de estaciones utilizadas para el balance h´ıdrico Para cada estaci´on, se estim´o el valor del par´ametro Rm´ax. Para ello, se determin´o la textura dominante (superficial y subsuperficial) en cada estaci´on a partir del Atlas de Suelos de la Rep´ublica Argentina (INTA, 1995), a escala 1:500.000. En base a la textura del suelo en cada estaci´on, se obtuvieron de tablas los l´ımites inferior y superior del rango posible de Rm´ax, para una profundidad est´andar de 1 m (Grassi, 1998). El valor de Rm´ax en cada estaci´on, para la situaci´on actual de los suelos, se adopt´o igual al promedio entre los l´ımites inferior y superior sugeridos. Tabla 1 muestra un listado de las estaciones consideradas en el estudio, los valores de Rm´ax adoptados y las fuentes de datos de precipitaci´on y temperatura medias mensuales. Los valores medios anuales de precipitaci´on y temperatura y los valores de exceso y d´eficit anuales calculados con el modelo se muestran en la Tabla 2. A modo comparativo, tambi´en se aplic´o el modelo con los datos sint´eticos de precipitaci´on y temperatura medias mensuales generados con el modelo CIMA (CIMA, 2005) para el mismo per´ıodo. Estaci´on Provincia Rm´ax [mm] Fuente de datos de P y T Bella Vista Corrientes 190 Bella Vista INTA (*1). Ceres Santa Fe 190 Ceres Aero SMN (*1). Concordia Entre R´ıos 80 Concordia Aero SMN (*1). Corrientes Corrientes 190 Corrientes Aero SMN (*1). El Colorado Formosa 190 P: DPAyS, T: interpoladas entre Formosa y Pres. R. S. Pe˜na. Formosa Formosa 190 Formosa Aero SMN (*1). Gob. Virasoro Corrientes 190 P: Est. Las Mar´ıas, T: interpoladas entre Posadas y Bella Vista. Gualeguaych´u Entre R´ıos 230 Gualeguaych´u Aero SMN (*1). Ituzaing´o Corrientes 190 P: EVARSA, T: interpoladas entre Posadas y Bella Vista. J. J. Castelli Chaco 190 P: APA, T: interpoladas entre Las Lomitas y Pres. R. S. Pe˜na. Las Bre˜nas Chaco 190 Las Bre˜nas INTA (*1). Las Lomitas Formosa 80 Las Lomitas SMN (*1). Paran´a Entre R´ıos 230 Paran´a INTA (*1). Paso de los Libres Corrientes 230 Paso de los Libres Aero SMN (*1). Posadas Misiones 190 Posadas Aero SMN (*1). Pres. R. S. Pe˜na Chaco 190 Pres. R.S. Pe˜na INTA (*1). Rafaela Santa Fe 230 Rafaela INTA (*1). Reconquista Santa Fe 186 Reconquista Aero SMN (*1). Resistencia Chaco 230 Resistencia Aero SMN (*1). Rosario Santa Fe 190 Rosario Aero SMN (*1). Santa Fe Santa Fe 230 P: FICH, T: Sauce Viejo Aero SMN (*1). Villa ´ Angela Chaco 190 Villa Angela SMN (*1). P: precipitaciones medias mensuales, T: temperaturas medias mensuales, INTA: Instituto Nacional de Tecnolog´ıa Agropecuaria, SMN: Servicio Meteorol´ogico Nacional, DPAyS: Direcci´on Provincial de Agua y Saneamiento de Formosa, EVARSA: Evaluaci´on de Recursos S.A., APA: Administraci´on Provincial del Agua del Chaco, FICH: Facultad de Ingenier´ıa y Ciencias H´ıdricas (*1): Estad´ısticas Climatol´ogicas 1981/1990 (SMN) Tabla 1. Listado de estaciones, valores del par´ametro Rm´ax y fuentes de datos hidrometeorol´ogicos utilizados en el balance h´ıdrico Efectos del cambio clim´atico sobre la disponibilidad de recursos h´ıdricos en la regi´on litoralmesopotamia (Argentina) 107 Estaci´on Balance con datos P y T Balance con datos P y T Errores relativos Observados Modelo CIMA P T Exc D´ef P T Exc D´ef P T Exc D´ef [mm] [ C] [mm] [mm] [mm] [ C] [mm] [mm] [ %] [ %] [ %] [ %] Bella Vista 1465.4 20.8 419.0 0.0 1418.5 20.8 362.6 0.0 -3.2 0.2 -13.5 0.0 Ceres 941.0 19.0 0.0 2.9 965.3 19.7 0.0 23.5 2.6 4.0 0.0 718.5 Concordia 1308.0 18.9 388.9 21.0 1350.8 18.8 426.3 6.6 3.3 -0.5 9.6 -68.8 Corrientes 1645.6 21.2 565.9 0.0 1500.3 21.5 396.1 0.0 -8.8 1.4 -30.0 0.0 El Colorado 1282.5 21.9 143.8 0.0 1352.3 21.9 222.3 0.0 5.4 -0.4 54.5 0.0 Formosa 1498.4 22.1 342.3 0.0 1522.2 22.1 376.0 0.0 1.6 -0.2 9.8 0.0 Gob. Virasoro 1935.1 21.1 857.9 0.0 1990.4 21.1 920.7 0.0 2.9 -0.2 7.3 0.0 Gualeguaych´u 1195.5 17.8 309.2 0.0 1171.2 17.8 288.8 0.0 -2.0 0.3 -6.6 0.0 Ituzaing´o 1854.9 21.3 766.8 0.0 1902.1 21.4 805.4 0.0 2.5 0.6 5.0 0.0 J. J. Castelli 1065.5 21.6 0.0 45.6 1128.5 22.0 0.0 12.0 5.9 1.5 0.0 -73.6 Las Bre˜nas 1079.4 20.8 31.7 0.0 1196.0 21.3 98.9 0.0 10.8 2.7 211.8 0.0 Las Lomitas 941.9 22.6 0.0 268.1 932.0 22.7 0.0 284.5 -1.1 0.6 0.0 6.1 Paran´a 1075.0 18.2 178.4 0.0 1060.6 18.7 127.8 0.0 -1.3 3.0 -28.3 0.0 P. de los Libres 1611.3 19.7 617.6 0.0 1650.7 19.8 662.2 0.0 2.4 0.3 7.2 0.0 Posadas 1947.7 21.5 847.9 0.0 1934.3 21.0 888.0 0.0 -0.7 -2.5 4.7 0.0 Pres. R. S. Pe˜na 1189.0 20.7 155.4 0.0 1270.8 21.5 167.9 0.0 6.9 3.8 8.1 0.0 Rafaela 1018.7 17.8 141.7 0.0 993.9 19.0 48.9 0.0 -2.4 6.6 -65.5 0.0 Reconquista 1380.1 19.9 379.5 0.0 1276.8 20.3 251.1 0.0 -7.5 1.8 -33.8 0.0 Resistencia 1556.7 21.0 487.4 0.0 1489.0 21.5 382.6 0.0 -4.3 2.3 -21.5 0.0 Rosario 933.2 17.3 68.8 0.0 985.9 17.8 100.4 0.0 5.6 3.1 45.9 0.0 Santa Fe 1051.8 18.9 105.3 0.0 1053.9 18.8 118.0 0.0 0.2 -0.5 12.1 0.0 Villa ´ Angela 1170.5 20.9 103.8 0.0 1255.4 21.0 185.4 0.0 7.2 0.6 78.6 0.0 Promedio 1.2 1.3 11.6 26.5 P: precipitaci´on anual, T: temperatura media anual, Exc: exceso anual, D´ef: d´eficit anual Tabla 2. Valores de excesos y d´eficit anuales calculados con el modelo de balance para el per´ıodo 1981/1990 Los valores anuales de precipitaci´on y temperatura simulados con el modelo CIMA se aproximan a los respectivos valores observados, con errores relativos de hasta ±10 %. Sin embargo, los valores de exceso y d´eficit anuales calculados con el modelo a partir de los datos sint´eticos tienen errores relativos mayores al ±20 % en algunas estaciones, con respecto a los valores obtenidos con el mismo modelo pero en base a datos observados (por ej. en Las Bre˜nas, Villa ´ Angela, Rafaela, El Colorado, Rosario, Reconquista, Corrientes y Resistencia). Esto es debido a que el modelo de balance es no lineal; esto es, no existe proporcionalidad entre las variables de entrada y de salida. Se encontr´o que peque˜nos errores en las variables de entrada originan errores mayores en las variables de salida; es decir, se amplifican los errores. Teniendo en cuenta esto, los resultados obtenidos con el modelo en base a los datos observados del per´ıodo de referencia se adoptaron como representativos de la situaci´on actual y fueron utilizados en posteriores comparaciones con los escenarios futuros. Las Figuras 4 y 5 muestran las isol´ıneas de excesos y d´eficit medios anuales, respectivamente, calculados con el modelo para el per´ıodo 1981/1990. Los resultados obtenidos muestran la existencia de excesos h´ıdricos anuales en la mayor parte de la regi´on en estudio, con un gradiente de aumento hacia el Noreste, en concordancia con la variaci´on areal de la precipitaci´on anual media. Tambi´en se observa un d´eficit anual nulo en la mayor parte de la regi´on, con excepci´on de los sectores ubicados al Oeste de las provincias de Formosa, Chaco y Santa Fe. Los excesos medios anuales obtenidos con el modelo son compatibles con los escurrimientos medios anuales medidos en cursos de la regi´on (Subsecretar´ıa de Recursos H´ıdricos, 2004). APLICACI´ ON DEL MODELO AL PER´ IODO 2081/2090 – ESCENARIO A2 Para cada estaci´on considerada, se calcularon las diferencias entre los valores medios mensuales de precipitaci´on simulados con el modelo CIMA para los per´ıodos 2081/2090 - Escenario A2 y 1981/1990. De manera similar, se calcularon las diferencias entre los valores medios mensuales de temperatura del aire (Tablas 3 y 4). 108 R. Pedraza, M. del Valle, G. Pusineri y M. D’El´ıa Figura 4. Excesos medios anuales calculados para el per´ıodo 1981/1990 [mm] Figura 5. D´eficit medios anuales calculados para el per´ıodo 1981/1990 [mm] Estaci´on E F M A M J J A S O N D Anual Bella Vista 13.3 4.8 -4.3 12.6 25.1 -46.5 -18.8 -14.9 -7.8 -24.5 -29.1 -28.8 -118.9 Ceres 28.7 2.3 3.4 24.6 17.2 -13.7 1.5 -11.0 -2.3 -12.3 -30.4 -22.8 -14.9 Concordia 10.5 4.5 -36.7 27.9 37.7 -42.3 -10.7 2.8 -17.4 -14.4 -17.1 -18.9 -74.1 Corrientes -16.3 1.3 7.4 12.5 16.6 -57.9 -19.6 -13.1 -3.9 -20.6 -5.7 -25.1 -124.5 El Colorado -21.8 7.7 15.9 19.1 10.3 -38.4 -10.6 -8.7 -12.9 -9.9 -0.2 1.2 -48.4 Formosa -26.2 8.5 36.0 8.5 12.7 -40.5 -6.6 -12.1 -26.9 -0.4 -8.2 3.0 -52.2 Gob. Virasoro -51.6 -14.3 2.3 15.5 24.0 -61.3 -14.8 -16.9 -12.6 -16.7 -22.0 -21.8 -190.3 Gualeguaych´u -6.2 8.4 5.8 22.2 35.0 -18.2 1.5 -8.5 0.7 -20.5 -9.6 -8.4 2.4 Ituzaing´o -42.0 -20.2 8.1 15.6 28.1 -64.8 -11.3 -17.6 -12.4 -25.4 -18.6 -21.7 -182.1 J. J. Castelli -26.0 30.7 -11.3 15.5 10.7 -22.1 -10.3 -6.7 -4.1 -6.5 4.0 3.1 -23.1 Las Bre˜nas -33.7 13.8 -14.6 20.2 9.4 -22.8 -9.5 -8.1 -4.2 -25.1 -6.1 -6.3 -86.9 Las Lomitas -33.2 13.3 41.2 4.6 -10.6 -20.2 -9.5 -8.9 -13.2 8.8 11.4 29.0 12.6 Paran´a 12.2 7.3 4.9 27.5 21.9 -15.4 4.1 -4.9 -5.4 -12.4 -6.3 8.3 41.9 P.d. los Libres -5.4 20.8 -9.0 -8.3 12.6 -71.5 -21.6 -11.0 -16.5 -24.5 -16.3 -27.3 -177.9 Posadas -30.6 -14.0 13.6 11.6 32.2 -47.4 -10.5 -17.2 -18.4 -20.8 -17.9 -8.8 -128.2 P. R. S. Pe˜na -40.3 31.3 5.9 17.9 2.5 -32.4 -13.6 -8.5 -3.8 -25.1 -2.4 -18.0 -86.7 Rafaela 29.2 6.8 -6.2 18.7 20.8 -14.1 4.8 -6.5 -5.8 -9.3 -19.1 2.8 22.0 Reconquista 7.1 3.9 -15.4 26.9 24.6 -31.4 -8.5 -2.6 -7.8 -16.2 -24.5 -36.8 -80.7 Resistencia -13.3 -3.0 6.3 12.8 15.8 -57.6 -18.3 -12.5 -3.9 -20.8 -14.5 -24.0 -132.9 Rosario 9.6 15.7 -9.4 40.2 15.2 -6.2 9.5 -4.7 -8.3 -4.6 -6.0 2.1 53.0 Santa Fe 13.5 4.8 1.9 28.4 21.0 -15.1 4.0 -4.6 -4.6 -12.2 -7.0 9.0 39.0 Villa ´ Angela -20.1 -1.4 -1.4 29.1 14.0 -25.4 -8.6 -9.2 -4.6 -22.9 -10.0 -22.5 -83.2 Tabla 3. Diferencias entre precipitaciones medias mensuales simuladas con el modelo CIMA para los per´ıodos 2081/2090 - Escenario A2 y 1981/1990 [mm] La Figura 6 muestra la variaci´on entre las precipitaciones medias anuales simuladas con el modelo CIMA para los per´ıodos 2081/2090 - Escenario A2 y 1981/1990. Efectos del cambio clim´atico sobre la disponibilidad de recursos h´ıdricos en la regi´on litoralmesopotamia (Argentina) 109 Estaci´on E F M A M J J A S O N D Anual Bella Vista 13.3 4.8 -4.3 12.6 25.1 -46.5 -18.8 -14.9 -7.8 -24.5 -29.1 -28.8 -118.9 Ceres 28.7 2.3 3.4 24.6 17.2 -13.7 1.5 -11.0 -2.3 -12.3 -30.4 -22.8 -14.9 Concordia 10.5 4.5 -36.7 27.9 37.7 -42.3 -10.7 2.8 -17.4 -14.4 -17.1 -18.9 -74.1 Corrientes -16.3 1.3 7.4 12.5 16.6 -57.9 -19.6 -13.1 -3.9 -20.6 -5.7 -25.1 -124.5 El Colorado -21.8 7.7 15.9 19.1 10.3 -38.4 -10.6 -8.7 -12.9 -9.9 -0.2 1.2 -48.4 Formosa -26.2 8.5 36.0 8.5 12.7 -40.5 -6.6 -12.1 -26.9 -0.4 -8.2 3.0 -52.2 Gob. Virasoro -51.6 -14.3 2.3 15.5 24.0 -61.3 -14.8 -16.9 -12.6 -16.7 -22.0 -21.8 -190.3 Gualeguaych´u -6.2 8.4 5.8 22.2 35.0 -18.2 1.5 -8.5 0.7 -20.5 -9.6 -8.4 2.4 Ituzaing´o -42.0 -20.2 8.1 15.6 28.1 -64.8 -11.3 -17.6 -12.4 -25.4 -18.6 -21.7 -182.1 J. J. Castelli -26.0 30.7 -11.3 15.5 10.7 -22.1 -10.3 -6.7 -4.1 -6.5 4.0 3.1 -23.1 Las Bre˜nas -33.7 13.8 -14.6 20.2 9.4 -22.8 -9.5 -8.1 -4.2 -25.1 -6.1 -6.3 -86.9 Las Lomitas -33.2 13.3 41.2 4.6 -10.6 -20.2 -9.5 -8.9 -13.2 8.8 11.4 29.0 12.6 Paran´a 12.2 7.3 4.9 27.5 21.9 -15.4 4.1 -4.9 -5.4 -12.4 -6.3 8.3 41.9 P.d. los Libres -5.4 20.8 -9.0 -8.3 12.6 -71.5 -21.6 -11.0 -16.5 -24.5 -16.3 -27.3 -177.9 Posadas -30.6 -14.0 13.6 11.6 32.2 -47.4 -10.5 -17.2 -18.4 -20.8 -17.9 -8.8 -128.2 P. R. S. Pe˜na -40.3 31.3 5.9 17.9 2.5 -32.4 -13.6 -8.5 -3.8 -25.1 -2.4 -18.0 -86.7 Rafaela 29.2 6.8 -6.2 18.7 20.8 -14.1 4.8 -6.5 -5.8 -9.3 -19.1 2.8 22.0 Reconquista 7.1 3.9 -15.4 26.9 24.6 -31.4 -8.5 -2.6 -7.8 -16.2 -24.5 -36.8 -80.7 Resistencia -13.3 -3.0 6.3 12.8 15.8 -57.6 -18.3 -12.5 -3.9 -20.8 -14.5 -24.0 -132.9 Rosario 9.6 15.7 -9.4 40.2 15.2 -6.2 9.5 -4.7 -8.3 -4.6 -6.0 2.1 53.0 Santa Fe 13.5 4.8 1.9 28.4 21.0 -15.1 4.0 -4.6 -4.6 -12.2 -7.0 9.0 39.0 Villa ´ Angela -20.1 -1.4 -1.4 29.1 14.0 -25.4 -8.6 -9.2 -4.6 -22.9 -10.0 -22.5 -83.2 Tabla 4. Diferencias entre precipitaciones medias mensuales simuladas con el modelo CIMA para los per´ıodos 2081/2090 - Escenario A2 y 1981/1990 [mm] Figura 6. Variaci´on entre las precipitaciones medias anuales del per´ıodo 2081/2090 - A2 y del per´ıodo 1981/1990 [mm] Se observa que para el per´ıodo 2081/2090 - Escenario A2 se espera una leve disminuci´on de las lluvias en gran parte de la regi´on, en relaci´on al per´ıodo de referencia. Las principales reducciones se ubican al Este de Corrientes, de hasta 190 mm/a˜no. Son excepciones a esta reducci´on, los sectores Centro y Sur de la Prov. de Santa Fe y Sureste de Entre R´ıos. En estos sectores, las precipitaciones anuales se mantendr´ıan sin cambios o aumentar´ıan ligeramente respecto a las actuales, hasta en 50 mm/a˜no. La Figura 7 muestra la variaci´on entre las temperaturas medias anuales simuladas con el modelo CIMA para los per´ıodos 2081/2090 - Escenario A2 y 1981/1990. Se observa un aumento generalizado significativo de la temperatura media anual del aire, con incrementos de 3.0 ◦C al Sur y de 4.1 ◦C al Norte de la regi´on. Para cada estaci´on, se calcularon las precipitaciones medias mensuales del per´ıodo 2081/2090 - Escenario A2 como suma de las precipitaciones medias mensuales observadas en el per´ıodo de referencia m´as las diferencias mensuales de la Tabla 3. De manera similar, para cada estaci´on se calcularon las temperaturas medias mensuales del per´ıodo 2081/2090 - Escenario A2 como suma de las temperaturas medias mensuales observadas en el per´ıodo de referencia m´as las diferencias de la Tabla 4. Se aplic´o el modelo de balance a cada estaci´on, considerando como datos de entrada los valores de precipitaci´on y temperatura media mensual del per´ıodo 2081/2090 - Escenario A2 obtenidos como se ha descrito en el p´arrafo anterior. 110 R. Pedraza, M. del Valle, G. Pusineri y M. D’El´ıa Estaci´on Rm´ax Exc D´ef ∆Exc ∆Exc ∆D´ef [mm] [mm] [mm] [ %] [mm] [mm] Bella Vista 160 20.9 243.8 -95 -398.1 243.8 Ceres 160 0.0 354.2 0 0.0 351.3 Concordia 60 189.5 238.9 -51 -199.5 217.9 Corrientes 160 145.4 276.8 -74 -420.5 276.8 El Colorado 160 11.8 545.8 -92 -132.0 545.8 Formosa 160 11.3 402.0 -97 -331.0 402.0 Gob. Virasoro 160 419.7 332.2 -51 -438.2 332.2 Gualeguaych´u 200 79.3 4.3 -74 -229.9 4.3 Ituzaing´o 160 297.7 326.9 -61 -469.1 326.9 J. J. Castelli 160 0.0 654.2 0 0.0 608.6 Las Bre˜nas 160 0.0 556.8 -100 -31.7 556.8 Las Lomitas 60 0.0 961.5 0 0.0 693.4 Paran´a 200 0.0 30.2 -100 -178.4 30.2 Paso de los Libres 200 147.5 152.3 -76 -470.1 152.3 Posadas 160 351.0 254.2 -59 -496.9 254.2 Pres. R. S. Pe˜na 160 0.0 429.2 -100 -155.4 429.2 Rafaela 200 0.0 79.0 -100 -141.7 79.0 Reconquista 156 8.2 153.4 -98 -371.2 153.4 Resistencia 200 76.4 277.1 -84 -411.0 277.1 Rosario 160 0.0 79.7 -100 -68.8 79.7 Santa Fe 200 0.0 148.3 -100 -105.3 148.3 Villa ´ Angela 160 0.0 498.5 -100 -103.8 498.5 Promedio -73 -234.2 302.8 Exc: excesos medios anuales, D´ef: d´eficit medios anuales, ∆Exc y ∆D´ef: variaci´on de los excesos y d´eficit medios anuales del per´ıodo 2081/2090 - Escenario A2 respecto del per´ıodo 1981/1990 Tabla 5. Resultados obtenidos para el per´ıodo 2081/2090 - Escenario A2 Figura 7. Variaci´on entre las temperaturas medias anuales del per´ıodo 2081/2090 - A2 y del per´ıodo 1981/1990 [mm] Para cada estaci´on, se tuvo en cuenta una disminuci´on del par´ametro Rm´ax respecto del valor actual. El escenario futuro A2 est´a asociado con intensas pr´acticas agr´ıcolas, las que originar´an una degradaci´on de la estructura del suelo, por disminuci´on de materia org´anica y compactaci´on. Esta degradaci´on producir´a una menor velocidad de infiltraci´on en el horizonte superficial, una menor conductividad hidr´aulica de los horizontes subsuperficiales y una disminuci´on del agua ´util. Por esta raz´on, para cada estaci´on se adopt´o el l´ımite inferior del rango posible de Rm´ax recomendado en la literatura. Los resultados obtenidos se muestran en la Tabla 5. Las Figuras 8 y 9 muestran las isol´ıneas de excesos y d´eficit medios anuales, respectivamente, calculados con el modelo de balance h´ıdrico para el per´ıodo 2081/2090 - Escenario A2. Se observa una dr´astica disminuci´on de los excesos h´ıdricos en la regi´on. Si bien los excesos medios anuales calculados conservan un gradiente de aumento hacia el Noreste, similar al per´ıodo de referencia, el ´area con excesos se reduce y los montos anuales son significativamente menores a los de aqu´el per´ıodo. El exceso anual m´aximo es de 400 mm/a˜no y se ubica en la Prov. de Misiones, mientras que los sectores al Oeste y Centro de las provincias de Formosa, Chaco y Santa Fe y sector Este de Entre R´ıos tienen un exceso nulo. Efectos del cambio clim´atico sobre la disponibilidad de recursos h´ıdricos en la regi´on litoralmesopotamia (Argentina) 111 Figura 8. Excesos medios anuales calculados para el per´ıodo 2081/2090 - A2 [mm] Figura 9. D´eficit medios anuales calculados para el per´ıodo 2081/2090 - A2 [mm] Figura 10. Variaci´on entre los excesos medios anuales del per´ıodo 2081/2090 - A2 y del per´ıodo 1981/1990 [mm] Tambi´en se observa que pr´acticamente toda la regi´on pasa a tener d´eficit h´ıdrico, con un gradiente de aumento hacia el Noroeste. A diferencia de la situaci´on actual, donde los sectores con d´eficit h´ıdricos se limitan a la zona Oeste de las provincias de Formosa, Chaco y Santa Fe, para el escenario A2 a´un los sectores ubicados en el Centro y Este de la regi´on pasan a tener d´eficits estacionales, durante meses de primaveraverano. La Figura 10 muestra las isol´ıneas de variaci´on (negativa) entre los excesos medios anuales del per´ıodo 2081/2090 - Escenario A2 y del per´ıodo de referencia. La Figura 11 muestra las isol´ıneas de variaci´on (positiva) de d´eficit medios anuales entre el per´ıodo 2081/2090 - Escenario A2 y el per´ıodo de referencia. La disminuci´on de excesos y el aumento de d´eficit se producir´ıan fundamentalmente por aumento de la evapotranspiraci´on, debido al aumento de la temperatura del aire. Los cambios implicar´ıan una reducci´on tanto de la escorrent´ıa directa como de la recarga subterr´anea. Los cursos aut´octonos, cuya escorrent´ıa se genera por lluvias ca´ıdas dentro de la regi´on en estudio, ver´ıan reducidos sus caudales medios, lo que implicar´ıa un cambio en la ten-