Observed and Projected Changes in Temperature and Precipitation in the Core Crop Region of the Humid Pampa, Argentina

The core crop region of the Humid Pampa is one of the most productive agricultural lands around the world and depends highly on climate conditions. This study assesses climate variability, climate extremes, and observed and projected climate changes there, using 1911–2019 observations and CMIP5 model simulations. Since 1970, the annual mean temperature has risen by 1 °C and the mean annual minimum and maximum temperatures by 2 and 0.5 °C, respectively. The frequency of warm days and nights increased, and cold days and nights decreased. Heatwaves became longer and more intense, and cold waves decreased with less frost events. Annual precipitation increased by 10% from 1911, mainly in summer, and years with excess precipitation outnumbered those with a deficit. Both intense precipitation events and consecutive dry days grew, suggesting more annual precipitation falling on fewer days. Projections show a warming of 1 °C by 2035, regardless of the scenario. From then on until 2100, mean temperature will increase by 2 and 3–3.5 °C in the RCP4.5 and RCP8.5 scenarios, respectively. Annual precipitation will grow 8 and 16% from current values by 2100 in the RCP4.5 and RCP8.5 scenarios, respectively. No major precipitation changes are projected in the RCP2.6 scenario.


Introduction
Global warming brings increased occurrences of extreme weather and climate worldwide [1]. Significant intensification of heavy precipitation can be seen on global and continental scales [2,3]. Consistently, southeastern South America (SSA) tends towards warmer and moister conditions, with a noticeable increase in the magnitude and frequency of extreme temperatures and heavy rainfall (e.g., [4,5]). In particular, central-northeastern Argentina has been experiencing these significant climate changes since the 1960s [6,7]. The continuous increase in frequency and in intensity of heavy rainfall has caused great losses in agro-industrial production, severely damaged water management, and affected the urban centers in the region [8,9]. However, human activities are also impacted by increased occurrence of consecutive dry days, which lead to more seasonal droughts [10,11] as well as by extreme-maximum and minimum-temperature events, which have grown in frequency, intensity, and duration [12].
The Humid Pampa extends over central-northeastern Argentina, and its core crop region is centered in the southern half of the province of Santa Fe (see Figure 1). The core crop region is acknowledged as one of the most productive agricultural lands in the world, thanks to the combination of its climate and the characteristics of soil sediments. Lovino et al. [13] have reported enhanced annual precipitation since the 1960s and Extreme precipitation may lead to floods and droughts that have strong negative impacts on human activities and the environment. Floods occurred in the core crop region in 2018-2019. Immediately following the flooding period, a severe drought started in 2019 and lasted during 2020. This drought took place over a large part of SSA and caused extraordinarily low water levels in the Parana River [15,16]. The Rosario Stock Exchange (BCR) estimated a loss of USD 244 million in the agro-export sector of Rosario in the first quarter of 2020 [17,18]. The sector had already faced losses of USD 230 million as a consequence of the 2018-2019 flood [19]. Moreover, low water in the Paraná River also affected the river ecosystem and other essential activities such as fishing and water purification [17].
This study examines climate variability, climate extremes as well as observed and projected climate changes in the core crop region of the Humid Pampa in the 20th and 21st centuries. The resulting information might contribute to estimate the impacts of climate variability and change aiming at reducing disaster risks and implementing appropriate management policies and adaptation actions in the affected sectors. The final goal is to make science-based information available to policy and decision makers in the public and private sectors for the development of policies, actions, and plans aimed at reducing the vulnerability of the region to extreme hydroclimatic events. To fulfill this objective, we analyzed the variability and changes in temperature and precipitation using high-quality meteorological records and future climate projections for the core crop region of the Humid Pampa. The paper is organized as follows; Section 2 describes the methods and the data used. The observed changes in temperature and precipitation from the beginning of the 20th century to present are described in Section 3.1. For the same period, we identified and examined climate extremes using different indices, mainly focused on intense precipitation and extreme temperature events, which are in Section 3.2. Future climate projections in the short-, medium-, and long-term for the 21st century are presented in Section 3.3. The conclusions of the study are in Section 4.

