atmosphere

Article
Using the Weather Research and Forecasting (WRF)
Model for Precipitation Forecasting in an Andean
Region with Complex Topography
Gonzalo Yáñez-Morroni 1 ID
, Jorge Gironás 1,2 * ID
, Marta Caneo 3 , Rodrigo Delgado 3
and René Garreaud 4
1 Departamento de Ingeniería Hidráulica y Ambiental, Pontificia Universidad Católica de Chile,
Vicuña Mackenna 4860, Macul, Santiago 7820436, Chile; [email protected]
2 Centro de Investigación para la Gestión Integrada de Desastres Naturales (CIGIDEN), Centro de Desarrollo
Urbano Sustentable (CEDEUS), Centro Interdisciplinario de Cambio Global UC (CCG-UC),
Vicuña Mackenna 4860, Macul, Santiago 7820436, Chile
3 Dirección Meteorológica de Chile (DMC), Portales 3450, Estación Central, Santiago 9170018, Chile;
[email protected] (M.C.); [email protected] (R.D.)
4 Departamento de Geofísica, Centro de Ciencia del Clima y la Resiliencia, Universidad de Chile,
Blanco Encalada 2002, Santiago Centro, Santiago 8370449, Chile; [email protected]
* Correspondence: [email protected]; Tel.: +562-2-3544227




Received: 14 June 2018; Accepted: 30 July 2018; Published: 2 August 2018


Abstract: The Weather Research and Forecasting (WRF) model has been successfully used in weather
prediction, but its ability to simulate precipitation over areas with complex topography is not optimal.
Consequently, WRF has problems forecasting rainfall events over Chilean mountainous terrain and
foothills, where some of the main cities are located, and where intense rainfall occurs due to cutoff
lows. This work analyzes an ensemble of microphysics schemes to enhance initial forecasts made
by the Chilean Weather Agency in the front range of Santiago. We first tested different vertical
levels resolution, land use and land surface models, as well as meteorological forcing (GFS/FNL).
The final ensemble configuration considered three microphysics schemes and lead times over three
rainfall events between 2015 and 2017. Cutoff low complex meteorological characteristics impede the
temporal simulation of rainfall properties. With three days of lead time, WRF properly forecasts the
rainiest N-hours and temperatures during the event, although more accuracy is obtained when the
rainfall is caused by a meteorological frontal system. Finally, the WSM6 microphysics option had the
best performance, although further analysis using other storms and locations in the area are needed
to strengthen this result.

Keywords: WRF forecast; rainfall; complex topography; Andean watershed; flash floods

1. Introduction
Natural disasters with hydrological, meteorological or climatological origin are a severe and
worldwide problem which causes loss of life and property damage [1]. Global warming appears
to be positively correlated with future flood risk at global scale: an increase of 4 ◦ C will drastically
increase flood risk in several countries representing ∼70% of the world population [2]. Unevenness in
future rainfall projections over South America doesn’t facilitate future flood risk projections, although,
in the past two decades, several flood events took place in Andean regions [3]. For example, Chile’s
Atacama Desert—the driest desert in the world—recorded 65 mm of rainfall in just three days in March
2015 [4,5]. This rainfall unleashed a torrential flood due to the combined effect of high temperatures,
a steep topography and erosion favored by the soil granulometry and its infiltration capacity [6].

Atmosphere 2018, 9, 304; doi:10.3390/atmos9000304 www.mdpi.com/journal/atmosphere

Atmosphere 2018, 9, 304 2 of 23

The event caused several human lives to be lost, the interruption of water and electricity supplies for
weeks and an estimated economic cost of ∼US$1.5 billion [7].
Numerical weather prediction (NWP) models are broadly and successfully used for weather
prediction and research. However, the simulation of rain events over complex orography
(i.e., mountaneous regions) is still a challenge, as NWP models may not resolve the underlying
topography within a high resolution [8,9]. Furthermore, complex topography affects the meteorology
by modifying the surface heat flux, albedo, wind speed and direction [10–13]. A description of
this complex topography helps to understand the behavior of the planet boundary layer (PBL) and
microphysics (MP) schemes, allowing to study the slope wind-flow and valley wind-flow. However,
none of the PBL parameterizations accurately predict the abrupt wind speeds and temperature
profiles near the surface [14–16]. Increasing the vertical resolution of the computation grid near
the ground-surface has been done to improve the modeling results in places with complex topography
[17]. The finer grid provides more details in complex orographic zones, and a better performance
depending on the regional climatology [18]. Nevertheless, a finer domain grid resolution has a limited
impact on the traditional verification scores, and it doesn’t always improve the rainfall forecast.
A method for improving the performance of NWP models is the use of an ensemble
approach [19–22]. This approach consists of repeating simulations under the same initial and boundary
conditions but varying the physics scheme parametrization each time. This parametrization typically
includes: radiative transfer, vegetation and soil characteristics, microphysics, and flux interaction
of heat, moisture and momentum in the soil/atmosphere interphase, among others. Ensembles’
performances have been widely used by different authors [23–27] for capturing climate projections
uncertainties, or achieving an optimal physics scheme parametrization. On the other hand, forecasting
models are in general very dependent on the local condition (i.e., topography, hydrology, time of
year, etc.). To overcome this issue, an iterative testing exploring different parameters’ values is
a suitable tool for knowing which parameters need to be specified and investigated with greater
accuracy [28]. Furthermore, a realistic model performance cannot be attributed to the achievement of
a single scheme parametrization. Due to the complexity and nonlinearity of the atmospheric equations,
the totality of the schemes’ parametrizations is involved in the output performance.
Another issue commonly tested for rainfall prediction over big watersheds is the forecast lead-time,
or the time prior to the forecast date. An appropriate lead-time can considerably enhance the simulation
accuracy. Rainfall forecast performance is sensitive to the temporal and spatial scale, becoming worse
with lead-times longer than five to six days, and eventually meaningless after nine days [29]. Advances
in NWP have made lead times of two weeks feasible [30], but forecast is more reliable for shorter
lead-times, i.e., three days, for lighter rainfalls within bigger basins [29]. Hence, a rainfall forecast in
a small watershed for heavy rainfalls is a challenge, even when shorter lead-times are considered.
On the other hand, models also require 6–12 h for spin up, so if the initialization is too close in time to
the actual storm onset, the model may miss some early-stage precipitation.
Even though many studies using the Weather Research and Forecasting (WRF) model [31] have
already embraced its limitations over complex topography [32–36], few studies have focused on
the South American Andean mountainous region, where the Andes cordillera reaches more than
5000 m.a.s.l. within 200 km from the Pacific sea border. WRF was used over the Nahuelbuta Mountains,
in coastal southern Chile (37◦ S–38◦ S), to successfully simulate observed seasonal and daily mean
rainfall distributions [37]. WRF was also used to study the direct effect of the Andean topography
on wind speed and direction over the Argentinian foothill (Mendoza) [38]. Finally, WRF was used to
forecast urban PM10 and PM2.5 pollution events over Santiago’s foothill, but their physical scheme
combination did not comprise any representation of rainfall events [39]. They only briefly discussed
the modeling of the coastal lows development into Santiago’s valley. On the other hand, the Chilean
Meteorological Agency (DMC in Spanish) performs daily WRF simulations with a physics scheme
calibrated to predict weather in Chile Central Valley, a 1200 km length region, which barely captures
Atmosphere 2018, 9, 304 3 of 23

the Andean mountain and foothill topography. These simulations are made using a 6 km horizontal
resolution and 50 vertical levels, more densified in the top and bottom boundaries.
The objectives of this work are to study and test different WRF configurations and physic
processes parameterizations for the forecasting of rainfall over mountainous Central Chile, as well
as to test the effects of different lead-times. The performance of the WRF simulations is evaluated
according to their ability to represent relevant characteristics of the temperature dynamics and rainfall
events for prediction of floods (i.e., length, peak intensity, concurrent freezing level and the rainiest
N-consecutive hours) over the Andean front range in central Chile, particularly the Quebrada de
Ramón catchment, a small basin located in the area. In our analysis, we used three storm events with
different characteristics that took place between 2015 and 2017.

2. Data and Methods

2.1. Field Data

