Search

Various mechanisms are proposed to explain landslide recovery time in the time following major earthquakes. However, research on prescribing possible recovery times following an earthquake is still relatively new. This paper provides an insight into factors governing landslide recovery time, which could be considered as a step forward in predictive modeling for landslide recovery time. To accomplish this, we examined 11 earthquake-affected areas based on the characteristics of both landslide events and landslide sites associated with diverse morphologic and climatic conditions. Our analyses indicate that the dominant characteristics of post-seismic landslide mechanisms determine the recovery time. The characteristics can be identified based on: (i) the fraction of area affected by landslides (%), (ii) mean relief and its standard deviation (m), (iii) average daily accumulated precipitation (mm) and (iv) rainfall seasonality index. If there are not enough co-seismic landslide deposits or not enough relief to trigger large deposits on hillslopes, then the recovery processes are mostly controlled by new landslides caused by a strength reduction of hillslope materials. In most of the cases, this brings a relatively quick recovery process in which the majority of post-seismic landslides may happen within a year or even in a month if sufficient intense rainfalls occur soon after the earthquake. If the predisposing factors create large co-seismic landslide deposits, then remobilization of material takes the role of the dominant mechanism and recovery may take years. Overall, our analyses show that the recovery takes relatively longer if a large amount of co-seismic landslide material is deposited within a high-relief mountainous environment where precipitation rate is low and not persistent.

We present the first detailed Quaternary landform map of the Büyük Menderes Graben System, located in western Turkey which is one of the most active extensional domains in the world. The main map was produced with a combination of TanDEM-X (12.5 m resolution), Red Relief Image Map, unmanned aerial vehicle, Google Earth images, and multiple fieldworks. The main map is presented at a scale of 1:160,000 although landforms were mapped at a scale of 1:15,000. The ten Quaternary landforms were defined considering their surface morphology and depositional environment. The alluvial fans, river terraces, and floodplains are the most common landforms. The spatial pattern of the alluvial fans and river terraces showed the variable rates of tectonics and surface processes along the strike of the graben system. The distribution of the meander cut-offs and paleo-shorelines imply that Büyük Menderes Graben System is vulnerable to geohazards like flash-flood sedimentation and flood events.

Badland areas are present in all continents, excluding Antarctica, and play a critical role in establishing local erosion and sedimentation rates. The presence of unconsolidated rocks (e.g., marls, sandstone, mudstone etc.) is a major driver controlling the distribution of badlands, which together with other environmental components, such as climate, tectonics, vegetation, and topography, determine their forms and processes. The mutual interaction of controlling factors in badlands areas provides a basis for developing a holistic approach to clarify their distribution patterns. Turkey's geodynamic evolution has led to the emergence of marine sedimentary rocks, pyroclastics, and continental clastics, especially in line with the uplift of the Anatolian Plateau and volcanism during the last 8 Ma.This study aims to explore the country-scale distribution of badlands and the controlling factors of this badland distribution in Turkey. Remarkably wide badlands landscapes (4494 km2) have been visually inspected using Google Earth ProTM to further digitize and extract geomorphological units by applying high-resolution multispectral images provided by WorldView-4/Maxar Technology and CNES/Airbus. To obtain exact boundaries, we eliminated contiguous flat areas surrounding the identified badlands by using red relief image map (RRIM) mosaics that express surface concavity and convexity combined with topographic slope derived from a digital elevation model of 5-m spatial resolution. Last, to determine the controlling factors of badlands distribution, we have compiled a global data set comprising 1-km resolution layers of mean annual precipitation, temperature and precipitation seasonality, aridity, NDVI, rainfall erosivity factor, elevation, and majority values of regional lithology in sub-catchments units. The enhanced investigation of the complex relationship that expresses the controlling factors of badlands distribution, has been conducted by K-means unsupervised cluster analysis.Our comprehensive regional analyses exploring the distribution and environmental attributes of major Turkish badlands identified five different groups or clusters of badlands that display spatial coherence with climatic and tectonic settings. We argue that Turkey's climatic and topographic transition zones, varying from Mediterranean climate dominated areas to the more arid Central Anatolian Plateau, and tectonically‑induced topographic barriers play a relevant role in discriminating these groups of badlands. Moreover, the Anatolian diversity of sedimentary rocks, which consists of Neogene and Paleogene continental clastics, volcano clastics & pyroclastics, and lacustrine deposits, makes an essential contribution to the identified, extensive badland distribution.This study has been produced benefiting from the 2232 International Fellowship for Outstanding Researchers Program of the Scientific and Technological Research Council of Turkey (TUBITAK) through grant 118C329. The financial support received from TUBITAK does not mean that the content of the publication is approved in a scientific sense by TUBITAK.

Along the western flank of the northern margin (Central Pontides) of the Central Anatolian Plateau, the humidity from the Black Sea is much higher than the central and eastern flanks and creates a complex relationship between surface and tectonic processes by triggering intense mass wasting activity and aggradation within narrow valleys. We identified three incised fill terrace levels and used Optically Stimulated Luminescence (OSL) dating to calculate fluvial sediment ages and cosmogenic 36Cl exposure dating to calculate limestone boulders exposure ages across the terrace surface. Stratigraphical interpretations and OSL ages of the lowest levels revealed that a fluvial fill terrace formed in the main valley at 275.6 ± 12.8 ka and was overlain by a main river-tributary junction alluvial fan that was abandoned at 39.5 ± 3.5 ka. The results collectively show the influence of climate, topography, hillslope processes, and lithology on aggradation-incision patterns of main rivers. Prolonged aggradation can prevent the channel equilibrium required to calculate rock uplift rates while also causing a new base-level and aggradation upstream. This effect can be exacerbated in uplifting mountainous regions with limited depositional areas. Bedrock incision rates based on the fluvial terrace age were between 0.15 and 0.2 mm/a since 39.5 ± 3.5 ka. However, the high aggradation within this segment of the main valley prevented incision of the channel bedrock for long periods, causing a potential underestimation of the rock uplift rate calculation. Our local period of aggradation appears to be related to increased aggradation and decreased bedrock incision rates measured 14 km upstream that were previously assumed to be the result of decreased tectonic uplift rates. This demonstrates the importance of corroborating strath terrace incision rate estimations with ages and incision rates of downstream fill terraces, if present, to check for potential interference with the tectonic signal. This work was supported by the European Commission within the Marie Curie-ITN ALErT Project [grant number FP7-PEOPLE-2013-ITN , number 607996] and the Istanbul Technical University BAP Project [grant number TDK-2017-40776 ]. We would also like to thank the CNRS-ASTER Laboratory in Aix en Provence France for AMS measurements and AcmeLabs, Canada for measurements of major and trace element concentrations. Reviews by the Geomorphology editor, Sarah Boulton (University of Plymouth), Tim Barrows (University of Wollongong), and an anonymous reviewer have greatly helped to improve the content of this work. This work was supported by the European Commission within the Marie Curie-ITN ALErT Project [grant number FP7-PEOPLE-2013-ITN, number 607996] and the Istanbul Technical University BAP Project [grant number TDK-2017-40776]. We would also like to thank the CNRS-ASTER Laboratory in Aix en Provence France for AMS measurements and AcmeLabs, Canada for measurements of major and trace element concentrations. Reviews by the Geomorphology editor, Sarah Boulton (University of Plymouth), Tim Barrows (University of Wollongong), and an anonymous reviewer have greatly helped to improve the content of this work. Publisher's Version

