Introducción

Para calcular el inicio y el término de la época lluviosa se han utilizado numerosas técnicas y diferentes variables meteorológicas. Por ejemplo, Zuo et al. ⁽¹⁾ utilizaron patrones predecibles de la precipitación monzónica en el verano del Hemisfério Norte, comparando los modos principales de anomalías de precipitación observada en los monzones y luego los compararon con los del pronóstico actualizado del modelo Climate Forecast System Version 2 ⁽²⁾ encontraron correlaciones significativas con la variabilidad de la temperatura superficial del mar (TSM) del Atlántico y Pacífico tropical. Enfield y Alfaro ⁽³⁾ utilizan dos enfoques independientes para el análisis de las relaciones interoceánicas con la lluvia en América Central. El primero es un análisis de descomposición de valores singulares (DVS), aplicado a conjuntos de datos de rejilla de TSM y precipitación acumulada (PCPA) mensual para detectar diversas configuraciones de la TSM del Pacífico y del Atlántico y sus asociaciones con la precipitación, donde encontraron modos que relacionan la TSM y la PCPA. El segundo, la relación de los índices de TSM con las fechas de inicio de la estación lluviosa (IELL) y término de la estación lluviosa (TELL) en América Central donde uno de sus resultados más importante es la relación que un Atlántico Tropical Norte (ATN) cálido (frío) favorece una expansión (contracción) de la temporada de lluvias en ambos extremos en el sur de América Central.

En la región Mesoamericana, Díaz-Alcántara ⁽⁴⁾ propone que la actividad convectiva de la Zona de Convergencia Intertropical (ZCIT) induce circulaciones relacionadas con ondas ecuatoriales que llevan humedad y provocan las lluvias. Utiliza datos de un modelo numérico y de estaciones meteorológicas de la costa del Pacífico Mesoamericano, como criterio para el cálculo de la fecha del IELL, realiza pendientes de PCPA utilizando el método de mínimos cuadrados con una ventana de diez días. Para realizar el cálculo de la fecha del IELL usó datos del modelo Global Precipitation Climatology Project 1-Degree Daily Combination 1 ⁽⁵⁾, y define como el IELL al cambio de pendiente en la PCPA para cada año.

Alfaro et al. ⁽⁶⁾ utilizaron las distribuciones de la TSM en los océanos Pacífico y Atlántico para explicar las variaciones de las fechas del IELL y el TELL sobre América Central, donde sus resultados muestran correlaciones significativas entre el índice del ATN con el IELL, el Índice de Oscilación del Sur (IOS) con el TELL y el Niño 3 también con el TELL. Alfaro y Cid ⁽⁶⁾, usaron datos de precipitación pentadal (precipitación acumulada durante 5 días), utilizando una versión modificada de la metodología propuesta por ⁽⁸⁾. Por otro lado, Hernández y Fernández ⁽⁹⁾ relacionan datos de precipitación y evapotranspiración diaria para obtener umbrales del inicio y término de la estación lluviosa utilizando el método FAO Penman-Monteith o FAO-56 de Allen et al. ⁽¹⁰⁾. El método del análisis del inicio y término de la época lluviosa (ITELL) utilizando péntadas de precipitación acumulada también fue utilizado por Nakaegawa et al. ⁽¹¹⁾ en Panamá, y Rivera-Ramos et al. ⁽¹²⁾ en la zona sur de Guatemala, los últimos con las mismas pentadas aplicaron un segundo criterio que consiste en la diferencia entre el valor promedio anual de precipitación pluvial y el valor de precipitación de cada pentada propuesto por Liebmann y Marengo ⁽¹³⁾.

En la presente investigación se utilizará un método elaborado para capturar un cambio estacional en el régimen de precipitación en América Central adaptado para ser usado con modelos de pronóstico y análisis retrospectivo. Este método, que solo usa datos de precipitación, fue diseñado por Bombardi y Carvalho ⁽¹⁴⁾, basado en el método propuesto por Liebmann y Marengo ⁽¹³⁾. El objetivo de este trabajo es comparar dos metodologías para calcular el inicio y término de la estación lluviosa de estaciones meteorológicas con ciclo anual de precipitación bimodal en la región centroamericana ⁽¹⁵⁾ utilizando primero el análisis de las péntadas de precipitación diaria acumulada y segundo el cálculo de la desviación de la precipitación acumulada (DPA), ambas metodologías sujetas a diferentes criterios de detección.

Materiales y métodos

