Search

Cliff bands are common in drylands and their evolution is often influenced by hydrogeomorphic processes. It has been previously suggested that cliff retreat patterns and morphology are affected by the properties and frequency-magnitude relations of rainstorms. However, basic questions on this topic persist because landscape evolution models typically do not account for the surface processes like runoff generation and sediment transport that occur under short-duration (sub-hourly) intense rainfall. Here we test the hypothesis that changes in rainstorm properties can systematically alter cliff retreat patterns and morphology. We developed a novel numerical model that simulates the response of cliffs and associated sub-cliff slopes to various rainstorm regimes to (1) identify dominant cliff morphologies, and (2) examine if extreme rainstorm properties are encoded in the topography. The new model utilizes the Landlab modeling toolkit and includes an explicit novel representation of surface processes that occur during short-duration rainstorms, including cliff-weathering, infiltration, runoff generation, clast fragmentation, and size-dependent sediment transport. Using a suite of numerical experiments, we vary model parameters and rainfall types and simulate changes in cliff retreat patterns and morphology. Our model results agree well with analytical predictions for cliff morphology under a control case of no transport on the sub-cliff slope, indicating a good representation of processes. Furthermore, sensitivity analyses on cases where sediment transport is explicitly included show that cliff evolution is highly dependent on both the grain size of sediment derived from the cliff and the rainfall intensities. These two factors can alter retreat patterns and determine whether and how fast the cliff can be buried under its own sediment. Numerical experiments based on rainfall and field measurements from the central Negev desert (eastern Mediterranean) demonstrate that including the dynamics of high-intensity rainfall and sediment grain size can help explain observed topographic trends. In addition, for a given imposed storm depth, we find that the rainstorm intensities pattern strongly influences both the cliff retreat and its morphology. Short rainstorms with higher intensities are much more erosive than longer storms with lower intensities. This latter case frequently triggers cliff burying. Taken as a whole, our results demonstrate that cliff evolution and morphology are significantly affected by storm-scale sediment transport dynamics and thus highlight the importance of incorporating high-resolution rainfall forcing into landscape evolution models of dryland landforms.

Deciphering aspect-related hillslope asymmetry can enhance our understanding of the influence of climate on Earth’s surface morphology and the linkage between topographic morphology and erosion processes. Although hillslope asymmetry is documented worldwide, the role of microclimatic factors in the evolution of dryland cliffs has received little attention. Here, we address this gap by quantifying aspect-dependent bedrock weathering, slope-rill morphology, and sub-cliff clast transport rates in the hyperarid Negev desert, Israel, based on light detection and ranging (LiDAR)-derived topography, clast-size measurements, and cosmogenic 10Be concentrations. Cliff retreat rates were evaluated using extrapolated profiles from dated talus flatirons and 10Be measurements from the cliff face and sub-cliff sediments. We document systematic, aspect-dependent patterns of south-facing slopes being less steep and finer-grained relative to east- and north-facing aspects. In addition, cliff retreat and clast transport rates on slopes of the south-facing aspect are faster compared to the other aspects. Our data demonstrate that bedrock weathering of the cliff face and the corresponding grain size of cliff-derived clasts delivered to the slopes constitute a first-order control on cliff retreat and sediment transport rates. We demonstrate that the morphology of the cliff and the pattern of bedrock weathering co-vary with the solar radiation flux and hence that cliff evolution in hyperarid regions is modulated by aspect-dependent solar radiation. These results help to better understand interactions between climate and dryland surface processes.

Figures S1–S8, Text S1, Table S1.

