Next Article in Journal
Effect of Water Deficit on Morphoagronomic and Physiological Traits of Common Bean Genotypes with Contrasting Drought Tolerance
Next Article in Special Issue
Transformation of the Flow Regime of a Large Allochthonous River in Central Europe—An Example of the Vistula River in Poland
Previous Article in Journal
Physics-Based Simulation of Ocean Scenes in Marine Simulator Visual System
Previous Article in Special Issue
Relationship between Water Temperature of Polish Rivers and Large-Scale Atmospheric Circulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Streamflow Intensification Driven by the Atlantic Multidecadal Oscillation (AMO) in the Atrato River Basin, Northwestern Colombia

1
Department of Geography, Faculty of Humanities, Universidad del Valle, Calle 13 #100-00, Cali 25360, Colombia
2
Postgraduation Program CLIAMB, Instituto Nacional de Pesquisas da Amazônia (INPA)—Universidade do Estado do Amazonas (UEA), Ave. André Araújo 2936, Manaus 69060-001, Brazil
3
Instituto Nacional de Pesquisas Espaciais, Centro de Previsão de Tempo e Estudos Climáticos, Divisão de Modelagem e Desenvolvimento. Av. dos Astronautas 1758, São José dos Campos 12227-010, Brazil
4
Escola Superior de Tecnologia, Universidade do Estado do Amazonas, Av. Darcy Vargas 1200, Parque 10 de Novembro, Manaus 69065-020, Brazil
5
Department of Agricultural Engineering, Federal University of Viçosa, Viçosa 36570-900, Brazil
6
Research Group in Water Resources Engineering and Soil (IREHISA), School of natural Resources and Environmental Engineering, Universidad del Valle, Calle 13 #100-00, Cali 25360, Colombia
7
Research Group of Hydrological and Environmental Modelling (GIMHA), Research Institute of Water and Environmental Engineering (IIAMA), Universitat Politècnica de València, 46022 Valencia, Spain
8
Research Group of Perception and Intelligent Systems (PSI), School of Electrical and Electronic Engineering, Universidad del Valle, Calle 13 #100-00, Cali 25360, Colombia
9
Engineering School, Natural and Environmental Resources Engineering School (EIDENAR), Universidad del Valle, Calle 13 #100-00, Cali 25360, Colombia
*
Author to whom correspondence should be addressed.
Water 2020, 12(1), 216; https://doi.org/10.3390/w12010216
Submission received: 4 November 2019 / Revised: 17 December 2019 / Accepted: 23 December 2019 / Published: 13 January 2020
(This article belongs to the Special Issue Effects of Oceanic-Atmospheric Oscillations on Rivers)

Abstract

:
The impact of the Atlantic Multidecadal Oscillation (AMO) on the variations in the streamflow in the Atrato River Basin (ARB) during the 1965–2016 period was analyzed here by considering the cold (1965–1994) and warm (1995–2015) phases of this oscillation. The mean streamflow increased after 1994 (AMO phase change). This increase is related to the strengthening of the zonal gradients of the sea surface temperature (SST) and sea level pressure (SLP) between the tropical central Pacific and the tropical Atlantic after 1994 (warm AMO phase). These gradients contributed to strengthen the Walker cell related upward movement over northern and northwestern South America, in particular during November-December (ND). Consistently, the frequency (R20 mm) and intensity (SDII) of extreme daily rainfall events increased during the 1995–2015 period. Our results show a connection between the AMO and the increase in the streamflow in the ARB during the last five decades. These results contribute to the studies of resilience and climate adaptation in the region.

Graphical Abstract

1. Introduction

Colombia is vulnerable to hydrometeorological risks (e.g., river overflows, floods, and landslides) caused by several phenomena of different space-temporal scales [1,2]. The hydroclimate regimes in Colombia are mainly driven by the latitudinal migration of the Intertropical Convergence Zone. However, other factors may also have a significant contribution. Some of them are the orography of the Andes, thermal driven local circulations [3,4,5,6,7,8,9], earth-atmosphere feedback processes [10], atmospheric rivers [7], the Chocó Jet and the associated mesoscale convective systems [11,12], and ocean-atmospheric processes that occur in the Atlantic and Pacific Oceans [13,14,15,16,17].
The El Niño Southern Oscillation (ENSO) in the interannual time scales is the main ocean-atmospheric coupled phenomenon that affects the rainfall distribution in tropical regions and causes droughts and/or floods in the watersheds involved [15,18,19]. Several studies have documented the ENSO effects on rainfall and streamflow in northern and northwestern South America (SA) [1,18,20,21,22,23,24,25] with consequent negative impacts on the economy, biodiversity and human health in Colombia [6,25,26]. The ENSO warm phase or El Niño (EN) is associated with negative anomalies of the rainfall, soil moisture, evapotranspiration, and streamflow, mainly in central, northern, and western Colombia [1,20,25,27,28], which are caused by variations in both the Walker and Hadley cells [29,30,31]. The La Niña (LN) events are related to the opposite anomaly patterns of these variables. At interdecadal scales, the average discharges of Colombian rivers and precipitation correlate significantly and negatively with the North Atlantic Oscillation (NAO) index and the Pacific Decadal Oscillation (PDO) index and the rainfall in the North Caribbean stations correlate significantly and positively with the sea surface temperature (SST) in the tropical Atlantic [13,15,32].
However, the ARB streamflow might present low-frequency variability which is modulated by the SST low-frequency variability mode. One of these modes in known as the Atlantic Multidecadal Oscillation (AMO) in the North Atlantic [33,34]. This mode is characterized by anomalies of the same signal in the North Atlantic with two centers of maximum anomalies at 15° N and 55° N and a periodicity between 65 and 80 years. The warm phase (positive) occurred during the 1860–1885 and 1930–1964 periods and the cold phase (negative) during the 1895–1924 and 1965–1994 periods [35].
Several studies have provided evidence of the AMO influence on the interannual climate variability in many regions in SA [36,37,38,39,40,41]. Chiessi et al. [40] indicated the significant effect of the AMO on the SA monsoon system, such that during the warm AMO (WAMO) phase the main area with summer rainfall is displaced northward, affecting the rainfall and runoff in the region. Dong et al. [41] showed that the amplitude of the Niño 3 index during WAMO is reduced in relation to the CAMO, so that the AMO can modulate ENSO variability. Additionally, Kayano and Capistrano [36] demonstrated that EN events during the cold AMO (CAMO) phase were generally stronger than those during WAMO phase and that the EN-related rainfall anomalies were more organized and extensive during the CAMO than in the other AMO phase. They also showed that LN events were stronger during the WAMO than during the CAMO. Timmermann et al. [39] found that the strong (weak) ENSO variability during the CAMO (WAMO) phase reduces (increases) the magnitude of the SST annual cycle in the tropical Eastern Pacific, and modulates the ENSO effect on the South American rainfall. In addition, Kayano et al. [38] indicated that the negative multidecadal relation between AMO and the PDO, during the WAMO reinforces the anomalous cooling of the tropical Pacific and favors the settlement of LN events. This leads to a configuration of the dipolar pattern with positive rainfall anomalies in northern SA and negative anomalies in central and eastern SA, especially from winter to summer, which is modulated by the anomalous Walker and Hadley cells. In turn, during the CAMO and warm PDO phase, the warming in the tropical Pacific is reinforced and the LN pattern and its effects on SA rainfall are weakened.
Floods, droughts, and other extreme hydrometeorological events have become recurrent natural disasters over the past several decades, and are expected that these types of hydrometeorological threats will intensify in the future [42,43]. The 2010–2012 LN event caused extreme rainfall and flooding in Colombia with disasters in the Magdalena and Cauca river valleys [44,45] and in the Atrato River valley, a region associated with a high degree of social vulnerability in the country [46]. The ARB is located in extreme northwestern Colombia known as the Chocó Biogeographic (Figure 1) and is one of the richest regions in the world in terms of natural and cultural diversity [47,48,49]. It is a region with some of the highest rainfall levels on earth, with an annual mean rainfall greater than 12,500 millimeters [7,50,51]. Therefore, understanding this region’s extreme hydrometeorological events is a priority for the conservation and planning of water resources of a region that shows strong contrasts.
In this sense, a better understanding of the influence of the low-frequency changes associated with the AMO in the ARB streamflow is aimed in this study, with the main focus on the peaks during May-June (MJ) and November-December (ND). The years with extreme ARB streamflow were stratified according to the AMO phases during the last five decades (1965–2015). For these years, the related anomalous oceanic and atmospheric patterns and the anomalies in the total rainfall and extreme rainfall indices were examined. Section 2 describes the data and methodology. The mean anomaly fields for SSTs, sea level pressures (SLP), vertical velocities at 500 hPa and rainfall are presented and discussed in Section 3. Finally, the conclusions are presented in Section 4.

2. Materials and Methods

2.1. Study Area

The ARB is located between 5° and 9° N and between 75° and 80° W covering an area of approximately 37,730 km2 and an altitude range of 1–3820 meters above the sea level (Figure 1a). The Atrato river is born in the Western Cordillera, running 600 km before delivering nearly 4140 m3·s−1 each year to the Urabá Gulf in the Caribbean Sea [52,53,54,55]; one of the longest navigable rivers in Colombia, only exceeded by the Magdalena and the Cauca Rivers [56]. Moreover, its specific streamflow of 0.109 m3·s−1km−2, is among the world’s highest runoff rates [7,25,53,54].
The precipitation annual cycle in the ARB (Figure 1b) is dominated by a uniform regime with rainfall throughout the year, especially between April and December, that exceeds the 450 mm·month−1 [57,58]. According to Velasquez-Restrepo et al. [53] and the Institute of Hydrology, Meteorology and Environmental Studies (IDEAM) [54], the upper basin (i.e., upstream Bellavista stream gauge station) has a drainage area of 15,358 km2, a mean annual streamflow of 2500 m3·s−1 for the 1965–2015 period, two notable peaks during May–June (MJ, 2640 m3·s−1) and November–December (ND, 2820 m3·s−1), and the lowest streamflow (1770 m3·s−1) in March (Figure 1c). So, considering the streamflow annual cycle the analyses focus on peaks during MJ and ND.

