Variations of Rainfall Rhythm in Alto Pardo Watershed , Brazil : Analysis of Two Specific Years , a Wet and a Dry One , and Their Relation with the River Flow

This research aims to understand the variability and rhythm of rainfall for two specific standard-years, and their relation with the river flow of the Alto Pardo watershed, located in southeastern Brazil, and thus identify atmospheric systems that can cause extreme events, and which may be reflected in heavy rainfall, floods, or drought episodes. Therefore, the research chose to investigate the years 1983 and 1984, rainy and dry standard-years respectively in the study area, where rainfall was described and spatialized through the geostatistical method of kriging at the monthly level and the rhythmic analysis technique was applied in order to identify what weather types are usual and extreme in the area. The results indicate that a high involvement of the frontal system in the year 1983 was responsible for the episodes of greater rainfall and peak water flow, especially in stationary front episodes. The year 1984 presented low rainfall in summer, a meteorological drought during the year, and the predominance of tropical air masses in relation to the frontal systems. The comparison between the two extreme years, a wet and a dry one, made it possible to understand the frequency and the chaining of the atmospheric systems during this period for the Alto Pardo watershed.


Introduction
Rainfall variability is really important at the regional and local scales, due to its ability to affect human society.At this level of detail, although the atmospheric conditions are usually expected most of the time, extreme events can cause weather conditions that the resident population is unprepared for and can cause severe rains [1,2], floods [3,4], or episodes of drought [5][6][7].Therefore, in knowing various impacts caused by the rain, there is little doubt that society as a whole has become more vulnerable to extreme weather conditions [8,9].
Such extreme conditions of weather and climate occur throughout the entire world, however, in recent decades, the intense concern about the imminence of a possible increase in the occurrence of extreme events has grown [10,11].The possibility of affecting more people on the planet has motivated researchers to understand climate variability and dynamics and the climate elements which can be responsible for these impacts, and their frequency and intensity in the atmosphere.For example, we can mention the studies of [12][13][14][15].
Atmospheric circulation is central to understanding the variability of the regional climate, because it presents a dynamic linked to the succession of different atmospheric systems on this scale.Monteiro [16], seeking to connect the dynamics of the types of weather to the other geographical has contributed to changes in the regional landscape by agricultural activity and has demanded a large amount of water supply for irrigation.

Study Period, Data, and General Analysis
Based on the principle of this study, which seeks to investigate through the rhythmic analysis the weather types that can result in extreme rainfalls, floods, or droughts, the first step was to identify and choose a period with extreme weather conditions.Therefore, we analyzed previous studies that had identified and classified years of extreme conditions based on the standard-year method [45,46].These surveys identified a great contrast between the years 1983 and 1984, the first one being designated as an extremely wet yet and the second classified as a dry year in the Pardo watershed.Another factor that emphasizes the choice of this period (1983)(1984) was the occurrence of a very strong El Niño event in southern hemisphere in 1983 [47], which caused an extremely wet year in the South and Southeast of Brazil [48].
To support the choice of the years 1983 and 1984, a general description of the area was carried out to verify the behavior of such years from a long series of data.For this, we collected daily rainfall and natural water flow data of Caconde city, located in the center of the basin.The pluviometry station belongs to the Department of Water and Electrical Energy [49], a São Paulo state agency, and the fluviometric data were generated by the National Electric System Operator [50].It is important to emphasize that we chose to use the natural flow data because the agency service provided model reconstructs the original water flow of the Pardo River, identifying the data that would occur if there were no human action, as the water stored in the reservoir and its use in agriculture and urban water supply [51].This type of data has made it possible to analyze the daily flow of the river without having to take into account the expansion of the use of water resources in this stage of the research.
The general analysis of the data was made using the historical series between 1966 and 2015, a 50 year period.The choice of this time series was based on the data availability, the date of construction of the hydroelectric plant, and the occurrence of possible gaps in the daily rainfall series.In the case of gaps in the pluviometry data, we calculated the average of the daily data of the three nearest stations [52].The bar graphs were prepared using Excel software and the boxplot graphs were generated by the statistics software R version 3.2.2..

Spatial Interpolation of Rainfall
For the mapping of the geostatistical rainfall interpolation, we use the Kriging technique [53] through the Surfer 13 software.This procedure allows for the continuity of the phenomena in space, even though gaps were left between the sampling points.The monthly rainfall data used in this study were collected through the Hidroweb tool, administered by the National Water Agency (ANA) of Brazil.In total, 18 pluviometric stations were used to perform the rainfall mapping,

