Landscapes on the edge: examining the role of climatic interactions in shaping coastal watersheds using a coastal–terrestrial landscape evolution model

Incised coastal gullies (ICGs) are dynamic features found at the terrestrial‐coastal interface. Their geomorphic evolution is driven by the interactions between processes of fluvial knickpoint migration and coastal cliff erosion. Under scenarios of future climate change the frequency and magnitude of the climatological drivers of both terrestrial (fluvial and hillslope) and coastal (cliff erosion) processes are likely to change, with an adjunct impact on these types of coastal features. Here we explore the response of an incised coastal gully to changes in both terrestrial and coastal climate in order to elucidate the key process interactions which drive ICG evolution. We modify an extant landscape evolution model, CHILD, to incorporate processes of soft‐cliff erosion. This modified version, termed the Coastal‐Terrestrial‐CHILD (CT‐CHILD) model, is then employed to explore the interactions between changing terrestrial and coastal driving forces on the future evolution of an ICG found on the south‐west Isle of Wight, UK. It was found that the magnitude and frequency of storm events will play a key role in determining the future trajectory of ICGs, highlighting a need to understand the role of event sequencing in future projections of landscape evolution. Furthermore, synergistic (positive) and antagonistic (negative) interactions were identified between coastal and terrestrial parameters, such as wave height intensity and precipitation duration, which act to modulate the impact of changes in any one parameter. Of note was the role played by wave height intensity in driving coastal erosion, which was found to play a more important role than sea‐level rise in determining rates of coastal erosion. This highlights the need for a greater focus on wave height in studies of soft‐cliff erosion. Copyright © 2014 John Wiley & Sons, Ltd.     

Emergent characteristics of rockfall inventories captured at a regional scale

High-resolution rockfall inventories captured at a regional scale are scarce. This is partly owing to difficulties in measuring the range of possible rockfall volumes with sufficient accuracy and completeness, and at a scale exceeding the influence of localized controls. This paucity of data restricts our ability to abstract patterns of erosion, identify long-term changes in behaviour and assess how rockfalls respond to changes in rock mass structural and environmental conditions. We have addressed this by developing a workflow that is tailored to monitoring rockfalls and the resulting cliff retreat continuously (in space), in three-dimensional (3D) and over large spatial scales (>10⁴ m). We tested our approach by analysing rockfall activity along 20.5 km of coastal cliffs in North Yorkshire (UK), in what we understand to be the first multi-temporal detection of rockfalls at a regional scale. We show that rockfall magnitude–frequency relationships, which often underpin predictive models of erosion, are highly sensitive to the spatial extent of monitoring. Variations in rockfall shape with volume also imply a systemic shift in the underlying mechanisms of detachment with scale, leading us to question the validity of applying a single probabilistic model to the full range of rockfalls observed here. Finally, our data emphasize the importance of cliff retreat as an episodic process. Going forwards, there will a pressing need to understand and model the erosional response of such coastlines to rising global sea levels as well as projected changes to winds, tides, wave climates, precipitation and storm events. The methodologies and data presented here are fundamental to achieving this, marking a step-change in our ability to understand the competing effects of different processes in determining the magnitude and frequency of rockfall activity and ultimately meaning that we are better placed to investigate relationships between process and form/erosion at critical, regional scales. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Fluvial thermal erosion investigations along a rapidly eroding river bank: application to the Lena River (central Siberia)

In Central Yakutia, frozen river banks are affected by a combination of thermal and mechanical erosion. Exceptional bank retreat of up to 40 m per year is observed. This results from ground thawing produced by heat transfer from the ?ow of water through the frozen ground, followed by mechanical transport of the thawed sediments. A one‐dimensional model is proposed to estimate the thermal erosion ef?ciency. A test of this model is a comparison of results obtained from experiments carried out in a cold room. A hydraulic channel allows measurements of the thaw front propagation, as well as the thermal erosion rate, in simulated ground ice that is subjected to warm water ?ow. Various laboratory simulations demonstrate the validity of the mathematical model for the range of laboratory conditions. A hierarchy of parameters (Reynolds number, water and ground ice temperatures) is proposed to explain the present ef?ciency of thermal erosion along the Siberian rivers. From the characteristics of the Lena River (geometry, temperature and discharge) during the ?ood season, the erosion of banks with different ice content predicted by the model is in agreement with ?eld observations. Copyright © 2003 John Wiley & Sons, Ltd.     