The extensive gypsum karst of Sivas, Turkey is one of the most outstanding examples of bare gypsum karst in the world. It displays a number of remarkable geomorphic features, including: (1) two stepped planation surfaces cut-across folded gypsum developed during an initial phase of slow base level deepening punctuated by periods of stability; (2) unusual deeply entrenched gypsum canyons related to a subsequent phase of rapid fluvial incision and water table lowering; (3) a polygonal karst of superlative quality mainly developed in the upper surface; (4) relict valleys disrupted by sinkholes in the lower erosional surface; (5) a large number of bedrock collapse sinkholes mostly associated with the lower surface; and (6) numerous cover subsidence sinkholes developed in the valley floors. This work analyses the spatial distribution, characteristics and evolution of the sinkholes within the broad Plio-Quaternary geomorphological and paleohydrological evolution of the epigene karst system dominated by autogenic recharge. A cartographic sinkhole inventory has been produced in an area covering 2820 km(2) with morphometric data and including 295 bedrock collapse sinkholes and 302 cover subsidence sinkholes. The different sinkhole types show a general spatial zonation controlled by the hydrogeological functioning of the different sectors: (1) solution sinkholes (polygonal karst) in the upper recharge area; (2) bedrock collapse sinkholes in the lower denudation surface and close to the base level, where well developed caves are inferred; and (3) cover subsidence sinkholes, with high densities probably associated with areas of preferred groundwater discharge. The morphology of the bedrock collapse sinkholes, varying from small cylindrical holes to large and deep tronco-conical depressions with gentle slopes reflect to geomorphic evolution of these sinkholes that reach exceptionally large hectometre-scale diameters. Their evolution, involving substantial enlargement and deepening, is attributed to the solutional removal as solute load of large volumes of gypsum by downward vadose flow. This type of morphological evolution with significant post-collapse solutional denudation differs from that observed in carbonate rocks characterised by lower solubility and erodibility. The analysis of historical imagery reveals that bedrock collapse sinkholes currently have a very low probability of occurrence and that buried cover subsidence sinkholes are used for urban development creating risk situations. (C) 2021 The Author(s). Published by Elsevier B.V.

The main purpose of the study is to determine the general characteristics of the landslide sizes observed in Cretaceous and Eocene aged flysch assemblages at the Western Black Sea region of Turkey by using magnitude and frequency relations. For this purpose, the magnitude and frequency relations were investigated by considering power-law scaling characteristics of the geological formations. The probability distributions were also examined by considering Double Pareto and Inverse Gamma distribution models. According to the power-law relations, the rollover effects were observed at 0.047 and 0.048 km2, and the fractal dimensions of the distributions were obtained as −1.97 and −1.41 for Cretaceous and Eocene flysch assemblages, respectively. Considering the probability distributions, the best-fits were acquired from the models Double Pareto with three parameters estimated by maximum likelihood estimation for Cretaceous flysch and Kernel density estimation for Eocene flysch. When we compared these results with the results of a study carried out in the same flysch but in another sub-catchment, it is concluded that rollover effects and fractal dimensions may not be generalized, that means, the parameters may differ site to site depending on not only spatial resolution but also morphological, climatic, and anthropogenic features of the region in concern, and conversely, landslide size distributions fit Double Pareto distribution models in general.

Various mechanisms are proposed to explain landslide recovery time in the time following major earthquakes. However, research on prescribing possible recovery times following an earthquake is still relatively new. This paper provides an insight into factors governing landslide recovery time, which could be considered as a step forward in predictive modeling for landslide recovery time. To accomplish this, we examined 11 earthquake-affected areas based on the characteristics of both landslide events and landslide sites associated with diverse morphologic and climatic conditions. Our analyses indicate that the dominant characteristics of post-seismic landslide mechanisms determine the recovery time. The characteristics can be identified based on: (i) the fraction of area affected by landslides (%), (ii) mean relief and its standard deviation (m), (iii) average daily accumulated precipitation (mm) and (iv) rainfall seasonality index. If there are not enough co-seismic landslide deposits or not enough relief to trigger large deposits on hillslopes, then the recovery processes are mostly controlled by new landslides caused by a strength reduction of hillslope materials. In most of the cases, this brings a relatively quick recovery process in which the majority of post-seismic landslides may happen within a year or even in a month if sufficient intense rainfalls occur soon after the earthquake. If the predisposing factors create large co-seismic landslide deposits, then remobilization of material takes the role of the dominant mechanism and recovery may take years. Overall, our analyses show that the recovery takes relatively longer if a large amount of co-seismic landslide material is deposited within a high-relief mountainous environment where precipitation rate is low and not persistent.

