Integrating Geotechnical Investigation with Hydrological Modeling for Mitigation of Expansive Soil Hazards in Tabuk City , Saudi Arabia

Due to global climatic changes, flash floods are followed as a yearly disaster with high magnitude of influence. During the years 1981, 1988, 2010, 2012 and in January 2013, Tabuk city, northwest of Saudi Arabia suffered huge flash floods. These are major factors affecting on the swelling behavior of expansive Tabuk shale. The examined geotechnical properties of the surface and subsurface lithology of the sedimentary deposits distinguished Tabuk city into three zones. The expansive zone is spread in the middle and the non-expansive zones are distributed in the east and west of the city. The Watershed Modeling System (WMS) and Hydrologic Engineering Center (HEC-1) models were used to delineate and identify the drainage system and basin morphometry, where flash floods and accumulation of water might take place. Integration between geotechnical distribution maps of the expansion soil and surface hydrological data in terms of runoff maps was done. It has been identified the whereabouts the soils which have expansion characteristics and areas prone to flooding and surface runoff. They are helpful in defining the hazard zones map. Based on this map, it can be avoided constructions on the risk neighborhoods such as Al Qadsiyah, Al Maseif, Arrwdah, Al Nakhil and Al Rajhi. Also, it can suggest that the western side of Tabuk city is suitable for future urban extension. These results will help planners and citizens to create alternative development scenarios and determine their impact on the future urbanization patterns. Moreover, the direction of surface runoff flow or storm water discharge should be away from the expansion soil areas. Therefore, constructing dams on the outlet of the high-risk basins, south of Tabuk city is an important solution to control flash flood events, as well as increase groundwater recharge.

Tabuk City, Watershed Modeling System

Introduction
Soils and soft rock that tends to swell or shrink due to changes in moisture content are commonly known as expansive soils [1].Damages to infrastructure caused by expansi-ve soil movement have been widely reported in many countries such as Australia, China, United States, India and South Arica [2].In America, the economic loss of the expansive soil question approximately amounts to above $15 billion every year [3] [4] [5] more than twice the damage from natural disasters.
One of the major geotechnical problems in Kingdom of Saudi Arabia (KSA), as has recently been noticed, is the swelling potential of the soils upon change in moisture content.
Tabuk city, northwest KSA, is underlain by a shale; it is the primary source of expansive soils [6] [7] [8] [9] [10].The problem in Tabuk is a terrifying for the administrative, municipality officials and legal authorities; luxury villas, modern built, commercial buildings, schools, hospital, mosques and roads are subjected to upheaval forces.It is resulted in serious cracks, tilting, twisting and sticking doors and windows.Damages seen in some districts are estimated in the hundreds of millions US dollars [11] [12] [13].Conducted comprehensive studies covering greater parts of the Tabukcity are included in [6] [10] [12] [14] [15] [16].
The hazard posed by expansive soils is greatest in regions with pronounced wet and dry seasons [17].Therefore, building and infrastructure built in arid to semi-arid regions will be predisposed to expansive soil problems when compared with humid regions that maintain a rather uniform soil moisture condition throughout the year [18].In years of extreme temperatures or rainfall, the damage is the most severe [19].Factors that control the behavior of expansive soils have been thoroughly investigated by many researchers [1] [20] [21] [22].The primary factor is a change in water content.Pronounced changes in water content in soils seasonally occur in many areas of the world.The potential change in water content is generally attributed to the environmental conditions which influence volume change of soil is seasonal moisture variations due to rainfall, runoff and groundwater table.
Water is introduced to the sub-soil through infiltration of rainfall and runoff, rise in the groundwater level, and human activities (leaking from subsurface utilities i.e., water or sewage systems) could adversely affect ambient moisture conditions [23] [24].Surface drainage leads to moisture accumulation or ponding, which can provide a source of moisture for expansive soil by infiltration.
Tabuk city is hardship surface water scarcity compensated by flash flooding, which has been affecting the entire city for many years during the winter season.The catchment experienced major floods during the years 1981, 1988, 2010, 2012 in January 2013, causing unfortunately the death of many people and serious damage to residential buildings [25] [26].Hydrological processes must play a dominant role in the swelling of the soil and rocks.Surface hydrology deals with movement of water along the watershed because of precipitation.Drainage basin morphometry plays a significant role in delineating the areas prone to runoff or flooding [27].Poor surface drainage is a frequent problem associated with constructions on expansive soils.It leads to water accumulation or ponding, which can provide a source of moisture for expansive soil by infiltration.The extent of the infiltration is a function of the transverse and longitudinal gradients in the drainage tributaries, lag time of runoff, volume and peak discharge of runoff and the hydrologic soil groups.
The main objective of the current study is to investigate the runoff hazards on the expanded shale of Tabuk residential neighborhoods and to identify basin tributaries; surface runoff, flash floods location, and depressions allow accumulation or standing water, which activates the expansive shale in a localized manner leading to differential heaving of expansive soil.In addition, it will be carried an integration of geotechnical investigation and surface hydrology to construct a hazard expansive soil map of Tabuk city.It is helpful in determining the suitable future extension and development of urban areas.