Rates of surface processes on slopes,slope retreat and denudation

Results taken from 270 publications on rates are summarized, and collated with those from 149 publications reviewed previously (Young, 1969, 1974). The data are classified by major climatic zone, normal or steep relief, and consolidated or unconsolidated rocks. Representative rates and their ranges are given for soil creep, solifluction, surface wash, solution (chemical denudation), rock weathering, slope retreat, cliff (free face) retreat, marine cliff retreat, and denudation, the last being compared with representative rates of uplift. Solifluction is of the order of 10 times faster than soil creep, but both cause only very slow ground loss. Solution is an important cause of ground loss for siliceous rocks, on which it may be half as rapid as on limestones. Total denudation, brought about mainly by surface wash, reaches a maximum in the semi-arid and probably also the tropical savanna zones. Acceleration of natural erosion rates by human activities ranges from 2–3 times with moderately intense land use to about 10 times with intensive land use (and considerably higher still where there is recognized accelerated soil erosion). Where there is active uplift, typical rates are of the order of 10 times faster than denudation, although in some high, steep mountain ranges these may approach equality.     

Power law statistics of cliff failures,scaling and percolation

Large cliff failures have often been characterized by their size in terms of the horizontal eroded area at the cliff top, and the maximum local retreat of the coastline. Field studies suggest that the frequencies of these two quantities decrease as power laws of their respective magnitudes, defining two different decay exponents. Moreover, the horizontal area increases as a power law of the maximum local retreat, identifying a third exponent. This observation indicates that the geometry of cliff failures are similar for different magnitudes. It is shown that, under a ‘scaling hypothesis’, the exponents satisfy a precise mathematical relation that relates the exponent of the distributions of magnitudes of events and their geometry – a relation that could be useful to validate the analysis of field catalogs. In order to understand the observed scenario, a numerical model of rocky coast erosion is developed. Despite its minimal character, it reproduces the observed cliff failure statistics. Strikingly, the model precisely reproduces the geometric similarity of cliff failures observed in the field. These results suggest that the theory of percolation, which underlies the erosion model, can possibly be used as a guide to decipher the physics of rocky coast erosion and could provide precise predictions of the statistics of cliff collapses. Copyright © 2015 John Wiley & Sons, Ltd.     

Global distribution of coastal cliffs

Previous studies have estimated that coastal cliffs exist on about 80% of the global shoreline, but have not been validated on a global scale. This study uses two approaches to capture information on the worldwide existence and erosion of coastal cliffs: a detailed literature survey and imagery search, and a GIS-based global mapping analysis. The literature and imagery review show coastal cliffs exist in 93% of the combined recognized independent coastal states and non-independent coastal regions worldwide (total of 213 geographic units). Additionally, cliff retreat rates have been quantified in at least one location within 33% of independent coastal states and 15% of non-independent regions. The GIS-based mapping used the near-global Shuttle Radar Topography Mission 3 arc second digital elevation model and Arctic Coastal Dynamics Database to obtain near-global backshore coastal elevations at 1 km alongshore intervals comprising about 1,340,000 locations (81% of the world vector shoreline). Backshore coastal elevations were compared with the mapped distribution of European coastal cliffs to produce a model training set, and this relationship was extended globally to map the likelihood of coastal cliff locations. About 21% of the transects (17% of the world vector shoreline) were identified as mangroves and eliminated as potential cliff locations. The results were combined with estimates of cliff percentages for Greenland and Antarctica from the literature, extending the global coverage to estimate cliff occurrence across 89% of the world vector shoreline. The results suggest coastal cliffs likely exist on about 52% of the global shoreline. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

The role of beach morphology on coastal cliff erosion under extreme waves