Study Period, Data, and General Analysis
Based on the principle of this study, which seeks to investigate through the rhythmic analysis the weather types that can result in extreme rainfalls, floods, or droughts, the first step was to identify and choose a period with extreme weather conditions.Therefore, we analyzed previous studies that had identified and classified years of extreme conditions based on the standard-year method [45,46].These surveys identified a great contrast between the years 1983 and 1984, the first one being designated as an extremely wet yet and the second classified as a dry year in the Pardo watershed.Another factor that emphasizes the choice of this period (1983)(1984) was the occurrence of a very strong El Niño event in southern hemisphere in 1983 [47], which caused an extremely wet year in the South and Southeast of Brazil [48].
To support the choice of the years 1983 and 1984, a general description of the area was carried out to verify the behavior of such years from a long series of data.For this, we collected daily rainfall and natural water flow data of Caconde city, located in the center of the basin.The pluviometry station belongs to the Department of Water and Electrical Energy [49], a São Paulo state agency, and the fluviometric data were generated by the National Electric System Operator [50].It is important to emphasize that we chose to use the natural flow data because the agency service provided model reconstructs the original water flow of the Pardo River, identifying the data that would occur if there were no human action, as the water stored in the reservoir and its use in agriculture and urban water supply [51].This type of data has made it possible to analyze the daily flow of the river without having to take into account the expansion of the use of water resources in this stage of the research.
The general analysis of the data was made using the historical series between 1966 and 2015, a 50 year period.The choice of this time series was based on the data availability, the date of construction of the hydroelectric plant, and the occurrence of possible gaps in the daily rainfall series.In the case of gaps in the pluviometry data, we calculated the average of the daily data of the three nearest stations [52].The bar graphs were prepared using Excel software and the boxplot graphs were generated by the statistics software R version 3.2.2..

Spatial Interpolation of Rainfall
For the mapping of the geostatistical rainfall interpolation, we use the Kriging technique [53] through the Surfer 13 software.This procedure allows for the continuity of the phenomena in space, even though gaps were left between the sampling points.The monthly rainfall data used in this study were collected through the Hidroweb tool, administered by the National Water Agency (ANA) of Brazil.In total, 18 pluviometric stations were used to perform the rainfall mapping, covering the Alto Pardo watershed and surrounding areas.Table 1 presents in detail the individual information on each of the stations.The choice of the pluviometric stations was made through the existence of a rain gauge and the rainfall data availability for the period.The spatial coverage of the seasons can be observed in Figure 2, highlighting the 18 distribution points throughout the study area, associated with the regional orography.In the figure, it is possible to notice a good overall distribution of the rainfall stations, but it shows a lower density of stations in the southeast region of the basin, coinciding with a high altitude area and a low urbanization level, factors that contributed to the minor presence of the rain gauges.However, the existence of stations near the study area aims to compensate for, even partially, the absence of data for the interpolation process.Thus, the spatial interpolation of the rainfall monthly level during the years 1983-1984 was made seeking to provide the spatial rainfall distribution of the area, which during the qualitative analysis of the results was related to the frequency of the atmospheric systems that have passed through the region at the monthly level.The choice of the pluviometric stations was made through the existence of a rain gauge and the rainfall data availability for the period.The spatial coverage of the seasons can be observed in Figure 2, highlighting the 18 distribution points throughout the study area, associated with the regional orography.In the figure, it is possible to notice a good overall distribution of the rainfall stations, but it shows a lower density of stations in the southeast region of the basin, coinciding with a high altitude area and a low urbanization level, factors that contributed to the minor presence of the rain gauges.However, the existence of stations near the study area aims to compensate for, even partially, the absence of data for the interpolation process.Thus, the spatial interpolation of the rainfall monthly level during the years 1983-1984 was made seeking to provide the spatial rainfall distribution of the area, which during the qualitative analysis of the results was related to the frequency of the atmospheric systems that have passed through the region at the monthly level.

