WRF Sensitivity for Seasonal Climate Simulations of Precipitation Fields on the CORDEX South America Domain

: Dynamic numerical models of the atmosphere are the main tools used for weather and climate forecasting as well as climate projections. Thus, this work evaluated the systematic errors and areas with large uncertainties in precipitation over the South American continent (SAC) based on regional climate simulations with the weather research and forecasting (WRF) model. Ten simulations using different convective, radiation, and microphysical schemes, and an ensemble mean among them, were performed with a resolution of 50 km, covering the CORDEX-South America domain. First, the seasonal precipitation variability and its differences were discussed. Then, its annual cycle was investigated through nine sub-domains on the SAC (AMZN, AMZS, NEBN, NEBS, SE, SURU, CHAC, PEQU, and TOTL). The Taylor Diagrams were used to assess the sensitivity of the model to different parameterizations and its ability to reproduce the simulated precipitation patterns. The results showed that the WRF simulations were better than the ERA-interim (ERAI) reanalysis when compared to the TRMM, showing the added value of dynamic downscaling. For all sub-domains the best result was obtained with the ensemble compared to the satellite TRMM. The largest errors were observed in the SURU and CHAC regions, and with the greatest dispersion of members during the rainy season. On the other hand, the best results were found in the AMZS, NEBS, and TOTL regions. This study evaluated the systematic errors and areas with large uncertainties in a simulated precipitation over the South American continent. Different convective, radiation, and microphysical schemes of the WRF model were tested. The analyses through Taylor diagrams showed that, in all sub-domains, the ensemble of all experiments was the best option, with results consistent with those from the satellite data obtained from TRMM. The results also showed that the WRF simulations were better than those of the ERAI reanalysis, which was used as initial and boundary conditions for the WRF model, showing the added value of dynamic downscaling.