Se utilizan datos diarios de precipitación acumulada de 1968 al 2012 de las 101 estaciones mostradas en la figura 1. Estas estaciones en su mayoría son de tipo bimodal con una estación lluviosa y seca muy marcada ⁽¹⁵⁾, escogidas según el criterio de Alfaro e Hidalgo ⁽¹⁶⁾. Para tener series de tiempo completas se rellenan los datos con la metodología de rellenado de datos geofísicos expuesta por Alfaro y Soley ⁽¹⁷⁾ donde como primera aproximación se usa un rellenado por medio de un filtro predictivo AR(p) y luego utilizando el método de componentes principales ⁽¹⁸⁾. Para estos datos se utilizaron 5 componentes principales.

Figura 1 Mapa de la distribución de estaciones utilizadas en América Central. 

Cálculo de las fechas de inicio y final de la época lluviosa

Primero, para el análisis de las péntadas de precipitación diaria acumulada, basado en la metodología desarrollada en Alfaro ⁽¹⁹⁾, y también aplicada en ⁽¹⁶, ²⁰, ²¹, ²²⁾, las series de tiempo de precipitación diaria se filtraron utilizando un promedio móvil de peso triangular con una ventana de 31 días ⁽²³⁾, para evitar o minimizar las interrupciones del inicio debido a posibles fortalecimientos de corta duración de los vientos alisios. Se aplicó un algoritmo para identificar sistemáticamente el inicio a la serie de tiempo de precipitación diaria filtrada. La metodología de Gramzow y Henry ⁽¹⁴⁾, que fue posteriormente adaptada en ⁽³, ⁶, ⁷, ²⁴⁾, consiste en que para un año en particular y en un registro de estación en particular, el inicio de las lluvias (en días julianos) se asocia con aquel día en que la precipitación suavizada es mayor a 5 mm día−1 durante al menos 10 días (o dos pentadas seguidas) y si en los siguientes 5 días (por ejemplo, la siguiente pentada), la precipitación observada es superior a 0.1 mm día−1. El final de la estación lluviosa es calculado de una manera similar invirtiendo la serie e iniciando del último día del año hacia el inicio.

Segundo, para calcular las fechas de inicio y final de la época lluviosa, basado en el cálculo de la desviación de la precipitación acumulada (DPA), se utilizó para cada año la derivada en el tiempo de la desviación de la precipitación acumulada S(n) (1) del promedio n expuesta por Bombardi et al. ²⁵, donde P(i) es la precipitación diaria el día i, PC es el promedio diario anual de la precipitación (razón climatológica de la precipitación) y una variante de esta. La primera utiliza una media móvil triangular de tres valores para suavizar la curva de S(n). Como criterio de selección desde el inicio de cada año, el primer día donde la derivada pasa de valores negativos a positivos es considerado el inicio de la estación lluviosa, y se mantiene mientras los valores positivos persisten durante tres días.

El final de la estación lluviosa es calculado de una manera similar invirtiendo la serie e iniciando del último día del año hacia el inicio.

La variante propuesta en este trabajo como tercera opción, utiliza S(n), con la diferencia que para suavizar esta curva se aplica un promedio móvil de peso triangular de treinta y un valores. Además, el criterio de selección del inicio de la estación lluviosa se escoge el primer día donde la derivada pasa de valores negativos a positivos, y se mantiene mientras los valores positivos persisten durante diez días. Para el término de la estación lluviosa se invierte la serie y se aplica la metodología anterior. A este criterio se le llamará Bombardi modificado.

Significancia estadística de las diferencias del promedio

Se calcula la duración de la época lluviosa con la diferencia entre los promedios del TELL y el IELL para el conjunto de datos de Alfaro e Hidalgo ⁽¹⁶⁾, la metodología de Bombardi et al. ⁽²⁵⁾ y la variante propuesta en la sección anterior. Se realizan dos pruebas de hipótesis de diferencia de medias utilizando el promedio de los conjuntos de datos, la prueba paramétrica t de dos muestras y la prueba no paramétrica de suma de rangos Wilcoxon ⁽²⁶⁾, entre los datos del IELL, TELL y duración de la época lluviosa (definida como TELL-IELL+1).

Resultados

En la figura 2 se muestra el IELL para las estaciones de la figura 1. Se observa que los inicios tempranos de la ELL se dan desde el 28 de febrero y los inicios tardíos hasta el 28 de junio, entre los mapas se pueden notar similitudes y diferencias entre los resultados. La parte norte de América Central tiene inicios más tardíos que la parte sur (Costa Rica, Panamá), excepto el sur de Honduras que tiene inicios tempranos de la ELL. En las tres metodologías hay una estación en la Península de Azuero de Panamá, en la región conocida como el Arco Seco, con inicio más tardío que las de su alrededor, por otro lado, el sur de Costa Rica y este de Panamá tiene el inicio más temprano de toda la región entre los resultados de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾ y la metodología propuesta.

