Complexities of landscape evolution during incision through layered stratigraphy with contrasts in rock strength

Variation in the erodibility of rock units has long been recognized as an important determinant of landscape evolution but has been little studied in landscape evolution models. We use a modified version of the Channel‐Hillslope Integrated Landscape Development (CHILD) model, which explicitly allows for variations in rock strength, to reveal and explore the remarkably rich, complex behavior induced by rock erodibility variations in even very simple geologic settings with invariant climate and tectonics. We study the importance of relative contrasts in erodibility between just two units, the order of these units (whether hard rocks overlie soft or soft rocks overlie hard) and the orientation of the contact between the two units. We emphasize the spatial and temporal evolution of erosion rates, which have important implications for basin analysis, detrital mineral records, and the interpretation of cosmogenic isotope concentrations in detrital samples. Results of the landscape evolution modeling indicate that the stratigraphic order of units in terms of erodibility, the gross orientation of the contact (i.e. dipping away or toward the outlet of the landscape) and the contact dip angle all have measurable effects on landscape evolution, including significant spatial and temporal variations in erosion rates. Steady‐state denudation conditions are unlikely to develop in landscapes with significant contrasts in rock strength in horizontal to moderately tilted rock layers, at least at the scale of the entire landscape. Additionally, our results demonstrate that there is no general relation between rock erodibility and erosion rates in natural settings. Although rock erodibility directly controls the erosion rate constant in our models, it is not uncommon for higher erosion rates to occur in the harder, less erodible rock. Indeed erosion rates may be either greater or less than the rock uplift rate (invariant in time and space in our models) in both hard and soft rocks, depending on the local geology, topography, and the pattern of landscape evolution. Copyright © 2016 John Wiley & Sons, Ltd.     

Response of transient rock uplift and base level knickpoints to erosional efficiency contrasts in bedrock streams

Knickpoints in bedrock streams are often interpreted as transient features generated by a change in boundary conditions. It is typically assumed that knickpoints propagate upstream with constant vertical velocities, though this relies on a stream being in erosional steady state (erosion rate equals rock uplift rate) prior to the knickpoint's formation. Recent modeling and field studies suggest that along-stream contrasts in rock erodibility perturb streams from erosional steady state. To evaluate how contrasts in rock erodibility might impact knickpoint interpretations, we test parameter space (rock erodibility, rock contact dip angle, change in rock uplift rate) in a one-dimensional (1D) bedrock stream model that has variable rock erodibility and produces a knickpoint with a sudden change in rock uplift rate. Upstream of a rock contact, the vertical velocity of a knickpoint generated by a change in rock uplift rate is strongly correlated with how the rock contact has previously perturbed erosion rates. These knickpoints increase vertical velocity upon propagating upstream of a hard over soft contact and decrease vertical velocity upon propagating upstream of a soft over hard contact. However, interactions with other transient perturbations produced by rock contacts make for nuances in knickpoint behavior. Rock contacts also influence the geometry of knickpoints, which can become particularly difficult to identify upstream of soft over hard rock contacts. Using our simulations, we demonstrate how a contact's along-stream horizontal migration rate and cross-contact change in rock strength control how much an upstream reach is perturbed from erosional steady state. When simulations include multiple contacts, the knickpoint is particularly prone to colliding with other transient perturbations and can even disappear altogether if rock contact dips are sufficiently shallow. Caution should be taken when analyzing stream profiles in areas with significant changes in rock strength, especially when rock contact dip angles are near the stream's slope.     

Analysis of a database for anthropogenic sinkhole triggering and zonation in the Naples hinterland (Southern Italy)

Anthropogenic sinkholes are generally due to the collapse of man-made underground caves and represent a major threat, especially in urban contexts, where they could cause damage to people, buildings and lifelines. The hinterland of Naples (Campania, Southern Italy) is one of the most urbanized areas of Italy and is historically affected by frequent sinkhole phenomena. This study analyzes a database of both man-made caves and sinkholes collected by the authors over several years. The aim is to comprehend the predisposing and triggering factors of sinkholes in order to attempt a zonation of occurrence. The predisposing factor resulted to be the presence of a preexisting network of caves within the Campanian Ignimbrite tuff and their geometric features. Generally, the shallower the depth of the chamber roof and the lower the thickness of the vault, the higher is the frequency of sinkholes. Furthermore, an intrinsic fragility is represented by the access shafts of vertical wells, usually filled in and abandoned. Meanwhile, the main triggering mechanism identified consists in saturation of the subsoil, due to water leaks coming from buried sewage and water pipelines. The macrozonation of sinkhole occurrence shows that the highest class is achieved where the tuff is shallow and both man-made caves and historical sinkholes are present.

     

Considerations concerning the non-rigid Earth nutation theory

This paper presents the reflections of the Working Group of which the tasks were to examine the non-rigid Earth nutation theory. To this aim, six different levels have been identified: Level 1 concerns the input model (giving profiles of the Earth's density and theological properties) for the calculation of the Earth's transfer function of Level 2; Level 2 concerns the integration inside the Earth in order to obtain the Earth's transfer function for the nutations at different frequencies; Level 3 concerns the rigid Earth nutations; Level 4 examines the convolution (products in the frequency domain) between the Earth's nutation transfer function obtained in Level 2, and the rigid Earth nutation (obtained in Level 3). This is for an Earth without ocean and atmosphere; Level 5 concerns the effects of the atmosphere and the oceans on the precession, obliquity rate, and nutations; Level 6 concerns the comparison with the VLBI observations, of the theoretical results obtained in Level 4, corrected for the effects obtained in Level 5.Each level is discussed at the state of the art of the developments.