The Relevance of the Core Crop Region
The Humid Pampa is an extensive region in central-northeastern Argentina. Its core crop region includes the center and south of Santa Fe province, east of Córdoba province, and north of Buenos Aires province (see Figure 1). The Humid Pampa presents a flat landscape with dominance of Quaternary aeolian sediments. Its prairie-like ecosystem develops in humid to sub-humid temperate to subtropical climates (North Pampa, [20]). The southeast of the Santa Fe province extends over the Pampean loessic belt [21,22], which is the most extensive in the Southern Hemisphere and one of the thickest loess deposits in the world [22]. The characteristics of the Pampean loess and the present climate make soils of the region highly productive and known worldwide for their fertility [21,22].
The Paraná River, which drains the second largest basin (2,600,000 km 2 ) in South America, flows across the Humid Pampa. The city of Rosario is located in the southeast of Santa Fe, on the right margin of the Paraná River, i.e., the core crop region of the Humid Pampa (see Figure 1). Rosario rests over a complex floodplain extending 600 km in the heart of the inland lowlands of the continent [22]. Tributaries to the Paraná River in the area of Rosario are the Carcarañá and the Saladillo Rivers (see Figure 1). Rosario has a privileged location both from the geographic and environmental perspectives. It is the most important city in the core crop region and one of the most important in Argentina, as it is a major industrial, trade, and financial center [23], with 43% of the province population. It is also the main agro-industrial and port hub in the region, and in recent years, it became the largest soybean and by-product export node in the world, with 78% of Argentina's installed capacity [24]. The port of Rosario is connected to the north with Mercosur countries through the Parana-Paraguay waterway, which is the primary communication corridor with southern Brazil and Paraguay [25].

Observed Data
Temperature and precipitation were examined using a long-term record of observations from Rosario Aero meteorological station (32.92 • S, 60.78 • W) of the National Weather Service of Argentina. This station, located in the core crop region of the Humid Pampa, is the only weather station in the region with long, complete, and reliable records of observations. These observations started in the first decades of the 20th century, thus allowing an analysis of long-term variability and historical changes. As a note of caution, the long-term variability as well as the historical and projected trends here reported are localized over the core crop region and do not represent the entire Humid Pampa region.
Daily temperature observations cover the period from 1935 to 2019, daily precipitation from 1936 to 2019, and monthly precipitation from 1911 to 2019. The temperature and precipitation time series meet the quality control requirements established by the World Meteorological Organization [26]. The time series of temperature and precipitation have very few missing data (2.8% for daily temperature, 0.25% for daily precipitation, and 3.1% for monthly precipitation). Data gaps were filled using the normal ratio method [27] and linear regression [28]. We mainly used the normal ratio method, while linear regression was applied when there were good correlations between station time series, or when data were available from only one neighboring station.

Detecting Variability and Changes in Time Series
A singular spectrum analysis (SSA, [29,30]) was applied to study the temporal variability of temperature and precipitation time series, including extreme events. The method determines the spatio-temporal structures of trends and quasi-oscillatory modes in time series (see details in [10]). SSA makes it possible to identify nonlinear long-term trends that describe the changes in temperature, precipitation, and extreme climate indices. Trend significance at the 95% confidence level was tested against a red noise null hypothesis using Monte Carlo method [31] with an ensemble of 1000 independent realizations.
For the long period (1911-2019) of monthly precipitation, we identified the years in which precipitation was above (below) the 90th (10th) percentile in summer (DJF), autumn (MAM), winter (JJA), and spring (SON), for the entire period of analysis. The percentile threshold selected enabled us to obtain a representative sample of excess/deficit events in each season of the year. Precipitation series were detrended applying linear regression. The relationship of these events with ENSO phases was analyzed using the classification by Trenberth [32], based on the Oceanic Niño index (ONI), i.e., the 3-month mean SST anomaly for the Niño 3.4 region [33]. Since ONI starts in 1950, for the previous years, we adopted an ENSO classification presented in Berri et al. [34].

Extreme Climate Indices
The Expert Team on Climate Change Detection and Indices (ETCCDI) defined a set of descriptive indices for temperature and precipitation extremes, which allows analyzing their frequency, intensity, and duration [35]. ETCCDI indices enable a uniform perspective of observed changes in climate extremes. In this study, a subset of ETCCDI indices was computed for Rosario Aero data (see Table 1). Our analysis includes fourteen temperaturebased and five precipitation-based indices as the most relevant for extreme events having greatest impact in northeastern Argentina (see [10]), where the study region is located.  10) are calendar days in the 90th (10th percentile) of daily maximum (minimum) temperature calculated for a 5-day window centered on each calendar day in the base period 1981-2010 (adapted from Donat et al. [36]).