Se observa que el criterio de Bombardi en Belice y norte de Guatemala no concuerda con los resultados de Alfaro e Hidalgo ⁽¹⁶⁾ y Bombardi modificado, mostrando inicios tempranos cuando en las dos metodologías anteriores muestran inicios tardíos. Para el IELL los datos muestran un rango del inicio desde el 28 de febrero al 28 de junio.

Los diagramas de caja del conjunto de datos de la figura 2 son mostrados en la figura 3. Se observa que los resultados de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾ tienen valores atípicos correspondientes a inicios tempranos de la ELL, para esta serie la mediana se da el 13 de mayo, el primer cuartil el 30 de abril y el tercer cuartil el 20 de mayo. Los resultados de la metodología modificada propuesta presenta un valor atípico correspondiente a un inicio tardío de la ELL, y tienen la mediana el 10 de mayo, el primer cuartil el 4 de mayo y el tercer cuartil el 17 de mayo. Los resultados de la metodología propuesta por Bombardi et al. ⁽²⁵⁾ muestran valores atípicos correspondientes a inicios tempranos de la ELL, una mediana el 19 de abril, el primer cuartil el 5 de abril y el tercer cuartil el 28 de abril.

Figura 2 Promedio del número de día juliano y fechas de inicio de la época lluviosa para las 101 estaciones entre 1968 al 2012 utilizando: (a) resultados de la metodología de Alfaro e Hidalgo ⁽²⁵⁾, (b) resultados de la metodología modificada propuesta, (c) resultados de la metodología utilizada por Bombardi et al. ⁽¹⁶⁾. 

Figura 3 Diagramas de caja del inicio de la ELL para los resultados obtenidos de la figura 2. IELL: resultados de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾, vi_mod: resultados de aplicar la metodología modificada propuesta, vi_Bmb: resultados de aplicar la metodología propuesta por Bombardi et al. ⁽²⁵⁾. Se muestran las fechas del primer cuartil, la mediana y tercer cuartil. 

Para el TELL (figura 4) se puede observar a partir de los resultados de las tres metodologías que, del centro de Costa Rica hacia el sur, junto con Providencia y San Andrés, unas estaciones al norte de Honduras, norte de Guatemala y en el centro y este de Belice, se tiene el término de la estación lluviosa más tardío (más colores verdes). Se observa que, en el centro de Costa Rica hacia el norte, la parte central y Pacífico de Nicaragua, Honduras, El Salvador, y Guatemala, se tiene el término de la estación lluviosa más temprano (más colores cafés).

La metodología de Bombardi et al. ⁽²⁵⁾ muestra el final de la estación lluviosa más larga al noroeste de Honduras que en las otras dos metodologías donde lo muestran con final más temprano (tonos más cafés). El rango del término de la estación lluviosa como resultado de los tres criterios aplicados para el conjunto de datos es del 28 de agosto al 26 de diciembre.

Los diagramas de caja para el conjunto de datos de TELL son mostrados en la figura 5. Para estas series no se presentaron valores atípicos. Los resultados de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾ muestran la mediana el 13 de noviembre, el primer cuartil el 20 de octubre y el tercer cuartil el 6 de diciembre, los resultados de la metodología modificada muestran la mediana el 9 de noviembre, el primer cuartil el 25 de octubre y el tercer cuartil el 30 de noviembre, por otro lado, los resultados de aplicar la metodología de Bombardi et al. ⁽²⁵⁾ muestran la mediana el 4 de diciembre, el primer cuartil el 10 de noviembre y el tercer cuartil el 12 de diciembre.

Figura 4 Similar a la figura 2 para el término de la época lluviosa. 

Figura 5 Diagramas de caja que del término de la estación lluviosa del promedio obtenido para todas las estaciones de la figura 1: (TELL) Resultados de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾ (vf_mod) resultados de aplicar la metodología modificada propuesta, (vf_Bmb) resultados de aplicar la metodología propuesta por Bombardi et al. ⁽²⁵⁾. Se muestran la fecha del primer cuartil, la mediana y tercer cuartil. 

Los promedios de la duración de la época lluviosa se muestran en la figura 6, la cual está definida para cada estación como TELL − IELL + 1. Se obtuvieron resultados en su mayoría una duración prolongada de la estación lluviosa (colores más verdes) desde el centro de Costa Rica hacia el Sur y en la isla de San Andrés en las tres metodologías. Los conjuntos de datos de las diferentes metodologías muestran una duración de la estación lluviosa más reducida en promedio en estaciones en el Pacífico Norte de Costa Rica, Nicaragua, parte sur y central de Honduras, El Salvador, sur de Guatemala, y en el norte de Belice. Contrasta en la metodología utilizada por Bombardi et al. ⁽²⁵⁾ estaciones en el centro y este de Belice con una duración más prolongada de la ELL. En general para estos resultados no se muestra una duración corta de la estación lluviosa (colores más cafés).

