Rainfall variability in the Brazilian northeast biomes and their interactions with meteorological systems and ENSO via CHELSA product

ABSTRACT Brazilian biomes are home to a significant portion of the world’s biodiversity, with a total of 14% of existing species and still concentrate 20% of the world’s water resources. However, changes in biomes have a direct impact on rainfall patterns and water recycling. Based on this, the objective was to evaluate the variability of rainfall in the four existing biomes in the Northeast Brazil (NEB) and their interaction with the ENSO climate variability mode and regional scale meteorological systems via CHELSA product. For this, monthly rainfall data were used from 1979 to 2013, with a spatial resolution of 1 km × 1 km of the CHELSA product, and seasonal and annual rainfall patterns were extracted via boxplot. It was found that the rainy season in the Amazon, Caatinga and Cerrado biomes occurred between January and April, with varying intensities, except for the Atlantic Forest. Such seasonality patterns are associated with the NEB meteorological systems, with emphasis on ITCZ (all Biomes), UTCV (Amazon, Caatinga and Cerrado), Frontal Systems (extreme south of Caatinga, Cerrado and Atlantic Forest) and EWD/TWD in the (Atlantic Forest). In the inter-annual scale, the remarkable influence of ENSO was verified, mainly in the years 1983, 1985, 1989, 1993, 1998, 2009 and 2012. It is noteworthy that 1985 was the wettest year of the period, with a surplus in all biomes, while the driest year differs between the Amazon (1983), Atlantic Forest and Caatinga (1993) and Cerrado (2012) biomes. The study via orbital product in NEB showed that anthropogenic processes and natural variability interfere with the forms of rain interception in the biomes and hence in rainfall patterns and water recycling in NEB.


Introduction
Brazilian biomes are shelter to a significant portion of the world's biodiversity (Agostinho, Thomaz, & Gomes, 2005), holding about 14% of all existing species (Lewinsohn & Prado, 2002). In the Atlantic Forest biome, the concentration of water resources is equivalent to 20% of the world's existing water resources, being one of the important centers of biodiversity due to the combination of high levels of wealth and endemism (Agostinho et al., 2005;Clemente, Oliveira Júnior, & Louzada, 2017;Sano, Rosa, Brito, & Ferreira, 2010;Scarano & Ceotto, 2015), however little explored (Agostinho et al., 2005;Scarano & Ceotto, 2015). In the last century, Holdridge (1947) proposed a classification of the world's plant formations using climate data from a triangular diagram. It showed the relationship between potential evapotranspiration (PET) and rainfall, annual precipitation and biotemperature. Therefore, the origin and evolution of the biome concept goes through the relationship between phytophysiognomy and climate variables, especially rainfall. (Coutinho, 2006).
One of the main challenges for maintaining biomes is related to recent changes in land use and occupation, especially the replacement of natural landscape areas by agriculture, livestock and urban development ( Roesch et al., 2009;Santos et al., 2018;Scarano & Ceotto, 2015). Another challenge is linked to the impacts of land use and occupation on climate change (Marengo, Williams, Alves, Soares, & Rodriguez, 2016a;Salazar, Nobre, & Oyama, 2007), which are a constant threat to the biodiversity preservation of Brazilian biomes, especially those predominant forestry and with greater species richness and endemism: Amazon and Atlantic Forest (Aleixo, Albernaz, Grelle, Vale, & Rangel, 2010;Clemente et al., 2017;Marengo et al., 2016a;Panisset et al., 2018;Salazar et al., 2007).
When it comes to the climate aspect, one of the factors responsible for the cyclical occurrence of these events is El Niño-Southern Oscillation (ENSO) variability mode, which influences the dynamics in the oceans and, consequently, alters large-scale circulation patterns (Kayano, Andreoli, & Souza, 2013;Kucharski, Polzin, & Hastenrath, 2008;Souza & Cavalcanti, 2009). This change in the circulation pattern affects the intensities and frequencies of regional and seasonal meteorological systems, and hence the rainfall on the biomes in northeastern Brazil (NEB) (Kayano et al., 2013;Kucharski et al., 2008;Marengo et al., 2016a).
To date, the assessment of rainfall dynamics in NEB biomes has been restricted to studies, considering only one biome (Barbosa et al., 2019;Martins & De Oliveira, 2011;Mataveli et al., 2017) or more than a biome at state level (Barros Correia Filho, Dos Santos, Diogo, & de Amorim, 2018). However, a regional study with all biomes and their relationship with rainfall via orbital product need to be investigated by the Brazilian and global scientific community. Therefore, the objective is to evaluate the variability of rainfall in the four existing NEB biomes and their interaction with the ENSO climate variability mode and regional meteorological systems via CHELSA product.