2.2. Hydrometeorological Data

The stream flow data of the 1965–2015 period at the stream gauge stations of Belén, San Antonio, and Bellavista (at 26, 18 and 12 m a.s.l., respectively) from the IDEAM, located in upper ARB were used in preliminary analyses (Appendix A) based on the principal component analysis (PCA) and correlations among the stations. These analyses indicated that the Bellavista streamflow time series reproduced most of the variations shown in the other two streamflow time series. So, the streamflow data at Bellavista was used here as representative of the study domain.
Monthly precipitation data were obtained from the observed meteorological network (from IDEAM). Only the most extensive records at Tagachi, Buchado, and Bellavista were considered (Figure 1a). Hence, the average precipitation of these three stations represent the average precipitation in the ARB. Furthermore, daily precipitation data from the Global Meteorological Forcing Dataset for Land Surface Modeling (GMFD) [59] at a horizontal resolution of 0.25° × 0.25°were also used. The GMFD data result from a combination of the global datasets and reanalysis from NCEP–NCAR reanalysis, Climatic Research Unit monthly climate variables, Global Precipitation Climatology Project (GPCP) daily precipitation, Tropical Rainfall Measuring Mission three-hourly precipitation and the NASA Langley monthly surface radiation budget. This dataset provides a significant improvement over the original reanalysis variables and has been used for a wide variety of applications and diagnostic studies in hydroclimatology [24,60,61]. This data was obtained in the area limited at 65° W, 85° W, 20° N and 5° S for the 1965–2015 period.

2.3. Atmospheric and Oceanic Data

Monthly reconstructed SST dataset (National Oceanic and Atmospheric Administration/Extended Reconstructed Sea Surface Temperature—NOAA/ERSST) [62], version 5, derived from the International Comprehensive Ocean-Atmosphere Dataset [63] was used. The monthly SLP and 500 hPa vertical velocity in pressure coordinate data were obtained from the National Center for Environmental Prediction and National Center for Atmospheric Research (NCEP/NCAR) [64] Reanalysis 1 project [65]. The SST data have a latitude-longitude resolution of 2° × 2° and the SLP and vertical velocity data, of 2.5° × 2.5°. These atmospheric and oceanic data were obtained in the domain limited at 120°E, Greenwich longitude, 40° N and 40° S during the 1965–2015 period. The monthly standardized anomalies of these variables were based on the 1965–2015 period. Additionally, the AMO index used in this study was accessed at (https://www.esrl.noaa.gov/psd/data/timeseries/AMO/) [35]. This index is based on the Kaplan SST data.

2.4. Extreme Rainfall Index

Daily rainfall data from GMFD were used to calculate the extreme rainfall indices proposed by the Expert Team on Climate Change Detection and Indices (ETCCDI) [66,67]. According to Folland et al. [68], changes in these extremes often have greater impacts than the changes in the mean values and can be climate risk indicators. Two indices related to extreme rainfall were calculated at the monthly scale. The frequency index (R20 mm) characterizes heavy rains and is the number of days in a month for which the daily rainfall is greater than or equal to 20 mm [69]. The simple daily intensity index (SDII) is the ratio between the total annual rainfall and the number of wet days (≥1 mm). These indices have been widely used to quantify severe rainfall and to detect possible relationships between hydrological threats and areas that are potentially vulnerable to catastrophic events [24,70,71,72,73,74,75].

2.5. Mean SST Anomaly Patterns During the AMO Phases

Following the methodology proposed by Kayano et al. [76,77], the reconstructed SST data of the 1870–2015 period were used only to support that the SST anomaly patterns during the 1965–1994 and 1995–2015 periods indeed contain, respectively, the CAMO and WAMO features [35,78]. In order to remove the climate change effects, the linear trend for the 1870–2015 of the time series at each grid point was eliminated using the method of least squares. Monthly standardized anomalies based on the 1965–2015 period were obtained in the domain limited at 60° N, 60° S, 120° E and the Greenwich longitude. Then, the bimonthly SST anomaly mean states were obtained for MJ and ND during the 1965–1994 and 1995–2015 periods. The statistical significance was assessed using the Student’s t-test with a confidence level of 95% and the degrees of freedom are the number of years in each period.

2.6. Composite Analysis

The extreme streamflow events in MJ and ND of the 1965–2015 period were selected and stratified according to AMO phases. Bimonthly streamflow time series for MJ and ND at Bellavista stream gauge station were obtained. These bimonthly series were then divided into two time series, one covering 30 years of the 1965–1994 period and the other, 21 years of the 1995–2015 period. These periods are referred to as the periods before and after 1994. Subsequently, the percentile time-series (R) were constructed by ranging the years from one to 30 (one to 21) and dividing them by 30 (21). Then, for each period, the extreme events above the percentile 75 were selected.
First, monthly standardized anomalies of precipitation, SST, SLP, and 500 hPa vertical velocity in pressure coordinate during the 1965–2015 period were obtained and used in the composite analyses. In this case, the means and standard deviations were based on the 1965–2015 period. Bimonthly composites of the anomalies of the rainfall, SST, SLP, and 500 hPa vertical velocity in pressure coordinate were obtained for MJ and ND. Furthermore, the composites for the extreme rainfall indices R20 mm and SDII were also obtained for MJ and ND. Composite analyses of the extreme events before and after 1994 were performed separately. The difference composites between these two periods are also presented.
The statistical significance of the composites was evaluated using the Student’s t-test at 95% confidence level, and the degrees of freedom, the number of events for each period. The significance test used for comparing two means [79] considers that a variable X with n values, standard deviation S and Student’s t distribution has a significant mean when the absolute value of the mean exceeds t α ( n 1 ) S / ( n 1 ) ; in which t α ( n 1 ) corresponds to the cumulative probability for (n − 1) degrees of freedom in Student’s table and α = 0.05 . The two-sample t-test of Levine and Wilks was used to evaluate the statistical significance of differences in the composites before and after 1994 [80]. This parametric test considers two variables, X 1 and X 2 , with values n 1 and n 2 , standard deviations S 1 and S 2 and composites indicated by X 1 ¯ and X 2 ¯ and assumes that the difference X 1 X 2 has a Student t distribution with n 1 + n 2 2 degrees of freedom and that the absolute values of X 1 ¯ X 2 ¯ greater than t α ( n 1 + n 2 2 ) ( n 1 1 ) S 1 2 + ( n 2 1 ) S 2 2   n 1 + n 2 n 1 n 2 ( n 1 + n 2 2 ) are statistically significant.

3. Results and Discussion

3.1. SST Mean States

Figure 2 shows the bimonthly SST anomaly patterns during the two selected periods. During 1965–1994, both seasons show significant negative anomalies in most of the North Atlantic and in the tropical West Pacific, and the negative ones in the extratropical South Atlantic and Southeast Pacific, with the negative anomalies in Tropical North Atlantic (TNA) being significant during ND and non-significant during MJ (Figure 2a,c). The maps during the 1995–2015 period show inverted signal patterns (Figure 2b,d). These maps show previously documented typical AMO features [35,39,76,81]. Thus, the two selected periods present contrasting SST mean states, under which the streamflow and rainfall in the ARB may have distinct characteristics which may be related to the low-frequency background associated with the AMO.

3.2. Difference of the Average Precipitation in the ARB and Streamflow in Bellavista

The differences between the two AMO phases of the bimonthly streamflow at Bellavista and precipitation in the ARB in MJ and ND were evaluated using the Student’s t-test (Table 1)). The t-test has been used to compare the differences between climate and environmental time series [82,83]. The normality of the series was evaluated using the Shapiro and Wilk [84] and Anderson and Darling [85] tests. All time series satisfied the normality test. The differences in the means were analyzed at a significance level of 5%. Table 1 shows the t-test comparisons of streamflow and precipitation time series before and after 1994 in MJ and ND. In both seasons, the mean streamflow at Bellavista shows a significant difference before and after 1994 and a statistically significant increase (p < 0.05) in the mean discharge from the period before to after 1994. Also, there was a statistically significant increase in the average precipitation in the ARB from the period before to after 1994 in ND.
Figure 3 shows the ARB streamflow and AMO index, both filtered with a 10-year moving mean filter. In both seasons, the streamflow presents a significant (p < 0.05) increase in the mean values of the period before and after 1994, which accompanies the change of the AMO phase. The mean streamflow values before 1994 and after were of 2530 m3·s−1 and 2800 m3·s−1 in MJ, and 2615 m3·s−1 and 3106 m3·s−1 in ND (Figure 3). These results are consistent with the significant correlations of 0.68 (p < 0.05) between the filtered MJ streamflow and AMO index time series, and 0.75 between the filtered ND streamflow and AMO index time series. Similar graphics were constructed for the average precipitation in the ARB (Figure 4). The ND precipitation time series shows a significant (p < 0.05) increase with the means values before and after 1994 of 472 mm and 579 mm. The filtered streamflow and precipitation time series show higher correlation in ND (r = 0.50) than in MJ (r = 0.43). Therefore, there were estimated increases from the period before to after 1994 of 12% and 19% in the streamflow in MJ and ND, and of 21% in the precipitation in ND.