Climatic impact on landscape morphology was previously demonstrated under pronounced gradients in average climatic properties such as mean annual precipitation or temperature. However, in arid areas, where both meteorological observations and rainfall measurements are scarce and the latter is meager, short-term and highly variable in space and time, the determination of meaningful “average climatic” conditions and their variability is challenging. Although it is generally acknowledged that surface processes in arid landscapes should be effected by short-duration rainfall intensities and their extremes, the topographic sensitivity to storm-scale properties were rarely quantified. Here, we attempted to bridge this gap by documenting systematic precipitation variations along a 40 km arid escarpment (Ramon crater) in the central Negev desert (Israel) and their associated topographic signature.We used 0.5 m pixel-1 LiDAR-derived topographic data coupled with field measurements to characterize the morphology of cliffs and slopes along the entire Ramon crater. Sub-hourly rainfall intensities were characterized using an 8-year record of high-resolution, convection-permitting, numerical weather model prediction (NWP). Frequency analyses of rainfall intensity and its spatial variation were conducted using a novel statistical method and used to determine runoff and sediment transport along sub-cliff slopes, through grid-based hydrological simulations of synthetic rainstorms with different frequencies.Our results indicate that due to a pronounced decreasing gradient in the number of rain storms per year, the mean annual rainfall decreases from ~100 mm in the southwest (SW) cliff segment to ~40 mm in the northeast (NE) segment. However, in the drier NE cliff segment, extreme rainfall intensities such as the ones occurring during a storm with a 100-year return period are higher. Topographic cliff gradients and the percentage of exposed bedrock over the cliffs increase toward the drier NE cliff section. Sub-cliff slopes in the NE are systematically straighter, shorter, and associated with a smaller clast sizes relative to the wetter (SW) part of the escarpment. Hydrological simulations reveal that under extreme storms, sediment is mobilized by sheetwash on the NE slopes but is less mobile on the wetter SW slopes. In addition, incised gullies and disconnected talus-flatirons are more frequent in the NE and correlate with the higher erosion efficiency of extreme rainstorms in this zone. Our results indicate that significant morphologic differences can be imprinted in arid landforms due to spatial gradients in the properties of extreme rainstorms.

Talus-pediment slopes are a common morphologic feature in arid areas and constitute a prominent runoff and sediment source at the watershed and channel scales. The evolution of talus-pediment sequences (talus flatirons) was often linked to climatic cycles, although the physical processes that may account for such a link remained obscure. Our approach is to integrate field measurements, high-resolution radar rainfall data and numerical modeling to link the frequency of storms and the resulted hillslope runoff and sediment transport. We present a quantitative hydrometeorological analysis of rainstorms and their geomorphic impact, potentially involved in the evolution of arid talus-pediment slopes in the Negev desert (Israel). Artificial rainstorms were designed based on intensity-duration-frequency curves and simulated in the field using a rainfall simulator. Then, the obtained experimental results were up-scaled to the entire slope length using a fully distributed hydrological model. In addition, natural storms and their hydro-geomorphic impacts were monitored using X-band radar and time-lapse cameras.These integrated analyses constrain the rainfall threshold for local runoff generation at rain intensity of 14-22 mm h-1 for a duration of 5 min for the study area conditions. We characterized small-scale runoff-generating convective rain cells using an X-band radar and found that small convective cells (~30 km2), having extremely high internal spatial gradients in rainfall intensity and low velocity (-1), have the potential to generate local hillslope runoff. The frequency of local runoff-producing rainstorms is ~1-3 per year, but most of these storms activate only small parts of the hillslope. Modeling results indicate that a full extent hillslope runoff occurs under much rarer rainstorms of at least 100-years return interval (1% or less). During such rainstorms, the shear stress produced by the runoff flow (sheetwash) is capable of transporting surface clasts at a distance of ~80 m downslope. However, transport of coarse clasts in the upper parts of the slopes is most probably gravitationally controlled. The erosion efficiency of discrete rare events (1% or less) on the lower part of the slopes highlights their potential to trigger incision and lead to cliff dissection. This study results support the hypothesis that a climatic shift in terms of the properties and frequency of extreme rainstorms, rather than the common views of it as changing precipitation means, can play an important role in shaping and in transforming landscapes in such arid setting.