The study area corresponds to the Quebrada de Ramón basin, a 36 km2 mountainous watershed
located in central Chile whose highest elevation is 3250 m. Its complex and steep topography has
a maximum elevation gradient of 220 m per 1 km and a high average slope (∼38%) [40]. Its outlet is
located at an elevation of 800 m, in the east part of Santiago, the capital of Chile. The urban dynamics
in this area neither considered high return period flows nor detrital floods. A clear example of the
consequences took place on 3 May 1993, when a warm storm produced large floods and landslides in
this and other front range basins nearby [41]. The flood injured 3500 people and caused 50 human
losses. The situation was aggravated by the presence of houses located in the Piedmont and the
floodplains [42].
Santiago (33◦ 300 S; 70◦ 420 W) exhibits a semi-arid Mediterranean climate, with a mean annual
precipitation of 310 mm and about 25 rainy days during the austral winter (May–September) and warm,
dry summers [43]. This general weather pattern has remained constant over time, since no trends were
found on monthly rainfall, mean and maximum annual precipitation for the period 1950–2018 using
the Mann–Kendall test [44,45], with a 5% significance level. Nevertheless, central Chile—and Santiago
in particular—have experienced a prolonged drought since 2010 onwards that seems to herald drier
conditions in the future [46].
For this study, data were collected from the Apoquindo meteorological station, located near
Quebrada de Ramón’s centroid (33◦ 260 S; 70◦ 280 W) at an elevation of 1625 m.a.s.l. Six meteorological
variables are registered every 10 min: temperature, relative humidity, solar radiation, wind speed and
direction, and rainfall (Figure 1). Moreover, hourly and daily precipitation and temperature data from
13 meteorological stations located in Santiago and its surroundings were used to verify the orographic
gradient and to analyze the return period of the storms events (Table 1).
The freezing level location (i.e., a elevation representative of the snow line that controls the spatial
occurrence of rain and snow) was obtained from the Santo Domingo atmospheric sounding station
record [47]. Soundings are made twice a day (00Z and 12Z) at 105 km from the city of Santiago
(33◦ 650 S; 71◦ 610 W). These measurements are considered to describe Santiago atmosphere’s features,
as the city is located within the area of measurements’ representativeness (defined by a 150 km radius),
based on the negligible effects of topography and near coast interactions in the free troposphere [48].
Atmosphere 2018, 9, 304 4 of 23

Figure 1. Quebrada de Ramón’s basin (outlined in red), in Santiago de Chile’s Piedmont, and the
location of Apoquindo meteorological station and the other 13 meteorological stations in Santiago’s
valley and foothills. White crosses are latitude/longitude coordinates.

Table 1. Meteorological stations names, coordinates and recorded variables: T (air temperature) and P
(precipitation). Stations located inside Santiago are indicated with *.

Name ID Latitude Longitude Elevation (m.a.s.l.) Variables

San José Guayacán SJ 33◦ 370 S 70◦ 210 W 928 Hourly T, P
Apoquindo AP 33◦ 270 S 70◦ 280 W 1625 Hourly T, P
Quebrada de Macul * QM 33◦ 300 S 70◦ 310 W 950 Daily T, P
Antupirén * AN 33◦ 300 S 70◦ 310 W 904 Daily P
Cerro Calán * CC 33◦ 240 S 70◦ 320 W 904 Daily T, P
Tobalaba * TO 33◦ 270 S 70◦ 330 W 650 Daily T, P
La Platina * PL 33◦ 340 S 70◦ 380 W 630 Hourly T, P
Quinta Normal * QN 33◦ 270 S 70◦ 410 W 534 Daily T, P
Lo Pinto * PI 33◦ 160 S 70◦ 440 W 512 Hourly T, P
San Pablo * SP 33◦ 270 S 70◦ 450 W 490 Hourly T, P
Pudahuel PU 33◦ 240 S 70◦ 480 W 482 Daily T, P
Rinconada de Maipú RM 33◦ 300 S 70◦ 510 W 462 Hourly T, P
Hacienda Lampa HL 33◦ 170 S 70◦ 510 W 493 Hourly T, P
El Paico PA 33◦ 420 S 70◦ 600 W 275 Hourly T, P

2.2. Meteorological Characterization of Rainfall Events

For the the scheme parameterization and calibration of WRF, we considered three rain events:
a frontal system rainfall (19 October 2015, hereafter referred to as OCT15), a cutoff low episode (17
Atmosphere 2018, 9, 304 5 of 23

April 2016, hereafter referred to as APR16) and a hybrid between cutoff low and frontal system rainfall
(11 May 2017, hereafter referred to as MAY17). Cutoff lows are meteorological phenomena that take
place in medium to high latitudes (20◦ S–40◦ S), more likely during autumn and winter (i.e., from April
to September) [49,50]. They can cause extreme cold weather and precipitation in high elevations above
1000 m, coupled with strong winds and occasional thunderstorms [51]. These three rainfall events are
considered to be representative of the possible range of rainfall events in this region.
Meteorological forcing data were obtained from NCEP/NCAR Reanalysis. In this study, we
use the period 1981–2010 for computing the mean anomalies during the studied rainfall events.
An anomaly shows the difference of any meteorological variable from its long-period mean value for
a given location.

2.2.1. 19 October 2015 Rainfall Event (OCT15)

Between 16–20 October 2015 (largest intensities between 19 October 2015 11:00 a.m. UTC and 20
October 2015 12:00 a.m. UTC), Chilean central valleys were affected by an event that left a relatively
high precipitation over the mountainous region (i.e., 46 mm in the Apoquindo gauge vs. 21 mm in
Quinta Normal station, Table 1), whose rainiest day had a return period of 2.3 years according to the
Quinta Normal records. This event was caused by a cutoff low developed over the Southeast Pacific
(80◦ –90◦ W and 30◦ S), which produced a nucleus of negative anomalies in the 500 hPa geopotential
level, with a magnitude of −150 m. The event was intensified by the jet stream at 250 hPa, with its
nucleus of largest intensities (30 m/s positive anomalies) located just above the Chilean central region
(30◦ S) (Figure 2, upper row).

2.2.2. 17 April 2016 Rainfall Event (APR16)

In 15–17 April 2016 (15 April 2016 9:00 a.m. UTC to 17 April 2016 8:00 p.m. UTC), an extreme
rainfall event affected Chile’s central region. On 17 April, 108 mm of rainfall were measured at
the Apoquindo meteorological station in 24 h. According to Quinta Normal historical records,
the maximum hourly rainfall of this event had a 46-year return period, even though the accumulated
rainfall of 14 h within the storm peak had a return period of 10 years [52].
The frontal system moved from north to south, with negative anomalies in the atmospheric
pressure and geopotential height at 500 hPa (−150 geopotential meters). At the same time,
an intensification of both subtropical and polar jet streams took place at 250 hPa, allowing higher
wind intensities in the system nucleus (positive anomalies of 35 m/s in the wind vectors), favoring the
cloudiness development over central Chile (Figure 2, central row).
Satellite images showed a significant moisture contribution, advected from lower latitudes, which
incorporated more precipitable water into the system. This quasi-stationary weather front led to
a warm winter storm-like disturbance [53], allowing more rainfall to precipitate over central Chile.
During most parts of a warm storm, precipitation is caused by mechanical uplift of moist air over the
Andes, developing its maximum precipitation in mid-elevations inlands. This behavior is consistent
for several warm storms in other world regions, where heavy rainfall and flood events happen [53].
Additionally, the freezing level was above 3000 m, causing a large contribution of liquid
precipitation at relatively high temperatures. In fact, during 16 April, the Mapocho river in Santiago’s
basin, flooded the central part of the city, where a riverbed modification due to a construction took
place. Because of its large duration, maximum intensity and total amount, this rainfall event was
used for initial WRF modifications over the grid discretization, land use and initial and boundary
condition data.
Atmosphere 2018, 9, 304 6 of 23

Figure 2. NCEP/NCAR Reanalysis composite anomaly for OCT15 (upper row), APR16 (middle row)
and MAY17 (bottom row) events, for 500 mb geopotential height (m) (left column) and 250 mb vector
wind (m/s) (right column).
Atmosphere 2018, 9, 304 7 of 23

2.2.3. 11 May 2017 Rainfall Event (MAY17)

With 53 mm of rainfall recorded in Apoquindo station in two days, the event of 11 May 2017 4:00
a.m. UTC to 12 May 2017 1:00 a.m. UTC affected a large portion of Chile. This event corresponds
to a return period of three years in Quinta Normal hourly records. In the previous days, the system
presented a trough with closed circulation at 500 hPa, a typical cutoff low characteristic. Additionally,
the jet stream at 250 hPa had a strong zonal direction in the trough posterior part (Figure 2, bottom row).
During the rainfall event, a low pressure center on the surface associated with the weather front
was observed. At higher elevations, the typical cutoff low pattern was not visible, as there was no clear
closed-circulation and the jet stream was coupled to westerlies. At 500 hPa, the trough contributed
with additional divergence, generating more cloudiness over the affected area.

2.3. Meteorological Forcing Data

