Sensitivity Study on the Influence of Parameterization Schemes in WRF _ ARW Model on Short-and Medium-Range Precipitation Forecasts in the Central Andes of Peru

A sensitivity study of the performance of the Weather Research and Forecasting regional model (WRF, version 3.7) to the use of different microphysics, cumulus, and boundary layer parameterizations for shortand medium-term precipitation forecast is conducted in the Central Andes of Peru. Lin-Purdué, *ompson, and Morrison microphysics schemes were tested, as well as the Grell–Freitas, Grell 3d, and Betts–Miller–Janjic cumulus parameterizations. *e tested boundary layer schemes were the Yonsei University and Mellor–Yamada–Janjic. A control configuration was defined, using the *ompson, Grell–Freitas, and Yonsei University schemes, and a set of numerical experiments is made, using different combinations of parameterizations. Data from 19 local meteorological stations and regional and global gridded were used for verification. It was concluded that all the configurations overestimate precipitation, but the one using theMorrison microphysical scheme had the best performance, based on the indicators of bias (B) and root mean square error (RMSE). It is recommended not to use the Betts–Miller–Janjic scheme in this region for low resolution domains. Categorical forecast verification of the occurrence of rainfall as a binary variable showed detection rates higher than 85%. According to this criterion, the best performing configuration was the combination of Betts–Miller–Janjic and Morrison. Spatial verification showed that, even if all the configurations overestimated precipitation in some degree, spatial patterns of rainfall match the TRMM and PISCO rainfall data. Morrison’s microphysics scheme shows the best results, and consequently, this configuration is recommended for shortand medium-term rainfall forecasting tasks in the Central Andes of Peru and particularly in the Mantaro basin. *e results of a special sensitivity experiment showed that the activation or not of cumulus parametrization for the domain of 3 km resolution is not relevant for the precipitation forecast in the study region.


Introduction
Mesoscale meteorological models are a powerful tool, both for operational simulation and for atmospheric investigations [1].In this sense, the Weather Research and Forecasting (WRF) model [2] is currently one of the most used in the world for these purposes, since it supports the use of very high-resolution grids for domains in any region of the planet and allows changes in its physics schemes configuration to tune it for regional conditions without having to compile it every time.Its development began in the second half of the 1990s as a collaborative partnership, mainly between the National Center for Atmospheric Research (NCAR), the National Oceanic and Atmospheric Administration (NOAA), the Air Force Time Agency (AFWA), the Naval Research Laboratory, the University of Oklahoma, and the Federal Aviation Administration (FAA), all of the United States of America.
In the present work, a research is carried out to evaluate how the WRF model simulates rainfall in the Peruvian Central Andes.For validation, observational data from 19 stations in the region are used, as well as gridded observational data.Most of the rainfall in Peru is concentrated in the period between the months of September and April [3], defining a marked seasonality, with a dry season between May and August [4][5][6].In this region, rainfall plays an important economic role, since 71% of the arable land in the basin depends on it for crops [7].
Regional weather numerical models, including WRF, have been previously used to carry out studies on different mountainous regions, including part of the Andes [8][9][10][11].However, in Peruvian Central Andes, numerical modeling studies are scarce.In a recent paper [12], the influence of orography on the diurnal cycle of rainfall in the Central Andes has been investigated using WRF.However, there are no studies evaluating the sensitivity of the model to changes in its configuration; for example, the use of certain parameterization schemes of physical processes or the characteristics of the domains, which is of great interest due to the orographic complexity of the region, may be relevant to the forecasting skill of the model for operational applications.In the present investigation, WRF sensitivity to changes in its cumulus parameterization, microphysics, and boundary layer schemes is evaluated, relative to the results of its short-and medium-term precipitation forecast in the study region.
In the specialized bibliography, as well as in the user manual, there are dissimilar parametrization schemes of the different physical processes that WRF solves, among them those of boundary layer, clusters, and microphysics, which would be the objects of verification of this work.us, among others, for the boundary layer, there are the Mellor-Yamada-Janjic scheme (MYJ) [13] and the Yonsei University scheme (YSU) [14].For convection, the parameterizations of Betts-Miller-Janjic (BMJ) [15], Grell 3D (GRELL3), and Grell-Freitas [16], were used, and for microphysics, the parameterizations of ompson [17], Morrison [18], and Lin-Purdué [19] were used.
High-resolution numerical modeling has specific issues to deal with, as the choice to include or not cumulus parameterization schemes (CPS), and explicit microphysics.
e WRF user manual (v.3.7) suggests the explicit solution of convection for grids with these resolutions; however, in the scientific literature, there does not seem to be a total consensus on this aspect.For example, in [20], it is proposed that an explicit treatment of convection for a resolution, for example, of 4 km, shows a more accurate description of the physics of convective systems but not necessarily provides a better forecast.us, in [21], although it is stated that "Idealized modeling demonstrates that, at a fine resolution such as 2 or 4 km, classic super-cellular convective features including a reflectivity hook, midlevel rotation, and storm splitting can be resolved explicitly, and CPS is not needed," but it is concluded that it is not appropriate to assume that a simulation that uses a grid of small spacing (<5 km) will not need a CPS. is will depend on the synoptic forcing and the time of the year.
Consequently, and bearing in mind that the objective of this research is to verify the sensitivity of WRF to the use of different parametrization schemes for punctual precipitation forecast in this complex region, it was decided to include the convection parameterization schemes also for the domain of 3 km resolution.However, a special numerical experiment was made, using the best performing configuration, deactivating the cumulus parameterization for the domain of 3 km resolution to investigate the sensitivity of the forecast to this procedure.

