A channel head stability criterion for first‐order catchments

The present analysis derives a stability criterion for long‐term equilibrium channel heads. The concept of finite perturbation analysis is presented, during which the surface is subjected to perturbations of a finite amplitude and resulting changes in flow path structure and slope are computed. Based on these quantities the analysis predicts whether the perturbed location is going to erode, be filled in or remain steady. The channel head is defined geometrically as the focus point of converging flow lines at the bottom of hollows. It is demonstrated that stability at the channel head grows out of the competition between the rate of flow path convergence and the degree of profile concavity. Analytical functions are derived to compute channel head‐contributing area and ‐slope, flow path convergence and profile concavity as a function of perturbation depth, distance from the crest and the initial slope. In a numerical model these quantities point to the long‐term equilibrium channel head position, which is shown to depend also on the width to length ratio of hollows. It is also demonstrated that the equilibrium channel head position is sensitive to the base‐level lowering/non‐dimensional slope length ratio and to the slope of the initial topography. Morphometrical measurements both in the field and on simulated topographies were used to test the theoretical predictions.     

The role of aridification in constraining the elevation range of Holocene solifluction processes and associated landforms in the periglacial belt of the Sierra Nevada (southern Spain)

Climate variability during the Mid‐Late Holocene has influenced the activity of geomorphic processes in the current periglacial belt of the Sierra Nevada. We studied two types of sedimentary records that reveal a synchronous timing for slope instability in this high semi‐arid massif: solifluction landforms and mountain lake sediments. Lithological and sedimentological properties of both records have recorded numerous cycles of different magnitude of slope processes in the massif. Solifluction deposits record seven phases of solifluction activity and soil development during the last 7 ka bp and lake sediments show evidence of eight periods with increased geomorphic activity in the catchments over the last 6 ka bp . Although present‐day climate conditions do not promote active solifluction processes in the Sierra Nevada, colder and wetter periods during the Holocene triggered solifluction and transported coarse‐grained sediments into the lakes. By contrast, warm phases favoured soil formation and spread an incipient vegetation cover over the headwaters of the highest valleys, diminishing the grain size of the particles reaching the lakes. Lake sediments record an aridification trend in the massif intensifying since 4·2 ka bp that has conditioned solifluction activity to shift gradually to higher elevations. During major cooler phases such as the Little Ice Age active solifluction was recorded back down to 2500 m altitude. Copyright © 2010 John Wiley & Sons, Ltd.     

Variability in erosion rates related to the state of landscape transience in the semi‐arid Chilean Andes

We quantify erosion rates in the higher sectors of the Huasco Valley, in the Main Cordillera of the semi‐arid Andes of Chile, using elevation differences between three successive geomorphic markers (pediments and paleo‐valleys) and the present day valley. Available Ar‐Ar ages of Neogene pediments are used to estimate mean erosion rates for the three periods (16 to 13 My, 13 to 8 My, and following 8 My). The landscape of the Huasco Valley is in a transient state, as indicated by well‐preserved pediment surfaces in interfluves, valleys deeply incised by fluvial and glacial erosion and scarped head‐valleys that represent the current knickzones. Higher erosion rates (45–75 m/My) are calculated for the more recent period (< 8 My) during which deep incision developed compared to previous periods (6–31 m/My). Quantitative data indicate that glaciers had a much higher erosional capability than fluvial activity in the higher sectors of the Main Cordillera. Comparison with erosion rates calculated in other drainage basins of the Chilean Andes suggests that the variability of erosion rates depends on the landscape's transient erosive state. The landscape's geomorphologic response to the uplift of the Main Cordillera results in the retreat of a knickzone, for which retreat velocity depends on precipitation rate pattern and glacial erosion intensity. Copyright © 2011 John Wiley & Sons, Ltd.     

Evaluating restoration of man‐made slopes: a threshold approach balancing vegetation and rill erosion