ARW-WRFv 3.5 [31], a state-of-the-art mesoscale NWP model was used in this work. The model is
suitable for a wide range of applications, such as weather routine forecast, research simulations,
and evaluation of parameters in simulated systems. WRF solves the scalar conservation and
compressible non-hydrostatic Euler equations, through vertical coordinates (η levels) with a variable
grid density. WRF system is supported and maintained by NCAR [54].
Initial and boundary meteorological conditions according to the NCEP operational Global Forecast
System (GFS) and the NCEP Final Operational Global Analysis data (FNL) were initially compared.
GFS is composed of four blocks (atmosphere, ocean, land/soil and sea/ice), providing a 0.5◦ resolution
grid with the forecast of atmospheric and land–soil variables every 3 h for modeling periods after
2006 [55]. Operational forecasters, such as the DMC in Chile, use GFS as boundary condition to
simulate precipitation at finer resolution via NWP. On the other hand, FNL is a 1◦ resolution grid
map produced every 6 h, prepared approximately an hour after GFS started, which allows more
observational data to be used (∼10%) in the upgrade of initial and boundary conditions. FNL data are
available for periods after 1999 [56]. Finally, sea surface temperature (SST) at a 0.5◦ resolution from the
NCEP SST analyses database were also provided to the WRF model.
An initial testing over GFS/FNL meteorological forcing considered three different lead times
for APR16 event (i.e., 120, 96 and 72 h): although both datasets performed similarly, the rainfall total
length provided by GFS was more feasible. Thus, only the GFS dataset was used in the simulations.

2.4. WRF Model and Physics Schemes

Initial testing over the APR16 event were made focused over (1) two vertical level resolutions,
(2) two land use and three land surface models. Results from these initial simulations allowed the
definition of the final schemes parameterization and grid resolution, which was used in a final set
of simulations to test different MP schemes and lead times. The forecasts were carried out with a
32-core machine (2.30 GHz Intel Xeon E5-2698 v3 in Santiago, Chile), taking less than 5 h to retrieve
operational daily forecasts at a resolution of 6 km. This time considers 1 h in data preparation and
downloading of initial and boundary conditions, (GFS 0.5◦ for DMC forecasts), 2.5 h in running WRF
and 45–60 min in post-processing.

2.4.1. Simulation Domains and Topography Complexity

The domain of the WRF simulations was composed of three nested grids with 54, 18 and 6 km of
resolution. These were based on a Lambert Projection centered at Santiago de Chile, interacting with
each other through a two-way nesting strategy (Figure 3). The biggest or parent domain (d01, 54 km)
covered South America’s western region, the Pacific and part of the Atlantic Ocean. As recommended
by the DMC, the objective of its extent is to embrace South America and a good portion of the Pacific
Ocean, which allows capturing phenomena approaching the continent beforehand. The first nested
domain (d02, 18 km) embraced Chilean and Argentinian central regions and the Pacific Ocean near the
Atmosphere 2018, 9, 304 8 of 23

shore. The innermost domain (d03, 6 km), the main focus of this study, covered a reduced portion of
the Pacific Ocean, central Chile and the Andes mountains.

Figure 3. WRF model domains for all the simulations. Domain d01 correspond to the whole plot,
and the white point is over Santiago, Chile.

WRF can numerically diverge with high elevation gradients, resulting in anomalous vertical
wind speeds, with gradients of ∼300 m per 1 km [7]. Moreover, due to the model grid resolution,
the complex topography produced a considerable overestimation of Apoquindo station’s elevation
(+314 m). Following the approach by Carvalho et al. [57], the Apoquindo station location in WRF was
moved within a 5 km radius. The best output location to represent the station was 5 km further north
of its real location, which produced a height underestimation of 200 m, but a better performance of the
rainfall temporal distribution.
The prediction of the microphysics and PBL processes, and thus the rainfall, is expected to
improve with a finer vertical grid resolution. However, doubling the number of vertical levels from 31
to 62 did not enhance quantitative precipitation forecasts in the central US [58]. Based on previous
configurations [59,60], Saide et al. [39] proposed an optimal density of 39 vertical levels in the study
region, with a first layer at 10 m and six levels below 100 m, which allowed the best forecasting of wind
speed, temperature and chemical compound concentrations. This discretization was compared against
one with 50 vertical levels currently used by DMC. This finer discretization reproduced temperature
profiles and the saturated atmosphere near the surface more accurately, and thus it was chosen as the
default discretization for the following simulations.

2.4.2. Physics Schemes and Land Use

Short and long wave radiation schemes determine radiative fluxes in the simulation. The Dudhia
short wave scheme [61] and the RRTM long wave radiation scheme [62] were chosen. On the other
hand, cumulus schemes simulate the sub-grid processes related to convective clouds. Given the WRF
capacity to explicitly solve this, a parameterization is not needed in the inner domain. In the external
domains, the Grell 3D Ensemble [63] was used because of its acceptable performance when simulating
convective rainfall over Santiago [39].
Atmosphere 2018, 9, 304 9 of 23

WRF’s Land Use (LU) for terrain characterization uses different categories to depict landscapes,
crops, vegetation, forests, and urban areas. A comparison among LU datasets shows that the USGS
(U.S. Geological Survey) underestimates the extension of Santiago area, probably due to the use of old
maps . Hence, MODIS data were chosen for further simulations.
Land Surface models (LSM) calculate heat and moisture fluxes above the land, sea and ice cover.
The simplest LSM physics option that considers a five-layer model for thermal diffusion was not
considered, as it neither includes vegetation effects, nor the changes in snow cover, or soil moisture
over time. Soil moisture is the most significant part of flux exchanges between the surface and the
first level of the model, and also a crucial factor that affects near surface temperature and wind [15].
Noah-MP LSM was finally chosen given the improvements to the simulation of surface fluxes, timing
of snow water equivalent and runoff peaks [64].
The PBL scheme determines surface heat and moisture fluxes due to eddy transports in
the remaining volume of the atmospheric column over the terrain. In this study, the Mellor–
Yamada–Nakanishi–Niino (MYNN) Level 2.5 parameterization was used as a local approach with
total kinetic energy (TKE) closure [65], providing a good representation of good vertical profiles,
which in turn produces a closer agreement with the planet boundary layer shape, magnitude and
maximum values. In the Santiago area, the MYNN scheme represented accurately the diurnal PBL
height [39]. PBL was not tested in this study because WRF simulations of precipitation have been
shown to weakly depend on land surface, PBL and radiation parameterizations for winter storms over
complex topography [66].
The final WRF configuration is summarized in Table 2, for which 50 vertical levels of vertical
resolution are applied together with the GFS dataset.

Table 2. Final WRF simulation schemes.

Physical Scheme Parametrization

Short-wave radiation Dudhia
Long-wave radiation RRTM
Cumulus Grell 3D Ensemble
Planet Boundary Layer MMYN 2.5
Soil Layer MMYN
Land Surface Model Noah-MP

2.4.3. Microphysics
Microphysics (MP) schemes explicitly resolve processes of water, cloud and precipitation, and
their mixed-phases (i.e., ice–water interaction). The Lin et al. scheme (hereafter referred to as LIN)
employs six forms of water (water vapor, cloud water, cloud ice, rain, snow and graupel). Moreover,
this scheme allows the explicit inclusion of snow, and the correct simulation of changes from cloud
ice to snow and then to graupel. In general, LIN better describes the dynamics of the clouds inner
processes [67], being appropriated for high resolution simulations.
The WRF-Single-Moment-Microphysics scheme (WSM) varies according to the class,
i.e., the number of prognostic water substance variables. The WSM 3-Class (hereafter referred to
as WSM3), used in the DMC initial forecast, contains water vapor, cloud water/ice, and rain/snow.
In this scheme, rain and water occur above the freezing level, and snow and ice below it. On the other
hand, the WSM 6-Class (hereafter referred to as WSM6) contains mixing ratios of water vapor, cloud
water, cloud ice, snow, rain and graupel, and thus is more appropriate for high resolution simulations.
Although both WSM schemes have similar behaviors when simulating extreme rainfalls at low grid
resolutions (∼45 km), total rainfall and maximum intensity are better simulated in WSM6 at finer
resolutions (∼5 km) [68]. Moreover, for extreme rainfalls, the combined effect of microphysics and ice
sedimentation (available in LIN and WSM6) provides a better representation of cloud covered areas,
Atmosphere 2018, 9, 304 10 of 23

mean temperatures in the upper troposphere and surface rainfall amounts [69]. Nonetheless, in this
study, we tested the LIN, WSM3 and WSM6 schemes.

2.4.4. Lead Time

Forecast accuracy depends not only on the parameterization of physical schemes and grid
resolution, but also on the lead time. Commonly, forecasts get worse with increasing lead time
because synoptic characteristics are constantly changing. A smaller outer domain and shorter lead
times allow more accurate rainfall simulations, whereas larger domain sizes and lead times can increase
inner variabilities in WRF, displacing the spatial rainfall band and affecting rainfall forecasting [70].
Although acceptable precipitation forecasting can be obtained with lead times up to six days, values
of three days or less produce much better results [29]. In this study, lead times of 3, 4 and 5 days
(72, 96 and 120 h) previous to the day with the maximum rainfall intensity were analyzed.
In summary, the final WRF configuration depicted in Table 2 will be used to test three MP
parameterizations and three lead times, as described in Table 3. Due to computational resources
constraints, all the MP parameterizations are considered for the APR16 event, but only the LIN and
WSM6 parameterizations will be used for the OCT15 and MAY17 events. Because the more detailed
depiction of water phases is associated with these last two parameterizations, their comparison fulfills
the objective of finding the best MP scheme to forecast the selected rainfall events.

