1. Introduction

Even though demand on water utilization has been increasing as a result of the increase in world population, industrialization and urbanization, the water resources potential stays constant. Furthermore, increase in agricultural activities, expansion of industrial and metropolitan urban areas and the effects of the climate change have severe adverse effects, decreasing land resources and putting pressure on water resources. Therefore, effective management of land and water resources becomes important. One of the key elements to achieve an effective water management is setting up technological infrastructures that provide gathering and storing of all relevant data on water resources. It is not easy to find appropriate data every time for a particular basin at the same quality and resolution. As such, remote sensing and GIS technologies gain importance for gathering data in the implementation of Integrated Water Resources Management (IWRM) plans. Implementation of IWRM effectively requires knowledge in basin dynamics, surface and groundwater resources potential of the basin and daily sectorial water consumption in the basin. Stakeholders in the basin should have an access to all those information related to the basin. In this context, many countries have collected basin- wide water data and established National Water Information Systems to store this data and share it with the public. European Environmental Agency together with the member countries is conducting a web based study on Water Information System for Europe within the context of the Water Framework Directive.

Implementation of IWRM effectively also requires a systematic definition of the watershed together with all lower level components. Basin, sub-basin definition should be considered within this context, and basin topologies should be established hierarchically. It is very important to define a unique code for each branch of the river network draining the basin. Even though studies related to basin and sub basin drainage system coding have been performed in many European countries, European Union still conducts studies known as “European Rivers and Catchments-ERICA” and “European Catchments and Rivers Network Systems-ECRINS” to provide a common coding base for determining accurately the basin and sub basin river network systems. In Turkey, determination of basin-sub basin drainage systems has been done by using a classical approach. With hydrologic perspective, Turkey was divided into 25 main basins, and some of the sub basins were determined during the water resources development pro- jects in these basins. The recent development of computer technologies and Geographic Information Systems (GIS) and use of these tools in water sources management pro- jects in those basins in Turkey helped engineers and water scientists. Various organizations in Turkey who are working in water sector, are using GIS technologies actively, especially in determining the study area’s hydrologic drainage characteristic for water resources project development. However, none of these studies systematically approached the determination of basin and sub basing drainage system coding.

Impacts of climate change affected the hydrology of watersheds adversely, and watershed authorities have been facing serious challenges in terms of reduction of water potential. Thus, water allocation among agriculture, hydropower, industry and domestic usage sectors became very important. Significant changes in magnitude and intensity of precipitation in basins caused problems in flood risk management. These challenges cause problems in water transmission between basins not only in one country but also for transboundary watersheds located in arid and semiarid regions. Thus engineers and water resources scientists have been working on these problems by using satellite information and GIS technologies to obtain solutions for these problems. Thus effective and accurate solution of water resources problems at basin scale requires the determination of watershed boundary, shape and size by determining river drainage network accurately and studying the drainage characteristics of the basin by using digital elevation models. Numerous studies may be found in literature for stream ordering [1] [2] . Horton provided the theoretical basis for quantitative analysis of drainage network structures [3] . He defined the main branch of a river discharging to a sea, lake or ocean as the first degree stream and all the branches joining to this main branch as the second degree streams, and the branches joining to the second degree streams as the third order streams and so on. In this definition the highest degree streams in the basin receive only the overland flow. The method that was defined by Horton for stream ordering has two stages [4] . In the first stage he defined the streams receiving only overland flow as the first degree streams. When two first degree streams join, they make up a second degree stream and when two second degree streams join, they make up a third degree stream. This ordering continues to the outlet of the basin.

In the method given by Strahler [5] - [7] a stream, which takes no other branch but overland flow as shown by number 1 in Figure 1, is called a first order stream. When two first order streams join, they make up a second order stream as shown by number 2 in Figure 1, a third order stream is produced when two second order streams are joined and can also take in other first and second order streams as shown by number 3 in Figure 1. This procedure is continued until the outlet of the basin and the order of the basin is equal to the order of the main stream (defined by number 4 in Figure 1) which is the longest branch in the basin. Shreve also defined a method for stream and basin ordering. Shreve’s definition of the first order branches is the same as Strahler’s method, but from then on branches take an order value, which is the summation of the orders of joining branches [1] . Strahler-Horton system is widely used among the methods