Location of Study Area
Tabukcity lies in the northwestern part of KSA, between latitudes 28˚20'N and 28˚28'N and longitudes 36˚26'E and 36˚36'E.It is situated on an alluvial plateau with an elevation of about 760 m [28].Tabuk area is almost flat; it gently slopes downward to the north and northeast with gradually moderate relief.The plateau is frequently cut with several wadis, which drain northwards for some 20 km to the Qa'a Sharawa dry lakebed (Figure 1).This area has the lowest relief and forms an internal drainage for the wadis.

Geotechnical Properties
The soil samples (disturbed and undisturbed) were collected from the drilled 72 boreholes to depth varied from 10m to 25 m distributed in the Tabuk city (Figure 1).Laboratory tests are carried out to attest the potential expansiveness of the soils by grain size distribution, soil moisture content, Atterberg limits, and swell tests.All of these tests were performed in conformity with the American Standards, ASTM [29].

Hydrological Methodology
The used rainfall data is collected from the Hydrology Division of the Ministry of Water and Electricity [30].The available rainfall records cover a period of 49 years (1965 to 2004).HyFrAn PLUS (Hydrological Frequency Analysis) program [31] had been used for the available rainfall data to obtain the corresponding rainfall depth values of the probabilities distribution (Gumbel, General Extreme value GEV, Pearson Type III and Log-Pearson Type III) for different return periods.The Chi-square test comparing computed values with observed values was carried out to identify the best-fit method.
Delineating drainage networks is essential to predict flow of surface water runoff, and consequently the flash flood hazards zones [32] [33].Basin delineation was condu-cted using the Shuttle Radar Topography Mission (SRTM v.3, 1 arc) 30-m digital elevation model "DEM" (Projected Coordinate System: WGS_1984_UTM_Zone_37 N) in the Watershed Modeling System (WMS, v.10).WMS and Global Mapper (v.17) were used to provide most of the morphometric parameters of the studied basins and their stream attributes [34].
The time of concentration constitutes an important factor in flood assessment studies, since it is the time required by runoff to travel from the most distant point to the basin's outlet point and lag time as the time from the center of mass of excess rainfall to the hydrograph peak [35].For ungagged watersheds, the Soil Conservation Service (SCS) suggests that the UH lag time (t L ) may be estimated as 0.6 times the t c (t L = 0.6 t c ).
A number of detailed hydrologic models have been developed to estimate the peak discharge and runoff hydrograph for a given rainfall distribution, such as the Hydrologic Engineering Center model (HEC-1), the SCS's Technical Release-20 model (TR-20), the SCS's Technical Release-55 model (TR-55).Sonbol [36] used HEC-1 with SCS Unit Hydrograph option to simulate peak discharges in Sinai, Egypt and concluded that the model reasonably predicted observed peak flow.The SCS method is suited for humid, semi-arid, and arid conditions [37].The SCS method has been successfully applied to several ephemeral watersheds in arid environments similar to Saudi Arabia [38] [39] [40] [41] [42].HEC-1 demonstrated acceptable performance in rainfall runoff modeling for different case studies in arid regions [43] [44] [45].
In the present study, the selected computer Rainfall-Runoff model was the Hydrologic Engineering Center (HEC-1) [32] under the Watershed Modeling System (WMS v. 10) software as it was successfully applied for hydrologic modeling of watercourses and flash floods of the daring system that flow into Tabuk city.
The peak flood discharge and volume of runoff of eight drainage basins were estimated by applying the Soil Conservation Services (SCS) method using the Hydrologic Engineering Center (HEC-1) model under WMS software program.The HEC-1 was run for 7 return periods (2, 3, 5, 10, 20, 50 and 100 years) using the soil conservation service curve number (SCS CN) method to estimate the infiltration losses with an average antecedent moisture condition (AMC II) [37].The SCS CN is a function of the hydrologic soil group, land use, land cover and antecedent moisture conditions [46].
Based on the geological, soil, land use and land cover characteristics, the weighed curve number (CN) is 80 for the study area catchments.

