Heavy Rainfall Simulation over Sinai Peninsula Using the Weather Research and Forecasting Model

Heavy rainfall is one ofmajor severeweather over Sinai Peninsula and causesmany flash floods over the region.The good forecasting of rainfall is very much necessary for providing early warning before the flash flood events to avoid or minimize disasters. In the present study using the Weather Research and Forecasting (WRF) Model, heavy rainfall events that occurred over Sinai Peninsula and caused flash flood have been investigated. The flash flood that occurred on January 18, 2010, over different parts of Sinai Peninsula has been predicted and analyzed using the AdvancedWeather Research and Forecast (WRF-ARW)Model.The predicted rainfall in four dimensions (space and time) has been calibrated with the measurements recorded at rain gauge stations.The results show that the WRF model was able to capture the heavy rainfall events over different regions of Sinai. It is also observed that WRF model was able to predict rainfall in a significant consistency with real measurements. In this study, several synoptic characteristics of the depressions that developed during the course of study have been investigated. Also, several dynamic characteristics during the evolution of the depressions were studied: relative vorticity, thermal advection, and geopotential height.


Introduction
Heavy rainfall is one of the major severe weather in Egypt particularly in arid and semiarid regions especially if it is steep and mountainous regions such as Sinai Peninsula and the Eastern desert of Egypt.Short duration of heavy rainfall over a relatively small drainage area can lead to devastating flash flood, consequently causing a number of fatalities and tremendous damages.
It can destabilize soils along mountain slopes, resulting in landslides and mudslides that cause severe damage to nearby villages.
Heavy rainfall is usually resulting from individual mesoscale storms or mesoscale convective systems embedded in synoptic-scale disturbances [1].High-resolution observations and numerical modeling technique are required to better predict heavy rainfall events, where the forecasting of heavy rainfall is very important for many decision makers who are sensitive to the occurrence of precipitation.An accurate quantitative precipitation forecast can identify the potential for heavy precipitation and possible associated flash flooding, as well as provide information for hydrological interests.
Since heavy rainfall and flash flood can lead to severe damage and losses for both life and infrastructure, the need to warn people in advance is, thus, an important goal, but heavy rainfall and flash floods forecasting is considered a difficult task, particularly in arid and semiarid regions, since they take place in very short time interval.Hence, rainfall forecasting is very much necessary for providing an early warning before the flash flood events to avoid or minimize disasters.An early warning system (EWS) for flash floods has been developed for part of the Sinai Peninsula of Egypt, a hyperarid area confronted with limited availability of field data, limited understanding of the response of the valley (wadi) to rainfall, and a lack of correspondence between rainfall data and observed flash flood events [2].
Sinai Peninsula contains many complex terrains; so rainfall forecasting for complex terrains is challenging and requires numerical simulations at very high resolution.
In arid areas such as Sinai in Egypt, rainfall is mainly caused by squall line and convective cloud mechanisms and by low-intensity frontal rain, causing storm floods [3].The complexity of terrain influences the weather in a variety of ways.Under stable atmospheric conditions, the terrain generates internal gravity waves that distribute momentum over wider areas.These processes may be related to strong winds and turbulence that influence the air traffic.Under unstable conditions, convective clouds and precipitation are generated over complex terrain which can grow into severe thunderstorms.All the mentioned processes occur on spatial length scales smaller than 100 km, usually even smaller than 10 km.
Because of the above-mentioned challenges, the mesoscale meteorological model named the Advanced Research Weather Research and Forecasting (WRF-ARW) is selected for rainfall forecasting in this research.WRF is developed mainly by the US National Centre for Atmospheric Research (NCAR) in collaboration with many other research centres and universities.WRF allows forecasting weather in complex terrains such as the one in Sinai and in the same time is considering orographic features.It is suitable for a broad spectrum of applications across scales ranging from meters to thousands of kilometres.
The main objectives of this study can be summarized as follows: (a) the understanding of the processes leading to the flood event at Sinai using WRF?(b) the investigation of the capability of WRF to simulate successfully the rainfall amount during a flood event?