Soil erosion affects agricultural landscapes worldwide, threatening food security and ecosystem viability. In arable environments, soil loss is primarily caused by short, intense rainstorms, typically characterized by high spatiotemporal variability. The complexity of erosive events challenges modeling efforts and explicit inclusion of extreme events in long-term risk assessment is missing. This study is intended to bridge this gap by quantifying the discrete and cumulative impacts of rainstorms on plot-scale soil erosion and providing storm-scale erosion risk analyses for a cropland region in northern Israel. Central to our analyses is the coupling of (1) a stochastic rainfall generator able to reproduce extremes down to 5-minute temporal resolutions; (2) a processes-based event-scale cropland erosion model (Dynamic WEPP, DWEPP); and, (3) a state-of-the-art frequency analysis method that explicitly accounts for rainstorms occurrence and properties. To our knowledge, this is the first study in which DWEPP runoff and soil loss are calibrated at the plot-scale on cropland (NSE is 0.82 and 0.79 for event runoff and sediment, respectively). We generated 300-year stochastic simulations of event runoff and sediment yield based on synthetic precipitation time series. Based on this data, the mean annual soil erosion in the study site is 0.1 kg m-2 [1.1 t ha-1]. Results of the risk analysis indicate that individual extreme rainstorms (>50 return period), characterized by high rainfall intensities (30-minute maximal intensity > ~60 mm h-1) and high rainfall depth (>~200 mm), can trigger soil losses even one order of magnitude higher than the annual mean. The erosion efficiency of these rainstorms is mainly controlled by the short-duration (≤30 min) maximal intensities. The results demonstrate the importance of incorporating the impact of extreme events into soil conservation and management tools. We expect our methodology to be valuable for investigating future changes in soil erosion with changing climate.

Water volume estimates of shallow desert lakes are the basis for water balance calculations, important both for water resource management and paleohydrology/climatology. Water volumes are typically...

Many of the world’s drylands are characterized by interior drainage systems that terminate at shallow desert lakes or playas. Except for episodic flooding these largely ephemeral water bodies, remain mostly dry. Surveying and mapping their respective floor topography in a suitable resolution for calculating water balance is a difficult and laborious task. As this is crucial for water resources management and reconstructing paleoenvironmental conditions, diverse methods and efforts were applied. However, detailed, high-quality bathymetric surveys in such environments are rare and have only been conducted in a few such lakes. This is primarily due to their shallow, low-relief, large areas, and often remote characteristics, which complicate application of conventional topographic surveying techniques. Therefore, satellite-based remote sensing is an essential complementary approach for deriving bathymetry of such lakes.Global digital elevation models, such as NASA’s Shuttle Radar Topography Mission (SRTM) or ASTER’s GDEM, are unsuitable for accurate measurements of these ephemeral lakes, mainly because of their impenetrability to water and their high signal-to-noise ratio in the low-relief environments. Recent studies addressed these complications by combining remote sensing data with local calibrations of in-situ measurements, or alternatively, by relating shoreline altitudes with precise altimetry. This approach requires a spatial interpolation of individual measurements. Therefore, it is prone to errors that demand intensive efforts to be reduced; even then the errors may remain larger than the actual depth of a flooded lake. Moreover, such methods are hard to apply in complex lakes with multiple sub-basins.To tackle these problems, we developed a simple methodology to derive a high-resolution (30 m per pixel) bathymetry of shallow desert lakes. In this new method we combine two sources of data: a >30-yr record of Landsat-derived surface water occurrence data; and accurate high-resolution elevation data, acquired by the NASA’s recently launched ICESat-2 satellite (Ice, Cloud, and Land Elevation Satellite-2). We test the proposed new method over three ephemeral lakes around the world: Lake Eyre, Australia, with its complex shallow lake system, consisting of a few sub-basins; Sabkhat Al-Mellah, Algeria; Lago Coipasa, Bolivia. The accuracy of the resulted bathymetric maps was evaluated using cross-validation of ICESat-2 scans, yielding low RMSD values of ~0.3 m, versus ~2.5 m of the SRTM data (validated through other ICESat-2 scans). At Lago Coipasa, we show that bathymetry was effectively determined even when the lake was full of water (up to a few meters depth). This high-resolution, low-error bathymetry mapping complement other globally available topographic data.