Geological and Geotechnical Setting
Geology of the study area was based on the geology of Tabuk Quadrangle mapped (GM-137 Sheet No. 28B) by Janjou [47] (Figure 2), field investigation and description and analyses of the drilling boreholes samples; in addition to the general geology and subsurface hydrology [48]- [56].

Surface Geology
Tabuk area is located within the sedimentary cover in the northwestern part of the Arabian Peninsula Shelf.The outcrop rock units in the Tabuk city belong to Paleozoic and Quaternary sediments (Figure 2).
Qasim and Sarah Formations represent the exposed Paleozoic sediments in Tabuk city, where it outcrops in several hills within the town limits (Figure 3).Qasim and Sarah Formations are belonging to Tayma and Tabuk groups, respectively [47].It composed of a thick sequence of shale, siltstone, and sandstone units that belong to marine and continental facies.
It rests conformably but with a sharp contact on the sandstones of the Saq Formation.Age of Qasim Formation is Early to Late Ordovician.Its maximum thickness is about 1070 m but diminishes gradually in the southeast direction.The Formation is subdivided into four members, from bottom to top: Hanadir, Kahfah, Ra'an and Quwarah.In Tabuk city, Ra'an and Quwarah members are present.Ra'an Member consists of gray and dark-red, finely laminated silty clay to greenish gypsiferous clay (Figure 2 and Figure 3(a)).Quwarah Member is overlies Ra'an member.It is late Ordovician age Figure 2. General geological map of Tabukcity (Modified from Janjou [47]).[57].Quwarah Member consists of medium to fine-grained micaceous sandstone and silty clay (Figure 2).Quaternary sedimentary deposits overly Paleozoic sedimentary rocks, covered most of the study area (Figure 2).It represented by alluvial deposits (active and inactive terraces), khabra and sabkhah deposits (Figure 3(c) & Figure 3(d)).They consist mainly of silt, clay and blown sand with variable proportions of gravel and rock fragments of basement complex and vesicular basalt (scoria), its thickness ranges from 0.5 m to 6 m.

Subsurface Lithology
Through the nine lithological cross sections (Figure 1 & Figure 4), it can be divide Tabuk city according to the subsurface lithology into three zones; these are western, central, and eastern zones.
The western zone consists of yellowish to brownish color, very fine to medium grain size, moderately sorted sandstone rock.It is found from surface to the extending depth of the borings (Figure 4).It is interbedded with thin lamination of shale.In some loca-

