Assessment of the RegCM4 Performance in Simulating the Surface Radiation Budget and Hydrologic Balance Variables in South America

The ability of the Regional Climate Model v4 (RegCM4) to simulate the surface radiation budget and hydrological balance variables over South America have been evaluated. For this purpose, a 34-year long simulation was carried out with the regional climate model RegCM4 over South America on the CORDEX domain. The model is forcing by ERA-Interim reanalysis. The results show that RegCM4 simulates the main patterns of the variables associated with the surface radiation budget and hydrological balance in the four seasons of the year compared to the observations (CLARA2 and CRU/PERSIANN). However, the cloudiness and surface radiation budget variables: Cloud Fraction Cover (CFC), net shortwave (SW) and longwave (LW) radiation at surface are overestimated, mainly over the oceans. This is associated with the errors in the CFC due to the deficiency of the model in representing the low-level clouds. Some differences are also noted in the hydrological balance. The intensity and temporal evolution of precipitation, especially in the central and southern Amazon, may be associated with the selected domain, which fails to adequately represent the influence of the adjoining oceans. In northern and northeast parts, the differences are associated with deficiencies of RegCM4 in representing precipitation rates. Although the deficiencies, taking into account that the model is capable to reproduce the general pattern of some important variables of the surface radiation budget and hydrological cycle, it may be a useful tool for climate studies.


Introduction
The surface radiative balance and hydrologic cycle are responsible for various atmospheric processes such as precipitation, evapotranspiration, wind intensity and temperature, in the regional and global scales (Wallace and Hobbs 2006). They are particularly important in tropical regions since the amount of radiation and the rates of precipitation are higher. South America (SA) region (between the latitudes 56° S-12 °N) is influenced by several processes that occur in the middle and tropical latitudes, as well as their interactions. SA is strongly affected by the adjoining oceans (Pacific and Atlantic) where the sea surface temperature (SST) variability controls the seasonal migration and interannual variation of the Intertropical Convergence Zone (ITCZ) and South Atlantic Convergence Zone (SACZ) (Garreaud et al. 2009;Reboita et al. 2010). These large-scale characteristics strongly affect the SA tropical and subtropical precipitation regimes. Also, the complex topography of the Andes mountains along the west of the continent, the vegetation of the Amazon basin and the large amount of radiation received by the tropical region maintain the SA weather and climate characteristics (Garreaud et al. 2009;Nobre et al. 2009;Reboita et al. 2010).
Some variables of the hydrological cycle is directly linked to the changes in the radiation balance and atmospheric temperature (Inglezakis et al. 2016). The northern part of SA, where the Amazon basin is located, is characterized by the large quantity of solar radiation received. The solar radiation received at the top of the atmosphere between 5° N and 10° S varies from the maximum of 424.8 Wm −2 in December-January to the minimum of 355.3 Wm −2 in June-July (Salati and Marques 1984). At the surface, the incident solar radiation is around 185.2-208.3 Wm −2 (Pereira et al. 2017). The seasonal cycle of incident solar radiation in the central 1 3 Published in partnership with CECCR at King Abdulaziz University Amazon shows maximum values in September-October and minimum values in December-February (Culf et al. 1996;Pereira et al. 2017). This temporal distribution is in the major part controlled by the nebulosity associated with Amazon convection (Horel et al. 1989).
The net surface shortwave (SW) and longwave radiation (LW) budget are important in climate modeling because they regulate the energy balance of the Earth and control the diurnal and annual cycles (Kothe et al. 2011). In the context of the CLARIS-LPB Project (Europe-South America Network for Climate Change Assessment and Impact Studies in La Plata Basin), Pessacg et al. (2014) evaluated the net surface shortwave and longwave radiation in the SA using seven regional models with a spatial resolution of around 50 km for the period from 1990 to 2008. Their results showed that in the most of the regional models the net SW is overestimated in the tropical regions of SA and La Plata basin, and underestimated over the ocean during the austral summer and winter seasons. The net LW did not show a defined pattern in the analyzed models, each one of them showing different results. On the other hand, the evaluations done by Lange et al. (2015) and Brisson et al. (2016) showed that parameterization of the convection and subgrid-scale clouds can improve the deficient representation of cloud cover, precipitation, and surface radiation budget, thus reducing biases. In the hydrological cycle, the precipitation is the meteorological variable more analyzed because it is the result of several physical processes including the thermodynamic, the changes in the water phase, and net radiation among others. In the Amazon, in the rainiest region of SA, the precipitation rates are maintained by the evapotranspiration, which is recycled by local sources (Salati et al. 1979;Zemp et al. 2017), the moisture carried from the Atlantic tropical ocean by the trade winds (Drumond et al. 2008;Duran-Quesada et al. 2009;Zemp et al. 2017) and the water vapor flux convergence associated to the regional and large-scale circulations (Chen 1985;Nobre et al. 2009). On the other hand, evapotranspiration is a dominant way in which the arriving net radiation and the precipitation are transferred from the surface to the atmosphere (Nepstad et al. 2004;da Rocha et al. 2009). On the regional scale, the moisture transport inside and outside from Amazon is extremely important for the rainfall regime, particularly during the rainy season. The moisture exported outside from Amazon is transported to the east of the Andes mountains by the SA Low-Level Jets (LLJs), with a maximum wind speed of 15 ms −1 at 850 hPa, near 17° S and 62° W. This contributes to the precipitation in the La Plata basin and southeast Brazil (Marengo 2004(Marengo , 2006Arraut and Satyamurty 2009;Martinez and Dominguez 2014;Montini et al. 2019).
Although Regional Climate Models (RCMs) are good tools for studying and evaluating the various processes that occur in the different components of the climate system, they still have some deficiencies in simulating the surface radiation balance and hydrological cycle processes in SA (Pessacg et al. 2014;Chou et al. 2012Chou et al. , 2014Reboita et al. 2014;Franchito et al. 2014;Lange et al. 2015;Llopart et al. 2017Llopart et al. , 2020Erfanian and Wang 2018). In other regions of the world RegCM4 also shows some deficiencies to represent the atmospheric variables Chiacchio et al. 2015;Almazroui et al. 2016;Nogherotto et al. 2016). Some studies showed that over Europe the RegCM overestimates the surface net SW radiation by 5 and 60 Wm −2 , respectively, in the boreal winter (December-January and February) and summer (June, July and August), while in the north Africa it is underestimated by − 30 to − 40 Wm −2 . On the other hand, the net LW radiation is underestimated by − 40 and − 30 Wm −2 , respectively, in most part of the model domain and with values 0 Wm −2 in the north Africa (Chiacchio et al. 2015). In a similar study, over the north Africa, Panitz et al. (2014) found that the COSMO model underestimates the net SW radiation until − 150 Wm −2 over the oceanic regions, in January-March and July-September. Over the continent, the values are underestimated by around − 50 Wm −2 . In the areas where the net SW flux was overestimated the values are around 30 Wm −2 . Otherwise, the net LW flux showed underestimated values by around − 30 and − 40 Wm −2 , and overestimated values with the same range of variation.
Particularly, for SA, there are few studies that have been conducted using regional models devoted to the concomitant validation of the surface radiation processes and hydrological cycle. For that purpose, the objective of the present work is to make a simultaneous analysis of these processes using the Regional Climate Model RegCM4 integrated over the Coordinated Regional Climate Downscaling Experiment (CORDEX) domain (Giorgi et al. 2009;Giorgi and Gutowski 2015). This study aims to contribute to a better knowledge of the processes associated with the surface radiation budget and hydrological cycle, especially the precipitation, over SA. Section 2 shows the RegCM4 description, the initial and boundary conditions for the model integration and the observed data used to examine the model performance in simulating the surface radiation balance and hydrological cycle over SA; Sect. 3 contains the results, and the discussions and conclusions are presented in Sects. 4 and 5, respectively. This paper becomes different from earlier studies since the role of the surface SW and LW radiation in the precipitation and potential evapotranspiration is emphasized.

