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.     

Modelling considerations in probabilistic performance‐based seismic evaluation: case study of the I‐880 viaduct

Limitations associated with deterministic methods to quantify demands and develop rational acceptance criteria have led to the emergence of probabilistic procedures in performance‐based seismic engineering. The Pacific Earthquake Engineering Research performance‐based methodology is one such approach. In this paper, the impact of certain modelling decisions made at different stages of the evaluation process on the performance assessment of a typical multi‐bent viaduct is examined. Modelling, in the context of this paper, covers hazard modelling, structural modelling and loss modelling. The specific application considered in this study is a section of an existing viaduct in California: the I‐880 interstate highway. Several simulation models of the viaduct are developed, a series of nonlinear time‐history analyses are carried out to predict demands, measures of damage are evaluated and the probability of closure of the viaduct is estimated using the specified hazard for the site. It is concluded that the methodology offers several advantages over existing deterministic performance‐based procedures. Results of the investigation indicate that the assessment methodology is particularly sensitive to the reliability of decisions made by bridge inspectors following a seismic event, and to the dispersion in the demand estimation, which in turn is influenced by several factors including soil–structure interaction effects and ground motion scaling procedures. Copyright © 2005 John Wiley & Sons, Ltd.     

On the effect of the rigid sidewall on the dynamic stiffness of embedded circular footings

The effect of the rigid sidewall, which is usually combined with embedded footings, on the dynamic stiffness of the footings is considered. An efficient numerical technique is used to calculate the static and dynamic stiffness of circular footings embedded in a stratum. The results show that the increase in static stiffness with increasing height of the sidewall is most significant in the case of rocking. The dynamic stiffness coefficients change considerably, if the sidewall extends higher than about half the depth of embedment. The damping coefficients corresponding to vertical vibrations and rocking are likewise affected by the height of the sidewall. The damping coefficients corresponding to torsional and horizontal vibrations increase considerably with increasing height of the sidewall.     

Geotechnical characterisation of pyroclastic soils involved in huge flowslides

Landslides of the flow type involving granular geo-materials frequently result in casualties and damage to property because of the long travel distance and the high velocities that these may attain. This was true for the events that took place in Campania Region (Southern Italy) in May 1998, involving pyroclastic soils originating from explosive activities of the Somma-Vesuvius volcano. Although these phenomena have frequently affected various areas of the Campania region over the last few centuries, there were no useful geological and geotechnical references available in the aftermath of the May 1998 events. For this reason Salerno University, which was involved in the scientific management of the emergency, addressed the issue of acquiring data on the geological, geomorphological and hydrogeological features of the slopes where the landslides had taken place. The information acquired made it possible to set up a slope evolution model that is able to interpret, from a geological point of view, past and more recent landslides that had occurred in the same area. As preliminary geotechnical analyses had already validated the above model, more detailed investigations were performed both on the pore pressure regimen of the covers still in place as well as on the physical and mechanical properties of pyroclastic soils, in saturated and unsaturated conditions. The present paper begins by discussing the data acquired during the .rst phase of the studies and then goes on to illustrate the laboratory results so far obtained with the aid of approximate procedures. These help advance our knowledge of pyroclastic soils within a reasonable time frame, thus improving landslide triggering analysis.     

Identification and isotopic studies of early Cambrian magmatism (El Carancho Igneous Complex) at the boundary between Pampia terrane and the Río de la Plata craton,La Pampa province,Argentina

