1. Introduction
Since the last decades of the twentieth century, Remote Sensing (RS) has become a popular approach for detecting LULC changes. RS imagery has been universally applied in urban expansion and environmental analysis, with a large number of functional and sophisticated strategies being developed [
1]. LULC classification plays an essential role in research on the physical and social sciences concerning urbanization and soil mapping [
2]. LULC is an essential parameter for detecting human activities and natural changes due to climate change. LULC change detection is valuable for a broad variety of subjects, such as urban development, coastal dynamics, deforestation monitoring, shoreline changes and river transportation using multitemporal satellite images [
3,
4,
5]. LULC changes are affected by human beings seeking to meet different needs, for instance, residential, agricultural, industrial, mining and other infrastructural facilities, and as such are main concerns connected with the economic and sustainable development of a region [
6,
7]. The timely and precise detection and identification of changes in the Earth’s surface provide a foundation for a deeper understanding of human-nature interactions and correlations [
8,
9,
10,
11,
12,
13]. Formerly, in situ data had been used for studying LULC, although such kinds of data are inadequate in a spatial and temporal extent. It is also hard, if not impossible, to collect data from unreachable areas. Thus, data collection is drawn-out and often expensive. Therefore, such data are not dependable for effective planning and decision-making [
14,
15]. Satellite images are a valuable source of information on LULC as they provide swift, cyclic and precise data acquisition from the RS system [
16]. Multitemporal RS data, mainly time-series Multispectral Scanner (MSS), Thematic Mapper (TM), Enhanced Thematic Mapper Plus (ETM+) and Landsat Data Continuity Mission (LDCM) images, have been extensively applied for detecting LULC changes since the launch of the Landsat satellite in the early 1970s [
17,
18,
19]. The reason for the wide use of Landsat images by researchers is because they are relatively inexpensive, more accessible and have a long history.
The expansion of built-up areas is widely known as a means of establishing an unacceptable climate in terms of air quality and controlling natural resources. This is because urban sprawl is turning the natural land cover into man-made land cover to meet human needs [
20]. Investigating spatial and temporal forms of urban development has become one of the most critical challenges in monitoring and analyzing present and future directions of urban growth issues [
21]. Over the last 300 years, the global trend of LULC has been characterized by deforestation and the expansion of agricultural land [
22]. However, during the last decades, developing countries have experienced a reduction in rural land and an expansion of urban areas through urbanization [
23,
24,
25,
26]. In many respects, urban growth phenomena have become unsustainable, in the majority of cities globally [
27]. In fact, urbanization itself is a hot topic of concern globally, where people are leaving rural areas and settling in big cities [
28].
The increase in impermeable areas in a catchment is known to boost flood risk [
29]. Urbanization is considered to be the major cause of transformations in hydrologic and hydraulic processes, decreasing drainage capacity and increasing the number of flash floods in urban areas [
30,
31]. Consequently, it increases the total runoff volume and peak discharge of storm runoff events [
23,
32,
33,
34,
35,
36]. Therefore, it is essential and useful to understand to what degree, for what purpose, and at what place LULC changes occur. Such radical changes in LULC have drawn the attention of researchers more and more to determine the extreme effects of these transformations on several directions of urban planning and the environment, including strategic land management, air quality requirements, and the probability of decreasing flash floods [
37,
38,
39,
40,
41,
42]. Many papers have investigated the effects of LULC changes on urban flash floods. Sanyal et al. [
34] regarding the Konar catchment in eastern India, showed that the increase in peak discharge for different LULC scenarios, like 1976 to 2004, was 1023.3 to 1194.7 m
3/s with a shorter time to peak, at one hour and 10 min. Al-Aysen, Riyadh, in Saudi Arabia, after a 15% increase in urbanization in the mentioned watershed, caused nearly a 30% step-up in the peak discharge for the five-year return period event. An increase of 115% appeared when the basin became entirely urbanized for the five-year return period event [
43]. Wałęga et al. [
29] determined the effect of catchment impermeability on flow variability in the Sudół Dominikański stream, Cracow, Poland and their investigation showed an increment in maximum flows over time connected to an increase in catchment imperviousness.
The KRI, from 2003 to present, witnessed its most significant development and expansion. The LULC patterns of Erbil, the capital city of the KRI, have changed extensively in the last two decades. These changes are because of the growth of the economy, and the region being comparatively safer than middle and southern Iraq. There is still a high demand for buying properties. Moreover, the Kurdistan Regional Government (KRG) encourages local and international investors to build residential and industrial areas. At the same time, the probability of the occurrence of urban flash floods in the center of Erbil Province is also increasing accordingly. Hameed, in his work that concerned the central district part of Erbil Province, showed that between 1984 and 2014, built-up areas and bare land expanded, while areas of vegetation and farmlands decreased dramatically [
44]. Since 2010, the number of urban flash floods has increased, and as a consequence, many sewer pipes overflowed, which led to sewage escaping and mixing with stormwater. Flooding also caused more power outages as many local and small power generators were flooded [
45].
During the last decades, the LULC of Erbil Province has significantly changed. Thus, the resulting changes need to be mapped urgently in order to recognize the direction of LULC transformation in the different parts of the city to assist in any possible land development planning, flood control and resource management. Assessing, identifying and classifying LULC changes in Erbil Province, and particularly how rapid urban development affects the center of the city, would contribute to managing the hydro-structures effectively and monitoring environmental changes. Additionally, it would help authorities in Erbil municipality and the General Directorate of Water and Sewerage to identify priorities, as well as helping them in making the right decisions. In this paper, we examined changes in LULC in the studied area in five phases. The time-steps started in 1984, 1994, 2004, 2014 and 2019. The studied area is located in Erbil-KRI (North of Iraq). We attempted to find answers to the four research questions concerning LULC changes and their influence on the runoff coefficient (the runoff coefficient is a unitless determinant that is used to convert the amounts of rainfall to runoff, or the ratio of the peak runoff rate to the rainfall intensity [
46]) in the studied area from 1984 to 2019:
(1) How was the LULC changed during this period, after every ten years?
(2) How did the economic, social and political situation affect the LULC changes?
(3) How did this transformation affect the hydrological response in the studied area?
(4) Is the probability of urban flash floods increasing?
4. Discussion
RS data, by employing the maximum likelihood method for image classification, were used to provide detailed information to describe the LULC dynamics around the central district of Erbil Province. Regarding the percentage of the accuracy assessment of the classified images (85.60%), this could be considered to be an acceptable classification. A Kappa statistics value of between 0 and 1 for the analysis of classified maps has become the standard for the agreement of maps in practice [
65]. The Kappa values for the results of our study and the validation of the visual interpretation classification of five classified images were over 0.79, which means that all maps agree well (a Kappa value between 0.61 and 0.8 is considered to correspond well [
66,
67]).
Urbanization in Erbil city has witnessed both dramatic changes and shifts in different periods. These changes principally corresponded to legislation factors, in parallel with the political situation, and economic and social determinants over four subsequent periods: (1920–1958), (1958–1991), (1991–2003), (2003–present) [
68]. Between 1984 and 2003, Iraq generally and the Kurdistan Region, which was established in 1992, were in a period of slow economic growth, due to the number of wars that happened in this area, such as the Iraq–Iran war (22 September 1980–20 August 1988), Gulf War (2 August 1990–28 February 1991), and Kurdish conflict (1994–1998). However, the Kurdistan Region, and especially Erbil as the capital of the region, witnessed considerable growth and development after 2004. Generally, since 1991, the KRI has been a de facto independent entity making its own regulations and policies, and according to the new constitution, which was prepared in 2005, Iraq is a federal country, with the Kurdistan Region as the only federal state within its borders [
69]. Further analysis showed that the increase in built-up areas was about 25 km
2 between 1984 and 2004, which nearly equals 1.25 km
2 per year. Furthermore, only between 2004 and 2014, the increase in built-up areas was about 63 km
2, which is the equivalent of 6.3 km
2 per year, which can be considered as very high, indicating rapid urban growth in the city with a spatial and temporal change in land use over only ten years. Oil production and real estate investment caused an economic boom. The population of Erbil Province grew from approximately 95,000 inhabitants in 1965 to about 1.3 million in 2010 [
70], recently reaching more than 2 million.
Figure 8 shows the evolution of built-up areas from 1984 to 2019. The current study, supported by previous research in this area, demonstrates that increasing urban development in flood plain zones will boost peak discharge, reduce the time to peak, and increase the runoff volume [
2,
29,
71,
72,
73,
74].
The findings of this study indicate that urbanization resulted in an increased risk of urban flash floods. Another cause of this risk is also the changing characteristics of precipitation as a result of global climate change [
47,
63]. Such floods are dangerous to people and destroy the urban infrastructure due to the overflow of the stormwater drainage system [
75] and increase erosion of the land surface [
76]. Generally, urban development will not only boost peak discharge but will also elevate the runoff volume and the extent of the flooded zone within a catchment. Because Erbil city center is more or less a flat territory, the flooding will last for a longer time compared to steeper areas, as well as due to the circularity of the main roads around Erbil Citadel, such as Barzani circular street, Kurdistan circular street, Peshawa circular road, 120 m ring road and the under-construction 150 m ring road. All these mentioned roads cross rivers and streams passing through the city. As a consequence, the routes of most of the rivers and streams have been either blocked or changed to smaller sewer culverts or pipes. This means for urbanized areas, streets, underpasses and low elevation areas, the possibility to be flooded is greater compared to other areas. Erbil center district has suffered from urban flash floods since the 1930s. In the last two decades, the number of flash floods due to heavy rainfalls has increased. Recently, on February 20, 2020, some neighborhoods in the Erbil center district were flooded. According to the meteorological station of the Erbil Directorate of Irrigation, the amount of 24.3 mm rainfall was measured within six hours. Moreover, the highest rainfall intensity recorded in this episode was 9.8 mm/hour in the middle of the event. Water accumulated in the expected areas within a short time. The inundation level was more than 25 cm in some areas, as shown in
Figure 9a,b.
5. Conclusions
This study investigated both the dynamic changes in LULC and evolution in the hydrological response in the center of Erbil Province, KRI, from 1984 to 2019 using multispectral satellite data. To sum up our research questions from the introduction section, our outcomes are presented as follows:
The results revealed that between 1984 and 2004, bare land and built-up areas steadily increased, while between 2004 and 2019, built-up areas soared by almost 245% (from just 55.823 to 136.658 km2). Permeable areas like agricultural land and vegetation steadily decreased from 1984 to 2019. Only between 2014 and 2019 did vegetation areas increase, especially in the downstream area of the city. In our opinion, these changes are related to farming types in these areas, and the remote sensing data showed the class as vegetation;
Three factors (economic, social and political) influenced the LULC changes in the studied area. For instance, after 2004, the KRI generally, and especially Erbil, enjoyed a decade of prosperity and development due to many factors such as stable security, local and international investment, and the successful implementation of some of the American free market-based developmental policies;
The influence of urbanization on the flood peak discharge, runoff volume and runoff coefficient were investigated for the depth of storm rainfall at the probability distribution of 10%, which is equal to 71.16 mm under different urbanization scenarios. The simulations showed that urban development could considerably aggravate flooding caused by a given storm due to the hindrance of natural drainage and decreasing permeability;
As it is obvious in the analysis, parallel to urban development, the runoff volume and peak discharge increased as well. Therefore, the probability of more flash floods in the city is increasing likewise. If the authorities do not take steps related to adopting a new strategy in order to stop urban flash floods, such flash floods will target areas that are denser in urban development and population.
The outcome of this study shows how important remote sensing is for the detection of previous, present and future human activities, in a relatively simple, cheap and easy way. Moreover, without using RS data and GIS techniques, the analysis of LULC in third world countries is somehow difficult or impossible. The results achieved from the study indicate that there is an increase in the runoff peak and volume, which is due to the hydrological response related to urbanization in the city. As a consequence of this, the possibility of flash floods is increased and it becomes a risk to life quality in the studied area, especially in the center of the city. The findings of the study are also important for policymakers and local authorities when deciding how to re-adapt the city to the current situation and plan for future development. In addition, the findings and the prepared LULC scenarios of this study will be useful in the research field on the flash floods mitigation in the studied area.