described above for ordering the streams and basins (Figure 1).

In order to design the structure of a river network topologically, the river network and its drainage area must be coupled so that the interaction between them can be established for coding. In general, Pfafstetter Coding System is used to achieve this work. This coding system gives information on natural channels, manmade irrigation and drainage canals, natural lakes, and reservoirs, shorelines, islands, and boundaries of the basin. In this system, branches of a river draining at most four major subareas are numbered evenly with single digit numbers, starting from downstream toward upstream. Inter-basins drained by the main branch are also numbered oddly with a single digit number starting from downstream toward upstream. The drainage area of the upstream branch numbered with 9 is always greater than the drainage area of the side branch numbered with 8 as seen in Figure 2 below. If there is a closed basin and/or a sub-basin, it is numbered with “0”. If the main river channel discharges to the sea, it is numbered with “1” [8] [9] .

Several countries derived their own river coding by using the Strahler-Horton system. United States Geological Survey developed “Hydrological Unit Code-HUC” for river basin coding. 1/250,000 scaled topographic maps were used during the first application of this coding. In this first application study, basins were categorized as Region, Sub-region, Accounting, and Catalog Units and each of these hydraulic units received HUC code consisting of 8 characters as “rr ss aa cc”. In this coding system “Region” is the first level hydraulic unit that includes main rivers or river systems discharging to a bay. “Sub-Region” is the second level hydraulic unit that includes main rivers, closed basins, and rivers discharging to the coastal zones. “Subregions” in the third level are divided into “accounting-units”. “Catalog Units” are the fourth level hydraulic units that include water collecting areas known as watersheds. This coding system was revised in 2010, and the new revised code consists of 12 characters that include six levels. In Germany, basin coding studies started in 1993 by a group named “Landerarbeitsgemeinschaft Wasser-LAWA” [10] . In this study, water contribution zones, consistent with hierarchical numerical numbers in the basin, were taken into consideration. As

a result of this study Germany was divided into six basins, namely the Donau with coding number 1, the Rheine with coding number 2, the Ems with coding number 3, the Wesser with coding number 4, the Elbe with coding number 5, and the Oder with coding number 6. As in the Pfafstetter coding system, 4 main river branches in each basin were considered in the coding studies. However, the direction of the coding study in LAWA system started from upstream toward downstream in contrast to pfafstetter coding system [10] . REGINE coding system is used in Norway. This system was developed by Norwegian Water Resources and Energy Administration in the context of Norwegian Water Information System. Rivers longer than 10 km, lakes greater than 1 km² surface area, reservoir areas and tunnel entrance locations were all considered in the REGINE system. In this system alphanumeric combination of coding was used for the determination of regions of the hydrologic drainage areas. Those regions were converted to main basins with single outlets [10] . Portuguese Institute for Water Resources also performed a study for river basin coding for those watersheds larger than 1 km². Portuguese hydrologists studied basin coding of various European basins and decided to use pfafstetter coding system for their river basins. This coding system may also provide a basis data for the European Water Framework Directive studies [11] [12] . Romania used Gravelius method for coding the Danube River which is the main river of Romania [13] . European Environmental Agency (EEA) conducted a feasibility study to develop a database of the boundaries of the transboundary rivers that are of common interest to various nations of Europe. Based on this feasibility study of EEA, low resolution maps with 1/100,000 scale and medium resolution maps with 1/25,000 scale were used in a pilot study to form a basin database to provide for the urgent data needs of the European nations. A river basin coding system known as ERICA-CS was developed within the framework of this project. This coding system was developed by using 1/25,000 scale river network that is based on pfafstetter coding system. This coding system is named as modified pfafstetter coding system and it is defined as “MM BBB N1 N2 N3 N4 A”, where MM is the marine code, BBB is the marine border code, N1 N2 N3 N4 is the nested catchment code, and A is the area band [14] . The modified pfafstetter coding system was modified in 2012 by including a country code at the beginning of the coding and also by increasing the nested catchment code from 6 digits to 22 digits so that each of the stream levels are incorporated into the coding [15] .

