The ability of the WRF-ARW (Weather Research and Forecasting-Advanced Research WRF) model to forecast extreme rainfall in the Central Andes of Peru is evaluated in this study, using observations from stations located in the Mantaro basin and GOES (Geostationary Operational Environmental Satellite) images. The evaluation analyzes the synoptic conditions averaged over 40 extreme event cases, and considers model simulations organized in 4 nested domains. We first establish that atypical events in the region are those with more than 27 mm of rainfall per day when averaging over all the stations. More than 50% of the selected cases occurred during January, February, and April, with the most extreme occurring during February. The average synoptic conditions show negative geopotential anomalies and positive humidity anomalies in 700 and 500 hPa. At 200 hPa, the subtropical upper ridge or “Bolivian high” was present, with its northern divergent flank over the Mantaro basin. Simulation results show that the Weather Research and Forecasting (WRF) model underestimates rainfall totals in approximately 50–60% of cases, mainly in the south of the basin and in the extreme west along the mountain range. The analysis of two case studies shows that the underestimation by the model is probably due to three reasons: inability to generate convection in the upstream Amazon during early morning hours, apparently related to processes of larger scales; limitations on describing mesoscale processes that lead to vertical movements capable of producing extreme rainfall; and limitations on the microphysics scheme to generate heavy rainfall.

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 short- and medium-term precipitation forecast is conducted in the Central Andes of Peru. Lin-Purdué, Thompson, and Morrison microphysics schemes were tested, as well as the Grell–Freitas, Grell 3d, and Betts–Miller–Janjic cumulus parameterizations. The tested boundary layer schemes were the Yonsei University and Mellor–Yamada–Janjic. A control configuration was defined, using the Thompson, 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 the Morrison 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 short- and medium-term rainfall forecasting tasks in the Central Andes of Peru and particularly in the Mantaro basin. The 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.

Information on the vertical structure of rain, especially near the surface is important for accurate quantitative precipitation estimation from weather and space-borne radars. In the present study, the rainfall characteristics, from a vertically pointed profile Radar in the Mantaro basin (Huancayo, Peru) are observed. In summary, diurnal variation of near-surface rainfall and bright band height, average vertical profiles of the drop size distribution (DSD), rain rate, radar reflectivity (Z_e) and liquid water content (LWC) are investigated to derive the rainfall characteristics. Diurnal variation of rain rate and bright band height show the bimodal distribution, where frequent and higher rain rate occurred during the afternoon and nighttime, and more than 70% bright band height found between 4.3–4.7 km. The average vertical profiles of Z_e show the opposite characteristics above and below the melting level (ML) and depend on the near-surface rain rate. For example, the average Z_e profiles have a negative gradient above the ML, whereas below, the ML, the gradient depends on the near-surface rain rate. The rain rate and LWC show the opposite behavior, and both consist of a positive (negative) gradient below (above) the ML. The vertical growth of DSD parameters depend on the near-surface rain rate, and a higher concentration of large-sized of droplets are observed for higher near surface rain rate, however, the dominant modes of droplets are <1 mm throughout the vertical column. However, the most significant variation in DSD growth is observed for near-surface rain rate ≥20 mm/h. These findings suggest using different retrieval techniques for near surface rain estimation than the rest of the vertical profile and high rain rate events. The improved understanding of the tropical Andes precipitation would be very important for assessing climate variability and to forecast the precipitation using the numerical models.

The aim of the present study is to analyze the triggering mechanisms of three thunderstorms (TSs) associated with severe rainfall, hail and lightening in the tropical central Andes of Peru, specifically above the Huancayo observatory (12.04 ${}^{\circ }$ S, 75.32 ${}^{\circ }$ W, 3313 m a.s.l.) located in the Mantaro valley during the spring-summer season (2015–2016). For this purpose, we used a set of in-situ pluviometric observations, satellite remote sensing data, the Compact Meteorological Ka-Band Cloud Radar (MIRA-35C), the Boundary Layer Tropospheric Radar and downscaling model simulations with the Weather Research and Forecasting (WRF) Model (resolutions: 18 km, 6 km and 2 km), and the Advance Regional Prediction System (ARPS) (resolution: 0.5 km) models in order to analyze the dynamic of the atmosphere in the synoptic, meso and local scales processes that control the occurrence of the three TS events. The results show that at synoptic scale, the TSs are characterized by the southern displacement of the South-east Pacific Subtropical Anticyclone up to latitudes higher than 35 ${}^{\circ }$ S, by the weakening and south-eastern displacement of the Bolivian high–North east low system and by the intrusion of westerly winds along the west side of the central Andes at upper and medium levels of the atmosphere. At meso-scale, apparently, two important moisture fluxes from opposite directions are filtered through the passes along the Andes: one from the north-west and the other from the south-east directions converge and trigger the deep convection into the Mantaro valley. These moisture fluxes are generated by the intrusion of the sea-breeze from the Pacific ocean along the west of the Andes coupling with upper and middle westerly winds and by the thermally induced moisture fluxes coming from the South American low level jet at the east side of the Andes. At the local scale, there is a low-level conditional instability in the previous hours as well as during the occurrence of the TSs above the Huancayo observatory. In addition, the simulation results indicated the possibility of generation of inertial gravity waves in the Amazon basin, associated with geostrophic adjustment which transports energy and moisture into the central Andes plateau and consequently intensifies the thunderstorms above the Mantaro valley.

The present study presents a detailed analysis of the diurnal and monthly cycles the surface boundary layer and of surface energy balance in a sparse natural vegetation canopy on Huancayo observatory (12.04 ${}^{\circ }$ S, 75.32 ${}^{\circ }$ W, 3313 m ASL), which is located in the central Andes of Perú (Mantaro Valley) during an entire year (May 2018–April 2019). We used a set of meteorological sensors (temperature, relative humidity, wind) installed in a gradient tower 30 m high, a set of radiative sensors to measure all irradiance components, and a set of tensiometers and heat flux plate to measure soil moisture, soil temperatures and soil heat flux. To estimate turbulent energy fluxes (sensible and latent), two flux–gradient methods: the aerodynamic method and the Bowen-ratio energy-balance method were used. The ground heat flux at surface was estimated using a molecular heat transfer equation. The results show minimum mean monthly temperatures and more stable conditions were observed in June and July before sunrise, while maximum mean monthly temperatures in October and November and more unstable conditions in February and March. From May to August inverted water vapor profiles near the surface were observed (more intense in July) at night hours, which indicate a transfer of water vapor as dewfall on the surface. The patterns of wind direction indicate well-defined mountain–valley circulation from south-east to south-west especially in fall–winter months (April–August). The maximum mean monthly sensible heat fluxes were found in June and September while minimum in February and March. Maximum mean monthly latent heat fluxes were found in February and March while minimum in June and July. The surface albedo and the Bowen ratio indicate semi-arid conditions in wet summer months and extreme arid conditions in dry winter months. The comparisons between sensible heat flux ( $Q_H$ ) and latent heat flux ( $Q_E$ ), estimated by the two methods show a good agreement (R ${}^2$ above 0.8). The comparison between available energy and the sum of $Q_E$ and $Q_H$ fluxes shows a good level of agreement (R ${}^2$ = 0.86) with important imbalance contributions after sunrise and around noon, probably by advection processes generated by heterogeneities on the surface around the Huancayo observatory and intensified by the mountain–valley circulation.