Initial and Border Conditions.
e experiments were carried out using model WRF v. 3.7 (Weather Research and Forecasting, version 3.7).e initial and boundary conditions were taken from the "Global Operational Analysis" of the National Center of Environmental Prediction (NCEP), final analysis FNL (https://rda.ucar.edu/datasets/ds083.2/), every 6 hours, with horizontal resolution of 1 °× 1 °. e analysis contains surface data and 26 higher air levels and includes the following surface variables: sea level pressure, higher air variables (geopotential height, air temperature, and zonal and meridional components of the wind and vertical air velocity), sea surface temperature, and soil parameters.Even though there are currently Global Forecast System (GFS) outputs with resolution of 0.5 °and 0.25 °available, in the present research, border conditions of FNL with coarser resolution were used, with the objective that the physical WRF parameterizations determine the behavior of the forecast in the largest possible measure.On the contrary, for a sensitivity study, analysis data are preferred over prognostic output, which would have introduced implicit additional forecast uncertainty.

Topography Data.
e topographic data Global 30 Arc-Second Elevation (GTOPO30) of the Geological Survey of the United States (USGS), which is a default option in the model, were replaced by the digital elevation model of the Shuttle Radar Topography Mission (SRTM; https://dds.cr.usgs.gov/srtm/version2_1/)[22,23], with a resolution of 90 m, which for South America has an average horizontal error of 9.0 m and an average absolute vertical error of 6.2 m and improves about 10 times (both in spatial resolution and in vertical accuracy) the continental digital elevation model GTOPO30 [24].

Model Configuration.
e simulations are performed for three domains (Figure 1), whose characteristics are specified in Table 1.e unidirectional nesting technique has been applied in the simulations, running all three domains at the same time.e model is initialized at 12 UTC in all cases.
e model was configured with 28 vertical levels, up to approximately 16,000 meters, for which the variable "e_vert" was used in the configuration of the "namelist.input."In this case, 57% of the levels were framed below 6000 meters.
A "control" simulation was defined, and the cumulus, boundary layer, and microphysical schemes were subsequently changed in the rest of the numerical experiments (only one at a time).
To carry out the sensitivity study, we have previously considered previous work results in similar conditions, 2 Advances in Meteorology focusing on orographic complexity and/or tropical condition.us, the parametrization schemes used in the control simulation (CTR in Table 2) and in the rest of experiments were mostly based on the suggestions of the WRF v. 3.7 "User Guide" and in previous results of sensitivity studies and applications of the model in the Andes [12] and in domains with complex topography in tropical and subtropical regions, such as the Himalayas [25,26], the south of India [27], and the Caribbean [28].
For the control experiment, the following parameterization schemes were selected: for microphysics, the ompson scheme, the soil model is the so-called Uni ed  Noah Land Surface Model [29]; for convection, the Grell-Freitas parameterization, it explicitly represents updraft and downdraft and includes the "detrainment" of cloud and ice.
e radiation model was the RRTMG (rapid radiative transfer model for general circulations models) [35], which is a relatively recent version of the RRTM (rapid radiative transfer model) [36], with a better representation of cloudiness, unresolved by the mesh of the model.It applies the k-correlated method to implement an algorithm that is characterized by its speed [37].
Six con gurations of the model were tested, which have been summarized in Table 2: the control con guration (CTR); C_BMJ and C_GRELL3, with Betts-Miller-Janjic and Grell 3D as cumulus parameterizations; MP_MR and MP_LP, where Morrison and Lin-Purdué parameterizations have been used cumulus; and BL_MYJ, where the Mellor-Yamada-Janjic boundary layer scheme was applied.Notice that the acronyms of the con gurations are composed by the type of the scheme that substitutes the one in the control con guration, CTR (C, convection; MP, microphysics; BL, boundary layer), and two or three letters referring to the applied scheme.

Study Period and Veri cation of Results.
e simulations were carried out with a forecast horizon of up to 10 days, for which 9 tens (periods of ten consecutive days) of the months of December, January, and February of the years 2007, 2009, 2010, 2011, and 2012 were selected.All the selected dates belong to the rainy period of the year.e particular tens were chosen among the rainiest for each month.Figure 2 shows the average of 24-hour rainfall for each forecast period for the 19 weather stations considered for punctual veri cation during the 9 selected tens.
e stations considered belong to the observation network of the National Meteorology and Hydrology Service of Peru (SENAMHI), whose data are subject to veri cation.e tens considered were the following: ( e coordinates and altitude of the stations considered for the study are shown in Table 3. e veri cation of the results is carried out in three ways: punctual veri cation the forecast of 24-hour precipitation values, using numerical descriptive statistical measures, the second, punctual veri cation of the occurrence of precipitation as a categorical binary variable, and the third, veri cation of the ability of the model to reproduce patterns of spatial distribution of precipitation.Model output was interpolated using the Cressman method [38] to the stations grid points and compared to the stations grid point with the "in situ" measured data on the basis of the following statistics: bias (B), root mean squared error (RMSE), and mean absolute error (MAE) are calculated as follows: where B k , RMSE k , and MAE k are the bias, root mean square error, and absolute mean error for the k forecast term (k 1, 10), respectively.Pf ij and Po ij are the simulated (forecast) and observed precipitation values, respectively, for the ith station and the jth period of ten days (i 1, 19; j 1, 9)."B" estimates the di erence of the mean value of the numerical forecast relative to the observations; "RMSE" is a measure of the magnitude of the mean error of the simulation with respect to the observations, while "MAE" constitutes the mean absolute deviation of the simulated variable with respect to the observations.e binary categorical veri cation was carried out with the help of the so-called "contingency table," used in [39], but in our case, the fuzzy veri cation technique was not used, but direct point-to-point veri cation was used, which is based on the point-to-point comparison of the predicted and observed elds and the calculation of several indices, such as the probability of detection (POD), the proportion of 4 Advances in Meteorology false alarms (FAR), the bias (B), and the "negative correct" (NC).As a criterion of the precipitation event, a 0.5 mm threshold was used.
In this case, a modi cation to the ACC (accuracy) index is proposed, and the Weighted E ectiveness Index (WEI) is introduced (2), which increases the relative weight of the misses, based on the criterion that a precipitation event that has occurred without being predicted, will have a greater negative impact on society and the economy than a false alarm: WEI hits + negative corrects hits + negative corrects + falses alarms + 1.25 * misses , ( where hits the event was forecast and it happened; false alarms the event was forecast and it did not happen; negative corrects the event was not forecast and it did not happen; and misses the event was not forecast and it happened.e spatial veri cation was carried out using the gridded data of the product 3B42 of the "Tropical Rainfall Measuring Mission" (TRMM 3B42), which is composed of a set of global "multisatellite" precipitation analysis data with a horizontal resolution of 0.25 °and time resolution of 3 h, available between 50 °S and 50 °N [40,41] and of the "Peruvian Interpolate data of the SENAMHI's Climatological and Hydrological Observations" (PISCO) [42].In this sense, it is important to point out that PISCO data have little information in the area of the Peruvian Amazonia (Figure 3), which could in uence the generated precipitation patterns [42].To obtain the spatial distribution of bias, the TRMM grid was adjusted (by interpolation) to the output grid of WRF and subsequently the di erence of the predicted grid and that of TRMM was calculated.

Results of Punctual Veri cation: 24-Hour Cumulative
Precipitation.Table 4 shows the results of the simulations in each domain in terms of forecast quality measures.

Advances in Meteorology
It can be seen that the model overestimated precipitation in all domains regardless of the con guration.
is is consistent with [12], who found that WRF shows a signicant positive bias over the Tropical Andes, when comparing the results of the diurnal cycle produced by the model with reference to the TRMM 3B42 product.
For the 18 km domain, the best results were obtained with the C_GRELL3 scheme (6.63 mm/day), which at the same time showed an "RMSE" of 12.41 mm/day and an "MAE" of 9.37 mm/day.e poorer result was obtained when the conguration C_BMJ was applied, which showed a signi cantly higher "B" than the rest of the experiments (12.08 mm/day).
For the 6 km domain, the MP_MR con guration gave the best results, with a bias of 4.66, "RMSE" of 10.62, and "MAE" of 8.07 mm/day.For this domain, all the con gurations showed a remarkable improvement of their indicators in relation to the domain of 18 km.In the particular case of C_BMJ, it showed bias of 6.71 and "RMSE" of 11.67, very similar to the rest of the con gurations, so it can be partially concluded that this con guration can be used for high-resolution domains, for example, 6 km, but not for lowresolution domains, this is 18 km in this case.
For the 3 km domain, the MP_MR con guration also showed the best indicators with bias of 3.68, "RMSE" of 10.02, and "MAE" of 7.32.e veri cation indicators did not show signi cant improvement for this domain in relation to those achieved in the 6 km domain.
So in general, the MP_MR con guration gave the best results, mainly for the ner domains.
Figure 4 shows the evolution of the average "RMSE" for the 10 days of the simulation horizon, between domains 2 and 3 with each con guration used.It is noted that the "RMSE" increases as the simulation period increases.e con guration C_BMJ and MP_MR shows relatively low RMSE values in all forecast days, growing slowly from day 1 to day 5, but from the sixth day of forecast on, C_BMJ continues increasing the error, while MP_MR curve shows almost no slope until day 10.Here, the MP_MR scheme showed the best behavior regarding RMSE during the whole period.
Figure 5 shows similar results as Figure 4, but for the B indicator.As can be seen, the lowest "B" was shown by the MP_MR con guration, producing the less rainfall overestimation.Figures 4 and 5 indicate that the MP_MR was not only the con guration showing the best results in a general way but also had the most stable behavior throughout the forecast horizon.Based on the point verication, it is further shown that all the tested con gurations overestimate precipitation in the study area.5 shows the statistics of categorical binary veri cation for the di erent con gurations.In a general sense, it can be noticed that the levels of detection of the di erent con gurations are greater than 85%.However, the B was greater than 1, indicating that the model overestimates the presence of precipitations in the region.In the table, it is observed that, for the 6 km domain, the highest POD index (94.39%)was obtained with the C_BMJ con guration, while the highest "negative correct" of 65.31% was obtained with the MP_MR con guration.e highest WEI was that of the C_BMJ scheme.e false alarms and B indices show close values for all the con gurations.
For the 3 km domain, also the C_BMJ con guration showed the best POD (93.07%), with an NC of 62.74; however, the highest NC (62.8%), and WEI (74.18) were shown by the control con guration.In general, for the two domains, the CTR, C_BMJ, and MP_MR con gurations were the most e ective.
Additional veri cations were performed, only for the stations located in the northern, central, and southern    regions of the basin, both quantitatively and qualitatively.However, the results obtained were similar to those shown considering all the basin stations, and we do not show them in this paper.
As we had mentioned in the introduction, from the results of the experiments with the di erent parameterizations, a new experiment was carried out using the MP_MR con guration (which showed the best B and RMSE indicators for the entire forecast period) with the objective to verify the need or not to use CPS for the 3 km domain in the study region.In this sense, the results are shown in Table 6.In the table, the number of the tens corresponds to the order indicated in the Data and Methodology, and the best results for each ten are indicated in bold type.As can be seen, when looking at all the tens, the best results were obtained in some cases when CPS was used, and in others, when the cumulus scheme was deactivated.So, this result con rms the idea that it is not appropriate to assume that, for domains of less than 5 km, the parameterization of clusters must be deactivated, as was concluded in [21].

Results of the Veri cation of the Spatial Precipitation Field
Forecast.Figure 6 shows the average spatial distribution of the precipitation from the TRMM and PISCO databases and the model output for the di erent con gurations used for the studied periods, corresponding to the 6 km domain.Regarding the data of TRMM and PISCO (Figures 6(a) and 6(b)), it can be seen that both databases show very similar patterns of precipitation.In this sense, the greater di erence is that TRMM shows a maximum of rainfall around 12 °S and 73 °W, which is less notorious in PISCO, due to the fact that in that region there are scarce meteorological stations.Two maximum precipitation regions have been highlighted in red circles in Figures 6 and 7 in the zone of Amazonia, to the east of the mountain range.Another precipitation maximum of smaller area (also surrounded by a red circle in the gures) is located at the southeast of the Mantaro Basin, which contour line is represented by a thick black line in the gure.
An elongated rainy area which starts from the south of the basin and extends towards the southeast has been indicated with a thick red curve in each gure.Another signi cant feature is the relative spatial rainfall minimum shown by PISCO along part of the eastern boundary of the basin, mostly in its northern half, indicated by black circles in Figures 6(a  It can be noticed that, in general, all con gurations reproduced correctly the precipitation patterns shown by PISCO and TRMM but overestimating the reference accumulations, shown in Figures 6(a) and 6(b), mainly on the eastern slope of the mountain range (hot spots).However, it must be mentioned that it has been found [43] that in this sector of the Peruvian mountain range, TRMM underestimates precipitation by at least 40%, so that the overestimation of WRF relative to the "real precipitation value" in that case could be much lower.e con guration which best reproduced the maximum in 12 °S-72 °W, relative to TRMM was C_BMJ.
It must be pointed out that none of the con gurations reproduced the less rainy strip shown by PISCO and, to a lesser extent, also by TRMM, along the eastern limit of the basin in its northern half (black circles in Figures 6(a) and 6(b)).In fact, WRF produces more precipitation in the eastern sector of the basin than in the west, which matches the observations only for the southern half of the basin.
Figure 7 shows the zoomed-in view of TRMM and PISCO precipitation patterns and WRF output for the 3 km domain.In TRMM and PISCO (Figures 7(a) and 7(b)), local maxima and minimum of precipitation are indicated as in Figure 6, using red circles fox maxima and with a black circle for the minimum of precipitation in the eastern limit of the northern half of the basin, practically covering the valley (which is represented by a small closed black contour near the western boundary of the basin).In a general sense, it is observed that the basin in its southern half is rainier towards the east, while in the northern half it is rainier towards the west, which can be seen more clearly in the PISCO data.It can also be noticed, more clearly in TRMM, that the south of the valley is drier than the north.
As shown in Figures 7(c)-7(h), the model reproduces the precipitation nuclei to the east of the basin, although in a general sense it overestimates accumulated rainfall.It can be seen that toward the southern half of the basin, the MP_MR con guration reproduced more clearly the fact that the eastern sector is rainier than the western.e CTR con guration also reproduces this pattern but much more extended to the west.e rest of the con gurations were less precise in this regard.In general, the MP_MR con guration was the one reproducing best the spatial distribution of rainfall in the south of the basin.For the northern half of the basin, all the con gurations re ected the eastern sector as rainier than the western, which does not reproduce the patterns in TRMM and PISCO.
From the above considerations, it can be concluded that for the 6 km domain, the precipitation pattern shown by In order to quantify the spatial distribution of the bias of the model, the bias (B) maps of both domains were made for di erent forecast horizons.However, in this case, bias will be calculated relative to TRMM data, considering that PISCO re ects a signi cant de cit of rainfall to the southeast of the basin.Figure 8 shows the spatial B distribution for the second day of forecast of the 6 km domain, which con rms that in general all the con gurations overestimate the precipitation; however, it is noticed that it underestimates both for the western sector of the mountain range and for the coast, although this sector of the Andes is much less rainy than its eastern sector.In this case, it can be seen that the MP_MR con guration produces the least overestimation, as we had observed in the general picture.
Figures 9 and 10 show the results for the sixth and tenth days of forecast, where it is visually noticed that the MP_MR con guration has the lowest bias.In this case, it is signi cant that, for the sixth day the con guration, C_BMJ shows a noticeable improvement in relation to the second, specically in the Amazonia zone, with insigni cant overestimations.
Figures 11-13 are analogous to Figures 8-10, but for the 3 km domain.As for the 6 km domain, in this case, the model overestimates precipitation, and even in the interior of the basin, some of the con gurations overestimate more than in the 6 km domain.
For this domain, in all cases, it is also noticed that MP_MR showed the best results, but with the particularity that outside the basin, B increased with forecast horizon, while the behavior inside the basin was more stable.us, from a spatial point of view, the model clearly shows an overestimation in the Mantaro basin and, in general, in the entire eastern sector of the ridge, while in its western sector, underestimation of precipitation was observed.From the spatial point of view, it is con rmed that MP_MR was the con guration with the lowest overestimation, behaving in a stable way from the second to the tenth forecast day, with a moderate increase in RMSE as the forecast horizon increased.From the above arguments, it can be concluded that the sensitivity to the microphysical scheme was higher than to the convective scheme, which indicates the importance of the microphysics parameterization in WRF for precipitation forecasting tasks in the study region.As an example, Figure 14 shows the vertical pro le of rain water mixing ratio (Qrain) produced by the 3 microphysical schemes used for the Huayao station and, in parentheses, the predicted rainfall in each case, averaged for all tens used in study.It can be noticed the correspondence between "Qrain" and predicted precipitation.
Similar to the results obtained in this research, in [25], a sensitivity study was made of simulated monsoon precipitation to the cloud microphysics schemes in WRF for the summer period in the valley of Langtang, Himalaya.It was obtained that, in a general sense, the model underestimated the accumulated precipitation in 10 days, and also the Morrison microphysical scheme showed the best results.
Our results for the Andes are consistent with the nding of [26], where it was concluded that, in general, the doublemoment Morrison's scheme represents more correctly the microphysics processes in the Himalaya, possibly because of the better representation of the warm and cold processes related with the formation and evolution of the parameters of particle size distributions, in comparison with onemoment parameterizations.

Summary and Conclusions
Simulations were developed with the regional model WRF, for domains of 18, 6, and 3 km of spatial resolution, with the main objective of determining the ability of several congurations of the model to forecast the eld of precipitation in the short and medium terms in the complex orography conditions of the Central Andes of Peru.
e initial and boundary conditions were generated by the NCEP "Global Operational Analysis," nal analysis FNL, every 6 hours, with horizontal resolution of 1 °× 1 °. e behavior of the numerical forecasts for the rst 10 days of simulation was analyzed, concluding that the mean square error for the higher resolution domains increases with the simulation term, in general until the tenth day.e lowest mean square error was obtained for the MP_MR con guration (Morrison for microphysics and Grell-Freitas for convection) during the whole period, remaining with little variation after the sixth forecast day.
is con guration also provided the lowest positive bias, although in all cases a clear overestimation of rainfall was observed.12 Advances in Meteorology e point veri cation showed that, in general, all the tested con gurations overestimate rainfall in the region.In this sense, the MP_MR con guration performed better, according to the bias and RMSE indicators, which alludes not only to the entire forecast period but also to all the forecast terms.It is important to notice that from the sixth day on, the RMSE for MP_MR did not increase more, which does not apply to the rest of the schemes used.Another important result in this case is that the C_BMJ con guration ( ompson microphysics and Betts-Miller-Janjic convection) showed poor results for the 18 km domain, so it is recommended not to use it in this region for low-resolution domains.
e categorical binary veri cation showed detection rates of precipitation above 85% in all cases.e B was all positive, which indicates that the model overestimated the presence of rainfall in the region.In general, the con gurations showing the best indicators were C_BMJ and MP_MR.
e spatial veri cation showed that all the schemes produced rainfall distribution patterns quite similar to those of TRMM and PISCO, although in all cases with positive bias in precipitation accumulations.In this sense, the con guration that less overestimated precipitation was MP_MR.It also reproduced better precipitation pattern in the southern half of the basin.Regarding the northern half of the basin, all the con gurations showed the eastern sector as rainier than the west, which does not agree with TRMM and PISCO.
e bias spatial distribution generally showed positive values in the basin and in general throughout the eastern sector of the Andes and the Amazonia, while towards the west of the Andes and the coast, it showed values close to zero or negative.
e con guration with lower bias was MP_MR, which also showed the lowest values of RMSE, in both cases for the two domains analyzed.e C_BMJ con guration, although with higher bias and RMSE, showed improvement for the greater terms, which is an indicator to be taken into account for its use in longer-term forecast.
Based on the results, it is concluded that the MP_MR and C_BMJ con gurations were the ones with the best results in a general sense.Consequently, the MP_MR con guration is recommended for short-and medium-term rainfall forecasting tasks in the Central Andes of Peru and particularly in the Mantaro basin.
e Betts-Miller-Janjic cumulus parameterization will be further investigated by the authors for the region in combination with other physical schemes to explore its possible application potential in medium range precipitation forecasting, using multicon guration ensembles.It was also concluded that the activation or not of cumulus parametrization for the domain of 3 km resolution is not relevant for the precipitation forecast in the study region.Advances in Meteorology

Figure 1 :
Figure 1: Domains used in the simulations and distribution of the stations considered for the veri cation of the results.

Figure 2 :
Figure 2: Mean of accumulated 24-hour precipitation for each forecast period in the 19 meteorological stations considered during the 9 selected tens.

Figure 3 :
Figure 3: Network of stations used for the elaboration of PISCO (taken from [42]), with percentages of data per station, considering the period 1981-2016.

Figure 4 : 3 . 2 .
Figure4: Evolution of the "RMSE" average of domains 2 and 3 of all stations with each con guration scheme used, for the 10 days of simulation horizon.

Figure 5 :Figure 7 :
Figure 5: Evolution of the average "B" of domains 2 and 3 with each con guration, for the 10 days of the simulation horizon.
) and 6(b).In TRMM, this zone is less extended to the north.e maps in Figures6(c)-6(h) show the output of WRF for the di erent con gurations used in the experiments.

Table 1 :
Main characteristics of domains, initial, and boundary conditions.

Table 2 :
Parameterization schemes used in the experiments.

Table 3 :
Stations used to validate the simulations.

Table 4 :
Statistics of the point veri cation of precipitation forecast for the 3 simulated domains and 6 combinations of parametrizations, including the control experiment, average of all the stations considered.

Table 5 :
Statistics of the point veri cation for the domains 2, 3, and 6 combinations of parametrizations, including the control experiment, average of all the stations considered.
e smaller errors are given in bold.

Table 6 :
Statistics of the point quantitative veri cation for every ten, simulated with CPS and without CPS in domain 3.