3.3. Composite Analysis

Composites considering the streamflow extreme events were calculated in the periods before and after 1994. The extreme events selected for the CAMO and WAMO phases are listed in Table 2.
Figure 5 and Figure 6 show the SST, SLP, and 500 hPa vertical velocity anomaly patterns during the periods before and after 1994 and their differences. The SST composite during the CAMO shows significant negative anomalies in the subtropical North Atlantic, the Tropical South Atlantic (TSA), the equatorial Western and subtropical South Pacific during MJ and in TNA, and most of the tropical Pacific during ND (Figure 5a,g). Additionally, small areas of anomalous warming occurred in the tropical central Pacific in MJ and in the subtropical areas of the Pacific and TSA in ND. During the WAMO, the SST anomaly patterns were characterized by an anomalous warming in a boomerang-shaped area in the Australasian region in both bimonthly periods and in most of the subtropical North Atlantic in MJ and in the TNA in ND, while negative anomalies were depicted in the central tropical Pacific in MJ (not significant) and ND (significant) (Figure 5b,h). This analysis therefore shows that the greatest differences in SST composites between the two AMO phases occurred in the TNA, subtropical North Atlantic, and central tropical Pacific, which defined the SST variations associated with the low-frequency background shown in Section 3.1 (Figure 5c,i). Furthermore, the SST anomalies exhibited higher magnitudes and greater spatial extents in the tropical sector in ND than in MJ.
Since 1994, anomalous cooling in the tropical Eastern Pacific and warming in the central North Pacific in ND indicate the establishment of cold PDO phase (CPDO). This result agrees with Poveda et al.’s [1] findings on the negative relation between the PDO index and the streamflow in Colombia. During the same two-month period, positive SST anomalies in the TNA and equatorial Atlantic characterize the WAMO. The positive SST anomalies in these Atlantic sectors, together with the negative SST anomalies in the tropical Pacific characterize an Atlantic-Pacific inter-basin positive SST gradient mode (Figure 5b,h), which creates a positive feedback modulating the climate variability in these Oceanic basins. Previous studies argued that the WAMO/CPDO low-frequency mean state strengthens the cooling in the tropical Pacific and, therefore, leads to a higher percentage of La Niña (LN) events [38,86]. This process explains why the LN-related SST anomaly patterns in the tropical Pacific are stronger after 1994 (Figure 5h). These results reinforce the fundamental role of the Atlantic in determining the intensities of the LN events and, consequently, the rainfall in SA [36,38,87,88,89]. Coherently, all events selected in ND after 1994 coincided with the LN event occurrences (Table 2). Choi et al. [90] showed that the ENSO decadal oscillation presents a stronger relationship with the SST mean state during its mature stage (November–January). Li et al. [91] and Barichivich et al. [92] observed, since the beginning of the 1990s, a strong warming trend in the tropical Atlantic and Indo-Western Pacific Oceans and a cooling trend in the Eastern Pacific, which intensifies the LN response in the tropical Pacific. The results in the present study confirm their findings (Figure 5c,i). Consistent with the SST composites, the SLP anomaly patterns after 1994 show an abnormally high pressure center in the central and eastern tropical Pacific, and an anomalous low pressure center in the TNA, both centers more intense in ND than in MJ (Figure 5e,k). These patterns lead to a stronger inter-basin SLP gradient, especially during ND (Figure 5k).
Consistent with the SST and SLP anomalies, the anomalous Walker cells before 1994 featured downward motion in the tropical Eastern Pacific and northern and northwestern SA and upward motion in the central tropical Pacific in MJ, and downward motion over most of the tropical Pacific between the South American coast and 180° W and upward motion on the northern coast of SA and northeastern Brazil in ND (Figure 6a,d). On the other hand, after 1994, the SST and SLP gradients between the tropical central Pacific and the tropical Atlantic strengthened the Walker cell with an upward motion to northern and northwestern SA and downward in the tropical central Pacific and the Caribbean Sea in both bimonthly periods (Figure 6b,e). The composite differences represent the intensification (weakening) of the SLP in the Pacific (tropical Atlantic and northwestern SA) associated with anomalous cooling (warming) in the tropical Pacific (TNA and subtropical North Atlantic) and the intensification of Walker circulation with downward (upward) motion in the tropical central Pacific (northern and northwestern SA) after 1994 (Figure 6c,f). Since the Walker circulation connects the tropical oceans, and its upward branch extends to northern and northwestern SA, intensification of the atmospheric circulation and the convergence of humidity on these regions are expected, increasing the rainfall and the streamflow in the ARB, with the highest intensity in ND, when the pattern is better established. The results show a physical connection between the ARB streamflow and the recent SST and SLP inter-basin gradients, boosted by the Atlantic warming during the WAMO, as observed in other studies [91,93,94,95].
Our results are in agreement with the finding that the atmospheric response to the tropical Atlantic warming can lead to an intensification of the upward branch of the Walker circulation in northern SA [96,97,98]. This process is accompanied by a downward motion on the central equatorial Pacific and upward in the Indo-Pacific region. These anomalous motions increase the zonal pressure gradient between the tropical Pacific and Atlantic and increase the divergence at the surface, generating a cooling in the central equatorial Pacific, thus favoring the development of LN events. This is in agreement with the above analyzed SST anomaly patterns with large negative SST anomalies in the tropical eastern Pacific after 1994 in ND (Figure 5h).
Consistent with the anomalous SST and SLP patterns with cooling in the tropical Pacific, especially during WAMO, positive rainfall anomalies were observed in Colombia in both bimonthly periods (left column in Figure 7 and Figure 8). After 1994, the positive rainfall anomalies increased mainly in the Caribbean region and in the Andes region adjacent to the ARB in MJ (left column in Figure 7). After 1994, significant positive anomalies with larger magnitudes in ND than in MJ covered a large part of the Colombian territory, including the ARB (left column in Figure 8). The rainfall increases are associated with the strengthened zonal SST and SLP gradients between the tropical central Pacific and Atlantic in ND (Figure 5h,k). These gradients during ND strengthen the anomalous upward motion in northern and northwestern SA, increasing the rainfall in the ARB and the Colombian Caribbean. An examination of the anomalous rainfall characteristics using the monthly database of the GPCC [99,100] is consistent with the GMFD results, with the greatest rainfall increase in the Colombian Caribbean during MJ (see Appendix B, left column in Figure A2) and more extensive anomalies in ND, covering the northern, central, and western parts of the country (right column in Figure A2).
Severe local rainfall events and environmental changes can often increase the incidence of hydrological risks, such as floods and landslides, especially in tropical regions, where an increase in the intensity and frequency of rainfall extremes in the coming decades is expected [24,75,101,102,103,104]. Here, rainfall extremes were analyzed using the R20 mm index and SDII.
The differences in composite in MJ for the R20 mm show no significant increases in the number of very rainy days in the ARB (middle column in Figure 7), or increases associated with the SDII (Figure 7). This result agrees with the lack of significant differences found in the t-test of the means before and after 1994 for the mean precipitation in the ARB shown in Section 3.2. However, in ND, both extreme rainfall indices describe significant increases in most of Colombia, especially in the Caribbean, Pacific, and Andean regions (middle and right columns in Figure 8). These results show that the upward change of the mean streamflow in the ARB at Bellavista in ND coincides with the increases in intensity and frequency of extreme rainfall events during the WAMO.

4. Conclusions

We evaluated how the low-frequency background associated with the AMO modifies the anomaly patterns of some atmospheric and oceanic variables associated with extreme streamflow events in the Atrato River basin (ARB), particularly during the two bimonthly streamflow peaks in May–June (MJ) and November–December (ND). Through a Student’ t-test, the mean rainfall and mean streamflow were compared considering the cold AMO (CAMO) and warm AMO (WAMO), which refer to the 1965–1994 and 1995–2015 periods, respectively. Our results showed a significant increase in the streamflow means at Bellavista in MJ and ND and in the rainfall means in ND from the first to the second period (Table 1, Figure 4).
Composite analyses of the SST and SLP anomalies based on the streamflow extreme events in MJ and ND during each AMO phase reproduce the main anomaly patterns of these variables previously documented to be related to the AMO [38,39,41,86]. The rainfall anomaly composites show a more defined pattern, with significant positive anomalies in extensive areas during the WAMO phase, which are indeed modulated by the LN events in the tropical Pacific, especially during ND. In other words, rainfall anomalies related to the cooling of the tropical Pacific during WAMO (1995–2015) are more pronounced and occupy more extensive areas than those during the CAMO phase (1965–1994).
Our findings show that the Atlantic Ocean plays an important role in modulating the rainfall and streamflow variability in the ARB. At the low frequency time scale, the SST and SLP anomalies in this Oceanic sector create the east-west SST and SLP gradients between the tropical Atlantic and eastern Pacific Oceans, which during the WAMO phase favors the La Niña event establishment [77]. Physical explanation is that the changes in the tropical Atlantic SST modify the inter-basin (tropical Atlantic and tropical eastern Pacific) pressure gradient, creating an atmospheric bridge that transmits the SST information between the basins [94,105]. The warm tropical Atlantic during the WAMO intensifies the Walker circulation [106] and confirms the previously discussed AMO forcing during the last two decades (1995–2015) [93,107,108,109]. McGregor et al. [94] indicated that the tropical Pacific Ocean shows changes in the average state at a multidecadal scale resulting from the AMO, particularly in the tropics, along with interannual and decadal scales [94,105,109].
The strengthening of the Walker circulation after 1994 played a fundamental role by causing upward motion in northern and northwestern South America (SA) (Figure 6), which contributed to the positive rainfall anomalies in the ARB. Significant increases in intensities (SDII) and frequencies (R20 mm) of extreme rainfall events in ND were found in the ARB during the WAMO (Figure 8). However, in MJ there was no significant increase in the extreme rainfall in the basin (middle and right columns in Figure 7).
Our analyses suggest that the low-frequency background associated with the AMO phases modulated the rainfall and streamflow in the ARB during the last five decades. However, the changes in the streamflow and rainfall in the ARB might also be partially related to the global warming in the 1965–2015 period. The results here improved our understanding on the hydrometeorological variability in the ARB. We believe that our results might contribute to the assessment of climate resilience and the monitoring of hydrometeorological threats. Changes in land use may alter the basin responses to rainfall, however, this issue is outside the scope of this study, and future studies are needed to evaluate possible changes.