Table 3. Final WRF simulations with different lead times and MP schemes.

Rainfall Event MP Schemes Simulation Beginning (00:00) Lead Times (h) Simulations
OCT15 LIN & WSM6 15, 16 & 17 October 2015 72, 96 & 120 6
APR16 LIN, WSM3 & WSM6 12, 13 & 14 April 2016 72, 96 & 120 9
MAY17 LIN & WSM6 6, 7 & 8 May 2017 72, 96 & 120 6

2.5. Model Validation

Observed data from the Apoquindo and other stations located in Santiago valley and Piedmont
(Figure 1) were used to assess the quality of the meteorological simulations. Note that this assessment
excludes the first day of all the simulations, used for warming up the model (spin-up period). Given the
difference in the starting date of the simulations (lead time), the valid temporal range was used to
compare them with observed data. Thus, we try to remove the bias to be able to compare each of the
simulation’s performances. Simulated temperature and rainfall time series were assessed via the mean
absolute error (MAE), a robust metric that prevents single events from having a large impact on the
statistic [27].

N
1
MAE =
N ∑ | S t − Ot | , (1)
t =1

where N is the total number of data, and St and Ot are the simulated and observed meteorological data.
In addition, we also use the concept of error, which corresponds to the difference between St and Ot at
any time t. Finally, the simulations’ performances were also characterized through the Nash–Sutcliffe
Efficiency coefficient (NSE):
2
∑tN=1 (St − Ot )
NSE = 1 − 2
. (2)
∑tN=1 (Ot − Ōt )
NSE has a range [−∞, 1], where values larger than 0 imply a prediction better than the average of
the observations (Ō). Since the NSE is extremely sensitive to outlier data, a modified version of the
Index of Agreement (IoA) [71] is also used. IoA is a less sensitive to outliers metric that ranges [−1, 1].
Values bigger than 0.8 are considered to indicate a good performance of the model:
Atmosphere 2018, 9, 304 11 of 23

∑tN=1 |St − Ot |

if ∑tN=1 |St − Ot | ≤ 2 ∑tN=1 Ot − Ōt ,
 
1 −


2 ∑tN=1 Ot − Ōt 
  
IoA = (3)
 2 ∑tN=1 Ot − Ōt 
− 1, in other case.


 N
∑ t =1 | S t − Ot |
Finally, we adopted the tolerance criteria for absolute error, defined in Table 4, used by DMC.

Table 4. Absolute difference tolerance criteria for simulated variables used by DMC.

Meteorological Variable Absolute Difference Tolerance Criteria

Dew temperature 1 ◦ C in surface, 2 ◦ C in atmosphere
Temperature 1 ◦ C in surface, 2 ◦ C in atmosphere
Wind speed 2.57 m/s (∼5 knots) for all data
Wind direction 20◦ in surface, 15◦ in pressure levels above the 850 hPa

Non-zero hourly simulated rainfall pulses were assumed to be those larger than 1 mm, while an
inter-event arrival time (IEAT) (i.e., the minimum dry time between two independent rainfall events)
of 30 h recommended for Santiago [72] was adopted. To assess the temporal distribution of the most
intense portion of the rainfall event, the rainiest N-consecutive hours (NRH) were computed, where N
is an integer. This approach implies finding the N-consecutive hourly pulses with the largest amount
of rainfall, with N integer from 1 (the rainiest hour) to 10 (the 10 consecutive hours with the largest
amount of cumulative precipitation). Thus, it is possible to assess the simulation in terms of temporal
distribution of the most intense rainfall pulses, which largely affect the basin hydrologic response.
Temperature time series associated with each NRH (T NRH ), were analyzed via MAE, emphasizing
both precipitation and temperature when N = 5 h (i.e., NRH(5) and T NRH (5) ). This time corresponds
to the estimated time of concentration of Quebrada de Ramón, a metric that is representative of the
hydrological response time of the catchment.

3. Results

3.1. Local Conditions

To assess the local performance of WRF over the rainfall events, we first characterized the overall
rainfall orographical gradient using total rainfall amounts registered in all the meteorological stations
(Table 1). Figure 4 compares total simulated rainfall amounts (considered to be the mean of the values
simulated using the three lead times for every MP option) against the observed ones. Different markers
according to the elevation are used in the figure. The left panels compare in every meteorological
station simulated and observed data, and the dashed line corresponds to the 1:1 ratio. The right panels’
plot rainfall against height, the scatter in the plot represents simulated values for several locations
in Santiago city and the west hill slope of the Andean foothills, to give a spatial context for the WRF
rainfall distribution.
For the OCT15 rainfall event (Figure 4, upper row), WRF overestimate rainfall amounts,
particularly at higher elevations. The orographical gradient inside Santiago city is accurately simulated,
but the scatter plot shows more dispersion at elevations higher than 700 m.a.s.l. (i.e., the mountainous
area of Santiago). The LIN MP option produces a bigger bias than the WSM6 option for all the points
belonging to the Santiago basin, although this bias is reduced for lower elevations. In general, for the
OCT15 event, WRF predicts slightly more rainfall as the latitude increases southwards, for the points
inside Santiago, excluding the mountainous area. Overall, the WSM6 scheme provides better results
(i.e., less biased) with a 5% significance level (Mann–Kendall test [44,45]).
Atmosphere 2018, 9, 304 12 of 23

Figure 4. Total observed and simulated precipitation according to height (left) and orographical
gradient (right) for OCT15 (upper row), APR16 (central row) and MAY17 (lower row) events. The
mean value of all lead times (72, 96 and 120 h) and three MP parameterizations (LIN, WSM6 and WSM3
when used) are considered. In right panels, observed data are plotted in black squares.

In the APR16 event (Figure 4, middle row), WRF overestimate rainfall amounts for all the stations
located above 800 m.a.s.l., corresponding to Santiago’s foothills. Predictions for lower elevations
are much better with small underestimations of the observed precipitation. WRF simulates a strong
correlation between total precipitation and height, although the observed rainfall ranges between
80–120 mm. In the orographical gradient, rainfall was understimated in the west part of Santiago
(lower elevation) and overestimated in the foothills. The scattered pattern appears once more, with
more dispersion (and thus bigger biases) for the foothill stations located in the southern area of Santiago
(AN, QM, and SJ stations). This also happened for the Apoquindo station, where a ∼100 mm bias
was simulated. The bias reduces for gauges located at northern latitudes below 1000 m.a.s.l. (gauges
Atmosphere 2018, 9, 304 13 of 23

CC and TO). Overall, the LIN scheme produced larger biases, whereas WSM3 and WSM6 perform
similarly depending on the station. In fact, WSM3 performs better in the mountainous zone of Santiago
(CC, AN, QM and AP), while WSM6 performs better in QN, PL and SJ (southwest Santiago valley and
the mountainous area in the East, respectively). The biases are considerably larger in the mountainous
zone, so no real advantage can be addressed over the WSM3 scheme.
Finally, for the MAY17 event (Figure 4, bottom row) rainfall is overestimated particularly at higher
elevations, although the rainfall in the highest two stations was understimated by WRF. The smooth
orographical gradient in the recorded data was captured in the WRF simulation. In Apoquindo station,
both MP options (LIN and WSM6) produced similar bias (∼20 mm). Once more, the WSM6 scheme
performed generally better as it produced smaller biases.

3.2. Rainfall and Temperature Simulation