Index
Index Name Index Definition Unit The fifth phase of the Coupled Models Inter-Comparison Project (CMIP5) used in the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC) provides historical and future climate simulations. One type of the CMIP5 experiments consisted of long-term historical simulations on the centennial time scale [37]. Longterm historical simulations allow for assessing the model ability to represent historical and present climate together with observed changes. In addition, they provide initial conditions for climate projections under representative concentration pathway scenarios (RCPs, [38], defined in Section 2.3.2 herein). Those long-term historical simulations were carried out for the period 1850-2005-i.e., from the beginning of the industrial revolution era to the near present. Lovino et al. [39] evaluated the ability of 25 CMIP5 models to represent historical and present climate in central-northeastern Argentina during the 1901-2005 period ( Figure 1). They selected the best-performing models, which include 9 models for temperature and 7 models for precipitation (see Table A1).
The General Circulation Model (GCM) simulations used by Lovino et al. [39] were corrected using a scaling method [40] to remove systematic errors in the simulated variables by fitting them to observed historical values from the Rosario Aero series. The scaling method (also called "direct approach") used in this study has been applied successfully in different previous precipitation and temperature studies (e.g., [41][42][43]). The method uses a scaling factor defined as the ratio between observed and simulated mean monthly variables for the period estimated in the historical simulations. Thus, Equation (1) shows that the corrected variable y  (1) In summary, the scaling method was used to correct systematic errors of the historical simulations and also adjust the trajectory of future simulated variables to that of historical observations. Bias correction was used to build the area-averaged time series of historical simulations and future projections. We also calculated the spatial distribution of future changes with non-corrected multi-model ensembles, in both historical and future simulations (see Section 2.3.2). Thus, the estimated changes are not affected by the bias correction method.

Future Simulations
Different climate change scenarios were used to estimate future climate trends in precipitation and temperature. Such a modeling exercise required assessing the historical variability and extreme events of these variables. To do so, we used the multi-model ensembles of GCM future projections that best represented historical and current climate in the region. The IPCC representative concentration pathway scenarios (RCPs, [38]) used here were: mitigation (RCP2.6), stabilization (RCP4.5), and high emissions (RCP8.5) were applied to the models listed in Table A1 to simulate monthly temperature and precipitation. No RCP scenario projections were available for three of the models used by Lovino et al. [39], i.e., CESM1-FASTCHEM, MIROCh4, and EC-EARTH. Thus, we used projections from 7 GCMs for temperature and from 5 GCMs for precipitation. Projections for RCP4.5 and RCP8.5 were available from all these models. However, RCP2.6 was available only for 3 models for each variable. Projected time series of temperature and precipitation from each GCM and the average GCM ensemble-called multi-model ensemble-were estimated for the closest grid point to Rosario Aero in the short (2006-2035), medium (2036-2065), and long term (2066-2100). Mean changes were analyzed in each model projection and time period.
Moreover, mean spatial changes were also analyzed for central-northeastern Argentina ( Figure 1) using the multi-model ensemble. The spatial distribution of mean expected changes was examined for the three climate change scenarios in the periods 2006-2035, 2036-2065, and 2066-2100. This required re-gridding all the GCMs to a common grid with horizontal resolution 1 • longitude-latitude, and time averages were calculated within the region bounded by 36 • S-26 • S latitude and 65 • W-58 • W in longitude ( Figure 1). Changes were estimated for each GCM, and the ensemble of changes was generated with all the models at every grid point.

Results and Discussion
3.1. Variability and Trends 3.1.1. Temperature Figure 2 shows the annual averages of the mean, maximum, and minimum temperatures at Rosario Aero station in the 1936-2019 period, together with their significant nonlinear long-term trend. A slight increase is observed in the mean temperature ( Figure 2a) in the first decades of the 20th century. This warming reached 1 • C from the 1970s to present, from an average value of 17.3 to 18.3 • C. Annual average minimum temperature ( Figure 2b) showed a sustained rise throughout the period, although it was more marked in recent decades, reaching a warming of 2 • C (from 10.5 to 12.5 • C). Annual average maximum temperature ( Figure 2c) decreased very slightly from the beginning of the analyzed period to the end of the 1970s; but from then on, it progressively increased by about 0.5 • C. Our results suggest that the minimum temperature recorded a fourfold increase compared to maximum temperature. Consistently, minimum temperature increased globally at rates more pronounced than maximum temperature (e.g., [1,36]). The joint analysis of the three temperatures shows that it was minimum temperature that drove the rise in mean temperature, rather than maximum temperature (2 and only 0.5 • C, respectively). Interestingly, Rosario Aero station is located at an airport, in an area that has undergone increasing urbanization over the years. As a consequence, the urban heat island effect could also be causing an increase in minimum temperature, mainly because of limited nighttime cooling.
Seasonal time series of mean, maximum, and minimum temperatures were also analyzed for the same period, along with their resulting nonlinear trends (figures not shown). No significant trends were observed in autumn and winter mean temperatures. On the other hand, we found significant trends in summer and spring, when the largest increase in average temperature-of the order of 2 • C-occurred. In summer, this increase would be explained mainly by an increase of 3 • C in minimum temperature, since the maximum temperature showed almost no variation. In spring, warmer minimum temperature is the major contributor to mean temperature increase, although there is also some contribution from maximum temperature.
The difference between maximum and minimum temperature, that is, the diurnal temperature range (DTR), remains practically unchanged during autumn and winter throughout the analyzed period. Spring and particularly summer DTR has been decreasing over the years. Summer DTR decreased from 15 to 12 • C. This result agrees with the findings by Lovino et al. [10], who reported a DTR decrease in northeastern Argentina during summer. Figure 3 presents the annual cycle of monthly precipitation as well as the average and maximum values of monthly maximum 1-day precipitation (RX1day) at Rosario Aero station in the 1911-2019 period. The annual cycle of precipitation displays significant intra-annual variability. According to Doyle [44], Rosario has a precipitation regime with a short period of low precipitation in winter and fairly uniform precipitation in the rainiest months. Figure 3 shows that the rainy season runs from October to April, when it rains 78% of the annual average. The main forcing factors during the rainy season include the South American Low-Level Jet (SALLJ) east of the Andes [45,46] and the Mesoscale Convective Systems (MCSs). The SALLJ is a key factor that activates convection and favors precipitation in the subtropical plains of South America [47]. MCSs are frequent and account for a large percentage of total precipitation during the warm rainy season (October-April) over northeastern Argentina [48]. March is the rainiest month, with an RX1day average of 59 mm (Figure 3, light blue dot) and an absolute maximum of 179 mm-the same as the maximum RX1day for December ( Figure 3, red dots). Summer months (DJF) exhibit precipitation values close to 110 mm/month and average RX1day of near 43 mm along the study period (Figure 3, light blue dots). Winter months (JJA) registered monthly values of 30 mm and an average RX1day of nearly 19 mm. Low winter precipitation is a common feature throughout the region [13,44]. In this season, the northward displacement of the South Atlantic anticyclone restricts water vapor inflow over the south-central region of the continent, which, added to lower solar radiation, results in the weakening of convective processes and heavy rains [49,50].