Study area
The evaluated biomes are located in the NEB region. The NEB is located between 01°02ʹ30" and 18°20ʹ07" S and 34°47ʹ30" and 48°45ʹ24" W. It is limited to the north (N) and east (E) with the Atlantic Ocean, to the south (S) with the states of Minas Gerais (MG) and Espírito Santo (ES), and to west (W) with the states of Pará (PA), Tocantins (TO) and Goiás (GO). The NEB is formed by nine states: Maranhão (MA), Piauí (PI), Ceará (CE), Rio Grande do Norte (RN), Paraíba (PB), Pernambuco (PE), Alagoas (AL), Sergipe (SE) and Bahia (BA). The NEB has an extension of 1,554 million km 2 , and an estimated population of 56.8 million (IBGE, 2019). In NEB are four biomes, they are: (a) Amazon; (b) Caatinga; (c) Cerrado; and (d) Atlantic Forest, illustrated in Figure 1.

NEB biomes
The Amazon biome is located in the state of Maranhão and has an estimated extension of 112.5 thousand km 2 , which represents 7% of the total NEB area (IBGE, 2002), at coordinates 0°47ʹ33"-05°37ʹ02" S and 43°37ʹ54"-48°53ʹ05" W. Amazon biome has great biodiversity and the existence of several conservation units (CUs), such as Gurupi Biological Reserve, Lençóis Maranhenses National Park, Mirador State Park and the newly created Chapada das Mesas National Park (Martins & De Oliveira, 2011).
The CUs are territorial spaces and their environmental resources, including jurisdictional waters, with relevant natural characteristics, legally established by the government aiming at conservation and with defined limits under special management regime, to which adequate guarantees of protection are applied (Brasil, 2000). Caatinga is located in the central NEB region, with an estimated area of 831.2 thousand km 2 representing approximately 53% of NEB and 7% of Brazilian territory. In Caatinga there are several economic activities geared towards farming and industrial purposes followed by a wealth of environments and species, with 932 species of plants, 178 mammals and 590 birds (MMA, 2019).
The Cerrado is the second largest biome in South America, located in central Brazil with an estimated area of 2.03 million km 2 , corresponding to 22% of the national territory. It is located in the states of Maranhão, Piauí and Bahia with an estimated area of 465.5 thousand km 2 (about 30% of the NEB area) (MMA -Ministério do Meio Ambiente, 2019). The Cerrado is recognized as savannah, being the most biodiverse biome in the world, sheltering 11,627 native plant species, of which approximately 4,400 are endemic species (existing only in this biome). In addition, it interacts with other biomes and has a phytophysiognomy (vegetation aspect of a region) diverse and unique (MMA -Ministério do Meio Ambiente, 2019).
The Atlantic Forest over NEB covers an area of approximately 164.7 thousand km 2 (about 10%). However, of this total, the biome occupies an area of approximately 19.4 thousand km 2 corresponding to 2.21% of its territory (Tabarelli, Melo, & Lira, 2006

Rainfall data -CHELSA
To evaluate the rainfall dynamics in the NEB biomes, we used rainfall data from the Climatologies at high resolution for the earth's land surface areas (CHELSA) product (Karger et al., 2017). CHELSA data have a spatial resolution of 1 km × 1 km, with monthly timescale. The study period was from 1979 to 2013.
CHELSA product consists of data from weather stations arising from the Global Precipitation Climate Centre (GPCC), with a temporal resolution of 50 km × 50 km. This data set is subjected to the downscaling statistical technique for grid refinement from 50 km to 25 km, using the B-Spline interpolation method with the help of Era-Interim precipitation data. In the case of rainfall, we assessed the effects of the topography, wind pattern and Atmospheric Boundary Layer (ABL) height, which are refined from 25 km to 1 km (Karger et al., 2017).
In addition, this procedure includes the wind effect and ABL height for a grid adjustment from 25 km to 1 km, whose data are available at: <http://chelsa-climate.org/timeseries. They are stored in Geotiff extension files. After obtaining these files, rainfall data were pre-processed for the boundary regions of the four existing NEB biomes. After selecting each biome, the files were submitted to the following steps: (1) organization; (2) manipulation; (3) extraction of monthly rainfall values; (4) obtaining descriptive and exploratory statistics for the biomes.
Based on the annual CV%, the Amazon ranged from 66.8% to 116.0%, the Caatinga ranged from 48.6% to 125.0%, the Cerrado ranged from 76.8% to 129.3%, and the Atlantic Forest ranged from 30.5% to 71.2%. Highlighting the Caatinga and Cerrado biomes, with CV% values greater than 100% in 10 of the 35 years (28.5% of the total), this indicated greater dispersion of rainfall data, and although the time series is heterogeneous (Wilks, 1995). It is noteworthy that the standard deviation was higher than the average, this result showed high discrepancy in the time series, which in turn is associated with two factors: i) data from orbital products and ii) high spatiotemporal variability of rain in the NEB, already verified by previous studies (Lyra, Oliveira-Júnior, & Zeri, 2014;Rao, Hada, & Herdies, 1995). While the Atlantic Forest showed the lowest variability in rainfall, with values below 50%, in 20 of the 35 years of the time series (57.14% of the total) (Correia Filho, Lucio, & Spyrides, 2016;Lyra et al., 2014).
In the monthly scale (Table 3), the rainy season occurrence of each of the existing biomes in the NEB stands out. In the Amazon it occurred from December to May, with Table 1. Moderate and strong intensity ENSO (aTSM ± 0.5°C) from 1979-2017 (Huang et al., 2015).
Regarding seasonality (Table 4), the wettest quarter is identified, the months of March-April-May (MAM) in all NEB biomes, except for the Cerrado, whose quarter was December-January-February (DJF). The average rainfall for Amazonia, Caatinga and Atlantic Forest were 302.8 mm ± 128.4 mm, 107.6 mm ± 65.4 mm and 144.4 mm ± 56.1 mm, respectively. In the case of the Cerrado, it was 134.3 mm ± 82.2 mm. CV% in Amazonia was 42.4% for the MAM quarter, followed by 43.6% in Caatinga and Floresta Atlantica was 58.6% and 34.7% in the June-July-August (JJA) quarter, respectively. In the Cerrado, it was 43.6% in the quarter DJF.