In the OCT15 storm, simulations with 72 and 96 lead time could not simulate the rain pulse
observed during 17 October. As a countermeasure, a base amount of precipitation (10 mm), equivalent
to the registered precipitation until 17 October, was considered in Figure 5, to establish comparisons
on even ground. The frontal system conditions of the OCT15 event allowed a more accurate prediction
of the temporal distribution of the rainfall (IoA = 0.64 and for the mean of LIN simulations, and 0.67
for the mean of WSM6 simulations), since most MP schemes captured the initial rainfall (19 October).
Simulations with 120 h of lead time also captured a previous pulse (17 October), as seen in Figure 5
(upper row). For shorter lead times, both MP schemes not only predicted the occurrence of another
precipitation pulse (20 October), which is inexistent in observed data, but compensates the underrated
rainfall pulse of 19 October. Given this additional pulse, the simulations with lead time of 72 h have
a bias of less than 2 mm in the final hours of the rainfall event.
None of the simulations could accurately predict the APR16 rainfall event (Figure 5, left central
plot). There was an initial ∼48 h lag between the observed and simulated rainfall peaks, as WRF
predicts rainfall earlier. Thus, the rainfall beginning is considered to be 16 April, ignoring the WRF
first rainfall pulse of 15 April by subtracting the cumulative rainfall until that date. This improved the
model performance, as a smaller MAE (∼1 mm) and a larger mean IoA value are obtained for all MP
options. Nevertheless, some values of IoA are still negative. Furthermore, for simulations with shorter
lead times, the rainfall pulse was correctly predicted (17 April), especially with a lead time of 72 h.
Overall, the hourly forecast of the last storm hours was slightly more accurate (∼−20 mm bias). As the
OCT15 and MAY17 simulated events don’t show these errors, they are not explained by the complex
topography, but by the rainfall properties which seem to be maximized in the initial GFS dataset. This
error in the rainfall amount was observed in all the meteorological stations of Table 1 (not shown).
Indeed, the APR16 event had unique meteorological characteristics and a large magnitude, which
could partially explain the lag errors in the prediction.
In the MAY17 storm, the 72 h lead time simulations start after the rain pulse registered during
7 May. Therefore, the registered precipitation until that date was added as a base amount in Figure 5.
The same as before, this was made to establish even comparisons. For the MAY17 rainfall simulations,
there is an uneven ensemble performance for both MP options (LIN and WSM6 IoA ∼0.58), as the
total rainfall amount is underestimated for all lead times (Figure 5, bottom row). However, most
simulations capture the rainfall pulse of 7 May, with the exception of LIN72 (and WSM672 to a
lesser extent) probably because its initialization was about 9 h before the initial pulse, within the
model spin-up time period. Although all simulations capture the 10 May rainfall pulse, errors of
∼−40 mm obtained with lead times of 96 are not produced by the simulation with lead time of 120 or
72 h. For these two simulations, WRF simulated intermittent rainfall with lower intensities during
10–12 May for 120 h lead time, and 10–13 May for 72 h lead time. The observed intensities were higher
during the first half of 10 May, but the WRF model could not capture them. The LIN scheme produces
the closest steep slope in agreement with the observed data (the last hours of 9 May) for a lead time of
120 h, but the WSM6 scheme produces the same agreement with a 120 h or 72 h lead time.
Atmosphere 2018, 9, 304 14 of 23

For temperature time series, there is an accurate performance for the whole ensemble, since the
temperature dynamics was captured for the OCT15 and MAY17 events with minor errors. The mean
performance between both schemes was very similar, with IoA = 0.77 for the OCT15 event and IoA
= 0.69 for the MAY17 event. Steep temperature increments are well captured by all simulations,
but overlapping and disagreement appear just before the rain started (∼1 day, from the last hours
of 18 October for OCT15, and from 10 May for MAY17, and ∼3 days, from 14 April for APR16),
where a general temperature performance decrease takes place. Temperatures are overestimated in
OCT15 and MAY17 events, just before the storm peak (19 April for APR16 and 10 May for MAY17).
A constant underestimation between 14–17 April can be seen in APR16, where all schemes had a
similar performance (IoA ∼0.47), probably due to the complex meteorological configuration of the
event (Figure 5, right panels).

Figure 5. Observed and simulated precipitation (left) and temperature (right) time series for OCT15
(upper row), APR16 (central row) and MAY17 (lower row) events. Three lead times (72, 96 and 120 h)
and three MP parameterizations (LIN, WSM6 and WSM3 when used) are considered. Graphical
modifications to the rainfall time series are indicated in the previous paragraphs.

3.3. Rainfall Forecast Performance

In addition to the total rainfall amount and dynamics, another relevant attribute of precipitation
events is the timing of the rainiest pulses. A good forecasting of this dynamics is essential for modeling
and predicting peak flows in a catchment of a given size. To visualize the performance of the forecasting
on this regard, we plot the observed and simulated N rainiest hours (NRH) for all the different MP
parameterizations and lead times, with N ranging between 1 and 10 h (Figure 6). Overall, the wide
range of simulated curves contains the observed NRH curve, although major differences between the
observed and simulated curves for the APR16 event were observed. NRH(5) values are overstimated
by all the ensemble. Unfortunately, the parameterization that best reproduces this value (WSM3120 ) has
Atmosphere 2018, 9, 304 15 of 23

the largest associated MAE value (Figure 6, central row). Generally, the MAE of the T NRH ensemble
values were above the tolerance criterion shown in Table 4.

Figure 6. Observed and simulated N rainiest hours (NRH, left) and T NRH MAE (right) from 1 to 10 h,
for OCT15 (upper row), APR16 (central row) and MAY17 (lower row) simulations. The vertical black
dotted line represents the Quebrada de Ramón time of concentration. Three lead times (72, 96 and
120 h) and three MP parameterizations (LIN, WSM6 and WSM3 when used) are considered.

For the OCT15 rainfall event, a wide range of NRH(5) values is simulated, with all of them
underestimating the observed value (22 mm). Nevertheless, the values provided by the LIN scheme
are closer to the observed data than those from the WSM6 scheme. In general, a smaller lead time
provides values closer to observations. As an ensemble, the average NRH(5) value associated with
all the simulations underrates in ∼10 mm the observed data. For the MAY17 event, all of the NRH(5)
values underestimate the observation, except the value from the LIN scheme with a 120 h lead time.
Overall, mean WSM6 values for all the lead times are closer to the observed data ( ∼−10 mm bias).
For the LIN scheme, smaller lead times are related with worst performance.
Among the used MP schemes, WSM6 shows more consistency, regardless of the lead time,
and the errors are within the tolerance criterion for temperature (Table 4). LIN performs similarly,
but results are less consistent and become worse with shorter lead times. Overall, the average WSM3
parameterization is the worst in reproducing NRH(5), while its relation with the lead time is less
consistent than the LIN and WSM6 options for NRH(5) and T NRH (5) MAE (Figure 6).
Total rainfall is overestimated for the APR16 event, whereas, for the other rainfall events, this
amount is within the range of simulated values. These values tend to be better with shorter lead times,
although no clear trend is observed. Furthermore, the average of the WSM6 simulations is worse than
Atmosphere 2018, 9, 304 16 of 23

that of the LIN scheme in reproducing the observed total amount. On the other hand, the rainfall
length is reasonably predicted by the ensemble for all the events, with WSM6 performing better than
the LIN parameterization, although real values are not within the 25th–75th percentiles. Overall, and
because it predicts the NRH(5) values, WSM6 is considered to be a good MP option to choose for
rainfall forecasting in the study area, and eventually other front-range catchments nearby.

3.4. Freezing Level Height

The freezing level height is a relevant factor in the hydrological response of the catchment because
it defines the area which receives rainfall instead of snowfall. We assessed the prediction of the freezing
level height for different MP schemes and lead times, in order to find any advantage in a certain
combination that could be useful under a flash flood warning scenario.
Freezing level height is quite well simulated by all the ensemble for the OCT15 and MAY17
events. Less differences between simulated and observed data are obtained during freezing level
height declines. On the other hand, when the freezing level height increases, the ensembles are more
discrepant. Regardless of the behavior of this level, during the rainfall peak, the error tends to be
negligible (Figure 7). For APR16, the decline of the freezing level is overestimated, especially for longer
lead times. This illustrates how the atmosphere stability development during cutoff low events is
better captured using small lead times.

Figure 7. Freezing level height simulated with lead times (LT) of five (left), four (center) and three
days (right), for OCT15 (upper row), APR16 (central row) and MAY17 (lower row) events.
Atmosphere 2018, 9, 304 17 of 23

Shorter lead times improve the forecast of the OCT15 and MAY17 events, better capturing the
freezing level increase ∼2 days before the peak intensity (Figure 7, right panels). Longer lead times
could not accurately simulate this peak value. A lead time of 72 h (3 days) was optimal for capturing
the freezing level height development. Even for APR16, where there’s no real agreement between
simulated and real data, this lead time produces relatively small error values. No clear advantage was
found in relation to the MP schemes, which is expected given the complex interactions that define
the freezing level height, more related to the PBL schemes and the vertical grid resolution. In that
sense, a broader study incorporating a full test of the PBL scheme and several vertical densities are
recommended for the future.

3.5. Ensemble Performance