Erosion of hard‐rock coastal cliffs is understood to be caused by a combination of both marine and sub‐aerial processes. Beach morphology, tidal elevation and significant wave heights, especially under extreme storm conditions, can lead to variability in wave energy flux to the cliff‐toe. Wave and water level measurements in the nearshore under energetic conditions are difficult to obtain and in situ observations are rare. Here we use monthly cliff‐face volume changes detected using terrestrial laser scanning alongside beach morphological changes and modelled nearshore hydrodynamics to examine how exposed cliffs respond to changes in extreme wave conditions and beach morphology. The measurements cover the North Atlantic storms of 2013 to 2014 and consider two exposed stretches of coastline (Porthleven and Godrevy, UK) with contrasting beach morphology fronting the cliffs; a flat dissipative sandy beach at Godrevy and a steep reflective gravel beach at Porthleven. Beach slope and the elevation of the beach–cliff junction were found to influence the frequency of cliff inundation and the power of wave–cliff impacts. Numerical modelling (XBeach‐G) showed that under highly energetic wave conditions, i.e. those that occurred in the North Atlantic during winter 2013–2014, with H_s = 5.5 m (dissipative site) and 8 m (reflective site), the combination of greater wave height and steeper beach at the reflective site led to amplified wave run‐up, subjecting these cliffs to waves over four times as powerful as those impacting the cliffs at the dissipative site (39 kWm^‐1 compared with 9 kWm^‐1). This study highlighted the sensitivity of cliff erosion to extreme wave conditions, where the majority (over 90% of the annual value) of cliff‐face erosion ensued during the winter. The significance of these short‐term erosion rates in the context of long‐term retreat illustrates the importance of incorporating short‐term beach and wave dynamics into geomorphological studies of coastal cliff change. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Sea stack formation and the role of abrasion on beach‐mantled headlands

Sea stacks are common and striking coastal landforms, but few details are known about how, how quickly, and under what conditions they form. We present numerical and analytical models of sea stack formation due to preferential erosion along a pre‐existing headland to address these basic questions. On sediment‐rich rocky coasts, as sea cliffs erode and retreat, they produce beach sediment that is distributed by alongshore sediment transport and controls future sea cliff retreat rates. Depending on their width, beaches can encourage or discourage sea cliff erosion by acting either as an abrasive tool or a protective cover that dissipates wave energy seaward of the cliff. Along the flanks of rocky headlands where pocket beaches are often curved and narrow due to wave field variability, abrasion can accelerate alongshore‐directed sea cliff erosion. Eventually, abrasion‐induced preferential erosion can cut a channel through a headland, separating it from the mainland to become a sea stack. Under a symmetrical wave climate (i.e. equal influence of waves approaching the coastline from the right and from the left), numerical and analytical model results suggest that sea stack formation time and plan‐view size are proportional to preferential erosion intensity (caused by, for example, abrasion and/or local rock weakness from joints, faults, or fractures) and initial headland aspect ratio, and that sea stack formation is discouraged when the sediment input from sea cliff retreat is too high (i.e. sea cliffs retreat quickly or are sand‐rich). When initial headland aspect ratio is too small, and the headland is ‘rounded’ (much wider in the alongshore direction at its base than at its seaward apex), the headland is less conducive to sea stack formation. On top of these geomorphic and morphologic controls, a highly asymmetrical wave climate decreases sea stack size and discourages stack formation through rock–sediment interactions. Copyright © 2014 John Wiley & Sons, Ltd.     

Environmental controls on soil frost activity in the Western Cape mountains,South Africa

This study aims to analyse the environmental controls on soil frost processes in the Western Cape mountains of South Africa. Two microclimatic monitoring stations were established on different substrates at about 1900 m a.s.l. recording air and soil temperature, soil moisture and precipitation over periods of five and two years respectively. Other data available are snow cover estimations and soil textural data. Results show the region to experience surficial diurnal frost only. The frequency of effective frost days in the sandstone areas is extremely limited due to insulation by snow cover and vegetation, effectiveness of the zero-curtain effect and high albedo values of the surface. Irrespective of climatic controls, sandstone-derived sediments are found to be too coarse to develop segregation ice. These strata underlie over 90 per cent of the Western Cape mountains over 1000 m a.s.l. Monitoring on shales indicates 12 and 16 diurnal frost cycles for needle-ice growth for 1993 and 1994, respectively. © 1998 John Wiley & Sons, Ltd.     