Author Contributions

Conceptualization, W.L.C., R.V.A., M.T.K., T.C., and F.F.; methodology, W.L.C., A.A., M.T.K., R.V.A., F.F. and I.A.R.; software was managed by W.L.C., A.A., T.C., and R.A.F.d.S.; validation, W.L.C., M.T.K., and R.V.A.; formal analysis, W.L.C., A.A., M.T.K., and R.V.A.; investigation, W.L.C., A.A., M.T.K. and R.V.A.; data curation, W.L.C., T.C., F.F., and A.A.; original draft preparation, W.L.C., R.V.A., I.A.R., and A.A.; reviewing and editing, W.L.C., A.A., M.T.K., R.V.A., T.C., Y.C.-E., W.A.-M., and I.A.R.; visualization, W.L.C., A.A., T.C., W.A.-M., and I.A.R.; supervision, R.V.A., M.T.K., F.F., R.A.F.d.S., Y.C.-E., and W.A.-M. All authors have read and agreed to the published version of the manuscript.

Funding

The first author was supported by the Doctoral Scholarship of the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brazil (CAPES), Finance Code 001 and Universidad del Valle (Cali-Colombia). The Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) of Brazil partially supported the second and third authors under grants 302322/2017-5 and 305942/2015-8, respectively. The fourth author was supported by the CAPES-Brazil and Universidade Federal de Viçosa for doctoral studies. The fifth author was supported by the Program for Strengthening Regional Capacities in Research, Technological Development and Innovation in the department of Nariño and the CEIBA Foundation for doctoral studies. The seventh author was supported by the CAPES-Brazil. This research was done partially during an internship of the first and fifth author at the Research Institute of Water and Environmental Engineering of the Universitat Politècnica de València. This work was also supported by the Spanish Ministry of Science and Innovation through Project TETISCHANGE (RTI2018-093717-B-I00).

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design, analysis and interpretation of data, in the writing manuscript, or in the decision to publish the results.

Appendix A

The streamflow data in Bellavista, San Antonio and Belén for the 1965–2015 period were provided by the Institute of Hydrology, Meteorology and Environmental Studies (IDEAM) (Figure 1a). These data were subjected to the principal component analysis (PCA) and the correlations among them. For PCA, the methodology proposed by Jolliffe and Cadima [110] was used. For the monthly analysis, the first mode explains 92% of the total monthly streamflow variance. For the bimonthly analyses, the MJ first mode and ND first mode explain 79% and 89% of the corresponding total bimonthly streamflow variances (Figure A1). These results indicate that the three stations are highly correlated with each other and that the first mode of each analysis incorporates the largest amount of information from the original variables. In addition, the monthly correlations between Bellavista and the other stations were greater than 0.76 (p < 0.05) at monthly scale (Table A1). The analyses here ensure that the streamflow time series at Bellavista can be used as representative of the other two stations in the ARB.
Figure A1. First mode of variability for: (a) monthly scale, (b) May–June (MJ), and (c) November–December (ND).
Figure A1. First mode of variability for: (a) monthly scale, (b) May–June (MJ), and (c) November–December (ND).
Water 12 00216 g0a1
Table A1. Stream gauge station significant correlations at the 0.05 level.
Table A1. Stream gauge station significant correlations at the 0.05 level.
Stream gauge stationsMonthlyMJND
Bellavista vs Belén0.720.380.72
Bellavista vs San Antonio0.890.790.83
San Antonio vs Belén0.770.450.80

Appendix B

Figure A2. Patterns of standardized precipitation anomalies associated with extreme flow events at Bellavista in May–June (MJ in left column), and November–December (ND in right column), (CAMO (upper line), WAMO (middle line) and difference (WAMO minus CAMO) (lower line)). The contour range is 0.2 standard deviation. The continuous line (dashed) refers to significant positive (negative) values at a 95% confidence level. The data were acquired from the Global Precipitation Climatology Centre (GPCC 2019) Full Data Reanalysis V2018, at a horizontal spatial resolution of 0.5° × 0.5° [100]. Standardized monthly anomalies were calculated using historical data for the 1965–2015 period.
Figure A2. Patterns of standardized precipitation anomalies associated with extreme flow events at Bellavista in May–June (MJ in left column), and November–December (ND in right column), (CAMO (upper line), WAMO (middle line) and difference (WAMO minus CAMO) (lower line)). The contour range is 0.2 standard deviation. The continuous line (dashed) refers to significant positive (negative) values at a 95% confidence level. The data were acquired from the Global Precipitation Climatology Centre (GPCC 2019) Full Data Reanalysis V2018, at a horizontal spatial resolution of 0.5° × 0.5° [100]. Standardized monthly anomalies were calculated using historical data for the 1965–2015 period.
Water 12 00216 g0a2