Study Area
The Sinai Peninsula (Figure 1), study area, is located in the far northeast of Egypt.It represents about 6% of Egypt's area which is about 61,000 km 2 .The Sinai Peninsula comprises a wedge-shaped block of territory with its base along the Mediterranean Sea coast to the north and its apex bounded by the Gulfs of Suez to the west and Aqaba to the east.
Its southern portion consists of rugged, sharply serrated mountains.These reach elevations of more than 2,400 meters; among them is Mount Catherine, Egypt's highest mountain, which has an elevation of 2,642 meters.The central area of Sinai consists of two plateaus, Al-Tih and Al-Ajmah, both deeply indented and dipping northward toward Wadi (Valley) El Arish.
Toward the Mediterranean Sea, the northward plateau slope is broken by dome-shaped hills; between them and the coast are long, parallel lines of dunes, some of which are more than 100 meters high.
It is dissected by the largest wadi in Sinai, Wadi El Arish, which emerges from elevated gravelly plains and terraces in the south to a distance of about 20 km till the coast Mediterranean Sea to the north.The rainfall and floods are the only sources of renewable water resources in Sinai Peninsula.

Climate of Sinai Peninsula. It is characterized by the
Mediterranean climate in the northernmost of Sinai to be close to the desert and semidesert climate to the south.
Much of the Sinai is hot, or very hot with a higher temperature inland, but there is a more temperate region near the north coast and over the mountains.During the period May/June to September/October, the mean daily maximum is 28 ∘ C to 37 ∘ C in the North, 31 ∘ C to 42 ∘ C near the south coast, and 35 ∘ C to 41 ∘ C inland.Minimum temperatures average between 20 ∘ C and 25 ∘ C in the summer.
The winter season is a little less harsh with day high's in the mid-teens and possible 20's, and evenings often falling to around 6 ∘ C-10 ∘ C, and may drop below 0 ∘ C.
The amount of rainfall in Sinai decreases from the northeast towards the southwest.The greatest amount of the annual rainfall was found at Rafah station (304 mm) in the northeast.The annual rainfall average is about 120 mm along the Mediterranean coast.It decreases in the uplands to the south to about 32 mm.The annual average of rainfall all over Sinai is about 40 millimeters, 27 millimeters from it is estimated to come from one storm that may provide 10 millimeters at a time.
Along the Mediterranean Coast, 60% of the rain occurs in the winter, while 40% falls during the transitional seasons.
The annual rainfall in the southern region is significantly less than the northern one, where it reaches 20 mm in the coastal areas over the Gulfs of Aqaba and Suez.Its amount increases to 70 millimeters over the mountain regions.
The rainfall in autumn and winter seasons is ranging from medium to heavy especially at some high terrain areas.The rainfall is completely nonsexist in summer season.

Data and Methodology
WRF model is used in this research with its nesting capability to simulate the heavy rainfall with fine grid spacing.
During the storm mesoscale, convective systems are highly interacted with synoptic-scale environment.
During this study we conducted multinested experiment for four domains with different horizontal resolutions of 81 km (DM1), 27 km (DM2), 9 km (DM3), and 3 km (DM4), respectively, as shown in Figure 2. The second domain (DM2) will be used to simulate the synoptic situation over Egypt on 18 January, 2010, while the last domain (DM4) will be used to investigate the ability of the WRF model to simulate and predict rainfall to compare it with rain gauges at Sinai.
The model output is adopted to produce its output every hour to be compatible with the rain gauges observations and two-way nesting domain starting from 17 to 19 January, 2010.The 27 vertical layers with the model top of 50 hPa are used in this study.
The initial and lateral boundary meteorological data which used to run the model has been downloaded from the National Centers for Environmental Prediction (NCEP), global final analyses on 1 ∘ ×1 ∘ degree from Global Forecasting System (GFS), and it is updated every six hours.
The hourly rainfall data, measured in millimeter (mm), was obtained from the rain gauges stations installed by the Water Resources Research Institute (WRRI); WRRI is the Egyptian governmental research institute with the mandate for flash flood management.
WRF offers multiple physics options that can be combined in any way.The options typically range from simple and efficient to sophisticated and more computationally costly and from newly developed schemes to well-tried schemes such as those in current operational models.Table 1 represents the selected physics of WRF model.

