Results from a new model of river basin evolution

This paper briefly describes a model of the erosional development of catchments and their channel networks. The model differentiates between the dominant transport processes in hillslope and channels. The development of channels and hillslopes occurs in an integrated manner as a function of physically observable mechanisms. The growth of a river basin is qualitatively described. The model concepts are used to study the basin during periods of growth (transient periods), as well as during dynamic equilibrium. This leads to hypotheses about the relationship between slopes, relief, tectonic uplift, erodability, runoff, and catchment area. It is shown that the model leads to very reasonable and desirable behaviour of hillslopes during retreat and degradation.     

Model prediction capacity of ephemeral gully evolution in conservation tillage systems

Ephemeral gully (EG) erosion has an important impact on agricultural soil losses and increases field surface hydrology connectivity and transport of pollutants to nearby water bodies. Watershed models including an EG component are scarce and not yet properly evaluated. The objective of this study is to evaluate the capacity of one such tool, AnnAGNPS, to simulate the evolution of two EG formed in a conservation tillage system. The dataset for model testing included runoff measurements and EG morphological characteristics during 3 years. Model evaluation focused on EG evolution of volume, width, and length model outputs, and included calibration and testing phases and a global sensitivity analysis (GSA). While the model did not fully reproduce width and length, the model efficiency to simulate EG volume was satisfactory for both calibration and testing phases, supporting the watershed management objectives of the model. GSA revealed that the most sensitive factors were EG depth, critical shear stress, headcut detachment exponent coefficient b, and headcut detachment leading coefficient a. For EG outputs the model was additive, showing low sensitivity to interactions between the inputs. Prediction of EG spatial evolution on conservation tillage systems requires improved development of gully erosion components, since many of the processes were developed originally for traditional tillage practices or larger channel systems. Our results identify the need for future research when EG form within conservation tillage systems, in particular to study gully headcut, soil erodibility, and width functions specific to these practices.     

Modeling flocculation in a hypertidal estuary

When fine particles are involved, cohesive properties of sediment can result in flocculation and significantly complicate sediment process studies. We combine data from field observations and state-of-the-art modeling to investigate and predict flocculation processes within a hypertidal estuary. The study site is the Welsh Channel located at the entrance of the Dee Estuary in Liverpool Bay. Field data consist of measurements from a fixed site deployment during 12–22 February 2008. Grain size, suspended sediment volume concentration, and current velocity were obtained hourly from moored instruments at 1.5 m above bed. Near-bottom water samples taken every hour from a research vessel are used to convert volume concentrations to mass concentrations for the moored measurements. We use the hydrodynamic model Proudman Oceanographic Laboratory Coastal Ocean Modelling System (POLCOMS) coupled with the turbulence model General Ocean Turbulence Model (GOTM) and a sediment module to obtain three-dimensional distributions of suspended particulate matter (SPM). Flocculation is identified by changes in grain size. Small flocs were found during flood and ebb periods—and correlate with strong currents—due to breakup, while coarse flocs were present during slack waters because of aggregation. A fractal number of 2.4 is found for the study site. Turbulent stresses and particle settling velocities are estimated and are found to be related via an exponential function. The result is a simple semiempirical formulation for the fall velocity of the particles solely depending on turbulent stresses. The formula is implemented in the full three-dimensional model to represent changes in particle size due to flocculation processes. Predictions from the model are in agreement with observations for both settling velocity and SPM. The SPM fortnight variability was reproduced by the model and the concentration peaks are almost in phase with those from field data.     

An audio-magnetotelluric investigation in Terceira Island (Azores)

Ten audio-magnetotelluric soundings have been carried out along a profile crossing the Serra do Cume caldera in the eastern part of the Terceira Island (Azores). The main objectives of this investigation were to detect geoelectrical features related with tectonic structures and to characterize regional hydrological and hydrothermal aspects mainly those related to geothermal fluid dynamics.Three-dimensional numerical investigation showed that the data acquired at periods shorter than 1 s are not significantly affected by ocean effect. The data was analysed using the Smith's decomposition method in order to investigate possible distortions caused by superficial structures and to estimate a global regional strike. The results suggest that in general the soundings were not distorted. A regional N55°W strike was chosen for the two-dimensional data inversion.The low-resistivity zones (10–30 ohm-m) displayed in the central part of the 2-D geoelectrical model have been interpreted as caused by hydrothermal circulation. The low-resistivity anomalies at the ends of the profile might be attributed to alteration zones with interaction of seawater intrusion. High-resistivity (> 300 ohm-m) values have been related with less permeable zones in the SW of Cinco Picos and Guilherme Moniz caldera walls.     