Data
Data from the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim Reanalysis (EIN15) are used as the initial and boundary conditions in the numerical experiment. This reanalysis is available from January 1979 to August 2019 and has 1.5° × 1.5° (latitude, longitude) spatial resolution, 6 h temporal frequency and 37 vertical levels (Dee et al. 2011). An inter-comparison on reanalysis data made by de Almeida et al. (2018), indicated that in South America the EIN15 shows a better performance in relation to the NCEP-DOE Reanalysis 2 and Climate Forecast System Reanalysis (CFSR) version 1.
The hydrological cycle variables and the water balance (precipitation and potential evapotranspiration) are evaluated using the Climate Research Unit data (CRU, version 4.05), which have a spatial resolution of 0.5° × 0.5° (latitude, longitude) and monthly temporal frequency (Harris et al. 2014(Harris et al. , 2020. This dataset is available from 1901 up to the present time. Over the oceans, precipitation data from the Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks-Climate Data Record (PERSIANN-CDR, hereafter called as PERSIANN) develped by the Center for Hydrometeorology and Remote Sensing (CHRS) of University of California, Irvine (UCI) are used. This data set comprises the globe between the latitudes 60° S-60° N, with 0.25° × 0.25° (latitude × longitude) spatial resolution and is available from January 1983 up to the present time (Ashouri et al. 2015). To have the same spatial resolution the PERSIANN data was interpolated to 50 km using a bilinear method.
The radiation balance is evaluated using data from the Climate Monitoring Satellite Application Facility (CM-SAF), Cloud, Albedo, Radiation dataset, AVHRR-based, Edition 2 (CLARA2), derived from the AVHRR (Advanced Very High-Resolution Radiometer), measures the National Oceanic and Atmospheric Administration (NOAA) satellites and Meteorological Operational satellite program (Metop) in polar orbit, as performed by the European Organization for Meteorological Satellites Exploration (EUMETSAT) (Karlsson et al. 2017). The second version is an improvement of the first version, where original visible radiation is intercalibrated and homogenized using the Moderate Resolution Imaging Spectroradiometer (MODIS) data taken as a reference before the application of several parameters' retrievals. The data covers a 34-year period  in the daily and monthly scales and have a 0.25° × 0.25° (latitude, longitude) spatial resolution (Karlsson et al. 2017). The variables used for the analysis are CFC, albedo a net SW and LW radiation. It is available https:// doi. org/ 10. 5676/ EUM_ SAF_ CM/ CLARA_ AVHRR/ V002.
The equations for calculating the net SW and LW radiation are: where: SIS = incoming solar surface radiation, SAL = surface Albedo, SDL = surface downwelling longwave radiation, SOL = surface outgoing longwave radiation. Additionally, temperature at 2 m from the CRU dataset is used.

Methodology
The Regional Climate Model used in the present study is the RegCM4 version 4.5.10 from the International Centre for Theoretical Physics (ICTP) . The model dynamics is based on the Fifth Generation Mesoscale Model (MM5) of the National Center for Atmospheric Research (NCAR)-Pennsylvania State University (PSU) (Grell et al. 1994). RegCM4 is a limited area model for a compressible atmosphere, with Arakawa B-grid and sigma vertical coordinate. It can be used in any region of the world in a wide range of studies (study of processes, climate projections and even paleoclimate). Since its first version in 1989, the RegCM has gone through several updates with the implementation of new schemes of convection, surface, microphysics, etc. ). RegCM4 has two surface schemes. In the present study the Biosphere-Atmosphere Transfer Scheme (BATS) is used. Based on Reboita et al. (2014), in the present study the configuration of BATS with MIXED1 convection [Grell (1993) over land and Emanuel and Zivkovic-Rothman (1999) over ocean] is used. In the Conclusions of that paper, the authors recommended the use of this configuration of RegCM4.3 to carry out numerical experiments in the CORDEX domain over SA.
The RegCM4 is continuously running for a 34-year period  using EIN15 data in the CORDEX domain over SA, with the three first years as spin-up. To improve the model simulations, some parameters associated with vegetation type 6 (Amazon forest) in the BATS scheme are modified according to da Rocha et al. (2012) to obtain a better representation of the Amazon soil moisture. These modifications allow us to improve the evapotranspiration rates due to the greater availability of water in the soil column, as described in da Rocha et al. (2012). The modified values are given in Table 1.
The configuration of the RegCM4 used in this study is given in Table 2. Figure 1 shows the model domain of RegCM4, with a spatial resolution of 50 km. To obtain consistency of the analysis, the period of both observations and simulations is the same. So, in this work, the precipitation, potential evapotranspiration, CFC, albedo, temperature at 2 m and net surface LW radiation are analyzed from the period 1982-2012 [31 years (1979-1981 is the spin-up)]. While, the net surface SW radiation is considered from the period 1992-2012 (21 years), both for CLARA2 data and RegCM4 simulation. This is because of the lack of surface incoming SW radiation (the net surface SW radiation is calculated directly from the surface incoming SW radiation and albedo) in most part of tropical and subtropical regions during the period 1982-1991. The potential evapotranspiration calculation in RegCM4 is based on modified version of Penman-Monteith method, which is used in the preparation of CRU data (Harris et al. 2014). Figure 2 shows the comparison between seasonal variation of the horizontal winds simulated by the RegCM4 and the EIN15 at low (850 hPa, Fig. 2a-f) and high levels (200 hPa, Fig. 2g-l). In the austral summer (DJF), at low levels, near Equator, the RegCM4 simulates very well the easterlies trade winds over the Atlantic ocean. The trade winds are responsible for the advection and transport of water vapor from the Atlantic Ocean to the SA continent (Garreaud et al. 2009;Duran-Quesada et al. 2009;Zemp et al. 2017;Zanin and Satyamurty 2020). The Andes cordillera, located in the west SA, acts as an obstacle to the displacement of trade winds so that there is a shift southward along the east region of the mountains (Montini et al. 2019;Zanin and Satyamurty 2020). In the regions of LLJ and northwest Amazon the simulated values are underestimated up to − 4 ms −1 . In the northeast and north SA the RegCM4 overestimates the wind values up to 4 ms −1 . In the austral winter (JJA), the areas where the winds are underestimated and overestimated are the same as in the austral summer, but with low values. The northward and southward displacement of the Pacific and south Atlantic subtropical high is well reproduced by the RegCM4. There is a good agreement between the location of these subtropical highs simulated by the model and their climatological position (Sun et al. 2017;Fahad et al. 2021). At high levels, in the austral summer, the RegCM4 is capable for simulating the main synoptic systems, such as the Bolivian high and the associated cyclonic vortex, as well as the westerlies (Fig. 2h). The RegCM4 shifts westward the Bolivian high compared to that observed in EIN15 (Fig. 2g). Similar deficiencies in the Bolivian high location were found by Llopart et al. (2019). This deficiencies in the simulation of the Bolivian high provoke an anticyclonic circulation in the southwest side of Bolivia, with underestimated values of winds in the north of the circulation and overestimated values in the south of the circulation (Fig. 2i). Erfanian and Wang (2018) also showed that the deficiencies in the representation of the Bolivian high location causes an anticyclonic circulation at high levels in south Bolivia in December. In the austral winter, the westerlies are well simulated by the model. The position of strong westerlies (35° S-25° S) is well reproduced by the model. However, in the continental region of SA the winds are overestimated up to 6 ms −1 , while in oceanic regions the overestimated   In general, the wind pattern simulated by the RegCM4 is similar to that obtained in earlier studies (Fernandez et al. 2006;Seth et al. 2007;Lange et al. 2015;Reboita et al. 2018a), although the RegCM4 presents higher or lower performance in some regions. Figure 3 shows, respectively, the seasonal mean  of CFC for the CLARA2 data, RegCM4 simulation and the   for ERA-Interim (left column), RegCM4 (middle column) and difference between them (right column) for low levels 850 hPa (a-f) and high levels 200 hPa (g-l) for the austral summer (DJF) and winter (JJA). The red lines at the low levels represent the sea level pressure. The unit of wind is ms −1 difference between them. In the austral summer, which is the rainy season, a higher concentration of CFC is noted in the regions of the ITCZ and the SACZ (Fig. 3a, b). In general, there is a good agreement between the RegCM4 simulation and the distribution of cloud cover observed in CLARA2, mainly over the continent, although over the oceans the agreement is lower (Fig. 3c). In the central SA and parts of northeast Brazil, the RegCM4 underestimates the cloud cover (around − 20%), while in the regions near the eastern littoral of Brazil, parts of central Andes mountains, south and southeast Brazil and Paraguay, the cloud cover is overestimated (values around + 20%).

Surface Radiation Balance
Over the oceans, the RegCM4 underestimates the CFC with values up to − 30% over the north Atlantic ocean and Published in partnership with CECCR at King Abdulaziz University the Pacific ocean (near Peru and Chile). Simulation for the austral summer is better than the results presented in Lange et al. (2015) because the simulated values are lower (underestimated/overestimated), mainly in the continental region. In the austral autumn, the ITCZ is shifting northwards. From Fig. 3d, e, the region of higher CFC is located toward the north of around 15° S, although the RegCM4 values are underestimated in western Brazil and northwestern SA (values up to − 20% in some parts) and overestimated in the eastern and southern part of Brazil (values around + 20%) (Fig. 3f). Over the oceans, the spatial distribution of the CFC is similar to that noted in the austral summer. However, the RegCM4 simulation is underestimated (values lower than − 30% over the Pacific ocean and up to − 20% in the Atlantic ocean).
In the austral winter, the areas with higher CFC are located northward from 3° N (Fig. 3g, h) due to the northward shifting of the ITCZ as a response to the increase of the SST. In the continent, the higher CFC are seen northward from 5° S. Southward from 25° S the CFC is maintained mainly by the frontal systems that occur with more frequency in this season (Andrade and Cavalcanti 2004;Cardozo et al. 2015). The RegCM4 underestimates the CFC in the ITCZ regions (values lower and up to − 30% over the Pacific and Atlantic oceans, respectively). On the other hand, in most part of the Andes mountains region, central and east the RegCM4 simulation overestimates the CFC by + 5 to + 20% (Fig. 3i). For the austral winter, the results of Lange et al. (2015) are better than those of the RegCM4. This may be due to the types of RCMs and the periods of analysis used in these studies. In the austral spring, the ITCZ starts moving southwards so that the CFC configuration is similar to that observed in the austral summer, mainly over the continent (Fig. 3j). There is a good agreement between the RegCM4 simulation and the CLARA2 data, as seen in Fig. 3j, k. However, the RegCM4 values are underestimated in the central and northwest SA and parts of Argentina and Paraguay (around − 5 to − 20%), while in some parts of the northeast and south of SA, they are overestimated by + 5% to + 20% (Fig. 3l). Again, the CFC is underestimated by the RegCM4 over the oceans. Figure 4 shows the mean seasonal net SW at surface for the CLARA2 data, RegCM4 simulations and the difference between them, for the period 1992-2012. As can be seen in Fig. 4a-c, the RegCM4 is able to reproduce the spatial distribution of SW in the austral summer. As noted in Fig. 4c, northward from 20° S over the continent, the SW fluxes are lower compared to the rest of the continent because the simulated CFC (and precipitation) are higher in this season. Otherwise, over the oceans the simulated SW fluxes are overestimated in the regions where the CFC is underestimated. The areas where the simulated values are overestimated agree with the spatial pattern of the observed CFC. In some parts of the Pacific ocean, mainly in the regions near the littoral of Peru and Chile as well as some regions of the Atlantic ocean, the overestimated values show areas where SW are higher than + 40 Wm −2 . This may be due to the fact that the RegCM4 is not able to simulate adequately the lowlevel cloud cover which is dominant in the regions. Garreaud et al. (2007) showed that in the coastal regions of Peru and Chile this kind of cloud is due to the thermal inversion which occurs between 1 and 2 km height.
In addition, in these regions there is wind subsidence caused by the South Pacific Subtropical High (SPSH). As showed in Fig. 7.6 of Boucher et al. (2013), the low-level clouds are mainly located in the oceanic regions. In the austral autumn and spring, the surface net SW in the continent is underestimated in the entire region of Andes mountain, while in central SA it is overestimated (Fig. 4f, i), with values close to those observed in the austral summer (Fig. 4c). In the austral winter, the RegCM4 simulation of SW is better than in the other seasons. In most of the continent, the values of SW are slightly overestimated (Fig. 4g-i). In the Andes region, the simulated values are underestimated by around − 30 Wm −2 . Similar to the other seasons, in some parts of the ocean regions, the simulated SW values are overestimated due to the fact that the model is not able to simulate adequately the low-level clouds. Figure 5 shows the seasonal spatial distribution of the net LW at the surface (CLARA2, RegCM4 simulations and the difference between them). It can be noted, in all the seasons the RegCM4 simulations of the net LW shows lower differences compared to that observed in CLARA2 data than in the case of the simulations of the net SW (Fig. 4). In the austral summer (Fig. 5a-c), the net LW is underestimated by the RegCM4 in the central SA (around − 10 to − 20 Wm −2 ), while southward of 20° S, over the Andes mountains, the underestimated values are up to − 40 Wm −2 . Regions of underestimated values of the net LW up to − 40 Wm −2 are also seen in the eastern Pacific near the littoral of Peru and Chile. Otherwise, over some regions of the Atlantic Ocean, the net LW is overestimated, but with smaller values compared to the underestimated regions. In the austral autumn and spring, the spatial distribution of the positive and negative biases in the continent (Fig. 5f, l) is similar to that obtained in the austral summer (Fig. 5c), although the values of the biases are lower. These results are better than previous version of RegCM4 (Pessacg et al. 2014), especially in the summer and winter seasons.
Over the oceans, the spatial distribution of the positive biases in the austral autumn is similar to that noted in the austral summer, while in the austral spring the spatial pattern of the biases and the austral winter are alike (Fig. 5i). As can be seen in Fig. 5g-i, in this season, the RegCM4 shows better performance in reproducing the observed LW. However, the observed values of LW are higher (Fig. 5g) due to the lower CFC (Fig. 3a) during the dry season. In the center-west of SA the simulated values of LW are overestimated, with some areas showing values up to + 20 Wm −2 . In the Andes region, the simulated values are underestimated (around − 10 Wm −2 ). In the Pacific ocean, near the littoral of Chile and Peru, the region of overestimated values reaches + 40 Wm −2 , while in the Atlantic Ocean the underestimated values (around − 10 Wm −2 ) are predominant. Figure 6 shows the mean seasonal albedo for the observations (CLARA2), the RegCM4 model and the differences between them. As can be seen in Fig. 6c, f, i, l, the modeled albedo is overestimated in all the seasons over the eastern and southern parts of SA and oceans. Over the oceans, the overestimated values are higher (lower) than 0.25 in the austral autumn and winter (austral summer and spring). The RegCM4 performs better to represent the albedo in Fig. 4 As in Fig. 3, but for the net shortwave radiation at surface (Wm −2 ) for 1992-2012 period the continental regions, although some parts of the Andes region where the Atacama Desert and other desertic areas (Paracas, Sechura, etc.) are located, the modeled values are overestimated (higher than 0.25). A feature of the desertic regions is that the albedo varies from 0.3 to 0.45 (Rechid et al. 2009;Riihela et al. 2013). It can be seen in Fig. 6a that the CLARA2 data underestimated the albedo values in these areas. This may be due to the deficiency of the CLARA2 instruments that incorrectly calculate the heterogeneity of the near-infrared reflectance near the pixel scales of CLARA2 (Karlsson et al. 2017).
On the other hand, the RegCM4 overestimated values in these desertic areas due probably to the fact that the BATS scheme was calibrated using values of albedo of North America (Dickinson et al. 1993). It may be taken into account that the RegCM4 albedo is prescribed considering  Fig. 3, but for the net longwave radiation at surface (Wm −2 ) the type of soil cover so that in the continental region the albedo values show a better distribution by the kind of vegetation cover than in CLARA2. In the Amazon basin the RegCM4 simulates well the characteristic albedo of the rainforest during the year (values of 0.09-0.16). In the austral autumn, winter and spring, the modeled values are overestimated in the east SA (values higher than 0.15 in the austral winter). This may be due to the higher RegCM4 reflectance of the vegetation cover in this region. Figure 7 shows the temperature at 2 m above the surface (T 2m ) over SA for 1982-2012 from CRU data. RegCM4 simulations and the difference between them. T 2m is the result of the net of radiation (SW and LW). Thus, it is a very important climate variable due to its influence on the atmospheric Fig. 6 As in Fig. 3, but for the surface albedo processes. According to its latitudinal and seasonal variation, the earth's surface absorbs SW radiation for about 9-12 h. As can be seen in Fig. 7a, in the austral summer, northward from 30° S, in the regions where the altitude is lower than 750 m, both CRU and RegCM4 simulation show higher values of T 2m . However, the RegCM4 underestimated T 2m in almost the entire SA continent (values between − 0.5 and − 4.0 °C). The lower values are located in the central and south Andes. Otherwise, in the north Argentina T 2m is overestimated, with values higher than 2 °C (Fig. 7c). In this region, which are part of the La Plata basin, several models show the same bias (Reboita et al. 2014;Chou et al. 2014).
In the austral winter, the RegCM4 model shows the same deficiency to represent the spatial distribution of T 2m as in the case of the austral summer. The values of T 2m are underestimated in most of SA, with the higher values in the Andes Fig. 7 As in Fig. 3, but for 2 m air temperature (°C) Published in partnership with CECCR at King Abdulaziz University region. On the other hand, there is an overestimation of T 2m in the central SA (values between + 0.5 and + 2.0 °C). This is a common characteristic of the different versions of the RegCM4. In the austral autumn and spring, the RegCM4 simulations show a better performance compared to the austral summer and winter. However, in the north Argentina, the overestimated values (lower than + 1 °C) are noted, while in some parts of SA there are underestimated values from − 0.5 to − 4.0 °C. The underestimated values in most of SA may be associated with the higher values of albedo and CFC simulated by the RegCM4 in the different seasons, as seen in Figs. 3 and 6, respectively. Figure 8 shows the mean seasonal precipitation from CRU/ PERSIANN, RegCM4 simulations, and the difference between them for the period 1982-2012. In the austral summer, the position of the ITCZ is well simulated by the RegCM4 compared with CRU data, despite in some areas overestimated and underestimated values are noted (Fig. 8c). The model is also capable to simulate the SACZ, although its position is shifted southwards. This may be due to the fact that modeled low-level winds are stronger than the observations so that the SACZ is displaced southwards due to the more intense winds in the northeast of the system, as shown in Fig. 2. Otherwise, the SA region is strongly influenced by the adjacent oceans that modulate the different atmospheric interactions, especially those associated with the large scale (Yoon and Zeng 2010;Erfanian and Wang 2018;Wang 2019). Erfanian and Wang (2018) pointed out that the COR-DEX domain is not able to adequately represent large-scale processes (for example, teleconnections) in the SA and the adjoining oceans (Pacific and Atlantic). On the other hand, they showed that a larger domain than CORDEX allow us a better representation of the location and intensity of the rainy and dry regions in AS.

Hydrological Cycle
In the other seasons, it can be seen that the RegCM4 is capable to reproduce the northward and southward movement of the ITCZ according to the season, following the seasonal heating of the oceans in each season. However, the branch of the ITCZ located in the region of Pacific Ocean shows underestimated values of precipitation in all the seasons. This may be associated with the RegCM4 dynamics where the stronger than observed winds displace the moisture flux that favors the precipitation in this region (Fig. 2c,  f). In austral autumn, the ITCZ starts to return northwards. Consequently, the higher precipitation rates in the SA continent are located northward of 12° S. While in austral winter, as the ICTZ reaches its position northwards the high precipitation rates in SA continent are placed northward from 6° S (see Fig. 8b). The position of the precipitation maxima in SA simulated by the RegCM4 is similar to those in earlier studies (Garreaud et al. 2009;Reboita et al. 2010).
In the austral spring, the regions of SA with higher precipitation rates are placed northward of 10° S because of the displacement of the ITCZ southwards. However, the RegCM4 has a tendency to underestimate the values of the precipitation rates besides showing a larger extension of the ITCZ. This is related to the fact that in the CRU data the precipitation maxima are located in the northwest SA, while in the RegCM4 simulation they cover the entire region. In the southern La Plata basin, a higher frequency of frontal systems occurs in the austral autumn, winter, and spring ( Fig. 8d, g, j), where the austral winter is the season with highest number of cold fronts (Andrade and Cavalcanti 2004;Cardozo et al. 2015). This contributes to the high precipitation rates in the region. As shown by de Jesus et al.
(2016), the frontal system contributes with more than 84% of precipitation rates in this region. Nevertheless, the RegCM4 underestimates the values when compared to CRU data (Fig. 8f, i, l). The underestimated values exhibited in these regions are in agreement with the results of de Jesus et al. (2016), which showed that the inaccuracy of precipitation by RegCM4 varies between − 25% in summer up to − 36% in spring. This is a problem noted in the different available configurations of the RegCM4, which show deficiencies in these areas (da Rocha et al. 2012;Reboita et al. 2014;de Jesus et al. 2016).
The difference between the precipitation and the potential evapotranspiration (P − PET) is the net water flux from the atmosphere to the earth's surface (Swenson and Wahr 2006). So, it is taken as a measure of the water availability of the terrestrial hydrological cycle (Greve and Seneviratne 2015), providing a measure of the hydric stress. Particularly, in the Amazon forest the evapotranspiration plays an important role in the surface-atmosphere interaction processes since it is the mechanism by which large water amounts are transported to the atmosphere, thus contributing to the precipitation rates in the region.
As shown in Fig. 9a, in the austral summer the CRU data present negative maximum values of P − PET over the Andes mountains between 55° S-20° S, with magnitude up to − 6 mm day −1 in the southward part and with lower values in the west between 20° S-4° S. A similar distribution pattern of P − PET in the austral spring is observed. This is in agreement with Figs. 3 and Fig. 4, which show that these regions receive large amounts of SW due to the low CFC, allowing that the surface absorbs a higher energy amount to heat the several soil layers. Thus, the values of precipitation are low so that PET is dominant. As can be seen in Fig. 9b, the RegCM4 is capable to simulate the spatial pattern of P − PET in these areas, although with lower spatial extension and weaker negative values. This is due to fact that the RegCM4 simulations present higher precipitation values (Fig. 8b, c).
Moreover, in these areas T 2m is underestimated (Fig. 7c). As can be seen in Fig. 9a, in the austral summer, the north SA region also presents negative P − PET values, showing that PET is dominant due to the low CFC and high SW. In the RegCM4 simulation a similar pattern is obtained, though with higher PET values. Another region with dominant PET values is observed in the semiarid northeast Brazil. From   Fig. 9a, b, d, e, it can be noted that in the austral summer and autumn the extension of the negative P − PET is smaller in RegCM4 simulations than in CRU data, while the opposite occurs in the austral spring (Fig. 9j, k).
The continental dry season occurs in the austral winter, particularly in the central SA, where negative P − PET values are observed as a response to the low CFC and high T 2m (Figs. 3, 7, respectively). The negative values of P − PET are more intense in the east side of central SA in the RegCM4 simulations, as shown in Fig. 9c. This is because the RegCM4 overestimates the observed T 2m (Fig. 7c). In the austral winter, the distribution pattern of P − PET is similar to the other seasons. However, the P − PET values are positive between 50° S and 30° S in the west side of the Andes (Fig. 9g). This is due to the higher frequency of frontal systems reaching the region in the wintertime (Carvalho et al. 2004) which, in crossing the mountain, cause the increase of precipitation in the west side of the cordillera. Moreover, in this season the south Pacific high is located close to the continent contributing to the occurrence of precipitation. The RegCM4 simulation presents a larger (smaller) territorial and latitudinal extension of the positive (negative) values of P − PET in the region of the Andes mountains. Similar to Fig. 9 The same as in Fig. 8, but for the precipitation minus potential evapotranspiration (P-PET) (mm day −1 ) the other seasons, in the Andes region the RegCM4 shows underestimated and overestimated values of P − PET, however the RegCM4 errors may be due to SW and P biases in this area.
The observed positive P − PET values over SA follow the seasonal distribution of precipitation, where in the austral summer the ITCZ and SACZ are the principal systems that contribute to the high precipitation rates (Fig. 9a). Even though the SW is high, the high values of CFC significantly reduce the SW that reaches the surface and thus the SW values do not play a significant role in the PET processes. In these regions, the main contribution to the convective processes is due to the moisture transported by the tropical Atlantic trade winds (Drumond et al. 2008;Duran-Quesada et al. 2009;Zemp et al. 2017). In the austral autumn, the ITCZ displaces northwards so that the positive values of P − PET are located to the north of 20° S (Fig. 9d). However, the RegCM4 simulated values are higher than the observed 4 mm day −1 , as shown in Fig. 9e, f.
In the austral winter, there are negative values of P-PET in the central SA due to the movement of the ITCZ northwards and areas with positive values located to the north of 5° S (Fig. 9g). As shown in Fig. 9h, the simulation of the RegCM4 is in a good agreement with the observations, although overestimated and underestimated values can be seen. In the austral spring, as the ITCZ starts to return southwards so that the positive values of P-PET are located in the northwest and parts of the central SA (Fig. 9j), while in the RegCM4 simulation the positive values appear in the entire north and some parts of SA (Fig. 9k). Figure 10 shows the box-plot for the annual cycle of precipitation, surface net radiation (SW minus LW) and T 2m to period of analysis. As can be seen, precipitation simulated by the RegCM4 is overestimated in the regions of Amazon and northeast Brazil, mainly in the austral summer and spring months, while in the other months the simulated values are close to those in CRU data. On the other hand, it is noted that in the Amazon region the variability of precipitation simulated by the RegCM4 shows an interannual distribution similar to that in CRU data, confirming the model good performance in this region (Fig. 8). In the region of the northeast Brazil, the interannual distribution of precipitation is also well represented by the RegCM4, with little variability in July and September and higher variability in the other months. The high overestimated values (with precipitation rates similar to those in the Amazon region) is due to the fact that the RegCM4 is not able to simulate the precipitation rates in the regions with ocean-surface contrast, as can be seen in the third column of Fig. 8. This is a problem of regional models (Giorgi and Gutowski 2015). In the regions of the southeast Brazil and La Plata basin, the RegCM4 performs better the interannual variability of precipitation, showing simulated values in good agreement with those in CRU data. In the four regions, the RegCM4 shows better performance in the months from May to September. The interannual variability of surface net radiation is in a good agreement with CLARA2. Also, the biases are lower than those in Figs. 3 and 4. It can be noted that in La Plata basin region there is a larger amplitude: the maxima (around 185 Wm −2 ) and minima (around 30 Wm −2 ). In the other three regions, the surface net radiation varies between 110 and 190 Wm −2 , each region showing own maxima and minima because of their different meteorological conditions. At last, it is noted that in the regions of Amazon, northeast Brazil and southeast Brazil T 2m is underestimated along the year, even thus the RegCM4 reproduces well the interannual variability observed in CRU data. In addition, the amplitude of the interannual variability in the regions of Amazon and northeast Brazil is similar, varying between 22 and 28 °C, while in the southeast Brazil it varies between 16 and 26 °C. The better performance of the model in simulating the interannual variability is for La Plata basin region, showing amplitudes from 10 to 27 °C. The results presented for the Amazon region is similar to those in Lange et al. (2015), however the simulation in the present study shows lower biases.

Discussion
The analysis carried out on the surface radiative balance and the hydrological cycle indicate that CFC plays a very important role in both analyzed components. It is noted that the regions with higher values of CFC are higher precipitation rates areas, as well as with lower SW at the surface and LW emitted to space occur. While in the regions with lower CFC values the opposite occurs (Figs. 2, 3, 4, 8). Undoubtedly, clouds are one of the most relevant components of the atmosphere system (Stephens 2005;Henderson et al. 2013;Matus and L'Ecuyer 2017), providing one of the strong climate feedback mechanisms (Stephens 2005;Soden and Vecchi 2011). The biases in CFC during the months of austral summer and winter are lower than those in Pessacg et al. (2014) and Lange et al. (2015), with values lower than 30%. As well as, the simulation in the oceanic regions is better than that in Pessacg et al. (2014) and Lange et al. (2015), although the values are higher than 30%. Moreover, the cloud microphysics needs to be improved, mainly in the oceanic regions, where the RegCM4 shows poor performance. Regarding the microphysical mechanisms for cloud feedback, significant uncertainties still remain (Gettelman and Sherwood 2016). To reduce the biases in CFC simulation in the RegCM4 the methodologies of Lange et al. (2015) and Brisson et al. (2016) may be used. This may be subject for a forthcoming study. Pessacg et al. (2014) and Chiacchio et al. (2015) showed that the SW and LW are strongly influenced by CFC, albedo and T 2m . This is in agreement with the results presented in Figs. 2, 5 and 6, where is noted that in the regions with lower biases of these variables the SW and LW are better represented. On the other hand, the albedo directly affects the radiation flux that is absorbed by different types of surface, where the vegetation and oceans have low albedo. The RegCM4 simulations show albedo values relatively higher (lower) than the observed to the south (north) of 20° S, with higher values of + (−) 0.01. These albedo biases are mainly associated with the fact that the albedo is parameterized considering each type of surface and does not adequately represent their annual variability. The overestimated albedo in the RegCM4 causes lower values of T 2m northward of 20° S over the SA continent. This is because of the higher albedo reflects more SW, so that there is lower amount of energy to heat the surface. Also, the underestimated T 2m southward of 20° S and in other regions of SA can be associated with the BATS scheme used in the present study, which calculates the soil temperature by force-restore method ) that splits the ground into two layers to approximate the surface energy budget. Over the oceans, higher values of albedo are also noted. However, the surface energy deficit may be compensated by a high amount of SW reaching the surface due to the lower CFC. However, the errors in T 2m in regional models have lower impact (lower than 10%) than the errors in LW (Pessacg et al. (2014). Pessacg et al. (2014) and Kothe et al. (2011) showed that the errors in SW and LW are mainly due to the uncertainties in CFC. The errors in SW and LW presented in Figs. 3 and 4 are relatively lower than, or in some cases equal to, those in earlier studies over SA, Europe, and north Africa (Pessacg et al. 2014;Kothe et al. 2011;Chiacchio et al. 2015;Panitz et al. 2014). The main sources of uncertainties may be associated with the configurations of the models and the geographic influence. Many dynamical, physical, and biologic processes that occur in the surface are a direct response to the SW and LW radiative fluxes. Thus, a good representation of SW and LW fluxes in regional and global models is essential.
The hydrological cycle is an essential geophysical process that provides continuous distribution of water between the different components of the climate system. Therefore, a good representation of precipitation is important to better understand the different processes that occur in the diverse phenomena of the atmospheric system. In most of SA, it is noted that in the regions with higher values of CFC (Fig. 3) there is an occurrence of high precipitation rates (Fig. 8). Thus, a good representation of CFC, as well as the cloud microphysics processes, will improve the simulation of the precipitation rates. In addition, a good representation of the dynamic of the regional circulation is important due to the fact that it permits a good simulation of the atmospheric systems in the different parts of SA. On the other hand, the RegCM4 deficiencies in intensity and temporal evolution of precipitation, mainly in the central and southern Amazon, are associated with the selected domain (Erfanian and Wang 2018), which fails to properly represent the adjoining oceans teleconnection processes. In northern and northeast parts, the differences are associated with deficiencies of RegCM4 in representing precipitation rates. It must be taken into account that the global and regional models have yet deficiencies to represent the nonlinear relationship between SST and precipitation over tropical regions so that overestimated or underestimated modeled precipitation rates may be noted Flato et al. 2014;Zhang et al. 2015). In addition, the deficiencies in representing the precipitation in the southeastern SA (see Fig. 8) would be associated with the poor representation of the cyclogenesis by RegCM model (Reboita et al. 2018a, b). The bias of precipitation rates over the Andes region may be associated with the deficiencies of regional models to represent satisfactorily the intensity of the winds in mountain slopes due to their abrupt changes (see Fig. 2). As a consequence, the microphysics processes are affected and, in general, lead to the overestimation of precipitation in the summer months (the winds carry moisture from Amazon to the Andes). In addition, recent studies have shown that regional models have some deficiencies in representing the diurnal cycle (Bechtold et al. 2004;Lange et al. 2015;Ou et al. 2020). This is another important factor that may influence the performance of the model simulation of precipitation. Evidently, improvements in cloud representation will allow us to better depict the spatial and temporal precipitation rates in regional and global simulations.

Conclusion
In this paper, we conducted a 34-long year transient simulation to assess the ability of the RegCM4 to represent the surface radiative balance and the hydrological cycle over the CORDEX SA domain.
The RegCM4, in general, was able to reproduce the main patterns of the observed seasonal variables associated with the surface radiative balance and hydrological cycle. However, some differences were noted. The RegCM4 overestimated the variables associated with the radiative balance (CFC, SW, LW and albedo), mainly over the oceans, in all the seasons. The net surface SW and LW errors were associated to the errors in the CFC due to the inaccuracy of the model in representing the lowlevel clouds, which is predominant in the regions of the Pacific Ocean near to Peru and Chile as well as most part of the south Atlantic ocean (Boucher et al. 2013). It must be pointed out that the CFC plays an important role in modulating the net surface radiation. Besides, the albedo was not well simulated by the model which contributes to the errors in SW and LW.
The RegCM4 simulations of the hydrological cycle variables also presented differences compared with the observations. Although the position of the ITCZ in the austral summer was well simulated by the RegCM4, in some areas, overestimated and underestimated values were noted. In the other seasons, the RegCM4 was capable of reproducing the northward and southward movement of the ITCZ, according to the seasonal heating of the oceans. However, over the Pacific Ocean the seasonal warming was underestimated, which may be associated with the RegCM4 dynamics where stronger than observed winds shift the position of the moisture flux that sustain the precipitation in this region. Some deficiencies in the model simulations may be associated with the CORDEX domain, which fails to adequately represent the oceanic processes that influence the SA climate (Erfanian and Wang, 2018).
It is clear that simulation of surface temperature has to be improved to better represent evapotranspiration, which is an important component of the hydrological cycle since it transfers soil moisture to the atmosphere. In addition, the cloud cover fraction also needs to be improved to reduce the biases in precipitation, PET, and the radiative surface fluxes SW and LW. To reduce the deficiencies in the model simulations, new physical schemes have been implemented in RegCM4. For example, the scheme of the University of Washington which gives a better simulation of the low cloudiness and the complex microphysics scheme (hydrometeors) (Nogherotto et al. 2016) is recently implemented into RegCM4. Also, the Community Land Model (CLM) surface scheme may be used to obtain a better simulation of temperature and surface radiative balance, as shown by Reboita et al. (2014). The use of these schemes may be the subject of forthcoming study.