Rhythmic Analysis
To generate the rhythmic analysis, a traditional technique of Brazilian Geographic Climatology [54,55], we used weather element graphs, synoptic charts, and satellite images to understand the atmospheric circulation and weather types on a regional scale.Thus, we collected daily and hourly data of the Graminha weather station, belonging to the Department of Water and Electrical Energy [49], located at 21 • 34 south latitude, 46 • 37 west longitude, and at 880 m altitude in Caconde city.The data collected by the weather station refers to the following climatic elements: precipitation, air temperature, atmospheric pressure, air humidity, duration of sunshine, wind direction and intensity, and cloudiness.Satellite images were collected in daily and hourly levels through the International Satellite Cloud Climatology Project (ISCCP) of the Global ISCCP B1 Browse System (GIBBS) [56], and the synoptic charts were provided by the Brazilian Navy [57].
The rhythmic analysis was performed to identify the atmospheric circulation through the manual method, and the synoptic weather classification in daily and hourly data was applied to identify the atmospheric systems that are occurring in the period and its circulation over time.The classification of the atmospheric system was made in two daily schedules in Brazil: 9 a.m.(12 GMT) and 9 p.m. (00 GMT).It is noteworthy that a manual method, despite the criticism it received in the 21st century due to the high degree of subjectivity involved [58], brings the researcher closer to understanding the atmospheric circulation and dynamics in the region during the research, so it is still valid for use [59], even with the introduction of the automated synoptic classification and statistically sophisticated methods [60].Regardless of the choice of method adopted in the study, all the classifications are inherently subjective to some degree [61], and these influences require the expertise of the researcher.In the case of South America, the non-existence of automated methods widely accepted for classifying "weather types" reinforces the use of the manual method in this study.
The identification of atmospheric systems in the area, at the regional level, adopted a classification system based on the relative positions of the air masses and fronts in South America [62,63], through the atmospheric flows of the anticyclones, cyclones, and frontal passages, as well as the response of climatic elements on a smaller-scale.In order to perform this procedure, we constructed the rhythmic analysis charts through the generation of line graphs of all the climatic variables collected at the Graminha weather station at the daily level, which have been aligned and organized in sequence.The graphs facilitate the analysis of variability of the data and can be observed in [41].After that, the synoptic charts and satellite images were analyzed for each day of the year, always associating the day-to-day synoptic weather conditions with the variations of the climate data in the region, to identify the dominant atmospheric systems in the two daily schedules.This genetic classification adopted the same methodology in several other studies in Brazil [64][65][66].The classes of the atmospheric systems and the basis of the weather types are as follows: 1.
Air masses with tropical characteristics- The systems identified above are only those that occur frequently in southeastern Brazil.During the years 1983-1984, there were no cases of Equatorial Continental mass (mEQ), Occlude Polar Atlantic Front (OCL), and Atlantic Intertropical Convergence Zone in the Alto Pardo watershed.Regarding the origin of the systems, their synthesis [65] can be classified as follows:
The air masses that can be observed in the southeastern region of Brazil, especially in the region of the study area, can be classified both by the tropical source region (mTA, mTAC, and mTC) and the polar source region (mPA, mPV, and mPVC).Table 2 presents a brief description of the air masses' region of origin and explains how these atmospheric systems are classified, based on criteria to be observed and used in the synoptic charts and the variability of the climatic data for the classification of the atmospheric systems.Figure 3a illustrates the distribution of air masses by the source region, and their main incursions into the study area.The presence of isobars of the high pressure sector influencing the study area, through the air flows of east and northeast observed in the synoptic chart, reinforced by the tropical characteristics of the climatic variables, such as high temperatures and relative humidity.

mTAC
This system is formed from the mTA, when it remains over the continent for a few days and although it still has a high pressure, it loses its original properties, mainly decreasing the relative humidity.
Always preceded by mTA and presenting a pressure pattern similar to this air mass, it is distinguished by the significant change in the climatic variables that can be observed in the data of the study area, such as the great decrease in relative humidity, higher maximum temperature, clear sky, and high insolation.

mTC
Low pressure system of Chaco, in the Tropic region east of the Andes, with cyclonic circulation of the surface and anticyclonic at the upper levels, as a consequence of the intermittent thermal-orographic depression.The characteristics of this air mass is hot and dry.
The presence of isobars of the low pressure sector of the center of South America in the study area, coming from the air flows to the west observed in the synoptic chart.A large decrease in relative humidity, elevation of maximum temperature, clear sky, and high insolation are characteristics of this air mass.

mPA
Anticyclone with the source region in the Atlantic Ocean in southern South America, presenting cold temperatures as the main characteristic.Because of the high pressure, this air mass tends to flow towards the lower latitudes.
It always occurs after the passage of a cold front, when the incursion of this anticyclonic system happens.In addition to observing the high pressure influencing the region in the synoptic charts, another factor that should be noted is the weather data, mainly the presence of lower temperatures and high atmospheric pressure for the area.

mPV
This system is formed from the mPA, when it remains outside the source region for a few days and loses its original properties, mainly raising its temperature.
In the synoptic chart one must observe the incursion and the influence of the polar anticyclone on the study area.However, it differs clearly from mPA by the sensible decrease of the original atmospheric pressure, in addition to the higher temperature when acquiring the tropical characteristics.

