SARchaeology—Detecting Palaeochannels Based on High Resolution Radar Data and Their Impact of Changes in the Settlement Pattern in Cilicia (Turkey)

: The fertile alluvial plain of Cilicia is bordered by the Taurus and Amanus mountain ranges to the west, north and east and the Mediterranean Sea to the south. Since the Neolithic Period, Plain Cilicia was an important interface between Anatolia and the Levant. The alluvial plain is dominated by three rivers: Tarsus, Seyhan and Ceyhan. The avulsion history of the lower course of the rivers Seyhan and Ceyhan during the Holocene remains an unresolved issue. The knowledge about how former river courses have changed is essential for the identiﬁcation of ancient toponyms with archaeological sites. The analysis of silted up riverbeds based on high-resolution digital elevation models (TanDEM-X) and historic satellite imagery (CORONA) in this paper provide the ﬁrst indications for the reconstruction of former river channels. Further evidence is given by the evaluation of the settlement patterns from 3rd to 1st millennium BC.


Introduction
The changes of river courses, which have had a strong influence on the development of sites and is important for the localization of sites has not been investigated in Cilicia in detail so far.Only the development of river deltas and the changes of the river Tarsus have been studied in detail [1][2][3][4].For the identification of toponyms known from written sources with archaeological sites, precise knowledge about the changes of river courses is essential.The aim of this paper is a first look at the development of the Seyhan and Ceyhan rivers during the late Holocene.Due to the lack of available coring data, this will be done by the analysis of remote sensing data in combination with a study of changes of the settlement cluster in Plain Cilicia.In southern Mesopotamia, studies about ancient watercourses on which all communication moved and along which all settlement ranged itself were undertaken already in the 1960s and 1970s on the basis of archaeological data and remote sensing data (LANDSAT) [5].Since then, much more data has become available e.g., CORONA and TanDEM-X.In this study, these archives will be analyzed for the region of Plain Cilicia.
Plain Cilicia (gr.Kilikia Pedias, lat.Cilicia Campestris) covers an area of approximately 8000 km 2 in the northeastern corner of the Mediterranean Sea.Like almost no other landscape in Asia Minor, this plain is defined by strong topographic contrasts: a large fertile plain surrounded by steep mountain ranges.The Middle Taurus Range (turk.Orta Toroslar), which exceeds 3700 m in elevation, borders the landscape to the west and north.In the east, the settlement cluster is confined by the Amanos (turk.Nur Da gları) that separates Plain Cilicia from inner Syria and reaches elevations of up to 2200 m.The Gulf of Iskenderun and the Gulf of Mersin form its southern boundary (Figure 1).Natural passes through the mountains still form the main routes of access to the lowlands.The "Cilician Gates" (turk.Gülek Bo gazı) at an elevation of 1290 m connect Anatolia with the plain and serve as the main passageway for railroad tracks and major vehicular roads up until this day.In the westernmost part of Plain Cilicia, forming the transition to Rough Cilicia, the Göksu Valley stretches north-south from the Konya Plain to coastal Silifke and represented an important connection to Anatolia in antiquity [6] (p. 88).The Bahçe Pass ("Amanian Gates") in the East leads over the Amanus Range and serves as a passage between Cilicia and the Islahiye Plain.Further south, the Orontes Valley can be reached over the Belen Pass ("Syrian Gates").Access from the north is facilitated by several tributary valleys of the Seyhan and Ceyhan rivers.Late Bronze and Iron Age rock reliefs as well as inscriptions like Hanyeri, ˙Imamkulu, Taşçı and Fraktın [7], mark their significance in antiquity.The plain itself is divided into a western (Çukurova) and an eastern (Yukarıova) settlement cluster by a natural border formed by the foothills of the Middle Taurus Range and the Misis Mountains.The approximately 760 m high Misis Mountains [8] (pp.305-307) extend northwards from the western end of the Gulf of Iskenderun (Figure 1).The average height of the Çukurova is about 15 m above sea level.The height of the Yukarıova varies from 20-100 m above sea level.
elevations of up to 2200 m.The Gulf of Iskenderun and the Gulf of Mersin form its southern boundary (Figure 1).
Natural passes through the mountains still form the main routes of access to the lowlands.The "Cilician Gates" (turk.Gülek Boğazı) at an elevation of 1290 m connect Anatolia with the plain and serve as the main passageway for railroad tracks and major vehicular roads up until this day.In the westernmost part of Plain Cilicia, forming the transition to Rough Cilicia, the Göksu Valley stretches north-south from the Konya Plain to coastal Silifke and represented an important connection to Anatolia in antiquity [6] (p. 88).The Bahçe Pass ("Amanian Gates") in the East leads over the Amanus Range and serves as a passage between Cilicia and the Islahiye Plain.Further south, the Orontes Valley can be reached over the Belen Pass ("Syrian Gates").Access from the north is facilitated by several tributary valleys of the Seyhan and Ceyhan rivers.Late Bronze and Iron Age rock reliefs as well as inscriptions like Hanyeri, İmamkulu, Taşçı and Fraktın [7], mark their significance in antiquity.The plain itself is divided into a western (Çukurova) and an eastern (Yukarıova) settlement cluster by a natural border formed by the foothills of the Middle Taurus Range and the Misis Mountains.The approximately 760 m high Misis Mountains [8] (pp.305-307) extend northwards from the western end of the Gulf of Iskenderun (Figure 1).The average height of the Çukurova is about 15 m above sea level.The height of the Yukarıova varies from 20-100 m above sea level.Three main rivers, originating in the Taurus Mountains, cross the lowlands and discharge into the Mediterranean Sea: Berdan or Tarsus Çayı (gr.Kydnos), Seyhan (hitt.Šamri/Sapara, gr.Saros) and Ceyhan (hitt.Puruna, gr.Pyramos) [10] (p. 6).These rivers form the deltaic complex that occupies the Three main rivers, originating in the Taurus Mountains, cross the lowlands and discharge into the Mediterranean Sea: Berdan or Tarsus Çayı (gr.Kydnos), Seyhan (hitt.Šamri/Sapara, gr.Saros) and Ceyhan (hitt.Puruna, gr.Pyramos) [10] (p. 6).These rivers form the deltaic complex that occupies the Adana Basin (Figure 2).The Adana Basin represents the onshore extension of the marine Cilicia Basin [11] (p. 123).The 6000 m thick sedimentary succession of the Adana Basin is spanning Miocene to Recent [12] (p.190) and is built by alluvial variable grainsize deposits [13] (p. 135).All rivers are braided in their upper reaches with channel gradients ranging from 1:40 to 1:400, while the channel gradients decrease to >1:2000 in the lower 35-110 km of their courses, causing them to meander [14] (p.57).The major tributaries of the 560 km long Seyhan are the Göksu, Zamantı, Çakıt and Körkün Rivers [15] (p.647).Prior to construction of the first dam in 1954 the Seyhan River Delta expanded its area on average by 2.8 ha per year [16] (p.64).The dam greatly reduced sedimentation in the delta [17] (p.5).The 510 km long Ceyhan River drains today into the Gulf of Iskenderun and formed a more than 100 km 2 large river delta.Paleoenvironmental studies of the Iskenderun Bay show a muddy sedimentation, which coincides with a major delta progradation of the Ceyhan around 2140 BP [18] (p.20).Therefore, a drainage into the Mediterranean Sea at Kap Karataş during the Iron Age is generally accepted [18] (p.20).
Geosciences 2017, 7, 109 3 of 20 Adana Basin (Figure 2).The Adana Basin represents the onshore extension of the marine Cilicia Basin [11] (p. 123).The 6000 m thick sedimentary succession of the Adana Basin is spanning Miocene to Recent [12] (p.190) and is built by alluvial variable grainsize deposits [13] (p. 135).All rivers are braided in their upper reaches with channel gradients ranging from 1:40 to 1:400, while the channel gradients decrease to >1:2000 in the lower 35-110 km of their courses, causing them to meander [14] (p.57).The major tributaries of the 560 km long Seyhan are the Göksu, Zamantı, Çakıt and Körkün Rivers [15] (p.647).Prior to construction of the first dam in 1954 the Seyhan River Delta expanded its area on average by 2.8 ha per year [16] (p.64).The dam greatly reduced sedimentation in the delta [17] (p.5).The 510 km long Ceyhan River drains today into the Gulf of Iskenderun and formed a more than 100 km 2 large river delta.Paleoenvironmental studies of the Iskenderun Bay show a muddy sedimentation, which coincides with a major delta progradation of the Ceyhan around 2140 BP [18] (p.20).Therefore, a drainage into the Mediterranean Sea at Kap Karataş during the Iron Age is generally accepted [18] (p.20).The course of the coastline underwent major changes.There are a few studies for a reconstruction of the palaeoshorelines available [2,4,14,17,18,[20][21][22][23].A sea level high stand around 6000-5000 BP is generally accepted, followed by a significant regression of the sea level of several meters as well as a slight rise to the modern datum [24].A mid-Holocene coastline about 5-10 km inland from the present-day coastline is assumed [1,4].
Also from a historical perspective, Plain Cilicia is defined by contrasts between independence and foreign influence up to annexation.During the Chalcolithic Period, there is strong evidence for connections to the Amuq (Figure 1) and other Levantine and north Mesopotamian regions [25].At the beginning of the Bronze Age, end of the 4th millennium BC, an alteration is visible by a change in the ceramic assemblage and a first orientation to the Anatolian plateau in the north [26] (pp.141-142), [27].However, there are still parallels to the Amuq region and the Northern Levant.At the beginning of the 2nd millennium BC, Syro-Cilician Pottery is predominant [28] (p.11).When Cilicia became the The course of the coastline underwent major changes.There are a few studies for a reconstruction of the palaeoshorelines available [2,4,14,17,18,[20][21][22][23].A sea level high stand around 6000-5000 BP is generally accepted, followed by a significant regression of the sea level of several meters as well as a slight rise to the modern datum [24].A mid-Holocene coastline about 5-10 km inland from the present-day coastline is assumed [1,4].
Also from a historical perspective, Plain Cilicia is defined by contrasts between independence and foreign influence up to annexation.During the Chalcolithic Period, there is strong evidence for connections to the Amuq (Figure 1) and other Levantine and north Mesopotamian regions [25].At the beginning of the Bronze Age, end of the 4th millennium BC, an alteration is visible by a change in the ceramic assemblage and a first orientation to the Anatolian plateau in the north [26] (pp.141-142), [27].However, there are still parallels to the Amuq region and the Northern Levant.At the beginning of the 2nd millennium BC, Syro-Cilician Pottery is predominant [28] (p.11).When Cilicia became the Kingdom Kizzuwatna independent, mid 2nd millennium BC, there is a sudden change in the material evidence, with a strong orientation to the Anatolian plateau visible.From mid-15th to the beginning of the 12th century BC, Kizzuwatna stood under Hittite or Mittani supremacy [28] (p.36).At the transition from the Late Bronze to the Early Iron Age Plain Cilicia became again independent as the kingdom Geosciences 2017, 7, 109 4 of 20 Kingdom Kizzuwatna independent, mid 2nd millennium BC, there is a sudden change in the material evidence, with a strong orientation to the Anatolian plateau visible.From mid-15th to the beginning of the 12th century BC, Kizzuwatna stood under Hittite or Mittani supremacy [28] (p.36).At the transition from the Late Bronze to the Early Iron Age Plain Cilicia became again independent as the kingdom Ḫiyawa/Que [29].During the 8th century BC Plain Cilicia was integrated into the Assyrian Empire [30] (pp.[7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24].The knowledge of the material culture of Kizzuwatna and Ḫiyawa/Que is still fragmentary.Most information come from excavations of the main Bronze and Iron Age sites in the plain.These are from West to East: Mersin-Yumuktepe, Tarsus-Gözlükule, Adana-Tepebağ, Sirkeli Höyük, Tatarlı Höyük and Kinet Höyük (see Figure 3).Since no archive has been discovered in Cilicia so far, only little is known about the ruling dynasties and the political organization for the Late Bronze and Early Iron Age.Today, the Çukurova is an important agriculture center with high productivity.A dozen archaeological surveys give some evidence about the distribution of settlements in the Çukurova and the Yukarıova.The settlement pattern has not changed fundamentally over the time.Mainly during the Bronze Age and in the Roman Period there is a consolidation of settlement visible.A lack of evidence exist in the foothills of the Taurus Mountains north of Adana as well in an area between the rivers Tarsus and Seyhan (Figure 4).
The knowledge of the material culture of Kizzuwatna and Geosciences 2017, 7, 109 Kingdom Kizzuwatna independent, mid 2nd millennium BC, there is a material evidence, with a strong orientation to the Anatolian plateau visib beginning of the 12th century BC, Kizzuwatna stood under Hittite or Mittan At the transition from the Late Bronze to the Early Iron Age Plain Cilicia be as the kingdom Ḫiyawa/Que [29].During the 8th century BC Plain Cilicia Assyrian Empire [30] (pp.[7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24].The knowledge of the material culture of Kizzuwatna and Ḫiyawa/Q Most information come from excavations of the main Bronze and Iron Age are from West to East: Mersin-Yumuktepe, Tarsus-Gözlükule, Adana-T Tatarlı Höyük and Kinet Höyük (see Figure 3).Since no archive has been discovered in Cilicia so far, only little is dynasties and the political organization for the Late Bronze and Earl Çukurova is an important agriculture center with high productivity.A doze give some evidence about the distribution of settlements in the Çukurova settlement pattern has not changed fundamentally over the time.Mainly du in the Roman Period there is a consolidation of settlement visible.A lack foothills of the Taurus Mountains north of Adana as well in an area betwe Seyhan (Figure 4).
Geosciences 2017, 7, 109 4 of 20 Kingdom Kizzuwatna independent, mid 2nd millennium BC, there is a sudden change in the material evidence, with a strong orientation to the Anatolian plateau visible.From mid-15th to the beginning of the 12th century BC, Kizzuwatna stood under Hittite or Mittani supremacy [28] (p.36).At the transition from the Late Bronze to the Early Iron Age Plain Cilicia became again independent as the kingdom Ḫiyawa/Que [29].During the 8th century BC Plain Cilicia was integrated into the Assyrian Empire [30] (pp.[7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24]. The knowledge of the material culture of Kizzuwatna and Ḫiyawa/Que is still fragmentary.Most information come from excavations of the main Bronze and Iron Age sites in the plain.These are from West to East: Mersin-Yumuktepe, Tarsus-Gözlükule, Adana-Tepebağ, Sirkeli Höyük, Tatarlı Höyük and Kinet Höyük (see Figure 3).Since no archive has been discovered in Cilicia so far, only little is known about the ruling dynasties and the political organization for the Late Bronze and Early Iron Age.Today, the Çukurova is an important agriculture center with high productivity.A dozen archaeological surveys give some evidence about the distribution of settlements in the Çukurova and the Yukarıova.The settlement pattern has not changed fundamentally over the time.Mainly during the Bronze Age and in the Roman Period there is a consolidation of settlement visible.A lack of evidence exist in the foothills of the Taurus Mountains north of Adana as well in an area between the rivers Tarsus and Seyhan (Figure 4).
The historian and geographer Abū l-Fidāʾ (1273-1331) reports [31] (p.465) that the river Seyhan and Ceyhan join after Misis: "Le Syhan, au-dessous d'Adana et de Messyssa, se joint au Djyhan; tous deux ne forment plus qu'un seul fleuve, qui se jette dans la mer entre Ayas et Tharse" [32] (p.62f).Since no archive has been discovered in Cilicia so far, only little is known about the ruling dynasties and the political organization for the Late Bronze and Early Iron Age.Today, the Çukurova is an important agriculture center with high productivity.A dozen archaeological surveys give some evidence about the distribution of settlements in the Çukurova and the Yukarıova.The settlement pattern has not changed fundamentally over the time.Mainly during the Bronze Age and in the Roman Period there is a consolidation of settlement visible.A lack of evidence exist in the foothills of the Taurus Mountains north of Adana as well in an area between the rivers Tarsus and Seyhan (Figure 4).
The localization of the important classical site Mallos is not clear.It is known that the site must have been close to the river mouth of the Pyramos, the ancient Ceyhan.Forlanini suggest that the site was known in the Late Bronze Age as the harbor town MLWM.He identifies the site at Domuztepe (Seton-Williams No. 74 [26]), situated east of the river Ceyhan, 20 km northeast of the Kap Karataş [10] (p.12).
The localization of the important classical site Mallos is not clear.It is known that the site must have been close to the river mouth of the Pyramos, the ancient Ceyhan.Forlanini suggest that the site was known in the Late Bronze Age as the harbor town MLWM.He identifies the site at Domuztepe (Seton-Williams No. 74 [26]), situated east of the river Ceyhan, 20 km northeast of the Kap Karataş [10] (p.12).

Settlement Database
To analyze changes in the settlement cluster of Cilicia, a database with the focus on Bronze and Iron Age sites was built.Each site has been documented with its coordinates and its occupation phases (spanning from Neolithic to Ottoman Period).Because less than 3% of the sites have been investigated by excavation it was not possible to provide precise information regarding settlement size during any particular period.Most information for the settlement database was provided by a dozen surveys conducted in the area (Figure 4).

Settlement Database
To analyze changes in the settlement cluster of Cilicia, a database with the focus on Bronze and Iron Age sites was built.Each site has been documented with its coordinates and its occupation phases (spanning from Neolithic to Ottoman Period).Because less than 3% of the sites have been investigated by excavation it was not possible to provide precise information regarding settlement size during any particular period.Most information for the settlement database was provided by a dozen surveys conducted in the area (Figure 4).
Tepeba g Höyük, in the center of the present city of Adana, was first excavated with a small sondage in 1936 by the Museum Adana.Recently, excavations have been resumed first by the Çukurova University and later by the Museum of Adana.It can be identified with hitt.Adaniya and classical Antioch on the Saros [29] (pp.191-195), [34] (p.154).To gain information about routes through the Taurus Helmuth Theodor Bossert and Bahadir Alkım conducted a survey in the Upper Ceyhan valley from 1945-1952 and in 1955 [35].They discovered the late Hittite fortress Karatepe-Arslantaş (hitt.Azatiwataya) and started to excavate there in 1946 [36][37][38].The first and still most comprehensive survey of Plain Cilicia was directed by Marjory Veronica Seton-Williams in 1951 [26].An important site at the eastern border of the Çukurova, Yakapınar-Misis, was excavated from 1955-1959 by Ludwig Budde and Helmuth Theodor Bossert.In 2000, the work was resumed by Giovanni Salmeri and Ana Lucia D'Agata after conducting a survey in the area.The site might be identified with Zunna  [43].
Information from written sources have available.

CORONA Satellite Imagery
A total of 19% of the sites in the settlem anymore and 76% are damaged.However historical satellite imagery.Declassified in the United States from 1958-1972 [44] have b palaeochannels in Plain Cilicia.The CORO before the intensification of the agriculture destroyed due to the erecting of dams [45 systems KH-4A and KH-4B were analyzed ( of 1.8 m [47] (p.119).For images of the missi of the Middle East were used (Center Arkansas/U.S. Geological Survey; [48]).To missions (DS1043, DS1101, DS1103, DS110 Scollar was used [49].The georeferencing wa identification of archaeological sites and geo visually. ara mentioned in Hittite sources [10] (p. 6).In classical time its name was Mopsuhestia and for a short time Seleucia ad Pyramum [34].Over 400 sites mainly dating from the Classical to the Byzantine Period were investigated by Friedrich Hild and Hansgerd Hellenkemper and published in volume five of the Tabula Imperii Byzantini [34].Since 1990 Mustafa H. Sayar is conducting an epigraphic survey in the entire region.Marie-Henriette Gates in 1991 and later Ann Killebrew and Gunnar Lehmann (2006Lehmann ( -2009) ) were conducting a survey (Mopsos Survey) along the eastern coast of the Gulf of Iskenderun [39,40].From 1992-2011 excavations at the harbor site Kinet Höyük under the direction of Marie-Henriette Gates were undertaken.The site can be identified with Hittite Izziya and classical Issos [41].Only 15 km upstream of Yakapınar-Misis is the site Sirkeli Höyük located.Excavations started in 1936-1937 by John Garstang as part of the Neilson Expedition.Later Excavations were directed by Barthel Hrouda 1992-1996 and Horst Ehringhaus (1997).The work was resumed in 2006 by Mirko Novák, Ekin Kozal and Deniz Yaşin Meier [42].There is strong probability to identify Sirkeli Höyük with Kummanni known from Hittite religious texts [10].A neighboring cult city to Kummanni was the site Lawazantiya/akkad.Lusanda.The site of Tatarlı Höyük is a strong candidate for its identification.Excavations have started in 2007 by K. Serdar Girginer [43].
Information from written sources have been added to the database with up to 1000 sites when available.

CORONA Satellite Imagery
A total of 19% of the sites in the settlement database are not visible on modern satellite imagery anymore and 76% are damaged.However, many of these archaeological sites are still visible on historical satellite imagery.Declassified in 1995, CORONA satellite images (Figure 5) acquired by the United States from 1958-1972 [44] have been used in addition to Digital Elevation Models to map palaeochannels in Plain Cilicia.The CORONA images are valuable because they were acquired before the intensification of the agriculture and the growth of urban areas.Some of the sites were destroyed due to the erecting of dams [45,46].Images from different missions with the camera systems KH-4A and KH-4B were analyzed (Figure 5).The images have an average spatial resolution of 1.8 m [47] (p.119).For images of the mission DS1107, orthorectified images of the CORONA Atlas of the Middle East were used (Center for Advanced Spatial Technologies, University of Arkansas/U.S. Geological Survey; [48]).To remove geometric distortion of the images from further missions (DS1043, DS1101, DS1103, DS1105, DS1109) the CORONA software written by Irwin Scollar was used [49].The georeferencing was done with ArcGIS v. 10.2 (3rd Order Polynomial).The identification of archaeological sites and geomorphic features on CORONA imagery was performed visually.

TanDEM-X: TerraSAR-X Add-on for Digital Elevation Measurement
The TanDEM-X system is a space borne radar interferometer that is based on two satellites (TerraSAR-X 3) orbiting in close formation [50].The high resolution X-band SARs (9.65 GHz) can be operated in different imaging modes (High Resolution Spotlight, Stripmap, ScanSAR).For this study, two TanDEM-X acquisitions were made possible in two different imaging modes at horizontal polarization.In Stripmap (SM) mode the course of the Ceyhan River was covered at an azimuth resolution of 2.54 m.The High Resolution Spotlight (HS) mode enabled the registration of a smaller area of interest (approx.5 × 5 km 2 ) at a nominal resolution of 1.13 m in azimuth (see [51] for a detailed description of image processing).In addition to TanDEM-X, other freely available DEM products from the shuttle radar topography mission (SRTM C-band version 4.1 with 90 m resolution, SRTM X-band with 30 m resolution) as well as the ASTER Global DEM (V.2) were used to generate a high resolution DEM-mosaic (Figure 6).The global spatial resolution is 5 m.55% of the area are covered by TanDEM-X data divided into 8 spotlight and 22 stripmap scenes.The accuracy of the height value is 4.65 m on average with a standard deviation of 3.82 m [45].The TanDEM-X DEM mosaic provided the basis for the extraction of cross sections of height values (profiles) for the detailed analysis of the topography across riverbeds and paleochannels (Figure 6).

TanDEM-X: TerraSAR-X Add-on for Digital Elevation Measurement
The TanDEM-X system is a space borne radar interferometer that is based on two satellites (TerraSAR-X 3) orbiting in close formation [50].The high resolution X-band SARs (9.65 GHz) can be operated in different imaging modes (High Resolution Spotlight, Stripmap, ScanSAR).For this study, two TanDEM-X acquisitions were made possible in two different imaging modes at horizontal polarization.In Stripmap (SM) mode the course of the Ceyhan River was covered at an azimuth resolution of 2.54 m.The High Resolution Spotlight (HS) mode enabled the registration of a smaller area of interest (approx.5 × 5 km 2 ) at a nominal resolution of 1.13 m in azimuth (see [51] for a detailed description of image processing).In addition to TanDEM-X, other freely available DEM products from the shuttle radar topography mission (SRTM C-band version 4.1 with 90 m resolution, SRTM X-band with 30 m resolution) as well as the ASTER Global DEM (V.2) were used to generate a high resolution DEM-mosaic (Figure 6).The global spatial resolution is 5 m.55% of the area are covered by TanDEM-X data divided into 8 spotlight and 22 stripmap scenes.The accuracy of the height value is 4.65 m on average with a standard deviation of 3.82 m [45].The TanDEM-X DEM mosaic provided the basis for the extraction of cross sections of height values (profiles) for the detailed analysis of the topography across riverbeds and paleochannels (Figure 6).In addition, VHR (very high resolution) satellite images in GoogleEarth and ArcGIS were visually analyzed.

Palaeochannels and Relict Canals
Palaeochannels are sections of abandoned channels [52].They become fully or partly infilled by abandoned channel accretion [53] (pp.136).Palaeochannels are most commonly preserved as exposed cross-sections, abandoned surface forms, or buried channels [54].They were mapped based on the analysis of geomorphic features visible in CORONA satellite imagery from the 1960s, high resolution DEM-data and freely available VHR satellite images (Figure 7, palaeochannels in blue and green color).
Relict canals are abandoned canals for irrigation.They can be traced by their pattern in modern fields and the remains of ancient levees [55] (p.289) (Figure 7, relict canals in pink color).
Since most of the palaeochannels were found between today's river course of Seyhan and Ceyhan, this area will be discussed here.In addition, VHR (very high resolution) satellite images in GoogleEarth and ArcGIS were visually analyzed.

Palaeochannels and Relict Canals
Palaeochannels are sections of abandoned channels [52].They become fully or partly infilled by abandoned channel accretion [53] (p.136).Palaeochannels are most commonly preserved as exposed cross-sections, abandoned surface forms, or buried channels [54].They were mapped based on the analysis of geomorphic features visible in CORONA satellite imagery from the 1960s, high resolution DEM-data and freely available VHR satellite images (Figure 7, palaeochannels in blue and green color).
Relict canals are abandoned canals for irrigation.They can be traced by their pattern in modern fields and the remains of ancient levees [55] (p.289) (Figure 7, relict canals in pink color).
Since most of the palaeochannels were found between today's river course of Seyhan and Ceyhan, this area will be discussed here.Not all palaeochannels are visible on aerial or satellite imagery.Older channels are buried beneath more recent sediments.However, the mapped palaeochannels show different morphologies.The idealized river course corresponds to a sine curve [56] (p.147) (Figure 8).Not all palaeochannels are visible on aerial or satellite imagery.Older channels are buried beneath more recent sediments.However, the mapped palaeochannels show different morphologies.The idealized river course corresponds to a sine curve [56] (p.147) (Figure 8).Not all palaeochannels are visible on aerial or satellite imagery.Older channels are buried beneath more recent sediments.However, the mapped palaeochannels show different morphologies.The idealized river course corresponds to a sine curve [56] (p.147) (Figure 8).The sinuosity ratio gives an indication how much a river meanders by measuring the length of a channel reach and dividing this by the straight line distance [53] (p.138).Rivers in a flat plain follow a winding path where the meander amplitude (A) is higher than the meander wavelength (λ) (Figure 8).This can be seen exemplary in the Çukurova about one kilometer east and west of the rivers Seyhan and Ceyhan.The blue palaeochannels (Figure 7) show the most recent abandoned channels.The former Seyhan river mouth, active until 1000 AD, is clearly visible with a high sinuosity ratio (Figure 7, No. 1) [1,21].A more recent alteration of the Ceyhan River seems to have occurred at the beginning of the 20th century.From the Middle Ages until 1935, the Ceyhan River mouth diverted eastwards [2] (Figure 7, No. 2), thus its present-day location is situated west of where it had been in medieval times [15] (p.656).
Less meandering structures with a lower sinuosity ratio show the green palaeochannels (Figure 7).They are located west of the river Ceyhan but belong to a former river course of the Seyhan as DEM-data show.A main channel of this former river course can be traced up to 7 km before the coastline northeast of the Akyatan Lagoon (Figure 9).The sinuosity ratio gives an indication how much a river meanders by measuring the length of a channel reach and dividing this by the straight line distance [53] (p.138).Rivers in a flat plain follow a winding path where the meander amplitude (A) is higher than the meander wavelength (λ) (Figure 8).This can be seen exemplary in the Çukurova about one kilometer east and west of the rivers Seyhan and Ceyhan.The blue palaeochannels (Figure 7) show the most recent abandoned channels.The former Seyhan river mouth, active until 1000 AD, is clearly visible with a high sinuosity ratio (Figure 7, No. 1) [1,21].A more recent alteration of the Ceyhan River seems to have occurred at the beginning of the 20th century.From the Middle Ages until 1935, the Ceyhan River mouth diverted eastwards [2] (Figure 7, No. 2), thus its present-day location is situated west of where it had been in medieval times [15] (p.656).
Less meandering structures with a lower sinuosity ratio show the green palaeochannels (Figure 7).They are located west of the river Ceyhan but belong to a former river course of the Seyhan as DEM-data show.A main channel of this former river course can be traced up to 7 km before the coastline northeast of the Akyatan Lagoon (Figure 9).The data of the high resolution TanDEM-X DEM allow the generation of cross sections (Figure 10).The large palaeochannel is visible as a positive height anomaly, which is formed by natural levees.Deposition of fine sand, silt and clay alongside and within the channels results in the aggradation of levees which gradually raise the rivers until they flow several meters above the plain on low ridges.As cross sections of the modern Seyhan and Ceyhan rivers show, their riverbed is located 2-4 m above the plain as well (Figure 10).The data of the high resolution TanDEM-X DEM allow the generation of cross sections (Figure 10).The large palaeochannel is visible as a positive height anomaly, which is formed by natural levees.Deposition of fine sand, silt and clay alongside and within the channels results in the aggradation of levees which gradually raise the rivers until they flow several meters above the plain on low ridges.As cross sections of the modern Seyhan and Ceyhan rivers show, their riverbed is located 2-4 m above the plain as well (Figure 10).
The construction of canals for irrigation results in levees as well.In addition, channel excavation and cleaning operations result in the up-cast of clean-out banks alongside the canal [57] (p.415), [58] (p.339).Compared to palaeochannels, (relict) canals show a less meandering course with a low sinuosity ratio.The course of presumed canals was mapped in pink (Figure 7).The spoil banks are visible on optical imagery as CORONA representing white lines bordered the canal (Figure 11).Relict canals itself can be identified as darkened depressions, filled with moist sediments [59].Water flowing in active canals has a dark signature as well (Figure 11).The construction of canals for irrigation results in levees as well.In addition, channel excavation and cleaning operations result in the up-cast of clean-out banks alongside the canal [57] (p.415), [58] (p.339).Compared to palaeochannels, (relict) canals show a less meandering course with a low sinuosity ratio.The course of presumed canals was mapped in pink (Figure 7).The spoil banks are visible on optical imagery as CORONA representing white lines bordered the canal

Distribution of Archaeological Sites in the Area
Even though for precise dating of palaeochannels and relict canals coring and the analysis of sediment is necessary, the location of archaeological sites give some evidence about the dating of these relict canals and palaeochannels (Figure 12).
First, a relation between the age of an archaeological site and its distance to the modern coastline can be observed.The older a site the farther it is situated from the sea.In this area, the closest chalcolithic sites are in distance of 18-20 km to the shoreline.These are Dervisli (Seton-Williams No. 135 [26]) and Domuz III (Seton-Williams No. 102 [26]).While the closest Bronze Age site Gavur Köy (Seton-Williams No. 119 [26]) is only 13 km from today's coastline located.The known Iron Age sites are located in a similar distance to the shoreline.Sites dating to the Roman Period are located right behind the Akyatan Laggon and at the Kap Karataş (Magarsos [60]).These settlement distributions give an idea about shifts of the coastline.
In the area of the large palaeochannel, there is only one site (Domuz III, Seton-Williams No. 102 [26]) dating to the Chalcolithic period known.However, several sites dating to Hellenistic and Roman period are situated close to the palaeochannels (Figure 12, mapped in green).We assume the danger of flooding during the Chalcolithic Period affected this distribution.The accumulation of sediments and the ideal slope for irrigation were the preconditions for sites in later periods.Sites on natural levee ridges were first occupied after these landforms were already developed and no longer subject to flooding [61].
Several sites dating mainly from Bronze Age to the Roman Period are situated along relict canals (mapped in pink).Therefore, the canals have to be in use since the Bronze Age.

Distribution of Archaeological Sites in the Area
Even though for precise dating of palaeochannels and relict canals coring and the analysis of sediment is necessary, the location of archaeological sites give some evidence about the dating of these relict canals and palaeochannels (Figure 12).
First, a relation between the age of an archaeological site and its distance to the modern coastline can be observed.The older a site the farther it is situated from the sea.In this area, the closest chalcolithic sites are in distance of 18-20 km to the shoreline.These are Dervisli (Seton-Williams No. 135 [26]) and Domuz III (Seton-Williams No. 102 [26]).While the closest Bronze Age site Gavur Köy (Seton-Williams No. 119 [26]) is only 13 km from today's coastline located.The known Iron Age sites are located in a similar distance to the shoreline.Sites dating to the Roman Period are located right behind the Akyatan Laggon and at the Kap Karataş (Magarsos [60]).These settlement distributions give an idea about shifts of the coastline.
In the area of the large palaeochannel, there is only one site (Domuz III, Seton-Williams No. 102 [26]) dating to the Chalcolithic period known.However, several sites dating to Hellenistic and Roman period are situated close to the palaeochannels (Figure 12, mapped in green).We assume the danger of flooding during the Chalcolithic Period affected this distribution.The accumulation of sediments and the ideal slope for irrigation were the preconditions for sites in later periods.Sites on natural levee ridges were first occupied after these landforms were already developed and no longer subject to flooding [61].
Several sites dating mainly from Bronze Age to the Roman Period are situated along relict canals (mapped in pink).Therefore, the canals have to be in use since the Bronze Age.

Dating of Palaeochannels and the Relict Canals
A soil survey of the Çukurova University [62] show that the palaeochannels mapped in Figure 7 correspond to late Holocene river terraces while the surrounding plain was built during the middle Holocene (Figure 13).The late Holocene river terraces spread in the area between today's river course of Seyhan and Ceyhan.
Additionally, the mapping of former river delta soils give some evidence of former location of the coastline (Figure 13).Fluvisols and Vertisols of the delta reach up to 20 km inland.
Based on the data of the soil survey of the Çukurova University (Figure 13) and the presence of Bronze Age sites along relict canals (Figure 12) we propose that the river Seyhan followed already in the Bronze Age (3rd and 2nd millennium BC) a western course like today (Figure 14).The intake of the relict canals must have been located west of the supposed former course of the large palaeochannel (mapped in green, Figure 12).Therefore, the large palaeochannel must date before the Bronze Age, possibly in the Neolithic (9th-7th millennium BC) or Chalcolithic Period (6th-4th millennium BC).However, a shared river mouth in previous times of former Seyhan and Ceyhan rivers-as authors in Antiquity like Pseudo-Scylax or in Medieval times like Abū l-Fidāʾ reported [31] (p.465)-is conceivable for the time when the large palaeochannel was active.
The Late Bronze Age harbor town MLWM could have been a river harbor situated east of the river Ceyhan likely identified with the archaeological site Domuztepe [10] (p.12). Figure 14 shows a preliminary reconstruction of the palaeo-landscape from the mid to the late Holocene.

Dating of Palaeochannels and the Relict Canals
A soil survey of the Çukurova University [62] show that the palaeochannels mapped in Figure 7 correspond to late Holocene river terraces while the surrounding plain was built during the middle Holocene (Figure 13).The late Holocene river terraces spread in the area between today's river course of Seyhan and Ceyhan.
Additionally, the mapping of former river delta soils give some evidence of former location of the coastline (Figure 13).Fluvisols and Vertisols of the delta reach up to 20 km inland.
Based on the data of the soil survey of the Çukurova University (Figure 13) and the presence of Bronze Age sites along relict canals (Figure 12) we propose that the river Seyhan followed already in the Bronze Age (3rd and 2nd millennium BC) a western course like today (Figure 14).The intake of the relict canals must have been located west of the supposed former course of the large palaeochannel (mapped in green, Figure 12).Therefore, the large palaeochannel must date before the Bronze Age, possibly in the Neolithic (9th-7th millennium BC) or Chalcolithic Period (6th-4th millennium BC).However, a shared river mouth in previous times of former Seyhan and Ceyhan rivers-as authors in Antiquity like Pseudo-Scylax or in Medieval times like Ab ū l-Fidā' reported [31] (p.465)-is conceivable for the time when the large palaeochannel was active.
The Late Bronze Age harbor town MLWM could have been a river harbor situated east of the river Ceyhan likely identified with the archaeological site Domuztepe [10] (p.12). Figure 14 shows a preliminary reconstruction of the palaeo-landscape from the mid to the late Holocene.Today, several modern sites are located along a modern canal following the western relict canal (Figure 15, No. 1).However, there are no modern sites located along the eastern relict canal (Figure 15, No. 2).Therefore, the eastern relict canal must have been in use during previous periods.Since several sites dating to the Bronze Age are located along these relict canals it is possible that they were in use from the Bronze Age to Roman Times.Today, several modern sites are located along a modern canal following the western relict canal (Figure 15, No. 1).However, there are no modern sites located along the eastern relict canal (Figure 15, No. 2).Therefore, the eastern relict canal must have been in use during previous periods.Since several sites dating to the Bronze Age are located along these relict canals it is possible that they were in use from the Bronze Age to Roman Times.

Irrigation System
The average annual rainfall in the Çukurova is over 750 mm/year [63].This rate is high enough for rainfed agriculture.Supplementary irrigation systems increase the productivity of rainfed areas.Due to its fertile alluvial soils and mild climate, the Çukurova produce high crop yields twice a year.This makes the region an important contributor to the gross value of Turkeys agricultural crop production of 5% in 1995 [63].In antiquity the region was famous for agricultural products [34] (p.109-111).
A large limestone statue of a male bearded figure standing on a basalt basis showing a chariot pulled by two bulls was found 1997 in Çineköy (Figure 3), a village in the middle of the Çukurova [64].The male figure represents a god, who must be identified with the Storm-god Tarhunzas.A bilingual inscription in Luwian and Phoenician is engraved in the monument [64][65][66].The more than 2.50 m high monument is shown today in the Museum Adana.The inscription is formed of two main parts: an introduction dedicated to the self-presentation of the king and a longer section describing the king's achievements [66].The king praise himself: "I caused the plain of Hiyawa to prosper" [66]

Irrigation System
The average annual rainfall in the Çukurova is over 750 mm/year [63].This rate is high enough for rainfed agriculture.Supplementary irrigation systems increase the productivity of rainfed areas.Due to its fertile alluvial soils and mild climate, the Çukurova produce high crop yields twice a year.This makes the region an important contributor to the gross value of Turkeys agricultural crop production of 5% in 1995 [63].In antiquity the region was famous for agricultural products [34] (pp.109-111).
A large limestone statue of a male bearded figure standing on a basalt basis showing a chariot pulled by two bulls was found 1997 in Çineköy (Figure 3), a village in the middle of the Çukurova [64].The male figure represents a god, who must be identified with the Storm-god Tarhunzas.A bilingual inscription in Luwian and Phoenician is engraved in the monument [64][65][66].The more than 2.50 m high monument is shown today in the Museum Adana.The inscription is formed of two main parts: an introduction dedicated to the self-presentation of the king and a longer section describing the king's achievements [66].The king praise himself: "I caused the plain of Hiyawa to prosper" [66] (p.187).He named the plain also as "land of the rivers" [64] (p.972).There is no large Late Bronze Age site in vicinity known which would explain the location of such monument.However, the find spot is located at a relict canal.Therefore, it might by possible that the monument was erected at an artificial canal.Additional irrigation must have used already during the Late Bronze and Iron Age.
The overlay of the modern canal network for irrigation with mapped relict canals and palaeochannels show that the course of modern canals corresponds to relict canals and palaeochannels (Figure 16).Today, most canals in the Çukurova are constructed with concrete slabs.Many canals have a meandering course even though a straight channel network would be more effective and reduce the number of necessary channels.Thus, the ideal slope of the palaeochannels was followed when constructing modern canals.(p.187).He named the plain also as "land of the rivers" [64] (p.972).There is no large Late Bronze Age site in vicinity known which would explain the location of such monument.However, the find spot is located at a relict canal.Therefore, it might by possible that the monument was erected at an artificial canal.Additional irrigation must have used already during the Late Bronze and Iron Age.The overlay of the modern canal network for irrigation with mapped relict canals and palaeochannels show that the course of modern canals corresponds to relict canals and palaeochannels (Figure 16).Today, most canals in the Çukurova are constructed with concrete slabs.Many canals have a meandering course even though a straight channel network would be more effective and reduce the number of necessary channels.Thus, the ideal slope of the palaeochannels was followed when constructing modern canals.

The Advantage of High Resolution DEM Data for Archaeology
Only a few archaeological sites are published with their coordinates.Therefore, we needed additional information for the precise localization of the sites for the settlement database.The CORONA imagery was an important source of information.Still, small sites were not always visible on this historic satellite imagery.For such cases, the TanDEM-X based height model was an important addition.A low elevated hill can be still visible in a height model while it is not visible anymore on optical data.
Large palaeochannels can be detected easily and even smaller palaeochannels can be uncovered with a stretched color ramp.High resolution DEM-data helps not only to detect features.It is also possible to take precise measurements.The data allows generating cross section to make small differences in height of 1 m visible.Nevertheless, it does not replace ground truthing.

The Advantage of High Resolution DEM Data for Archaeology
Only a few archaeological sites are published with their coordinates.Therefore, we needed additional information for the precise localization of the sites for the settlement database.The CORONA imagery was an important source of information.Still, small sites were not always visible on this historic satellite imagery.For such cases, the TanDEM-X based height model was an important addition.A low elevated hill can be still visible in a height model while it is not visible anymore on optical data.
Large palaeochannels can be detected easily and even smaller palaeochannels can be uncovered with a stretched color ramp.High resolution DEM-data helps not only to detect features.It is also possible to take precise measurements.The data allows generating cross section to make small differences in height of 1 m visible.Nevertheless, it does not replace ground truthing.

Conclusions
The combined use of remote sensing data, soil surveys and archaeological data are first indications for the reconstruction of ancient river courses, when there is only limited coring available.Based on the available data we assume a more eastern course of the river Seyhan during the early or middle Holocene/the Neolithic Period.A major avulsion took place during the late Neolithic or Chalcolithic Period.This shift of the river Seyhan to the west must have an influence to the distribution of sites.In comparison to southern Mesopotamia, where major shifts of rivers were leading to depopulation of large areas, the consequences in the Çukurova might have been less drastic.While irrigation was essential for agriculture in southern Mesopotamia farmers in the Çukurova could still rely on rainfed agriculture for their subsistence.The influence of channel avulsion on settlement must be further investigated.
Today, the landscape of the Çukurova is structured from a bird's eye view by the agricultural area as well as modern river beds, irrigation canals and channel avulsion.Despite the meandering character of palaeochannels modern irrigation canals still follow their course.The location of archaeological sites along both relict canals and modern irrigation canals demonstrate that irrigation canals are in use-maybe not continuous-over a very long period.The finding of such a big monument as the Çineköy statue close to a relict canal gives further evidence for the use of irrigation in this area at least since the Late Bronze Age.
However, the limitations of this study can only be overcome by further investigations.The analysis of drilling cores along transects to get information about the past shifts of the coastline and the sedimentation history of the Ceyhan and Seyhan rivers is necessary.

Figure 2 .
Figure 2. Morphological and simplified tectonic map of the region with the locations of major Neogene basins and the main fault zones (FZ) (after [19] Figure 1), © Susanne Rutishauser.

Figure 2 .
Figure 2. Morphological and simplified tectonic map of the region with the locations of major Neogene basins and the main fault zones (FZ) (after [19] Figure 1), © Susanne Rutishauser.

Figure 3 .
Figure 3. Map of Plain Cilicia and its neighboring regions with the most important Bronze and Iron Age sites as well as modern cities.Topographic data based on SRTM [9] © Susanne Rutishauser.

Figure 3 .
Figure 3. Map of Plain Cilicia and its neighboring regions with the most impo Age sites as well as modern cities.Topographic data based on SRTM [9] © Sus

Figure 3 .
Figure 3. Map of Plain Cilicia and its neighboring regions with the most important Bronze and Iron Age sites as well as modern cities.Topographic data based on SRTM [9] © Susanne Rutishauser.

Figure 3 .
Figure 3. Map of Plain Cilicia and its neighboring regions with the most important Bronze and Iron Age sites as well as modern cities.Topographic data based on SRTM [9] © Susanne Rutishauser.

Figure 4 .
Figure 4. Distribution of archaeological sites investigated during different surveys in the Cilician Plain © Susanne Rutishauser.

Figure 4 .
Figure 4. Distribution of archaeological sites investigated during different surveys in the Cilician Plain © Susanne Rutishauser.

Figure 6 .
Figure 6.DEM-Mosaic based on TanDEM-X, Aster GDEM and SRTM data with modern river courses © Susanne Rutishauser.

Figure 6 .
Figure 6.DEM-Mosaic based on TanDEM-X, Aster GDEM and SRTM data with modern river courses © Susanne Rutishauser.

Figure 8 .
Figure 8.The dotted line represents an idealized river course, the continuous line a meandering river course.λ: meander wavelength, A: meander amplitude, r: meander radius, AG: meander belt © Susanne Rutishauser.

Figure 8 .
Figure 8.The dotted line represents an idealized river course, the continuous line a meandering river course.λ: meander wavelength, A: meander amplitude, r: meander radius, AG: meander belt © Susanne Rutishauser.

Figure 8 .
Figure 8.The dotted line represents an idealized river course, the continuous line a meandering river course.λ: meander wavelength, A: meander amplitude, r: meander radius, A G : meander belt © Susanne Rutishauser.

Figure 9 .
Figure 9. High-resolution digital elevation model based on TanDEM-X data (Stripmap) of the Çukurova.The grey arrows point to a large palaeochannel.The location of the cross sections is represented by black lines.The dots mark a distance of 2000 m © Susanne Rutishauser.

Figure 9 .
Figure 9. High-resolution digital elevation model based on TanDEM-X data (Stripmap) of the Çukurova.The grey arrows point to a large palaeochannel.The location of the cross sections is represented by black lines.The dots mark a distance of 2000 m © Susanne Rutishauser.

Geosciences 2017, 7 , 109 11 of 20 Figure 10 .
Figure 10.Cross sections of height values from the TanDEM-X DEM for the rivers Seyhan and Ceyhan and the large palaeochannel.For the location, see Figure 9.The blue rectangles mark the location of the riverbed resp.palaeochannel © Susanne Rutishauser.

Figure 10 .
Figure 10.Cross sections of height values from the TanDEM-X DEM for the rivers Seyhan and Ceyhan and the large palaeochannel.For the location, see Figure 9.The blue rectangles mark the location of the riverbed resp.palaeochannel © Susanne Rutishauser.

Figure 11 .
Figure 11.Detail of a CORONA scene (DS1105-2267, DA003, 20.11.1968).The course of the large palaeochannel visible on the DEM is marked with red arrows.A modern canal is following its course.The blue arrows points to a recent canal from the 1960s.The green arrow point to a relict canal © Susanne Rutishauser.

Figure 11 .
Figure 11.Detail of a CORONA scene (DS1105-2267, DA003, 20.11.1968).The course of the large palaeochannel visible on the DEM is marked with red arrows.A modern canal is following its course.The blue arrows points to a recent canal from the 1960s.The green arrow point to a relict canal © Susanne Rutishauser.

Figure 12 .
Figure 12.Sites dating from the Chalcolithic to the Roman Period in the Çukurova as well as palaeochannels and relict canals © Susanne Rutishauser.

Figure 12 .
Figure 12.Sites dating from the Chalcolithic to the Roman Period in the Çukurova as well as palaeochannels and relict canals © Susanne Rutishauser.

Figure 14 .
Figure 14.Preliminary reconstruction of the palaeo-landscape of the eastern Çukurova.The schematic black lines represent former river courses of Seyhan and Ceyhan, the pink lines show relict canals, © Susanne Rutishauser.

Figure 14 .
Figure 14.Preliminary reconstruction of the palaeo-landscape of the eastern Çukurova.The schematic black lines represent former river courses of Seyhan and Ceyhan, the pink lines show relict canals, © Susanne Rutishauser.

Figure 14 .
Figure 14.Preliminary reconstruction of the palaeo-landscape of the eastern Çukurova.The schematic black lines represent former river courses of Seyhan and Ceyhan, the pink lines show relict canals, © Susanne Rutishauser.

Figure 15 .
Figure 15.Location of archaeological and modern sites in the Çukurova.Modern river courses (grey), palaeochannels (blue, green) and relict canals (pink) in the Çukurova © Susanne Rutishauser.

Figure 15 .
Figure 15.Location of archaeological and modern sites in the Çukurova.Modern river courses (grey), palaeochannels (blue, green) and relict canals (pink) in the Çukurova © Susanne Rutishauser.

Figure 16 .
Figure 16.Modern river courses (grey), palaeochannels (blue, green) and relict canals (pink) in the Çukurova.The network of thin grey lines shows the modern network for canals for irrigation (mapped based on freely available VHR data in ArcGIS) © Susanne Rutishauser.

Figure 16 .
Figure 16.Modern river courses (grey), palaeochannels (blue, green) and relict canals (pink) in the Çukurova.The network of thin grey lines shows the modern network for canals for irrigation (mapped based on freely available VHR data in ArcGIS) © Susanne Rutishauser.
Geosciences 2017, 7, 109 direction of John Garstang started to excava in its first year of the survey.The work was site might be identified with Hittite Ellipra, N Tepebağ Höyük, in the center of the p sondage in 1936 by the Museum Adana.Çukurova University and later by the Muse classical Antioch on the Saros [29] (pp.191 through the Taurus Helmuth Theodor Bosse Ceyhan valley from 1945-1952 and in Karatepe-Arslantaş (hitt.Azatiwataya) and st most comprehensive survey of Plain Cilici 1951 [26].An important site at the eastern b from 1955-1959 by Ludwig Budde and Helm Giovanni Salmeri and Ana Lucia D'Agata a identified with Zunnaḫara mentioned in H Mopsuhestia and for a short time Seleucia ad Classical to the Byzantine Period were inve and published in volume five of the Tabula conducting an epigraphic survey in the ent Killebrew and Gunnar Lehmann (2006-2009 eastern coast of the Gulf of Iskenderun [39,4 Höyük under the direction of Marie-Henri with Hittite Izziya and classical Issos [41].On Höyük located.Excavations started in 1936-Later Excavations were directed by Barthe work was resumed in 2006 by Mirko Novák probability to identify Sirkeli Höyük with neighboring cult city to Kummanni was t Höyük is a strong candidate for its identif Girginer