Trends in Extreme Indices and Seasonal Analysis of Precipitation and Temperature in the Northwest Region of Rio Grande do Sul, Brazil ()
1. Introduction
Probably the most important environmental challenge of this century is to adapt to climate change and develop strategies to minimize its effects. The major climate changes predicted for this century include an increase in average temperature, elevated CO2 concentration in the atmosphere, changes in rainfall spatial distribution, increase in evapotranspiration during summer and greater occurrence of extreme events such as windstorms, hail and heavy rainfall [1] .
There is already a consensus that, even if all precautions to reduce emissions and carbon recovery are taken, an increase of 2˚C in the average temperature of the Earth till the end of the XXI century (in a more optimistic scenario) seems inevitable [1] [2] .
As a consequence, the global average surface temperatures registered during the last decade showed the highest increases, probably through the last millennium, being more expressive in the northern high latitudes. On the other hand, it is more difficult to characterize changes in precipitation regimes, as this is strongly influenced by regional climate phenomena such as El Niño and La Niña in South America. About this fact, the El Niño of 1998 was one of the most severe occurrences of the last century [3] .
The meteorological variables exhibit both spatial and temporal variability and, allied to this fact, it is known that the climate models are still very inaccurate. It is possible that, even using different models to evaluate the uncertainties on climate projections, future climate conditions may not be covered by any of them.
These projections are generated by general or regional circulation models (GCMs and RCMs, respectively) as a result of future scenarios of climate forcings caused by the emission of greenhouse gases and aerosols [4] . Such models comprise partial differential equations systems based on physical laws of fluid motion and chemical substances [5] , ranging from simple approaches of punctual energy balance to three-dimensional general circulation models, which attempt to model all complexities of the climate system of the Earth [3] . Since it is extremely difficult to properly simulate the climate, climate models often contain some simplifications, repre- sented by different parameterizations. The use of different sets of parameters results in different output simulations or members, representing how sensitive to the parameters the model is. Thus, it is usual to considerer a “control” member associated with other disturbances such as a “high” member refers. The term “high” refers to a large disturbance on outputs associated to a certain set o parameters, whereas a “low” member indicates that the model is not that much sensitive, for example.
One can cite numerous papers on climate change which have used such models globally [6] [7] ; some continental, for South America [8] - [10] ; and for Brazil and/or some of its regions [11] - [18] . There is a consensus in most of the studies on climate change performed for South America that there will be an increase (decrease) in the frequency of warm (cold) nights and a higher frequency of days with heavy rainfall [14] .
The study conducted for Brazil and presented by [19] included six GCMs and considered only the most pessimistic (A2) and the most optimistic (B2) scenarios. Climate scenarios are representations of future climate based on assumptions about changes in the composition of the atmosphere, as consequence of human activities. The projections for future climate in the country were compared in terms of anomalies with respect to the baseline scenario represented by the period of 1961-1990. The trends indicate an increase in temperature of up to 4˚C in 2080, while larger differences are likely to occur from June to October, i.e., winters will become warmer compared to the summers. On the other hand, for changes in average daily precipitation in Brazil, there will be a decrease in the monthly amounts compared to the baseline period, but only few differences between future periods of short (2025s), medium (2055s) and long-term (2085s) are projected.
Their results also reveal that the south region of Brazil is likely to suffer with the largest increases in temperature especially during the spring months, while the projection for 2080 according to the A2 scenario is 5.5˚C above the baseline average. Furthermore, it is estimated that the highest rainfalls for Brazil occur during the summer (40%) and autumn (29%) months. The precipitation anomalies observed in the South region showed a decrease throughout all the seasons.
It should be noted that the forecasts made for Brazil by [19] were obtained by the averaged projections of six GCMs and contradictory results on precipitation (increase rather than decrease) were found by other authors. For instance, [20] predicted an increase in annual precipitation for the entire south region of Brazil, while [21] found no significant changes or only a small increase up to the end of this century. These results serve to demonstrate that, among other aspects, the uncertainties of the models remain large and difficult to be assessed, justifying more researches, especially at regional scale with the use of RCMs.
It is known that 80% of the cultivated areas in the world are non-irrigated [20] . Thus, agricultural productivity becomes one of the most vulnerable sectors in developing countries, such as Brazil. In this context, Brazil is responsible for 27% of global soybean production [22] , and the State of Rio Grande do Sul (RS) is the third largest producer in this country. Yet, the north and northwest regions are the largest producers of soybean in this state, becoming evident that the biggest challenge faced by agriculture over this century will be to deal with possible climate changes.
The main agronomic effects of climate change are related to the respiratory cycle of the plant, effects of changes in temperature and precipitation and rising sea level with a consequent reduction of cultivated areas due to an increase in soil salinity [23] . Another expected consequence of this century is that the elevated temperatures will cause a delay in the sub-periods of crop sowing and flowering, as evidenced by some authors for various crops in future scenarios [23] - [25] .
In general, with respect to precipitation, it is expected that the rainy season will become wetter, while dry seasons will be even drier. In locations where there is an increase in water deficit for agriculture during the dry months, irrigation will require larger quantities of water, intensifying conflicts over its use [26] . It is known that continuous periods of time without precipitation lasting for more than 10 days during the rainy seasons negatively affect agricultural productivity [22] .
It is possible that the soybean cultivation in south Brazil becomes increasingly difficult. The occurrence of heavy precipitation and strong winds may limit the use of machinery in some agricultural areas [27] , besides the fact that hailstorms can destroy crops almost instantly. It is therefore extremely important to determine the frequency and intensity of extreme events, especially on precipitation, and which are their preferred locations of occurrence.
As mentioned by [8] , the development of climate change scenarios at a regional scale can reveal pronounced variability, which may suggest a geographic distribution of impacts. In their studies, an increase of approximately 3˚C - 4˚C for the south region of Brazil is expected, and for this region, most of the models of the Intergovernmental Panel on Climate Change (IPCC) fourth Special Report agree to an increase in precipitation under a warmer climate.
Thus, given the agronomic importance of the State of Rio Grande do Sul in the Brazilian economy and the lack of studies on climate change specifically for the most productive region of this State, the present study aims to conduct an investigation to detect changes in temperature and precipitation in the northwest region of this state with the use of different GCMs and RCMs climate models projections.
2. Materials and Methods
2.1. Study Area
The study area comprises the northwest region of the state of Rio Grande do Sul (Figure 1). The geographical location of this state inserts it under the influence of atmospheric systems such as cold fronts and mesoscale convective systems, which may cause natural disaster and severe extreme events. In addition, the climatic phenomena El Niño and La Niña also determine the occurrence of floods and droughts in southern Brazil [28] .
A recent study presented by [28] and representative of the 2003-2009 period showed that 70% of drought events recorded in the Rio Grande do Sul correspond to the northern half of the State. Moreover, higher frequency of dry events is registered during the summer months, so that crop yield during this period decreases because the crop cultivation period coincides with the most common dry months, from November to March [29] .
The land use in the northwest region is essentially agricultural and, in most of the cases, crops are non- irrigated. In addition, this region is the largest soybean producer in the state, strongly contributing to its economic development. These reasons justify the choice of this region to analyze the vulnerability of this area under a changing climate.
2.2. Meteorological Data
Points in Figure 1 indicate the seven locations for which there are climate projections of the meteorological variables in the region. These projections include precipitation and temperature for the baseline period (1961- 1990), near term (2011-2040), midterm (2041-2070) and far term (2071-2100). Future periods will be repre- sented in this work by the central year on each (2025s, 2055s and 2085s). The information for each location is shown in Table 1.
Ten different projections for the A1B greenhouse gas emission scenario were considered at each location: five GCMs and five projections derived from the ETA regional climate model. The A1B scenario belongs to the
Table 1. Meteorological data locations in the northwest region of Rio Grande do Sul.
family of climate scenarios A1, which suggests a future world of very rapid economic growth, global population that peaks mid-century and declines thereafter, and rapid introduction of new and more efficient technologies. The A1B scenario refers to a balance across fossil-intensive e nonfossil energy sources [30] .
Global models description is shown in Table 2. The ETA model is descendant of the HIBU model (Hydrometeorological Institute and Belgrade University), previously developed by [31] . It is a regional climate model that couples the HadCM3 general circulation model as lateral boundary condition for mesoscale simulations. A more detailed description of this method can be found in [32] . The coupling of this global model considers a set of regional members (CTRL, LOW, MID and HIGH), which represent the climate sensitivity. The horizontal resolutions of the ETA model applied in this study were 20 and 40 km, then providing five projections: ETA 20-CTRL, ETA 40-CTRL, ETA 40-LOW, ETA 40-MID, ETA-40 HIGH. Thus, for each of the 7 locations shown in Table 1, there are 10 sets of meteorological data derived from these different models. Additionally, only the member CTRL of all GCMs was applied.
Table 2. General circulation models applied in this study.
*Abbreviation adopt in this study. Adapted from [33] .
Daily meteorological data refer to the mean temperature (T, ˚C) and precipitation (P, mm). As agriculture in the northwest region is moved primarily by soybean cultivation, it is necessary to determine whether changes in precipitation and temperature may occur during the months in which soybeans are cultivated. Therefore, we first carried out a monthly analysis of T and P data in order to investigate this point. In a second moment, we analyzed the occurrence of extreme rainfall events by the use of climate extreme indices (Table 3). This analysis could not be performed for the temperature, since only mean temperature data were available for the study area.
In Table 3, a wet day occurs when the daily rainfall rate is greater or equal to 1 mm, and a dry day is when this rate is lower than 1 mm. The same climate extreme indices (or some of them) were recently applied by other authors such as [13] and [14] . These indices will be of great importance, since they help detect the occurrence of extreme events related to precipitation. We are also interested in analyzing spatial distribution of these extreme events, considering seven different locations in the northwest region.
In order to detect trends in annual indices, the Mann-Kendall test [34] [35] was performed. It consists of a non-parametric test for assessing if there is a monotonic upward or downward trend of the variable of interest over time. Thus, it can be seen as a nonparametric test for zero slope of the linear regression of time-ordered data versus time [36] . The test was conducted at a 5% significance level.
3. Results and Discussions
3.1. Annual Analysis of Temperature and Precipitation
A brief analysis of annual meteorological data is presented in Table 4 and Table 5. Table 4 provides the average annual temperature and the anomalies for all different locations in the study area. The anomaly refers to the difference between a future time slice and the baseline period. It is clear that the projections indicate that the temperature will increase in all locations, when compared to the baseline period (1961-1990). The largest anomalies were predicted for the end of the century, and a temperature diminution in the short term (2025s) is expected only in location 6. Since some locations are close to each other, similar projections were found probably because they belong to neighbor resolution cells of the climate model. They were maintained in the study because different precipitation projections were provided by the models even for these locations.
The smallest anomalies were found for locations 4, 5 and 6, situated in the north, central and south portion of the region, respectively. Analyzing the long term period, we can note mean annual temperatures expected to be from 20.6˚C to 22.7˚C, which may suggest a spatial distribution of the changes of temperature and that there may be an increase of almost 3˚C, in average, till the end of this century. Similar result were found by [8] for the Southern Brazil for A2 scenario. They found that the warming is projected by about 3˚C - 4˚C given by regional climate models.
Analyzing the information presented in Table 5, one may observe that, as well as for temperature, it is expected an increase in the annual amount of rainfall, and it is more expressive for locations 2, 3, 4 and 5. Again, only in location 6, the models projected a reduction in precipitation in the near term. Indeed, for this location, the smallest anomalies were found, suggesting that probably it will not be affected by climate changes, or only small climate disturbances (not significant) might occur. The worst scenarios were projected for the same locations abovementioned, situated in higher altitudes. Increase in precipitation for Southern Brazil was also found by [8] . This is probably due to an increase in dry days, whereas the consecutive wet days are unchangeable (see Section 3.3).
Table 3. Climate extreme indices related to precipitation of climate projections for the northwest region [14] .
Table 4. Mean annual temperature and anomalies for each location in the northwest region of RS. The values correspond to the average of all ten projections (five global and five regional projections).
Table 5. Mean annual precipitation and anomalies for each location in the northwest region of RS. The values correspond to the average of all ten projections (five global and five regional projections).
Based on the data presented in these tables, we may conclude, a priori, that it could be expected positive impacts of these future scenarios. However, these values refer to the averaged annual temperature and precipitation for all models and they show nothing about the seasonal cycles of these variables. It is possible that higher amounts of precipitation are concentrated in months which are not interesting for agriculture in the northwest region of Rio Grande do Sul. Thus, to reasonably identify climate changes that might have impacts on agriculture, it is crucial to understand how these changes are distributed throughout the year. It is also important to point out that, in this study, we are not interested to evaluate the impacts on agriculture themselves, but only to detect whether they are expected to occur or not. How intense or frequent they are probably to be is out of the scope of this study. So, the analysis presented here is an attempt to identify climate changes in the region, based only on the temperature and precipitations projections.
3.2. Seasonal Analysis of Temperature and Precipitation
An initial analysis of temperature and precipitation projections was performed by taking the average of the monthly values of these variables given by all models in all locations and then calculating the anomalies compared to the baseline period (1961-1990) (Figure 2). The anomalies refer to the differences between each future period of 30 years and the baseline period, also represented by a series of 30 years of data.
It should be noted that the largest increases of temperature are projected for winter and early spring, from June to October. Additionally, the seasonal variability is similar for all future time slices. The average annual anomalies for 2025s, 2055s and 2085s are 1.2˚C, 2.1˚C and 2.9˚C, respectively. The months with higher temperature anomalies are May and October (2025s), and August and September (2055s and 2085s). These months do not coincide with the summer months of the main crop cultivation (soybean) in the region, which is practiced from November to March. Indeed, these are the months when minor anomalies were obtained. However, it is projected that such anomalies are in the range of 1.0˚C to 2.9˚C, which may not be neglected, especially at the end of the century.
The lowest mean anomalies provided for precipitation are expected during the winter months. For August, a reduction in precipitation of nearly 10 mm for the periods of 2025s and 2055s is expected, but any change for 2085s. For the remaining months of the year, all the anomalies are positive, i.e., increase in precipitation. This increase is more significant from April to May and from October to December. These results are coherent with those presented by [8] , who also found an increase in precipitation for fall months and a reduction during the winter at the end of the century. The largest increases in precipitation coincide with winter crops (wheat and oats, sown in May) and the sowing period of the summer crop (soybean, sown in November). Thus, one may conclude that it is expected that the impacts will be positive with respect to precipitation. Furthermore, this increase in precipitation during the months of the summer crop also coincides to minor increases in temperature during this period.
However, higher amounts of precipitation does not necessarily indicate a positive impact, because higher amounts of rain concentrated in short periods of time can be associated to a greater occurrence of extreme events such as heavy rains and floods. If this is the case, then it is possible that the period of sowing of the summer crop (more important economically) is more strongly affected by climate change (discussed in the next section).
In order to evaluate the differences given by each climate model for each location, Table 6 presents the results of monthly temperature projections, divided into: averaged for general circulation models (AVG) and averaged for regional climate model projections (AVR). That is, each value in this table is the average of five projections of the ETA RCM or five projections of the GCMs. In the case of precipitation, it was decided to make a representation by graphs (Figures 3(a)-(g)), since the extent of the differences are larger and more noticeable than for temperature.
Figure 2. Seasonal cycles of temperature and annual precipitation anomalies. The curves correspond to the variable averaged for all climate projections and locations.
Table 6. Average monthly temperature projected by different climate models in all locations.
aAVG: monthly mean temperature averaged for GCMs projections; bAVR: monthly mean temperature averaged for RCMs projections.
In most cases, the results demonstrate that higher temperatures are projected by ETA model with its various projections. Only for April (all locations), higher temperatures were provided by the projections of GCMs. Moreover, higher temperatures are expected in locations 1, 2, 6 and 7 in the east and central portions of the region. December, January and February are the warmest months, while June and July are the coldest months, with amplitude of more than 10 degrees among these months. All projections indicate an increase in temperature in all locations, and it is possible that winters will become warmer compared to the summers. The anomalies in 2025s have lower fluctuations than for subsequent periods and it is expected that the major differences will be felt between 2025s and 2055s than between 2055s and 2085s.
The study conducted by [37] showed similar results on the analysis of average and minimum temperatures in southern Brazil during the 1960-2002 period. They are consistent with the trends provided by the projections used in this study. The authors point out to a warming in the southern region of the country with positive trends in maximum and minimum temperatures. They add also that warming will be more intense during the winter than in the summer, probably due to the greater number of warm days in the winter.
Regarding precipitation, GCMs agree with each other in different locations (lower variability). In some cases, however, they provide projections contrary to the ETA model. They also provide lower seasonal variability with respect to the baseline period. About this fact, it is observed that for the months of December, January and February, the GCMs indicate that there will be little or no change in comparison to 1961-1990. Hence, this entails that it is not expected changes during most of the months of the summer crop (November to March). Additionally, smaller increases in temperature for these months may indicate that the impacts are likely to be less severe. Regarding precipitation, GCMs agree with each other in different locations (lower variability). In some cases, however, they provide projections contrary to the ETA model. They also provide lower seasonal variability with respect to the baseline period. About this fact, it is observed that for the months of December, January and February, the GCMs indicate that there will be little or no change in comparison to 1961-1990. Hence, this entails that it is not expected changes during most of the months of the summer crop (November to March). Additionally, smaller increases in temperature for these months may indicate that the impacts are likely to be less severe.
Only the months of March, May and November had strictly positive anomalies for all future periods in all locations. The months of July, August and September had the largest reductions, and almost all anomalies were negative in August, generated by both GCMs and the projections of the ETA model. In these months, negative anomalies projected by ETA model were obtained only for the months of February and December in locations 1 and 6, respectively, in all future periods. The projections for August provided the highest reductions, reaching anomalies of up to 20 mm.
Comparing the anomalies between AVG and AVR, it is observed that all anomalies calculated for the ETA model are much larger than those provided by GCMs. Hence, the regional climate model projects a situation that is more favorable associated to an increase in precipitation amounts. This increase, in general, is greater the further it gets from the baseline period for each month, except for the months of February, July and September, when we observe the largest increases in the middle of the century (2055s). The major positive anomalies were obtained for May and November, which are the months when the winter and summer crop are sown, respectively.
In general, it can be said that the models agree with each other that during the months of summer crop there will be an increase in precipitation, and it is possible that agriculture will not be negatively affected by this variable. With regard to the anomalies, it was found that the projections of positive anomalies have greater magnitude than the negative anomalies.
3.3. Extreme Rainfall Indices
To facilitate the visualization of the results, only the mean values of the selected indices for rainfall projections given by GCMs and the RCM projections (AVG and AVR, respectively) are presented in Table 7, for each location. However, the Mann-Kendall test was performed for each annual series of the indices considered in each location and for each model.
The first index, R95p, refers to the total rainfall of the days when precipitation is greater than the 95th percentile. Thus, if there is an increase trend of this index, this may indicate that the increase of precipitation discussed in the previous section is partly due to the occurrence of heavy rain (more concentrated in time). In all locations, this index was higher for GCMs than for the ETA model and is greater the further it gets from the baseline period (Table 7). However, the hypothesis test reveals increasing trends in all locations only for the ETA 20 model
Table 7. Mean annual precipitation indices for all locations and time slices, based on averaged regional and global climate models.
and any trend for HADCM3 and NRCCCSM models in any location. It was found that the locations 6 and 7 showed no trends by none of the models. This fact indicates a spatial distribution of impacts, suggesting that the climate near these locations is different from the rest of the region and will be less affected by a changing climate. Statistical evidence of trend was found for almost all regional projections, while for the global model only GFCM21 presented trend in all locations, except in location 7.
Although average precipitation given by GCMs was higher, this result is not in contradiction with the hypothesis test, which demonstrated higher trends for the ETA model. This can be better understood in Figure 4. This figure presents the trends that have been tested for ETA 20 and NRCCCSM models. Note that, although the NRCCCSM model has a higher average R95p index, its trend line is less steep than for ETA 20 model. Thus, the difference in average values refers only to variability among models, while the Mann-Kendall test identifies trends over time.
All locations showed increasing trends for at least 5 of the 10 climate models, except for the location 7. In this location, trend was detected only for the ETA 20 and ETA 40 HIGH models for the level of significance adopted in this study. The graphs also show the large interannual variability of this index.
The analysis of R99p index follows the same logic as above, but considers the sum of the days with even more intense rainfall, above the 99th percentile. However, trend for this index was detected for only few projections, and in any location that has happened for more than 3 climate models. The R95p and R99p indices also indicate how much of the total annual precipitation is probably to occur in more intense or extreme events. For example, the total rainfall in 2050 in location 1 for the ETA 40 CTRL model is 2596 mm. Of this total, is expected that 532 mm (approximately 20%) precipitate in rainy days with precipitation above 95th percentile and 213 mm (approximately 8%) occur in days with precipitation above the 99th percentile. In agriculture, these may be the worst precipitation events, because the crops can be quickly destroyed. For models that registered positive trends may be due to either the occurrence of heavy rains or more frequent occurrence of these extreme rainfall events. To answer this question, the trend of the number of days necessary to generate a certain amount of precipitation should have been analyzed, but this is out of the scope of this work.
For the PRCPTOT index, which corresponds to the total annual precipitation of wet days (P ≥ 1 mm), the highest values were provided by the ETA model and its members. Except for locations 1 and 6, the ETA model predicted up to 600 mm more than the baseline period for this index in 2085s. With respect to the GCMs, they suggested an increase in precipitation of 160 mm in the whole region for the same period. The Mann-Kendall test for this index detected an upward trend only for ETA model and rejects this assumption for all GCMs.
Similar results were obtained for R30 index (days), which counts the number of days in the year with precipitation greater than 30 mm. This index also helps in the detection of extreme events, since they can be associated with the previous two indices. For example, the 95th percentile of ETA 20 in 2022 and 2033 (not shown) were 30 mm, approximately. If we associate the indices R95p and R30, we may notice that in 2022, only 5 days (R30, in 2022) were responsible for a rain of 260 mm (R95p, 2022); and in 2033, 7 days (R30, in 2033) accounted for a rain of 360 mm (R95p, in 2033). Although this analysis is not presented here for all projections and locations, it is very useful to understand how many days are needed to generate a certain amount of precipitation.
For GCMs, it is noted from Table 7 that the R30 index has no changes with respect to the period 1961-1990, and it is at most 3 days more than in 2085s. However, the rates calculated for the regional projections reveal up to 9 days with precipitation greater than 30 mm in 2085s. In this case, only the ETA 40 HIGH showed trend in locations 6 and 7 and any trend in location 5. In other locations, no more than 3 trends were detected for the other climate projections.
Figure 4. Differences among models and R95 index increasing trends (mm) for ETA 20 and NRCCCSM models, location 1 (y1 and y2 are trend lines for ETA 20 and NRCCCSM models, respectively).
For all other indices, any trend was detected neither for ETA model, nor for the GCMs. Comparing the RX1 and RX5 indices, one may observe that the maximum precipitation in a day (RX1) is proportionately much more extreme than the total rainfall in five consecutive wet days (RX5). Although there is no upward trend, these average values are of concern because they represent days with very intense rainfall. Finally, CDD and CWD indices varied between 16 and 21 days and between 7 and 10 days, respectively. There is no statistical evidence of trend for these indices. Locations 1 and 4 provided lower consecutive dry days and higher consecutive wet days.
4. Conclusions
This study aimed to detect evidence of climate change in the northwest region of the State of Rio Grande do Sul, Brazil. This region is the third largest producer of soybean in the state, contributing significantly to the national productivity. Ten projections of general and regional circulation models were considered in seven different locations in the region. The analysis was carried out for the period from 2011 to 2100 and compared with the baseline period (1960-1990).
In most cases, both the ETA model and GCMs agreed about an increase in temperature and precipitation in all seven locations. Projections indicate an increase of mean annual temperature of almost 3˚C till the end of the century, a little lower than found by other authors. Additionally, we expect that the locations 6 and 7 are likely to suffer with the highest average temperatures. Regarding annual cumulative precipitation, only for the location 6 is predicted a reduction in rainfall amounts. It can be concluded that it is possible that the climate in these locations (more severe towards the center of the state) is different from the rest of the region and that the spatial distribution of the impacts of climate change on agriculture will not be uniform.
The seasonal cycles of monthly temperature and precipitation show that the months of the summer crop will be the least affected by the increase in these variables. However, this fact does not mean that the changes in T and P are not significant. Regarding this fact, one may observe a great variability between the projections of ETA model and its members and the GCMs applied in this study. This can be confirmed by analyzing the trends in extreme precipitation indices, which showed that most of upward trends were detected for the indices provided by the regional model projections. Since RCMs are constructed for smaller areas with a higher resolution, their projections are believed to be closer to reality than global models. Thus, based on the abovementioned, it could be assumed that it is possible that future climate will be characterized by more frequent extreme events of precipitation in the northwest region.
Additionally, the results presented here are useful to demonstrate the influence of the downscaling on the results. Furthermore, the average precipitations projected by the GCMs were higher than for the ETA model, but increasing trends were not detected for most indices calculated from global projections.
Furthermore, the differences between the results found in this study and the others mentioned in this paper are in part due to different greenhouse gas emission scenarios considered. Therefore, our results are based on projections for a scenario (A1B) which is not as pessimistic as A2. For example, the work presented by [19] indicates temperature increases up to 5.5˚C in 2085s for A2 scenario.
Since there is no other research applying regional climate models for the northwest region of Rio Grande do Sul, we believe that this study contributed to improving the findings on climate change for the south of Brazil. Even if indirectly, this work also helps demonstrate that, for the same model, there is a great variability. This variability increases the uncertainty about the conclusions and exposes the difficulty in understanding the extent to which climate change may manifest in a particular location.
Acknowledgements
We are very thankful to the National Council of Scientific and Technological Development (CNPq) for sponsoring this research and acknowledge all anonymous reviewers for their suggestions and improvements in this paper.