Influence of lithostratigraphy on the retreat of an unconsolidated sedimentary coastal cliff (St. Lawrence estuary,eastern Canada)

This paper investigates the processes involved in unconsolidated cliff recession using LiDAR surveys (2005, 2010 and 2013) and aerial photographs (1964–2012) at Pointe‐au‐Bouleau, on the north shore of the St. Lawrence estuary, in eastern Canada. The high lithostratigraphic variability of the sediments allowed for the identification of stratigraphic and lithological variables that explain the evolution of coastal cliffs. Space‐for‐time substitution was also used to assess how lithostratigraphy controls the evolution of emerged glaciomarine coastal cliffs over decadal to centennial timescales. This case study presents new quantitative data that contributes to a better understanding of the role of sediment architecture, stratigraphy and geomorphology on coastal evolution. The methodological approach includes the development of a new conceptual model suitable for identifying erosion on cliff coastlines. The high spatial resolution methodology (<5 cm) used herein demonstrates the need for further research using LiDAR data in order to quantify the processes involved in the evolution of coastal cliffs. Copyright © 2015 John Wiley & Sons, Ltd.     

Climate change and soil erosion in a small alpine lake basin on the Loess Plateau,China

Multi‐proxy indices retrieved from sediments in Lake Chaonaqiu, an alpine lake on the western Loess Plateau (LP) of China, were used to reconstruct a precipitation history over the last ~300 years. The results correlate well with records from tree rings and historical documents in neighboring regions. We show that the lake oscillated between two states, i.e. wetter climatic conditions, which favored denser vegetation cover, and promoted weaker catchment soil erosion; and drier climatic conditions, which lead to less vegetation coverage, correlate with stronger surface soil erosion. Several intensive soil erosion events were identified in the sediment cores, and most of these occurred during decadal/multi‐decadal dry periods, and correlate well with flood events documented in historical literature. The results of this study show that soil erosion by flood events is particularly intense during dry periods, and further highlights the role of vegetation cover in the conservation of water and soil in small lake basins on the Chinese LP. Copyright © 2016 John Wiley & Sons, Ltd.     

Influence of rock mass strength on the erosion rate of alpine cliffs

Collapse of cliff faces by rockfall is a primary mode of bedrock erosion in alpine environments and exerts a first‐order control on the morphologic development of these landscapes. In this work we investigate the influence of rock mass strength on the retreat rate of alpine cliffs. To quantify rockwall competence we employed the Slope Mass Rating (SMR) geomechanical strength index, a metric that combines numerous factors contributing to the strength of a rock mass. The magnitude of cliff retreat was calculated by estimating the volume of talus at the toe of each rockwall and projecting that material back on to the cliff face, while accounting for the loss of production area as talus buries the base of the wall. Selecting sites within basins swept clean by advancing Last Glacial Maximum (LGM) glaciers allowed us to estimate the time period over which talus accumulation occurred (i.e. the production time). Dividing the magnitude of normal cliff retreat by the production time, we calculated recession rates for each site. Our study area included a portion of the Sierra Nevada between Yosemite National Park and Lake Tahoe. Rockwall recession rates determined for 40 alpine cliffs in this region range from 0·02 to 1·22 mm/year, with an average value of 0·28 mm/year. We found good correlation between rockwall recession rate and SMR which is best characterized by an exponential decrease in erosion rate with increasing rock mass strength. Analysis of the individual components of the SMR reveals that joint orientation (with respect to the cliff face) is the most important parameter affecting the rockwall erosion rate. The complete SMR score, however, best synthesizes the lithologic variables that contribute to the strength and erodibility of these rock slopes. Our data reveal no strong independent correlations between rockwall retreat rate and topographic attributes such as elevation, aspect, or slope angle. Copyright © 2009 John Wiley & Sons, Ltd.     