Results and Discussions
Egypt is subjected to thunderstorms and heavy rains started at the north coast, the Red Sea, and the Sinai Peninsula during 18 January, 2010.This extreme weather and intensive rainfall led to flash flood events over Sinai Peninsula.
The development of intensive weather events that invade Egypt during January 18, 2010, were characterized by "exceptional and extremely heavy rainfall, " which affected a wide part of Egypt, including Sinai Peninsula and fatal to some Bedouin tribes located in its path.This heavy rainfall is accompanied by strong winds and thunderstorms.
The Mediterranean region is considered to be one of the most cyclogenetic areas in the world, usually favoring the development of weak low-pressure systems.Occasionally, these systems develop into deep cyclogenesis that cause a series of severe weather events as they cross the Mediterranean.
Generally, the winter season in Middle East is known to be associated with the development of low-pressure systems (Cyprus lows) over the eastern Mediterranean Sea [4].Such systems have been observed at upper levels, particularly at 500 hPa.Rainfall has been observed to be heavy over Sinai and extended to southern parts of Egypt during our case study in winter season, January 18, 2010.Kahana et al. [5] found that Mediterranean Sea cyclones deepen in general, in association with upper level troughs.Upper level troughs are regarded as the key factors in activity of midlatitudes, such as the Mediterranean Sea system, particularly in their front sides, where positive vorticity advection and enhanced convection take place [6].Ferraris et al. [7], however, found that the polar continental air from central Asia moving toward the eastern Mediterranean, and interacting with the relatively warm sea surface temperature (SST), produces enhanced lower level instability.
Because of the importance of the Mediterranean cyclogenesis, many studies on weather systems in the region have been conducted.The synoptic situation started by the development of the subtropical jet stream in the upper troposphere.The maximum wind speed ranged between 120 and 135 knots as shown in Figure 3, with the southern extension of the polar jet to the south of the Mediterranean.
Figure 3 shows the wind speed of the subtropical branch of the jet stream at 200 hPa for 4 consecutive times per day for our case study, with the southern extension of the polar jet to the east of the Mediterranean Sea.The synoptic situation can also be analyzed by referring to the midtroposphere.Figure 4 shows the geopotential height at 500 hPa level for the same day; the synoptic system at 500 hPa level also follows the subtropical jet development.

International Journal of Atmospheric Sciences
The axis of the trough which is located over east Europe shifts southeast.This trough line orientation indicates that the affected area has extended and engulfed the whole eastern Mediterranean area.This development agrees with Krichak et al. [8,9] and Krichak and Alpert [10] who showed that the meridional orientation of the upper trough line represents the major synoptic process of the development of the East Mediterranean cyclone.As the subtropical jet stream shifts south and splits, the system with high geopotential heights moves eastward and allows the development of a low system in the eastern Mediterranean.
The relative vorticity advection is considered as one of the most important parameter that leads to occurrence of heavy rainfall and flash flood events in Sinai Peninsula.As shown in Figure 5, the positive relative vorticity advections leads to move the upper atmospheric trough from west to east (steering of trough) and reached to Sinai Peninsula.
From Figure 6, we can also follow the temperature advection that calculated using centered finite difference scheme of second order.In the stage of rapid development, there is a cooperative interaction between the upper level and surface flows; strong cold advection is seen to occur west of the trough at the surface, with warm advection to the east.This pattern of thermal advection is a direct consequence of the fact that the trough at 500 hPa lags (lies to the west of) the surface trough.These are consistent with the results of Charney [11], Holton [6,12], and Bluestein [13], where the growth of midlatitude cyclone is associated with a meridional gradient of temperature.We can expect that the cold advection on the west of the trough is responsible for the decreasing of temperature over Egypt in that period.
The mean sea level pressure (MSLP) of the case study is shown in Figure 7 and follows the system development at the 850 hPa level as shown by the heat advection (Figure 6).The intensification of Azores high, to the west of the region, leads the low system to be developed in the eastern of the Mediterranean region.In addition, this produces a larger air mass temperature difference between the east and west of the low trough line and allows the Red Sea trough to extend northward.The Red Sea trough extension and the meridional orientation of the upper trough permit favorable conditions for the formation and development of the eastern Mediterranean cyclone.
Figure 8 shows the total precipitation associated with the east Mediterranean storm.The rainfall amount increases as the low-pressure system deepens and the precipitation area extends eastward following the surface low.