We have identified late Early Cambrian metaigneous rocks very poorly exposed at the Estancia El Carancho, in central La Pampa province, Argentina. They comprise calc–alkaline metadiorite and metagranite, and tholeiitic metapyroxenite and metagabbro. They are jointly referred to as the El Carancho Igneous Complex, and regarded to pertain to the Pampean magmatic arc and backarc, respectively. Titanite U–Pb SHRIMP dating of the metapyroxenite yielded 528 ± 5 Ma, and zircon U–Pb SHRIMP dating of the metadiorite yielded 520 ± 1.4 Ma. Hafnium isotope determinations on the dated zircons show ¹⁷⁶Hf/¹⁷⁷Hf ratios corresponding to positive εHf values from + 7.18 to + 9.37; Hf model ages of the Cambrian zircons yielded 884 Ma. It is interpreted that the metadiorites of the Complex crystallized from an Early Neoproterozoic (Tonian) juvenile source. We argue that the inferred occurrence of juvenile Tonian magmatic rocks in the (otherwise, mostly Paleo-Mesoproterozoic) substratum of the southern Pampia terrane could indicate a zone of thinned basement possibly associated with the early stage of Rodinia's breakup. In addition, the studied segment of the Pampean magmatic arc is contaminated by also juvenile, Late Mesoproterozoic crust, as evidenced by the presence of xenocrystic cores of 1140–1194 Ma – T_DM-Hf 1720 Ma and εHf values of + 3.24 to + 4.85 – in the Cambrian zircons, hence suggesting that the studied segment of the Pampean magmatic arc was intruded into juvenile Late Mesoproterozoic magmatic arc rocks. The El Carancho Igneous Complex would be located at the tectonic boundary between the Pampia terrane and the Río de la Plata craton. This boundary stands out in the aeromagnetic data as a change in the structural orientation about a roughly N-S line located approximately at 65° W and representing the suture zone between the Pampia terrane and the Río de la Plata craton. Our geotectonic model envisages westward dipping subduction of oceanic crust beneath the Pampia terrane; the El Carancho Igneous Complex would, therefore, have been originated on the Pampia side (upper plate) of the suture. Slivers of the arc- and backarc-type rocks would have been tectonically imbricated in the suture zone during the Pampean orogeny.     

Tectono-sequence stratigraphy and U–Pb zircon ages of the Rincón Blanco Depocenter,northern Cuyo Rift,Argentina

Extensional processes that followed the Gondwanan Orogeny rise to the development of a series of rift basins along the continental margin over older accreted Eopaleozoic terranes. Stratigraphic, structural, paleontological, and isotopic studies are presented in this work in order to constrain the ages of the sedimentary infilling and to analyze the tectosedimentary evolution of one of the Cuyo basin depocenters, known as Rincón Blanco. This asymmetrical half-graben was filled by continental sediments under a strong tectonic control. The infilling was strongly controlled by tectonics which in term produced distinctive features along the whole sedimentary sequence. Using a combination of lithological and structural data the infilling was subdivided into packages of genetically linked units bounded by regional extended surfaces. Several tuffs and acid volcanic rocks have been collected across the whole section of the Rincon Blanco sub-basin for SHRIMP and LA-MC-ICPMS U–Pb zircon dating. The ages obtained range from 246.4 ± 1.1 Ma to 230.3 ± 1.5 Ma which is the time elapsed for the deposition of three tectono-sequence units separated by regional unconformities and mainly constrained to the Middle Triassic. They are interpreted as a result of a reactivation of the extensional system that has evolved along strike as segments of faults that linked together and/or as laterally propagating faults. Regional correlation with coeval rift basins permits to establish north-south propagation in the extensional regime along the western margin of SW Gondwana. This trending started in the lowermost Triassic and extended until the latest Triassic. Two of them were precisely correlated with Cerro Puntudo and Cacheuta half-graben systems. The new data indicate that the three sequences were mostly deposited during the Middle Triassic (246 to 230 Ma), with no evidence of sedimentation during Norian and Rhaetian, which is in conflict with some previous biostratigraphic studies.     

Modelling size effect on rock aggregates strength using a DEM bonded-cell model