mPVC
Similar to mPV, this system is formed from mPA when it loses its original properties; however, it was modified by the continental trajectory of the anticyclone migration, mainly decreasing the relative humidity.
It presents the same criteria of mPV, however it is classified as mPVC when the anticyclone trajectory occurs inside the continent, resulting in very low relative humidity.
The incursion of mid-latitude polar air towards lower latitudes is linked to the cold fronts, whose trajectory and movement exhibit a different behavior throughout the year.In general, the origin of the frontal systems in the region of study is found in the advances of the mPA towards the Brazilian Southeast, generally associated with cloudy weather, high winds, and precipitation.During autumn and winter, the penetration of polar air is more intense in the Brazilian southeast, while in spring and summer, this incursion presents greater resistance by the tropical air masses [67].The different movements of these systems are identified by the synoptic charts through their symbols, and the criteria for classification can be observed in Table 3.In addition, the table provides information on the Tropical Instability Line (LI), identified as troughs, and is associated with the air masses, not being individually classified as an atmospheric system.Figure 3b

Type
Criteria for Classification1

FPA
The system is represented in the synoptic chart as the limits of mPA's advance toward the study area.
A drop in temperature and atmospheric pressure should be noted in the climatic variables, in addition to an increase in the relative humidity and cloudiness.

FPR
This system occurs when, after the polar air has advanced over the region of study, a trough is observed in this anticyclone forming a squall line, which can be identified in the synoptic chart by the symbol of frontolysis in the continent and the symbol of frontogenesis in the ocean.Increased humidity and cloudiness can be observed in the climate data.

DIS
It is represented in the synoptic chart as a frontolysis over the study area.In the climatic variables, an increase of the sunshine and decrease of the humidity and cloudiness can be observed during the day.

REP
This system is classified when an approximation of the FPA is observed through the synoptic chart, but it does not occur directly over the study area.Even so, a significant change has already been observed in the climatic variables such as increased cloudiness and humidity, decreased thermal amplitude and, in most cases, precipitation.

EST
This system has a specific symbol for stationary front in the synoptic chart, and can be classified when a frontal system is moving very slowly or is stalled for a few days.Mostly cloudiness and precipitation are observed.

QTE
This system can be identified as a warm front in the synoptic chart.It is usually associated with increased cloudiness and temperature, and the occurrence of precipitation.

LI
Lines of instability may appear on air masses, especially mTA and mTC, identified in the synoptic chart as a trough line.These lines of pressure intensify the convective movement.Rainfall and increase of the humidity can usually be observed.
After identifying all the atmospheric systems and weather types at a daily and hourly level, the next step is to analyze the data generated through a qualitative point of view.It is important to observe the frequency of the atmospheric systems over the months and seasons for the two standard-years (a wet and a dry one), and to compare them.However, the technique of rhythmic analysis consists of the analysis of the sequence of weather types, making it possible to identify the most frequent atmospheric systems, standard, or distinguished atmospheric systems in a continuous sequence.To observe the rhythm, which is more or less a regular return to the same atmospheric conditions and weather patterns, year after year, we analyzed and identify the chain of atmospheric sequences that normally occur in the region, in which the resident population is accustomed to the resulting weather types, and identified the atmospheric sequences that cause extreme episodes in the region.

General Analysis of Rainfall Variability
The monthly variability of the rainfall volume at the Caconde station was high, as can be seen in Figure 4.A preliminary analysis of the figures enables us to decipher the measures that are present in the graph, showing the monthly volume of rain that usually occurs in the region, between the first (25th percentile) and the third quartile (75th percentile).
The seasonality of the rains observed in Figure 3a was expected, representing a regime of high precipitation in spring and summer in contrast to autumn and winter, according to the meteorological calendar of the four seasons in the southern hemisphere.However, when we observe the year 1983 in Figure 3b, only August is found to be below the historical average, and in nine months the rainfall volume was higher than the third quartile, and in some cases could be classified as an outlier, demonstrating the exceptional nature of this period.In the year 1984, seven months appear below the average, four of which are below the first quartile, but three of them (January, February, and March), are registered as dry periods in usually wet months, which gives evidence of a possible meteorological drought [68], during the months of low precipitation until the month of December.To observe the time variable of the analyzed data in this general analysis, Figure 5 relates the sequential monthly data of precipitation to the water flow in the Pardo River.As observed in the previous figure, the seasonal rainfall variability was highlighted, evidencing the region's rainfall regime.When associated to the natural flow of the Pardo River, there is a clear correspondence between the increase in the volume of the rainy months and the growth in the water flow, an ordinary situation observed in small and medium-sized hydrographic basins.One point of the figure that stands out is the recorded flood peaks, especially those exceeding 400 m 3 /s, values that are four times higher than the monthly average for the month of the greatest average flow (January-98.8m 3 /s), and which can cause large impacts along the drainage of the region.In total, this number is exceeded six times in the historical series, and although the maximum peak was recorded in the last days of the year 1966, the year of 1983 presents two of these flood peaks, in the months of January and December.

