Received 30 July 2016; accepted 22 August 2016; published 25 August 2016
The stability of rock slopes is significantly influenced by the structural setting of the rock in which the slope is excavated which refers to naturally occurring breaks in the rock such as bedding planes, joints and faults which are naturally called discontinuities. The properties of discontinuities relative to stability include orientation, persistence, roughness and infilling  .
The main aims of the present study are identifying the regional stresses led to these discontinuities, assessment of the effect of rock discontinuities (orientation, density, spacing, persistence, and aperture) on rock slope stability and the influence of rock slope stability on the proposed Abo-Hederya quarry that affects the safety of the quarry.
2. The Study Area
The study area is located in Abo-Hederya (longitude 47˚19'42.9''E, and latitude 23˚26'28.2''N) proposed quarry in Missan Governorate 95 km east of Amara city within Mukdadiya Formation (Figure 1).
The exposed formation in the study area is Mukdadiya formation. The name Mukdadiya formation was used by  to replace the Lower Bakhtiari formation  . Mukdadiya (Lower Bakhtiari) Formation comprises up to 2000 meters of fining upwards cycles of gravely sandstone, sandstone and red mudstone. The Formation was deposited in fluvial environment in a rapidly subsiding foredeep basin. The age based on the vertebrate fauna is Late Miocene  . The lower contact is recognized lithologically by the first occurrence of the pebbly sands. It is distributed mostly in the Foothill Zone. The upper contact of the formation that overlain Mukdadiya Formation is diachronous  ; it is recognized lithologically by the appearance of the first thick (more than 1 meter) conglomerate horizon  . In the study area the upper contact of Al-Mukdadiya formation is an erosional surface (angular unconformity surface), with the overlain quaternary deposits.
The study area is located within Himreen subzone, foothill zone, which is a part of unstable shelf according to tectonic map of Iraq (1996) (Figure 2)   . It is in an extension of asymmetrical Himreen anticlines, with
Figure 1. Location and geological map of study area  .
Figure 2. Tectonic map of Iraq (modified after  ).
fold axis trending NW-SE, and bounded between two thrust faults with NW-SE direction (Figure 1), as a result of the collision between the Arabian and Iranian plates. The dip of the beds has an average value ranging from 20 to 8˚ because of their location near the hinge area of the fold as in Figure 1.
Field work was carried out for 8 days, from 7th to 14th Feb., 2015. It covered the discontinuity, lithology, and geomorphologic features survey. It was carried out in 10 stations. The area of each station was not less than 2.5 × 2.5 m, while the location and distribution of these stations were mainly based on the availability of exposures and the importance of engineering and structural position. At each station notes were collected for beds (lithology, thickness and attitude) and discontinuities (orientation, density, spacing, persistence, and aperture).The office work of the present study included calculating, representing and analyzing of the field data.
4. Geomorphologic Features of Study Area
The most common geomorphologic features found in the studied area were as follows:
Strike Ridges: Strike Ridges are Common within Mukdadiya as in Figure 3(a).
Rill Erosion: This feature is more effective and influential on the soft rocks that cover large areas (Figure 3(b)).
Gully Erosion: This feature is effective on solid rocks like Sandstone (Figure 3(c)).
Badland Erosion: This type of sculpture is common in Mukdadiya Formation composed of complex valleys and irregular lands (Figure 3(d)).
Mesa: Mesa formed when the beds are horizontal or near to horizontal (or sub horizontal) as in Figure 3(e).
Sand Sheets: They are scattered among different part of the studied area as in Figure 3(f).
Discontinuity is a general term for any mechanical break that has zero or relatively low tensile strength, it is a collective term that includes feature such as joint (a barren closed discontinuity on which no measurable slip is observed), vein (a discontinuity in which visible mineral infill took place), fissure (a discontinuity with a void space), fault (a discontinuity along which a detectable slip has occurred), and bedding plane (a discontinuity that separates the adjacent rock beds)  -  . The discontinuities strongly effect the mechanical and hydrological properties of a rock mass, in term of its strength, deformability, stability, porosity, and permeability  , accordingly, the discontinuities play a major role in the design and safety of civil works  and the occurrence and movement of ground water and hydrocarbons  (H. There are many field and laboratory parameters which can be used to evaluate any site under study. The type of these parameters depends on the purpose of the study  . In the present study only five field parameters were used. These are: orientation, density, spacing, persistence, and aperture. The parameters have direct relations to the geotechnical aspects of the study area.
5.1. Discontinuity Orientation
Discontinuity orientation is the attitude of discontinuity in space. It could show grouping into sets (group of parallel discontinuities that maintain their orientation over considerable area) and system (applied to two intersecting discontinuity sets with relatively fixed angle of intersection)  .
The discontinuity orientation can be described by the dip of the line of steepest inclination measured   ,
Figure 3. Geomorphologic Features of study area ((a) Strike Ridges; (b) Rill Erosion; (c) Gully Erosion; (d) Badland Erosion; (e) Mesa; (f) Sand sheet).
from horizontal, and by the dip direction measured clock wise from true north. Example: dip direction/dip angle (025˚/45˚). The average readings of the discontinuity orientations in the stations of the study area were designated in Table 1 as 1; 2; 3; 4; 5; 6; 7; 8; 9; 10. Discontinuity (fracture) analysis according to the principal tectonic axes was done to detect the geometrical and genetic relationship with the stresses in the region. Fractures were representative in the form of shear fracture sets (okl, hko, hol and hkl) and extension fracture set (ac).
The common discontinuity sets are persistently developed and are characterized by being normal to the bedding, while the less common discontinuity sets are characterized by being inclined at various degrees to the bedding (Figures 4(a)-(i)). Stereograms (stereographic projection) of orientations of discontinuities groups using stereonet ver. 8.9 software were done as shown in (Figures 5(a)-(j)). From stereographic projection, indicates
Figure 4. (a) The dominant fractures in station; (b) The dominate fractures in station; (c) The dominant fractures in station; (d) The dominant fractures in station; (e) The dominant fractures in station; (f) The dominant fractures in station; (g) The dominant fractures in station; (h) The dominant fractures in station; (i) The dominant fractures in station; (j) The dominant fractures in station.
that there are no conjugate systems of fractures in the same station to make use of it in finding the maximum principal stress. Therefore, tension set fractures (ac) were considered, which always accompanist to conjugate system and resulting from the same maximum stress, to find the maximum principal stress in the study area which represent the regional stress in the area. It was found that the maximum principal stress direction in the study area was between N15E to N40E.
Figure 5. Stereogram of orientations of discontinuities (fractures) groups in stations (1, 2, 3, 4, 5, 6, 7, 8, 9, and 10).
5.2. Discontinuity Density (d)
It is the cumulative trace length of discontinuities in a given area regardless of their orientation   (Figure 6).
In the present study, it was measured within an area about 2.5 × 2.5 m of each field station using the following equation:
d = discontinuity density
A = area of discontinuity density measurement (m2)
i = discontinuity number
N = whole number of discontinuities within A
L = discontinuity length (m).
It showed wide variation (Table 1) and it was used to evaluate the porosity and permeability of the exposed rocks of the study area.
Figure 6. Schematic block diagram shows the discontinuity spacing, persistence and aperture modified after   .
Table 1. The average readings of the discontinuities in the study area.
5.3. Discontinuity Spacing
Discontinuity spacing is the perpendicular distance between adjacent discontinuities  . Discontinuity spacing is one of the most important parameters because of its direct influence on rock slope stability and intact rock block size   . In the present study, the average spacing of each discontinuity set was measured at each station, and then the average spacing of discontinuities was computed for each station (Table 1).
5.4. Discontinuities Persistence
Discontinuities persistence is the trace length of discontinuities as observed on the surface of exposures and it is one of the most important rock mass parameters, because of its direct influence on rock slope stability and quarry. In the present study, the average persistence of discontinuity was measured for each station (Table 1).
5.5. Discontinuities Aperture
Discontinuities Aperture is the perpendicular distance between adjacent rock walls discontinuities    . The discontinuities aperture was measured for each discontinuities set. Next, the averages of apertures were found for all the observed discontinuities at each station (Table 1). The values of discontinuities aperture were increasing toward the center and decreasing on the top of the hills between the valleys.
Abo-Hederya is an arid area where wind is considered as principal geomorphological factor influencing rock slope geomorphology. Aeolian differential erosion acting on weak strata is forming overhanging unstable slopes.
The most common fractures have high dip angle (sub vertical) therefore they act as back release surfaces (BRS), lateral release surfaces (LRS) or composite back release surfaces (CBRS) but not as sliding surfaces. Their spacing and orientation also determine the size and shape of the failed rock blocks.
The high or low values of the discontinuities parameters have direct relationship with the structural position of the study area.
The slope stability analysis suggests that the different discontinuities sets have direct influence on slope failure in the study area. Nevertheless slope failure is limited in size distribution and has no much hazard on the quarry itself.
Fracture analysis reveals the geometrical and genetic relationship with the principal tectonic stresses, and it indicates that the maximum principal stress causing these fractures is directed N15 - 40E which is compatible with the regional stress resulting from Arabian plate motion towards NE. It has been found that tension fractures are common in the study area, while shear fractures are relatively less. This might be related to the structural position where the area is on the extension of the fold axis area.
The main stress in study area is compression and there might be some local stresses which have been repre- sented in the form of tension fractures set (ac) and shear fractures sets (hko, hol and hkl).
 ISRM (International Society for Rock Mechanics) (1978) Suggested Methods for the Quantitative Description of Discontinuities in Rock Masses. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 15, 319-368.
 Sens, Z. and Kazi, A. (1984) Discontinuities Spacing and RQD Estimates from Finite Length Scanline. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 3, 203-212.
 Henning, P.H., Olson, J.E. and Thompson, L.B. (2001) Combining Outcrop Data and Three Dimensional Structure Models to Characterize Fractured Reservoirs (and Example from Wyoming). American Association of Petroleum Geologists AAPG Bulletin, 84, 830-849.
 Cruden, D.M. and Hu, X.Q. (1994) Toppling on Underdip Slope in the Higwood Pass Alberta, Canda. Quarterly Journal of Engineering Geology, 27, 57-68.
 Evans, R.S. (1981) An Analysis of Secondary Toppling Rock Failures the Stress Redistribution Method. Quarterly Journal of Engineering Geology, 14, 77-86.