A viscoelastic shear zone model of compressional and extensional plate boundaries

A model is proposed describing the mechanical evolution of a shear zone along compressional and extensional plate boundaries, subject to constant strain rate. The shear zones are assumed as viscoelastic with Maxwell rheology and with elastic and rheological parameters depending on temperature and petrology. Stress and strain are computed as functions of time and depth. For both kinds of boundaries the model reproduces the existence of a shallow seismogenic zone, characterized by a stress concentration. The thickness of the seismogenic layer is evaluated considering the variations of shear stress and frictional strength on faults embedded in the shear zone. Assuming that a fault dislocation takes place, the brittle-ductile transition is assumed to occur at the depth at which the time derivative of total shear stress changes from positive to negative values. The effects of different strain rates and geothermal gradients on the depth of the brittle-ductile transition are studied. The model predictions are consistent with values inferred from seismicity data of different boundary zones.     

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.     

Rock grade control structures and stepped gabion weirs: Scour analysis and flow features

Rock and stepped gabion weirs are peculiar hydraulic structures that received relatively little attention in technical literature. Nevertheless, they can be successfully used for river restoration instead of traditional hydraulic structures. They have the advantage of being elastic structures and to preserve the natural environment. They can easily adapt to the in situ conditions and can be effortlessly modified according to the different hydraulic or geometric conditions which can occur in a natural river. The present study aims to analyze the effects of their presence on flow pattern and on the scour hole occurring downstream. The analysis involved scour processes, hydraulic jump types, stilling basin morphology and flow patterns. Two different hydraulic jump types were distinguished and classified. It was shown that the flow regime deeply influences the scour process, which evolves much more rapidly when a Skimming Flow regime takes place. Empirical relationships are proposed to evaluate maximum scour depth, maximum axial length, and non dimensional axial profiles.     

Low Visibility Formation and Forecasting on the Northern Coast of Brazil

Visibility analysis and forecast at the Maceio International Airport in the Brazilian Northeast (NEB) was the principal goal of this investigation. Surface meteorological data of the Maceio International Airport were used for low visibility frequency study. Low visibility in NEB was provoked more frequently by light fog (LF) formation (1,098 or 92 h month^?1 on average). Haze and fog were very rare (81 h and one event per year, respectively on average). Light fog with a visibility less than 2 km usually was detected together with rain or drizzle. Low visibility was observed more frequently at night and during the rainy season. Applications of the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model for light fog forecast were tested. Thermodynamic processes were studied by vertical profile, elaborated by: (1) National Centers for Environmental Prediction (NCEP) reanalysis data for Maceio (because of some radiosonde absence) and (2) forecast vertical temperature and humidity profiles were produced, using Air Parcels Trajectories of the HYSPLIT model at the pattern levels. The synoptic situations before and during low visibility phenomena were analyzed using different products of NCEP reanalysis, the high resolution (10 km) ETA model and infrared satellite images. Wave disturbance in the trade winds field, localized on the northwest periphery of the South Atlantic subtropical High, usually accompanied the phenomena. A humidity advection, weak ascendant movement and thermal inversion absence at the low levels were created by these waves. The middle level’s descendent movement provoked the humidity accumulation at levels below. Satisfactory results of the HYSPLIT model applications for light fog forecast were obtained with 12 h antecedence. In particular, stable level forecast by the ETA model was forecast satisfactorily with 12 h antecedence; vertical movements were predicted better with up to 48 h antecedence. The PSU/NCAR mesoscale model (MM5) and PAFOG models were tested for analysis and forecast of an intensive fog event. Intensive fog provoked a fatal accident of a small airplane near the Maceio Airport in 2007. These fog formation processes were studied by NCEP reanalysis data, the high resolution regional model MM5, and satellite and radar data. Fog formation was simulated by PAFOG model and satisfactory results were obtained with 10 h antecedence.     

On the SAR derived alert in the detection of oil spills according to the analysis of the EGEMP

Satellite services that deliver information about possible oil spills at sea currently use different labels of “confidence” to describe the detections based on radar image processing. A common approach is to use a classification differentiating between low, medium and high levels of confidence. There is an ongoing discussion on the suitability of the existing classification systems of possible oil spills detected by radar satellite images with regard to the relevant significance and correspondence to user requirements. This paper contains a basic analysis of user requirements, current technical possibilities of satellite services as well as proposals for a redesign of the classification system as an evolution towards a more structured alert system. This research work offers a first review of implemented methodologies for the categorisation of detected oil spills, together with the proposal of explorative ideas evaluated by the European Group of Experts on satellite Monitoring of sea-based oil Pollution (EGEMP).     

