Numerical modelling of the effects of foundation scour on the response of a bridge pier

Foundation scour can have a detrimental effect on the performance of bridge piers, inducing a significant reduction of the lateral capacity of the footing and accumulation of permanent settlement and rotation. Although the hydraulic processes responsible for foundation scour are nowadays well known, predicting their mechanical consequences is still challenging. Indeed, its impact on the failure mechanisms developing around the foundation has not been fully investigated. In this paper, numerical simulations are performed to study the vertical and lateral response of a scoured bridge pier founded on a cylindrical caisson foundation embedded in a layer of dense sand. The sand stress–strain behaviour is reproduced by employing the Severn-Trent model. The constitutive model is firstly calibrated on a set of soil element tests, including drained and undrained monotonic triaxial tests and resonant column tests. The calibration procedure is implemented considering the stress and strain nonuniformities within the samples, by simulating the laboratory tests as boundary value problems. The numerical model is then validated against the results of centrifuge tests. The results of the simulations are in good agreement with the experimental results in terms of foundation capacity and settlement accumulation. Moreover, the model can predict the effects of local and general scour. The numerical analyses also highlight the impact of scouring on the failure mechanisms, revealing that the soil resistance depends on the hydraulic scenario.

     

Variability in drift macroalgal abundance in relation to biotic and abiotic factors in two seagrass dominated estuaries in the Western Gulf of Mexico

We examined the spatial and temporal variability in drift macroalgal abundance in two seagrass dominated estuarine systems on the Texas coast: Redfish Bay (in the Copano-Aransas Estuary) and Lower Laguna Madre. Measurements of benthic macroalgal variability were made in conjunction with a suite of biotic (seagrass biomass, percent cover, blade width and length, shoot density, epiphyte biomass, seagrass blade C:N ratios, and drift macroalgal abundance and composition) and abiotic (inorganic nitrogen and phosphorus concentrations, chlorophylla, total suspended solids, light attenuation, salinity, temperature, total organic carbon and porewater NH₄ ⁺) indicators. All parameters were measured at 30 sites within each estuary semiannually from July 2002 to February 2004. Principal components analysis (PCA) was used to examine relationships between drift macroalgal abundance and biotic and abiotic parameters. In both Redfish Bay and Lower Laguna Madre, drift macroalgal distribution was widespread, and during three of four sampling periods, abundance was equal to abovegro und biomass ofThalassia testudinum, the dominant seagrass. Drift macro algal abundance was highly variable within sites, between sites, and between seasons in both estuaries. No significant differences in drift macroalgal abundance were found between Redfish Bay and Lower Laguna Madre. In Redfish Bay, drift macroalgae (90.1±10.2 gm⁻²) tended to accumulate in bare patches within seagrass beds. In Lower Laguna Madre, drift macroalgae (72.7±10.7 gm⁻²) tended to accumulate in areas of dense seagrass cover rather than in bare areas. We found no relationship between drift macroalgal abundance and low (<2μM) water column nutrient concentrations, and although several of our measured parameters were related to drift macroalgal abundance, none alone sufficiently explained the variability in abundance noted between the two estuarine systems. The contrasting patterns of macroalgal accumulation between Redrish Bay and Lower Laguna Madre likely reflect differences in water circulation characteristics between the two regions as dictated by local physiography, in cluding the shape and orientation of the lagoons, with seasonal variations in macroalgal abundance related to changes in freshwater inflow and nutrient loading.     

Performance of masonry buildings during the Emilia 2012 earthquake