In 1999, Joint Research Center of the European Commission launched a study on a database for basin boundaries and river networks of Europe within the context of “Euro Landscape” Project, and by using 1/250,000 - 1/500,000 scaled maps. Findings of this study were published under the title “Catchment Characterization and Modelling” (CCM). Later on, Space Shuttle Topography Mission-SRTM data were combined with the Catchment Characterization and Modelling database to obtain the Digital Elevation Maps of Europe, and CCM study was revised as CCM2 [16] . “European Catchments and Rivers (ECRINS)” was developed by incorporating Corrine Land Cover and Water Framework Directive reports into CCM2 in 2012 [16] - [18] .

A study was also performed to develop hydrological river basin coding for Turkish river catchments based on the evaluation of the ECRINS coding system, briefly discussed above. The developed hydrological river coding is flexible and versatile for all types of basins in Turkey. Using this hydrological river basin coding will define river basin boundaries accurately, prevent conflicts in sectorial water allocation caused by uncertainty in locations of water source, produce a base for determining water potential and estimating extreme hydrological events of the basin, produce a basis for the prevention of water disputes among stakeholders within the basin, and help implementation of IWRM at basin level.

2. Methodology for the Derivation of the River Basin Coding in Turkey

Hydrological river basin coding developed in this study is based on the European Catchments and River Network System-ECRINS, and is defined in the Table 1 given below.

The priority matters for the development of basin and sub basin studies for rivers and basins in Turkey are as follows: Hydrological River Basin coding developed for Turkish river basins should meet the needs of the reports requested by the European Union, should easily be adoptable to other coding studies, and should protect areas and hydrological structures of the existing 25 river basins and sub basins of Turkey. The general form of the code that covers all these points is given in Table 1. In this table the country code MS is defined and published by “ISO 3166-1-ALPHA-2”. For example TR is used for Turkey. In Table 1 code BD defines the river basin district code shown in Figure 3, which is a special code for Turkey especially from the management point.

Even though Turkey is divided into 25 river basins, this division does not show the features of actual hydrological river basins as defined under the Strahler-Horton system. Those 25 basins may not only cover the major rivers but may also cover small basins in the coastal zones whose creeks discharge directly to the seas. Gediz, Buyuk Menderes, Kızılırmak, Yesilırmak, Sakarya, and Ceyhan river basins are this type of basins. In order to include the small basins which is actually do not contribute their discharge to the major river of the basin, River Basin District (RBD) definition is adopted (Figure 3) and included in the coding. Thus, small creeks that discharge directly to the Marmara Sea, the Eastern and Western Black Sea, and the Mediterranean sea regions are grouped to form one basin in the respective regions. This approach is also adopted in closed basins such as the Konya, the Lake Van, the Akarcay, and the Burdur-Beyse- hir Lakes closed basins. This definition is also consistent with the European Water Framework Directive.

On the other hand, Turkey has transboundary basins where the river originates in the Turkish territory and discharges to the seas in the neighboring countries. The Euphrates, the Tigris, the Aras, and the Çoruh rivers originate from Turkey, and then cross from the boundary of Turkey and discharge to the sea within the neighbor country. The Maritza river (Meriç) and the Orontes river (Asi) originate from the neighbor country and discharge to the sea in Turkey. Portions of these transboundary river basins that are located within the Turkish territory are considered and adopted to the River Basin District (RBD) definition (Figure 4).

In Table 1, the code H defines Hydrological System Number. This number defines whether the major river of the basin is discharging to the Mediterranean Sea, Black Sea, Persian Gulf, Caspian Sea or inland lakes within Turkey. Alphabetic capital characters are used for defining hydrological system code. For example, code “M” is used for basins that contribute their discharge to the Mediterranean Region. Studies done for the development of water resources projects for the 25 Basins in Turkey are all evaluated, and those 25 basins are named as River Basin Districts in Turkey is shown in Figure 4 and Table 2.