Evaluation of force-based and displacement-based out-of-plane seismic assessment methods for unreinforced masonry walls through refined model simulations

The performance of force-based and displacement-based seismic assessment methods for the life-safety limit state check of out-of-plane loaded unreinforced masonry walls is evaluated on the basis of refined numerical simulations. For this purpose, a discrete element model of a vertically spanning wall is built and validated against experimental results from static and dynamic test conditions. The model is then analysed for a large range of wall configurations. For each configuration, a static pushover analysis and a series of incremental dynamic analyses are run, the latter permitting to determine the capacity of the wall under dynamic loading. The accuracy of the assessment methods in predicting the acceleration at which the walls collapse is evaluated. It is found that the displacement-based method is more accurate, robust, and safe than the force-based method. The comparison also shows that for walls characterised by a relatively high ratio of axial load to Euler's critical load, both assessment methods lead to an overestimation of the wall capacity. As a remedy, a modification to the methods based on a recently developed mechanical model is put forward and tested. For the force-based method, it is additionally suggested to set for walls with relatively high overburden ratios the behaviour factor equal to 1. To ensure reproducibility of this study, all input and output files of the numerical simulations are made publicly available.     

On the correction of spectra by a displacement reduction factor in direct displacement-based seismic design and assessment

The application of some design and assessment approaches, such as the direct displacement-based design (DDBD) and the capacity spectrum methods, requires the modification of elastic design spectra by some displacement reduction factor, to account for the appropriate energy dissipation capacity of different structures. While several equations to correlate dissipation and hysteresis cycles are available, once the displacement reduction factor has been obtained, the correction of the spectra is operated reducing the displacement demand accordingly and conserving the period of vibration at each point. This procedure is here discussed and proved to be conceptually inappropriate, because the spectral acceleration rather than the period should be kept at each point. The application of this alternative procedure may result in increased shear strength demand in design and in larger required displacement capacity for the same level of strength in assessment, if all other factors are not modified. However, the calibration of the reduction factors applied in DDBD has been extensive, and the method has proved to be effective in predicting displacement demands consistent with those resulting from refined nonlinear time history analysis; therefore, a possible introduction of the proposed correction will require equally extensive studies and possibly compensating corrections in the calculation of the equivalent damping. On the contrary, an appropriate calibration of the factors to be used in the application of the “capacity spectrum” method is still being developed, and the consideration of a constant acceleration may facilitate the derivation of effective equations.     

Addressing non-uniqueness in linearized multichannel surface wave inversion

The multichannel analysis of the surface waves method is based on the inversion of observed Rayleigh-wave phase-velocity dispersion curves to estimate the shear-wave velocity profile of the site under investigation. This inverse problem is nonlinear and it is often solved using 'local' or linearized inversion strategies. Among linearized inversion algorithms, least-squares methods are widely used in research and prevailing in commercial software; the main drawback of this class of methods is their limited capability to explore the model parameter space. The possibility for the estimated solution to be trapped in local minima of the objective function strongly depends on the degree of nonuniqueness of the problem, which can be reduced by an adequate model parameterization and/or imposing constraints on the solution.
In this article, a linearized algorithm based on inequality constraints is introduced for the inversion of observed dispersion curves; this provides a flexible way to insert a priori information as well as physical constraints into the inversion process. As linearized inversion methods are strongly dependent on the choice of the initial model and on the accuracy of partial derivative calculations, these factors are carefully reviewed. Attention is also focused on the appraisal of the inverted solution, using resolution analysis and uncertainty estimation together with a posteriori effective-velocity modelling. Efficiency and stability of the proposed approach are demonstrated using both synthetic and real data; in the latter case, cross-hole S-wave velocity measurements are blind-compared with the results of the inversion process.     

Transient processes in Stromboli’s shallow basaltic system inferred from dolerite and magmatic breccia blocks erupted during the 5 April 2003 paroxysm