Precipitation
The connection was examined between years of absolute maximum daily precipitation in the whole series (1936-2019) and the occurrence of ENSO (Figure 3    The large inter-annual variability observed in Figure 4, including water deficits and excesses, is characteristic of the central-northeastern region of Argentina and has cycles of 2-6 and 9 years [7]. As shown by Lovino et al. ( [10] and references therein), most excesses were related to EN years, while deficits were related to LN years. However, they reported that excess precipitation was observed to be more localized, rather than covering the entire study area. On the contrary, precipitation deficit caused droughts with large spatial extension ( [11] and references therein). To assess the connection between precipitation and ENSO at Rosario Aero, the years of excesses and deficits (90th and 10th percentiles, respectively) were identified according to the seasons of the year, for the entire observational period available (1911-2019). The ONI index from 1950 was applied, according to Trenberth [32], to classify years into EN, LN, or Ne. Years of excess/deficit before and after 1950 are presented in Table 2.  1943 1913, 1914, 1926, 1933, 1938, 1947/ 1914, 1922, 1940, 1949/1916, 1920, 1933, 1935, 1937, 1944, 1948 1911, 1919, 1933, 1939/1916, 1917, 1949 After 1950 EN 1977/ 1966, 1969/ 1965, 1972/ 1972/ LN 1956, 1971, 1976, 1984/2018/ 1973/1974, 1999Ne 19611959, 2007/ 1959/1995, 20081978, 1990, 2001/1966 The analysis of the whole observation period shows that precipitation excess/deficit years were, respectively, 11/3 in summer, 11/0 in autumn, 11/10 in winter, and 11/6 in spring. In other words, except for winter, there were much more years with excess precipitation than with deficit. In particular, no precipitation deficits were observed in the autumns of the 109 years analyzed. Table 2 shows that there were much more excess periods than deficit periods. In addition, the period with water excess has grown in the last two decades with the relationship 7/3 years of excess/deficit. Deficit was more common in winters and springs, although less common than excess, in agreement with Sgroi et al. [11]. Table 2 shows a predominance of years with summer precipitation excess under LN conditions. On the other hand, winter excess occurred under EN conditions, and spring excess was mostly related to Ne years. No remarkable relation was identified between autumn precipitation and large-scale forcing. Contrary to literature findings for the region (e.g., [51,52]), no summer deficits were recorded under LN, but mostly excesses. This result would indicate that mesoscale processes are the drivers of excess precipitation in summer (as noted above). In contrast, winters with excess precipitation occurred during EN conditions, that is, they were modulated by large-scale forcing, in line with what has been found for central-northeastern Argentina [53]. No winter deficit years were observed during EN and LN, and the few deficit years identified were Ne. In the transition seasons, no years were recorded with deficit in autumn; spring deficits occurred in LN and Ne years ( Table 2). Note that ENSO-related SST anomalies generally begin in the first half of a calendar year, they peak by the end of that year, and start to vanish in the following year [32]. Therefore, a complete ENSO cycle includes two consecutive years, and its effects are felt at different moments of the event cycle to a variable degree in the regions affected [34].

Frequency of Temperature Extremes
Figure 5a-d presents the changes in occurrence frequency of temperature extremes, i.e., warm days (TX90p, Figure 5a), cold days (TX10p, Figure 5b), warm nights (TN90p, Figure 5c), and cold nights (TN10p, Figure 5d). The frequency of warm days (TX90p, Figure 5a) remained approximately constant during the first decades. Starting in the 1990s, the average frequency of TX90p has grown by approximately 5%, from 10 to 15%, warm days per year in the present. A sustained positive trend in warm night frequency (TN90p, Figure 5c) was observed except at the beginning of the period, since the annual average of warm night frequency increased from 7 to 12% between 1960 and 2019. The frequency of occurrence of cold days (TX10p, Figure 5b) has been decreasing throughout the whole period, while the frequency of occurrence of cold nights showed a noticeable drop (from 18 to 8% annually, see Figure 5d), particularly from the 1960s onwards. Such a remarkable drop could be influenced by the gradual increase in minimum temperatures (Figure 2b).  Table 1 for index definition. All nonlinear trends detected with the SSA method are significant at the 95% confidence level.
These changes in frequency of temperature extremes are consistent with global trends [35,54]. At regional scales, the changes are also consistent with the results in Rusticucci et al. [55], who analyzed temperature variability in central-eastern Argentina. They reported a decrease in the frequency of cold extremes-TX10 and TN10-and a significant increase in warm extreme occurrence (TX90) in the period 1981-2010. Other authors suggested that temperature-related extremes were influenced by interannual variability (Ref [12] and references therein). Collazo et al. [56] reported interannual cycles in warm-day occurrence with periodicities of 2-4 years and 8 years. They also found that ENSO and the Pacific Decadal Variability (PDO) were key predictors of warm days in Argentina. Moreover, Müller et al. [49] and Collazo et al. [57] suggested that El Niño favored the occurrence of less cold extremes. In particular, the occurrence of El Niño during autumn favored the occurrence of less extreme cold events in the following winter throughout the region, as reported in Collazo et al. [57].  Figure 6a,b show the temporal evolution of the highest and lowest annual temperatures in the time series of maximum temperature. Figure 6c,d show the same for minimum temperature. The annual warmest maximum daily temperature (TXx, Figure 6a) showed a decrease from the beginning of the period until the 1980s, when it began a slight increase that continues to the present. From end to end, maximum annual temperature decreased by 2 • C throughout the analyzed period. The coldest maximum daily temperatures (TXn, Figure 6b) showed a minor oscillation with no resulting change throughout the study period. On the other hand, the warmest minimum daily temperature (TNx, Figure 6c) presented a sustained increase of over 2 • C that accelerated in the last decades. Finally, the coldest minimum daily temperature (TNn, Figure 6d) exhibited a slight increase of almost 1 • C. The range of warmest daily temperatures (Figure 6a,c) decreases from 16 • C in 1936 to 12 • C in 2016; but the amplitude of the coldest daily temperatures remains basically unchanged around 14 • C throughout the whole period (Figure 6b,d).
Changes in duration of extreme temperature values are presented in Figure 6e for warm periods and in Figure 6f for cold periods. The duration of warm spells (Figure 6e, WSDI) showed a tendency to almost double the WSDI in the long term, from an average of 5 days to almost 10 days in duration. Conversely, the duration of cold spells (Figure 6f, CSDI) decreased systematically, and in recent decades, it barely reached 1 to 2 days on average. Alexander et al. [59] also analyzed warm spell duration in Argentina in the period 1951-2003. They found positive and negative trends over different regions in the country, with no significant linear trends of WSDI. More recently, this result was reinforced by two studies, which used global data indices (HadEX2). The first of these studies, by Donat et al. [60] found the same spatial extent and trend inconsistency that Alexander et al. [59] found in the same region. The other study by Lu et al. [61] analyzed the duration of warm and cold spells globally, also using HadEX2 data from 1958-2010, and determined weaker trends for both indices in South America: positive for WSDI and negative for CSDI. Finally, based on observations, Lovino et al. [10] found a significant trend toward longer warm spells and shorter cold spells in the 1963-2013 period in northeastern Argentina.  Table 1 for index definition. All nonlinear trends detected with the SSA method are significant at the 95% confidence level.

Intensity and Duration of Precipitation Extremes
Figure 7a-c shows the time series of the indices of precipitation intensity extremes, i.e., the annual maximum 1-day precipitation (Figure 7a, RX1day), the annual maximum 5-consecutive day precipitation (Figure 7b, RX5day), and the simple daily intensity index (SDII), defined as the total annual precipitation divided by the number of rainy days in the year. The RX1day (Figure 7a) remained approximately constant for most part of the period, but in recent decades, it suggested a moderate increase. The RX5day (Figure 7b) showed an increase at the beginning of the study period, it then stabilized until recent decades when it increased again. The SDII (Figure 7c) showed a significant increase during the last three decades, in agreement with the increase observed in annual precipitation in recent decades (Figure 4). This result was also found in the analysis of Table 2, where 7/3 years of excess/deficit corresponded to the last two decades. Similarly, SDII trend values for Rosario Aero station were coincident with the slightly upward trends calculated/projected in southeastern South America [56,59]. In general, the whole of South America is becoming wetter, which is related to the occurrence of heavy rainfall events concentrated in few days [62].  Table 1 for index definition. All nonlinear trends detected with the SSA method are significant at the 95% confidence level. Figure 7d,e shows the duration of precipitation extremes characterized by the annual maximum number of consecutive dry days (Figure 7d, CDD) and the annual maximum number of consecutive days with precipitation greater than 1 mm (Figure 7e, CWD). Except for the short decrease during the first two decades, the number of consecutive dry days (Figure 7d) has increased steadily from the 1950s, from 30 to 40 days. In contrast, consecutive wet days (Figure 7e) remained rather constant throughout the whole study period, with a mean value slightly above 4 consecutive days. A recent study considering consecutive days with precipitation greater than 0.1 mm in central-northern Argentina showed that precipitation events of 6 consecutive days were not frequent [63].

Future Climate Projections
Future climate projections for temperature and precipitation were analyzed using GCM models (Table A1) and their multi-model ensembles, for different time horizons in each of the RCP2.6, RCP4.5, and RCP8.5 scenarios.
3.3.1. Temperature Figure 8 presents the time series of the multi-model ensembles that provide a graphic view of expected temperature variations in future climate change scenarios. In the RCP2.6 mitigation scenario (Figure 8a), the multi-model ensemble showed that there would be a small Tmean variation, with average values very close to current mean values in 2035, 2065, and 2100 (vertical lines in the figure). The inter-model range (shading in Figure 8a) in future projections is very similar to that of the historical period and relatively small, which makes the projection of the multi-model ensemble more robust.  In the scenario of high emissions RCP8.5 (Figure 8c), Tmean would grow steadily, until an average maximum value of 21.5 • C at the end of the 21st century. In the short term, until 2035, an average Tmean close to 18.5 • C would be expected, a similar value to that expected in the stabilization scenario, with an average of 19.5 • C in 2065. The inter-model range is small (Figure 8c), suggesting a coherent projection of the different models, which gives robustness to the results. Figure A1 presents expected average temperature changes in the short-, medium-, and long-term future periods, according to model projections for different scenarios. In the RCP2.6 scenario ( Figure A1a), positive changes of less than 1 • C are expected for the short (to 2035), medium (to 2065), and long term (to 2100). Only one of the analyzed GCMs projected negative changes of less than 0.5 • C. The multi-model ensemble showed changes of 0.4 • C by 2035 and up to 0.6 • C by 2100. In summary, in the mitigation scenario, the expected Tmean increase would be less than 1 • C for the whole 21st century. In the RCP4.5 scenario, all models and their multi-model ensemble showed a consistent increase of 0.6 • C by 2035, of approximately 1.2 • C by 2065, and between 1.5-2 • C by 2100 ( Figure A1b). There is high consistency among the average changes projected by each model, projecting robust changes for this scenario in each period. In the RCP8.5 scenario (Figure A1c), all models and their ensemble also showed very consistent positive changes among the different GCMs. For the 2006-2035 period, an increase is expected of about 0.6-1 • C, of 1.7-2 • C for 2036-2065, and of up to 3-3.5 • C for the 2066-2100 period. These results suggest that average annual temperature would increase significantly during the 21st century in the worst-case emission scenario. Figure 9 shows the spatial distribution of the mean changes in mean annual temperature estimated with the multi-model ensemble in central-northeastern Argentina. The greatest increases were observed in the north of the region in the three scenarios, with decreasing magnitude towards the south. This rise could reach 4 • C by the end of the century in the high greenhouse gas emissions scenario (RCP8.5). In particular, a warming of 0.5-1 • C was projected under the three climate-change scenarios in the core crop region in the 2006-2035 period (first column, Figure 9). Thus, the region would experience a warming between 0.5 and 1 • C by 2035, regardless of the scenario. However, in the following 2036-2065 period, mean temperature would rise 1-1.5 • C in the RCP2.6 and RCP4.5 scenarios, and between 1.5 and 2 • C in the RCP8.5 scenario (second column, Figure 9). By the end of the century-2066-2100-(third column, Figure 9), models projected changes of 1-1.5 • C in the mitigation scenario, 1.5-2 • C in the stabilization scenario, and 3-3.5 • C in the high-emission scenario.  Figure 10 presents the multi-model ensemble projections for annual precipitation. Precipitation would remain almost the same in the mitigation scenario (Figure 10a), with no significant differences in the short-, medium-, and long-term periods. The inter-model range is low, which gives greater reliability to the ensemble projection. In the stabilization scenario, precipitation might experience a less significant increase (Figure 10b). The intermodel range for the GCMs studied fluctuates around the average values, especially after 2065. In this scenario, annual precipitation is expected to increase to a mean value of approximately 1100 mm by 2100, while the mean value for the period 1971-2005 is 1020 mm. In the high-emission scenario, an increasing trend of annual precipitation is observed, although not as significant as that for temperature. In this scenario, annual precipitation could reach 1150 mm by 2100. The range of variation among GCMs was considerable, although it did not exceed the variability of the observed historical time series. These results suggest a gradual increase in precipitation along the 21st century in the RCP8.5 scenario. Our results show that none of the models predicted a trend in precipitation as noticeable as in temperature. Changes in precipitation would not be expected to affect annual amounts significantly, although they would likely impact extreme events in southeastern South America, as suggested by Seneviratne et al. [64], Donat et al. [65], and Pfahl et al. [66]. Therefore, the increasing trend observed in intense precipitation events (see Section 3.2.3) would likely continue in the coming decades. In the three emission scenarios, annual precipitation would not undergo significant changes in the study region. In the most unfavorable emission scenarios (RCP4.5 and RCP8.5), there would be slight increases in precipitation, and the variability would remain practically unchanged. Figure A2 summarizes the annual precipitation changes in future scenarios for the different studied periods. Precipitation changes simulated by the GCMs were not as consistent as those of temperature. In the mitigation scenario ( Figure A2a), variations in annual precipitation would be rather small, since expected changes would vary around 20-30 mm in the 2006-2035 period and around 20-60 mm in 2036-2065 and 2066-2100. These variations represent a 2 to 6% increase in annual precipitation, which is not significant given the great interannual variability of precipitation in the region. In the RCP4.5 scenario ( Figure A2b), precipitation would increase between 40 and 90 mm from 2006 to 2065. In the period 2066-2100, a greater increase, between 70 and 130 mm, would be expected, which represents 10% of current values. The only GCM model of the group that projects a minimum and slightly negative change is NorESM1-M. The most significant increases would be expected in the high-emission scenario RCP8.5 ( Figure A2c). In the 2006-2065 period, the range of variation among models oscillates between 40 and 140 mm (except NorESM1-M). In the 2066-2100 period, GCMs showed consistent increasing projections of 80-160 mm. This increase represents from 8 to 16% of average annual precipitation. Figure 11 presents the spatial distribution of mean changes in annual precipitation in central-northeastern Argentina. The distribution of mean changes is heterogeneous, a characteristic feature for this variable. Changes were mostly positive, indicating an increase in precipitation. However, under scenario RCP2.6 (first row of Figure 11), some regions in the northwest and the south will undergo negative changes-i.e., precipitation will slightly decrease throughout the century. In the rest of the region, increases will not exceed 50 mm. Under scenarios RCP4.5 and RCP8.5 (second and third rows of Figure 11), projected precipitation increases were greater than those under RCP2.6. Mean changes in precipitation will be greater in the west, center, and south of the region, with peak values of 100 to 150 mm by the end of the century in the high-emission scenario RCP8.5. In particular, in the core crop region around Rosario (Figure 11), projections under the mitigation scenario show a precipitation increase of up to 25 mm in the period 2006-2035 and from 25 to 50 mm for the rest of the 21st century. These values are not significant given the large interannual variability of precipitation (Figure 4a). In the stabilization scenario, precipitation will increase between 25 and 75 mm in the 21st century. The high-emission scenario projects the greatest increase, between 50 to 100 mm, in the first half of the century, and between 100 and 150 by the end of the 21st century. These projections would represent an increase of from 8 to 16% of current values in mean annual precipitation.

Precipitation
Several studies have shown that globally averaged precipitation increases with warming [2,3,67,68]. However, at regional scales, the situation may be different, as warming may result in more precipitation in one region and less in another. This lack of connection has been noted in other studies that attempted to evaluate such links in southern South America, where precipitation is projected to increase [2,67,69,70].
These results suggest that changes in precipitation would not be expected to affect annual amounts significantly. However, it would be highly likely that the frequency and intensity of extreme precipitation events would be affected [9,64]. The analysis of historical records showed increasing trends of precipitation intensity (see Figure 7). Should this trend continue in the coming decades as expected, intense precipitation events would continue to be more frequent in the region.

Conclusions
The core crop region of the Humid Pampa is one of the richest agricultural lands around the world and the main productive region in Argentina. It concentrates most of the country's agro-industrial production. The analysis of the climate changes in the 20th and 21st centuries is crucial, because besides its economic importance, the region is densely populated and is home to one of the most important cities in Argentina. This paper assesses climate variability and historical climate extremes from the beginning of the 20th century to the beginning of the 21st century and historical changes in temperature and precipitation. It also includes future simulations of temperature and precipitation under different emission scenarios in the 21st century, for the short, medium, and long term.
The analysis of temperature (1936-2019) revealed a 1 • C rise in annual mean temperature by from 1970 to present. Mean annual values of minimum and maximum temperature have risen by 2 and 0.5 • C, respectively. These changes occurred mainly in spring and summer, while no significant trends were seen in autumn and winter. From the beginning of the 20th century, annual precipitation (1911-2019) rose from an average of 950 mm in 1960 to 1050 mm between 1970 and 2017, which represents a 10% growth in average precipitation. The greatest change (almost 100 mm/season) occurred in summer (wet season), while winter precipitation (dry season) decreased slightly in recent years.
During the 109 years of monthly precipitation observations (1911-2019), 11 excess events were observed in each season of the year. In contrast, years with precipitation deficit showed notable differences depending on the season. There was no precipitation deficit in autumn, but deficits were identified in three summers, six springs, and up to 10 winters. These results reveal that the number of excess periods is much greater than the number of deficit periods. In addition, the number of excess periods has increased faster in the last two decades with a relationship 7/3 years of excess/deficit.
The effects of ENSO on precipitation depend on the season of the year. In summer, there is a predominance of years with excess precipitation in LN years, while winter excesses predominate in EN years. Absolute maximum daily precipitations were recorded in two LN summers in the entire 1936-2019 period. These results suggest a dominance of mesoscale processes in the summers with excess precipitation. Such dominance could be caused by LN inhibiting large-scale precipitation mechanisms ( [51,52]). In contrast, large-scale processes predominated in winters with excessive precipitation during EN, as it enhanced precipitation in the center-northeast of the country. In the transition seasons, ENSO seemed to have no influence on excess precipitation, while precipitation deficits during LN events were observed only in two springs and in one summer.
Temperature extremes showed significant trends towards warmer conditions, consistent with observed global trends ( [36,59]). The frequency of occurrence of warm days and nights increased considerably and cold days and nights decreased notably, with a significant drop in the number of frost events. Consistently, the duration and intensity of heatwaves increased. The opposite was seen for cold waves. Extremes of precipitation increased in frequency, with events lasting from 1 to 5 consecutive days, while the number of consecutive dry days also increased. These changes led to a concentration of annual precipitation in fewer days, which caused more intense events.
Future climate scenarios were estimated from a subset of GCMs selected in Lovino et al. [39]. Model projections indicated that temperature would continue to rise in central-northeastern Argentina during the 21st century. The growth rate is higher towards the north of the region and decreases towards the south. The warming would not exceed 1 • C in the first period analyzed (2006-2035), regardless of the scenario. From 2035, temperature would be expected to rise, reaching 2 • C at the end of the 21st century under the stabilization scenario. In the high-emission scenario, the rate of increase grows, with temperature changes reaching 3-3.5 • C.
Future average precipitation changes would have heterogeneous spatial distribution throughout central-northeastern Argentina. Future projections of annual precipitation showed gradual changes with a slight to moderate increase depending on the scenario. In the mitigation scenario, total annual precipitation would not undergo major variations during the 21st century. Average precipitation changes under the stabilization and highemission scenarios would be greater in the west, center, and south of the region, including the area of Rosario. Values could range from 100 to 150 mm by the end of the century in the high-emission scenario. In particular, in the core crop region, annual precipitation is projected to increase during the century in all the cases. The more critical the scenario, the greater increases in projected precipitation.
Our results of observed and projected changes in regional climate and the increase in extreme event occurrence show that new management strategies are needed for natural and human systems in order to prevent life and economic losses, damage to infrastructure, and natural resources. In addition to tending towards sustainability, this information would allow for improving the management of climate risk by adopting preventive measures to face the consequences of climate change and variability in the study region.

Data Availability Statement:
The CMIP5 multi-model ensemble [37] data presented in this study are openly available in https://esgf-node.llnl.gov/search/cmip5/ (accessed on 26 February 2019). Restrictions apply to the availability of Rosario Aero meteorological station data. Data was obtained from the National Weather Service of Argentina and are available upon request from https://www. smn.gob.ar/ (accessed on 26 February 2019). Table A1. Selection of best-performing CMIP5 models in northeastern Argentina (see text for details). Variables: monthly surface air temperature (T) and precipitation (P).