Modelling net erosion responses to enviroclimatic changes recorded upon multisecular timescales

Land use change has been recognized throughout the Earth as one of the most important factors influencing the occurrence of rainfall-driven geomorphological processes. However, relating the occurrence of historical soil erosion rates is difficult because of the lack of long-term research projects in river basins. Also, complex models are not adequate to reconstruct erosion rate changes because they require significant input data not always available on long timescales. Given the problems with assessing sediment yield using complex erosion models, the objective of this study is to explore a parsimonious scale-adapted erosion model (ADT) from the original Thornes and Douglas algorithms, which aims at reconstruction of annual net erosion (ANE) upon multisecular timescales. As a test site, the Calore River basin (3015 km² in southern Italy) provides a peculiar and unique opportunity for modelling erosion responses to climate and land cover changes, where input-data generation and interpretation results were also supported by documented hydrogeomorphological events that occurred before and after land deforestation. In this way, ANE_ADT-values were reconstructed for the period 1675–2004 by using precipitation indexes, complemented by recent instrumental records, and by using land cover statistics from documented agrarian sources. Pulses of natural sedimentation in the predeforestation period have been related to Vesuvius volcanic activity and changes in rainstorm frequency. After deforestation, the basin system became unstable with sudden fluctuations in the hydrogeomorphological regime contributing significantly to increased erosion and, in turn, sediment transport sequences via river drainage towards the Tyrrhenian coast.     

Structural and nutritional heterogeneity of riparian vegetation in Patagonia (Argentina) in relation to seasonal grazing by sheep

Riparian meadows of southern Patagonia are temporally and spatially heterogeneous habitats. They are thought to play a key role in regulating the dynamics of arid grazed steppes of the region. We conducted a 2-year study with sheep to determine the influence of two grazing conditions (deferment and spring grazing) on structural and nutritional parameters of three vegetation types apparently associated with a soil moisture gradient in a riparian meadow in southern Santa Cruz (Argentina). Spring deferment allowed forage accumulation in very moist and intermediate vegetation types by the beginning of the summer and had no detrimental impact on forage quality. Deferment had no effect on biomass accumulation of the drier vegetation type at our study site. Structural attributes of the vegetation appeared to affect spatial grazing patterns of sheep. Soil-related patchiness, rather than sheep grazing, appeared to control vegetation structure and nutritional value. The proportion of more mesic plant communities in riparian habitats of Patagonia could determine the outcome of plant–animal interactions in these systems. Practical recommendations of grazing capacity in paddocks, or more theoretical considerations of ecosystem dynamics of the Patagonian steppes need to explicitly consider the internal patchiness of riparian habitats.     

Crustal structure across the Grand Banks–Newfoundland Basin Continental Margin – II. Results from a seismic reflection profile

New multichannel seismic reflection data were collected over a 565 km transect covering the non-volcanic rifted margin of the central eastern Grand Banks and the Newfoundland Basin in the northwestern Atlantic. Three major crustal zones are interpreted from west to east over the seaward 350 km of the profile: (1) continental crust; (2) transitional basement and (3) oceanic crust. Continental crust thins over a wide zone (∼160 km) by forming a large rift basin (Carson Basin) and seaward fault block, together with a series of smaller fault blocks eastwards beneath the Salar and Newfoundland basins. Analysis of selected previous reflection profiles (Lithoprobe 85-4, 85-2 and Conrad NB-1) indicates that prominent landward-dipping reflections observed under the continental slope are a regional phenomenon. They define the landward edge of a deep serpentinized mantle layer, which underlies both extended continental crust and transitional basement. The 80-km-wide transitional basement is defined landwards by a basement high that may consist of serpentinized peridotite and seawards by a pair of basement highs of unknown crustal origin. Flat and unreflective transitional basement most likely is exhumed, serpentinized mantle, although our results do not exclude the possibility of anomalously thinned oceanic crust. A Moho reflection below interpreted oceanic crust is first observed landwards of magnetic anomaly M4, 230 km from the shelf break. Extrapolation of ages from chron M0 to the edge of interpreted oceanic crust suggests that the onset of seafloor spreading was ∼138 Ma (Valanginian) in the south (southern Newfoundland Basin) to ∼125 Ma (Barremian–Aptian boundary) in the north (Flemish Cap), comparable to those proposed for the conjugate margins.     

Interpretation of subhorizontal crustal reflections by metamorphic and rheologic effects in the eastern part of the Pannonian Basin

The geologic origin of subhorizontal reflections, often observed in crustal seismic sections, was investigated by establishing metamorphic facies and strength of rocks in depth, and correlating these properties to seismic reflection sections from eastern Hungary. Estimation of the depths of metamorphic mineral stability zones utilized the principles developed by Fyfe et al. and known geothermal data of the area. The strength versus depth profile was derived by relating local seismic P -wave interval velocities to Meissner et al. 's activation energy. The results show that the series of subhorizontal reflections, observed in the Pannonian Basin, are a consequence of combined metamorphic and rheologic changes in depths. The synthesis of the integrated data set suggests that the retrograde alteration of the pre-Tertiary basement above the percolation threshold was made possible by the softening effect of shear zones and their water-conducting capacity. The subhorizontal reflections of highest energy, of the consolidated crust below the percolation threshold, originate in the depths of greenschist, amphibolite and granulite metamorphic mineral facies, which were formed in geothermal and pressure conditions similar to those existing today. These results imply the overprint of earlier (Variscan) metamorphic sequences of the crust by more recent retrograde metamorphic processes.     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.