Introduction
The South America continent (SAC) spreads over a large latitudinal extension, ranging from equatorial areas to mid-latitudes, and it is influenced by several weather and climate patterns. The South American climate and its variability are affected by remote, regional, and local forcings [1]. Understanding the mechanisms acting on these three scales of variability, in particular how they impact the hydrological cycle, is critical for medium range and seasonal precipitation forecasts. In turn, this is important for water resources management, hydro-power generation planning, and agricultural activities in SAC [2].
The use of dynamic numerical models is the basis of weather and climate forecasts around the world. The general circulation models (GCMs) have shown good ability at pre-dicting and representing large-scale phenomena, also being employed for studies of future climate change scenarios [3]. However, the low spatial resolution hinders the accurate representation of mesoscale processes [3]. To improve the skill of global simulations, regional climate models (RCMs) with higher resolutions have been used for climate forecasting [4]. The higher accuracy of RCM simulations is partially due to better representation of the topography, coastlines, vegetation, and physical process parameterizations of meso and micro scales, which allow better representations of local influence factors on the climate of a particular region. For that reason, RCMs have become widely used to downscale low-resolution global climate simulations throughout different regions of the world. They provide a clear added value with respect to simulations with coarse resolution GCMs, particularly for variables, such as precipitation and near surface temperature [3,5,6].
COordinated Regional Climate Downscaling EXperiment (CORDEX)) [7] is a coordinated international effort whose goal is to improve the quality of downscaling in regional climate simulations with RCMs. Among other characteristics, CORDEX standardized reference areas for numerical studies. This strategy facilitates and adds quality to the intercomparison between simulation studies from different research groups. Within the several subgroups of CORDEX, there is the CORDEX-WRF component (http: //www.meteo.unican.es/wiki/cordexwrf/CordexWrfHome, accessed on 18 August 2021) that concentrates on evaluating the weather research and forecasting model (WRF) [8]. In this way, the simulations carried out in this work are inserted in the CORDEX-WRF context, with special focus to the modeling of precipitation.
In recent decades, RCMs have been broadly used to produce climate change projections over the SAC region in an attempt to better capture regional and local feedback processes, e.g., the CLARIS project [9]. In particular, these studies have evaluated different RCMs to understand present day climate and future climate scenarios [10][11][12][13][14][15][16][17][18][19][20]. For a review of regional climate modeling studies over SAC, reference [1] summarized the main progress made since the early efforts at the beginning of the 2000s. They report that the earlier research aimed to test the ability of a diversity of RCMs (e.g., ETA, MM5, RegCM, among others) in reproducing basic aspects of seasonal climate, while from 2002 to 2007 the focus was on evaluating model performance and sensitivity. Since 2008 most of the efforts have been focused on the production of regional climate change scenarios. Reference [21] have made a compilation of the main studies using RCMs for South America.
In the current work, we assess the ability of the WRF model to reproduce the observed precipitation variability in SAC under the CORDEX-SOUTH AMERICA framework. Our focus is on the investigation of the optimal combination between radiation, convection, and microphysics schemes, and on the detection of systematic errors and areas with large uncertainties. The radiation parameterization schemes aim to estimate the total radiative flux and provide the shortwave heating and long wave cooling rates in the cloud layer. Radiation plays a significant role in the surface water-energy budget at the land surface. The total radiative flux is essential to calculate net radiation, sensible heat flux, and latent heat flux at the land surface. These fluxes also depend on the soil and vegetation types. Recently, the impact of radiation on these fluxes over the Amazon basin has been discussed by [22]. The microphysics scheme handles the phase changes of water, including cloud and precipitation processes. Finally, the convection parameterization schemes lead with the subgrid-scale representation of physical processes related to a wide range of cloud types, such as shallow, deep, and stratiform clouds at a horizontal resolution more major than cloud scale (<3 km). Convection is one of the key processes in atmospheric models, considering its important role in the atmosphere's regional and global water and energy cycles by transporting heat, momentum, and moisture from the boundary layer to the free atmosphere. Moreover, convection is a significant driver of regional and global atmospheric circulations through the release of latent heat, as discussed by [23] on the South American summer upper-level circulation. In particular, in the tropics, where cumulus convection is extremely active, large amounts of precipitation and accompanied latent heat release play a key role in tropical atmospheric circulation and global energy balance [24]. As a consequence of this chain of mechanisms, different choices in radiation/convection/microphysics yield very different results in the simulated precipitation fields.
This study was conducted in a climate perspective (a 2-year simulation) instead of a weather forecast approach, which is the most frequent use of WRF. Despite this period being short compared to the usual climate timescale of 30-years, it does not prejudice the main goal of this study i.e., to find out which set of parameterizations is better suited to South America. In fact, longer simulations in the SAC using WRF are rare [25,26]. A similar study was made by [27], who applied the dynamical downscaling technique in the WRF model over the MED-CORDEX domain for a period of one year (2002), to investigate and validate the performance of different physics parameterizations. This study fills this gap and contributes toward the CORDEX effort over this region for subsequent climate modeling in future experiments. In addition, many users of the WRF model carry out research on different regions of the SAC and, therefore, our results can serve as a comparative reference.

Model Configuration and Numerical Experiments
Initial and lateral boundary conditions used to constrain the WRF model simulations are from the European Centre for Medium-Range Weather Forecasting (ECMWF) interim reanalysis (ERAI) [28]. The ERAI is the ECMWF global atmospheric reanalysis using the 4D-VAR data assimilation system, and it is a continuation of the previous ERA-40 atmospheric reanalysis dataset [28,29]. These data were available on a 0.75°× 0.75°grid from 1979 to 2019 at 6 h resolution for the 3D prognostic variables, which is a spatial resolution that agrees with our 50 km simulation. Despite the already available ERA-5 reanalysis (0.25°× 0.25°), its higher resolution is not suitable for driving 50-km regional model simulations, which several climate forecast centers in SAC still do, due to computer resource limitations.
Here, we use data from 2001 and 2002. The choice of these years had no special reason and they were chosen randomly from the dataset. In order to account for recent changes in land use (e.g., deforestation), updated datasets of land-surface albedo and greenness fraction at a high resolution were used instead of the standard WRF datasets (details in [30]).
Numerical simulations were conducted using the Advanced Research Weather Research and Forecasting model (WRF-ARW) version 3.9.1 [8]. The total simulated time is 16 months. It started at 00 UTC on 1 September, 2001, and finished on 21 UTC on 31 December 2002. The first four months of outputs were considered as the spin-up time, and were excluded from the final analysis. This is a more conservative approach than recommended by [4], and gives the model more time to equalize the surface hydrology with the external forcing, avoiding instabilities or drifts.
The simulations cover the CORDEX-South America domain ( Figure 1) with a horizontal resolution of 50 km and 38 vertical levels with the top of the model at 50 hPa. The domain has 150 by 200 west-east and north-south grid points, respectively, in the Mercator projection.
To examine the WRF model uncertainties in simulating the annual precipitation regime over South America, a sensitivity test was performed where we combined different physical parameterizations, as detailed in Figure 2. There are many different physics options in the WRF model for cumulus convection, turbulence, radiative transfer, gravity waves, microphysics, planetary boundary layer (PBL), land surface, etc. In the present study, we focus on cumulus parameterization (CU), radiation parameterization (RA) and microphysics (MP), while PBL (Mellor-Yamada-Janjic; MYJ) [31,32], surface layer (Monin-Obukhov Janjic; MOJ) [33,34] and land surface (Noah LSM) [35] schemes were kept fixed for all model configurations.
A total of 10 simulations were run, corresponding to the 10 different model configurations (see Figure 2). The ensemble mean was also evaluated. The ensemble mean experiment was generated from the mean of ten individual members. Three convection schemes were used, namely the Kain-Fritsch (KF) [36,37], Grell-Devenyi (GR) [38], and Betts-Miller-Janjic (BM) [39]. Three radiation schemes were used as well, the newer version of the rapid radiative transfer model (RG) [40], Dudhia scheme (DH) [41], and CAM scheme (C3) [42]. As for the microphysics, we evaluated the WRF single-moment six-class scheme (S6) [43] and WRF double-moment six-class scheme (D6) [44]. We considered the control simulation (CTRL) as the combination of KF cumulus scheme, RG radiation scheme, and S6 scheme for the cloud microphysics, which is the mostly employed model setup in operational forecasting in the SAC domain. For each model configuration, there is one or more model configurations that differ in just one of the parameterizations, allowing to evaluate the effect of changing each one in isolation.

Model Evaluation
The evaluation of the different model configurations was conducted in two ways. First, Tropical Rainfall Measuring Mission (TRMM) [45] version 3B42 V7 data were used as the "observational truth". TRMM has been widely used over South America [46][47][48], and recent studies [49] have shown a satisfactory agreement with ground observations by meteorological stations.
The evaluation of monthly precipitation variability and model systematic errors were performed for each of the CORDEX sub-domains shown in Figure 1  The simulated monthly average precipitation was assessed using Taylor diagrams [50]. These diagrams summarize the degree of similarity between each WRF experiment and the reference data with respect to the standard deviation (SD), root-mean-square error (RMSE), and Pearson correlation coefficient. These diagrams are especially useful in assessing multiple aspects of complex models or in gauging the relative skills of many different models simultaneously. In addition to Taylor diagrams, we used box-plot diagrams to evaluate whether errors were random or systematic in the annual cycle. For some of the analysis, we grouped the sub-domains in three representative regions: tropical (AMZN, AMZS, NEBN, NEBS, and PEQU; from 15°S to 5°N), subtropical (CHAC and SUDE; from 23°S to 15°S) and extratropical (SURU; from 35°S to 23°S). Figures 3-6 show the average daily precipitation field for the DJF/MAM/JJA/SON quarters, respectively. The first line of plots contains the observation reference data by TRMM; the original input data to the simulations, ERAI reanalysis; the mean of ten individual members (ENSE = ensemble); and the control simulation (CTRL). These four plots are hereafter referred to as "reference cases", and the individual simulation in the three remaining lines as "individual members", or simply "members".

Seasonal Analysis
In SAC, most of the precipitation volume falls on DJF, which is shown in Figure 3. This wet period, in some specific areas, extends to March. This is the case for the majority of the Amazon basin, southeast and northeast Brazil, Peru, and Central South America. Northern portions of Amazonia, which is partially in the northern hemisphere, have their wet season shifted forward to April and May. On the other hand, the winter months of JJA ( Figure 5) are the driest in the majority of the continent (except for the north hemisphere portion of SAC, subject to inter-tropical convergence zone (ITCZ) influence, as can be seen in Figure 3).  With respect to the wet trimester DJF (Figure 3), it is possible to see that TRMM and the reference cases agree, with respect to the spatial pattern of precipitation, but clearly disagree with respect to their precipitation volumes. While CTRL overestimates the observed values, ERAI underestimates it, and the ENSE of simulations is closer to TRMM than both CTRL and ERAI. This improvement of ENSE is an important result, and demonstrates that an ensemble based on physical parameterizations, plus the finer resolution was capable of providing more accurate results than the reanalysis fields from ERAI. Looking onto specific members (lines 2-4 of plots), it is possible to see that the general spatial pattern of precipitation is correct for most of them, except for GRCS and KFCD that largely underestimated rain in the center of SAC, and KFDS/KFRS, whose results went in the opposite direction i.e., largely overestimated precipitation. The results for the Kain-Fritsch parameterization are surprising, and they suggest that it is the combination of radiation and microphysics that determines precipitation fields when KF is the chosen cumulus parameterization. Given that microphysics at the resolution of 50 km is a weak forcing, and that convective storms dominate the scene in the wet season in SAC, it seems that the radiation choice is the critical parameter when KF is employed. Figure 5. As in Figure 3, but for June-July-August quarter. The following quarter MMA is shown in Figure 4. For the rainy equatorial latitudes, the same bias of TRMM, ERAI, ENSE, and CTRL was observed. Important differences can be seen in rain volume in South of Brazil, Uruguay, and Argentina. For these regions, the CTRL setting seems to better represent rain, while ENSE and ERAI underestimate precipitation. Dry regions (eastern portion of Brazil, central Argentina) were kept dry for the totality of members, and all simulations were capable to produce rain in the southern portion of Chile with little difference among them.
The dry season in SAC, the quarter JJA, was well reproduced in a broad sense by all members, CTRL, ENSE, and ERAI ( Figure 5. Some differences remain for the SURU region, where an underestimation of precipitation was observed in all cases (when directly compared to TRMM).
The quarter SON is the dry-to-wet transition in most of the SAC. As can be seen in Figure 6, important spatial differences were obtained for the reference cases ENSE, CTRL, and ERAI, as well for individual members. These differences (anomaly) also comprise the precipitation volumes. These differences can be better discussed by the analysis of the anomalies between reference cases and individual members to the observational reference of TRMM data. These anomalies are shown in Figures 7-10.
A common feature for all months, for both reference and member cases, is an underestimated precipitation in the SURU area. Apparently, there is an intrinsic difficulty to correct simulate precipitation in this region, which is partly due to its difficulty in simulating cold fronts, and which seems to be a known phenomenon already observed in previous studies [51][52][53][54].
With respect to the anomalies in DJF, one can see in Figure 7 that ERAI mostly underestimate precipitation in the SAC while ENSE anomalies are more evenly distributed throughout the continent. CTRL simulations presented more positive anomalies inside the continental area, except for a few areas in the extreme northeast coast of Brazil and north of Argentina. For individual members, all on them in the left column underestimated precipitation over continental area. In common, they employed the C3 radiation scheme. This spatially broad negative anomaly was observed for all convection parameterizations, which is a clear demonstration that C3 induces a negative bias during wet season. Looking into the other six members, a general result is the already mentioned underestimation in SURU area.
In the fall (MAM, Figure 8), among the reference cases, ENSE presented the best result with a smooth distribution of anomalies though the three reference cases except the underestimation in the SURU area. In the Equator region, CTRL overestimated precipitation in the south of the Amazon River, while ERAI showed the opposite sign but in the north side of the same Amazon River course. With respect to individual members, those using the C3 scheme (left column) still underestimate precipitation. Member with KF cumulus parameterization greatly differed again, showing the strong dependence on the radiation scheme.
For JJA (the driest quarter, Figure 9), there is a good agreement for the continental area. In fact, most positive anomalies are over the ocean, following the climatological position of ITCZ. The few areas with anomalies in the continent are over the west and central portions of Amazonia and SURU domain, with ENSE showing the best results.
Finally, in the dry-to-wet transition quarter of SON (Figure 10), it is noteworthy the ERAI and ENSE low deviations in continental areas (except for the underestimation in SURU). However, ENSE produced significant positive deviations in the north of the Atlantic ocean (ITCZ area). Either way, both got results far superior to the CTRL experiment, which resulted in significant portions of the Amazon basin with negative deviations. Among the individual members, the KFCD experiment presented a large area with negative deviation in the continental portion of the domain. Interestingly, this pattern was not followed by the other C3 schemes (on the left column-as observed in the DJF quarter; Figure 7). A possible explanation for this behavior is that the convection schemes in BMCS and GRCS were partially able to counterpart the bias induced by C3 in SON, but not in DJF. Further investigation is necessary to corroborate (or not) this hypothesis.

Evaluation of Sub-Domain Annual Cycle
In this section, the general behavior of the simulation ensemble performed in this study is evaluated with the observational results of the TRMM satellite and the reanalysis for each subdomain. This part of the analysis aims to verify the level of agreement and coherence between these datasets in a broader approach before discriminating by a numerical experiment. Figure 11 shows for each subdomain the time series of simulated precipitation (boxplots based on the 10 members) and the corresponding TRMM satellite monthly totals.

Amazon Basin (AMZN + AMZS)
For both sub-domains corresponding to the Amazon basin (Figure 11a,b) there is a good general agreement between simulations and TRMM, though a few differences are noteworthy to highlight. At the AMZN subdomain, during the rainy months (January-April), there is a negative anomaly for February that is compensated in April, with a positive anomaly of similar magnitude. This difference is not observed for AMZS, except for a small difference in January. In spite of this small differences, ENSE x TRMM has excellent agreement in the annual cycle.

Northeast of Brazil (NEBN + NEBS)
In general, ENSE and TRMM agree very well on the annual precipitation cycle in both sub-domains that represent the northeast region (NEBN and NEBS). The biggest discrepancy was for the NEBN region in the months of March/April, which is more influenced by ITCZ when compared to the rest of the northeast region as a whole (especially in the coastal region).
The northeast of Brazil is influenced by several meteorological systems, whose intensity varies according to the portion under analysis: ITCZ, cold/stationary fronts, land and sea breeze [55], and easterly wave disturbances [56,57]. The relative influence of each system varies significantly within the NEBN subdomain itself. Thus, the results of this work show that, on average, the models do not capture this relative modulation of meteorological components in a satisfactory way especially for the wettest months (March/April), and more detailed studies are necessary to understand the reason. Here, we hypothesize that the dispersion between ENSE/TRMM is probably due to spatial differences within the NEBN and NEBS boxes, such as between the coastal region and the semi-arid center. Another aspect is that, for March, the interquartile range (IQR) of the ENSE/NEBN is very low, which shows that the choice of parameterizations is indifferent with respect to making the model agree with the TRMM observation. Differently, IQR for April is significantly larger than March, which shows the difficulty in adjusting the real precipitation. Again, the likely explanation is a possible inhomogeneity in precipitation field within the NEBN domain, which could generate bias in the TRMM observational mean. Anyway, in an integrated bimonthly perspective (March/April), the ENSE approximately hits the total rainfall. A better accuracy in this type of forecast would be of great value, for example, for local agriculture, given the low total rainfall of the NE region as a whole.

South of Brazil and Uruguay (SURU)
As already discussed in Figures 3-10 (seasonal maps), both reference cases and individual members underestimated the precipitation observed by the TRMM satellite. Figure 11f shows that the monthly rainfall curve is systematically above all boxplots. In fact, most of the current operational regional and global atmospheric models used in numerical weather and climate prediction and CMIP models still show serious deficiencies in simulating the precipitation during the wet season over South America, and the La Plata Basin. In general, models tend to underestimate rainfall over these regions as shown in previous studies [51][52][53][54]. In fact, it is a region where models perform with low skill, mostly because the success of the simulation depends on a precise determination of a low level jet, cold fronts passage, and the generation of mesoscale convective systems (quite common in this region). Even small changes in the position of these meteorological systems induce huge changes in the precipitation fields. This is the case with our simulations in the SURU sub-domain, which is under the influence of all these mechanisms.
However, in spite of the systematic errors, one can see that the annual cycle is captured in a general way, with annual minimums in the winter months (JJA), and maximums in March/April and November/December.

Peru/Equator (PEQU), Chaco (CHAC), and Southeast of Brazil (SUDE)
For this sub-domain, the annual cycle was captured in an excellent way, in spite of some dispersion of members in the rainy season. For most of the year, the observed data are very close to the medians. At CHAC, simulations caused an earlier onset of the rainy season, resulting in maximum precipitation in January while observational maxima was in February. It is noteworthy to mention that the agreement in the PEQU is surprisingly good, especially considering that this is an area of rugged relief due to the Andes, and a 50 × 50 km resolution usually fails to capture the influence of topography on the intense orographic precipitation associated.

Average over All Domains (TOTL)
The last plot shows the area-weighted daily average precipitation for all simulated domains, and the corresponding TRMM observations. The agreement between simulations and observations is excellent, and the correct reproduction of the annual cycle is remarkable, except for March, when simulations underestimated the actual precipitation. Another point to highlight is the very low variability of simulations during the dry months. In fact, the interquartile Range for July-October is very low, which means that models not only agree with observations, but they also agree with each other.

Best Sub-Domain Settings-Taylor Diagrams Analysis
To assess the sensitivity of the model to different parameterizations and its ability to reproduce the simulated precipitation patterns, an analysis was performed using Taylor diagrams. These kinds of diagrams synthesize, in a single plot, the relevant parameters for choosing the best model from a statistical point of view. Thus, it was investigated how well each member (i.e., its precipitation field) agrees with the TRMM observations in terms of root mean square error (RMSE), correlation coefficient (R), and standard deviation (STD). Note that the model's ability to simulate the reference precipitation, as well as the impact of changing the parameterization schemes, resulted in different behaviors for each sub-domain.

Amazonia
In the AMZN region (Figure 12a), the ENSE obtained the most consistent result with respect to the reference data, with R = 0.91, RMSE = 0.4, and the same STD as the TRMM data. However, as the computational costs to generate a nine-member ENSE is prohibitive in most operational centers in SAC, we can attribute the BMRS experiment as the best alternative. The GRDS simulation obtained the same STD as the reference data, but with lower R and higher RMSE. With respect to the cumulus parameterization, BM presented the best and KF the worst result. Regarding the radiation scheme, a negative impact is evident when using option C3. Figure 12b shows that in the AMZS region, the model obtained more consistent agreement with TRMM data than observed in the AMZN region. We attribute this difference to the fact that the precipitation AMZN is modulated by the ITCZ [58], whose position in March 2002 was further north of its climatological position, and in April, further south (data not shown), which may have resulted in the model's difficulties in predicting such anomalies. In the AMZS region, the GR cumulus parameterization obtained the best results, unlike AMZN, where the GRDS and GRRS simulations presented the same STD as the TRMM. However, GRDS had slightly higher R and lower RMSE. The cumulus KF parameterization showed the worst results, with high RMSE and STD.
In both regions, the use of the D6 microphysics had a negative impact with the cumulus KF parameterization, with a strong reduction in the correlation with the reference data. The KF scheme uses low-level vertical movements as a trigger function and the removal of convective potential energy available as a closure. In fact, according to [59], in the Amazon region, the available convective potential energy (CAPE) varies with the intensity of convection, but does not explain convection variations, suggesting that a convective parameterization based on CAPE does not capture the observed behavior.

Northeast of Brazil
In the NEBN region (Figure 12c), the model's sensitivity to parameterization changes is small when compared to the AMZN and AMZS regions. The GR and BM parameterizations obtained very similar results, except for the BMRS experiment, with the worst performance for this region. GR had higher correlations compared to the TRMM and lower RMSE values, with the GRRS option being slightly more consistent. In general, the model was not very sensitive to changes in convective parameterization in this region, except when using the RRTMG scheme, which had a negative impact when combined with the BM and KF schemes. The best performance was obtained with ENSE, the same occurring in the NEBS region (Figure 12d). In this region, the KFDS and KFRD configurations presented results slightly more consistent with what was observed.

Andes
The Peru-Ecuador region (Figure 12e) presents a strong sensitivity of the model to the choice of parameterization schemes used. It is evident that the Grell parameterization obtained the best results, regardless of the radiation scheme, which in this case was insignificant. It is noteworthy that the GRDS and GRRS experiments had the same STD as the reference data, but GRDS was slightly superior and had lower RMSE. The BM and KF experiments showed low correlations with the TRMM (except the CTRL experiment), and high RMSE values. Again, ENSE obtained a result consistent with the reference data.

South of Brazilian and Uruguay
In the SURU region (Figure 12f), the GRRS experiment (KFCD) presented the result more consistent with what was observed. Similar to what happened in the NEBN region, this region is not very sensitive to changes in the parameterization schemes. A highlight of this region is that the use of the C3 radiation scheme had a negative impact, regardless of the cumulus parameterization.

Southeast Brazil
The best result for the SUDE region (Figure 12g) was obtained with the GRDS option, while the cumulus KF parameterization scheme presented the worst results, regardless of the radiation scheme used. This result is in agreement with the study by [60], who evaluated the ability of the MM5 model to simulate intraseasonal variability during the warm season over South America through sensitivity experiments, and their results showed that Grell's cumulus scheme was most suitable for subtropical latitudes.

Chaco
In the Chaco region (CHAC), as in other sub-domains, the ENSE option obtained the most consistent result compared to TRMM (Figure 12h). Individually, the BMCS experiment showed the highest correlation with the reference data and the lowest RMSE, being the most recommended. Regarding the experiments with the GR cumulus scheme, these showed low sensitivity when the radiation scheme was changed.

All Sub-Domains (TOTL)
Considering the total area ( Figure 12i) the most consistent results, compared to TRMM are the ENSE and the GRDS experiment, while the largest errors are from the experiments using KF, as expected, in view of the results presented and discussed in the individual sub-domains.
The microphysics scheme was changed only using KF. It was possible to observe that the impact is negative when using option D6 in the PEQU and CHAC regions, and slightly less significant in the NEBS, SURU, SUDE, and TOTL regions. In other areas, there was no impact.

Summary and Conclusions
This study evaluated the systematic errors and areas with large uncertainties in a simulated precipitation over the South American continent. Different convective, radiation, and microphysical schemes of the WRF model were tested. The analyses through Taylor diagrams showed that, in all sub-domains, the ensemble of all experiments was the best option, with results consistent with those from the satellite data obtained from TRMM. The results also showed that the WRF simulations were better than those of the ERAI reanalysis, which was used as initial and boundary conditions for the WRF model, showing the added value of dynamic downscaling.
Regarding the choice of the convective parameterization scheme, the GR parameterization was the best option in most areas (AMZS, NEBN, PEQU, SURU, and SUDE), and satisfactory in other areas. On the other hand, the worst results were obtained with KF, with greater errors and lower correlations. In terms of radiative parameterization scheme, the impact in the precipitation distribution and seasonality was less than the convective parameterization, as expected. The DH and RG options performed better, as well as the S6 microphysical scheme. The GRDS configuration turned out to be the best configuration for the total domain and KFRD the worst.
The largest errors were in the SURU and CHAC regions, and in some experiments in the PEQU region. The evaluation of the annual cycle showed that the model was able to reproduce the variations in precipitation in most regions, highlighting the SURU region, which is a region in which precipitation occurs throughout the year. This region stood out for also presenting greater systematic discrepancies. In addition to these facts, it was also noted that the greatest dispersion of members was concentrated in the rainy periods.
This study comes to fill this gap and to contribute to the CORDEX effort over the South American region for future climate model experiments and downscaling efforts. In addition, many users of the WRF model carry out research on different regions of the SAC and, therefore, our results can serve as a comparative reference. The important conclusion that can be drawn from this work is that WRF simulations accurately reproduces features at several time scales over the SAC and, additionally, it provides information at spatial scales not resolved by GCMs that can be extremely useful in the studies of climate change scenarios.
As a future activity, we intend to apply the Brazilian Global Atmospheric Model (BAM), which is the atmospheric module of the Brazilian earth system model (BESM), as initial and boundary conditions on the WRF model, aiming to predict the sub-seasonal to seasonal and climate variability on specific regions of South America and Brazil. Moreover, based on the new implementation of aerosol-cloud microphysical interactions and cloud processing of aerosols in BAM, we will investigate the effects of aerosols on seasonal rainfall in South America at the regional scale.

Data Availability Statement:
The datasets used in this study as initial and boundary conditions for the WRF model simulations and validation were obtained from the European Centre for Medium-Range Weather Forecasting (ECMWF) and Tropical Rainfall Measuring Mission (TRMM) version 3B42 V7, respectively. The links where to download the data were provided through https:// apps.ecmwf.int/datasets/data/interim-full-daily/levtype=sfc/ (accessed on 2 December 2021) and https://disc.gsfc.nasa.gov/datasets/TRMM_3B42_7/summary (accessed on 2 December 2021) to ECMWF and TRMM, respectively.