The earthquake sequence started on May \(20\) th 2012 in Emilia (Italy) affected a region where masonry constructions represent a large part of the existing building stock and the construction of new modern masonry buildings is a common practice. The paper is focused on the performance of common architectural configurations, typical for residential or business use. The large majority of old masonry buildings is made of fired clay bricks. The seismic performance of these buildings is particularly interesting since major past earthquakes in Italy affected areas with mainly stone masonry structures. Apart from examples showing systematic or peculiar structural deficiencies governing the vulnerability of several buildings, the overall seismic performance of these structures to repeated shaking, with PGA as large as 0.25–0.3 g was rather good, despite the major part of them were only conceived for carrying vertical loads. In fact, seismic design is mandatory in the area only since 2003. Modern low-rise masonry buildings erected after this date and incorporating seismic design and proper detailing resulted in most cases practically undamaged. The examples reported in the paper allow an evaluation of the superior performance of seismically designed modern masonry buildings in comparison to older ones.     

The added value of joint PP–PS inversion for reservoir characterization: A case study using Jubarte PRM seismic data (offshore Brazil)

The added value of the joint pre-stack inversion of PP (incident P-wave and reflected P-wave) and PS (incident P-wave and reflected S-wave) seismic data for the time-lapse application is shown. We focus on the application of this technique to the time-lapse (four-dimensional) multicomponent Jubarte field permanent reservoir monitoring seismic data. The joint inversion results are less sensitive to noise in the input data and show a better match with the rock physics models calibrated for the field. Further, joint inversion improves S-impedance estimates and provides a more robust quantitative interpretation, allowing enhanced differentiation between pore pressure and fluid saturation changes, which will be extremely useful for reservoir management. Small changes in reservoir properties are expected in the short time between the time-lapse seismic acquisitions used in the Jubarte project (only 1 year apart). The attempt to recover subtle fourth-dimensional effects via elastic inversion is recurrent in reservoir characterization projects, either due to the small sensitivity of the reservoirs to fluid and pressure changes or the short interval between the acquisitions. Therefore, looking for methodologies that minimize the uncertainty of fourth-dimensional inversion outputs is of fundamental importance. Here, we also show the differences between PP only and joint PP–PS inversion workflows and parameterizations that can be applied in other projects. We show the impact of using multicomponent data as input for elastic seismic inversions in the analysis of the time-lapse differences of the elastic properties. The larger investment in the acquisition and processing of multicomponent seismic data is shown to be justified by the improved results from the fourth-dimensional joint inversion.     

Flow time estimation with spatially variable hillslope velocity in ungauged basins

The automated spatial estimation of the hillslope runoff dynamics is used as a valuable tool for the estimation of the travel time distribution (flow time), a major factor for the hydrologic prediction in ungauged basins. In fact, while the flow time function is usually obtained by rescaling the flow paths with constant channel and hillslope velocities, in this work a spatially distributed kinematic component, as a function of terrain features and in particular slope and land use, is implemented and its influence on the hydrologic response is tested by means of the Width Function Instantaneous Unit Hydrograph (WFIUH) framework. Hillslope surface flow velocities are evaluated by applying different uniform flow formulas within an automated DEM-based (terrain analysis) algorithm. A comparison test of the performances of the Manning, Darcy, Maidment and Soil Conservation Service uniform flow equations is performed for several case studies in Italy pertaining to different climatic and geomorphic conditions. Results provide new insights for a better understanding of the flow time function also introducing a more parsimonious and physically-based calibration scheme of the WFIUH.     

Geophysical Finite-Element Simulation Tool (GeoFEST): Algorithms and Validation for Quasistatic Regional Faulted Crust Problems

GeoFEST (Geophysical Finite Element Simulation Tool) is a two- and three-dimensional finite element software package for the modeling of solid stress and strain in geophysical and other continuum domain applications. It is one of the featured high-performance applications of the NASA QuakeSim project. The program is targeted to be compiled and run on UNIX systems, and is running on diverse systems including sequential and message-passing parallel systems. Solution to the elliptical partial differential equations is obtained by finite element basis sampling, resulting in a sparse linear system primarily solved by conjugate gradient iteration to a tolerance level; on sequential systems a Crout factorization for the direct inversion of the linear system is also supported. The physics models supported include isotropic linear elasticity and both Newtonian and power-law viscoelasticity, via implicit quasi-static time stepping. In addition to triangular, quadrilateral, tetrahedral and hexahedral continuum elements, GeoFEST supports split-node faulting, body forces, and surface tractions. This software and related mesh refinement strategies have been validated on a variety of test cases with rigorous comparison to analytical solutions. These include a box-shaped domain with imposed motion on one surface, a pair of strike slip faults in stepover arrangement, and two community-agreed benchmark cases: a strike slip fault in an enclosing box, and a quarter-domain circular fault problem. Scientific applications of the code include the modeling of static and transient co- and post-seismic earth deformation, Earth response to glacial, atmospheric and hydrological loading, and other scenarios involving the bulk deformation of geologic media.     