Rainfall seasonality in NEB biomes
Seasonality of rainfall in NEB biomes via monthly boxplot showed high rainfall variability between biomes, especially in relation to median and interquartile range. Clearly, the Amazon (1847.9 mm ± 418.7 mm) (Figura 2(a)), Caatinga (773.1 mm ± 212.5 mm) (Figura 2(b)) and Cerrado (1196.6 mm ± 218.9 mm) (Figura 2(c)) biomes showed that the rainy season occurred between January and April, but with varying intensities. The exception was the Atlantic Forest biome (1363.8 mm ± 217.0 mm) (Figure 2(d)), with the rainy season between March and July. Outliers were remarkable in all biomes, especially in the rainy season, except for the Amazon biome.
It is noteworthy that in the Caatinga, Cerrado and Atlantic Forest biomes, the presence of outliers is associated with rainfall above the climatological average from CHELSA data. The months considered atypical were between January and March. When referring to the dry season in Amazon, Caatinga and Atlantic Forest biomes, it is between August and October, while in the Cerrado biome occurs between June and August (Nascimento, De Araújo, & Ferreira Júnior, 2011). During this dry season in NEB, mainly in the northeastern Amazon and Cerrado biomes, there is a significant increase in the number of fire outbreaks, mainly in the states of Maranhão. (Bezerra et al., 2018), Piauí and Bahia (Nascimento et al., 2011), from natural origin and/or anthropic activities (Bezerra et al., 2018;Caúla, Oliveira-Júnior, Lyra, Delgado, & Heilbron Filho, 2015;Mataveli et al., 2017;Sales et al., 2019).
Such a pattern of rainfall variability identified between biomes is directly associated with rainfall producing and inhibiting meteorological systems, which in turn act at distinct periods along the NEB (Lyra et al., 2014). It is noteworthy that NEB has a high interannual and intra-seasonal rainfall variability, with the thermal influence of the Pacific and Tropical Atlantic oceans, which control and modulate years of drought and flooding in the region (Alves et al., 2017;Correia Filho et al., 2016;Lyra et al., 2017), and therefore will drastically interfere with the rainfall dynamics of biomes.
Finally, Caatinga and Cerrado, as they are located in much of the NEB, both biomes have similar rainfall characteristics; however, with intensities different, but both biomes are influenced by the ITCZ in the northern portion, while the far south is due to forays from the frontal systems that detach from the South Atlantic Convergence Zone (SACZ) (Kousky, 1979;Lyra et al., 2014). In addition to SACZ and ITCZ, there is also the UTCV acting in this period Fedorova et al., 2017;Kousky & Alonso Gan, 1981).