To understand the ensemble performance, we plot the observed hyetograph in the Apoquindo
gauge (blue bars) and temperature series against the red region denoting the time at which more than
50% of the WRF ensemble predict rainfall (Figure 8). In the figure, we also plotted the 25th, 50th and
75th percentiles of the temperatures simulated by the ensemble. Note that the ensemble corresponds
to as all the simulations performed for each rainfall event, indistinctive from their lead time or MP
scheme, i.e., 6 simulations for OCT15, 9 for APR16 and 6 for MAY17 (Table 3).
The occurrence of the main pulses of the OCT15 event is very well predicted, although a final
pulse that did not happen in 20 October is produced. At the time of the rainfall peak, the ensemble
simulates a temperature ∼2 ◦ C higher than the observed one, due to the abrupt temperature decrement
prior to the rainfall beginning.
On the other hand, the WRF ensemble fails in identifying the temporal distribution of the more
intense APR16 rainfall event with a 48 h delay in the prediction of the occurrence (Figure 8, central
plot). However, a third rainfall pulse predicted by the ensemble matches the observed rainfall, which
can lead to wrong conclusions. As mentioned by Garreaud [53], the warm storms last from 12 to 36 h,
coupled with an air temperature drop (>3 ◦ C) and the highest precipitation within the first hours of
the event. This conditions indeed occurred during the APR16 event, particularly the pulse starting
on 17 April. The ensemble underestimates the observed temperature for ∼2 days before the rainfall
initiates, but, during its peak, this temperature was contained between the 25th and 75th percentiles.
The ensemble performance for the MAY17 event was more difficult to analyze, due to the frontal
characteristics of the event. Although it was unable to accurately predict the characteristics of the
rainfall occurrence, the ensemble clearly simulated three different instances where precipitation
happens. Rain is not properly predicted in the main period of intensive rainfall (11 May), with only
one of the three larger intensities being matched by a single red stripe. A final rainfall pulse (12 May)
is simulated by the model with almost half a day of anticipation. Temperature is accurately predicted
during the first rainfall event peak, but, for the main peak, there is a ∼2 ◦ C overestimation.
Overall, the WRF ability to simulate frontal rainfalls is acceptable with a proper MP scheme
parametrization, such as LIN or WSM6. In addition, the accurate prediction of rainfall is still too
complicated over the complex topography of Santiago.
Atmosphere 2018, 9, 304 18 of 23

Figure 8. Ensemble performance for OCT15 (upper plot), APR16 (central plot) and MAY17 (lower plot)
events. Red stripes indicate WRF rain occurrence prediction when the ensemble probability of rain
exceeds 0.5, and blue bars are the observed hourly precipitation in Apoquindo rain gauge. Temperature
time series includes 25 (blue), 50th (green) and 75 (magenta) percentiles for the ensemble forecast
against observed data (black squares).
Atmosphere 2018, 9, 304 19 of 23

4. Conclusions
Rainfall forecast over complex topography using WRF was studied through the simulation
of three events between 2015 and 2017 in Quebrada de Ramón, a 38 km2 mountainous Andean
watershed in central Chile. Three microphysics (MP) parametrizations (i.e., LIN, WSM6 and WSM3,
currently used by the Chilean Weather Agency) were tested to find an optimum model performance.
The simulations considered an horizontal resolution of 6 km and 50 vertical levels for improving
atmospheric temperature profiles, as well as 0.5◦ grid resolution GFS dataset. A realistic representation
of Santiago’s urban area was provided by MODIS, and Noah-MP model was used as the land surface
model. Finally, variable lead times of 72, 96 and 120 h before the rainfall start were also analyzed.
The uneven prediction of rainfall length and total amount by the LIN and WSM6
parameterizations tends to improve with shorter lead times. Both performed better than WSM3,
mainly for temperatures and rainfall intensities. The N-rainiest consecutive hours (NRH), a relevant
characteristic of storm events given the impact on the hydrologic response, could not be well predicted
by any of the MP options for cutoff low events. Frontal system events were better captured, mostly via
the WSM6 MP option, which also predicted well the concurrent temperatures with high intensities.
This is very relevant, as high flows and potential floods are typically produced by warm events when
high temperatures and rainfall take place at the same time [53].
The temporal evolution of precipitation, temperature and freezing level height were properly
predicted for shorter lead times, especially for frontal system events, while complex meteorological
cutoff low characteristics lead to poor forecasts. No clear trends in lead times were found, but shorter
values (72 h ahead rainfall event) tended to provide more accurate simulations. Further testing on PBL
schemes and vertical resolution is needed in this theme.
From our results, the WSM6 scheme resulted to be the best to simulate rainfall events in the
Andean watershed under study. Nevertheless, rainfall simulation in WRF over complex topography is
still a challenging issue, and its ability to accurately simulate rainfall, particularly non-frontal, events
over Chilean central mountainous and foothills areas, where some of the main cities are located, is far
from ideal. Further investigation should focus on simulating more rainfall events for which observed
data could be available, as well as testing additional microphysical schemes such as Thompson
scheme [73] or the Aerosol-aware Thompson scheme [74].
Finally, it would be of interest to improve rainfall forecasts by combining NWP tools with
geostatistical methods, which incorporate spatial statistics and local topography to describe, enhance
or predict rainfall in basins with complex topography [75–77]. Such geostatistical approach would
eventually offset to a certain extent WRF forecast problems over the Andean complex topography.

Author Contributions: G.Y.-M. and J.G. conceived the study, R.D. performed the WRF simulations; G.Y.-M., J.G.,
M.C. and R.G. analyzed the data; G.Y.-M. wrote the paper; and all the authors revised the paper.
Funding: The authors acknowledge funding from the following projects: ND-PUC Seed Fund, FONDECYT
1130378, FONDECYT 1161439, (CR)2/FONDAP 15110009, CIGIDEN/FONDAP 15110017 and CEDEUS/FONDAP
15110020.
Acknowledgments: The corresponding author is grateful for the Arturo Cousiño Lyon and Magíster Nacional
CONICYT scholarships, the Chilean Weather Agency computational resources and collaboration, and for the
constructive and valuable comments from Dr. Christian González, Ps. Mateo Zamorano, and two anonymous
reviewers.
Conflicts of Interest: The authors declare no conflict of interest.
Atmosphere 2018, 9, 304 20 of 23

Abbreviations
The following abbreviations are used in this manuscript:

ARW Advanced Research WRF

DMC Chilean Weather Agency
FNL Final Operational Global Analysis data
GFS Global Florecast System
LSM Land Surface Model
MODIS Moderate Resolution Imaging Spectroradiometer
MP Microphysics
NCAR National Center for Atmospheric Research
NCEP National Centers for Environmental Prediction
NRH N-rainiest consecutive hours
NWP Numerical Weather Prediction
PBL Planet Boundary Layer
RRTM Rapid Radiative Transfer Model
WRF Weather Research and Forecasting

References
1. EM-DAT. The International Disaster Database; Technical Report; Centre for Research on the Epidemiology of
Disasters (CRED), School of Public Health, Université Catolique de Louvain: Louvain, Bélgica, 2016.
2. Alfieri, L.; Bisselink, B.; Dottori, F.; Naumann, G.; Roo, A.; Salamon, P.; Wyser, K.; Feyen, L. Global projections
of river flood risk in a warmer world. Earth’s Future 2016, 5, 171–182. [CrossRef]
3. Kundzewicz, Z.W.; Kanae, S.; Seneviratne, S.I.; Handmer, J.; Nicholls, N.; Peduzzi, P.; Mechler, R.;
Bouwer, L.M.; Arnell, N.; Mach, K.; et al. Flood risk and climate change: Global and regional perspectives.
Hydrol. Sci. J. 2014, 59, 1–28. [CrossRef]
4. Bozkurt, D.; Rondanelli, R.; Garreaud, R.; Arriagada, A. Impact of warmer Eastern Tropical Pacific SST on
the March 2015 Atacama Floods. Mon. Weather Rev. 2016, 144, 4441–4460. [CrossRef]
5. Wilcox, A.C.; Escauriaza, C.; Agredano, R.; Mignot, E.; Zuazo, V.; Otárola, S.; Castro, L.; Gironás, J.;
Cienfuegos, R.; Mao, L. An integrated analysis of the March 2015 Atacama floods. Geophys. Res. Lett. 2016,
43, 8035–8043. [CrossRef]
6. Directors Guild of America. Caracterización de Suelos y Generación de Información Meteorológica Para
Prevención de Riesgos Hidrometeorológicos Cuencas Salado y Copiapó; Technical Report; Dirección General
de Aguas—Modelación Ambiental SPA; Directors Guild of America: Sunset Boulevard, CA, USA, 2016.
7. Otarola, S.; Dimitrova, R.; Leo, L.; Alcafuz, R.; Escauriaza, C.; Arroyo, R.; nez Morroni, G.Y.; Fernando, H.J.
On the role of the Andes on weather patterns and related environmental hazards in Chile. In Proceedings of
the 19th Joint Conference on the Applications of Air Pollution Meteorology with the AWMA, New Orleans,
LA, USA, 10–11 January 2016.
8. Arnold, D.; Morton, D.; Schicker, I.; Seibert, P.; Rotach, M.W.; Horvarth, K.; Dudhia, J.; Satomura, T.; Muller,
M.; Zangl, G.; Takemi, T.; Serafin, S.; Schmidli, J.; Schneider, S. Issues in high-resolution atmospheric
modeling in complex topography - The HiRCoT workshop. Croat. Meteorol. J. 2012, 47, 3–11.
9. Goger, B.; Rotach, M.W.; Gohm, A.; Stiperski, I.; Fuhrer, O. Current challenges for numerical weather
prediction in complex terrain: Topography representation and parameterizations. In Proceedings of the
2016 International Conference on High Performance Computing & Simulation (HPCS), Innsbruck, Austria,
18–22 July 2016; pp. 890–894. [CrossRef]
10. Bongioannini Cerlini, P.; Emanuel, K.A.; Todini, E. Orographic effects on convective precipitation and
space-time rainfall variability: Preliminary results. Hydrol. Earth Syst. Sci. 2005, 9, 285–299. [CrossRef]
11. Houze, R.A. Orographic effects on precipitating clouds. Rev. Geophys. 2012, 50. [CrossRef]
12. Jiménez, P.A.; Dudhia, J. Improving the representation of resolved and unresolved topographic effects on
surface wind in the WRF model. J. Appl. Meteorol. Climatol. 2012, 51, 300–316. [CrossRef]
13. Lorente-Plazas, R.; Jiménez, P.A.; Dudhia, J.; Montávez, J.P. Evaluating and improving the impact of
the atmospheric stability and orography on surface winds in the WRF mode. Mon. Weather Rev. 2016,
144, 2685–2693. [CrossRef]
Atmosphere 2018, 9, 304 21 of 23

