texto en 
Inglés (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
Introducción
Para una planificación y gestión adecuadas de los recursos hídricos se debe tener en cuenta la incertidumbre asociada con la variabilidad de la oferta y demanda de los recursos disponibles en una región. En ese aspecto, es clave conocer la evolución en tiempo y espacio de los recursos hídricos disponibles para la toma de decisiones e identificación de la fragilidad de los sistemas de una misma región.
Las situaciones de extremos hidrológicos, como periodos de excesos o escasez de agua, representan amenazas para la sociedad. En particular, las sequías, que son periodos de escasez, son fenómenos de relativa complejidad, que afectan de manera importante el desarrollo y aprovechamiento de los recursos hídricos en una región (Fernández, 1997). Avanzar en el conocimiento de sus características (frecuencia de ocurrencia, duración, extensión espacial, intensidad y magnitud) y su relación con los fenómenos macroclimáticos que intervienen en su formación, ayuda a la generación de modelos de predicción a mediano y largo plazos, que permitirían estar preparados para la gestión eficiente de los recursos hídricos durante los periodos de escasez.
En un trabajo previo realizado por los autores, se abordó la identificación y caracterización de las sequías en la República Argentina desde una concepción hidrológica (Díaz, Rodríguez, Dölling, Bertoni, & Smrekar, 2016). En este trabajo se identifican las periodicidades plurianuales tanto de las series de caudales de los sistemas fluviales de la República Argentina como de indicadores macroclimáticos, utilizando un análisis espectral.
Conocer las frecuencias de tiempo dominantes (caracterizando las periodicidades plurianuales) en las series de caudales permite comprender la evolución temporal de los procesos que intervienen en los ciclos hidrológicos y así avanzar en la generación de modelos conceptuales que expliquen dichos procesos, al igual que desarrollar modelos de pronóstico que ayuden en distintas escalas a prever los años de déficit hídricos, así como en la planificación y gestión de los recursos hídricos.
Región de estudio
La región de estudio incluye 14 cuencas de las regiones Centro, Norte y Cuyo de la República Argentina. En la Tabla 1 y la Figura 1 se describen sus principales características.
Tabla 1 Resumen de las características hidrológicas y geográficas de las cuencas analizadas.
Figura 1 Ubicación de las cuencas analizadas en el trabajo. Los números corresponden con el indicador mostrado en la Tabla 1.
Las series de caudales analizadas corresponden a cuencas de la región Centro, Norte, Noa y Cuyo de la República Argentina mostradas en la Figura 1. Los ríos presentan características diferentes en cuanto a localización de la cuenca, módulo, área de aporte y volumen de aporte anual (ver Tabla 1).
La cuenca del río San Juan cubre gran parte de la provincia de San Juan y abarca un pequeño sector del norte de Mendoza. Casi la totalidad del escurrimiento superficial proviene del derretimiento de la nieve o hielo de los glaciares en los Andes, y constituye la principal fuente de recarga de los acuíferos.
La cuenca del río Mendoza se ubica en el extremo noroeste de la provincia homónima y cubre una pequeña extensión del sur de la provincia de San Juan. Esta cuenca drena 90 km del frente de la cordillera de los Andes. Las aguas de este río provienen en su mayoría del deshielo. El río Atuel es el quinto afluente del desaguadero; es de régimen nival, aunque también recibe aportes pluviales, y tiene una longitud aproximada de 600 km. Su cuenca de alimentación está a más de 3 000 msnm en la cordillera de los Andes.
La cuenca del río Colorado (Subsecretaría de Recursos Hídricos, 2010) comprende cuatro ecorregiones con relieve y régimen de precipitaciones muy variado, desde 100 hasta 600 mm medios anuales. Es una cuenca pluvionival.
La cuenca del río Ctalamochita pertenece a la cuenca Carcaraña y es de carácter exorreico, terminando en la cuenca del Plata. Nace en las vertientes orientales de la Sierra Grande. En la actualidad, el río Tercero (Ctalamochita) está regulado por una cadena de reservorios artificiales.
El río Xanaes nace por la confluencia de los ríos Anisacate y de Los Molinos. El río Anisacate resulta de la unión de los ríos de la Suela y el San José. Atraviesa la Sierra Chica en una garganta estrecha como río retrogradante. El río Suquía posee una cuenca de unos 6 000 km2 cerca de Mar Chiquita y adquiere su nombre en el embalse San Roque. En el punto de estudio, el área de la cuenca es de 1 350 km2. Estos últimos son típicos ríos de montaña, donde las variaciones de caudal responden directamente a las variaciones de las precipitaciones en la cuenca. Los caudales por arriba del módulo suceden entre octubre y marzo; mientras que los estiajes son más acentuados y abarcan el periodo junio-septiembre.
La cuenca de aporte del río Dulce en el tramo estudiado abarca las provincias de Tucumán, Salta y Catamarca. Presenta una precipitación promedio anual de 800 mm y se considera de clima húmedo.
La cuenca del río Pasaje-Juramento tiene sus nacientes en las cumbres de los nevados de Cachi y Acay (4 895 msnm), en la cordillera oriental. La cuenca alta y media se desarrolla dentro de la provincia de Salta. El río Juramento está regulado a través del dique Embalse General Belgrano (Cabra Corral), de propósito múltiple. La cuenca del río Bermejo forma parte de la hoya hidrográfica del Plata y junto con el río Pilcomayo es el principal afluente del río Paraguay. La cuenca alta del río Bermejo se sitúa en el extremo noroeste de Argentina y extremo sursudeste de Bolivia. La precipitación anual media varía desde 200 mm en la parte oeste hasta 1 400 mm en el centro de la cuenca. El Pilcomayo es un río de montaña que tiene sus nacientes en los Andes bolivianos, a más de 5 000 m de alturas snm; pasa en dirección NO-SE por las sierras subandinas.
A lo largo de los tiempos geológicos, el río fue depositando gran parte de los sedimentos que transporta sobre la planicie chaqueña, construyendo así un gran abanico aluvial. La cuenca del río Paraná en el punto de análisis (estación Corrientes) tiene un área de 1 950 000 km2. Esta cuenca es parte de la cuenca del Plata, una de los mayores del mundo. Se nota una escasez pluvial en invierno (junio-agosto) en la mayor parte de la cuenca, con lluvias intensas en verano, siendo más reducidas en el oeste y en las regiones situadas al norte del trópico de capricornio, hasta el Planalto de Brasil; hacia el sur se caracteriza por abundantes precipitaciones. La precipitación media anual es de 1 200 mm (1961-1990) (Paoli & Schreider, 2000). El río Salado pertenece al complejo hídrico de la cuenca del Plata. Las fuentes del río Salado surgen en las estribaciones orientales de los Andes dentro de la provincia de Salta; su fuente principal se encuentra en la llamada Sierra de los Pasos Grandes, casi inmediatamente al sur del cerro nevado de Acay. Es un río de régimen pluvionival.
Metodología y resultados
Análisis espectral en las series de caudales medios anuales
El análisis espectral consistió en calcular, en primera instancia, los espectros de energía de las series históricas de caudales medios anuales, con el objetivo de identificar las frecuencias dominantes. En el análisis espectral se analiza la serie de caudales en el dominio de la frecuencia, utilizando la transformada rápida de Fourier. Después, el espectro de energía de la fluctuación se calculó con la siguiente ecuación:
Donde:
A continuación (Figura 2, Figura 3, Figura 4, Figura 5, Figura 6, Figura 7, Figura 8, Figura 9, Figura 10, Figura 11, Figura 12, Figura 13, Figura 14 y Figura 15) se presentan las series de caudales medios anuales con sus respectivos espectros adimensionales de energía de las fluctuaciones de los caudales (utilizando el valor máximo para su adimensionalización), resaltando las tres frecuencias consideradas dominantes (máximo valor) en cada espectro (Tabla 2).
Figura 2 Serie de caudales del río San Juan y espectro de energía de las fluctuaciones de los caudales.
Figura 3 Serie de caudales del río Atuel y espectro de energía de las fluctuaciones de los caudales.
Figura 4 Serie de caudales del río Dulce y espectro de energía de las fluctuaciones de los caudales.
Figura 5 Serie de caudales del río Suquía y espectro de energía de las fluctuaciones de los caudales.
Figura 6 Serie de caudales del río Mendoza y espectro de energía de las fluctuaciones de los caudales.
Figura 7 Serie de caudales del río Ctalanochita y espectro de energía de las fluctuaciones de los caudales.
Figura 8 Serie de caudales del río Anizacate y espectro de energía de las fluctuaciones de los caudales.
Figura 9 Serie de caudales del río Juramento y espectro de energía de las fluctuaciones de los caudales.
Figura 10 Serie de caudales del río Pilcomayo y espectro de energía de las fluctuaciones de los caudales.
Figura 11 Serie de caudales del río Salado y espectro de energía de las fluctuaciones de los caudales.
Figura 12 Serie de caudales del río Bermejo y espectro de energía de las fluctuaciones de los caudales.
Figura 13 Serie de caudales del río Xanaes y espectro de energía de las fluctuaciones de los caudales.
Figura 14 Serie de caudales del río Colorado y espectro de energía de las fluctuaciones de los caudales.
Figura 15 Serie de caudales del río Paraná y espectro de energía de las fluctuaciones de los caudales.
Tabla 2 Frecuencias dominantes y energía del espectro (adimensionalizada con su valor máximo) de cada cuenca estudiada.
En las figuras 6, 11, 12 y 15 se marcó con un rectángulo la primera frecuencia dominante del espectro. Estas frecuencias (mayores al largo de la serie) no se tuvieron en cuenta porque con el análisis de series sintéticas se demostró que son producto de tendencias en la serie o componentes de la baja frecuencia del largo de la serie.
Se puede observar en la Tabla 3 que los dos periodos dominantes se dan en los siguientes rangos:
< 3.7 años 2 casos
3.7 - 6.7 años 9 casos
8.5 - 10.7 años 4 casos
10.7 - 18.3 años 2 casos
18.3 - 32 años 11 casos
Tabla 3 Frecuencias dominantes de algunas de las cuencas estudiadas.
| Cuencas | Fourier | Wavelet | ||
|---|---|---|---|---|
| Periodo (años) | Periodo (años) | Periodo (años) | Periodo (años) | |
| Atuel | 4.7 | 10.7 | 4.1 | 23.4 |
| Ctalamochita | 4.3 | 25.6 | 4.1 | 23.4 |
| Anizacate | 5.3 | 16 | 5.8 | 16.5 |
| Suquía | 18.3 | 6.7 | 16.5 | 4.1 |
| Dulce | 21.3 | 8.5 | 8.3 | 23.4 |
| Paraná | 8.5 | 32 | 8.3 | 33.1 |
Se observa que para la mayoría de los ríos (excepto Colorado, Paraná y Dulce) se presentan fluctuaciones con periodos dominante entre 3 y 7 años. Periodos dominantes entre 7 y 11 años se observan en los ríos Colorado, Paraná, Atuel y Dulce; mientras que periodos dominantes entre 13 y 35 años se ve en todos los ríos, excepto Atuel, Salado y Pilcomayo. Existen resultados similares en trabajos de otros autores, como Vargas, Minetti y Poblete (2002), quienes detectaron fluctuaciones de baja frecuencia cuasi periódicas (de 22 a 26 años) en los caudales anuales del río Paraná; Compagnucci, Berman, Velasco-Herrera y Silvestri (2014) encontraron en la serie de caudales del río Paraná periodicidades de 3.5, 9 y 30 años, y en la serie de caudales del río Atuel, periodicidades de 4 a 5 años, 7, 11 y 22 años.
En la serie de caudales del río San Juan se apreciaron cinco bandas de frecuencia con amplitudes significativas, que se corresponden con periodos de unos 12.3, 7.4, 5.7 y 3.8 años (Correa & Guevara, 1992).
Análisis de coherencia Wavelet
Además de efectuar el análisis espectral con Fourier, se realizó un análisis de coherencia Wavelet en las series de caudales para validar los resultados obtenidos. Para las señales que presentan múltiples escalas de variabilidad y frecuencia en el espacio o tiempo puede ser de utilidad un análisis localizado de transformaciones con curvas tipo wavelets, a fin de discriminar las geometrías dominantes de las variables analizadas. La coincidencia que muestran los resultados de frecuencias dominantes obtenidos con las metodologías de análisis espectral con base en la transformada de Fourier y curvas tipo Wavelets son una muestra de que las señales analizadas no tienen un alto grado de variabilidad de frecuencias dominantes y, por lo tanto, con la aplicación del análisis espectral con base en Fourier se alcanzarían resultados satisfactorios.
Se realizó el análisis con la función Wavelet tipo Morlet (ko = 6).
Análisis de la contribución de distintos procesos a las fluctuaciones de las series de caudales
Después de realizar el análisis de frecuencias dominantes de las distintas series de caudales y encontrar que existen rangos de frecuencias que predominan en las series, se procedió a realizar un segundo análisis, que es el cálculo de la varianza. Dicho cálculo se llevó a cabo integrando el espectro de energía de las fluctuaciones en el rango de frecuencias analizado (Figura 16). Esto se realiza con el fin de evaluar el aporte porcentual que tienen a la varianza total y las varianzas de las fluctuaciones de cada rango de frecuencia característicos de los distintos procesos. Para efectuar esto, se procedió a aplicar filtros pasa-banda a los espectros entre los siguientes rangos de frecuencia:
El filtro pasa banda utilizado consiste en un filtro de Fourier (en primer lugar se verificó la herramienta con el análisis de una serie sintética, fluctuaciones de periodos conocidos), haciendo que la energía de las fluctuaciones sean iguales a cero para los rangos de frecuencias no incluidos en la banda del filtro pasa-banda utilizado. La Tabla 4 resume los resultados, donde se indica el porcentaje de aporte y la varianza que tienen las fluctuaciones de cada rango de frecuencia.
Tabla 4 Aporte a la varianza total de las fluctuaciones de las series de caudales en cada ancho de banda.
| Cuenca | Varianza total de las fluctuaciones de las series de caudales en cada ancho de banda | ||||||
|---|---|---|---|---|---|---|---|
| N a | N 2 | N 5 | N 11 | N bd | N vc | Suma | |
| San Juan | 45% | 3% | 27% | 5% | 4% | 7% | 87 |
| Mendoza | 40% | 0% | 9% | 3% | 20% | 3% | 75 |
| Atuel | 36% | 2% | 23% | 9% | 15% | 11% | 96 |
| Colorado | 50% | 8% | 12% | 11% | 10% | 1% | 92 |
| Ctalamochita | 50% | 0% | 19% | 8% | 19% | 0% | 96 |
| Xanaes | 43% | 0% | 9% | 8% | 15% | 9% | 84 |
| Anizacate | 19% | 2% | 38% | 12% | 5% | 0% | 76 |
| Suquía | 38% | 0% | 14% | 8% | 15% | 15% | 90 |
| Dulce | 28% | 6% | 4% | 16% | 25% | 16% | 95 |
| Juramento | 49% | 12% | 22% | 1% | 12% | 0% | 96 |
| Bermejo | 27% | 4% | 19% | 3% | 23% | 20% | 96 |
| Pilcomayo | 33% | 6% | 9% | 4% | 26% | 0% | 78 |
| Paraná | 31% | 3% | 16% | 21% | 9% | 16% | 96 |
| Salado | 30% | 5% | 36% | 2% | 2% | 2% | 77 |
La suma distinta al 100% se debe a la correlación entre las componentes.
En N a se consideran los procesos menores a 2.2 años y la parte aleatoria de todas las frecuencias. Correa y Guevara (1992) proponen que la varianza residual no explicada por las periodicidades significativas debe ser atribuida a la componente estocástica. Esta varianza aleatoria para 86% de los casos representa entre 30 y 50% de la varianza de las fluctuaciones, lo cual significa que la serie tiene un alto componente aleatorio. El resto puede explicarse en parte por distintos procesos con frecuencias características entre 2 y 3 años, 3 y 7 años, 7 y 13 años, 13 y 35 años y mayores a 35 años.
Se resalta en la Tabla 4 (sin considerar los procesos N a ) que para nueve cuencas de 14 estudiadas, los procesos de baja frecuencia dominan N bd y luego lo siguen en importancia los N 5. Para dos cuencas se resalta el proceso N v y para sólo una, el proceso N 11 (Paraná). Así, el análisis reveló que existe en las series históricas de caudales medios anuales escurridos fluctuaciones de alta, baja y media frecuencia, que aportan, con distinta significancia, a la variabilidad temporal de los caudales escurridos. Este resultado permite avanzar en el estudio de los procesos que se deben considerar en la explicación de fenómenos hidrológicos, como sequías y excesos hídricos, teniendo en cuenta que el proceso involucra fenómenos que tienen en el tiempo una evolución que incluye procesos a pequeña, media y gran escala, y aleatorios o estocásticos.
Relaciones entre la evolución temporal de sequías hidrológicas e indicadores macroclimáticos
En esta sección se analizan indicadores macroclimáticos y astronómicos, utilizando el mismo análisis espectral de las series temporales. Para conocer si existe una relación entre su evolución temporal con las sequías hidrológicas en las distintas escalas de tiempo.
En la Tabla 5 se describen los indicadores macroclimáticos analizados.
Tabla 5 Resumen de los indicadores macroclimáticos utilizados.
| Indicador | Descripción | Periodo |
|---|---|---|
| TNA | Índice del Atlántico norte tropical | 1948-2013 |
| TSA | Índice del Atlántico sur tropical | 1948-2013 |
| AMO | Oscilación multidecadal del Atlántico | 1861-2008 |
| AMM | Modo meridional del Atlántico | 1948-2001 |
| SOI | Índice de oscilación sur | 1951-2013 |
| PDO | Oscilación decadal del Pacífico | 1948-2013 |
| MS | Manchas solares | 1700-2013 |
| ONI | Índice de El Niño oceánico | 1950-2012 |
| Niño 3.4 | Temperatura promedio de la superficie del mar en las regiones El Niño 3 y El Niño 4 | 1950-2013 |
| Niño 1+2 | Temperatura promedio de la superficie del mar en las regiones 2+1 | 1950-2013 |
Los indicadores seleccionados son ONI, AMM, TSA, AMO, SOI, PDO, Niño 1+2, Niño 3.4 y TNA (ver Tabla 5). En la Tabla 5 y Tabla 6 se indican los periodos dominantes observados de las fluctuaciones y los aportes a la varianza (en porcentaje) de las fluctuaciones con los distintos anchos de bandas analizados anteriormente.
Las series de indicadores AMO, AMM, ONI, PDO, SOI, Niño 3.4, Niño 1+2, TSA y TNA se obtuvieron del sitio web del NOAA (2016), mientras que la serie de manchas solares se obtuvo de SILSO (2016).
En los índices PDO, TNA, AMM se encuentra un periodo de frecuencia dominante de 64 años-1 y en el índice TSA de 128 años-1; sin embargo, se descartan por presentar mucha incertidumbre, pues las series muestran longitudes inferiores a cien años. De este análisis es posible ver que los indicadores Niño 1+2, Niño 3.4 y ONI tienen una frecuencia dominante del orden de 4 a 5 años. Los indicadores manchas solares, AMM, SOI, TSA y TNA tienen frecuencias dominantes entre 8.5 y 13 años. Mientras que los indicadores PDO y AMO presentan una frecuencia dominante multidecadal.
En el análisis de varianza que se indica en la Tabla 7 se observa que los procesos N 11 aportan un 66% a la varianza del espectro de las fluctuaciones. Mientras que los procesos N 5 tiene un importante aporte en los indicadores ONI, Niño 3.4 y Niño 1+2.
Tabla 7 Aporte de la varianza en porcentaje a los distintos anchos de bandas analizados.
| MS | AMO | AMM | ONI | PDO | SOI | Niño 3.4 | Niño 1+2 | TSA | TNA | |
|---|---|---|---|---|---|---|---|---|---|---|
| N a | 2% | 1% | 56% | 52% | 21% | 44% | 49% | 64% | 43% | 52% |
| N 5 | 4% | 0% | 0% | 40% | 19% | 25% | 38% | 34% | 5% | 0% |
| N 11 | 66% | 1% | 18% | 6% | 14% | 18% | 6% | 0% | 12% | 14% |
| N bd | 5% | 10% | 1% | 0% | 9% | 2% | 0% | 0% | 7% | 5% |
| N v | 20% | 86% | 3% | 0% | 29% | 0% | 0% | 0% | 25% | 21% |
Análisis de correlaciones entre las fluctuaciones de indicadores macroclimáticos y series de caudales medios anuales
Conocer los indicadores macroclimáticos cuyas fluctuaciones puedan tener vinculación con las sequías en las distintas escalas de tiempo es una herramienta que ayuda a la predicción de tales fenómenos. En este apartado se analizan las correlaciones entre las fluctuaciones de indicadores macroclimáticos y de los caudales medios anuales, para distintos anchos de bandas. Las cuencas seleccionadas para este análisis son la de los ríos Suquía, Dulce, Paraná, San Juan y Atuel, por ser las series más extensas y representativas de cada región estudiada. Los filtros pasa-banda de Fourier utilizados fueron los siguientes: N bd : filtro pasa-banda entre 13 y 35 años; N 11: filtro pasa-banda entre 7-13 años; N 5: filtro pasa-banda entre 3 y 7 años. En este análisis se calculó el coeficiente de correlación entre las series temporales de cada serie filtrada de caudales de cada cuenca (N bd , N 11, N 5) y las series filtradas de los índices macroclimáticos en el mismo ancho de banda. Luego se graficaron las series que poseen mejor correlación. A continuación se presentan los resultados para los ríos analizados.
Río Dulce
En la Tabla 8 se muestra el coeficiente de correlación entre las series de caudales del río Dulce y las series de los distintos indicadores.
Tabla 8 Coeficiente de correlación entre las series de caudales del río Dulce y las series de los distintos indicadores.
A continuación, en las siguientes figuras, se muestran la evolución temporal de caudales y los indicadores macroclimáticos. Para graficar las series de caudales, éstas se normalizan restando la media y dividiendo por el desvío de cada serie.
En la Figura 17 se observa que las fases negativas del indicador AMM coincidieron con las sequías severas de la cuenca de mediados de la década de los años 40, 60 y fines de los 80, y la sequía de los años 2011-2013. Se destaca la alta correlación encontrada entre la serie de caudales y los indicadores MS, AMO, PDO, Niño 1+2 y TSA para el ancho de banda de 13 a 35 años. En la Figura 17b se observa que el indicador Niño 1+2 en su fase negativa sucedió luego de dos años, por periodos de sequías en la cuenca. Este comportamiento se observa para las sequías de mediados de la década de los años 50 y 60 (sequías severas por su permanencia en el tiempo), y para las sequías de los años 1988-1989, 1993-1995, 2002-2004 y 2011-2013. En la Figura 17e se observa una correspondencia entre valles y picos para la mayoría de los ciclos observados, excepto para el periodo 1994-2004. En la Figura 18 se resalta que las fases positivas del TSA coinciden con la sequía de 1962-1971, y anteceden las sequías 1988-1990 y 2005-2010. Cabe señalar que desde 1950, estas sequías son las más largas que se identificaron a la fecha.
Figura 17 a) Serie de caudales normalizados del río Dulce y del indicador AMM ambas filtradas con filtro pasa-banda entre 13 y 35 años. Caudales normalizados del río Dulce con filtro pasa-banda entre 7 y 13 años, e indicadores Niño 1+2 b) ONI c) y d) TNA; e) serie de caudales normalizados del río Dulce y del indicador TNA, ambas filtradas con filtro pasa-banda entre 3 y 7 años.
Figura 18 a) Serie de caudales del río Suquía y del indicador manchas solares; b) TSA, todos filtrados con filtro pasa-banda entre 13 y 35 años; c) caudales del río Suquía e indicador TSA; d) PDO, con filtro pasa-banda entre 7 y 13 años; e) serie de caudales del río Suquía y del indicador Niño 3.4, ambos filtrados con filtro pasa-banda entre 3 y 7 años.
En la Figura 19a se observa que la fase positiva del indicador TSA coincide con déficits hídricos (sequía 1966-1972), y periodos de caudales bajos (no identificados como sequías hidrológicas) en los años 1985-1988 y 1999-2008. En la Figura 19c se observa que las fases negativas de indicador Niño 1+2 coinciden con las sequías 1947-55 (en la etapa final), 1961-1963, 1977, 1995, y 2001-2005, y con la etapa de caudales bajos en 1986-1988 no identificadas como sequías hidrológicas. En la Figura 20a puede observarse que las fases negativas de indicador AMM concuerdan con los periodos de déficits hídricos en la cuenca, coincidiendo con las sequías identificadas en 1945-1952, 1966-1971, 1988-1996 y 2009-2013. Las fases negativas del indicador Niño 3.4 coinciden con los eventos de sequía de la cuenca en los años 1950-1951, 1960, 1962-1965, 1975, caudal bajo de 1985, 1988-1989, 1995-1996, 1999, 2004, 2007 (ver Figura 20c). En la Figura 21a se observa que las fases negativas del indicador AMM coinciden con los ciclos bidecadales de caudales bajos en la cuenca del río San Juan. Dentro de estas fases negativas se identificaron las sequías de 1945-53, 1966-1975, 1988-1996 y 2010-2013. En la Figura 21d se observa que las fases negativas del Niño 3.4 coinciden con la mayoría de los eventos de sequía identificados. Estos eventos son 1954-1957, 1960-1962, 1964, 1974-1975, 1979-1981, 1985, 1988-1989, 1993-1996, 1999, 2007-2008.
Figura 19 a) Serie de caudales del río Paraná y del indicador TSA; b) PDO, todos filtrados con filtro pasa-banda entre 13 y 35 años; c) serie de caudales del río Paraná y del indicador Niño 1+2, ambos filtrados con filtro pasa-banda entre 7 y 13 años; d) serie de caudales del río Paraná y del indicador ONI, y e) SOI, todos filtrados con filtro pasa-banda entre 3 y 7 años.
Figura 20 a) Serie de caudales del río Atuel y del indicador AMM, ambos filtrados con filtro pasa-banda entre 13 y 35 años; b) serie de caudales del río Atuel y del indicador SOI, ambos filtrados con filtro pasa-banda entre 7 y 13 años; c) serie de caudales del río Atuel y del indicador Niño 3.4, ambos filtrados con filtro pasa-banda entre 3 y 7 años.
Figura 21 a) Serie de caudales del río San Juan y del indicador AMM, ambos filtrados con filtro pasa-banda entre 13 y 35 años; b) serie de caudales del río San Juan y del indicador manchas solares, ambos filtrados con filtro pasa-banda entre 7 y 13 años; c) serie de caudales del río San Juan y del indicador SOI, ambos filtrados con filtro pasa-banda entre 3 y 7 años; d) serie de caudales del río San Juan y del indicador Niño 3.4, ambos con filtro pasa-banda entre 3 y 7 años.
Río Suquía
En la Tabla 9 se observa el coeficiente de correlación entre las series de caudales del río Suquía y las series de los distintos indicadores.
Tabla 9 Coeficiente de correlación entre las series de caudales del río Suquía y las series de los distintos indicadores.
Río Paraná
En la Tabla 10 se tiene el coeficiente de correlación entre las series de caudales del río Paraná y las series de los distintos indicadores.
Tabla 10 Coeficiente de correlación entre las series de caudales del río Paraná y las series de los distintos indicadores.
En la Figura 19 se presentan los más significativos.
Río Atuel
En la Tabla 11 se observa el coeficiente de correlación entre las series de caudales del río Atuel y las series de los distintos indicadores.
Tabla 11 Coeficiente de correlación entre las series de caudales del río Atuel y las series de los distintos indicadores.
En la Figura 20 y Figura 21 se presentan los más significativos.
Río San Juan
En la Tabla 12 se muestra el coeficiente de correlación entre las series de caudales del río San Juan y las series de los distintos indicadores.
Tabla 12 Coeficiente de correlación entre las series de caudales del río San Juan y las series de los distintos indicadores.
Conclusiones
En este trabajo se identificaron las periodicidades plurianuales en series de caudales escurridos en sistemas fluviales de la República Argentina. El análisis espectral realizado permite determinar que para 86% de las series de caudales analizadas, la componente aleatoria (alta frecuencia) explica entre 30 y 50% de las señales de caudales. Además se observó que existe una variabilidad en las series de alta, baja y media frecuencia que aportan a la fluctuación de los caudales en distintos porcentajes. Luego de la componente aleatoria, los procesos con predominio son los que tienen periodos dominantes entre 3 y 7 años (p. ej., ENSO), y los que tienen periodos entre 13 y 35 (oscilaciones bidecadales). Los resultados obtenidos en este trabajo permiten avanzar en el estudio de los procesos que se deben considerar en la explicación de fenómenos hidrológicos, como sequías y exceso de agua.
A nivel bidecadal se encontró una gran influencia del océano Atlántico sobre las series de caudales. En las cinco cuencas analizadas (ríos Suquía, Dulce, Paraná, Atuel y San Juan) se observó una alta correlación con el índice Atlantic Meridional Mode “AMM”, Dulce (0,82); Paraná (0,53); Atuel (0,87); San Juan (0,81) y Suquía (0,51). Siendo la más alta correlación encontrada en casi todas las cuencas, excepto para el Paraná y Suquía, que presentan una alta correlación negativa con el TSA Paraná (-0.69) y Suquía (-0.74). A nivel decadal se observan correlaciones significativas (mayores a 0.5) para el Paraná con el Niño 1+2 (0.68), el Atuel con el SOI (-0.62) y Suquía con TSA (-0.53).
En la escala plurianual se observan correlaciones negativas con el SOI, para las cuencas de los ríos Paraná (-0.59), Atuel (-0.50) y San Juan (-0.72). Fuertes correlaciones positivas con los indicadores ONI para el Paraná (0.60); San Juan y Atuel, con el Niño 3.4 (0.70; 0.46 respectivamente), y el Dulce con el TNA (0.54).
Esto no sólo confirma los resultados de trabajos anteriores (Compagnucci et al., 2014; Mauas, Buccino y Flamenco, 2011; Dölling, 2014; Vargas et al., 2002), además aporta nueva información con respecto a vínculos con indicadores del océano Atlántico para las distintas cuencas. Estos resultados son importantes para el entendimiento de la dinámica macroclimática de Argentina y de su pronóstico en el espacio y en el tiempo.