Spatial Analysis of Rainfall
In the previous procedure, the monthly precipitation volumes of the Caconde station gave an overview of the Alto Pardo basin, but described the variability from one specific rain gauge.To obtain a spatial view of the rain in the basin, Figure 6 distributes the monthly rain with contour lines, with 10 millimeters equidistant along the study area, through the kriging method for the years 1983-1984.It is important to emphasize that the precipitation maps were made to contribute to the   As observed in the previous figure, the seasonal rainfall variability was highlighted, evidencing the region's rainfall regime.When associated to the natural flow of the Pardo River, there is a clear correspondence between the increase in the volume of the rainy months and the growth in the water flow, an ordinary situation observed in small and medium-sized hydrographic basins.One point of the figure that stands out is the recorded flood peaks, especially those exceeding 400 m 3 /s, values that are four times higher than the monthly average for the month of the greatest average flow (January-98.8m 3 /s), and which can cause large impacts along the drainage of the region.In total, this number is exceeded six times in the historical series, and although the maximum peak was recorded in the last days of the year 1966, the year of 1983 presents two of these flood peaks, in the months of January and December.

Spatial Analysis of Rainfall
In the previous procedure, the monthly precipitation volumes of the Caconde station gave an overview of the Alto Pardo basin, but described the variability from one specific rain gauge.To obtain a spatial view of the rain in the basin, Figure 6 distributes the monthly rain with contour lines, with 10 millimeters equidistant along the study area, through the kriging method for the years 1983-1984.It is important to emphasize that the precipitation maps were made to contribute to the As observed in the previous figure, the seasonal rainfall variability was highlighted, evidencing the region's rainfall regime.When associated to the natural flow of the Pardo River, there is a clear correspondence between the increase in the volume of the rainy months and the growth in the water flow, an ordinary situation observed in small and medium-sized hydrographic basins.One point of the figure that stands out is the recorded flood peaks, especially those exceeding 400 m 3 /s, values that are four times higher than the monthly average for the month of the greatest average flow (January-98.8m 3 /s), and which can cause large impacts along the drainage of the region.In total, this number is exceeded six times in the historical series, and although the maximum peak was recorded in the last days of the year 1966, the year of 1983 presents two of these flood peaks, in the months of January and December.

Spatial Analysis of Rainfall
In the previous procedure, the monthly precipitation volumes of the Caconde station gave an overview of the Alto Pardo basin, but described the variability from one specific rain gauge.To obtain a spatial view of the rain in the basin, Figure 6 distributes the monthly rain with contour lines, with 10 millimeters equidistant along the study area, through the kriging method for the years 1983-1984.It is important to emphasize that the precipitation maps were made to contribute to the analysis of the monthly frequency of the atmospheric systems, since they facilitate the observation of the rainfall occurrence along several pluviometric stations located in the watershed.Even considering the spatial and temporal variability of the basin, it is possible to notice more rainfall in the high altitudes, especially in the wettest months of the period under examination.In the year 1983 it is possible to notice an extremely high volume of rainfall in the first three months, and there are cases that exceeded 400 mm monthly rainfall, and if you add up the first three months, the recorded volumes were above 1000 mm.Such volumes of rain observed contrast with the year of 1984, that presented little precipitation during the same period, such as in the case of the months of January, February, and March, which usually have more rain in the region.The normally dry months in Alto Pardo, between April and September, in the year 1983 registered a decrease in rainfall only between June and August, while in 1984 there was already a summer drought, which continued to show low precipitation by the end of November.At the end of the two years, the month of December was especially rainy for both periods.