High Anatolian orographic margins have large variations in terms of topographic relief, precipitation, and uplift rate. These variations lead to the dynamics of mass movements and surface runoff, which are the dominant geomorphological processes in ice-free mountain landscapes. There is growing recognition that large landslides are important agents of landscape evolution, resulting in massive slope failures, which can cause extensive and rapid topographic changes in many active orogenic belts. Unlike the cognatic orogenic plateau margins in the world, there are no studies available on the large landslides and their geomorphic impact at the margins of the Anatolian Plateau. This study presents results from a regional-scale inventory of 1290 large landslides (> 1 km2) that allow the characterization of spatial distribution and landslide-dominated landscapes in the northern margins of the Anatolian Plateau. The majority of large landslides are clustered in three main zones that correspond to the Western, Central, and Eastern Pontides, which is an east-west-trending orogenic belt that represents a coalesced tectonic entity in the northern section of Turkey. Nearly 80% of large landslides have occurred in a terrain with a mean hillslope relief of > 1000 m in those three landslide-dominated landscapes. The results of regional comparisons reveal that in addition to hillslope relief and steepness, lithotectonic differences largely control the abundance of landslides along the northern margins of the Anatolian Plateau. In this respect, the spatial distribution and abundance of large landslides imply a landscape in which lithological and tectonic controls on hillslope erosion are more significant than climate. The study further shows that the parallel or perpendicular position of the landslides with respect to the direction of the drainage network is effective as positive or negative feedback in response to fluvial dissection of the plateau margins. On the other hand, there is certainly a need for more comprehensive radiometric dating studies to understand the contribution of large landslides on the erosional decay rate of Anatolian Plateau margins. Furthermore, the presence of these large landslides and the derived deposits in this dynamic terrain provide a unique opportunity for deciphering the past climatic and seismic events.

In this study, the high-resolution palynological data derived from two cores taken from Mogan Lake in the Golbasi Basin of the Ankara region in central Anatolia (Turkey) are presented. The results provide the palaeobotanical, palaeoclimatological and palaeoenvironmental records of the last ca. 3100 years for the region which is characterised by rich natural resources and a long human occupation history. The major climatic events and the archaeological-historical periods, which are well established in previous works, along with textual sources, are related to vegetation history inferred from the palynological investigation in the study area. Five main vegetation phases are distinguished and interpreted. Our findings demonstrate that the phases of limited woody vegetation cover are chiefly coupled to several dry climate periods, with some exceptions. Steppe vegetation in an open landscape accompanied by more cultural indicators and lesser arboreal elements is inferred during the Near-East Aridification Phase at the LBA-Iron Age transition (c. 1200-850 BCE), Roman Warm Period before 700 CE, Medieval Warm Period at about 1020 CE, the 1500s and twentieth century. Only during the three dry episodes in the 1600s, 1800s and at the start of Modern Warming Period (in the twentieth century) partial arboreal recovery is deduced, and land-cover patterns are related to both natural and human-induced factors. The most extensive pine-dominated forest cover coincides mainly with the three wet climatic episodes, in Galatians times of the Hellenistic Period at c. 2080 BP, the 700s CE of the middle Byzantine era, the Seljuqs Dynasty and the earlier Ottoman Periods (c. 1165-1432/1435 CE). The former episode is represented by remarkable signs of human interference, while the others have marks of sporadically less cultivated landscape. It is concluded that a combination of various factors can be taken into account as the drivers of vegetation changes in the Ankara region, but climate appears to have been the primary agent in shaping the land-cover. (C) 2021 Elsevier Ltd. All rights reserved.

The Sinop Peninsula is located at the northernmost part of the Asia Minor (Anatolia). Its geographic position between the Central Pontide Mountains and the Black Sea together with the presence of young geological units and landforms provides favorable conditions for understanding onshore and offshore geological and geomorphic processes acting along the northern Anatolian coasts. Here, we focus on some landscapes that constitute one of the best examples along the Turkish Black Sea coast. These are inundated fluvial valleys, uplifted isthmus and marine terraces, and paleo- and active dunes.

We present a new, seismologically consistent expression for the total area and volume of populations of earthquake-triggered landslides. This model builds on a set of scaling relationships between key parameters, such as landslide spatial density, seismic ground acceleration, fault length, earthquake source depth, and seismic moment. To assess the model we have assembled and normalized a catalog of landslide inventories for 40 shallow, continental earthquakes. Low landscape steepness causes systematic overprediction of the total area and volume of landslides. When this effect is accounted for, the model predicts the total landslide volume of 63% of 40 cases to within a factor 2 of the volume estimated from observations (R2=0.76). The prediction of total landslide area is also sensitive to the landscape steepness, but less so than the total volume, and it appears to be sensitive to controls on the landslide size-frequency distribution, and possibly the shaking duration. Some outliers are likely associated with exceptionally strong rock mass in the epicentral area, while others may be related to seismic source complexities ignored by the model. However, the close match between prediction and estimate for about two thirds of cases in our database suggests that rock mass strength is similar in many cases and that our simple seismic model is often adequate, despite the variety of lithologies and tectonic settings covered. This makes our expression suitable for integration into landscape evolution models and application to the anticipation or rapid assessment of secondary hazards associated with earthquakes.

Morphological changes, caused by the erosion and deposition processes due to water discharge and sediment flux occur, in the banks along the river channels and in the estuaries. Flow rate is one of the most important factors that can change river morphology. The geometric shapes of the meanders and the river flow parameters are crucial components in the areas where erosion or deposition occurs in the meandering rivers. Extreme precipitation triggers erosion on the slopes, which causes significant morphological changes in large areas during and after the event. The flow and sediment amount observed in a river basin with extreme precipitation increases and exceeds the long-term average value. Hereby, erosion severity can be determined by performing spatial analyses on remotely sensed imagery acquired before and after an extreme precipitation event. Changes of erosion and deposition along the river channels and overspill channels can be examined by comparing multi-temporal Unmanned Aerial Vehicle (UAV) based Digital Surface Model (DSM) data. In this study, morphological changes in the Büyük Menderes River located in the western Turkey, were monitored with pre-flood (June 2018), during flood (January 2019), and post-flood (September 2019) UAV surveys, and the spatial and volumetric changes of eroded/deposited sediment were quantified. For this purpose, the DSAS (Digital Shoreline Analysis System) method and the DEM of Difference (DoD) method were used to determine the changes on the riverbank and to compare the periodic volumetric morphological changes. Hereby, Structure from Motion (SfM) photogrammetry technique was exploited to a low-cost UAV derived imagery to achieve riverbank, areal and volumetric changes following the extreme rainfall events extracted from the time series of Tropical Rainfall Measuring Mission (TRMM) satellite data. The change analyses were performed to figure out the periodic morphodynamic variations and the impact of the flood on the selected meandering structures. In conclusion, although the river water level increased by 0.4–5.9 meters with the flood occurred in January 2019, the sediment deposition areas reformed after the flood event, as the water level decreased. Two-year monitoring revealed that the sinuosity index (SI) values changed during the flood approached the pre-flood values over time. Moreover, it was observed that the amount of the deposited sediments in September 2019 approached that of June 2018.