Experimental solidification of anhydrous latitic and trachytic melts at different cooling rates: The role of nucleation kinetics

Two sets of cooling experiments were run at atmospheric conditions for two anhydrous starting latitic and trachytic melts: 1) five cooling rates (25, 12.5, 3, 0.5, and 0.125 °C/min) between 1300° and 800 °C, and 2) a 0.5 °C/min cooling rate from 1300 °C with quench temperatures at 1200°, 1100°, 1000° and 900 °C. Trachytic run-products are invariably glassy. Nucleation is also suppressed in the latitic run-products at the three highest cooling rates. Conversely, in the 0.5 and 0.125 °C/min runs, latites have a crystal content of  90 vol.%. The phases are: plagioclase, clinopyroxene, glass and iron-bearing oxide (in order of abundance). The variable quench temperatures, investigated by coupling experiments with Pt wire and Pt capsule sample containers in set 2, again did not produce crystallization of trachyte, whereas latitic samples are characterized by  10 vol.% of oxides, pyroxenes and plagioclase (in order of appearance), at temperature < 1000 °C. Effects of (preferential) heterogeneous nucleation on sample holders, of superheating degree, and chemical species loss during cooling are absent for both melt compositions. The difference of solidification paths between these two silicate melts can be ascribed only to their small chemical differences. In comparison with calculated equilibrium conditions all the experimental latitic and trachytic run-products revealed strong kinetic effects, interpretable in the light of the nucleation theory. The glass-forming ability (GFA) of trachyte is higher, whereas their critical cooling rate (Rc) is lower (< 0.125 °C/min), in comparison to latitic melts (Rc > 0.5 °C/min). The experimental results carried out in this study can be applied to lava flows and domes; trachytic lavas are able to flow for longer period with respect to latitic ones in a metastable condition. Glass-rich terrestrial lavas, i.e. obsidians, can be the result of sluggish nucleation kinetics due to the relative high polymerisation of evolved silicate melts.     

The influence of slope gradient and gully channel on the run-out behavior of rockslide-debris flow: an analysis on the Verghereto landslide in Italy

Landslides - Rockslide-debris flow may cause catastrophic damages because of its high speed and long run-out distance. The influences of slope gradient and gully channel on the entrainment,...     

The effect of crystallization on the rheology and dynamics of lava flows

The dynamics of a lava flow is studied by a two-dimensional model describing a viscous fluid with Bingham rheology, flowing down a slope. The temperature in the flow is calculated assuming that heat is transferred through the plug by conduction and is lost by radiation to the atmosphere at the top of the flow. Taken into account is that the increasing crystallization takes place in the flow as a consequence of cooling. The lava viscosity and yield stress are expressed as a function of crystallization degree as well as of temperature: in particular it is assumed that yield stress reaches a maximum value above the solidus temperature, according to experimental data. Dynamical variables, such as velocity and thickness of the flow, are calculated for different values of the maximum crystallization degree and the flow rate. The model shows how the lava flow dynamics is affected by cooling and crystallization. The cooling of the flow is controlled by the increase of yield stress, which produces a thicker plug and makes the heat loss slower. The increasing crystallization has two opposing effects on viscosity: it produces an increase of viscosity, but at the same time produces an increase of yield stress and hence reduces the heat loss and keeps the internal temperature high. As a consequence, lava flows are significantly affected by the dependence of yield stress on temperature and scarcely by the maximum crystallization degree.