Results of hydrological studies done in these basins showed that some of the streams and creeks considered within the Euphrates-Tigris River Basin are actually located in the Lake Orumiyeh Closed Basin and in the Orontes River Basin. Moreover, some of the creeks considered in the Lake Van Closed Basin are actually within the Aras River Basin. In order to be consistent with the previous studies done for those 25 basins, same names were used in the River Basin District definition developed in this study. However, hydrologists must be aware that these basins are not defined by Strahler-Horton systems, but are defined as River Basin Districts (RBD).

It is important to define whether a river is discharging to a sea or an inland water body (lake) for the coding system developed for Turkish river basins in this study. A two-digit code is used to be consistent with “ECRINS” coding system. According to the “ECRINS” coding system, rivers and streams do not necessarily discharge to the same sea or lake, as in the case of The Marmara Basin and The Konya Closed Basin. Thus rivers discharging to different water bodies are coded with the codes defined for those water bodies. Table 3 defines the codes of those water bodies used in this study.

Even though their codes are taken from “ECRINS” coding system, the codes of the Marmara Sea, the Akşehir Lake, the Burdur Lake, the Acıgöl, the Hotamış Lake, and the Akgöl which are not defined in the “ECRINS” Coding system, are integrated into the “ACRINS” Coding system in this study as shown in Figure 5 and Table 4. As shown in Table 4 rivers in the Marmara River Basin District are discharging to three different seas with different codes.

After the determination of Country Code, Basin District Code and Hydrological System Code, the remaining part of the coding consists of the coastal and sub basin codes. Coastal code is an important variable used in this study and the Pfafstetter Coding System was used to define this variable in the coding system. This coding for the coastal areas follows the procedure where four largest river basins on the Black Sea coast of Turkey, discharging to a sea, are taken and coded with even numbers as 2, 4, 6 and 8. Areas situated in-between those four large basins, known as inter basins, are coded with

odd numbers as 1, 3, 5, 7, and 9. Basins coded with even numbers are the main river catchments which can also be divided into smaller sub basins and coded by using Pfafstetter Coding System (Figure 6).

As can be seen from Figure 6, large river basin areas in the eastern Black Sea zone are coded as TR22.M5.2 and the inter basins are coded with TR22.M5.3. These areas can also be divided into subareas and coded by using Pfafstetter coding system. For example, in closed basins if the streams are discharging to a natural lake, the coding of the drainage areas is done in the clockwise direction according to Pfafstetter system (See Figure 7).

3. Application of the Method to Turkish River Basins

In order to show the flexibility and versatility of the coding system, as defined in Tables 1-4, it was applied to one of the major river basins in Turkey, the Yesilirmak River Basin, which has all the complexities in terms of small sub-basins that do not contribute their flows to the main river branch in the basin. This basin is located in the coastal region of the Black Sea in the Northern sector of Turkey (Figure 8). This basin takes its name from the main river which is called the Yeşilirmak River (Figure 9).

Figure 9 depicts the drainage network of the Yeşilırmak River basin at six different degrees. The degrees of the drainage system are shown by different colors that are determined by the Pfafstetter coding system. According to the Pfafstetter coding system, sub basins of this basin were determined as follows (please see Figure 10). The main drainage area of the Yeşilirmak River is defined by TR14.M5.2. The main subareas located at the outlet of the Yeşilirmak River Basin and discharging directly to the Black Sea, are coded by TR14.M5.1 at the right side of the Yeşilırmak River outlet and TR14.M5.3 at the left side of the Yeşilırmak River outlet. The main drainage area, coded by TR14.M5.2, is further divided into subareas starting from the downstream to the upstream of the Yeşilırmak River. Sub-basins that have the largest areas are determined as the main sub-basins and are coded by even numbers 2, 4, 6 and 8. The remaining sub areas depicting inter-basin characteristics are coded by odd numbers 1, 3, 5, 7 and 9 as shown in Figure 10. Based on the Pfafstetter coding system, the whole Yeşilırmak river basin is divided into 9 sub-basins.