Interactive effects of environmental variability and human impacts on the long-term dynamics of an Amazonian floodplain lake and a South Atlantic coastal lagoon

Human activities are exposing freshwater ecosystems to a wide range of stressors, whose direct and indirect effects can be alleviated or exacerbated through interactive effects with dynamic environmental drivers. This study used long-term data from two Neotropical lacustrine freshwater systems (Batata Lake, an Amazonian floodplain lake and Imboassica lagoon, an Atlantic coastal lagoon) subjected to different kinds of environmental fluctuations (i.e., flood pulse and sandbar opening) and anthropogenic impacts (i.e., siltation and eutrophication). Our objective was to determine whether the effects of human perturbations are contingent on modifications of important biotic and abiotic characteristics through environmental variability. For both ecosystems, environmental variability consistently interacted with anthropogenic perturbations to alter most of the variables analyzed, such as nutrient dynamics, chlorophyll-a concentration, zooplankton and benthic invertebrate species richness, and temporal community stability, which indicates that interactive effects between environmental variability and anthropogenic perturbations may impact a myriad of ecosystem properties. Furthermore, the nature of these interactive effects was highly dependent on the variable considered and on the ecosystem analyzed. For example, at Imboassica lagoon, sandbar openings interacted synergistically with trophic state to increase the phosphorus concentration in the water column. At Batata Lake, flooding generally alleviated the negative effects of siltation on species richness by both diluting inorganic suspended material concentration and by promoting local recruitment from the regional species pool. Such results indicate that our ability to understand and predict the outcome of anthropogenic impacts on inland aquatic systems can be hampered if we consider human stressors as “static” phenomena disconnected from dynamic interactions with major local environmental drivers.     

A quantitative evaluation of Playfair's law and its use in testing long‐term stream erosion models

Playfair's law (J. Playfair, illustrations of the Huttonian Theory of the Earth, 1802) requires any two tributaries in a river network to lower at the same rate near their junction. Although this law holds exactly at the junction, it is unclear how well it holds in the vicinity of the junction. This issue has practical importance because Playfair's law has been used to estimate parameters for detachment‐limited models of erosion. If the incision rate of a stream is modelled as βA^mSⁿ, where β is an erodibility parameter, A is the area drained by the stream, and S is the local gradient of the channel, then the ratio of the parameters m/n can be estimated from junctions by assuming that Playfair's law holds over the distance used to determine S for each tributary. In this paper, Playfair's law and associated m/n estimates are evaluated for simulated basins with constant and temporally varying uplift rates (or baselevel lowering rates). The results demonstrate that estimates of m/n may be biased for basins with upward‐concave stream profiles because the local slope must be approximated with an average upstream slope. In addition, when uplift rate varies temporally, knickpoints are shown to travel through the basins with constant vertical velocity. Because incision rates vary within the basin, Playfair's law only holds exactly at the junctions. These effects are more important when slopes are measured over longer distances. Finally, measurement techniques are presented which address these potential biases. Copyright © 2001 John Wiley & Sons, Ltd.     

Multiple scattering ofSH waves in 2-D media with many cavities