2018 Geoinformation for Disaster Management Conference, Gi4DM 2018 -- 18 March 2018 through 21 March 2018 -- -- 135177, The Black Sea Region is one of the most landslide prone area due to the high slope gradients, heavy rainfall and highly weathered hillslope material conditions in Turkey. The landslide occurrences in this region are mainly controlled by the hydro-climatic conditions and anthropogenic activities. Rapid regional landslide inventory mapping after a major event is main difficulties encountered in this densely vegetated region. However, landslide inventories are first step and necessary for susceptibility assessment since considering the principle that the past is the key to the future thus, future landslides will be more likely to occur under similar conditions, which have led to past and present instability. In this respect, it is important to apply rapid mapping techniques to create regional landslide inventory maps of the area. This study presents the preliminary results of the semi-automated mapping of landslides from unmanned aerial vehicles (UAV) with object-based image analysis (OBIA) approach. Within the scope of the study, ultra-high resolution aerial photographs were taken with fixed wing UAV system on Aug 17, 2017 in the landslide zones which are triggered by the prolonged heavy rainfall event on August 12-13, 2016 at Bartin Kurucaşile province. 10 cm resolution orthomosaic and Digital Surface Model (DSM) data of the area were produced by processing the obtained photographs. A test area was selected from the overall research area and semi-automatic landslide detection was performed by applying object-based image analysis. OBIA has been implemented in three steps: image segmentation, image object metric calculation and classification. The accuracy of the resulting maps is assessed by comparisons with expert based landslide inventory map of the area. As a result of the comparison, 80% of the 240 landslides in the area were detected correctly, Firat University Scientific Research Projects Management Unit: 1608F607, This study was supported by Anadolu University Scientific Research Projects Commission under the grant number 1608F607. The authors would like to thank İsmail Bulut from General Directorate of Combating Desertification and Erosion for his valuable suggestions on the site-selection.

Earthquake‐induced landslide (EQIL) inventories are essential tools to extend our knowledge of the relationship between earthquakes and the landslides they can trigger. Regrettably, such inventories are difficult to generate and therefore scarce, and the available ones differ in terms of their quality and level of completeness. Moreover, access to existing EQIL inventories is currently difficult because there is no centralized database. To address these issues, we compiled EQIL inventories from around the globe based on an extensive literature study. The database contains information on 363 landslide‐triggering earthquakes and includes 66 digital landslide inventories. To make these data openly available, we created a repository to host the digital inventories that we have permission to redistribute through the U.S. Geological Survey ScienceBase platform. It can grow over time as more authors contribute their inventories. We analyze the distribution of EQIL events by time period and location, more specifically breaking down the distribution by continent, country, and mountain region. Additionally, we analyze frequency distributions of EQIL characteristics, such as the approximate area affected by landslides, total number of landslides, maximum distance from fault rupture zone, and distance from epicenter when the fault plane location is unknown. For the available digital EQIL inventories, we examine the underlying characteristics of landslide size, topographic slope, roughness, local relief, distance to streams, peak ground acceleration, peak ground velocity, and Modified Mercalli Intensity. Also, we present an evaluation system to help users assess the suitability of the available inventories for different types of EQIL studies and model development.

WOS: 000424358600001, This paper presents the results of a geomorphological investigation of Mount Akdag landslide complex, located on the Western Taurus range, SW Turkey. The landslide, resulting in the collapse of a 5 km segment of Mount Akdag, covers an area of 9.8 km(2) and has a volume of about 3 x 108 m(3). The 1: 15,000 scale geomorphological map (Main Map) of the Akdag landslide presented here was produced from a combination of field studies and unmanned aerial vehicle-based DSM, including aerial-photo interpretation and geomorphometric analysis of landforms. The map shows contemporary and past geomorphological hillslope processes and landforms associated with the landslide complex. Our new map suggests that the Akdag landslide was a complex rockslide consisting of multiple landslide types that developed through various movements. We found that 43 highly developed and active mass movements on the main body and side slopes of the landslide made a major contribution to the sediment input. The sediment produced in the upstream severely threatens the touristic site of Saklikent Canyon in the lower catchment via an increasing magnitude of torrents. This map will provide important data for local and national administrators regarding monitoring up-to-date mass-wasting processes, understanding the landslide and its related sedimentary dynamics, and implementing risk assessment and necessary precautions., Scientific Research Projects Coordination Unit of Istanbul University [ONAP 33594 [ID734]]; Anadolu University Scientific Research Projects Commission [1210E156], The authors would like to thank the Scientific Research Projects Coordination Unit of Istanbul University (Project number: ONAP 33594 [ID734]) and Anadolu University Scientific Research Projects Commission for their financial support through grant no. 1210E156.

In areas where groundwater overexploitation occurs, land subsidence triggered by aquifer compaction is observed, resulting in high socio-economic impacts for the affected communities. In this paper, we focus on the Konya region, one of the leading economic centers in the agricultural and industrial sectors in Turkey. We present a multi-source data approach aimed at investigating the complex and fragile environment of this area which is heavily affected by groundwater drawdown and ground subsidence. In particular, in order to analyze the spatial and temporal pattern of the subsidence process we use the Small BAseline Subset DInSAR technique to process two datasets of ENVISAT SAR images spanning the 2002-2010 period. The produced ground deformation maps and associated time-series allow us to detect a wide land subsidence extending for about 1200 km2 and measure vertical displacements reaching up to 10 cm in the observed time interval. DInSAR results, complemented with climatic, stratigraphic and piezometric data as well as with land-cover changes information, allow us to give more insights on the impact of climate changes and human activities on groundwater resources depletion and land subsidence. © 2017, by the authors; licensee MDPI, Basel, Switzerland.