Los diagramas de caja de la figura 7 muestran que los resultados de aplicar las tres metodologías no generan valores atípicos. Para los resultados de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾ se muestra la duración de la ELL con una mediana de 186 días, el primer cuartil de 154 días y el tercer cuartil 217 días. Para la rutina modificada propuesta la mediana es de 185 días, el primer cuartil de 169 días y el tercer cuartil de 201 días. Los resultados de la metodología de Bombardi et al. ⁽²⁵⁾ muestran una mediana de 230 días, el primer cuartil de 204 días y el tercer cuartil de 250 días.

Los resultados de las 3 metodologías muestran una duración más prolongada de la época lluviosa del centro de Costa Rica hacia el sur, en estaciones en la costa este de Belice, norte de Guatemala y San Andrés. Los resultados de aplicar el criterio usado por Bombardi et al. ⁽²⁵⁾ muestra en Belice una duración mayor de la época lluviosa, este método de selección del criterio pudo ser sensible a frentes fríos, provocando falsas detecciones del término o inicio de la estación lluviosa Zárate-Hernández ⁽²⁷⁾. Se muestra una zona con menor duración de la época lluviosa en Nicaragua, Honduras, El Salvador y algunas estaciones al sur de Guatemala congruente con el área del Corredor Seco Centroamericano propuesta por Quesada-Hernández et al. ⁽²⁸⁾.

Figura 6 Promedio de la duración de la época lluviosa para las 101 estaciones entre 1968 al 2012 utilizando: (a) resultados de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾, (b) resultados de la metodología modificada propuesta, (c) resultados de la metodología utilizada por Bombardi et al. ⁽²⁵⁾. 

Figura 7 Diagramas de caja de la duración de la estación lluviosa del promedio obtenido para todas las estaciones de la Figura 1. (dur_ELL): resultados de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾, (dur_vmod): metodología modificada propuesta, (vf_Bmb): metodología propuesta por Bombardi et al. ⁽²⁵⁾. Se muestran las fechas del primer cuartil, la mediana y tercer cuartil. 

Las diferencias de medias para el IELL (figura 8) muestra diferencia hasta de 90 días para el inicio de la estación lluviosa. Entre los resultados se puede observar que la metodología modificada detecta un inicio más temprano desde Costa Rica hacia el sur, por otro lado, se observa inicios más tardíos en los resultados de Alfaro e Hidalgo ⁽¹⁶⁾ para estaciones en Nicaragua y alrededores. Las pruebas de diferencia de medias muestran que para la gran mayoría de las estaciones existe una diferencia significativa al 95% para la prueba t y para la prueba Wilcoxon entre los resultados de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾ y la metodología propuesta. Ambas pruebas de diferencia de medias arrojaron resultados de series con promedio estadísticamente igual entre metodología propuesta y la utilizada por Bombardi et al. ⁽²⁵⁾ (figura 8c y 8d).

Figura 8 Diferencia de los promedios del inicio de la época lluviosa entre: (a) y (b) los resultados de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾ y la metodología propuesta (c) y (d) resultados de la metodología de Bombardi et al. ⁽²⁵⁾ y resultados de la metodología propuesta. Se muestra la prueba de hipótesis de diferencia de medias: prueba t en (a) y (c), prueba de suma de rangos Wilcoxon en (b) y (d). Donde h = 0, no se puede rechazar la hipótesis nula (las series provienen de muestras con medias iguales) y h = 1, se rechaza la hipótesis nula (las series no provienen de muestras con medias iguales), ambas pruebas con nivel de significancia al 95%. 

Se obtuvieron diferencias de hasta 60 días para el TELL (figura 9). Diferencias cercanas a cero se muestran entre los resultados de Alfaro e Hidalgo ⁽¹⁶⁾ y la metodología modificada, en Panamá, el Pacífico Norte de Costa Rica, la costa Pacífico de Nicaragua, San Andrés, Honduras y en Belice. Los resultados de la prueba t y Wilcoxon para estas series de tiempo mostraron que la mayoría de las series tienen promedios significativos diferentes.

Figura 9 Igual a la figura 8 para el término de la estación lluviosa. 