Empirical evidence has shown that particle breakage affects the mechanical behaviour of granular materials. The source of this mechanism takes place at the particle scale, and the main consequence on the macromechanical behaviour is increasing compressibility. Due to the inverse correlation between particle size and particle crushing strength, coarse rockfill materials are particularly vulnerable to mechanical degradation due to particle breakage. However, such coarse materials do not fit in standard laboratory devices, and the alternative of large sample testing is usually unavailable or too expensive. Alternatively, recent works have proposed multi-scale approaches using the discrete element method (DEM) to carry out numerical testing of coarse crushable materials, although few studies have focused on size effects. This article presents the application of a DEM bonded-cell model to study particle size-strength correlation on angular rock aggregates. Each particle is modelled by a cluster of perfectly rigid polyhedral cells with Mohr–Coulomb contact law. Constant cell density within particles implies that the presence of potential fragmentation planes increases with size. Therefore, particle strength decreases with size. A comprehensive sensitivity analysis was carried out through 1477 particle crushing simulations in a given particle size. Based on published experimental data on calcareous rock aggregates, part of the simulations were used for calibration, and 97 additional simulations of a coarser size fraction were performed for validation. The results show a good agreement with the empirical data in terms of size effect and data scatter through Weibull statistics.

     

The stability of a remediated bed in Hamilton Harbour,Lake Ontario,Canada

ABSTRACT In situ measurements of lakebed sediment erodibility were made on three sites in Hamilton Harbour, Lake Ontario, using the benthic flume Sea Carousel. Three methods of estimating the surface erosion threshold (τ_c(0)) from a Carousel time series were evaluated: the first method fits measures of bed strength to eroded depth (the failure envelope) and evaluates threshold as the surface intercept; the second method regresses mean erosion rate (E_m) with bed shear stress and solves for the floc erosion rate (E_f) to derive the threshold for E_m = E_f = 1 × 10^?5 kg m^?2 s^?1; the third method extrapolates a regression of suspended sediment concentration (S) and fluid transmitted bed shear stress (τ₀) to ambient concentrations. The first field site was undisturbed (C) and acted as a control; the second (W) was disturbed through ploughing and water injection as part of lakebed treatment, whereas the third site (OIP) was disturbed and injected with an oxidant used for remediation of contaminated sediment. The main objectives of this study were: (1) to evaluate the three different methods of deriving erosion threshold; (2) to compare the physical behaviour of lacustrine sediments with their marine estuarine counterparts; and (3) to examine the effects of ploughing and chemical treatment of contaminated sediment on bed stability. Five deployments of Sea Carousel were carried out at the control site. Mean erosion thresholds for the three methods were: τ_c(0) = 0·5 (±0·06), 0·27 (±0·01) and 0·34 (±0·03) Pa respectively. Method 1 overpredicted bed strength as it was insensitive to effects in the surface 1–2 mm, and the fit of the failure envelope was also highly subjective. Method 2 exhibited a wide scatter in the data (low correlation coefficients), and definition of the baseline erosion rate (E_f) is largely arbitrary in the literature. Method 3 yielded stable (high correlation coefficients), reproducible and objective results and is thus recommended for evaluation of the erosion threshold. The results of this method correlated well with sediment bulk density and followed the same trend as marine counterparts from widely varying sites. Mass settling rates, expressed as a decay constant, k, of S(t), were strongly related to the maximum turbidity at the onset of settling (S_max) and were also in continuity with marine counterparts. Thus, it appears that differences in salinity had little effect on mass settling rates in the examples presented, and that biological activity dominated any effects normally attributable to changes in salinity. Bedload transport of eroded aggregates (2–4 mm in diameter) took place by rolling below a mean tangential flow velocity (U_y) of 0·32 ms^?1 and by saltation at higher velocities. Mass transport as bedload was a maximum at U_y = 0·4 ms^?1, although bedload never exceeded 1% of the suspended load. The proportion of material moving as bedload was greatest at the onset of erosion but decreased as flow competence increased. Given the low bulk density and strength of the lakebed sediment, the presence of a bedload component is notable. Bedload transport over eroding cohesive substrates should be greater in estuaries, where both sediment density and strength are usually higher. Significant differences between the ploughed and control sites were apparent in both the erosion rate and the friction coefficient (φ), and suggest that bed recovery after disruption is rapid (< 24 h). τ_c(0) increased linearly with time after ploughing and recovered to the control mean value within 3 days. The friction coefficient was reduced to zero by ploughing (diagnostic of fluidization), but increased linearly with time, regaining control values within 6 days. No long‐term reduction in bed strength due to remediation was apparent.