Atmospheric Systems and Rhythmic Analysis
The research procedure for the identification of atmospheric systems was applied to the years 1983 and 1984, and this stage of the study intended to synthesize the results in a clear and objective way.Tables 4 and 5 indicate the frequency of the atmospheric systems for each of the months and the synthesis of the total frequency of the frontal systems and air masses with tropical and polar characteristics.In addition, with the aim of better summarizing the results, bar graphs were generated from the data described in the tables, which allow a better comparison between the two years.Figure 7 illustrates the frequency of the atmospheric systems at the monthly level for the year 1983 (a) and 1984 (b).The results observed in the tables indicate a high frequency of occurrence of frontal systems and other disturbing influences in the south of Brazil in the year 1983, as compared to 1984.These results reinforce the idea that the performance of the polar front and its unfolding are the main factors that are responsible for the occurrence of intense rainfall in the state of São Paulo, as [69,70] suggests, and the observation of which will be discussed later in the paper.Another point that draws attention is the significant occurrence of the stationary front in the year 1983, while in the year 1984 it appears only in the first month of the year.
When comparing the data obtained in the monthly contribution of the atmospheric systems and rainfall maps of the previous sub-section, there is a relationship between the dry months with the predominance of air masses, both of polar or tropical origin.However, the predominance of air masses of tropical origin in the year 1984 during the first three months was numerically far superior to the year 1983, and is a fact that stands out when comparing the rain and weather types in those two years.
Nevertheless, understanding the dynamics and chaining of the atmospheric systems, and more broadly their rhythm in the daily data, is important for understanding the origin of extreme episodes and what kind of weather types are usually common in a region.Thus, we analyzed the atmospheric systems in sequence in order to understand the climatic rhythm.The systems were related to the  The results observed in the tables indicate a high frequency of occurrence of frontal systems and other disturbing influences in the south of Brazil in the year 1983, as compared to 1984.These results reinforce the idea that the performance of the polar front and its unfolding are the main factors that are responsible for the occurrence of intense rainfall in the state of São Paulo, as [69,70] suggests, and the observation of which will be discussed later in the paper.Another point that draws attention is the significant occurrence of the stationary front in the year 1983, while in the year 1984 it appears only in the first month of the year.
When comparing the data obtained in the monthly contribution of the atmospheric systems and rainfall maps of the previous sub-section, there is a relationship between the dry months with the predominance of air masses, both of polar or tropical origin.However, the predominance of air masses of tropical origin in the year 1984 during the first three months was numerically far superior to the year 1983, and is a fact that stands out when comparing the rain and weather types in those two years.
Nevertheless, understanding the dynamics and chaining of the atmospheric systems, and more broadly their rhythm in the daily data, is important for understanding the origin of extreme episodes and what kind of weather types are usually common in a region.Thus, we analyzed the atmospheric systems in sequence in order to understand the climatic rhythm.The systems were related to the rainfall and natural water flow at a daily level, thus giving a view of the chain of events and the repercussions on rainfall and water dynamics.Figure 8a,b shows the results of the atmospheric system associated to the described variables.
The atmospheric systems sequences observed through the rhythmic analysis chart indicate a high variability of weather types during the year 1983, especially in the rainy periods.An example of these intense atmospheric dynamics is the period from 12 to 21 of January that begins with the advance of a cold front over the region, and found resistance to the advance toward the northern region of the study area.This resistance resulted in a warm front, which stopped and remained in the region for a period of two days.On day 16, when the system threatened to move, it received more moisture from another low pressure system that remained over the Alto Pardo watershed and dissipated only on day 21.This sequence of chained frontal systems resulted in an extreme summer precipitation event, with prolonged rainfall and was largely distributed on a regional scale.As a result, the water flow reflected the heavy rains recorded along the drainage area for nine days and, on 21 January, reached the peak of 484 m 3 /s natural river flow.The episode described in this paragraph is just one example among several possibilities of observation through the daily analysis of the atmospheric systems.
rainfall and natural water flow at a daily level, thus giving a view of the chain of events and the repercussions on rainfall and water dynamics.Figure 8a,b shows the results of the atmospheric system associated to the described variables.
The atmospheric systems sequences observed through the rhythmic analysis chart indicate a high variability of weather types during the year 1983, especially in the rainy periods.An example of these intense atmospheric dynamics is the period from 12 to 21 of January that begins with the advance of a cold front over the region, and found resistance to the advance toward the northern region of the study area.This resistance resulted in a warm front, which stopped and remained in the region for a period of two days.On day 16, when the system threatened to move, it received more moisture from another low pressure system that remained over the Alto Pardo watershed and dissipated only on day 21.This sequence of chained frontal systems resulted in an extreme summer precipitation event, with prolonged rainfall and was largely distributed on a regional scale.As a result, the water flow reflected the heavy rains recorded along the drainage area for nine days and, on 21 January, reached the peak of 484 m 3 /s natural river flow.The episode described in this paragraph is just one example among several possibilities of observation through the daily analysis of the atmospheric systems.The year 1983 recorded several rainy episodes and, in addition to that discussed in the previous paragraph, other similar episodes were recorded that were related to the stationary front, which occurred between 30 January and 5 February, 4-8 March, 28 May and 4 June, 17-21 July, 22-25 September, 11-15 and 21-24 October, and 10-15 and 22-28 December.The natural flow of the Pardo River reflected such rainy episodes, and had a rise in the data during and after each of these episodes.Among them was the day 16 December, which after the dissipation of the frontal system and the migration of the old polar mass towards the region, recorded a river flow peak of 493 m 3 /s.
The year 1984 was characterized by high-pressure typical conditions from the air masses, in particular the Atlantic tropical mass and mTAC, generating little rain in the study area, clear skies, and hot and dry weather types.The only period of occurrence of the stationary front recorded in the rhythmic analysis graph was the second half of January; nevertheless, the frontal system presented little intensity and did not register a very high volume of precipitation, even dominating the regional The year 1983 recorded several rainy episodes and, in addition to that discussed in the previous paragraph, other similar episodes were recorded that were related to the stationary front, which occurred between 30 January and 5 February, 4-8 March, 28 May and 4 June, 17-21 July, 22-25 September, 11-15 and 21-24 October, and 10-15 and 22-28 December.The natural flow of the Pardo River reflected such rainy episodes, and had a rise in the data during and after each of these episodes.Among them was the day 16 December, which after the dissipation of the frontal system and the migration of the old polar mass towards the region, recorded a river flow peak of 493 m 3 /s.
The year 1984 was characterized by high-pressure typical conditions from the air masses, in particular the Atlantic tropical mass and mTAC, generating little rain in the study area, clear skies, and hot and dry weather types.The only period of occurrence of the stationary front recorded in the rhythmic analysis graph was the second half of January; nevertheless, the frontal system presented little intensity and did not register a very high volume of precipitation, even dominating the regional atmosphere for several days.The predominance of tropical air masses, few and low intensity frontal systems, in addition to low amounts of cold waves, was something that was significantly observed during the year through rhythmic analysis.Due to this reason, the natural water flow reduced gradually throughout the year, and did not show any peak flow, prolonging the drought and the scarcity of water resources in the region.It was only in the second half of November that it was possible to observe higher volumes of rain, brought about by the advance of frontal systems and lines of instabilities that formed the rainy month of December after a long period of drought.