Even though the dashed circle area, shown in Figure 10, is located at the left side of the Yeşilirmak River Outlet, and is part of the Yeşilırmak River Basin, it does not contribute its discharge to the Yeşilırmak River. This area discharges directly to the Black Sea as shown in Figure 11. The same approach, defined above for coding of the sub-basins, is adopted for determining the codes of the sub-basins and inter-basins of this area. Figure 11 also shows how this section of the Yeşilırmak River basin is coded so that even though it is considered as part of the Yeşilırmak River basin by the district code 14M5, the discharge of the Yeşilırmak River will not be influenced at all when a water resources development project is planned for this section of the basin. The remaining part of the Yeşilırmak River basin contributes all its discharge to the Yeşilirmak River flow.

The second example application of the developed coding system is adopted to the Akarcay closed basin shown in Figure 12. The Akarcay river is one of the important inland water resources of the Afyonkarahisar province. This river originates from the north of the city of Afyonkarahisar, then drains all the areas Afyonkarahisar province and discharges to the inland Eber Lake first. It then hydraulically connects the Eber

and the Akşehir inland Lakes as shown in Figure 12.

With the above-given definitions, the coding system developed for Turkish river basins works accurately and is applicable to all kinds of basins, such as the closed basins (Lake Van, and Konya), and Transboundary basins (Meriç (Maritza), Çoruh (Chorokhi) and Tigris-Euphrates).

4. Discussion

The developed coding system was applied to all 25 river basins in Turkey to determine whether it is applicable to the delineation of the watershed boundaries accurately provided that digital elevation map (DEM) information with high resolution is available. This coding, developed for Turkish river basins, is also appropriate for the management of Turkish water resources by State Hydraulic Works (DSI) which is the main water authority in Turkey. Turkey is divided into 25 hydrological river basins where at each basin Regional Water Authorities (DSI Regional Directorates) are established in order to prepare and implement basin-wide water resources development projects. To prevent conflict in sectorial water allocation and water disputes among the stakeholders in the basin, it is important to determine the watershed boundary along with its surficial drainage system. The coding system developed in this study provides a base for the solution of such problems. For example, if a project is developed for the Yeşilırmak River basin to supply water for agriculture, environment, industry, hydropower, and domestic use, after allocating all the water needs of the different sectors, baseflow of the river should always drain into the river channel. The stakeholders living at the very upstream of the Yeşilırmak River basin which is not contributing its flow to the Yeşilırmak River, may not be influenced by the impacts of the project on the Yeşilırmak River, because this part of the basin has its own water supply and drainage system even though it is considered to be within the Yeşilirmak River basin. The appropriate sub-basin delineation analysis can be done precisely by the help of the coding system described above. The coding described above not only prevents conflicts in sectorial water allocation caused by uncertainty in the locations of the water sources, but also provides a base for determining the water potential and estimating the extreme hydrological events of the basin. Similar analysis can also be done for closed basins and for transboundary basins easily.

5. Conclusions

GIS technologies have been extensively used in hydrological studies during the last two decades. Using GIS in the determination of watershed boundaries enables the water resources planners, hydrologists, engineers and academicians to determine water potentials, to estimate extreme hydrological events, to allocate water for agriculture, hydropower, industry and domestic uses, for water transmission between watersheds and for management of transboundary basins during floods and draughts.

A study was performed to develop hydrological river basin coding for Turkish river basins. In this development, firstly ECRINS and pfafstetter methods were combined and secondly River Basin District concept was incorporated into the coding so that the sub areas which are parts of the basin but do not contribute their discharges to the major river of the basin, are accounted within the system. The developed hydrological river coding system is flexible and versatile for all types of basins in Turkey. In order to demonstrate the flexibility and accuracy of the coding system, it was applied to three river basins in Turkey each of which have complex hydological and topographical features that are different from each other. These basins are the Yeşilirmak River basin which is one of the major river basins discharging to the Black Sea, the Maritza (Meriç) River basin which is a transboundary river basin, and the Akarçay Closed Basin that is drained by the Akarçay River that discharges to the inland Eber and Akşehir lakes. These applications showed that using the developed hydrological river basin coding defines river basin boundaries accurately. The accurate delineation of basin boundaries is helpful in the planning of a water resources development project, eliminates or minimizes uncertainty in sectorial water allocation caused by locations of water sources, produces a base for determining the water potential and estimating extreme hydrological events of the basin, produces a basis for the prevention of water disputes among stakeholders within the basin, and helps the implementation of IWRM at basin level.