14. Madala, S.; Satyanarayana, A.N.V.; Srinivas, C.V.; Boadh, R. Sensitivity of PBL Schemes of WRF-ARW Model in
Simulating Mesoscale Atmospheric Flow-Field Parameters Over a Complex Terrain; AGU Fall Meeting Abstracts;
American Geophysical Union: Washington, DC, USA, 2014; pp. A43A–3242.
15. Dimitrova, R.; Silver, Z.; Zsedrovits, T.; Hocut, C.M.; Leo, L.S.; Sabatino, S.D.; Fernando, H.J. Assessment
of planetary boundary-layer schemes in the Weather Research and Forecasting mesoscale model using
MATERHORN field data. Bound.-Layer Meteorol. 2016, 159, 589–609. [CrossRef]
16. Siuta, D.; West, G.; Stull, R. WRF hub-height wind forecast sensitivity to PBL scheme, grid length, and initial
condition choice in complex terrain. Weather Forecast. 2017, 32, 493–509. [CrossRef]
17. Pontoppidan, M.; Reuder, J.; Mayer, S.; Kolstad, E.W. Downscaling an intense precipitation event in complex
terrain: The importance of high grid resolution. Tellus A Dyn. Meteorol. Oceanogr. 2017, 69, 1271561.
[CrossRef]
18. Mass, C.F.; Ovens, D.; Westrick, K.; Colle, B.A. Does increasing horizontal resolution produce more skillful
forecasts? Bull. Am. Meteorol. Soc. 2002, 83, 407–430. [CrossRef]
19. Gneiting, T.; Raftery, A.E. Weather forecasting with ensemble methods. Science 2005, 310, 248–249. [CrossRef]
[PubMed]
20. Zhang, H.; Pu, Z. Beating the uncertainties: Ensemble forecasting and ensemble-based data assimilation in
modern Numerical Weather Prediction. Adv. Meteorol. 2010, 2010, 432160. [CrossRef]
21. Lee, J.A.; Kolczynski, W.C.; McCandless, T.C.; Haupt, S.E. An objective methodology for configuring and
down-selecting an NWP ensemble for low-level wind prediction. Mon. Weather Rev. 2012, 140, 2270–2286.
[CrossRef]
22. Bauer, P.; Thorpe, A.; Brunet, G. The quiet revolution of numerical weather prediction. Nature 2015,
525, 47–55. [CrossRef] [PubMed]
23. Ruiz, J.J.; Saulo, C.; Nogués-Paegle, J. WRF Model sensitivity to choice of parameterization over South
America: Validation against surface variables. Mon. Weather Rev. 2010, 138, 3342–3355. [CrossRef]
24. Evans, J.P.; Ekström, M.; Ji, F. Evaluating the performance of a WRF physics ensemble over South-East
Australia. Clim. Dyn. 2012, 39, 1241–1258. [CrossRef]
25. Kim, J.H.; Shin, D.B.; Kummerow, C. Impacts of a priori databases using six WRF microphysics schemes on
passive microwave rainfall retrievals. J. Atmos. Ocean. Technol. 2013, 30, 2367–2381. [CrossRef]
26. Katragkou, E.; García-Díez, M.; Vautard, R.; Sobolowski, S.; Zanis, P.; Alexandri, G.; Cardoso, R.M.;
Colette, A.; Fernandez, J.; Gobiet, A.; et al. Regional climate hindcast simulations within EURO-CORDEX:
Evaluation of a WRF multi-physics ensemble. Geosci. Model Dev. 2015, 8, 603–618. [CrossRef]
27. Ekström, M. Metrics to identify meaningful downscaling skill in WRF simulations of intense rainfall events.
Environ. Model. Softw. 2016, 79, 267–284. [CrossRef]
28. Hirabayashi, S.; Kroll, C.N.; Nowak, D.J. Component-based development and sensitivity analyses of an air
pollutant dry deposition model. Environ. Model. Softw. 2011, 26, 804–816. [CrossRef]
29. Siddique, R.; Mejia, A.; Brown, J.; Reed, S.; Ahnert, P. Verification of precipitation forecasts from two
numerical weather prediction models in the middle Atlantic region of the USA: A precursory analysis to
hydrologic forecasting. J. Hydrol. 2015, 529, 1390–1406. [CrossRef]
30. Buizza, R.; Leutbecher, M. The forecast skill horizon. Q. J. R. Meteorol. Soc. 2015, 141, 3366–3382. [CrossRef]
31. Skamarock, W.C.; Klemp, J.B.; Dudhia, J.; Gill, D.O.; Barkerand, D.M.; Duda, M.G.; Huang, X.; Wang, W.;
Powers, J.G. A Description of the Advanced Research WRF Version 3; NCAR Technical Note; NCAR: Boulder,
CO, USA, 2008.
32. Jiménez, P.A.; González-Rouco, J.F.; García-Bustamante, E.; Navarro, J.; Montávez, J.P.; de Arellano, J.V.G.;
Dudhia, J.; Muñoz-Roldan, A. Surface wind regionalization over complex terrain: Evaluation and analysis
of a high-resolution WRF simulation. J. Appl. Meteorol. Climatol. 2010, 49, 268–287. [CrossRef]
33. Zhang, H.; Pu, Z.; Zhang, X. Examination of errors in near-surface temperature and wind from WRF
numerical simulations in regions of complex terrain. Weather Forecast. 2013, 28, 893–914. [CrossRef]
34. Karki, R.; ul Hasson, S.; Gerlitz, L.; Schickhoff, U.; Scholten, T.; Böhner, J. Quantifying the added value of
convection-permitting climate simulations in complex terrain: A systematic evaluation of WRF over the
Himalayas. Earth Syst. Dyn. 2017, 8, 507–528. doi:10.5194/esd-8-507-2017. [CrossRef]
35. Soltanzadeh, I.; Bonnardot, V.; Sturman, A.; Quénol, H.; Zawar-Reza, P. Assessment of the ARW-WRF
model over complex terrain: The case of the Stellenbosch Wine of Origin district of South Africa.
Theor. Appl. Climatol. 2017, 129, 1407–1427. [CrossRef]
Atmosphere 2018, 9, 304 22 of 23