Geotechnical Properties of Expansive Soil
Shale is the primary source of expansive soils in Tabuk due to the presence of clay minerals (smectite and illite) that are derived from shale [6] [7] [8] [9] [58].Sabtan [10] found that, the degree of soil expansion in Tabuk ranges from low to very high.
From the geotechnical properties and mineralogical composition of Tabuk shale, it had low to mediums welling potential, according to the criteria proposed by [1] [59] [60] [61] [62].The mineralogical composition of clay constituents of the shale, investigated using X-ray diffraction analysis reveals that the predominant clay minerals are kaolinite with minor amounts of illite [62].The expansive lattice-type clay mineral (i.e., montmorillonite) cannot be traced in the shale.Both kaolinite and illite minerals are known for being low to moderately expansive materials.
The geotechnical properties of the surface and subsurface lithology of the sedimentary deposits and/orrocks showed that, shale layers of (Ra'an member) have expansive characteristic, whatever sandstone (Sarah Formation) and clayey sand with gravel and rock fragments (Quaternary alluvial deposits) layers are non-expansive.

Distribution of Expansive Soil
Based on the field geology, subsurface lithology, and geotechnical properties of the investigated soil samples of Tabuk city, it can be distinguished Tabuk city into three zones (Figure 5) from geotechnical view as follow.
The expansive zone is spread in the middle of the city, composed of shale.This zone is narrow in the south and become wider in the north of the city.The neighborhoods of Arrwadah, Al Farabi, Annahdhah, Al Nakhil, Rahil, Al Rajhi, Al Qadsiyah and Al Masef are situated in the expansive zone and suffered significant damage to constructions.
The non-expansive zones are distributed in the east and west of Tabuk city, represented by alluvial deposits (silt, clay and blown sand with variable proportions of gravel and rockfragments) and sandstone eastwards and westwards zones.The neighborhoods of the east, west, and southern parts of Tabuk city are underlain by non-expansive deposits and/or rocks and did not show problems to constructions.

Rainfall Analysis
Climate conditions play an important role in defining the hydraulic response of the watersheds existing in Tabuk city and it's environ.The most important factor affecting the hydraulic behavior of the wadi basins is rainfall.Its duration, intensity, distribution, and return periods are major influences.
In Tabuk city, the amount of rainfall is irregular through the years and through the months (Table 1, Figure 6).Monthly rainfall distributions are uneven, with little or no rainfall in the May, June, July, August and September.The maximum rainfall occurs in January, November, October and December (Figure 6).The daily maximum-recorded rainfall in the region was 96.2 mm on the 19 th of December 1985 [30].
The computed values and visual inspection revealed that Pearson Type III pdf produced the best fit in most of cases (Figure 7).The predicated maximal daily rainfall depth for 2, 3, 5, 10, 20, 50 and 100 return periods of Pearson Type 3 distribution is  presented in Table 2.

Delineation of Drainage System
Analysis DEM (Figure 8) of the area shows that the high lands occur towards the south, east, and the west of Tabuk area (Figure 9); the northern part occupies the relatively lower elevations towards Qa'a Sharawra.The general slope is towards the north (0.13 m/m).Flatness terrain present in the last eight kilometers of Tabuk city, starting from Arrawdah district to the northern border of Tabuk city, the slope is nearly 0.0002 m/m (Figure 9).

Peak Discharge
The results of the peak discharge and volume of runoff estimated for return periods of 2, 3, 5, 10, 20, 50 and 100 years from the run of the HEC-1 model of the drainage basins at the outlet of Tabuk area are represented in Table 4. From the hydrograph (Figure 11), it is noticed that the discharge rate of the drainage basins is highly variable from  4 and Figure 11).Qp: Peak discharge (m 3 /sec); b V: Volume of Runoff (m 3 ).

Discussion
The hazards posed by expansive soils have significant relation with water runoff.It leads to water accumulation or ponding, which can provide a source of moisture for expansive soil by infiltration.In Tabuk city, some parts of the watershed have runoff or flash floods on areas comprise expansive soil and hence the problem is duplicated due to first the flash flood and second the heaving of the soil.On the other hand, some locations have no flow, consequently no problems from this point of view even if the expansion soil is found.Volume of surface runoff can give singe of the potential swelling soil locations, where the expansion is linearly related to surface runoff.The drainage pattern represented in.Delineating the drainage pattern effect on Tabuk districts is helpful in determining the directions of running water and the possible locations for accumulation (Figure 12(a)).The flow direction of surface runoff or storm water discharge should be away from the expansion soil zones.Therefore, the integration was carried out between geotechnical distribution of the expanded soil and the drainage network over Tabuk neighborhoods (Figure 12(b)).Consequently, the riskiest zones in the Tabuk city have been identified (Figure 13).It found that, the main watercourse of Wadi Dab'an and Abu Nushayfah passes through area contains expansive shale.The residential neighborhoods built on the expansive shale in the middle part of Tabuk city affected by surface runoff flow of the tributaries of basins of Wadi Abu Nushayfah, Na'am and Atanah outlet, represent high-risk zones.These high-risk residential  neighborhoods are Al Qadsiyah, Al Maseif, Arrwdah, Al Nakhil and Al Rajhi (Figure 13).Flatness in the last eight kilometers of the basin (B-3) of Wadi Dab'an is starting from Arrawdah district to the north.It reduces surface runoff velocity and hence increases accumulation and stagnant water for a long time over the expansive soil, consequently maximizing the heaving of the expansive shale.The low-risk zones includes the western and the eastern neighborhoods (Figure 13), it is built on non-expansion soil or/rock (sandstone rocks in the west and silt, sand with boulders cobble rocks in the east).