Interannual rainfall variability in NEB biomes
In the annual analysis of rainfall in the NEB biomes (Figure 3), there is a highly variable and similar cycle in the annual rainfall pattern common to all biomes, but with marked differences in their intensities. Regarding the annual mean totals, the highest rainfall was found in the Amazon compared to the other biomes with an annual mean around 1847.9 mm year −1 . Annual rainfall totals identified in the Amazon biome are maintained by ITCZ in NNEB (Alves et al., 2017;Araújo et al., 2009;Reboita & Santos, 2015; and still associated with the process of evapotranspiration (ET) of the biome that contributes to the release of water vapor in the atmosphere and, consequently, causing the formation of rain-producing clouds (Fisch, Marengo, & Nobre, 1998;Salati, Dall'Olio, Matsui, & Gat, 1979;Xu, Agee, Wang, & Ivanov, 2019).
Regarding Caatinga, similar to the central pattern of the Northeastern Cerrado, UTCV influences the rainfall pattern in the region Fedorova et al., 2017;Kousky & Alonso Gan, 1981). However, the low annual totals recorded are associated with their geographical position, situated in the northeastern semiarid region which has an irregular rainfall regime mainly motivated by the orography (Correia Filho et al., 2018;Reboita et al., 2016).
In this case, this occurs due to the positioning of orographic barriers, such as the Borborema Plateau and Diamantina and Araripe Plateaus (Correia Filho et al., 2018;Reboita et al., 2016), which prevent wind flow and humidity transport from the Atlantic Ocean to the interior of the NEB, thus inhibiting rainfall to the leeward of these orographic barriers.
During the ENSO (El Niño) warm phase, there is a significant decrease in rainfall and prolonged drought regimes in NEB, especially in the northeastern semiarid region (Kayano et al., 2013;Marengo & Hastenrath, 1993;Marengo et al., 2016b;Rodrigues et al., 2011), where Caatinga is located. In addition, there is also a considerable increase in the number of fire records in the Amazon biome (Caúla et al., 2015;Fearnside, 2009). Regarding the Atlantic Forest, De Oliveira Souza et al. (2018) found a decrease in vegetation cover along the NEB during ENSO episodes (except in Alagoas and Pernambuco); in addition, there was a reduction associated with urban sprawl and areas susceptible to desertification. Filgueira, Neto, and Barbosa (2007) pointed out that the El Niño periods are marked by the record of low rainfall, preventing the outflow generation and the recovery of volumes stored in the dams, and mark the progressive decrease of the water mirrors, reaching the total emptying of the reservoirs.
During the ENSO (La Niña) cold phase, GISST, when associated with TAV negative in the southern Atlantic, affects the climatic positioning of the ITCZ, shifting it south of its climatic position (Reboita & Santos, 2015;Souza & Cavalcanti, 2009). This results in an increase in rainfall patterns and better distribution throughout the NEB (Kayano et al., 2013;Rodrigues et al., 2011). However, during the 2011-2012 La Niña occurrence, the negative Sea Surface Temperature Anomaly (SSTA) pattern in the combined Atlantic and Pacific Oceans resulted in the ITCZ displacement to north from its climatological position, which contributed to significantly decreased rainfall rates over the NEB (Rodrigues & McPhaden, 2014;Souza & Cavalcanti, 2009) which in turn especially affected the semiarid region and especially the Caatinga (Pereira, Justino, Malhado, Barbosa, & Marengo, 2014). Such effects have been devastating to NEB biomes, especially the economic sectors (Marengo et al., 2016b), and social and environmental problems (Marengo et al., 2018;Mariano et al., 2018;Pereira et al., 2014).
In addition, with the reduction of rainfall indices increased the biome ET, thus contributing to the reduction of soil moisture (Carrão, Russo, Sepulcre-Canto, & Barbosa, 2016) and hence aggravating the degradation of land use and occupation (Marengo & Bernasconi, 2015), and still contributing to the severity of the drought (Carrão et al., 2016;Marengo et al., 2018;Marengo & Bernasconi, 2015).

Actuation of ENSO phases on rainfall in NEB biomes
From the previous results (Figure 3), four episodes of La Niña (1985, 1989, 2009 and 2012) and two episodes of El Niño (1983, 1993 and 1998) were selected for an intra-annual and annual total analyses, as shown in Figure 4 and Table 2. Comparing the monthly rainfall series and the climatological mean, it is verified that the anomalous series presented rainfall above (below) the climatological mean between the months of October and December (January and August), except for the years 1985, 1989 and 2009, in the occurrence of La Niña. In relation to the annual totals, the years mentioned had excess rainfall which in turn contributed significantly to the maintenance of the biomes in NEB, especially for Caatinga in 1985, with a gain of 645.0 mm year −1 (83.4% above the climatological mean).
This situation was also verified in the Amazon and Cerrado biomes with a gain of 761.2 mm year −1 (41.2% above the climatological mean) and 712.6 mm year −1 (59.6% above the climatological mean). In the case of the Atlantic Forest this gain was modest, only 219.5 mm year −1 (16.1% above the climatological mean). These rainfalls above the climatological mean occurred in 1985 as a result of a La Niña categorized as weak according to the ONI index combined with positive TAV in the southern Atlantic. Both factors resulted in high NEB rainfall, with flooding records in several locations, mainly in the capitals of Piauí and Maranhão (Andrade et al., 2018;Chaves, Tavares, & de Andrade, 2017;De Oliveira, Guedes, & da Silva Costa, 2017;Pereira, Soares, Silva, de Assunção Montenegro, & de Souza, 2017). In 1983In , 1993 years there was an annual loss of rainfall according to the El Niño episodes which in turn resulted in various impacts on biomes and various economic sectors (Confalonieri, 2015;Marengo et al., 2018Marengo et al., , 2016bRao et al., 1995). Especially in 1983, where there was a significant reduction in total rainfall (Table 5) with a loss of 816.9 mm year −1 for the Amazon biome, with an annual mean of 1037 mm.year −1 . The largest percentage reduction was observed in Caatinga in 1993, with 49.5% below the climatological mean (Rao et al., 1995).
However, among the episodes with the highest annual rainfall, the largest was in 2012, considered atypical because it is La Niña year. This episode of La Niña was extremely severe, especially for Caatinga, where there was a loss of 359.9 mm year −1 , a 46.5% reduction in annual rainfall. This effect resulted in a huge $ 6 billion loss in the agricultural  1985,1989,1993,1998,2008  sector, and 997 of the 1794 NEB municipalities declared a state of emergency due to drought severity (Marengo et al., 2018).

Spatial distribution of rainfall in NEB biomes
Spatial distribution of rainfall on the dry and rainy period in the NEB biomes followed by the driest and wettest year of each biome in the period 1979-2013 were evaluated as shown in Figure 5-8. In common, all biomes presented the year as the wettest year, the Table 5. Climatological mean and annual accumulated in mm/year (percentages in %) of the years 1983, 1985, 1989, 1993, 1998, 2009 and 2012,   driest year for the Amazon biome ( Figure 5) was 1983, for the Caatinga (Figure 6) and the Atlantic Forest (Figure 8) was 1993, and for the Cerrado (Figure 7) was the year 2012. It is noteworthy that the biomes are strongly influenced by the rain-producing systems, previously evaluated. Moreover, there were marked differences in the drier years.
In the Amazon Biome, it is found that rainfall decreases significantly in the south portion, especially in the dry season. Regarding spatial climatology, the dry period ( Figure 5(a)) presented a well-defined homogeneous rainfall along the biome, ranging from 500 to 775 mm, while in the rainy season ( Figure 5(d)), the largest accumulations occur in the northern portion, with values between 1000 and 2040 mm. This spatial variability is due to the performance of the ITCZ and squall lines (Alves et al., 2017;Araújo et al., 2009;Reboita & Santos, 2015;. When analyzing the wettest year ( Figure 5(b,e)) and the driest year ( Figure 5(c,d)), again prevails the pattern defined by climatology, but with varying and distinct intensities, for example, in the wettest (dry) year, the northern portion was above the climatological mean with a maximum of 1193 mm (rainfall below the climatological mean between 112 and 521 mm) in the rainy season.
The Caatinga biome showed a rainy pattern (Figure 6(d)) similar to the Amazon biome with distinct intensities, where the largest accumulated occurred in the coastal range located in the northern portion, with values between 800 and 1487 mm. While in the dry season ( Figure 6(a)), the eastern portion is influenced by EWD/TWD with values greater than 700 mm, while most of the biome values do not exceed 300 mm. When the wettest year ( Figure 6(b,e)) and the dry year ( Figure 6(c,f)) are evaluated, respectively, the difference is quite distinct between them and the climatology. During the wettest year, the Araripe Plateau region recorded twice the climatological mean for the rainy season ( Figure 6(e)) with values higher than 2500 mm, and in the northwest with a maximum of 3420 mm. In the driest year (Figure 6(f)), the maximum occurred in the extreme northwest and north portions of the biome with maximum values reaching 1000 mm, while in much of the region, the values were below 350 mm and aggravated during the dry period with values below 200 mm. Such rainfall variability had a negative impact on the region, with significant agricultural and economic losses (Marengo et al., 2018(Marengo et al., , 2016b. Regarding the Cerrado biome (Figure 7), the differences in the climatic pattern between the dry (Figure 7(a)) and rainy (Figure 7(d)) periods are related to intensity and influence area. In the dry season (Figure 7(a)), the rainfall is concentrated on the northern end of the biome with values between 45 and 200 mm, while in the rainy period (Figure 7(d)) was well distributed with its maximum in the western portion of the biome, ranging from 1000 to 1757 mm. In relation to the wettest year (Figure 7(b,e)), some regions in this biome double the rainfall rates in the rainy season, especially in the north with values exceeding 2500 mm and a maximum of 3210 mm, and decreasing towards the south, with a minimum of 590 mm. In the driest year, in relation to the dry season (Figure 7(c)) did not reach the value of 400 mm while the rainy season (Figure 7(f)) did not exceed 1540 mm.
The Atlantic Forest biome ( Figure 8) showed a very distinct pattern between the dry and rainy periods (Figure 8(a,d)) with values between 141 and 1893 mm and 185 and 993 mm, respectively. These patterns are influenced by the influence of EWD/TWD and frontal systems, respectively (Alves et al., 2001;Gomes et al., 2019;Kouadio et al., 2012;Lyra et al., 2014). As for the wettest (Figure 8(b,e)) and driest (Figure 8(c,f)) years, the variations in the rainfall indices were significant, being more intense in the wettest year, especially in the coastal region, during the rainy season (Figure 8(e)) with values between 1350 and 2430 mm. While in the dry season (Figure 8(b)), the rainfall was concentrated on the extreme south portion, due to the frontal systems (Kousky, 1979) which in turn ranged from 400 to 974 mm.

Conclusions
Rainfall dynamics in the NEB via CHELSA product shows that the rainy season in the biomes occurs between January and April, with varying intensities, except for the Atlantic Forest where this occurs between April and August. Rainfall seasonality in NEB biomes is associated with spatiotemporal variability of meteorological rain-producing systems, for example, the ITCZ (all Biomes), UTCV (Amazon, Caatinga and Cerrado), Frontal Systems (far south of Caatinga, Cerrado and Atlantic Forest) and EWD/TWD in (Atlantic Forest). Regarding the inter-annual rainfall variability in the NEB biomes, the influence of the ENSO climate variability mode is clearly shown, being found mainly in specific years of the time series, such as 1983, 1985, 1989, 1993, 1998, 2009 and 2012. In relation to the years mentioned above, we highlight 1985 (the wettest of the time series), with a surplus in all evaluated biomes. The driest year stands out in the time series among the biomes, for example, on Amazon in 1983, Caatinga and Atlantic Forest in 1993, and Cerrado in 2012.
On the spatial scale, rainfall in each biome varies significantly as a function of the positioning and action of rainfall producing weather systems, followed by ENSO phases. The study via CHELSA product in NEB brings a discussion about rainfall dynamics in the existing NEB biomes. However, if constant changes in biomes by land use and occupation, forest fires and burns, desertification process and prolonged droughts continue along with natural variability, forms of rain interception in biomes will impact rainfall patterns and water recycling in NEB.