Discussion and Conclusions
This study examined the development of the 1983-1984 atmospheric systems and analysed the rhythm between these years.Through the methodological procedures adopted, it was identified that the frontal systems were mainly responsible for the origin of the rainfall episodes at the monthly and annual levels.The combination of these highly changeable weather conditions over the course of several days resulted in 88.5% of the total volume of rainfall in 1983 that precipitated in days under the passage of frontal systems, and the hasty volume in the year 1984 during the passage of the frontal systems was 76% of the annual total.These passages of cold fronts by the incursions of the Atlantic polar mass, whether through paths of easternmost longitude and damp or through more continental and dry routes, were more frequent in the year 1983, in some cases for several sequential days, and less frequent in the year 1984, with the predominance of the air mass systems.Another point was that, due to the greater number of polar advances in the year 1983, there were more days under the dominion of the polar air masses, mainly in the winter, with a relationship between the climatic rhythm with the number of days of cold weather.
The frequency of atmospheric systems during the months, as listed in the tables, is important because it is possible to interpret the weather types that have been most frequently in the Alto Pardo region during this period.Such data, when related to the spatial distribution of rainfall, allows an interpretation of this frequency at the regional level, especially for explaining the wettest months, which increased the frequency of frontal systems that resulted in months of higher occurrence of rains upstream of the Pardo River.However, the analysis cannot be limited only to the monthly level of the frequency of atmospheric systems, because the weather types interfere at the daily level, and in some cases at the hourly level, in the natural events and human beings.
These procedures approach the climatology of the geographical analysis, by enabling the correlation of this string of weather types to the geographical variables.For example, the month of February 1983 registered high rainfall in the north of the basin, although the south of the Alto Pardo watershed has a high altitude.With the daily analysis, it is possible to check the alternating situation between the warm front and the tropical continental mass with lines of instability during the days 6 to 13, which probably contributed to the high precipitation in the northwest basin.Another fact that the monthly analysis does not show is the predominance of the tropical masses after day 14, resulting in clear skies and hot and dry weather types for almost half of the month, even in the summer of a wet year.
The frequency of frontal systems in the study area was observed mainly in the year 1983, predominating in 34.1% of days, but even in the year 1984 the fronts were frequent, more specifically in 22.4% of days.The constant currents that the South displaced to the region of the Alto Pardo watershed, and that can interfere with the weather types, occurs due to the basin being located near the Tropic of Capricorn, more specifically to the north of the line, characterized for being the threshold between the tropical and temperate zones in the southern hemisphere of the planet.In a classical climate regionalization proposed by [71] and confirmed by [72,73], it is estimated that there is a large-scale limit of the climate transition near the tropic, however this does not have an accurate and static nature [26].That is, there is a great variability in the frequency of the frontal systems that move to the southeastern region of Brazil, and may have an impact on the source of regional rainfall.
A factor that must be considered for the planning of water resources in the region is the variation of the El Niño-Southern Oscillation (ENSO) phenomenon.Although there is no clear correlation between the rainy years and the occurrence of El Niño in Southeast Brazil, as in the southern region of Brazil [74], a very strong intensity of El Niño may have influenced the year 1983.In the analyzed region, the occurrence of the Stationary Polar Front (EST) was associated with most of the cases of intense rainfall on a timescale of two years.In addition, in some cases, this convective system is fed by moisture from the northern region of Brazil, configuring the South Atlantic Convergence Zone [75].Its occurrence is mainly in the spring and summer, and can cause intense rainfall, and has proven to be influential for the occurrence of floods in the region.
In the case of the meteorological droughts, high-pressure atmospheric systems such as the mTA and mTAC predominated during a period normally rainy as is the case in the summer, and the chaining of these weather types, in combination with mTC and the frontal systems in dissipating toward the ocean, contributed to the low rainfall observed on the spatial scale in the first months of the year 1984, that resulted in the reduction of the levels of natural water flow throughout the year, whereas autumn and winter feature a pattern that is usually dry in the region.A similar case in southeastern Brazil is described by [76] for the summer of 2014, that proposed the existence of a blockade for the passage of frontal systems due to the anomalous high pressure over the southeast region of Brazil.
The comparison between two extreme years, a wet and a dry one, with different atmospheric circulation patterns during the year, allowed us to obtain relevant results for the study area.However, the rhythmic analysis of different weather types, and their succession over time [77], cannot be limited only to these two years of study if we aim to achieve results even more relevant and applicable.For future research, we intend to apply the method of rhythm analysis to a historical series that includes at least 15 years of data, to gradually identify the succession pattern of weather types in the Alto Pardo watershed and Southeast Brazil.The long series of data can go through qualitative analysis, and identify the chaining of the atmospheric systems at the source of rain and drought episodes, beyond the statistical treatment to try to identify the relationship between the weather type and precipitation.We hope to apply the results of the rhythmic analysis to geographical variables in the region, as in the case of floods, air quality, agriculture, health, extreme events, and several others.Finally, we want the results of the manual method of identification of the atmospheric systems to contribute in the future to build an automatic method for South America, to identify the air masses and frontal systems through the synoptic data (charts) and meteorological data.