Los resultados de la diferencia de media entre las series muestran hasta 140 días entre la duración de la ELL (figura 10), los resultados entre la metodología de Alfaro e Hidalgo ⁽¹⁶⁾ y la metodología modificada muestran diferencias negativas al norte de América Central y diferencias más positivas desde el centro de Costa Rica hacia el sur. Entre la metodología modificada y la usada por Bombardi et al. ⁽²⁵⁾ se puede notar que en Belice se tienen diferencias más positivas. Diferencias cercanas a cero se obtuvieron cerca del Canal de Panamá, Pacífico Norte de Costa Rica y Nicaragua. Para estos mapas se observa una predominancia en diferencias significativas entre los promedios de sus series.

Figura 10 Diferencia del promedio de la duración de la época lluviosa entre: (a) y (b) los resultados de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾ y la metodología propuesta, (c) y (d) resultados de la metodología de Bombardi et al. ⁽²⁵⁾ y resultados de la metodología propuesta. Se muestra la prueba de hipótesis de diferencia de medias: prueba t en (a) y (c), prueba de suma de rangos Wilcoxon en (b) y (d). Donde h = 0, no se puede rechazar la hipótesis nula (las series provienen de muestras con medias iguales) y h = 1, se rechaza la hipótesis nula (las series no provienen de muestras con medias iguales), ambas pruebas con nivel de significancia al 95%. 

Conclusiones

Este estudio comparó los resultados del IELL, TELL y la duración de la ELL de tres metodologías donde se observa que para el IELL entre resultados de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾ y la metodología propuesta una mediana muy similar, siendo los resultados de Alfaro e Hidalgo ⁽¹⁶⁾ más dispersos en la detección de un IELL más temprano y más tardío para el mismo conjunto de estaciones (figuras 2a, 2b y 3). La metodología de Bombardi et al. ⁽²⁵⁾ detecta el IELL con una mediana casi 20 días y con inicios más tempranos que las otras dos metodologías. Esto indica que esta metodología es susceptible a frentes fríos detectando falsos inicios de la ELL (figuras. 4c y 5).

Al igual que en el IELL, los resultados para el TELL entre de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾ y la metodología propuesta poseen una mediana muy similar, los resultados de Alfaro e Hidalgo siendo más dispersos, con TELL más tempranos (figuras 2a, 2b y 3). Los resultados de Bombardi et al. ⁽²⁵⁾ arrojaron un TELL más tardío, con una mediana casi de 25 días después que las otras dos metodologías (figuras 4c y 5).

En cuanto a la duración de la ELL la mediana entre los resultados de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾ y la metodología propuesta es casi la misma solo con un día de diferencia (186 y 185 días respectivamente), sin embargo, los resultados de Alfaro e Hidalgo ⁽¹⁶⁾ son mucho más dispersos que las otras metodologías (figuras 6a, 6b y 7). La metodología de Bombardi et al ⁽²⁵⁾ muestra una duración de la ELL en general más larga que las otras dos metodologías, con una mediana de 230 días.

Por último, las pruebas de hipótesis realizadas muestran en su mayoría que las series del IELL, TELL entre los resultados de la metodología de Alfaro e Hidalgo ⁽¹⁶⁾ y la metodología propuesta son significativamente diferentes y con diferencia hasta de 90 días en algunas estaciones (figuras 8a, 8b, 9a y 9b). Los resultados de la duración de la ELL entre estas metodologías muestran diferencia de hasta 140 días en algunas estaciones. Las pruebas de hipótesis de diferencia de media entre los resultados de la metodología de Bombardi et al. ⁽²⁵⁾ y la metodología propuesta son estadísticamente similares solo para el IELL (figuras 8c y 8d ).

Tomando en cuenta las figuras 2, 4 y 6, el valor de tendencia central y la dispersión de las figuras 3, 5 y 7, se determina que la metodología propuesta genera resultados estadísticamente similares en tendencia central a los determinados por Alfaro e Hidalgo ⁽¹⁶⁾ y con menos variabilidad. Por otro lado, la metodología aplicada por Bombardi et al. ⁽²⁵⁾ detecta falsos inicios de la ELL, extienden el TELL y alarga la duración de la ELL en el sur de América Central, norte de Honduras y en Belice.

Agradecimientos

Al programa de Posgrado en Ciencias de la Atmósfera de la Universidad de Costa Rica (SP5905). La presente investigación se elaboró en el marco de varios proyectos inscritos en la Vicerrectoría de Investigación y Acción Social de la Universidad de Costa Rica: 805-C0- 610 (apoyado por el Fondo de Estímulo), EC-497 (apoyado por el fondo de extensión FEESCONARE), 805-C0-074, 808-C2-103 y 805-B9-454 (apoyado por el Fondo de Grupos).