An Overview of the Geomorphological Characteristics of the Pergamon Micro-Region (Bakırçay and Madra River Catchments, Aegean Region, West Turkey)

: Current landforms are the basis for understanding past geomorphodynamics and human activities. Based on multiple materials and methods, including geomorphometric analyses of a digital elevation model and visual interpretations of satellite images, different landscape units of the Bakırçay and Madra River catchments in the environs of ancient Pergamon are described. The area was mainly shaped by tectonics that formed a horst-and-graben structure; small Miocene horsts and NE-trending subgrabens in the Kozak and Yunt Da˘gı Mountains were separated by the NEE-trending Plio-Pleistocene Bergama Graben. The asymmetry in relief, drainage network, and sediment accumulation between the Kozak Mountains to the north and the Yunt Da˘gı Mountains to the south of the lower Bakırçay plain characterize the Pergamon Micro-Region. The regional relief characteristics, with wide ﬂat basin and plateau areas, are suitable for agriculture. Complemented by its richness in natural resources, the Pergamon Micro-Region became a preferred settlement area, with evidence of human impact since the Hellenistic–Roman period at the latest. As a consequence of settlement activities, several landscape engineering measures were implemented simultaneously with a parallel change in morphodynamics.

The current study focuses on one of these valleys, the Bakırçay valley, which is directly related to the Pergamon Micro-Region.This microregion has been settled at least since the Neolithic period [26].The site of Pergamon, on the hilltop northeast of modern Bergama (Figure 1A), had already become a polis by the fourth century BC and then became the capital of the Kingdom of Pergamon in 281-133 BC.Pergamon doubled its size prior to the beginning of the third century and subsequently the decline had already set in [27,28].The Hellenistic-Roman urban center of Pergamon (a UNESCO World Heritage site) and its micro-region have been intensively studied, predominantly from the archaeological point of view, since the end of the 19th century [27,[29][30][31][32][33][34].More recently, geoarchaeological research has been integrated, focusing on selected archaeological sites and related subbasins along the western lower Bakırçay plain [8][9][10][11][35][36][37][38][39][40][41][42] and the neighboring Madra River delta [43].(2) aqueduct: The Atlas Project of Roman Aqueducts (https://www.romaq.org/the-project/aqueducts/article/692#tab-details,last accessed 8 June 2021) [20] (p.150); (3) mine site [21,22]; (4) elevation data: TanDEM-X digital elevation model (12-m horizontal resolution) [23,24].(B) The location of the study area (red rectangle), highlighted with the major westward-trending rivers (blue lines) and their alluvial plains (green polygons) along the Aegean coastal region of west Turkey (modified from [1]), and the major tectonic structures (shown as black thin lines in the west of Anatolia and dashed lines linked to the sutures across the Aegean Sea) and tectonic zones (black thick line: the ˙Izmir-Ankara Suture Zone; BFZ: Bornova Flysch Zone) modified from [25]).(C) Geomorphological units and subunits of the Pergamon Micro-Region.
The current study focuses on one of these valleys, the Bakırçay valley, which is directly related to the Pergamon Micro-Region.This microregion has been settled at least since the Neolithic period [26].The site of Pergamon, on the hilltop northeast of modern Bergama (Figure 1A), had already become a polis by the fourth century BC and then became the capital of the Kingdom of Pergamon in 281-133 BC.Pergamon doubled its size prior to the beginning of the third century and subsequently the decline had already set in [27,28].The Hellenistic-Roman urban center of Pergamon (a UNESCO World Heritage site) and its micro-region have been intensively studied, predominantly from the archaeological point of view, since the end of the 19th century [27,[29][30][31][32][33][34].More recently, geoarchaeological research has been integrated, focusing on selected archaeological sites and related subbasins Land 2021, 10, 667 3 of 27 along the western lower Bakırçay plain [8][9][10][11][35][36][37][38][39][40][41][42] and the neighboring Madra River delta [43].
The majority of prehistoric and ancient settlements in the Pergamon Micro-Region were located in the transitional areas between the mountains and alluvial plains or along flat unconfined valleys (Figure 1A) [9][10][11][17][18][19].These areas appear to be particularly suitable for the foundation of human settlements and for local supplies and route networks.They are relatively flood-protected, adjacent to fertile plains for agriculture and to mountainous hinterland for horticulture and pastoralism [9, 17,18].Examples include the Greek polis Atarneus (from the sixth to the fourth century BC) in the southwestern Kozak Mountains [8,9] and the Late Chalcolithic/Byzantine settlement Güm4 in the northern Yunt Da gı Mountains (Figure 1A) [18,19].However, some of the original flood-safe areas became prone to flooding because of continuous aggradation of the lower Bakırçay plain and occasionally flood episodes, e.g., the site Teut114 after the 1st century CE [10] and the Late Chalcolithic-to-Bronze Age settlement Yeni Yelde girmentepe [11] in the western lower Bakırçay plain.In addition to the settlements adjacent to the main Bakırçay alluvial plain, prehistoric sites were located on or close to the coast during the Middle Holocene, when alluvial and colluvial conditions began to dominate the littoral areas [1].The settlements on the Madra River delta and around the Bakırçay River delta are prominent examples [36,38,43].Flourishing in a sheltered position, the Greco-Roman harbor city of Elaia (Figure 1A) avoided the presumably high sedimentation of the Bakırçay River [36,38].Several ancient human settlements and archaeological remains have been documented in the mountains adjoining the lower Bakırçay plain [17,31,44], which, however, lack systematic geoarchaeological research.
Previous geomorphological and sedimentological studies indicate not only the environmental potential for human activities, but also the Holocene local morphodynamics in the western lower Bakırçay plain and the adjacent piedmonts.Moreover, a meta-analysis of the 14 C-ages from site-related sediment sequences reconstructed several phases of increased geomorphodynamics during the Holocene, which are presumably related to multiple triggers, such as climate, population dynamics, and local settlement history [42].
The morphological complexity of a landscape in a catchment subsequently influences the sediment connectivity [52].Landscape configuration controls where and how long the sediments are stored [53,54].Wider valleys have larger accommodation spaces for sediments to be stored within the catchment and less in the basin outlet [54], which could promote landscape stability.Combined factors of neotectonics [48], cycles of climate [55] and sea level changes [1,39], vegetation degradation, and increased human activities [40,42] are able to increase sediment and landscape dynamics [53,54] and consequently shape the local geomorphology [56,57].
An identification of geomorphological characteristics forms the basis for a reconstruction of geomorphodynamics, which allows insights to be gained into the distribution of human activities in a landscape, including ancient human settlements and route networks.Despite the previously mentioned site-related geoarchaeological, regional, and supra-regional geological and paleoenvironmental studies, geomorphological characterizations of the Bakırçay and Madra River catchments in the environs of ancient Pergamon are currently lacking.Therefore, our main objectives are (a) to provide a catchment-scale overview and (sub)unit comparison of current landforms in the Bakırçay and Madra River catchments based on geomorphological mapping and geomorphometric analyses; and (b) to assess potential drivers of landscape dynamics and to infer the consequences for different landscape units.On this basis, the attractiveness and appropriateness of the Pergamon Micro-Region as a preferred settlement area since Hellenistic times are sketched.

Geographical Setting
The study area comprises the Bakırçay River catchment (ca.3382 km 2 ) and the neighboring Madra River catchment (ca.382 km 2 ), and the adjacent coastal areas along the Aegean Sea (Figure 1A).It includes a major part of the Pergamon Micro-Region [28], where important natural resources (such as stone as a building material, clay for ceramic production, soil for agricultural production, and water as a daily necessity) have been utilized to supply the city of Pergamon [17,43,58].Located between 38 • 54 -39 • 29 N and 26 • 48 -27 • 57 E, this area extends approximately 100 km in the E-W direction and 66 km in the N-S direction, covering an area of ca.4124 km 2 .
The Bakırçay River (the ancient Kaïkos) is one of the E-W oriented graben-related rivers of the Aegean Region in Turkey [5,[61][62][63] (Figure 1B).It originates from the Kocada g (in the Gölcük Da gları Mountains; Figure 1A).After flowing through the upper and the (eastern and western) lower Bakırçay alluvial plain, it enters into the Aegean Sea and develops a delta between Zeytinda g (west Yunt Da gı Mountains) and Çandarlı (south Kara Da gı Mountains) (Figure 1A).Its natural channel pattern changes from braided in the middle reaches to meandering in its lower reaches [11,61].The Madra River drains through the intramountainous basin of the Kozak interior plain and it incises the steep valley in the western Kozak Mountains (Figure 1A).After passing through the Altınova-Dikili littoral plain, deposits of the Madra River form a delta in the Aegean Sea [43], which is smaller than the Bakırçay delta.

Tectonics, Bedrock, and Soil
Geologically, the study area belongs to the westward-drifting Aegean-Anatolian microplate [25,64], comprised by the Sakarya Zone in the north and the Bornova Flysch Zone in the south (Figure 1B) [25].Among these continental zones is the NE-trending ˙Izmir-Ankara Suture Zone (IASZ) (Figure 1B) [25,48].NE-trending dextral strike slip faults of the North Anatolian Fault System and N-NW-trending extension of the Aegean Extensional System simultaneously controlled the tectonic structures in this area [48].Several horsts and grabens formed under the tectonic paleostress, e.g., the Kozak Horst, the Miocene Zeytinda g and Örenli-E giller grabens, and the Plio-Pleistocene Bergama and Ayvalik-Lesvos grabens (Figure 2) [48,65,66].The Bergama Graben continues to subside at an average rate less than 1 cm per 1000 years [39,61].The tectonic stress caused severe earthquakes during the last two millennia, for example, the Çandarlı earthquake in 105 CE [67] (p.54), the Dikili earthquake in 1939, and the Göçbeyli earthquake in 1919, the magnitudes of the latter two being ≥6.5 [68] (Figure 2).Within its western sector, several Miocene volcanic domes (Turkish tepe, hills) rise from the Bakırçay plain [50].Faults are mainly found in three zones, viz., the Altınova-Dikili fault zone (Figure 2: 1 ), the Zeytinda g-Bergama fault zone (Figure 2: 2 ), and the Soma-Kırka gaç fault zone (Figure 2: 3 ) [69][70][71][72].The latter two zones are dominated by normal faults that were active during the Holocene (Figure 2).Pre-Miocene granodioritic rocks form the core of the Kozak Mountains and are surrounded by Mesozoic metamorphic rocks, Miocene carbonate rocks, and volcanic formations [60] (Figure 2).The volcanic units are dominated by andesites, basalts, and pyroclastic rocks.The latter rock can also be found at the western Yunt Dağı Mountains, e.g., the Maruflar Horst (Figure 2).Miocene continental clastic rocks were deposited at the outer belt of the Kozak Horst, mostly covering the old carbonate rocks in the Gölcük Dağları Mountains and some at the depressions in the Yunt Dağı Mountains [60] (Figure 2).
Pre-Miocene granodioritic rocks form the core of the Kozak Mountains and are surrounded by Mesozoic metamorphic rocks, Miocene carbonate rocks, and volcanic formations [60] (Figure 2).The volcanic units are dominated by andesites, basalts, and pyroclastic rocks.The latter rock can also be found at the western Yunt Da gı Mountains, e.g., the Maruflar Horst (Figure 2).Miocene continental clastic rocks were deposited at the outer belt of the Kozak Horst, mostly covering the old carbonate rocks in the Gölcük Da gları Mountains and some at the depressions in the Yunt Da gı Mountains [60] (Figure 2).
The formation of different soils in this area mainly depends on lithology and vegetation cover [74].The predominant soil types are Cambisols-Leptosols in the mountains and Fluvisols in the alluvial plains [75,76].Beyond these, brown forest soils, which dominantly formed on granodiorite parent materials, are the most common soil type in the Kozak Mountains; Rendzinas are mainly distributed on calcareous-rich sedimentary bedrocks [74].In the area along Soma-Kırka gaç, Luvisols prevail [75,76].

Climate
Located in the Aegean Region, the study area has a typical Mediterranean subhumid climate [77].The area is characterized by dry and hot summers, and mild and wet winters ("Csa" climate according to the Köppen-Geiger classification [78]).Between 1981 and 2010, the annual temperature averaged 16.2 • C in Dikili and 14.7 • C in Bergama, and the annual precipitation averaged 711 mm in Dikili and 749 mm in Bergama (Figure 3).At both locations, most rain falls in December and January, reaching 140 mm per month; the dry period spans from May to October (Figure 3).The maximum precipitation recorded at the meteorological station in Dikili between 1939 and 1996 totaled 183 mm within 24 h [74].

Climate
Located in the Aegean Region, the study area has a typical Mediterranean subhumid climate [77].The area is characterized by dry and hot summers, and mild and wet winters ("Csa" climate according to the Köppen-Geiger classification [78]).Between 1981 and 2010, the annual temperature averaged 16.2 °C in Dikili and 14.7 °C in Bergama, and the annual precipitation averaged 711 mm in Dikili and 749 mm in Bergama (Figure 3).At both locations, most rain falls in December and January, reaching 140 mm per month; the dry period spans from May to October (Figure 3).The maximum precipitation recorded at the meteorological station in Dikili between 1939 and 1996 totaled 183 mm within 24 h [74].http://data.ceda.ac.uk/badc/cru/data/cru_ts/cru_ts_4.01/data, last accessed 8 November 2020) [79].
Between 9 and 3 ka BP, climate conditions in the eastern Mediterranean became gradually drier after the relatively wet early Holocene; the period between 6 and 3 ka BP appears to have been the driest period during the Holocene [55].More rapid (centennialscale) climate change events were interspersed around 6.7-6.3 ka BP, 4.5-4.3ka BP, and 3.5-3.3ka BP in the Mediterranean [55].During the past three millennia, the hydroclimate became slightly wetter [57].Climate changes, particularly rapid climate change events, Between 9 and 3 ka BP, climate conditions in the eastern Mediterranean became gradually drier after the relatively wet early Holocene; the period between 6 and 3 ka BP appears to have been the driest period during the Holocene [55].More rapid (centennial-scale) climate change events were interspersed around 6.7-6.3 ka BP, 4.5-4.3ka BP, and 3.5-3.3ka BP in the Mediterranean [55].During the past three millennia, the hydroclimate became slightly wetter [57].Climate changes, particularly rapid climate change events, are assumed to be in chronological agreement with major archaeological turnover events in the Aegean/Levantine region [80].

Vegetation
During the early Holocene (9500-8500 BP), the landscape ecology in the Mediterranean region remained stable, with the forest extent reaching a maximum [57].The gradual decline in forest between 7500 and 3000 BP was impacted both by climate change and human activities; subsequently, the landscape has continued to be subject to anthropogenic pressure [57].Now, the remnants of the natural vegetation cover have retreated to steep valley sites (Figure S6, Corine Land Cover data [81]), appearing predominantly as contracted vegetation [74].
Present-day vegetation cover in the study area is Thermo-Mediterranean [77] or Eu-Mediterranean-type [40].Forest communities and shrub formations (woody shrubs with evergreen leaves: maquis, and xerophyte bushes and thorns: garrigue) of the Mediterranean region are widespread in the mountainous areas [3].Dense forest (mainly coniferous forest) is dominant at elevations above 500 m a.s.l. at the outer belt of the Kozak intramountainous basin and the eastern Yunt Da gı Mountains (Figure S6).Stone pine forests (Pinus pinea) occur in the Kozak Mountains in elevations between 600 and 700 m a.s.l.[74].Shrub and herbaceous vegetation associations are widespread in the Kara Da gı Mountains at 200 and 500 m a.s.l.The western slopes of the Kozak Mountains and the Maruflar Horst (in Yunt Da gı Mountains) are intensively used for grazing [9,43] (Figure S6).The Zeytinda g Graben (the western Yunt Da gı Mountains) and the margins of the upper Bakırçay plain (Soma-Kırka gaç area) are largely cultivated with olive trees (Figure S6).The common vegetation of the Bakırçay and Madra alluvial plains comprises cash crops, mainly cotton and vegetables [3,74].Sporadically, crops are also cultivated in intramountainous basins and infilled valleys (Figure S6).

Survey of Geomorphological Features
Major geomorphological features were mapped using QGIS 3.8 [82], supported by a literature review, field observations, visual interpretations of freely available highresolution satellite images provided by Google Earth (<1 m resolution; DigitalGlobe, Inc., maps.google.com,last accessed 8 June 2020), and geomorphometric derivates from TanDEM-X data [23,24] (Section 3.2.).Their subdivisions with definition and classification rules are summarized in Table S1 [83].The materials used for interpretation and field photos of each landform unit are provided in Supplementary Materials Figures S1-S7.
Rivers were digitized from topographic and geological maps [69][70][71][72].Please note that the field observations were mainly conducted in the lower Bakırçay plain and its surrounding mountains, whereas the validation of the results from other areas is mostly based on the literature (e.g., [45,46,49]).

Extraction of Basic Geomorphometric Parameters
Based on a digital elevation model derived from TanDEM-X data (12-m horizontal resolution; please note that the elevations received from the TanDEM-X are ellipsoidal heights; the undulation of the geoid is around 40 m at the coast [23,24]), selected geomorphometric analyses [84] were conducted using various modules implemented in GRASS GIS [85] and SAGA GIS [86] in QGIS 3.8 [82].Computed geomorphometric derivates were shaded relief, slope, geomorphons, aspect, topographic position index and topographic wetness index (Table 1).Zonal statistical analyses of geomorphological subunits were processed using the Zonal statistics module in QGIS 3.8 [82].
Table 1.Geomorphometric derivates (based on TanDEM-X data [23,24]) used for the identification of landforms in this study.

Shaded relief
Hillshade Shaded relief representation of the terrain surface.Slope r.slope.aspect[87] The slope angle of a raster pixel.

Geomorphons r.geomorphon * [88]
A recognition approach for the classification and mapping of the ten most common landform elements.Aspect r.slope.aspect[87] The slope orientation of a raster pixel.
Topographic Position Index topographic position index (TPI) [89,90] Comparison of each raster pixel's elevation to the mean elevation of a specified neighborhood around that pixel.Topographic Wetness Index r.topidx [91][92][93] The capacity of a raster pixel to accumulate water.
* Regarding the parameters of r.geomorphon in this study, the outer search radius was set to be 1000 m, the inner search radius 50 m, and the flatness threshold 3 degrees.Geomorphon classes in the main text are shown in italics.Settings for other derivates remained at default values.
To obtain information on the major channel courses, river longitudinal profiles [94] (pp.294-312) were extracted for several major river courses from the TanDEM-X data [69][70][71][72].Noise in the river longitudinal profiles caused by reservoirs was manually excluded.Knickpoints were identified manually as inflexion points of the river longitudinal profiles.
To assess the shape of a mountain front, the mountain-front sinuosity index (Smf = Lmf /Ls; where Lmf is the planimetric length of a mountain front along the footslope which is at the topographic break in slope, and Ls is the length of a straight line at the former mountain front) was applied and calculated in QGIS 3.8 [82].This represents the straight mountain front produced by active tectonics when the Smf value is low, and sinuous or eroded front due to erosive processes and less active tectonics when the Smf value is high [95,96].

General Geomorphological Characteristics
The relief of the Bakırçay and Madra River catchments and their adjacent coastal areas can be grouped into three major landform types: (1) upland landforms including ridge, slope with debris cover and with bare rock, and valley with or without fluvial infill; (2) lowland landforms; and (3) anthropogenic landforms (Figure 4).
In the mountainous areas, upland landforms prevail, of which ridges and slopes with debris cover dominate.Outcropping bare rocks occur locally.Most of the perennial and intermittent river valleys are filled with fluvial deposits, whereas their tributary valleys in general have distinctly steeper inclining stream beds than the receiving streams and lack fluvial infills.Furthermore, anthropogenic landforms, including reservoirs, mines, and agricultural terraces, are widespread in the mountainous area.Alluvial landforms refer to the plains aggraded by the Bakırçay River; paleochannels in the lower Bakırçay plain document the anastomosing character of the Bakırçay River before straightening and embanking.To drain the lower Bakırçay plain, a wide network of canals was installed.The coastal areas are characterized by littoral plains and the deltas of Bakırçay and Madra Rivers.In the transition to the plain areas, the footslopes are fringed by colluvial deposits, locally interfingering with fan deposits from mountain rivers.S1 [83].Note that paleochannels and terraces are not exhaustive.Geomorphological (sub)units are indicated by numbers; for territorial extent of the geomorphological (sub)units see Figure 1C.
In the mountainous areas, upland landforms prevail, of which ridges and slopes with debris cover dominate.Outcropping bare rocks occur locally.Most of the perennial and intermittent river valleys are filled with fluvial deposits, whereas their tributary valleys in general have distinctly steeper inclining stream beds than the receiving streams and lack fluvial infills.Furthermore, anthropogenic landforms, including reservoirs, mines, and agricultural terraces, are widespread in the mountainous area.Alluvial landforms refer to the plains aggraded by the Bakırçay River; paleochannels in the lower Bakırçay plain document the anastomosing character of the Bakırçay River before straightening and embanking.To drain the lower Bakırçay plain, a wide network of canals was installed.The coastal areas are characterized by littoral plains and the deltas of Bakırçay and Madra Rivers.In the transition to the plain areas, the footslopes are fringed by colluvial deposits, locally interfingering with fan deposits from mountain rivers.
The entire area of the Bakırçay and Madra River catchments has a mean elevation of ca.377 m (σ = 233 m) with an elevation range of ca.1432 m (Table 2).Corresponding to the mountainous relief, the slope inclination averages around 12°, with steep slopes regularly exceeding 10° in the mountains and flat-gentle terrain less than 3° in the plain areas (Table 2; Figure 5A).Landform units as displayed by geomorphons are dominated by slopes and spurs in the mountainous areas whereas the geomorphon flat prevails in plain areas (Figure 5B).All of these geomorphometric parameters clearly show spatial diversities in different (sub)units (Figures 5 and S7).S1 [83].Note that paleochannels and terraces are not exhaustive.Geomorphological (sub)units are indicated by numbers; for territorial extent of the geomorphological (sub)units see Figure 1C.
The entire area of the Bakırçay and Madra River catchments has a mean elevation of ca.377 m (σ = 233 m) with an elevation range of ca.1432 m (Table 2).Corresponding to the mountainous relief, the slope inclination averages around 12 • , with steep slopes regularly exceeding 10 • in the mountains and flat-gentle terrain less than 3 • in the plain areas (Table 2; Figure 5A).Landform units as displayed by geomorphons are dominated by slopes and spurs in the mountainous areas whereas the geomorphon flat prevails in plain areas (Figure 5B).All of these geomorphometric parameters clearly show spatial diversities in different (sub)units (Figure 5 and Figure S7).
Table 2. Geomorphological parameters in the four geomorphological units and their subunits of the study area (see Figure 1C for the (sub)units).

Geomorphological (Sub)Units
Area (km  Our results show extensive asymmetrical features of the Bakırçay and Madra River catchments, as evident in their geomorphometric characteristics.In particular, the asymmetry between the Northern Highlands and the Yunt Da gı Mountains surrounding the (western and eastern) lower Bakırçay plain is pronounced (Figure 6).The swath profile characterizes the asymmetry between the mountains surrounding the western lower Bakırçay plain (Figure 6: profile A-A'), viz., the Kozak Mountains to the north, and the western Yunt Da gı Mountains to the south.The eastern lower Bakırçay plain, with the Northern Highlands adjoining to the north and the eastern Yunt Da gı Mountains adjoining to the south (Figure 6: profile B-B'), is also asymmetrical, although with less-distinct geomorphological features.* Elevation: ellipsoidal height, based on TanDEM-X data [23,24].
Our results show extensive asymmetrical features of the Bakırçay and Madra River catchments, as evident in their geomorphometric characteristics.In particular, the asymmetry between the Northern Highlands and the Yunt Dağı Mountains surrounding the (western and eastern) lower Bakırçay plain is pronounced (Figure 6).The swath profile characterizes the asymmetry between the mountains surrounding the western lower Bakırçay plain (Figure 6: profile A-A'), viz., the Kozak Mountains to the north, and the western Yunt Dağı Mountains to the south.The eastern lower Bakırçay plain, with the Northern Highlands adjoining to the north and the eastern Yunt Dağı Mountains adjoining to the south (Figure 6: profile B-B'), is also asymmetrical, although with less-distinct geomorphological features.

Major Geomorphological Units
Corresponding to the major lithological (Figure 2) and geomorphological characteristics (Figure 4), the study area was divided into the four major geomorphological units-1-Kara Dağı Mountains, 2-Northern Highlands, 3-Yunt Dağı Mountains, and 4-Plain areas-with a total of thirteen subunits (Figure 1C).

Major Geomorphological Units
Corresponding to the major lithological (Figure 2) and geomorphological characteristics (Figure 4), the study area was divided into the four major geomorphological units-1-Kara Da gı Mountains, 2-Northern Highlands, 3-Yunt Da gı Mountains, and 4-Plain areas-with a total of thirteen subunits (Figure 1C).

Geomorphological Unit 1-Kara Da gı Mountains
Geomorphological unit 1 corresponds to the Kara Da gı Mountains and covers an area of 198 km 2 (~5% of the study area; Table 2); this unit is not further subdivided.The Kara Da gı Mountains are dome-shaped, with the highest elevation at its central southwestern part (ca.860 m a.s.l.) and circumferentially gently decreasing slopes (Figure 1A).In its eastern part, slopes generally incline less than 16 • , whereas in its western part slopes reach inclinations up to 30 • (Figure 5A).In its eastern part, narrow elongated valleys less than 1 km wide cut through a rolling landscape of slopes and ridges; extended slope areas occur at the western part of the Kara Da gı Mountains at its dip to the Aegean (Figure 5B).Agricultural terraces occur locally along the elongated valley flanks; these valleys are mostly infilled with shallow fluvial and colluvial deposits (Figure 4).Bare rocks are dominantly exposed at the southeastern upslope areas, whereas the downslope connecting footslope is covered by colluvial deposits, interfingering with littoral deposits of the Bay of Elaia (Figure 4).

Geomorphological Unit 2-Northern Highlands
Geomorphological unit 2 corresponds to the Northern Highlands, which include the major Kozak Mountains (subunits 2.1 and 2.2), Korucu Da gları Mountains (subunit 2.3), Golcük Da gları Mountains (subunit 2.4), and their intramountainous basins (Figure 1A).This is the largest geomorphological unit in the study area, covering more than half of its area (Table 2).Subunits 2.1-2.3 were separated based on the major drainage divides, whereas subunit 2.4 was delineated according to differences in parent rocks with sed-imentary rocks characterizing subunit 2.4, and mainly volcanic and pyroclastic rocks characterizing subunit 2.3.The Madra River (subunit 2.1) directly flows through the Altınova-Dikili littoral plain into the Aegean Sea.
Geomorphological subunits 2.1 and 2.2 correspond to the Kozak Mountains, with the Madra River catchment draining the northern part of the Kozak mountains delineated as subunit 2.1.The Kozak Mountains have a mean elevation of around 530 m a.s.l. in subunit 2.1 and 440 m a.s.l. in subunit 2.2, with a total elevation range c. 1400 m (Table 2).Slope inclinations average 18 • (Table 2) and locally reach a maximum of 80 • (Figure 5A).Within the Kozak Mountains, the Kozak interior plain at 500-600 m a.s.l.altitude is located as a NE-striking intramountainous basin.The basin floor is filled by sediments, with the fluvial deposits of the Madra River forming overlying strata (Figure 4) [51].The relief of the Kozak interior plain is shallow (<3 • ), corresponding to its high topographic wetness index (TWI) values (Figure S7C).The predominating geomorphons are ridges, spurs, slopes, hollows, and valleys in the mountainous areas beyond the Kozak interior plain (Figures 5B  and 7).The high heterogeneity of the relief is also displayed by small scale changes in the values of the topographic position index (TPI) (Figure S7B).
Headwater areas of the Madra, Geyikli, Bergama, and Kestel Rivers rise from the Kozak mountains (Figure 4).Due to year-round spring discharge, these rivers sourcing from the Kozak Mountains are perennial.The river courses of the Madra River and the Kestel River show clear rectangular drainage patterns [43,51] (Figure 5A: a).The Madra River is dammed midstream at the transition of the Kozak Mountains into the coastal plain (Figure 4).Valley slopes are only rarely terraced, with agricultural terraces predominantly occurring at the Kozak Mountains' southern slopes (subunit 2.2) (Figure 4).Tributaries sourcing from the Kozak Mountains are short and intermittent, with valley infills predominantly found along their lower courses (Figure 4).Longitudinal profiles of the river courses (subunit 2.2) have a distinct concavity with steep upstream channel gradients that rapidly decrease at the mountain's foot zone and achieve a stretched characteristic in the downstream part (Figure 8), whereas flowing through the Kozak Mountains, repeatedly occurring knickpoints mark riffles and rapids (Figure 8).
Along the southwestern ridges of the Kozak Mountains, extended areas with outcropping bedrock occur (Figure 4).These locations frequently coincide with small extraction areas and small stone quarries, spread particularly across the northern slopes of the Kozak interior plain (subunit 2.1).Due to the plutonic rocks forming the northern and western part of the Kozak mountains, outcropping bedrock corresponds to granite boulders and their weathering products (Figure S4A).The Ovacik Gold Mine [21,98] is located at the southern footslopes of the Kozak Mountains at the transition to the western lower Bakırçay plain (Figure 2).Geomorphological subunits 2.3 and 2.4 exhibit a NE-striking depression-mountaindepression pattern, including the major morphotectonic landforms (from west to east): the Örenli-E giller Graben, Korucu Da gları Mountains, Savaştepe interior plain, Gölcük Da gları Mountains, and Gelenbe interior plain (Figures 1A and 2).The Örenli-E giller Graben (100-300 m a.s.l.) is characterized by a slightly rolling landscape with an average inclination of <6 • and an average TWI of seven; the graben floor is covered by fluvial and colluvial deposits with locally occurring landslides.Graben flanks, and the Kozak Mountains to the northwest and the Korucu Da gları Mountains to the east, have a strong relief, with c. 15 • average inclination (subunit 2.3; Figure 5A).The basin structures of the Savaştepe interior plain and the Gelenbe interior plain are located at an elevation of c. 250 m a.s.l.; like the Kozak interior plain, they are filled with fluvial and colluvial deposits, resulting in a shallow relief with an average TWI of 7.1 (Figure S7C).Outcropping bare rocks are widely scattered across the mountainous areas of subunits 2.3 and 2.4 (Figure 4).inantly occurring at the Kozak Mountains' southern slopes (subunit 2.2) (Figure 4).Tributaries sourcing from the Kozak Mountains are short and intermittent, with valley infills predominantly found along their lower courses (Figure 4).Longitudinal profiles of the river courses (subunit 2.2) have a distinct concavity with steep upstream channel gradients that rapidly decrease at the mountain's foot zone and achieve a stretched characteristic in the downstream part (Figure 8), whereas flowing through the Kozak Mountains, repeatedly occurring knickpoints mark riffles and rapids (Figure 8).Along the southwestern ridges of the Kozak Mountains, extended areas with outcropping bedrock occur (Figure 4).These locations frequently coincide with small extraction areas and small stone quarries, spread particularly across the northern slopes of the Kozak interior plain (subunit 2.1).Due to the plutonic rocks forming the northern and western part of the Kozak mountains, outcropping bedrock corresponds to granite boulders and their weathering products (Figure S4A).The Ovacik Gold Mine [21,98] is located at the southern footslopes of the Kozak Mountains at the transition to the western lower Bakırçay plain (Figure 2).Geomorphological subunits 2.3 and 2.4 exhibit a NE-striking depression-mountaindepression pattern, including the major morphotectonic landforms (from west to east): the Örenli-Eğiller Graben, Korucu Dağları Mountains, Savaştepe interior plain, Gölcük Dağları Mountains, and Gelenbe interior plain (Figures 1A and 2).The Örenli-Eğiller Graben (100-300 m a.s.l.) is characterized by a slightly rolling landscape with an average inclination of <6° and an average TWI of seven; the graben floor is covered by fluvial and colluvial deposits with locally occurring landslides.Graben flanks, and the Kozak Mountains to the northwest and the Korucu Dağları Mountains to the east, have a strong relief, with c. 15° average inclination (subunit 2.3; Figure 5A).The basin structures of the Savaştepe interior plain and the Gelenbe interior plain are located at an elevation of c. 250 m a.s.l.; like the Kozak interior plain, they are filled with fluvial and colluvial deposits, resulting in a shallow relief with an average TWI of 7.1 (Figure S7C).Outcropping bare Tributaries of Bakırçay River with sources in the Korucu Da gları Mountains (subunit 2.3) and the Gölcük Da gları Mountains (subunit 2.4) show concave longitudinal channel profiles, which are typical for mountainous tributaries with very steep upper courses repeatedly interrupted by knickpoints (Figure 8).A clear relationship between the degree of concavity and the ratio of the channel length to the altitude of the headwater area can be observed, indicating stronger morphodynamics for the Deliktas and Sitma Rivers than for the Ilica and Ya gcılı Rivers.
Reservoirs have been built at the middle-lower reaches of the major streams (Deliktaş, Sıtma, and Ya gcılı) (Figure 4), some of which collect runoff from drainage basin areas up to 450 km 2 in size, e.g., the Sevişler reservoir on the Ya gcılı River between the Korucu Da gları and Gölcük Da gları Mountains [99] (Figure 4).A large mining area is located next to the Soma-Kolin Thermal Power Plant in the Korucu Da gları Mountains (subunit 2.4; Figures 2 and 4) covering an area of c. 15 km 2 .Agricultural terraces mainly occur on the Gölcük Da gları Mountains (subunit 2.3; Figure 4).

Geomorphological Unit 3-Yunt Da gı Mountains
Geomorphological unit 3 corresponds to the northern part of the Yunt Da gı Mountains and incorporates about one quarter of the research area (Table 2).Including the Maruflar Horst, the ridges along the divide reach elevations >1200 m a.s.l.(Table 2) and are widely composed of pyroclastic rocks (88% of its area).The Zeytinda g and Yayaköy grabens (subunit 3.1; Figure 2) confine the Maruflar Horst to the NW and N, both being infilled by Quaternary fluvial deposits (Figure 4).In the NE of Maruflar Horst, De girmendere Graben and Kınık Horst (subunit 3.3) display the basin-range structure, building the transition into Bergama Graben.The Kınık Horst corresponds to a block of uplifted Miocene carbonates, whereas adjacent eastern Yunt Da gı Mountains (subunit 3.4) are mainly built out of continental clastic rocks (Figures 2 and 7).
Geomorphological subunit 3.1.The ground level of the Zeytinda g and Yayaköy grabens (subunit 3.1) is at 180-190 m a.s.l.altitude; Quaternary graben infills cause a moderate relief with inclinations <10 • (Table 2 and Figure 5A).Streams are generally short and intermittent, especially those close to Çalıbahçe and Zeytinda g (Figure 4).Additionally, agricultural terraces and mining areas are widely spread throughout the area (Figure 4).Geomorphological subunit 3.2.The Maruflar Horst (subunit 3.2) has a mean elevation of c. 516 m and forms the highest area of the Yunt Da gı Mountains.The horst is in the headwater areas characterized by small flats and footslopes (areas of less than 2 km 2 ), and wide and elongated ridges and slopes being up to c. 3 km wide (Figure 5B).In the headwater areas, valleys are saucer-shaped and infilled by shallow fluvial deposits.Valleys are deeply incised and still under active incision at the confluence into the Zeytinda g and Yayaköy grabens.At the lower courses of the Sınır and Gümüş Rivers, in particular, the valley floors widen and are infilled by fan deposits; fans gently incline and cover areas of c. 2 km 2 (Sınır) and 5 km 2 (Gümüş) (Figure 4).Agricultural terraces occur predominantly on the uphill graben shoulders of Zeytinda g and Yayaköy grabens.
Channel courses in the western Yunt Da gı Mountains (subunits 3.1 and 3.2) show distinct convex shapes in their longitudinal profiles in the upstream areas, tracing the transition of saucer-shaped valleys in plateau-like environments of the headwater areas to the graben flanks (Figure 8).Repeatedly occurring knickpoints mark fault lines crossed by the river courses (see Gümüş, Tekkedere, and Sınır Rivers in Figure 8).
Geomorphological subunit 3.3.The basin-range structure of Kınık Horst and De girmendere Graben (subunit 3.3) covers an elevation range of 1152 m; the 15.5 • average slope inclination documents the strong relief (Table 2).Landform units are characterized by short and narrow ridges, slopes, and valleys, whereas there is a lack of flats (Figure 7); overall, the landscape of Kınık Horst and De girmendere Graben is highly heterogeneous (Figure 5B).Corresponding to the strong relief, valleys are deeply incised, with fluvial deposits only occurring in basin-like valley openings, such as in the lower course of Yayaköy River (Figure 4).The size of most of these basins is less than 1 km 2 .Slopes are widely covered by debris and only locally bare rocks crop out (Figure 4); correspondingly, agricultural terraces are widely distributed.
Geomorphological subunit 3.4.The eastern Yunt Da gı Mountains (subunit 3.4) are characterized by steep slopes of more than 30 • inclination, particularly at the areas adjoining to the Soma valley and upper Bakırçay plain (Figure 5A).The slopes are widely covered by debris, whereas bare rocks only locally crop out.Soma valley, with its 1-2-km-wide valley floor, is also part of subunit 3.4, expressed by a distinct occurrence of flats (Figure 5B).Kısrakdere Coal Mine, with a quarrying area of more than 40 km 2 , is located in the central eastern Yunt Da gı Mountains (Figure S4E).Agricultural terraces occur only rarely (Figure 4).
Longitudinal profiles of channel courses in the eastern Yunt Da gı Mountains (subunit 3.3 and part of subunit 3.4) have a concave shape, locally interrupted by convex knickpoints, e.g., Kara and Işıkören Rivers (Figure 8).In particular, the longitudinal channel profile of the Işıkören River is marked by strong inclinations at the upper and middle courses, and a distinct knickpoint in the middle course.

Geomorphological Unit 4-Plains
The plain areas (geomorphological unit 4) are widely flat, with inclinations of less than 3 • (Figure 5A).Fluvial, colluvial, and littoral Quaternary deposits cover the plain areas (Figure 4).The headwater areas of Bakırçay River are located in the Gölcük Da gları Mountains from where the river course heads southward until reaching the basin of the upper Bakırçay plain (subunit 4.3).In the upper Bakırçay plain, the river flows westward before turning north when reaching the eastern spur of the Yunt Da gı Mountains.In particular, the Bakırçay River close to Soma flows along a meandering valley, which is predetermined by the adjoining mountain ranges and the local fault systems.Downstream of Soma, the valley opens to the lower Bakırçay plain.Bakırçay River flows in a WSW direction in the eastern lower Bakırçay plain (subunit 4.3) and turns to the SW after passing the Bergama fan (western lower Bakırçay River plain; subunit 4.2).Throughout the lower Bakırçay plain, the Bakırçay River flows almost along the footslopes of the Yunt Da gı Mountains.
At the transition of the western lower Bakırçay plain into the Aegean Sea, the delta of the Bakırçay River, with an area of 14 km 2 , is developed between the modern settlement of Çandarlı (subunit 1) and the Bay of Elaia, which is separated from the delta by the Bozyertepe ridge, a spur of the Yunt Da gı Mountains (subunit 3.1) [38].
Geomorphological subunit 4.1.The Altınova-Dikili littoral plain adjoins the footslopes of the western Kozak Mountains.Correspondingly, the area is composed of a 1-2-km-wide zone covered by coastal deposits interfingering with colluvial deposits along the western footslopes of the Kozak Mountains (Figure 4).Fluvial deposits from the Madra River aggrade and gradually form the Madra River delta downstream from the modern settlement of Altınova [43].The littoral plain also connects to the western lower Bakırçay plain (subunit 4.2) through a narrow corridor close to Dikili.
Geomorphological subunit 4.2.At the northern transition zone between the Kozak Mountains and the western lower Bakırçay plain, colluvial deposits and alluvial fans cover the foot zone of the Kozak Mountains, reaching a width of 2 km (Table 3).In contrast, at the southern boundary of the western lower Bakırçay plain, the colluvium along the footslope of the Yunt Da gı Mountains has a maximum width of 500 m (Figure 4).In addition, the mountain-front sinuosity index (Smf) of the Kozak Mountains and the Yunt Da gı Mountains differs, with a distinctly lower value for the Yunt Da gı Mountains (2.16) than for the Kozak Mountain front (2.73) (Table 3).In the northern part of the Bergama Graben, small dome-shaped volcanic hills occur (Figure 2), e.g., the Kalarga Tepe and E grigöl Tepe.The Bergama fan separates the western lower Bakırçay plain from the eastern lower Bakırçay plain.The city of Bergama was founded in the root-zone of the fan and provides its name.The Bergama fan covers an area of approximately 14 km 2 and receives sediments from the Kozak Mountains by means of the Bergama River and Kestel River, which have a combined drainage basin area of 156 km 2 .Other fans deposited in the western lower Bakırçay plain mostly cover areas of 1-2 km 2 and are deposited by Bakırçay tributaries with drainage basin areas <120 km 2 .In the western lower Bakırçay plain, widespread river regulation measures occur, including canals, stream straightening (e.g., the Do gramacı canal of the Geyikli River), and the construction of river embankments up to 5 m high [10].Abandoned channels and cut-off meanders occur predominantly close to the straightened and embanked river courses (Figure 4).The extensive drainage of wide areas of the lower Bakırçay plain influences local runoff regimes and water budgets.
Geomorphological subunit 4.3.In the eastern lower Bakırçay plain (ca.88 m; Table 2), colluvial deposits at the eastern transition zone from Göçbeyli to Kınık are up to 3 km wide, whereas in the western transition areas they are generally less than 0.5 km wide (Table 3).Alluvial fans with areas of 2-5 km 2 primarily occur in the east of the zone (Figure 4).The sinuosity value of the eastern Yunt Da gı Mountains' front (3.00) is slightly higher than that of the Northern Highlands (2.68) (Table 3).
Geomorphological subunit 4.4.The mean elevation of the upper Bakırçay plain (subunit 4.4) averages ca.213 m (Table 2).The slopes to the Akhisar plain of the Gediz River catchment are gentle; thus, they can be easily connected [17,30].At the footslopes along the upper Bakırçay plain, abundant colluvial sediments were deposited (Figure 4), building a two-kilometer-wide transition zone at the east (from Soma to Kırka gaç) and decreasing in extent to the west.The alluvial fans of the Bakırçay River tributaries have areas of less than 1 km 2 .Only a few canals are found in the upper Bakırçay plain.

Discussion
The distinct asymmetry of the topography in the Bakırçay and Madra River catchments (Figure 6) is an expression of the strong impact of tectonic movements.Locally varying rates of uplift and sinking due to tectonic compression and divergence resulted in a smallchambered relief with sharp borders between the different geomorphological units.These effects are amplified by the inclination of the rock strata forming the mountainous areas due to their locally highly variable characteristics of strike and dip.This diverse landscape promoted the area's suitability for usage since the Holocene.

Development of the Major Landform Units
The Kozak Mountains (geomorphological subunits 2.1 and 2.2) are characterized by high average elevation and steep slopes (Figure 6: A-A'), a NW-trending large-scale bare rock association, and rectangular drainage courses.The mountains comprise Miocene volcanic centers characterized by andesitic volcanism, which surround the plutonic core of the Kozak interior plain [100].In parallel, compressional forces from the westward drift of the Anatolian plate during the Miocene caused an uplift of the area [3,50].During the Miocene, Pliocene, and early Quaternary, the uplift of the Kozak Horst and subsidence of surrounding grabens shaped the Kozak Mountains (Figure 2) [45,51,101].The bare rock association dominating on the western ridges of the Kozak Mountains (subunit 2.1, Figure 4) is associated with the NW-SE-striking Altınova-Dikili fault zone [72] (Figure 2).Steep slopes are deeply fissured, resulting in a high infiltration capacity and shallow soils [74].The rectangular flow pattern of the Madra River is controlled by the local fault system (Figure 5A) and the Mio-Pliocene lava flows [51].In addition to lithological differences between pyroclastic and continental clastic rocks [60] (Figure 2), the rectangular drainage pattern of the Kestel River can also be attributed to the fault system [72] (Figure 5Ac).
A peculiarity of the Kozak Mountains is the wide and flat Kozak interior plain (geomorphological subunit 2.1), a large intramountainous basin filled by Quaternary sediments.Granodiorite boulders in the northern margin of the Kozak interior plain document spheroidal weathering (Figure S4A), building a typical granite landform [102].In addition, the slopes formed in the granodioritic bedrock result in frequently outcropping bare rock covered by a thin layer of weathering products.Comparable interior plains as large sediment sinks can be found at Savaştepe and Gelenbe, in geomorphological subunit 2.4 [103].The formation of the Savaştepe interior plain can be linked to the NE-trending ˙Izmir-Ankara Suture Zone (Figure 2) [25,48].
Compared to the Yunt Dagı Mountains (see below), the colluvial and alluvial fan deposits at the footslopes of the Kozak Mountains are largely extended, being directly related to the size and relief of the drainage basin, which indicate high erosion dynamics [94,[104][105][106].The upper parts of the Geyikli, Bergama, and Kestel Rivers' longitudinal profiles, which each have headwater areas located in the Kozak Mountains (Figure 8), clearly indicate high relief, resulting in strong incision processes [107,108].The knickpoints along the longitudinal profiles are indicative of faults if the bedrock remains unchanged [94,109].
In the western Yunt Da gı Mountains' hinterland, the Maruflar Horst (subunit 3.2) developed in a manner geologically similar to the Kozak Horst [45]; the bedrock of the Maruflar Horst is dominated by pyroclastic rocks formed by Early to Middle Miocene volcanism [45,60].
The Zeytinda g and Yayaköy grabens (subunit 3.1) form local sedimentary sinks.The NE-trending Zeytinda g Graben to the south of the Bakırçay plain-together with Örenli-E giller Graben to the north-developed under the Late Miocene east-west extension between the Kozak Horst and the Göçbeyli Horst [45].The formation of the orthogonal Bergama Graben during the Pliocene separated the Zeytinda g-Örenli-E giller Graben into two units [45], both having the typical topographical characteristics of graben structures with flat infilled graben floors and step faults along the graben shoulders [110] (p.74).The straight mountain fronts of the Yunt Da gı Mountains display regional tectonics, particularly the Holocene normal faults, e.g., the Sindel normal fault [70] (Figure 5Ab).
In comparison to the Kozak Mountains, the gentler relief in the western Yunt Da gı Mountains is expressed in less-pronounced channel gradients [111] (Figure 8).The convexupward shape of the Bakırçay tributaries draining the Yunt Da gı Mountains indicates the plateau-like relief in the headwater area, whereas the repeatedly superimposed concavities in their longitudinal profiles show their high rate of incision with rapid adjustment close to the Bakırçay plain [109].Fluvial sediments can only be found along the river courses in local sinks wherever there is a break in inclination corresponding to the system of sediment cascades [52][53][54].Overall, the Bakırçay tributaries flowing from the Yunt Da gı Mountains show an immature stage of channel development [94], with frequently occurring knickpoints, indicating Quaternary faults [69,72] (Figure 2).Although perennial streams such as the Sınır and Gümüş Rivers merge into the Bakırçay plain, the alluvial fans are only very slightly convex-shaped, due to the shallow inclination of the tributaries' river courses and the corresponding shallow relief of their drainage basins [18].
In the north of the eastern Lower Bakırçay Plain, the NE-trending Örenli-E giller Graben (Figure 2) exhibits a similar elevation distribution and gentle slopes to those of the Zeytinda g Graben (subunit 3.1) [45].The basis of the Örenli-E giller Graben forms a sink, preventing most of the sediment eroded from the mountains from reaching the eastern lower Bakırçay plain (Figure 4).The steep channels of the upper Deliktaş and Sıtma Rivers indicate the lasting uplift of the Kozak Horst [94] (p.306).The llıca and Ya gcılı Rivers show less steep longitudinal channel profiles than the Deliktaş and Sıtma Rivers, indicating lower uplift rates in their headwater areas than at the Kozak Horst [94].Knickpoints along llıca and Ya gcılı River course profiles imply fault lines of the ˙Izmir-Ankara Suture Zone [25].Due to the large drainage basins of the llıca, Ya gcılı, Deliktaş, and Sıtma Rivers, extended alluvial fans formed at their confluence into the eastern lower Bakırçay plain (Figure 4).
The Kınık Horst (subunit 3.3) in the south surmounts the eastern lower Bakırçay plain by approximately 500 m and steep slopes characterize its transition to the eastern lower Bakırçay plain (Figure 6: profile B-B').Its sinuous mountain front (Smf = 3; Table 3) indicates low to moderate tectonic activities [112].Only a narrow verge at the footslope of the Kınık Horst is covered by colluvial deposits (Figure 4).The two alluvial fans of the Yayaköy and Kara Rivers occur because their drainage basin sizes are distinctly larger than those deposited at the footslopes of the western Yunt Da gı Mountains.The rivers deeply incised into Kınık Horst; the lack of intermediate sediment storage areas resulted in large fan areas being slightly disproportionate to the size of their drainage basins.
The Kara Da gı Mountains to the west of the lower Bakırçay plain have a Miocene volcanic center [100] and form as the Kane peninsula reaching into the Aegean Sea.The N-S trending valleys and ridges in the Kara Da gı Mountains (Figure 4) correspond to Cenozoic NNW-SSE striking faults [49,69,71,73] (Figure 2).Extensive colluvial deposits in the southeast of the Kara Da gı Mountains in the transition zone to the western lower Bakırçay plain are indicative of high erosion activity, resulting in deeply incised valleys following faults [49].
The lower Bakırçay plain develops in the east-west trending Bergama Graben.Alluvial deposits of the Bakırçay River, in addition to colluvial and alluvial fan deposits from the adjoining mountains, are continuously aggrading the lower Bakırçay plain.Separating the upper and lower Bakırçay plain, the Soma valley, with a long and narrow shape, is controlled by a series of high-angle oblique-slip normal faults, active since the Pleistocene, with small lateral components that strike NE-SW, NW-E, N-S, and E-W [46,48,70] (Figure 2).The western extension of N-trending faults [70] between Soma and Kırka gaç forms a linear structure of steep slopes facing west.
The formation of the Bergama Graben started from the Early Pliocene [45,48].During the Quaternary, two NNE-to NE-trending dextral oblique-slip fault segments at Zeytinda g-Bergama and Kınık-Turgutalp deformed the former Bergama Graben into a transtensional pull-apart basin striking NE-SW with sediments up to 120 m thick [48].Due to the influence of the former fault zone, the Zeytinda g Graben and the surrounding ridges (subunit 3.1) were moved to the NE (towards the city of Bergama), causing the constriction between the eastern and western lower Bakırçay plain [45,48].Here, the straight channel pattern of the Bakırçay River (Figure 5A: b) displays the tectonic influence on the river course.Presumably due to the continuous uplift of the Kozak Mountains, the Bakırçay River flows close to the footslopes of the Yunt Da gı Mountains, contributing to the asymmetry of the plain [45].
An exceptional landform in the Plain Areas is the Bergama alluvial fan, which is formed from deposits of the Kestel and Bergama Rivers in the Northern Highlands.Mountainous relief [105] and seasonally high discharges allow high transport rates for the bedload deposited at the longitudinal gradient change at the transition from the Kozak Mountains to the Bakırçay Plain [105].
Widespread paleochannels in the lower Bakırçay plain indicate the formerly anastomosing channel pattern of the Bakırçay River, with strong deposition and redeposition dynamics during seasonal flood flows at least since ancient history [10,11].During the Early-Middle Holocene, the Bakırçay floodplain was assumed to be narrower than today [11].
The most prominent coastal landforms in the study area are the Altınova-Dikili littoral plain and the deltas of the Madra and Bakırçay Rivers.The littoral plain is part of the NW-trending Ayvalik-Lesvos Graben, which subsided during the Quaternary [101] (p.19).Similar to the southern Kozak Mountains, the tectonic movement caused steep slopes, which today merge into the coastal plain by colluvium deposited at the footslope [1,51].The Altınova-Dikili coastal area is dominated by Pliocene-Pleistocene alluvial fans based on torrential deposits that extend into the Aegean Sea [2,101].These coarse-detrital fan deposits were flooded by the Holocene marine transgression after 7000-6000 BP [101] and were subsequently overlain by alluvial delta deposits of the Madra River and the Karakoç River (the latter is located c. 5 km north to the Madra River) [39].
The corridor between the Altınova-Dikili littoral plain and the Bakırçay plain can be linked to the active tectonic activities at the Altınova-Dikili fault zone (Figure 2: 1 ) [113].Tectonic activities along the fault zone are also responsible for the thermal springs occurring in this area, e.g., the Kaynarca thermal spring that emerged after the 1939 earthquake [113].In addition to these effects of tectonics, the development of coastal areas was highly impacted by regional sea-level changes [1,114,115] and caused the deposition of interbedded coastal swamp-sandy marine sediments and terrestrial alluvial-deltaic deposits [38,43].

Usage Suitability and Human Impact during the Late Holocene
Prehistoric and ancient settlements in the Pergamon Micro-Region are dominantly found in the flat-gently-inclining transition areas between the mountains and plains (Figure 1A) [12,18].Those areas are easily accessible to diverse resources (including water, food, and building materials) and have geographical strategic significances (such as being flood-safe and providing through-traffic).Fertile plain areas provide the largest portion of the food supply of the micro-region [28].Throughout the Bakırçay plain, settlements prefer flood-safe places, such as the piedmonts along the footslopes of the adjoining mountain areas [10].In the western lower Bakırçay plain, some archaeological sites also occur close to the elevated volcanic cones, e.g., the Yeni Yelde girmentepe site (Figure 1A) [11].The slightly rolling and flood-safe piedmonts at the footslopes of the mountains are currently the preferred locations for pathway and road networks [17].Low-angle colluvial deposits attracted colonization and the foundation of cities, e.g., the modern cities of Soma (ancient Germe), Kırka gaç, and Kınık (Figure 1A) [17].However, the flood-safe areas are not always "safe".Fluvial aggradation during Holocene transgression may occasionally flood the ancient human occupations [10], agricultural fields, and roads, which may hinder the continuous usage of plain areas.
The strong relief of the Kozak Mountains is widely unsuitable for settlements, and only the gently inclined colluvial deposits adjoining the Bakırçay plain are favored for agriculture [8,9].An exception is the hill-top settlement of Kale A gılı (ancient polis Atarneus) in the lower Geyikli drainage basin (Figure 1A), which benefited from its prominent topographical location and advantages such as a strategic position in terms of visual control, easy access to timber from a mountainous hinterland, and fertile agricultural land in the environs [9].The foot-zone of the Kozak Mountains is a preferred area for the throughtraffic between the coast and Bergama, i.e., between the ancient polis of Atarneus and Pergamon (Figure 1A) [17].In addition, the flat relief and fertile fluvial sediments on the Kozak interior plain provide ideal conditions for agriculture (Figure S6).During antiquity, human settlements were found in the area, e.g., Perperene and Kytonion (Figure 1A) [17].Similarly, the low inclinations and infilled unconsolidated sediments make the Örenli-E giller Graben suitable for agriculture (Figure S6) and settlements (Figure 1A) [17].The Kozak Mountains have been incised by large rivers (the Madra and Bergama Rivers) that are integral parts of the road network, with connections among the Kozak interior plain (further to Adramytteion), the coastal areas, and ancient Pergamon (Figure 1A) already existing in ancient times [17].In addition, the llıca, Ya gcılı, Deliktaş, and Sıtma Rivers in the Northern Highlands formed deeply incised valleys that provided natural pathways connecting Pergamon with ancient settlements to the northeast and east, such as Ergasteria and Hadrianoutherai (Figure 1A) [17].
The Bergama fan is a terrestrial landform of high morphodynamics, characterized by frequent sheet-flow and rapid stream-channel shifts, which hampered settlement activities.However, after the incision of the Bergama River into its own deposits, the channel course stabilized, and the flooding risk was reduced.Consequently, natural advantages, including access to water (surface water and groundwater) for drinking and farming, fertile soils for food, and spaces for living, have attracted humans to settle in this area [105,116].Some scholars regard the major reasons for the founding of ancient Pergamon adjacent to the Bergama fan were the fan's physiographic endowment, good accessibility of natural resources, and good strategic location [30,117,118].In spite of the possible accessibility of water from the Bergama alluvial fan, gigantic aqueducts in ancient Pergamon were constructed for continuous water supply as a basic feature of civilized life during Hellenistic-Roman times, such as the aqueduct from the Kozak Mountains to the Acropolis of Pergamon (Figure 1A) built by Eumenes II (197-159 BC).This construction benefited not only from political, economic, and technical developments [20,58], but also from the geographical advantage of having the water source area in the topographically high Kozak Mountains.
For the Pergamon Micro-Region, most recent archaeological surveys in the gentle areas of the Yunt Dagı Mountains assume a large potential for human occupation [12,18,44,110].Olive cultivation in the western Zeytinda g Graben has been reported since at least ~2.8 ka BP [40].The local sediment sinks in the Yunt Da gı Mountains, including the fertile intramountainous basins and flat basin-like extended valleys at the river outlets, favor the development of small-scale settlements and agriculture, e.g., in the Tekkedere, Çalıbahçe, and Gümüş valleys (Figure 1A), in different archaeological periods at least since the Late Chalcolithic period [12,18,32,44,110].In the Yunt Da gı Mountains, small flats with infilled saucer-shaped valleys served for the founding of small hamlets.The soils formed on the fluvial deposits of the saucer-shaped valleys attracted farming; the hamlets were also able to be used as rest areas for travelers on the route from Pergamon to southern cities such as Apollonis and Magnesia (Figure 1A) [17], proving the importance of the areas adjoining the plain for traffic purposes.However, due to the active tectonics, such as earthquakes and the resulting landslides, the settlements in the west Yunt Da gı Mountains may be abandoned and relocated occasionally [110].
The Kara Da gı Mountains make the Bakırçay valley unique along the Anatolian Aegean coast because they restrict the connection between the Bakırçay plain and the Aegean Sea to two narrow corridors [119] (p.15).Settlement areas in the Kara Da gı Mountains predominantly occur along the coastal zone and close to small rivers [17,27,44].The fertile delta plains of Bakırçay and Madra Rivers favor agricultural activities [43], and their marshy ground also makes them widely unsuitable as settlement areas and for through-traffic.The site of the Hellenistic-Roman harbor city of Elaia [35,38] was chosen due to its protected location in the lee of the Kara Da gı Mountains.However, the shallow area of the Bay of Elaia was already silted during Roman times due to the high sediment load of the Bakırçay River and the adjacent slopes [38].
Currently, human activities in the Bakırçay and Madra River catchments and their adjacent coastal areas are not only affected by the characteristics of the natural environments, but also directly shape landforms and locally dominate the landscape.
Hydrologic constructions for water storage are common in the Northern Highlands.These constructions have affected discharge and resulted in reduced sediment transport and alluviation within the lower Bakırçay plain and coastal areas since the 20th century [1,8,10].Mechanized, industrial-scale agriculture is conducted throughout the entire eastern and western lower Bakırçay plain due to (a) river regulation and straightening, (b) drainage by a network of ditches, and (c) flood protection resulting from the construction of dams along the course of the Bakırçay River [18].
Rich mineral resources such as coal [46,120,121] and gold [21,122] are widely extracted [98,123,124], resulting in the formation of mining areas of varying sizes throughout the study area.Coal mining areas are located in the eastern part of the area (Figure 2), i.e., the Soma coal basin, which is one of the largest economic lignite-bearing alluvial basins in Turkey, with lignite and kerogen as the main products [46,124].A gold and silver mine is located at Ovacık-Narlica in the south of the Kozak Mountains (Figure 1A), where the precious metals formed during the Early Miocene [21,98].Mining activities in this area started around 600 BCE under the rule of King Croesus and were further exploited during the Roman period (c.400 CE) [21].The modern Ovacik Gold Mine has been active since ca.1990.Stone or sand quarries are widely distributed throughout the study area and are mostly small in size.The stone quarry at Asa gı Çuma (Figure 1A) has been used at least since the Hellenistic period [22,123].Mining of perlite, a weathering product emerging from the hydration of obsidian, is conducted in the western Yunt Da gı Mountains [18].
Agricultural terraces implemented for slope stabilization and water harvesting are widespread on the slopes of the intramountainous basins and along the slopes of the Bakırçay and Madra tributaries; today they are preferably used for olive groves [44,110] (Figure S6).Ancient agricultural terraces from the Aegean Region date to at least Greco-Roman times [125,126].

Conclusions
The Bakırçay and Madra River catchments and adjacent coastal areas form major parts of the Pergamon Micro-Region, corresponding to the environs of the UNESCO World Heritage site of Pergamon (modern Bergama).The landform characteristic of this area is highly influenced by tectonic processes and structural components.Overall, the Pergamon Micro-Region is an impressive example of a complex morphotectonic landscape, attributed to interfingering factors of tectonics, lithology, sea-level change, and human impact on erosion and sedimentation processes.Cenozoic tectonics and volcanism formed the basic horst and graben structures.Plio-Pleistocene tensional forces and crustal stretching formed the Bergama Graben and the adjoining horsts with their intramountainous basins.The Kozak Mountains to the north of the lower Bakırçay plain, and the Yunt Da gı Mountains to its south, show remarkably asymmetrical geomorphological features, displaying these tectonic influences.The continuous uplift of the Kozak Horst resulted locally in a mountainous relief.Extensive interior plains served as local sediment sinks, which vary in size within highland areas.Rich mineral resources, such as coal, gold and silver, and perlite, are present and have been continuously mined since antiquity, forming extensive quarrying areas.Overall, Quaternary tectonics resulted in a widely immature relief, which remains highly dynamic in the present day.
Due to the availability of several natural resources, the Pergamon Micro-Region was a preferred settlement area in early times.The infills of the intramountainous basins and graben structures provide suitable grounds for agriculture.Perennial streams sourced in the adjoining mountainous areas provide year-round water.Widespread andesitic rocks are an appropriate construction material.The forest areas in the adjoining mountains provided sufficient amounts of wood for construction needs and charcoal production.The topography of the Bergama Graben is characterized by two small seaside accesses and the potentially flood-safe footslopes of the adjoining graben flanks, thus providing a high degree of protection for the settlements on the Bakırçay plain.However, as early as during antiquity, unsuitable land use triggered an enhancement in morphodynamics, resulting in extended colluvial deposition and an increase in the major river sediment load, causing an aggradation of the alluvial plains and a progression of the delta areas, and finally the silting up of the Elaia harbor area.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/10.3390/land10070667/s1.Table S1: The definition and classification rules of landform units for the geomorphological map (modified from [83]).

Figure 3 .
Figure 3. Climate diagrams of Dikili and Bergama (reference period: 1981-2010), showing Mediterranean subhumid climate with arid conditions prevailing during summer months and extreme rainfall during winter months.Data: Climatic Research Unit (CRU) (database:

Figure 4 .
Figure 4. Geomorphological map of the Bakırçay and Madra River catchments according to the classification rules in TableS1[83].Note that paleochannels and terraces are not exhaustive.Geomorphological (sub)units are indicated by numbers; for territorial extent of the geomorphological (sub)units see Figure1C.

Figure 4 .
Figure 4. Geomorphological map of the Bakırçay and Madra River catchments according to the classification rules in TableS1[83].Note that paleochannels and terraces are not exhaustive.Geomorphological (sub)units are indicated by numbers; for territorial extent of the geomorphological (sub)units see Figure1C.

Figure 5 .
Figure 5. Morphometric characters of the Bakırçay and Madra River catchments.(A) Slope with perennial streams (streams were based on[69][70][71][72]) and faults (for fault legends see Figure2; the red dashed rectangles numbered with a, b, c and d denote the geomorphological features under the impacts of tectonics, see the Discussion section).(B) Geomorphons.See Figure1Cfor the geomorphological (sub)units.

Figure 6 .
Figure 6.Swath profiles (for location see Figure 1A), indicating distinctly asymmetrical geomorphological characteristics between the Northern Highlands and Yunt Dağı Mountains surrounding the western lower Bakırçay plain (A-A') and the eastern lower Bakırçay plain (B-B').Note: 20× vertical exaggeration.

Figure 6 .
Figure 6.Swath profiles (for location see Figure 1A), indicating distinctly asymmetrical geomorphological characteristics between the Northern Highlands and Yunt Da gı Mountains surrounding the western lower Bakırçay plain (A-A') and the eastern lower Bakırçay plain (B-B').Note: 20× vertical exaggeration.

Land 2021 , 27 Figure 7 .
Figure 7. Distribution of the morphometric and lithologic characteristics of the Bakırçay and Madra River catchments, outlining geomorphological subunit distribution of geomorphons and elevation; y-axis shows the total number of pixels per geomorphon and per altitude range.Pie charts show the distribution of lithology.

Figure 7 .
Figure 7. Distribution of the morphometric and lithologic characteristics of the Bakırçay and Madra River catchments, outlining geomorphological subunit distribution of geomorphons and elevation; y-axis shows the total number of pixels per geomorphon and per altitude range.Pie charts show the distribution of lithology.
Subunit 4.1 corresponds to the littoral area.Subunits 4.2-4.4display the alluvial plains of the Bakırçay River, separated by the Bergama fan and Soma valley into the western lower Bakırçay plain (subunit 4.2), eastern lower Bakırçay plain (subunit 4.3), and upper Bakırçay plain (subunit 4.4).The width of the lower Bakırçay plain in its central part (around the Bergama fan) is 2-3 km, whereas it increases to c. 10 km in the western lower Bakırçay plain (subunit 4.2) and c. 13 km in the eastern lower Bakırçay plain (subunit 4.3).

Figure S1 :
Çalıbahçe basin, indicating how we use the satellite images and geomorphons to generate the geomorphological map (ridge, valley without fluvial deposits, valley with fluvial deposits and slope with debris cover) and the validation in the field.(A) Data: Google Earth satellite images (WorldView-2, last accessed 8 June 2020).(B) Geomorphological map (legend refers to Figure 4, same below).(C) Landscape photograph.

Figure S3 :
Slope with bare rock in Tekkedere valley.(A) Data: Google Earth satellite images (WorldView-2, last accessed 8 June 2020).(B) Landscape photograph.Figure S4: (A) An example of mining fields and decomposed granite boulders at the northwestern slopes of the Kozak interior plain (subunit 2.1).Data: Google Earth satellite images (WorldView-2, last accessed 8 June 2020).Photographs of a canal on the Bergama Fan (B), slope with debris cover in the Tekkedere valley (C), artificial terraces in Çalıbahçe basin (D), Kısrakdere Coal Mine (E).

Figure S5 :
The western lower Bakırçay plain from the Sultantepe to the Yunt Da gı Mountains.(A) The hillshade of TanDEM-X data, indicating the features of alluvial plain and paleochannels which are highlighted in Figure S5: B. (B) Geomorphological map.(C) Landscape photograph.

Figure S6 :
The land cover map of the Bakırçay and Madra River catchments and the adjacent coastal areas.Data: Corine Land Cover (CLC) 2018 (Version 2020_20u1)[81].

Table 3 .
Mountain-front sinuosity index (Smf) and extension of colluvial deposits in the lower Bakırçay plain.