Figure 2 .
Figure 2. Elevation and pluviometry stations in the Alto Pardo watershed.

Figure 2 .
Figure 2. Elevation and pluviometry stations in the Alto Pardo watershed.

Figure 3 .
Figure 3. (a) Air masses by the source region and their main incursions into the Alto Pardo according to the data of [28,67] (b) Example of a synoptic chart provided by [57] (20 May 1983).

Figure 3 .
Figure 3. (a) Air masses by the source region and their main incursions into the Alto Pardo according to the data of [28,67] (b) Example of a synoptic chart provided by [57] (20 May 1983).

Figure 6 .Figure 6 .
Figure 6.Distribution of monthly rainfall in the Alto Pardo watershed in (a) 1983 year; (b) 1984 year.310 Figure 6.Distribution of monthly rainfall in the Alto Pardo watershed in the year (a) 1983; (b) 1984.

Figure 7 .
Figure 7. Atmospheric systems frequency at the monthly level in the tear (a) 1983; (b) 1984.

Figure 7 .
Figure 7. Atmospheric systems frequency at the monthly level in the tear (a) 1983; (b) 1984.

Figure 8 .
Figure 8. Rhythm of the atmospheric systems associated with rainfall and natural water flow in the year (a) 1983; (b) 1984.

Table 1 .
List of the 18 pluviometric stations used in the spatial interpolation.
1 Source: National Water Agency.

Table 1 .
List of the 18 pluviometric stations used in the spatial interpolation.

Table 2 .
Air masses: origin and criteria for classification.
mTASouth Atlantic Anticyclone high-pressure system over tropical latitudes, presenting high temperature and humidity characteristics.

Table 3 .
Frontal Systems and Individualized Systems: criteria for classification.

Table 4 .
Atmospheric systems frequency at the monthly level in the year 1983.

Table 5 .
Atmospheric systems frequency at the monthly level in the year 1984.

Table 5 .
Atmospheric systems frequency at the monthly level in the year 1984.