In this study we develop an empirical method to estimate the volume threshold for predicting coseismic landslide dam formation using landscape parameters obtained from digital elevation models (DEMs). We hypothesize that the potential runout and volume of landslides, together with river features, determine the likelihood of the formation of a landslide dam. To develop this method, a database was created by randomly selecting 140 damming and 200 non-damming landslides from 501 landslide dams and > 60 000 landslides induced by the Mw 7.9 2008 Wenchuan earthquake in China. We used this database to parameterize empirical runout models by stepwise multivariate regression. We find that factors controlling landslide runout are landslide initiation volume, landslide type, internal relief (H) and the H/L ratio (between H and landslide horizontal distance to river, L). In order to obtain a first volume threshold for a landslide to reach a river, the runout regression equations were converted into inverse volume equations by taking the runout to be the distance to river. A second volume threshold above which a landslide is predicted to block a river was determined by the correlation between river width and landslide volume of the known damming landslides. The larger of these two thresholds was taken as the final damming threshold. This method was applied to several landslide types over a fine geographic grid of assumed initiation points in a selected catchment. The overall prediction accuracy was 97.4% and 86.0% for non-damming and damming landslides, respectively. The model was further tested by predicting the damming landslides over the whole region, with promising results. We conclude that our method is robust and reliable for the Wenchuan event. In combination with pre-event landslide susceptibility and frequency–size assessments, it can be used to predict likely damming locations of future coseismic landslides, thereby helping to plan emergency response. Copyright © 2014 John Wiley & Sons, Ltd.

The 2002 M-w 7.9 Denali Fault earthquake, Alaska, provides an unparalleled opportunity to investigate in quantitative detail the regional hillslope mass-wasting response to strong seismic shaking in glacierized terrain. We present the first detailed inventory of similar to 1580 coseismic slope failures, out of which some 20% occurred above large valley glaciers, based on mapping from multi-temporal remote sensing data. We find that the Denali earthquake produced at least one order of magnitude fewer landslides in a much narrower corridor along the fault ruptures than empirical predictions for an M 8 earthquake would suggest, despite the availability of sufficiently steep and dissected mountainous topography prone to frequent slope failure. In order to explore potential controls on the reduced extent of regional coseismic landsliding we compare our data with inventories that we compiled for two recent earthquakes in periglacial and formerly glaciated terrain, i.e. at Yushu, Tibet (M-w 6.9, 2010), and Aysen Fjord, Chile (2007 M-w 6.2). Fault movement during these events was, similarly to that of the Denali earthquake, dominated by strike-slip offsets along near-vertical faults. Our comparison returns very similar coseismic landslide patterns that are consistent with the idea that fault type, geometry, and dynamic rupture process rather than widespread glacier cover were among the first-order controls on regional hillslope erosional response in these earthquakes. We conclude that estimating the amount of coseismic hillslope sediment input to the sediment cascade from earthquake magnitude alone remains highly problematic, particularly if glacierized terrain is involved. (C) 2014 Elsevier Ltd. All rights reserved.

The Black Sea Region of Turkey is one of the most landslide prone areas due to its high slope topography, heavy rainfall, and highly weathered hillslope material conditions. Preparation of landslide inventory maps is the first step in producing landslide susceptibility maps. Ground-based methods for mapping landslide occurrences are time-consuming and expensive. Additionally, landslide mapping based on satellite imageries and aerial photographs has some limitations, including climatic conditions, cost, and limited repetitive measurement capacity. Visual interpretation-based landslide mapping, which is based on satellite imageries and aerial photographs, is a time-consuming procedure that requires an experience-based expert opinion. Therefore, the data acquisition based on unmanned aerial vehicle (UAV) and landslide event inventory maps using an object-based classification approach can be superior to other methods in terms of speed and cost. In this study, we developed a semiautomatic model using object-based image analyses for rapid mapping of shallow landslides from the data obtained from UAVs after major landslide events in the Black Sea Region of Turkey. For this purpose, two test sites—Kurucasile (Bartin) and Cayeli (Rize)—were selected. Landslide mapping models were developed in the investigation sites, and the performance of the models was evaluated. The landslides' data obtained with the developed models were compared to the landslides' data produced by the experts. The comparison process revealed that landslides mapped by using UAV data have an accuracy rate higher than 86% according to the number of landslides and 83% according to the landslide area.

A precise, accurate and complete landslide inventory is indispensable for the establishment of reliable landslide susceptibility and hazard maps. In the preparation of landslide inventories, dense vegetation cover is the major obstacle that confounds the topographic signature of landslides. Today, the growing usage of light detection and ranging (LiDAR) technology in the field of geoscience illuminates the mystery of the landscapes cloaked in dense vegetation by providing a new visual acuity to researchers. Turkey, similar to many mid-latitude mountainous countries suffering from landslides, is still continuing landslide inventory mapping with conventional methods such as aerial photo interpretations (API), although 68% of landslide events occurred under dense forest cover in the north of the country. Despite the country-wide medium-scaled landslide catalogs, the number and the abundance of landslides covered by forests remain largely unknown regardless of the increasing availability of high-resolution remote-sensing data. From these motivating insights and drawbacks, the study's focus is to evaluate the capability of mapping landslides by visually analyzing airborne LiDAR DTM derivatives and compare the results with the 1:25,000 scaled API-based inventory to understand the potential contribution of LiDAR technology in Europe's deadliest country (Turkey) in terms of landslides. The landslide mapping results for a test area located in the densely forested Ulus Basin, Western Black Sea region, reveal that the extent and the number of the mapped landslides (n = 902) from the airborne LiDAR data are much higher than those of the available API-based landslide inventory, which includes 67 landslides. Comparative analysis on topographic signatures of landslides is also underpinned on a distinct discrepancy of slope failure diagnostic indicators captured by the two available DEMs derived from LiDAR and 1:25,000 scaled topographic maps, emphasizing the success rate of LiDAR-derived DEM in recognition and mapping landslides with higher precision and accuracy. Together with all analysis results, the landslide recognizability degree assessment based on the forest cover percentage and slope height differences of landslides highlights that the LiDAR data enable landslides to be defined up to 100 m2, while this value is 20 times lower than that of estimated values from API-based landslides. In addition to underlining the contribution of LiDAR technology to the recognition of landslides hidden under dense vegetation, these findings stress the importance of LiDAR data on complete and accurate landslide inventories for production of reliable susceptibility and hazard maps and further a better understanding of the landslide processes and reducing related losses.

Earthquake-triggered landslide dams are potentially dangerous disrupters of water and sediment flux in mountain rivers, and capable of releasing catastrophic outburst flows to downstream areas. We analyze an inventory of 828 landslide dams in the Longmen Shan mountains, China, triggered by the M-w 7.9 2008 Wenchuan earthquake. This database is unique in that it is the largest of its kind attributable to a single regional-scale triggering event: 501 of the spatially clustered landslides fully blocked rivers, while the remainder only partially obstructed or diverted channels in steep watersheds of the hanging wall of the Yingxiu-Beichuan Fault Zone. The size distributions of the earthquake-triggered landslides, landslide dams, and associated lakes (a) can be modeled by an inverse gamma distribution; (b) show that moderate-size slope failures caused the majority of blockages; and (c) allow a detailed assessment of seismically induced river-blockage effects on regional water and sediment storage. Monte Carlo simulations based on volumetric scaling relationships for soil and bedrock failures respectively indicate that 14% (18%) of the estimated total coseismic landslide volume of 6.4 (14.6) x 10(9) m(3) was contained in landslide dams, representing only 1.4% of the >60,000 slope failures attributed to the earthquake. These dams have created storage capacity of similar to 0.6x 10(9) m(3) for incoming water and sediment. About 25% of the dams containing 2% of the total river-blocking debris volume failed one week after the earthquake; these figures had risen to 60% (similar to 20%), and >90% (>90%) within one month, and one:year, respectively, thus also emptying similar to 92% of the total potential water and sediment storage behind these, dams within one year following the earthquake. Currently only similar to 0.08 x 10(9) m(3) remain available as natural reservoirs for storing water and sediment, while similar to 0.19 x 10(9) m(3), i.e. about a third of the total river-blocking debris volume, has been eroded by rivers. Dam volume and upstream catchment area control to first order the longevity of the barriers, and bivariate domain plots are consistent with the observation that most earthquake-triggered landslide dams were ephemeral. We conclude that the river-blocking portion of coseismic slope failures disproportionately modulates the post-seismic sediment flux in the Longmen Shan on annual to decadal timescales.

Onshore and offshore seismic and geologic-morphologic evidence from the wider region of the Canakkale Basin indicates that this area has been widely exposed to transpressional tectonism, which already commenced in the Pliocene. During this transpressional tectonism, the Gelibolu Fault and the Anafartalar Shear Zone on the Gelibolu Peninsula, as well as the Bozcaada-Biga Shear Zone on the Biga Peninsula were activated. As a consequence, the northern part of the Gelibolu Peninsula, and a broad zone between Bozcaada Island and the Karaburun Peninsula were uplifted to form the northern and southern boundaries of the Canakkale Basin, respectively. This remained a low-elevation intermontane basin between these two highlands. The original morphology of the Canakkale Basin may have developed as a coastal and shelf section of the large extensional Marmara Sea Basin at the end of the Late Miocene. This tectonic phase was followed in the Pliocene by the transpressional tectonism of the North Anatolian Fault Zone, which destroyed the initial morphology and formed the present V-shaped basin. The activity of the Gelibolu Fault and the Anafartalar Shear Zone along the northern boundary of the Canakkale Basin ended in the late Pleistocene with the initiation of the northern segment of the North Anatolian Fault Zone. The tectonism along the northern boundary of the Canakkale Basin thus shifted from a transpressional to a transtensional regime. Seismic data indicate that the Bozcaada-Biga Shear Zone continues to be active to the present day.

Projected scenarios of climate change involve general predictions about the likely changes to the magnitude and frequency of landslides, particularly as a consequence of altered precipitation and temperature regimes. Whether such landslide response to contemporary or past climate change may be captured in differing scaling statistics of landslide size distributions and the erosion rates derived thereof remains debated. We test this notion with simple Monte Carlo and bootstrap simulations of statistical models commonly used to characterize empirical landslide size distributions. Our results show that significant changes to total volumes contained in such inventories may be masked by statistically indistinguishable scaling parameters, critically depending on, among others, the size of the largest of landslides recorded. Conversely, comparable model parameter values may obscure significant, i.e. more than twofold, changes to landslide occurrence, and thus inferred rates of hillslope denudation and sediment delivery to drainage networks. A time series of some of Earth's largest mass movements reveals clustering near and partly before the last glacial-interglacial transition and a distinct step-over from white noise to temporal clustering around this period. However, elucidating whether this is a distinct signal of first-order climate-change impact on slope stability or simply coincides with a transition from short-term statistical noise to long-term steady-state conditions remains an important research challenge. Copyright © 2011 John Wiley & Sons, Ltd.

This paper presents the preliminary results of an extensive study of the mapping the distribution of landslides triggered by the Wenchuan earthquake in Sichuan Province, China, on 12 May 2008. An extensive landslide interpretation was carried out using a large set of optical high resolution satellite images (e.g. ASTER, ALOS, Cartosat-1, SPOT-5 and IKONOS) as well as air photos for both the pre- and post-earthquake situation. Landslide scarps were mapped as points using multi-temporal visual image interpretation taking into account shape, tone, texture, pattern, elevation and ridge and valley orientation. Nearly 60,000 individual landslide scarps were mapped. The landslide distribution map was compared with the inventory map that was prepared directly after the earthquake, which contains about 11,000 individual landslide points, through the calculation of normalized landslide isopleths maps. Remarkable differences were observed, as the earlier inventory mapping did not consider the pre-earthquake situation and did not consider all individual landslides. As part of the landslide inventory, landslides were identified that had blocked the drainage and had formed landslide dams. The landslide distribution was compared with a number of aspects, such as the seismic parameters (distance to epicenter, distance to fault rupture, co-seismic fault geometry and co-seismic slip distribution), and geology. The most remarkable correlation found was with the co-seismic slip distribution and the fault geometry. Landslide distribution in the section of the fault that had mainly a thrust component with low angle fault plane was found to be much higher than the sections that had steeper fault angles and a major strike slip component.

Seismic and multi-beam bathymetric data from the northern shelf and slope of the Cinarcik Basin, which is generated by the North Anatolian Fault Zone (NAFZ) located in the easternmost basin in the Marmara Sea, were re-interpreted to better understand the future sub-marine landslide susceptibility. Seismic data indicate that upper surface of the sub-marine extension of the Paleozoic rocks has an NNE–SSW oriented basin and a ridge type morphology controlled by the secondary faults of the NAFZ. Basins are fulfilled by Plio-Quaternary sediments, which are cut by strike-slip faults on the shelf and slope. The thickness of basin deposits reaches up to 130 m toward the linear northern slope of the Cinarcik Basin. A relatively recent sub-marine landslide, the Tuzla Landslide, cuts the slope of the Cinarcik Basin. The detailed morphological investigation indicates that the Tuzla Landslide is a deep-seated rotational landslide, which was likely triggered by activity of the NAFZ. Morphological analyses also indicate that the thick Plio-Quaternary deposits on the Paleozoic basement slid during the Tuzla Landslide event. This landslide is considered as a key event to understand the dynamics of the potential landslides on the northern shelf and slope of the Cinarcik Basin. Two areas locating on the eastern and the western sides of the Tuzla Landslide are considered as the potential areas for future sliding due to similarities of geological and geomorphological features with the Tuzla Landslide such as similar thick Plio-Quaternary deposits, similar slope morphology, and similar fault activity cutting the sediments. Considering this information, the purposes of the present study are to determine the dynamics of the possible landslide areas and to discuss their effects on the sub-marine morphology. In the light of the interpretations, the amounts of possible displaced material are obtained. Three different landslide scenarios due to possible slide surfaces for future landslides are developed and assessed. The first scenario is sliding of the sediments at the shelf break. The third scenario is a mass movement of almost whole basin deposits on the Paleozoic rocks. The latter one is evaluated as less important because of the volume of the displaced material, and the latter one is accepted as lowest possible event. Among the scenarios, the second scenario is accepted as the most critical and possible because of the amount of the slipped material and existence of faults rupture, which is considered as further sliding surfaces. These landslides will result in important changes in shelf, slope and basin floor in the study area.

In the international literature, although considerable amount of publications on the landslide susceptibility mapping exist, geomorphology as a conditioning factor is still used in limited number of studies. Considering this factor, the purpose of this article paper is to implement the geomorphologic parameters derived by reconstructed topography in landslide susceptibility mapping. According to the method employed in this study, terrain is generalized by the contours passed through the convex slopes of the valleys that were formed by fluvial erosion. Therefore, slope conditions before landsliding can be obtained. The reconstructed morphometric and geomorphologic units are taken into account as a conditioning parameter when assessing landslide susceptibility. Two different data, one of which is obtained from the reconstructed DEM, have been employed to produce two landslide susceptibility maps. The binary logistic regression is used to develop landslide susceptibility maps for the Melen Gorge in the Northwestern part of Turkey. Due to the high correct classification percentages and spatial effectiveness of the maps, the landslide susceptibility map comprised the reconstructed morphometric parameters exhibits a better performance than the other. Five different datasets are selected randomly to apply proper sampling strategy for training. As a consequence of the analyses, the most proper outcomes are obtained from the dataset of the reconstructed topographical parameters and geomorphologic units, and lithological variables that are implemented together. Correct classification percentage and root mean square error (RMSE) values of the validation dataset are calculated as 86.28% and 0.35, respectively. Prediction capacity of the different datasets reveal that the landslide susceptibility map obtained from the reconstructed parameters has a higher prediction capacity than the other. Moreover, the landslide susceptibility map obtained from the reconstructed parameters produces logical results.

Structural, mass-wasting and sedimentation processes along an active dextral shear zone beneath the Gulf of Saros and the NE Aegean Sea were investigated on the basis of new high-resolution swath bathymetric data and multi-channel seismics. A long history of dextral shearing operating since the Pliocene culminated in the formation of a NE-SW-trending, ca. 800-m-deep basin (the so-called inner basin) in this region, which is bordered by a broad shelf along its northern and eastern sides and a narrow shelf at the southern side. The western extension of the North Anatolian Fault Zone (the Ganos Fault) cuts the eastern shelf along a narrow deformation zone, and ends sharply at the toe of the slope, where the strain is taken up by two NE-SW-oriented fault zones. These two fault zones cut the basin floor along its central axis and generate a new, Riedel-type pull-apart basin (the so-called inner depression). According to the bathymetric and seismic data, these basin boundary fault zones are very recent features. The northern boundary of the inner depression is a through-going fault comprising several NE-SW- and E-W-oriented, overlapping fault segments. The southern boundary fault zone, on the other hand, consists of spectacular en-echelon fault systems aligned in NE–SW and WNW–ESE directions. These en-echelon faults accommodate both dextral and vertical motions, thereby generating block rotations along their horizontal axis. As the basin margins retreat, the basin widens continuously by mass-wasting of the slopes of the inner basin. The mass-wasting, triggered by active tectonics, occurs by intense landsliding and channel erosion. The eroded material is transported into the deep basin, where it is deposited in a series of deep-sea fans and slumps. The high sedimentation rate is reflected in an over 1,500-m-thick basin fill which has accumulated in Pliocene–Quaternary times.

Seismic profiling, bathymetric and physical oceanographic data collected from the Canakkale Strait revealed that the morphological evolution of the strait has been controlled by tectonic activity, and sediment erosion and deposition. Sediments in the strait have been sourced mostly by rivers draining the Biga Peninsula during lowstand periods. In highstand periods, by contrast, deposits in the strait were reworked by currents. The seafloor morphology of the Canakkale Strait is also controlled by a sequence of factors ranging from tectonics to current erosion and deposition. Channel deposits overlying the basement are being eroded at the narrower, meandering central section of the strait (the Nara Passage) due to high current velocities. The eroded sediments are deposited in the relatively linear and wider, northern and southern sectors of the strait exposed to low current velocities. As a result, the high-energy areas are more deeply incised due to the erosion, whereas deposition elevates the seafloor in the areas exposed to lower current energy. Three strike-slip faults, which possibly relate to the activity of the North Anatolian Fault Zone, are responsible for the irregular shape of the strait and this, in turn, controls the current velocity along the strait. The high-energy conditions probably commenced with the latest invasion of Mediterranean waters some 12 ka b.p., and have continued as a two-layered current system to the present day.

The Quaternary evolution of the Gulf of Izmit, situated on the tectonically active North Anatolian Fault Zone (NAFZ), was investigated using seismic reflection, paleontologic, and sediment textural data. On the basis of seismic stratigraphic and sedimentologic-paleontologic interpretations, four depositional units were distinguished within the Plio-Quaternary sequence of the Gulf of Izmit. According to these data, Plio-Quaternary deposits supplied from the northern terrestrial area started to accumulate during a progradational phase, in a south-facing half-graben. A coarse-grained sedimentary unit prograding into the gulf from the south since 200 ka indicates a dramatic variation in the evolution of the gulf, with the initiation of a new strike-slip fault of the NAFZ and a corresponding uplift of the Armutlu Peninsula in the south of the gulf. During the evolution of this fault from a wide shear zone consisting of right-stepped strike-slip faults and pull-apart basins to a localized principal fault zone, sediments were deposited under the influence of northerly prograding terrestrial and shallow-marine conditions due to relative sea-level fluctuations in the Marmara Sea. During this period, the Gulf of Izmit was invaded mainly by Mediterranean and partly by Black Sea waters. In the latest glacial period, shallow areas in the gulf became subaerially exposed, whereas the central and western sub-basins of the gulf turned into lakes. The present evolution of the Gulf of Izmit is controlled by the after effects of the new rupture of the NAFZ and the estuarine nature of the gulf environment. © Springer-Verlag 2007. Istanbul Üniversitesi Acknowledgements This paper is based on the MSc thesis of Erdal Dolu. We gratefully acknowledge the captains and crews as well as the scientists and technicians of the R/V Arar of Istanbul University, and the R/V TCG Çubuklu of the Turkish Navy, Department of Navigation, Hydrography and Oceanography, for their help during the bathymetric and seismic surveys. We thank Dr. Timur Ustaömer for his very valuable criticism. We also acknowledge the helpful comments of two anonymous referees and the editor of the journal, Dr. Burg Flemming, to improve the paper.

The Strait of Istanbul (SoI) (Bosphorus) is a narrow valley, which has evolved tectonically from a stream, and in which thick sediment deposits have accumulated in the course of its evolution. Detailed seismic and multi-beam bathymetric data have revealed that the upper parts of the deeper channel deposits consist of parallel strata, which have mostly been eroded subsequently to their deposition. The resulting erosion surface is represented by the present channel floor in the strait, the estimated volume of the eroded material being approximately 2×108 m3 . Erosion rate and seafloor morphology indicate that the flow direction was from the south to the north. This inner channel may have been formed by an abrupt flooding of the Black Sea by Mediterranean waters at the beginning of the latest connection between the Marmara and the Black seas. Subsequently, the Mediterranean bottom current of the modern two-way flow system, which was established at about 5–4 ka b.p., has given the latest shape to the strait floor.

The 12 January 2010 Mw 7.0 Haiti earthquake occurred in a complex deformation zone at the boundary between the North American and Caribbean plates. Combined geodetic, geological and seismological data posited that surface deformation was driven by rupture on the Leogâne blind thrust fault, while part of the rupture occurred as deep lateral slip on the Enriquillo–Plantain Garden Fault (EPGF). The earthquake triggered > 4490 landslides, mainly shallow, disrupted rock falls, debris-soil falls and slides, and a few lateral spreads, over an area of ~ 2150 km2. The regional distribution of these slope failures defies those of most similar earthquake-triggered landslide episodes reported previously. Most of the coseismic landslides did not proliferate in the hanging wall of the main rupture, but clustered instead at the junction of the blind Leogâne and EPGF ruptures, where topographic relief and hillslope steepness are above average. Also, low-relief areas subjected to high coseismic uplift were prone to lesser hanging wall slope instability than previous studies would suggest. We argue that a combined effect of complex rupture dynamics and topography primarily control this previously rarely documented landslide pattern. Compared to recent thrust fault-earthquakes of similar magnitudes elsewhere, we conclude that lower static stress drop, mean fault displacement, and blind ruptures of the 2010 Haiti earthquake resulted in fewer, smaller, and more symmetrically distributed landslides than previous studies would suggest. Our findings caution against overly relying on across-the-board models of slope stability response to seismic ground shaking.

Landslide dams caused by earthquakes are extremely hazardous disruptions of the flow of water and sediment in mountain rivers, capable of delivering large outburst floods that may devastate downstream areas. We analyzed a unique inventory of 828 landslide dams triggered by the M w 7.9 2008 Wenchuan tectonic earthquake, China, constituting ∼1.4% of the >60,000 coseismic slope failures mapped and attributed to this event. While 501 landslides blocked the rivers completely, the remainder caused only partial damming or channel diversion. The spatial distribution of landslide dams follows the same trend of that of the total landslide distribution, with landslide dams being most abundant in the steep watersheds of the hanging wall of the Yingxiu-Beichuan Thrust Fault, and in the northeastern part of the strike-slip fault near Qingchuan. Besides the co-seismic landslide density, the river width also played a key role in determining the landslide dam formation. Narrow rivers are more prone to be dammed than the wide rivers. The correlation between river width and landslide dam volume follows a linear relation, which can be used to roughly estimate the dam formation possibility. However, the applicability of this correlation needs to be validated in other regions. The decay (failure) rate of dams, defined here as the percentage of the number and area of landslide dams that have failed over time, shows that ∼25% of dams accounting for ∼30% of total landslide dam area failed one week after the earthquake. These percentages increased to ∼60% within 1 month, and to >90% within 1 year. The geomorphometric parameters were analyzed, revealing power-law relations between landslide area and dam width, landslide source area and dam area, as well as lake area and lake volume. The inventory presented in this study will enrich the worldwide earthquake-induced landslide dam database and will also contribute to a better understanding of the post-earthquake dam decay.

WOS: 000288802500012, The Cappadocia region has unique geomorphological features resulting from differential erosional processes which make it very attractive to tourists. Besides the fairy chimneys for which the area is best known, there are also impressive buttes and mesas. Buttes and mesas are formed in regions having flat-lying strata in which the uppermost levels are composed of well-cemented limestones and granular ignimbrites, whereas the lower parts and slopes consist of low-durability tuff and ignimbrites. This durability difference results in serious rockfall events. This study involves two-dimensional rockfall analyses in and near the Avanos, Zelve, and Cavusini areas, where volcano-sedimentary units of Neogene age outcrop, to provide a rockfall hazard map in which areas of tourism activity are also considered.

Seismic and bathymetric data from the Canakkale Strait and its extensions onto the shelves of the Marmara and Aegean seas indicate that the strait was formed mainly by an erosional event. Four seismic units are observed on seismic profiles. The lower two of these (units 4 and 3) constitute the basement of a regionally widespread erosional unconformity (ravinement), which developed during marine isotope stage 2 (MIS 2). The two upper units (units 2 and 1), which overlie the ravinement surface, form a higher-order sequence. Sequence stratigraphic analysis indicates that units 2 and 1 deposited as lowstand and highstand systems tracts respectively, since the end of MIS 2. The transgressive systems tract is represented by a major erosional event which occurred throughout the Canakkale sill area when the Mediterranean-Marmara Sea connection and, hence, the Canakkale Strait was formed. The existence of the erosive Sarkoy Canyon along the shelf edge of the southern Marmara Sea demonstrates that the flow direction causing the erosion was from south to north, thus proving that it was produced by Mediterranean water flowing over the sill into the Marmara Sea basin.