Conclusions
Tabuk city contains expanded shale, in the central part of the city.It expands when it absorbs more water and creates the potential to exert enough force on a building or other structures to cause damage.In addition, Tabuk city exposed too many flash floods, mainly from large drainage basin of Wadi Na'am, Al Baqqar and Dab'an which covered an area of about 2860 km 2 and had high peak discharge and volume of runoff 2201 m 3 /s and 151,574,352 m 3 , respectively and 41.19 m 3 /s and 2,930,891 m 3 for returned period 100 and 2 years, respectively.In the event of flash floods or surface runoff, water flows and collects in some places that contain the expanded shale, causing its expansion and thus the destruction of the buildings is built upon it.The direction of flow of surface runoff or storm water discharge should be away from the expansion soil zones.The quicker the water can flow out, the better it is.The integration was carried out between geotechnical distribution maps of the expansion soil and flow of runoff, and the riskiest places had been identified in Tabuk city.The resident neighborhoods built on the expansive shale in the middle of Tabuk city and affected by surface runoff of the tributaries of the basin of Wadi Abu Nishayfah, Na'am and Atanah outlet, represent high-risk zones.These high-risk neighborhoods are Al Qadsiyah, Al Maseif, Arrwdah, Al Nakhil and Al Rajhi.The low-risk zones include the western and the eastern districts; it is non-expansion soil.
This study showed the importance of characterizing geotechnical properties and hydrological elements as part of a comprehensive environmental analysis for future urban planning.Based on the hazard map of the expansive soils of Tabuk city, it found that the western side of the city lying in the upstream of the small basin and the subsurface layers is composed of sandstone.It can suggest that, the western side of Tabuk city is suitable for future urban extension.

Figure 1 .
Figure 1.Location map showing distribution of drilled boreholes and cross sections, Tabuk City, KSA.

Sarah
Formation is belonging to Tabuk Group.It is overlies Qasim Formation of Tayma Group.The age of Sarah Formation is of late Ordovician to early Silurian.Sarah Formation consists of medium to coarse-grained sandstone (Figure 2 and Figure 3(b)).
tion, sandstone overlies by thin sandy clay, clayey sand or sand blown layers, ranges in thickness from 0.5 m to 1 m.Sandstone represents the subsurface layers of the Al Aaen, Tabuk University, Al Quds Districts.

Figure 5 .
Figure 5. Distribution of the expansive and non-expansive soils in Tabuk city, KSA.

Figure 6 .
Figure 6.Statistical distribution of daily rainfall depth records over months of the year at Tabuk Weather Station 403,750.

Figure 7 .
Figure 7. Gumbel distribution fit for annual maximum daily rainfall at Tabuk weather station 403,750.

1 Figure 9 .
Figure 9. Cross section of Tabuk area from West to East (A), longitudinal section from South to North (B) Tabuk, KSA.
one basin to another, according to the basin area.Three basins (B-3, B-1 and B-5) have a high peak discharge and volume of runoff.The maximum peak discharge for the lowest return periods of 2 year of the Basins 3, 1 and 5 are 41.19, 20.39 and 15.5 m 3 /s, respectively, whatever for the highest return period 100 year are 2201, 1085 and 824 m 3 /s, respectively (Table

Figure 12 .
Figure 12.Drainage system overly on Tabuk districts (A) and integration between geotechnical distribution of the expanded soil and the drainage network (B) Tabuk city, KSA.
Two major strategies to mitigate hazards resulted from flash floods and runoff on the swelling soil zone (high-risk zone) in Tabuk districts.It is essentially based on prevent or reduce flow of surface runoff and flash floods from the offside Tabuk city and drainage storm water from inside Tabuk districts.The hazards can be mitigated by the following procedures: 1) Construction of water flow retardation structures such as dams, artificial ponds, and orifice divert the course of the wadi etc. in upstream along the main courses of Wadi Abu Nushayfah and Ghudayy, Na'am, Al Baqqar, Dab'an and Atanah and Ghuwaul.This would serve two goals; the first decreasing the flash floods flow to Tabuk city and thus reduces the chance of swelling soils; and the second increasing the contribution to the groundwater aquifer, which depends essentially on rainfall recharge.2) Nature wadis and artificial open channels that pass into expansion soil areas (hazard zones) within Tabuk city should be refine and lining, to prevent infiltration into the subsurface expansion soil.

3 )
Benefit from the surface water discharge by matching the paved street layout within the affected districts with the direction of flow, to reduce time of concentration and increase the speed of the discharge.In addition, the direction and slope of the main streets should be in the north direction toward Al Qa'a depression.4) Construction of storm water network in the hazard zones of Tabuk neighborhoods.

Table 1 .
Statistical analysis of daily rainfall data for the period from 1965 to 2004 at Tabuk Weather Station 403,750 Latitude: 28.38, Longitude: 36.63,Altitude: 778.

Table 2
. Estimated return periods for maximum daily rainfall depth (Pearson Type III) for Tabuk weather station 403,750.

Table 3 .
The main basins geometric parameters of Tabuk area, KSA.