References

  1. Poveda, G.; Vélez, J.I.; Mesa, O.; Hoyos, C.; Mejía, J.; Barco, O.J.; Correa, P.L. Influencia de fenómenos macroclimáticos sobre el ciclo anual de la hidrología Colombiana: Cuantificación lineal, no lineal y percentiles probabilísticos. Meteorol. Colomb. 2002, 6, 121–130. [Google Scholar]
  2. Sura, P. A general perspective of extreme events in weather and climate. Atmos. Res. 2011, 101, 1–21. [Google Scholar] [CrossRef]
  3. Mejía, J.; Mesa, O.J.; Poveda, G.; Vélez, J.I.; Hoyos, C.; Ricardo, M.; Barco, J.; Cuartas, A.; Montoya, M.; Botero, B. Distribución espacial y ciclos anual y semianual de la precipitación en Colombia. Dyna Medellin Colomb. 1999, 127, 7–26. [Google Scholar]
  4. Hastenrath, S. The Intertropical Convergence Zone of the Eastern Pacific revisited. Int. J. Climatol. 2002, 22, 347–356. [Google Scholar] [CrossRef]
  5. Poveda, G.; Waylen, P.R.; Pulwarty, R.S. Annual and inter-annual variability of the present climate in northern South America and southern Mesoamerica. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2006, 234, 3–27. [Google Scholar] [CrossRef]
  6. Hoyos, I.; Baquero-Bernal, A.; Jacob, D.; Rodríguez, B.A. Variability of extreme events in the Colombian Pacific and Caribbean catchment basins. Clim. Dyn. 2013, 40, 1985–2003. [Google Scholar] [CrossRef]
  7. Poveda, G.; Jaramillo, L.; Vallejo, L.F. Seasonal precipitation patterns along pathways of South American low-level jets and aerial rivers. Water Resour. Res. 2014, 50, 98–118. [Google Scholar] [CrossRef]
  8. Trojer, H. Meteorología y climatología de la vertiente del Pacífico colombiano. Rev. Acad. Colomb. Ciencias Exactas Físicas Nat. 2018, 41, 467. [Google Scholar] [CrossRef] [Green Version]
  9. Urrea, V.; Ochoa, A.; Mesa, O. Seasonality of Rainfall in Colombia. Water Resour. Res. 2019, 55, 4149–4162. [Google Scholar] [CrossRef]
  10. Poveda, G.; Mesa, O.J. Las fases extremas del fenómeno ENSO (El Niño y La Niña) y su influencia sobre hidrología de Colombia. Ing. Hidráulica en México 1996, XI, 21–37. [Google Scholar]
  11. Jaramillo, L.; Poveda, G.; Mejía, J.F. Mesoscale convective systems and other precipitation features over the tropical Americas and surrounding seas as seen by TRMM. Int. J. Climatol. 2017, 37, 380–397. [Google Scholar] [CrossRef]
  12. Yepes, J.; Poveda, G.; Mejía, J.F.; Moreno, L.; Rueda, C. Choco-jex: A research experiment focused on the Chocó low-level jet over the far eastern Pacific and western Colombia. Bull. Am. Meteorol. Soc. 2019, 100, 779–796. [Google Scholar] [CrossRef]
  13. Poveda, G. La hidroclimatología de Colombia: Una síntesis desde la escala inter-decadal hasta la escala diurna. Rev. Académica Colomb. Ciencias Tierra 2004, 28, 201–222. [Google Scholar]
  14. Zuluaga, I.C.M.; Poveda, G.; Mejia, J.F. Ciclo Diurno de la Precipitación sobre Colombia y el Pacífico Oriental durante 1998–2002 según la misión TRMM. In Proceedings of the XVI Seminario Nacional de Hidráulica e Hidrología, Armenia, Quindio, 29–31 October 2004. [Google Scholar]
  15. Poveda, G.; Álvarez, D.M.; Rueda, Ó.A. Hydro-climatic variability over the Andes of Colombia associated with ENSO: A review of climatic processes and their impact on one of the Earth’s most important biodiversity hotspots. Clim. Dyn. 2011, 36, 2233–2249. [Google Scholar] [CrossRef]
  16. Navarro-Monterroza, E.; Arias, P.A.; Vieira, S.C. El Niño-Oscilación del Sur, fase Modoki, y sus efectos en la variabilidad espacio-temporal de la precipitación en Colombia. Rev. Acad. Colomb. Ciencias Exactas Físicas Nat. 2019, 43, 120. [Google Scholar] [CrossRef]
  17. Poveda, G.; Álvarez, D.M. El colapso de la hipótesis de estacionariedad por cambio y variabilidad climática: Implicaciones para el diseño hidrológico en ingeniería. Rev. Ing. 2012, 36, 65–76. [Google Scholar]
  18. Córdoba-Machado, S.; Palomino-Lemus, R.; Gámiz-Fortis, S.R.; Castro-Díez, Y.; Esteban-Parra, M.J. Influence of tropical Pacific SST on seasonal precipitation in Colombia: Prediction using El Niño and El Niño Modoki. Clim. Dyn. 2015, 44, 1293–1310. [Google Scholar] [CrossRef]
  19. Philander, S.G. El Nino, La Nina, and the Southern Oscillation, 1st ed.; Philander, S., Ed.; Academic Press: New York, NY, USA, 1989; ISBN 0080570984. [Google Scholar]
  20. Pabón, J.D.; Montealegre, J.E. La Variabilidad Climatica Interanual Asociada Al Ciclo El Niño-La Niña–Oscilacion Del Sur Y Su Efecto En El Patron Pluviometrico De Colombia. Meteorol. Colomb. 2000, 2, 7–21. [Google Scholar]
  21. Kayano, M.T.; Andreoli, R.V.; Ferreira de Souza, R.A. Relations between ENSO and the South Atlantic SST modes and their effects on the South American rainfall. Int. J. Climatol. 2013, 33, 2008–2023. [Google Scholar] [CrossRef]
  22. Tedeschi, R.G.; Cavalcanti, I.F.A.; Grimm, A.M. Influences of two types of ENSO on South American precipitation. Int. J. Climatol. 2013, 33, 1382–1400. [Google Scholar] [CrossRef]
  23. Córdoba-Machado, S.; Palomino-Lemus, R.; Gámiz-Fortis, S.R.; Castro-Díez, Y.; Esteban-Parra, M.J. Assessing the impact of El Niño Modoki on seasonal precipitation in Colombia. Glob. Planet. Chang. 2015, 124, 41–61. [Google Scholar] [CrossRef]
  24. Ávila, Á.; Guerrero, F.C.; Escobar, Y.C.; Justino, F. Recent precipitation trends and floods in the Colombian Andes. Water 2019, 11, 379. [Google Scholar] [CrossRef] [Green Version]
  25. Poveda, G.; Jaramillo, A.; Gil, M.M.; Quiceno, N.; Mantilla, R.I. Seasonally in ENSO-related precipitation, river discharges, soil moisture, and vegetation index in Colombia. Water Resour. Res. 2001, 37, 2169–2178. [Google Scholar] [CrossRef]
  26. Córdoba-Machado, S.; Palomino-Lemus, R.; Gámiz-Fortis, S.R.; Castro-Díez, Y.; Esteban-Parra, M.J. Seasonal streamflow prediction in Colombia using atmospheric and oceanic patterns. J. Hydrol. 2016, 538, 1–12. [Google Scholar] [CrossRef]
  27. Puertas, O.L.O.; Carvajal-Escobar, Y. Incidencia de El Niño-Oscilación del Sur en la precipitación y la temperatura del aire en colombia, utilizando el Climate Explorer. Rev. Científica Ing. Desarro. 2011, 104–118. [Google Scholar]
  28. Ávila, Á.; Carvajal, Y.; Gutiérrez, S. El Niño and La Niña analysis influence in the monthly water supply at Cali River basin. Tecnura 2014, 18, 120–133. (In Spanish) [Google Scholar]
  29. Kousky, V.E.; Kagano, M.T.; Cavalcanti, I.F.A. A review of the Southern Oscillation: Oceanic-atmospheric circulation changes and related rainfall anomalies. Tellus A 1984, 36 A, 490–504. [Google Scholar] [CrossRef]
  30. Rasmusson, E.M.; Mo, K. Linkages between 200-mb tropical and extratropical circulation anomalies during the 1986–1989 ENSO cycle. J. Clim. 1993, 6, 595–616. [Google Scholar] [CrossRef] [Green Version]
  31. Zhou, J.; Lau, K.M. Principal modes of interannual and decadal variability of summer rainfall over South America. Int. J. Climatol. 2001, 21, 1623–1644. [Google Scholar] [CrossRef]
  32. Poveda, G.; Gil, M.M.; Quiceno, N. El ciclo anual de la hidrología de Colombia en relación con El ENSO y la NAO. Bull. Am. Meteorol. Soc. 1998, 27, 721–731. [Google Scholar]
  33. Sung, M.K.; An, S.-I.; Kim, B.M.; Kug, J.S. Asymmetric impact of Atlantic Multidecadal Oscillation on El Niño and la Niña characteristics. Geophys. Res. Lett. 2015, 42, 4998–5004. [Google Scholar] [CrossRef] [Green Version]
  34. Schlesinger, M.E.; Ramankutty, N. Low-frequency oscillation. Nature 1994, 367, 723–726. [Google Scholar] [CrossRef]
  35. Enfield, D.B.; Mestas-Nuñez, A.M.; Trimble, P.J. The Atlantic multidecadal oscillation and its relation to rainfall and river flows in the continental U.S. Geophys. Res. Lett. 2001, 28, 2077–2080. [Google Scholar] [CrossRef] [Green Version]
  36. Kayano, M.T.; Capistrano, V.B. How the Atlantic multidecadal oscillation (AMO) modifies the ENSO influence on the South American rainfall. Int. J. Climatol. 2014, 34, 162–178. [Google Scholar] [CrossRef]
  37. Martínez, J.I.; Obrochta, S.; Yokoyama, Y.; Battarbee, R.W. Atlantic Multidecadal Oscillation (AMO) forcing on the late Holocene Cauca paleolake dynamics, northern Andes of Colombia. Clim. Past Discuss. 2015, 11, 2649–2664. [Google Scholar] [CrossRef] [Green Version]
  38. Kayano, M.T.; Andreoli, R.V.; Souza, R.A.F. de El Niño–Southern Oscillation related teleconnections over South America under distinct Atlantic Multidecadal Oscillation and Pacific Interdecadal Oscillation backgrounds: La Niña. Int. J. Climatol. 2019, 39, 1359–1372. [Google Scholar] [CrossRef]
  39. Timmermann, A.; Okumura, Y.; An, S.I.; Clement, A.; Dong, B.; Guilyardi, E.; Hu, A.; Jungclaus, J.H.; Renold, M.; Stocker, T.F.; et al. The influence of a weakening of the Atlantic meridional overturning circulation on ENSO. J. Clim. 2007, 20, 4899–4919. [Google Scholar] [CrossRef] [Green Version]
  40. Chiessi, C.M.; Mulitza, S.; Pätzold, J.; Wefer, G.; Marengo, J.A. Possible impact of the Atlantic Multidecadal Oscillation on the South American summer monsoon. Geophys. Res. Lett. 2009, 36, 1–5. [Google Scholar] [CrossRef] [Green Version]
  41. Dong, B.; Sutton, R.T.; Scaife, A.A. Multidecadal modulation of El Niño-Southern Oscillation (ENSO) variance by Atlantic Ocean sea surface temperatures. Geophys. Res. Lett. 2006, 33, 2–5. [Google Scholar] [CrossRef]
  42. Dottori, F.; Szewczyk, W.; Ciscar, J.-C.; Zhao, F.; Alfieri, L.; Hirabayashi, Y.; Bianchi, A.; Mongelli, I.; Frieler, K.; Betts, R.A.; et al. Increased human and economic losses from river flooding with anthropogenic warming. Nat. Clim. Chang. 2018, 8, 781. [Google Scholar] [CrossRef]
  43. Mora, C.; Spirandelli, D.; Franklin, E.; Lynham, J.; Kantar, M.; Miles, W.; Smith, C.; Freel, K.; Moy, J.; Louis, L.; et al. Broad threat to humanity from cumulative climate hazards intensified by greenhouse gas emissions. Nat. Clim. Chang. 2018, 8, 1062–1071. [Google Scholar] [CrossRef]
  44. Arias, P.A.; Martínez, J.A.; Vieira, S.C. Moisture sources to the 2010–2012 anomalous wet season in northern South America. Clim. Dyn. 2015, 45, 2861–2884. [Google Scholar] [CrossRef]
  45. Enciso, A.; Carvajal Escobar, Y.; Sandoval, M. Hydrological analysis of historical floods in the upper valley of Cauca river. Ing. Compet. 2016, 18, 46. [Google Scholar]
  46. Hoyos, N.; Escobar, J.; Restrepo, J.C.; Arango, A.M.; Ortiz, J.C. Impact of the 2010-2011 La Niña phenomenon in Colombia, South America: The human toll of an extreme weather event. Appl. Geogr. 2013, 39, 16–25. [Google Scholar] [CrossRef]
  47. Marchese, C. Biodiversity hotspots: A shortcut for a more complicated concept. Glob. Ecol. Conserv. 2015, 3, 297–309. [Google Scholar] [CrossRef] [Green Version]
  48. Mittermeier, R.A.; Gil, P.R.; Hoffmann, M.; Pilgrim, J. Hotspots:Earth’s biologically richest and most threatened terrestrial ecoregions. Conserv. Int. 2004, 400. Available online: https://www.academia.edu/1438756/Hotspots_revisited_Earths_biologically_richest_and_most_endangered_terrestrial_ecoregions (accessed on 26 December 2019).
  49. Camargo, X.S. The Ecocentric Turn of Environmental Justice in Colombia. King’s Law J. 2019, 1757–8442. [Google Scholar] [CrossRef]
  50. Hurtado-Montoya, A.F.; Mesa-Sánchez, Ó.J. Reanalysis of monthly precipitation fields in Colombian territory. Dyna 2014, 81, 251. [Google Scholar] [CrossRef]
  51. Durán-Quesada, A.M.; Reboita, M.; Gimeno, L. Precipitation in tropical America and the associated sources of moisture: A short review. Hydrol. Sci. J. 2012, 57, 612–624. [Google Scholar] [CrossRef]
  52. Lasso, C.A.; Gutiérez, F.P.; Morales-Betancourt, M.A.; Agudelo, E.C.; Ramírez-Gil, H.; Ajiaco-Martínez, R.E. Pesquerias Continentales de Colombia: cuencas del Magdalena-Cauca, Sinú, Canalete, Atrato, Orinoco, Amazonas y vertiente del Pacífico; Instituto de Investigación de los Recursos Biológicos Alexander von Humboldt: Bogotá, Colombia, 2011; ISBN 9789588343624. Available online: http://repository.humboldt.org.co/handle/20.500.11761/9332 (accessed on 26 December 2019).
  53. Velásquez-Restrepo, M.; Poveda, G. Estimación del balance hídrico de la región Pacífica. Dyna 2019, 86, 297–306. [Google Scholar] [CrossRef]
  54. IDEAM. 2014-Estudio Nacional del Agua-ENA 2014; Instituto de Hidrología, Meteorología y Estudios Ambientales—IDEAM: Bogotá, Colombia, 2014; ISBN 9789588067704. Available online: http://documentacion.ideam.gov.co/openbiblio/bvirtual/023080/ENA_2014.pdf (accessed on 25 September 2019).
  55. Mosquera-Machado, S.; Ahmad, S. Flood hazard assessment of Atrato River in Colombia. Water Resour. Manag. 2007, 21, 591–609. [Google Scholar] [CrossRef]
  56. Palomino-Ángel, S.; Anaya-Acevedo, J.A.; Simard, M.; Liao, T.H.; Jaramillo, F. Analysis of floodplain dynamics in the Atrato River Colombia using SAR interferometry. Water 2019, 11, 875. [Google Scholar] [CrossRef] [Green Version]
  57. Estupiñan, A.R.C. Estudio de la variabilidad espacio temporal de la precipitación en Colombia. Doctoral Dissertation, Universidad Nacional de Colombia, Bogotá, Colombia, 2016. Available online: http://bdigital.unal.edu.co/54014/1/1110490004.2016.pdf (accessed on 5 October 2019).
  58. Guzmán, D.; Ruíz, J.F.; Cadena, M. Regionalización de Colombia Según la Estacionalidad de la Precipitación Media Mensual, a Través Análisis de Componentes Principales (ACP); IDEAM: Bogota, Colombia, 2014. [Google Scholar]
  59. Department of Civil and Environmental Engineering; Princenton-University Global Meteorological Forcing Dataset for Land Surface Modeling (GMFD). Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory. Available online: http://rda.ucar.edu/datasets/ds314.0/ (accessed on 3 May 2019).
  60. Sheffield, J.; Goteti, G.; Wood, E.F. Development of a 50-year high-resolution global dataset of meteorological forcings for land surface modeling. J. Clim. 2006, 19, 3088–3111. [Google Scholar] [CrossRef] [Green Version]
  61. Chaney, N.; Sheffield, J.; Villarini, G.; Wood, E.F. Development of a High-Resolution Gridded Daily Meteorological Dataset over Sub-Saharan Africa: Spatial Analysis of Trends in Climate Extremes. J. Clim. 2014, 27, 5815–5835. [Google Scholar] [CrossRef]
  62. National Oceanic and Atmospheric Administration; Extended Reconstructed Sea Surface Temperature (NOAA/ERSST) SST Data. Available online: https://www.esrl.noaa.gov/psd/data/gridded/ (accessed on 10 July 2018).
  63. Huang, B.; Thorne, P.W.; Banzon, V.F.; Boyer, T.; Chepurin, G.; Lawrimore, J.H.; Menne, M.J.; Smith, T.M.; Vose, R.S.; Zhang, H.M. Extended reconstructed Sea surface temperature, Version 5 (ERSSTv5): Upgrades, validations, and intercomparisons. J. Clim. 2017, 30, 8179–8205. [Google Scholar] [CrossRef]
  64. National Center for Environmental Prediction; National Center for Atmospheric Research (NCEP/NCAR) Monthly Gridded Atmospheric Data. NCEP/NCAR Reanalysis 1. Available online: https://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.html (accessed on 3 May 2019).
  65. Kalnay, E.; Kanamitsu, M.; Kistler, R.; Collins, W.; Deaven, D.; Gandin, L.; Iredell, M.; Saha, S.; White, G.; Woollen, J.; et al. The NCEP/ NCAR reanalysis project. Bull. Am. Meteorol. Soc. 1996, 77, 437–471. [Google Scholar] [CrossRef] [Green Version]
  66. Zhang, X.; Yang, F.; Canada, E. RClimDex (1.0) User Manual. Clim. Res. Branch Environ. Can. 2004, 2004, 1–23. [Google Scholar]
  67. Zhang, X.; Alexander, L.; Hegerl, G.C.; Jones, P.; Tank, A.K.; Peterson, T.C.; Trewin, B.; Zwiers, F.W. Indices for monitoring changes in extremes based on daily temperature and precipitation data. Wiley Interdiscip. Rev. Clim. Chang. 2011, 2, 851–870. [Google Scholar] [CrossRef]
  68. Folland, C.K.; Karl, T.R.; Jim Salinger, M. Observed climate variability and change. Weather 2002, 57, 269–278. [Google Scholar] [CrossRef]
  69. Sillmann, J.; Kharin, V.V.; Zwiers, F.W.; Zhang, X.; Bronaugh, D. Climate extremes indices in the CMIP5 multimodel ensemble: Part 2. Future climate projections. J. Geophys. Res. Atmos. 2013, 118, 2473–2493. [Google Scholar] [CrossRef]
  70. Aguilar, E.; Peterson, T.; Obando, P.R.; Frutos, R.; Retana, J.A.; Solera, M.; Soley, J.; García, I.G.; Araujo, R.M.; Santos, A.R.; et al. Changes in precipitation and temperature extremes in Central America and northern South America, 1961–2003. J. Geophys. Res. Atmos. 2005, 110, 1–15. [Google Scholar] [CrossRef]
  71. Santos, M.; Fragoso, M.; Santos, J.A. Regionalization and susceptibility assessment to daily precipitation extremes in mainland Portugal. Appl. Geogr. 2017, 86, 128–138. [Google Scholar] [CrossRef]
  72. Alexander, L.V.; Zhang, X.; Peterson, T.C.; Caesar, J.; Gleason, B.; Klein Tank, A.M.G.; Haylock, M.; Collins, D.; Trewin, B.; Rahimzadeh, F.; et al. Global observed changes in daily climate extremes of temperature and precipitation. J. Geophys. Res. Atmos. 2006, 111, 1–22. [Google Scholar] [CrossRef] [Green Version]
  73. Stephenson, T.S.; Vincent, L.A.; Allen, T.; Van Meerbeeck, C.J.; Mclean, N.; Peterson, T.C.; Taylor, M.A.; Aaron-Morrison, A.P.; Auguste, T.; Bernard, D.; et al. Changes in extreme temperature and precipitation in the Caribbean region, 1961–2010. Int. J. Climatol. 2014, 34, 2957–2971. [Google Scholar] [CrossRef]
  74. Debortoli, N.; Camarinha, P.; Marengo, J.; Rodrigues, R. An index of Brazil’s vulnerability to expected increases in natural flash flooding and landslide disasters in the context of climate change. Nat. Hazards 2017, 86, 557–582. [Google Scholar] [CrossRef]
  75. Ávila, A.; Justino, F.; Wilson, A.; Bromwich, D.; Amorim, M. Recent precipitation trends, flash floods and landslides in southern Brazil. Environ. Res. Lett. 2016, 11, 114029. [Google Scholar] [CrossRef]
  76. Kayano, M.T.; Rosa, M.B.; Rao, V.B.; Andreoli, R.V.; de Souza, R.A.F. Relations of the low-level extratropical cyclones in the southeast Pacific and South Atlantic to the Atlantic multidecadal oscillation. J. Clim. 2019, 32, 4167–4178. [Google Scholar] [CrossRef] [Green Version]
  77. Kayano, M.T.; Andreoli, R.V.; de Souza, R.A.F. Pacific and Atlantic multidecadal variability relations to the El Niño events and their effects on the South American rainfall. Int. J. Climatol. 2019, 1–18. Available online: https://doi.org/10.1002/joc.6326 (accessed on 26 December 2019). [CrossRef]
  78. Enfield, D.B.; Mestas-Nuñez, A.M. Multiscale variabilities in global sea surface temperatures and their relationships with tropospheric climate patterns. J. Clim. 1999, 12, 2734–2746. [Google Scholar] [CrossRef]
  79. Panofsky, H.G.; Brier, G.W. Some Applications of Statistics to Meteorology; Earth and Mineral Sciences Continuing Education, College of Earth and Mineral Sciences, Eds.; Pennsylvania State University: Pennsylvania, PA, USA, 1968; Available online: https://digital.libraries.psu.edu/digital/collection/digitalbks2/id/48274/ (accessed on 26 July 2019).
  80. Wilks, D.S. Statistical Methods in the Atmospheric Sciences, 3rd ed.; Wilks, D.S., Ed.; Academic Press: San Diego, CA, USA, 2011; ISBN 9780123850232. [Google Scholar]
  81. Jones, C.; Carvalho, L.M.V. The influence of the Atlantic multidecadal oscillation on the eastern Andes low-level jet and precipitation in South America. Npj Clim. Atmos. Sci. 2018, 1, 40. [Google Scholar] [CrossRef] [Green Version]
  82. Coe, M.T.; Latrubesse, E.M.; Ferreira, M.E.; Amsler, M.L. The effects of deforestation and climate variability on the streamflow of the Araguaia River, Brazil. Biogeochemistry 2011, 105, 119–131. [Google Scholar] [CrossRef]
  83. Costa, M.H.; Botta, A.; Cardille, J.A. Effects of large-scale changes in land cover on the discharge of the Tocantins River, Southeastern Amazonia. J. Hydrol. 2003, 283, 206–217. [Google Scholar] [CrossRef]
  84. Shapiro, A.S.S.; Wilk, M.B. An Analysis of Variance Test for Normality (Complete Samples). Biometrika 1965, 52, 591–611. [Google Scholar] [CrossRef]
  85. Anderson, T.W.; Darling, D.A. Asymptotic theory of certain goodness of fit criteria based on stochastic processes. Ann. Math. Stat. 1952, 23, 193–212. [Google Scholar] [CrossRef]
  86. Kucharski, F.; Parvin, A.; Rodriguez-Fonseca, B.; Farneti, R.; Martin-Rey, M.; Polo, I.; Mohino, E.; Losada, T.; Mechoso, C.R. The teleconnection of the tropical Atlantic to Indo-Pacific sea surface temperatures on inter-annual to centennial time scales: A review of recent findings. Atmos. 2016, 7, 29. [Google Scholar] [CrossRef] [Green Version]
  87. Araújo, R.G.; Andreoli, R.V.; de Souza, R.A.F.; Candido, L.A.; Kayano, M.T. Influence of EI Niño-Southern oscillation and equatorial Atlantic on rainfall over northern and northeastern regions of South America. Acta Amaz. 2013, 43, 469–480. [Google Scholar] [CrossRef] [Green Version]
  88. Kayano, M.T.; Oliveira, C.P.; Andreoli, R.V. Interannual relations between South American rainfall and tropical sea surface temperature anomalies before and after 1976. Int. J. Climatol. 2009, 29, 1439–1448. [Google Scholar] [CrossRef]
  89. Kayano, M.; Valéria Andreoli, R.; Ferreira de Souza, R.A. Evolving anomalous SST patterns leading to ENSO extremes: Relations between the tropical Pacific and Atlantic Oceans and the influence on the South American rainfall. Int. J. Climatol. 2011, 31, 1119–1134. [Google Scholar] [CrossRef]
  90. Choi, J.; An, S.-I.; Yeh, S.W. Decadal amplitude modulation of two types of ENSO and its relationship with the mean state. Clim. Dyn. 2012, 38, 2631–2644. [Google Scholar] [CrossRef]
  91. Li, X.; Xie, S.P.; Gille, S.T.; Yoo, C. Atlantic-induced pan-tropical climate change over the past three decades. Nat. Clim. Chang. 2016, 6, 275–279. [Google Scholar] [CrossRef]
  92. Barichivich, J.; Gloor, E.; Peylin, P.; Brienen, R.J.W.; Schöngart, J.; Espinoza, J.C.; Pattnayak, K.C. Recent intensification of Amazon flooding extremes driven by strengthened Walker circulation. Sci. Adv. 2018, 4, 1–7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. England, M.H.; McGregor, S.; Spence, P.; Meehl, G.A.; Timmermann, A.; Cai, W.; Gupta, A.S.; McPhaden, M.J.; Purich, A.; Santoso, A. Recent intensification of wind-driven circulation in the Pacific and the ongoing warming hiatus. Nat. Clim. Chang. 2014, 4, 222–227. [Google Scholar] [CrossRef]
  94. McGregor, S.; Timmermann, A.; Stuecker, M.F.; England, M.H.; Merrifield, M.; Jin, F.F.; Chikamoto, Y. Recent Walker circulation strengthening and pacific cooling amplified by Atlantic warming. Nat. Clim. Chang. 2014, 4, 888. [Google Scholar] [CrossRef] [Green Version]
  95. Espinoza, J.C.; Ronchail, J.; Marengo, J.A.; Segura, H. Contrasting North–South changes in Amazon wet-day and dry-day frequency and related atmospheric features (1981–2017). Clim. Dyn. 2019, 52, 5413–5430. [Google Scholar] [CrossRef]
  96. Losada, T.; Rodríguez-Fonseca, B.; Polo, I.; Janicot, S.; Gervois, S.; Chauvin, F.; Ruti, P. Tropical response to the Atlantic Equatorial mode: AGCM multimodel approach. Clim. Dyn. 2010, 35, 45–52. [Google Scholar] [CrossRef] [Green Version]
  97. Wang, C.; Enfield, D.B.; Lee, S.K.; Landsea, C.W. Influences of the Atlantic warm pool on western hemisphere summer rainfall and Atlantic hurricanes. J. Clim. 2006, 19, 3011–3028. [Google Scholar] [CrossRef]
  98. Wang, C.; Kucharski, F.; Barimalala, R.; Bracco, A. Teleconnections of the tropical Atlantic to the tropical Indian and Pacific Oceans: A review of recent findings. Meteorol. Z. 2009, 18, 445–454. [Google Scholar] [CrossRef] [Green Version]
  99. Global Precipitation Climatology Centre (GPCC) Monthly Gridded Precipitation Data. Available online: https://www.esrl.noaa.gov/psd/data/gridded/data.gpcc.html (accessed on 3 May 2019).
  100. Schneider, U.; Finger, P.; Meyer-Christoffer, A.; Rustemeier, E.; Ziese, M.; Becker, A. Evaluating the hydrological cycle over land using the newly-corrected precipitation climatology from the Global Precipitation Climatology Centre (GPCC). Atmos. 2017, 8, 52. [Google Scholar] [CrossRef] [Green Version]
  101. Marelle, L.; Myhre, G.; Hodnebrog, Ø.; Sillmann, J.; Samset, B.H. The Changing Seasonality of Extreme Daily Precipitation. Geophys. Res. Lett. 2018, 45, 352. [Google Scholar] [CrossRef] [Green Version]
  102. Betts, R.A.; Alfieri, L.; Caesar, J.; Feyen, L.; Gohar, L.; Koutroulis, A.; Lewis, K.; Morfopoulos, C.; Richardson, K.J.; Tsanis, I.; et al. Changes in climate extremes, fresh water availability and vulnerability to food insecurity projected at 1.5 °C and 2 °C global warming with a higher-resolution global climate model. Philos. Trans. A. Math. Phys. Eng. Sci. 2018, 2018, 1–27. [Google Scholar] [CrossRef]
  103. Tanoue, M.; Hirabayashi, Y.; Ikeuchi, H. Global-scale river flood vulnerability in the last 50 years. Sci. Rep. 2016, 6, 36021. [Google Scholar] [CrossRef] [PubMed]
  104. Marengo, J.A.; Valverde, M.C.; Obregon, G.O. Observed and projected changes in rainfall extremes in the Metropolitan Area of São Paulo. Clim. Res. 2013, 57, 61–72. [Google Scholar] [CrossRef] [Green Version]
  105. Li, X.; Gerber, E.P.; Holland, D.M.; Yoo, C. A Rossby wave bridge from the tropical Atlantic to West Antarctica. J. Clim. 2015, 28, 2256–2273. [Google Scholar] [CrossRef]
  106. Chikamoto, Y.; Timmermann, A.; Luo, J.J.; Mochizuki, T.; Kimoto, M.; Watanabe, M.; Ishii, M.; Xie, S.P.; Jin, F.F. Skilful multi-year predictions of tropical trans-basin climate variability. Nat. Commun. 2015, 6, 6869. [Google Scholar] [CrossRef] [Green Version]
  107. Sutton, R.T.; Hodson, D.L.R. Atlantic Ocean Forcing of North American and European Summer Climate. Science 2005, 309, 115–118. [Google Scholar] [CrossRef] [Green Version]
  108. Sutton, R.T.; Hodson, D.L.R. Climate response to basin-scale warming and cooling of the North Atlantic Ocean. J. Clim. 2007, 20, 891–907. [Google Scholar] [CrossRef]
  109. Levine, A.F.Z.; Frierson, D.M.W.; McPhaden, M.J. AMO forcing of multidecadal pacific ITCZ variability. J. Clim. 2018, 31, 5749–5764. [Google Scholar] [CrossRef]
  110. Jolliffe, I.T.; Cadima, J. Principal component analysis: A review and recent developments Subject Areas: Author for correspondence. Philos. Trans. R. Soc. A 2016, 374, 20150202. [Google Scholar] [CrossRef]