International Journal of Atmospheric Sciences
Evaluation of the performance of WRF model for prediction of heavy rainfall events over Sinai Peninsula will be done by comparison of the rainfall predicted by the WRF model and the rain gauges measurements; the complete information of the rain gauges is shown in Table 2.
One of the very important statistical tools is the root mean square error (RMSE) which gives a good overall measure of model performance.The weighting of (predictionobservation) by its square tends to inflate RMSE, particularly when extreme values are present.With respect to a perfect model, the root mean square error should approach zero.On the other hand, the mean bias (MB) is considered also to evaluate the model performance which represents the The equations for RMSE, MB and its percentages are given as follows: where  is number of observations,   and   are the predicted and observed values, respectively, and   is the average of observed values.These statistical measures are used extensively for evaluation of model forecasts of precipitation (e.g., [14][15][16]).The International Journal of Atmospheric Sciences Mean sea level pressure (hPa) at 00Z18JAN2010 quantitative validation of the simulated precipitation is done by calculating their differences with the corresponding rain gauges.
Figure 9 shows the comparisons between the predicted hourly rainfall and the measured ones in different regions over Sinai Peninsula during January 18, 2010.
From Figure 9, it can be noticed that the predicted rainfall by WRF model is in a good consistency and harmony with observed rainfall for all rain gauge stations.
In El-Gudairate, El-Rawafaa, El-Haithy2, and El-Themed Stations rainfall prediction by WRF is very close to the measured ones.But in the other stations such as El-Haithy1 and Ras-Shira, it gave small differences with the measured rainfall during January 18, 2010.
Table 3 shows the percentages of root mean square error and mean bias for predicted rainfall and the measured ones at different stations.
It can be deduced from this table that the maximum values of root mean square error and mean bias are for El-Haithy1 station, while the minimum values of root mean square error and mean bias are for El-Gudairate, El-Rawafaa, El-Haithy2, El-Themed, and Ras-Shira Stations.
This means the WRF model has high-performance rainfall prediction for all stations except El-Haithy1 station.The WRF model was able to simulate the synoptic situation which invaded Egypt and accompanied by heavy rains started at the north coast, the Red Sea, and the Sinai Peninsula during 18 January, 2010.
In the same time, the model simulation of the synoptic situation was in a good compatibility with the previous studies.
It may conclude that the performance of WRF model approximately gave a good and reasonable rainfall prediction compared to measurements at all stations over Sinai Peninsula.

Conclusions
The WRF model was able to simulate the synoptic situation which invaded Egypt and accompanied by heavy rains started at the north coast, the Red Sea, and the Sinai Peninsula during January 18, 2010.
In the same time, the model simulation of the synoptic situation was in a good compatibility and agreement with the previous studies over the Mediterranean region.
It may be concluded that the performance of WRF model approximately gave a good and reasonable rainfall prediction compared to measurements at all stations over Sinai Peninsula.

Figure 1 :
Figure 1: Study area and rain gauges location.

Figure 3 :
Figure 3: 200 hPa subtropical Jet stream (knots) at over northeastern Africa with the southern extension of the polar jet to the east of the Mediterranean at January 18, 2010.

Figure 8 :
Figure 8: Accumulated rainfall (mm) during January 18, 2010 plotted every six hours with the maximum volume of rainfall over Sinai Peninsula.

Table 1 :
Selected physics of the WRF model options.

Table 3 :
The percentages of root mean square error and mean bias for rainfall from WRF and the measurements at different stations.