The full waveform synthetic seismogram of multiple scatteredSH waves by many cylindrical cavities in two-dimensional homogeneous elastic media is computed. We used the so-called single-layer potential integral representation of the scattered field and a discretization scheme with line source distribution for each cavity. The total field is the sum of the incident wave plus the field radiated from all sources, each multiplied by an unknown complex constant representing its strength. These constants are determined by imposing the appropriate boundary conditions in the least-squares sense. Here we solve scattering problems involving one, two, four, twelve and fifty cavities regularly distributed in a half-space. The seismograms computed along the free-surface show regions where the incident wave is strongly attenuated, as well as the arrivals of all multiple scattered phases. The accuray of the method is estimated from the degree of agreement of our solution for one cavity with the corresponding analytical solution, and also from the magnitude of the residual tractions along the boundaries of two cavities separated at various distances. Finally we apply the method to compute the case of fifty cylindrical cavities, each of radiusa, randomly distributed in a region 80a wide by 30a deep in a half-space. The value of scattering loss is obtained from the amplitude decay of the primary wave with distance for wavelengths in the range from 1.7a to 13.3a, using the synthetic seismogram calculated for the same distribution of 50 cavities as above, but in full-space.     

Early-Morning Flow Transition in a Valley in Low-Mountain Terrain Under Clear-Sky Conditions

We investigate the evolution of the early-morning boundary layer in a low-mountain valley in south-western Germany during COPS (convective and orographically induced precipitation study) in summer 2007. The term low-mountain refers to a mountainous region with a relief of gentle slopes and with an absolute altitude that remains under a specified height (usually 1,500 m a.s.l.). A subset of 23 fair weather days from the campaign was selected to study the transition of the boundary-layer flow in the early morning. The typical valley atmosphere in the morning hours was characterized by a stable temperature stratification and a pronounced valley wind system. During the reversal period—called the low wind period—of the valley wind system (duration of 1–2 h), the horizontal flow was very weak and the conditions for free convection were fulfilled close to the ground. Ground-based sodar observations of the vertical wind show enhanced values of upward motion, and the corresponding statistical properties differ from those observed under windless convective conditions over flat terrain. Large-eddy simulations of the boundary-layer transition in the valley were conducted, and statistical properties of the simulated flow agree with the observed quantities. Spatially coherent turbulence structures are present in the temporal as well as in the ensemble mean analysis. Thus, the complex orography induces coherent convective structures at predictable, specific locations during the early-morning low wind situations. These coherent updrafts, found in both the sodar observations and the simulation, lead to a flux counter to the gradient of the stably stratified valley atmosphere and reach up to the heights of the surrounding ridges. Furthermore, the energy balance in the surface layer during the low wind periods is closed. However, it becomes unclosed after the onset of the valley wind. The partition into the sensible and the latent heat fluxes indicates that missing flux components of sensible heat are the main reason for the unclosed energy balance in the considered situations. This result supports previously published investigations on the energy balance closure.     

Assessment of hydrological alterations on wandering beetle assemblages (coleoptera: Carabidae and Tenebrionidae) in coastal wetlands of arid Mediterranean systems

Changes in wandering beetle assemblages (Carabidae and Tenebrionidae) of different habitats situated in coastal wetlands of a Mediterranean arid area (Mar Menor, SE Spain) were analysed in 1984, 1992 and 2003 by pitfall trapping. Over two decades, the increase in irrigated lands at watershed scale led to rising water tables in the Mar Menor wetlands, which affected their beetle communities. These hydrological changes caused an increase in the carabid population, particularly in the sites most affected by flooding, where halobionts and halophiles, which were practically absent in 1984 and 1992, had become dominant by 2003. In contrast, tenebrionid assemblages simplified with time and by 2003 were dominated by one or two generalist species.     

Shallow necking depth and differential denudation linked to post‐rift continental reactivation: The origin of the Cenozoic basins in southeastern Brazil

The southeastern Brazilian margin presents post‐breakup Cenozoic tectonism that created a series of grabens and small sedimentary basins, known as the Continental Rift of Southeastern Brazil. The formation of this rift occurred long after the South Atlantic ocean opening and has been attributed to different mechanisms like regional uplift induced by hotspot activity, pulses of Andean orogeny and reactivation of pre‐existing faults. However, the proposed models lack an analytical or numerical verification from a geodynamic point of view. Based on finite element modelling of the lithospheric stress field evolution we conclude that a shallow necking depth, consistent with the hyperextended southeastern Brazilian margin, combined with differential denudation of the continent, resulted in an extensional stress field in the upper crust that induced the observed Cenozoic tectonism.