Figure 1. (a) Atrato River basin, (b) annual cycle of the precipitation, and (c) annual cycle of the streamflow during 1965–2015.
Figure 1. (a) Atrato River basin, (b) annual cycle of the precipitation, and (c) annual cycle of the streamflow during 1965–2015.
Water 12 00216 g001
Figure 2. Sea surface temperature (SST) mean standardized anomaly patterns in May–June (MJ) during the (a) cold Atlantic Multidecadal Oscillation (CAMO) and (b) warm Atlantic Multidecadal Oscillation (WAMO) and in November–December (ND) during the (c) CAMO and (d) WAMO. The continuous (dashed) line indicates significant positive (negative) values at the 95% confidence level using the Student’s t-test. The contour range is 0.1 standard deviations with red (blue) shades for positive (negative) values. The continuous (dashed) line encompasses significant positive (negative) values at the confidence level of 95%.
Figure 2. Sea surface temperature (SST) mean standardized anomaly patterns in May–June (MJ) during the (a) cold Atlantic Multidecadal Oscillation (CAMO) and (b) warm Atlantic Multidecadal Oscillation (WAMO) and in November–December (ND) during the (c) CAMO and (d) WAMO. The continuous (dashed) line indicates significant positive (negative) values at the 95% confidence level using the Student’s t-test. The contour range is 0.1 standard deviations with red (blue) shades for positive (negative) values. The continuous (dashed) line encompasses significant positive (negative) values at the confidence level of 95%.
Water 12 00216 g002
Figure 3. Time series of streamflow (Q) in m3·s−1 and Atlantic Multidecadal Oscillation (AMO) index smoothed with a 10-year running mean filter for (a) May–June (MJ) and (b) November–December (ND) during the 1965–2015 period. The blue and purple horizontal lines indicate the average streamflow before and after 1994, respectively.
Figure 3. Time series of streamflow (Q) in m3·s−1 and Atlantic Multidecadal Oscillation (AMO) index smoothed with a 10-year running mean filter for (a) May–June (MJ) and (b) November–December (ND) during the 1965–2015 period. The blue and purple horizontal lines indicate the average streamflow before and after 1994, respectively.
Water 12 00216 g003
Figure 4. Time series of streamflow (Q) in m3·s−1 and precipitation (PRP) in mm smoothed with a 10-year running mean filter in (a) May–June (MJ) and (b) November–December (ND) during the 1965–2015 period. The horizontal blue and purple lines indicate the average precipitation before and after 1994, respectively.
Figure 4. Time series of streamflow (Q) in m3·s−1 and precipitation (PRP) in mm smoothed with a 10-year running mean filter in (a) May–June (MJ) and (b) November–December (ND) during the 1965–2015 period. The horizontal blue and purple lines indicate the average precipitation before and after 1994, respectively.
Water 12 00216 g004
Figure 5. MJ standardized anomaly patterns associated with extreme flow events at Bellavista (left column) for the SST during: (a) CAMO, (b) WAMO, and (c) WAMO minus CAMO and for the SLP during (d) CAMO, (e) WAMO, and (f) WAMO minus CAMO. ND standardized anomaly patterns associated with extreme flow events at Bellavista (right column) for the SST during (g) CAMO, (h) WAMO, and (i) WAMO minus CAMO and for the SLP during (j) CAMO, (k) WAMO, and (l) WAMO minus CAMO. The contour range is 0.2 standard deviations. Shaded areas in red (blue) indicate positive (negative) values. The continuous (dashed) line encompasses significant positive (negative) values at the confidence level of 95%.
Figure 5. MJ standardized anomaly patterns associated with extreme flow events at Bellavista (left column) for the SST during: (a) CAMO, (b) WAMO, and (c) WAMO minus CAMO and for the SLP during (d) CAMO, (e) WAMO, and (f) WAMO minus CAMO. ND standardized anomaly patterns associated with extreme flow events at Bellavista (right column) for the SST during (g) CAMO, (h) WAMO, and (i) WAMO minus CAMO and for the SLP during (j) CAMO, (k) WAMO, and (l) WAMO minus CAMO. The contour range is 0.2 standard deviations. Shaded areas in red (blue) indicate positive (negative) values. The continuous (dashed) line encompasses significant positive (negative) values at the confidence level of 95%.
Water 12 00216 g005
Figure 6. MJ standardized anomaly patterns associated with extreme flow events at Bellavista (left column) for the 500 hPa vertical velocity in pressure coordinates during (a) CAMO, (b) WAMO, and (c) WAMO minus CAMO. ND standardized anomaly patterns associated with extreme flow events at Bellavista (right column) for the 500 hPa vertical velocity in pressure coordinates during (d) CAMO, (e) WAMO, and (f) WAMO minus CAMO. The contour range is 0.2 standard deviations. Shaded areas in red (blue) indicate positive (negative) values. The continuous (dashed) line encompasses significant positive (negative) values at the confidence level of 95%.
Figure 6. MJ standardized anomaly patterns associated with extreme flow events at Bellavista (left column) for the 500 hPa vertical velocity in pressure coordinates during (a) CAMO, (b) WAMO, and (c) WAMO minus CAMO. ND standardized anomaly patterns associated with extreme flow events at Bellavista (right column) for the 500 hPa vertical velocity in pressure coordinates during (d) CAMO, (e) WAMO, and (f) WAMO minus CAMO. The contour range is 0.2 standard deviations. Shaded areas in red (blue) indicate positive (negative) values. The continuous (dashed) line encompasses significant positive (negative) values at the confidence level of 95%.
Water 12 00216 g006
Figure 7. MJ standardized anomaly patterns associated with extreme flow events at Bellavista during CAMO (upper line), WAMO (middle line), and the difference WAMO minus CAMO (lower line) for the precipitation in the ARB (left column), R20 mm (middle column) and SDII (right column). Contour ranges are of 0.2 standard deviations for the precipitation of two days and one day for the R20 mm and R20 mm difference, and of 5 mm.day−1 and 2 mm day−1 for the SDII and the SDII difference, respectively. The continuous (dashed) line refers to the significant positive (negative) values at a 95% confidence level.
Figure 7. MJ standardized anomaly patterns associated with extreme flow events at Bellavista during CAMO (upper line), WAMO (middle line), and the difference WAMO minus CAMO (lower line) for the precipitation in the ARB (left column), R20 mm (middle column) and SDII (right column). Contour ranges are of 0.2 standard deviations for the precipitation of two days and one day for the R20 mm and R20 mm difference, and of 5 mm.day−1 and 2 mm day−1 for the SDII and the SDII difference, respectively. The continuous (dashed) line refers to the significant positive (negative) values at a 95% confidence level.
Water 12 00216 g007
Figure 8. As in Figure 7 but for November–December (ND).
Figure 8. As in Figure 7 but for November–December (ND).
Water 12 00216 g008
Table 1. Paired t-test for mean streamflow and precipitation before and after 1994 in the high-flow seasons discharge for MJ and ND. Q stands for the streamflow and PRP for precipitation.
Table 1. Paired t-test for mean streamflow and precipitation before and after 1994 in the high-flow seasons discharge for MJ and ND. Q stands for the streamflow and PRP for precipitation.
Q BellavistaPRP ARB
(m3·s−1)(mm·day−1)
BeforeAfterBeforeAfter
May-June
Mean25302800552586
1st quartile22302563482517
3rd quartile28732992624632
SD407376118110
Observations30213021
t-Statistic−2.411 −1.026
Two-tail0.02 0.310
November-December
Mean26153106472579
1st quartile21232673415528
3rd quartile29703379551640
SD557544121143
Observations30213021
t-Statistic−3.126 −2.899
Two-tail0.003 0.006
Table 2. Years with extreme flow events at Bellavista in May–June (MJ) and November–December (ND) before and after 1994.
Table 2. Years with extreme flow events at Bellavista in May–June (MJ) and November–December (ND) before and after 1994.
May–June (MJ)November–December (ND)
BeforeENSO *AfterENSOBeforeENSOAfterENSO
1966S-EL2000S-LN1973S-LN1999S-LN
1967N2003M-EL1975S-LN2005W-LN
1970W-EL2006N1984W-LN2007S-LN
1978W-EL2007N1988S-LN2010S-LN
1991N2011S-LN1990N2011M-LN
1993N 1993N
1994N 1994M-EL
* ENSO Type: EL—El Niño. LN—La Niña. W—Weak. M—Moderate. S—Strong. N—Neutral. The classification of ENSO types was based on the El Niño Oceanic Index (ONI), available at https://origin.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ONI_v5.php.