The management of reclaimed slopes derived from industrial and civil activities (e.g. surface mining and road construction) requires the development of practical stability analysis approaches that integrate the processes and mechanisms that rule the dynamics of these ubiquitous emerging ecosystems. This work describes a new modelling approach focused on stability analysis of water‐limited reclaimed slopes, where interactive relationships between rill erosion and vegetation regulate ecosystem stability. Our framework reproduces two main groups of possible trends along the temporal evolution of reclaimed slopes: successful trends, characterized by widespread vegetation development and the effective control of rill erosion processes; and gullying trends, characterized by the progressive loss of vegetation and a sharp logistic increase in erosion rates. Furthermore, this analytical approach allows the determination of threshold values for the state variables (i.e. vegetation cover and rill erosion) that drive the system's stability, facilitating the identification of critical situations that require specific human intervention (e.g. revegetation or, in very problematic cases, revegetation combined with rill network destruction) to ensure the long‐term sustainability of the restored ecosystem. The application of our threshold analysis framework in Mediterranean‐dry reclaimed slopes derived from surface coal mining (the Teruel coalfield in central‐eastern Spain) showed a good field‐based performance. Therefore, we believe that this model is a valuable contribution for the management of water‐limited reclaimed systems, including those associated with rill erosion, as it provides a tool for the evaluation of restoration success and can play an important role in decision‐making during ecosystem restoration in severely disturbed landscapes. Copyright © 2011 John Wiley & Sons, Ltd.     

Impact of subsurface rock fragments on runoff and interrill soil loss from cultivated soils

Field and laboratory studies have indicated that rock fragments in the topsoil may have a large impact on soil properties, soil quality, hydraulic, hydrological and erosion processes. In most studies, the rock fragments investigated still remain visible at the soil surface and only properties of these visible rock fragments are used for predicting runoff and soil loss. However, there are indications that rock fragments completely incorporated in the topsoil could also significantly influence the percolation and water distribution in stony soils and therefore, also infiltration, runoff and soil loss rates. Therefore, in this study interrill laboratory experiments with simulated rainfall for 60 min were conducted to assess the influence of subsurface rock fragments incorporated in a disturbed silt loam soil at different depths below the soil surface (i.e. 0.001, 0.01, 0.05 and 0.10 m), on infiltration, surface runoff and interrill erosion processes for small and large rock fragment sizes (i.e. mean diameter 0.04 and 0.20 m, respectively). Although only small differences in infiltration rate and runoff volume are observed between the soil without rock fragments (control) and the one with subsurface rock fragments, considerable differences in total interrill soil loss are observed between the control treatment and both contrasting rock fragments sizes. This is explained by a rapid increase in soil moisture in the areas above the rock fragments and therefore a decrease in topsoil cohesion compared with the control soil profile. The observed differences in runoff volume and interrill soil loss between the control plots and those with subsurface rock fragments is largest after a cumulative rainfall (P_cum) of 11 mm and progressively decreases with increasing P_cum. The results highlight the impacts and complexity of subsurface rock fragments on the production of runoff volume and soil loss and requires their inclusion in process‐based runoff and erosion models. Copyright © 2011 John Wiley & Sons, Ltd.     

Response of the topside ionosphere over Arecibo to a moderate geomagnetic storm

We analyze the data obtained using the Arecibo incoherent scatter radar to examine the response of the topside ionosphere to a moderate geomagnetic storm that occurred during the period March 7–11, 2008. During this time period a magnetic storm with a non-monotonic main phase decrease in the Dst index occurred. The recovery phase also exhibited a secondary Dst decrease. During the initial phase of the storm, Te and Ti increased coincident with the arrival of the solar wind. The main phase registered an increase in proton concentration proportional to Ne while temperatures reached the lowest values. Variations in O⁺ concentration were not significant but a reduction in helium fraction was observed. Soon after the peak of the storm, the transition height between the topside ionosphere and the protonosphere, where H⁺ ions dominate composition, was lower than would be expected during quiet conditions and this behavior lasted for approximately 12 h.     

Adapting a texture synthesis algorithm for conditional multiple point geostatistical simulation

Computer vision provides several tools for analyzing and simulating textures. The principles of these techniques are similar to those in multiple-point geostatistics, namely, the reproduction of patterns and consistency in the results from a perceptual point of view, thus, ensuring the reproduction of long range connectivity. The only difference between these techniques and geostatistical simulation accounting for multiple-point statistics is that conditioning is not an issue in computer vision. We present a solution to the problem of conditioning simulated fields while simultaneously honoring multiple-point (pattern) statistics. The proposal is based on a texture synthesis algorithm where a fixed search (causal) pattern is used. Conditioning is achieved by adding a non-causal search neighborhood that modifies the conditional distribution from which the simulated category is drawn, depending on the conditioning information. Results show an excellent reproduction of the features from the training image, while respecting the conditioning information. Some issues related to the data structure and to the computer efficiency are discussed.