36. Jiménez-Esteve, B.; Udina, M.; Soler, M.; Pepin, N.; Miró, J. Land use and topography influence in a complex
terrain area: A high resolution mesoscale modelling study over the Eastern Pyrenees using the WRF model.
Atmos. Res. 2018, 202, 49–62. doi:10.1016/j.atmosres.2017.11.012. [CrossRef]
37. Garreaud, R.; Falvey, M.; Montecinos, A. Orographic precipitation in coastal southern Chile: Mean
distribution, temporal variability, and linear contribution. J. Hydrometeorol. 2016, 17, 1185–1202.
doi:10.1175/JHM-D-15-0170.1. [CrossRef]
38. Puliafito, S.; Allende, D.; Mulena, C.; Cremades, P.; Lakkis, S. Evaluation of the WRF model configuration for
Zonda wind events in a complex terrain. Atmos. Environ. 2015, 166, 24–32. doi:10.1016/j.atmosres.2015.06.011.
[CrossRef]
39. Saide, P.E.; Carmichael, G.R.; Spak, S.N.; Gallardo, L.; Osses, A.E.; Mena-Carrasco, M.A.; Pagowski,
M. Forecasting urban PM10 and PM2.5 pollution episodes in very stable nocturnal conditions
and complex terrain using WRF-Chem CO tracer model. Atmos. Environ. 2011, 45, 2769–2780.
doi:10.1016/j.atmosenv.2011.02.001. [CrossRef]
40. Catalán, M. Acodicionamiento de un Modelo Hidrológico Geomorfológico Urbano de Onda Cinemática
para su Aplicación en Cuencas Naturales. Ph.D. Thesis, Pontificia Universidad Católica de Chile, Santiago,
Chile, 2013.
41. Garreaud, R.; Rutllant, J. Precipitación estival en los Andes de Chile central: Aspectos climatológicos.
Atmósfera 1997, 10, 191–211.
42. Vargas, X. Corrientes de detritos en la Quebrada de Macul, Chile. Estudio de caudales máximos. Ingeniería
del Agua 1999, 6, 341–344. [CrossRef]
43. Dirección General de Aeronáutica Civil (DGAC). Dirección Meteorológica de Chile (DMC). Available online:
http://www.meteochile.cl (accessed on 22 March 2018).
44. Mann, H.B. Nonparametric tests against trend. Econometrica 1945, 13, 245–259. [CrossRef]
45. Kendall, M.G. Rank Correlation Methods, 4th ed.; Griffin, Charles: London, UK, 1975; p. 272.
46. Garreaud, R.D.; Alvarez-Garreton, C.; Barichivich, J.; Boisier, J.P.; Christie, D.; Galleguillos, M.; LeQuesne,
C.; McPhee, J.; Zambrano-Bigiarini, M. The 2010–2015 megadrought in central Chile: Impacts on regional
hydroclimate and vegetation. Hydrol. Earth Syst. Sci. 2017, 21, 6307–6327. [CrossRef]
47. University of Wyoming. College of Engineering. Soundings. Available online: http://weather.uwyo.edu/
upperair/sounding.html (accessed on 22 March 2018).
48. Roney, J.A. Statistical wind analysis for near-space applications. J. Atmos. Sol.-Terr. Phys. 2007, 69, 1485–1501.
[CrossRef]
49. Fuenzalida, H.; Sanchez, R.; Garreaud, R. A climatology of cutoff lows in the Southern Hemisphere.
J. Geophys. Res. Atmos. 2005, 110. [CrossRef]
50. Garreaud, R.; Fuenzalida, H. The Influence of the Andes on Cutoff Lows: A Modeling Study.
Mon. Weather Rev. 2007, 135, 1596–1613. [CrossRef]
51. Aguilar, X. Análisis de Eventos Extremos de Precipitación en División Los Bronces, Cordillera Central de
Chile. Ph.D. Thesis, Tesis Para Optar al Título Profesional de Meteorólogo, Departamento de Meteorología,
Universidad de Valparaíso, Valparaíso, Chile, 2010.
52. MOP. Manual de Drenaje Urbano: Guía para el Diseño, Construcción, Operación y Conservación de Obras de Drenaje
Urbano; Technical Report; Ministerio de Obras Públicas: Santiago, Chile, 2013.
53. Garreaud, R. Warm winter storms in central Chile. J. Hydrometeorol. 2013, 14, 1515–1534. [CrossRef]
54. National Center for Atmospheric Research (NCAR). WRF Users Page. Available online: http://www.mmm.
ucar.edu/wrf/users (accessed on 22 March 2018).
55. National Oceanic and Atmospheric Administration (NOAA). National Centers for Environmental
Information. ModelData. Available online: https://www.ncdc.noaa.gov/data-access/model-data/
(accessed on 22 March 2018).
56. National Centers for Environmental Prediction (NCEP). NCEP FNL Operational Model Global Tropospheric
Analyses, Continuing from July 1999. Available online: http://rda.ucar.edu/datasets/ds083.2/ (accessed
on 22 March 2018).
57. Carvalho, D.; Rocha, A.; Gomez-Gesteira, M.; Santos, C. A sensitivity study of the WRF model in wind
simulation for an area of high wind energy. Environ. Model. Softw. 2012, 33, 23–34. [CrossRef]
58. Aligo, E.A.; Gallus, W.A.; Segal, M. On the impact of WRF model vertical grid resolution on midwest
summer rainfall forecasts. Weather Forecast. 2009, 24, 575–594. [CrossRef]
Atmosphere 2018, 9, 304 23 of 23

59. Seaman, N.; Gaudet, B.; Zielonka, J.; Stauffer, D. Sensitivity of vertical structure in the stable boundary layer
to variations of the WRF model’s Mellor Yamada Janjic turbulence scheme. In Proceedings of the 9th WRF
Users Workshop, Boulder, CO, USA, 23–27 June 2008.
60. Rahn, D.A.; Garreaud, R. Marine boundary layer over the subtropical southeast Pacific during VOCALS-REx.
Part 1: Mean structure and diurnal cycle. Atmos. Chem. Phys. 2010, 10, 4491–4506. [CrossRef]
61. Dudhia, J. Numerical study of convection observed during the winter monsoon experiment using a
mesoscale two-dimensional model. J. Atmos. Sci. 1989, 46, 3077–3107. [CrossRef]
62. Mlawer, E.J.; Taubman, S.J.; Brown, P.; Iacono, M.J.; Clough, S.A. Radiative transfer for inhomogeneous
atmospheres: RRTM, a validated correlated-k model for the longwave. J. Geophys. Res. Atmos. 1997,
102, 16663–16682. [CrossRef]
63. Grell, G.A. Prognostic evaluation of assumptions used by cumulus parameterizations. Mon. Weather Rev.
1993, 121, 746–787. [CrossRef]
64. Niu, G.Y.; Yang, Z.L.; Mitchell, K.E.; Chen, F.; Ek, M.B.; Barlage, M.; Kumar, A.; Manning, K.; Niyogi, D.;
Rosero, E.; et al. The community Noah land surface model with multiparameterization options (Noah-MP):
1. Model description and evaluation with local-scale measurements. J. Geophys. Res. Atmos. 2011, 116, n/a.
[CrossRef]
65. Nakanishi, M.; Niino, H. An improved Mellor–Yamada level-3 model: Its numerical stability and application
to a regional prediction of advecting fog. Bound.-Layer Meteorol. 2006, 119, 397–407. [CrossRef]
66. Liu, C.; Ikeda, K.; Thompson, G.; Rasmussen, R.; Dudhia, J. High-Resolution Simulations of Wintertime
Precipitation in the Colorado Headwaters Region: Sensitivity to Physics Parameterizations. Mon. Weather
Rev. 2011, 139, 3533–3553. [CrossRef]
67. Lin, T.; Farley, R.D.; Orville, H.D. Bulk parameterization of the snow field in a cloud model. J. Clim. Appl.
Meteorol. 1983, 22, 1065–1092. [CrossRef]
68. Hong, S.Y.; Lim, J.O.J. The WRF single-moment 6-class microphysics scheme (WSM6). J. Korean Meteor. Soc.
2006, 42, 129–151.
69. Caneo, M. Sensibilidad a Diferentes Esquemas De Microfísica del WRF, en Chajnator-Chile. Ph.D. Thesis,
Tesis Para Optar al Título Profesional de Meteorólogo, Departamento de Meteorología, Universidad de
Valparaíso, Valparaíso, Chile, 2010.
70. Sikder, S.; Hossain, F. Assessment of the weather research and forecasting model generalized
parameterization schemes for advancement of precipitation forecasting in monsoon-driven river basins.
J. Adv. Model. Earth Syst. 2016, 8, 1210–1228. [CrossRef]
71. Willmott, C.J.; Robeson, S.M.; Matsuura, K. A refined index of model performance. Int. J. Climatol. 2012,
32, 2088–2094. [CrossRef]
72. Zegpi, M.; Fernández, B. Hydrological model for urban catchments—Analytical development using copulas
and numerical solution. Hydrol. Sci. J. 2010, 55, 1123–1136. [CrossRef]
73. Thompson, G.; Field, P.R.; Rasmussen, R.M.; Hall, W.D. Explicit forecasts of winter precipitation using an
improved bulk microphysics scheme. Part II: Implementation of a new snow parameterization. Mon. Weather
Rev. 2008, 136, 5095–5115. [CrossRef]
74. Thompson, G.; Eidhammer, T. A study of aerosol impacts on clouds and precipitation development in
a large winter cyclone. J. Atmos. Sci. 2014, 71, 3636–3658. [CrossRef]
75. Caseri, A.; Ramos, M.; Javelle, P.; Leblois, E. A space-time geostatistical approach for ensemble rainfall
nowcasting. In Proceedings of the 3rd European Conference on Flood Risk Management, Lyon, France,
17–21 October 2016. [CrossRef]
76. Kumari, M.; Singh, C.K.; Bakimchandra, O.; Basistha, A. DEM-based delineation for improving geostatistical
interpolation of rainfall in mountainous region of Central Himalayas, India. Theor. Appl. Climatol. 2017,
130, 51–58. [CrossRef]
77. Kumari, M.; Singh, C.K.; Basistha, A. Clustering Data and Incorporating Topographical Variables for
Improving Spatial Interpolation of Rainfall in Mountainous Region. Water Resour. Manag. 2017, 31, 425–442.
[CrossRef]

c 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access
article distributed under the terms and conditions of the Creative Commons Attribution
(CC BY) license (http://creativecommons.org/licenses/by/4.0/).