Submarine platform development by erosion of a Surtseyan cone at Capelinhos,Faial Island,Azores

Erosion of volcanic islands ultimately creates shallow banks and guyots, but the ways in which erosion proceeds to create them over time and how the coastline retreat rate relates to wave conditions, rock mass strength and other factors are unclear. The Capelinhos volcano was formed in 1957/58 during a Surtseyan and partly effusive eruption that added an ~2.5 km² tephra and lava promontory to the western end of Faial Island (Azores, central North Atlantic). Subsequent coastal and submarine erosion has reduced the subaerial area of the promontory and created a submarine platform. This study uses historical information, photos and marine geophysical data collected around the promontory to characterize how the submarine platform developed following the eruption. Historical coastline positions are supplemented with coastlines interpreted from 2004 and 2014 Google Earth images in order to work out the progression of coastline retreat rate and retreat distance for lava- and tephra-dominated cliffs. Data from swath mapping sonars are used to characterize the submarine geometry of the resulting platform (position of the platform edge, gradient and morphology of the platform surface). Photographs collected during SCUBA and ROV dives on the submarine platform reveal a rugged surface now covered with boulders. The results show that coastal retreat rates decreased rapidly with time after the eruption and approximately follow an inverse power-law relationship with coastal retreat distance. We develop a finite-difference model for wave attenuation over dipping surfaces to predict how increasing wave attenuation contributed to this trend. The model is verified by reproducing the wave height variation over dipping rock platforms in the UK (platform gradient 1.2° to 1.8°) and Ireland (1.8°). Applying the model to the dipping platform around Capelinhos, using a diversity of cliff resistance predicted from known lithologies, we are able to predict erosion rate trends for some sectors of the edifice. We also explore wider implications of these results, such as how erosion creates shallow banks and guyots in reef-less mid-oceanic archipelagos like the Azores. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.     

Defining the moment of erosion: the principle of thermal consonance timing

Geomorphological process research demands quantitative information on erosion and deposition event timing and magnitude, in relation to fluctuations in the suspected driving forces. This paper establishes a new measurement principle – thermal consonance timing (TCT) – which delivers clearer, more continuous and quantitative information on erosion and deposition event magnitude, timing and frequency, to assist understanding of the controlling mechanisms. TCT is based on monitoring the switch from characteristically strong temperature gradients in sediment, to weaker gradients in air or water, which reveals the moment of erosion. The paper (1) derives the TCT principle from soil micrometeorological theory; (2) illustrates initial concept operationalization for field and laboratory use; (3) presents experimental data for simple soil erosion simulations; and (4) discusses initial application of TCT and perifluvial micrometeorology principles in the delivery of timing solutions for two bank erosion events on the River Wharfe, UK, in relation to the hydrograph. River bank thermal regimes respond, as soil temperature and energy balance theory predicts, with strong horizontal thermal gradients (often >1 K cm^?1 over 6·8 cm). TCT fixed the timing of two erosion events, the first during inundation, the second 19 h after the discharge peak and 13 h after re‐emergence from the flow. This provides rare confirmation of delayed bank retreat, quantifies the time‐lag involved, and suggests mass failure processes rather than fluid entrainment. Erosion events can be virtually instantaneous, implying ‘catastrophic retreat’ rather than ‘progressive entrainment’. Considerable potential exists to employ TCT approaches for: validating process models in several geomorphological contexts; assisting process identification and improving discrimination of competing hypotheses of process dominance through high‐resolution, simultaneous analysis of erosion and deposition events and driving forces; defining shifting erodibility and erosion thresholds; refining dynamic linkages in event‐based sediment budget investigations; and deriving closer approximations to ‘true’ erosion and deposition rates, especially in self‐concealing scour‐and‐fill systems. Copyright © 2005 John Wiley & Sons, Ltd.     

Delineation of flood generating processes and their hydrological response