Share and Cite

MDPI and ACS Style

Loaiza Cerón, W.; Toshie Kayano, M.; Andreoli, R.V.; Avila, A.; Canchala, T.; Francés, F.; Ayes Rivera, I.; Alfonso-Morales, W.; Ferreira de Souza, R.A.; Carvajal-Escobar, Y. Streamflow Intensification Driven by the Atlantic Multidecadal Oscillation (AMO) in the Atrato River Basin, Northwestern Colombia. Water 2020, 12, 216. https://doi.org/10.3390/w12010216

AMA Style

Loaiza Cerón W, Toshie Kayano M, Andreoli RV, Avila A, Canchala T, Francés F, Ayes Rivera I, Alfonso-Morales W, Ferreira de Souza RA, Carvajal-Escobar Y. Streamflow Intensification Driven by the Atlantic Multidecadal Oscillation (AMO) in the Atrato River Basin, Northwestern Colombia. Water. 2020; 12(1):216. https://doi.org/10.3390/w12010216

Chicago/Turabian Style

Loaiza Cerón, Wilmar, Mary Toshie Kayano, Rita Valeria Andreoli, Alvaro Avila, Teresita Canchala, Félix Francés, Irma Ayes Rivera, Wilfredo Alfonso-Morales, Rodrigo Augusto Ferreira de Souza, and Yesid Carvajal-Escobar. 2020. "Streamflow Intensification Driven by the Atlantic Multidecadal Oscillation (AMO) in the Atrato River Basin, Northwestern Colombia" Water 12, no. 1: 216. https://doi.org/10.3390/w12010216

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop