Shaking‐table tests of flat‐bottom circular silos containing grain‐like material

According to Eurocode 8, the seismic design of flat‐bottom circular silos containing grain‐like material is based on a rough estimate of the inertial force imposed on the structure by the ensiled content during an earthquake: 80% of the mass of the content multiplied by the peak ground acceleration. A recent analytical consideration of the horizontal shear force mobilised within the ensiled material during an earthquake proposed by some of the authors has resulted in a radically reduced estimate of this load suggesting that, in practice, the effective mass of the content is significantly less than that specified. This paper describes a series of laboratory tests that featured shaking table and a silo model, which were conducted in order to obtain some experimental data to verify the proposed theoretical formulations and to compare with the established code provisions. Several tests have been performed with different heights of ensiled material – about 0.5 mm diameter Ballotini glass – and different magnitudes of grain–wall friction. The results indicate that in all cases, the effective mass is indeed lower than the Eurocode specification, suggesting that the specification is overly conservative, and that the wall–grain friction coefficient strongly affects the overturning moment at the silo base. At peak ground accelerations up to around 0.35 g, the proposed analytical formulation provides an improved estimate of the inertial force imposed on such structures by their contents. Copyright © 2015 John Wiley & Sons, Ltd.     

Sediment and particulate phosphorus characteristics in grassed waterways from row crop corn and alfalfa fields collected by manual University of Exeter samplers and automatic sampling

Phosphorus (P) export from agricultural lands above known threshold levels can result in adverse impacts to receiving water quality. Phosphorus loss occurs in dissolved and sediment‐bound, or particulate phosphorous (PP), forms, with the latter often dominating losses from row‐cropped systems. To target practices, land managers need good computer models and model developers need good monitoring data. Sediment monitoring data (e.g. radiometric finger printing and sediment P sorption capacity) can help identify sediment source areas and improve models, but require more sediment mass than is typically obtained by automatic sampling. This study compares a simple suspended sediment sampler developed at the University of Exeter (UE) with automatic sampling in intermittent channels draining corn and alfalfa fields. The corn field had a greater runoff coefficient (27%) than alfalfa (11%). No differences were found in enrichment ratios (sediment constituent/soil constituent) in PP (PP_ER) or percent loss on ignition (LOI_ER) between paired UE samplers on corn. The median LOI_ER for the UE samplers (1·9%) did not differ significantly (p > 0·13) from the automatic sampler (2·0%). The PP_ER from the UE samplers was on average 20% lower than the automatic samplers. A correlation (r² = 0·75) was found between sediment PP and % LOI from automatic samplers and UE samplers for particles < 50 µm, while for > 50 µm PP concentration did not change with changes in % LOI. Sediment ammonium‐oxalate extractable metals were similarly related to LOI, with the strongest correlation for iron (r² = 0·71) and magnesium (r² = 0·70). Copyright © 2011 John Wiley & Sons, Ltd.     

Representing Grounding Line Dynamics in Numerical Ice Sheet Models: Recent Advances and Outlook

Recent satellite observations of the Antarctic and Greenland ice sheets show accelerated ice flow and associated ice sheet thinning along coastal outlet glaciers in contact with the ocean. Both processes are the result of grounding line retreat due to melting at the grounding line (the grounding line is the contact of the ice sheet with the ocean, where it starts to float and forms an ice shelf or ice tongue). Such rapid ice loss is not yet included in large-scale ice sheet models used for IPCC projections, as most of the complex processes are poorly understood. Here we report on the state-of-the art of grounding line migration in marine ice sheets and address different ways in which grounding line migration can be attributed and represented in ice sheet models. Using one-dimensional ice flow models of the ice sheet/ice shelf system we carried out a number of sensitivity experiments with different spatial resolutions and stress approximations. These are verified with semi-analytical steady state solutions. Results show that, in large-scale finite-difference models, grounding line migration is dependent on the numerical treatment (e.g. staggered/non-staggered grid) and the level of physics involved (e.g. shallow-ice/shallow-shelf approximation).     

BET_VH: exploring the influence of natural uncertainties on long-term hazard from tephra fallout at Campi Flegrei (Italy)

In this paper, we explore the effects of the intrinsic uncertainties upon long-term volcanic hazard by analyzing tephra fall hazard at Campi Flegrei, Italy, using the BET_VH model described in Marzocchi et al. (Bull Volcanol, 2010). The results obtained show that volcanic hazard based on the weighted average of all possible eruptive settings (i.e. size classes and vent locations) is significantly different from an analysis based on a single reference setting, as commonly used in volcanic hazard practice. The long-term hazard map for tephra fall at Campi Flegrei obtained here accounts for a wide spectrum of uncertainties which are usually neglected, largely reducing the bias intrinsically introduced by the choice of a specific reference setting. We formally develop and apply a general method to recursively integrate simulations from different models which have different characteristics in terms of spatial coverage, resolution and physical details. This outcome of simulations will be eventually merged with field data through the use of the BET_VH model.     

The performance of churches in the 2012 Emilia earthquakes

In this study the damage suffered by churches during the 2012 Emilia seismic sequence in Italy is analysed, based on surveys and inspections carried out in the area. Similarly to what was observed after other Italian earthquakes, the damage to churches was severe. However, the Emilia churches present some characteristic features such as the use of unreinforced clay brick masonry. In order to appropriately address the performance of this class of buildings, typical architectural layouts and construction techniques are described. Such techniques are interpreted also in the light of the local seismic catalogue. Fifty churches are then selected and their damage is studied, with reference to typical local-collapse mechanisms of different macro-elements. The study highlights that the damage is often concentrated at the top section of the façade, in the clerestory walls, in the vaults and in the bell towers. Structural analyses are performed to explain some of the observations. The overturning of the top section of the façade is analytically addressed, modelling the friction interlocking. With reference to the case study of San Francesco in Mirandola, non-linear static and dynamic analyses allow us to correlate the directionality of damage to the higher seismic demand along the NS direction, to point out the negligible role of the large vertical component of ground motion and to emphasise the relevance of the buttresses for the seismic response of the façade.     

Injection and arrest of dykes: implications for volcanic hazards

Dykes are the principal channels through which magma reaches the surface in volcanic eruptions. For this reason dykes observed in the field are commonly assumed to be feeders to lava flows. The actual proportion of dykes reaching the surface is, however, poorly known. In order to develop models for the purpose of estimating volcanic hazard, this proportion must be known. This follows because such models should not only consider the probability of dykes being injected from magma chambers during periods of unrest in the associated volcanoes, but also the probability of the injected dykes being arrested. This paper presents field data on several thousand dykes from Iceland and Tenerife (Canary Islands) indicating that many, and probably most, dykes become arrested at various crustal levels and never reach the surface to feed eruptions. Using the results of analytical and numerical models, it is shown that, for common loading conditions, the stress field in the vicinity of a magma chamber may favour the injection and propagation of dykes while the stress field at a certain distance from the chamber favours dyke arrest. This means that many dykes that are injected from the chamber propagate only for a very limited distance from the chamber to the point where they become arrested. The implication is that during periods of unrest in volcanoes, the probability of volcanic eruption is only a small fraction of the probability of dyke injection from the source magma chamber.     

Discriminating the long distance dispersal of fine ash from sustained columns or near ground ash clouds: The example of the Pomici di Avellino eruption (Somma-Vesuvius,Italy)

Ash samples from tephra layers correlated with the Pomici di Avellino (Avellino Pumice) eruption of Somma-Vesuvius were collected in distal archives and their composition and particle morphology investigated in order to infer their behaviour of transportation and deposition. Differences in composition and particle morphologies were recognised for ash particles belonging to the magmatic Plinian and final phreatomagmatic phases of the eruption. The ash particles were dispersed in opposite directions during the two different phases of the eruption, and these directions are also different from that of coarse-grained fallout deposits. In particular, ash generated during magmatic phase and injected in the atmosphere to form a sustained column shows a prevailing SE dispersion, while ash particles generated during the final phreatomagmatic phase and carried by pyroclastic density currents show a general NW dispersion. These opposite dispersions indicate an ash dispersal influenced by both high and low atmosphere dynamics. In particular, the magmatic ash dispersal was first driven by stratospheric wind towards NE and then the falling particles encountered a variable wind field during their settling, which produced the observed preferential SE dispersal. The wind field encountered by the rising ash clouds that accompanied the pyroclastic density currents of the final phreatomagmatic phase was different with respect to that encountered by the magmatic ash, and produced a NW dispersal. These data demonstrate how ash transportation and deposition are greatly influenced by both high and low atmosphere dynamics. In particular, fine-grained particles transported in ash clouds of small-scale pyroclastic density currents may be dispersed over distances and cover areas comparable with those injected into the stratosphere by Plinian, sustained columns. This is a point not completely addressed by present day mitigation plans in case of renewal of activity at Somma-Vesuvius, and can yield important information also for other volcanoes potentially characterised by explosive activity.     

An experiment on reservoir representation schemes to improve hydrologic prediction: coupling the national water model with the HEC-ResSim

ABSTRACT

This study experiments with reservoir representation schemes to improve the ability to model active water management in the National Water Model (NWM). For this purpose, we developed an integrated water management model, NWM-ResSim, by coupling the NWM with HEC-ResSim, and two reservoir representation schemes are tested: simulation of reservoir operations and retrieval of scheduled operations. The experiments focus on a pilot reservoir domain in the Russian River basin – Lake Mendocino, California – and its contributing watershed. The evaluation results suggest that the NWM-ResSim improves the simulation performance of reservoir outflow from this managed reservoir over the NWM default level pool routing scheme. The degree of this improvement depends on the suitability of the operation guidance; the reservoir operations simulation scheme could have acceptable errors for the purposes of water resources management, but not for flood operations. Results of the retrieval scheme of scheduled operations demonstrated better performance for sub-daily flood operations.     

Spatio-temporal characteristics of small-scale wave–current ripples on the Ameland ebb-tidal delta

Ebb-tidal deltas are highly dynamic environments affected by both waves and currents that approach the coast under various angles. Among other bedforms of various scales, these hydrodynamics create small-scale bedforms (ripples), which increase the bed roughness and will therefore affect hydrodynamics and sediment transport. In morphodynamic models, sediment transport predictions depend on the roughness height, but the accuracy of these predictors has not been tested for field conditions with strongly mixed (wave–current dominated) forcing. In this study, small-scale bedforms were observed in the field with a 3D Profiling Sonar at five locations on the Ameland ebb-tidal delta, the Netherlands. Hydrodynamic conditions ranged from wave dominated to current dominated, but were mixed most of the time. Small-scale ripples were found on all studied parts of the delta, superimposed on megaripples. Even though a large range of hydrodynamic conditions was encountered, the spatio-temporal variations in small-scale ripple dimensions were relatively small (height 0.015 m, length 0.11 m). Also, the ripples were always highly three-dimensional. These small dimensions are probably caused by the fact that the bed consists of relatively fine sediment. Five bedform height predictors were tested, but they all overestimated the ripple heights, partly because they were not created for small grain sizes. Furthermore, the predictors all have a strong dependence on wave- and current-related velocities, whereas the ripple heights measured here were only related to the near-bed orbital velocity. Therefore, ripple heights and lengths in wave–current-dominated, fine-grained coastal areas ( mm) may be best estimated by constant values rather than values dependent on the hydrodynamics. In the case of the Ameland ebb-tidal delta, these values were found to be m and m. ©2019 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.