Knowledge about flood generating processes can be beneficial for numerous applications. Especially in the context of climate change impact assessment, daily patterns of meteorological and catchment state conditions leading to flood events (i.e., storylines) may be of value. Here, we propose an approach to identify storylines of flood generation using daily weather and snow cover observations. The approach is tested for and applied to a typical pre‐Alpine catchment in the period between 1961 and 2014. Five precipitation parameters were determined that describe temporal and spatial characteristics of the flood associated precipitation events. The catchment's snow coverage was derived using statistical relationships between a satellite‐derived snow cover climatology and station snow measurements. Moreover, (pre‐) event snow melt sums were estimated using a temperature‐index model. Based on the precipitation and catchment state parameters, 5 storylines were identified with a cluster analysis: These are (a) long duration, low intensity precipitation events with high precipitation depths, (b) long duration precipitation events with high precipitation depths and episodes of high intensities, (c) shorter duration events with high or (d) low precipitation intensity, respectively, and (e) rain‐on‐snow events. The event groups have distinct hydrological characteristics that can largely be explained by the storylines' respective properties. The long duration, high intensity storyline leads to the most adverse hydrological response, namely, a combination of high peak magnitudes, high volumes, and long durations of threshold exceedance. The results show that flood generating processes in mesoscale catchments can be distinguished on the basis of daily meteorological and catchment state parameters and that these process types can explain the hydrological flood properties in a qualitative way. Hydrological simulations of daily resolution can thus be analysed with respect to flood generating processes.     

A comparison of channel erosion in small urban and rural catchments,Armidale, New South Wales

Five years after the cessation of housing construction in a small urban catchment in Armidale, New South Wales, the stability of the channel bed and banks was investigated with the use of erosion pins. After eighteen months of monitoring, the rate of channel bank erosion was found to be 3.6 times greater than that of a nearby rural channel and the rate of knickpoint retreat was 2.4 times greater. There was, however, no evidence to suggest that the urban channel was inherently unstable, but that the increased rate of erosion was the product of changed runoff conditions associated with urban development. Both channels are considered to be in a state of equilibrium.     

Marine control over negative power law scaling of mass wasting events in chalk sea cliffs with implications for future recession under the UKCP09 medium emission scenario

Coastal cliff erosion represents a significant geohazard for people and infrastructure. Forecasting future erosion rates is therefore of critical importance to ensuring the resiliency of coastal communities. We use high precision monitoring of chalk cliffs at Telscombe, UK to generate monthly mass movement inventories between August 2016 and July 2017. Frequency–magnitude analysis of our inventories demonstrate negative power law scaling over 7 orders of magnitude and, for the first time, we report statistically significant correlations between significant wave height (H_s) and power law scaling coefficients (r² values of 0.497 and 0.590 for β and s respectively). Applying these relationships allows for a quantitative method to predict erosion at the site based on H_s probabilities and sea level forecasts derived from the UKCP09 medium emission climate model (A1B). Monte‐Carlo simulations indicate a range of possible erosion scenarios over 70 years (2020–2090) and we assess the impact these may have on the A259 coastal road which runs proximal to the cliffs. Results indicate a small acceleration in erosion compared with those based on current conditions with the most likely scenario at the site being 21.7 m of cliff recession by 2090. However, low‐probability events can result in recession an order of magnitude higher in some scenarios. In the absence of negative feedbacks, we estimate an ~11% chance that the A259 will be breached by coastal erosion by 2090. Copyright © 2018 John Wiley & Sons, Ltd.     

Sediment dynamics below retreating cliffs