We describe the mineralogy, geochemistry, and mesomicrostructure of fresh subvolcanic blocks erupted during the 5 April 2003 paroxysm of Stromboli (Aeolian Islands, Italy). These blocks represent ∼50 vol.% of the total erupted ejecta and consist of fine- to medium-grained basaltic lithotypes ranging from relatively homogeneous dolerites to strongly or poorly welded magmatic breccias. The breccia components are represented by angular fragments of dolerites entrapped in a matrix of vesiculated (lava-like to scoriae) crystal-rich (CR) basalt. All of the studied blocks are cognates with the CR basalt of the normal Strombolian activity or lavas and they are often coated by a few-centimeter thick layer of crystal-poor (CP) basaltic pumice erupted during the paroxysm. We suggest that they result from the rapid increase of pressure and related subvolcanic rock failure that occurred shortly before the 5 April 2003 explosion, when the uppermost portion of the edifice inflated and suffered brecciation as the result of the sudden rise of the gas-rich CP basalt that triggered the eruption. Dolerites and magmatic matrix of the breccias show major and trace element compositions that match those of the CR basalts erupted during normal Strombolian activity and effusive events at Stromboli volcano. Dolerites consist of (a) phenocrysts normally found in the CR basalts and (b) late-stage magmatic minerals such as sanidine, An_60-28 plagioclase, Fe–Mn-rich olivines (Fo_68-48), phlogopite, apatite, and opaque mineral pairs (magnetite and ilmenite), most of which are never found both in lava flows and scoriae erupted during the persistent explosive activity that characterizes typical Strombolian behavior. Subvolcanic crystallization of the Stromboli CR magma, leading to slowly cooled equivalents of basalts, could result from transient drainage of the magma from the summit craters to lower levels. Fingering and engulfing of the material that collapsed from the summit crater floor into the shallow basaltic system during the late evening of 28 December 2002 coupled with the short break in the summit persistent explosions between December 2002 and March 2003 permitted the CR magma pockets to solidify as dolerites, which were confined to the uppermost portion of the system and thus not involved in the ongoing flank effusive activity. Crystal size distribution of the basaltic blocks and crystallization of the finer-grained (<0.1 mm) mafic minerals of the dolerites over a time interval of ∼100 days closely agrees with the above interpretation. Vesicle filling (miarolitic cavities) locally found in some dolerites, with minerals deposited as vapor-phase crystallization is a result of continuous gas percolation through the rocks of the uppermost portion of the volcanic system. Poorly welded magmatic breccias formed during syn-eruptive processes of 5 April 2003, when the paroxysm strongly shattered the shallow subvolcanic system and many dolerite fragments were entrapped in the CR magma. In contrast, the high degree of welding between the dolerite clasts and the CR basaltic matrix in the strongly welded magmatic breccias provides a snapshot of subvolcanic intrusions of the CR basalt into the dolerite when, after a 2-month break in activity, CR magmas started to rise again to the summit craters. Blocks similar to these subvolcanic ejecta of 5 April 2003 were also erupted during previous paroxysms (e.g., 1930) suggesting that changes in the usual Strombolian activity (e.g., short breaks in the persistent mild explosions and/or flank effusive activity) lead to transient crystallization of dolerites in the shallow plumbing system.     

Long‐term sediment yield from a small catchment in southern Brazil affected by land use and soil management changes

Sediments produced from eroding cultivated land can cause on‐site and off‐site effects that cause considerable economic and social impacts. Despite the importance of soil conservation practices (SCP) for the control of soil erosion and improvements in soil hydrological functions, limited information is available regarding the effects of SCP on sediment yield (SY) at the catchment scale. This study aimed to investigate the long‐term relationships between SY and land use, soil management, and rainfall in a small catchment. To determine the effects of anthropogenic and climatic factors on SY, rainfall, streamflow, and suspended sediment concentration were monitored at 10‐min intervals for 14 years (2002–2016), and the land use and soil management changes were surveyed annually. Using a statistical procedure to separate the SY effects of climate, land use, and soil management, we observed pronounced temporal effects of land use and soil management changes on SY. During the first 2 years (2002–2004), the land was predominantly cultivated with tobacco under a traditional tillage system (no cover crops and ploughed soil) using animal traction. In that period, the SY reached approximately 400 t·km^?2·year^?1. From 2005 to 2009, a soil conservation programme introduced conservation tillage and winter cover crops in the catchment area, which lowered the SY to 50 t·km^?2·year^?1. In the final period (2010–2016), the SCP were partially abandoned by farmers, and reforested areas increased, resulting in an SY of 150 t·km^?2·year^?1. This study also discusses the factors associated with the failure to continue using SCP, including structural support and farmer attitudes.