The retreat of cliffs may constitute the dominant erosional response to base‐level fall in arid settings underlain by horizontally‐bedded sedimentary rock. These vertical cliffs typically loom above a relatively straight bedrock slope (‘plinth’) that is mantled with a thin layer of sediment and perched near the angle of repose. In detail, a plinth consists of a system of quasi‐parallel ridges and channels. We ask how the sediment supplied from a retreating cliff influences the erosion of the plinth hillslopes and channels, and how this affects the rate of cliff retreat. Motivated by field observations and high‐resolution topographic data from two sites in western Colorado, we develop a two‐dimensional (2D), rules‐based numerical model to simulate the erosion of channels draining a plinth and diffusive erosion of the intervening interfluves. In this model, retreat of a cliffband occurs when the height of the vertical cliff exceeds a threshold due to incision by channels on the plinth below. Debris derived from cliff retreat is distributed over the model plinth according to the local topography and distance from the source. This debris then weathers in place, and importantly can act to reduce local bedrock erosion rates, protecting both the plinth and ultimately the cliff from erosion. In this paper, we focus on two sets of numerical model experiments. In one suite, we regulate the rate of rockfall to limit the cliff retreat rate; in most cases, this results in complete loss of the plinth by erosion. In a second suite, we do not impose a limit on the cliff retreat rate, but instead vary the weathering rate of the rockfall debris. These runs result in temporally steady cliff‐plinth forms and retreat rates; both depend on the weathering rate of the debris. Copyright © 2011 John Wiley & Sons, Ltd.     

Coastal cliff retreat rates at Beit‐Yannay,Israel, in the 20th century

Research indicates that the aeolianite (Kurkar) cliffs along the Israeli Mediterranean coastline have continuously retreated eastward during the last few decades. There seems to be no dispute among Earth scientists regarding the general trend of cliff retreat. However the majority of papers displaying cliff retreat rates are based upon comparison of aerial photographs. Their lack of advanced geometric measurement methods causes a high margin of error. Public attention is focused upon the Beit‐Yannay coastal cliff since private homes are located along the southern section of the cliff crest. The current research compares the historic location of the cliff crest edge at Beit‐Yannay as observed in a series of aerial photographs taken during the period 1918–2000. Quantitative measurement methods included applications of satellite geodesy and digital photogrammetry and mapping. Research results offer quantitative, consecutive and highly accurate data regarding retreat rates over a relatively long period of 82 years. It is concluded that: 1. Annual average cliff retreat rates of the cliff crest is 20 cm/year. 2. Categorization of the study time span reveals periods displaying varying retreat rates such as 27 cm/year during 1918–1946, 21 cm/year during 1946–1973 and 10 cm/year during 1973–2000. 3. Maximum retreat distances of the cliff crest, over the study period were found to be approximately 25 m along the northern, lowest section of the cliff. Minimum distances of 11 m were identi?ed at the highest, southern section of the cliff. 4. The eolianite (Kurkar) cliffs along the Israeli Mediterranean coast throughout the 20th century have been an important source of sediment, contributing approximately 24 × 10⁶ m³ of sediments to the sediment balance of Israeli beaches. Copyright © 2004 John Wiley & Sons, Ltd.     

Lithological controls on soft cliff planshape evolution under high and low sediment availability

This paper addresses a series of geomorphic questions relating to large‐scale (> 1 km), long‐term (100 – 1,000 years) coastal planshape evolution. Previous research on soft‐cliff coasts has recognised the role of protective fronting beach volumes on reducing rates of cliff toe retreat. However, it is the maintenance of this critical threshold that ultimately determines two contrasting modes of shoreline behaviour: Mode A, in which there is little beach sediment and shoreline evolution is controlled by material strength; and, Mode B, when ample beach sediment means that shoreline evolution is controlled by longshore sediment transport. Here we use a numerical model (SCAPE) to investigate temporal and spatial changes in beach volume on a broader range of feedbacks than considered in previous models. The transition between Mode A and Mode B coasts is defined by relative sediment inputs to outputs and used to explore how these contrasting modes control the evolution of an initial linear frontage exhibiting longshore changes in cliff lithology (material resistance and the proportion of beach grade material in the eroded bedrock). Under Mode A, relative changes in material resistance result in long term heterogeneous rates of retreat, which result in the development of persistent headland and embayment features. However, under Mode B, feedbacks between coastal planshape, longshore sediment transport, beach volume and wave energy result in steady state retreat rates regardless of longshore variations in resistance. Results are compared and contrasted to previous simulations and site specific examples and a conceptual model of Mode A and Mode B interactions presented. Copyright © 2015 John Wiley & Sons, Ltd.