Elastic waves in nonhomogeneous media under 2D conditions: II. Numerical implementation

The purpose of this work is to present three methods of analysis for elastic waves propagating in two dimensional, elastic nonhomogeneous media. The first step, common to all methods, is a transformation of the governing equations of motion so that derivatives with respect to the material parameters no longer appear in the differential operator. This procedure, however, restricts analysis to a very specific class of nonhomogeneous media, namely those for which Poisson's ratio is equal to 0.25 and the elastic parameters are quadratic functions of position. Subsequently, fundamental solutions are evaluated by: (i) conformal mapping in conjunction with wave decomposition, which in principle allows for both vertical and lateral heterogeneities; (ii) wave decomposition into pseudo-dilatational and pseudo-rotational components, which results in an Euler-type equation for the transformed solution if medium heterogeneity is a function of one coordinate only; and (iii) Fourier transformation followed by a first order differential equation system solution, where the final step involving inverse transformation from the wavenumber domain is accomplished numerically. Finally, in the companion paper numerical examples serve to illustrate the above methodologies and to delineate their range of applicability.     

Scattering of seismic waves by cracks in multi-layered geological regions: I. Mechanical model

In this work, a hybrid boundary integral equation method (BIEM) is developed, based on both displacement and hypersingular traction formulations, for the analysis of time-harmonic seismic waves propagating through cracked, multi-layered geological regions with surface topography and under plane strain conditions. Specifically, the displacement-based BIEM is used for a multi-layered deposit with interface cracks, while the regularized, traction-based BIEM is used when internal cracks are present within the layers. The standard uni-dimensional boundary element with parabolic shape functions is employed for discretizing the free surface and the layer interfaces, while special discontinuous boundary elements are placed near the crack tips to model the asymptotic behaviour of both displacements and tractions. This formulation yields displacement amplitudes and phase angles on the free surface of a geological deposit, as well as stress intensity factors near the tips of the cracks. Finally, in the companion paper, numerical results are presented which show that both scattered wave and stress concentration fields are sensitive to the incidence seismic wave parameters and to specific site conditions such as surface topography, layering, the presence of cracks and crack interaction.     

Seismic response of lined tunnels in the half-plane with surface topography

In this work, we examine the seismic response of multiple tunnels reinforced with liners and buried within the elastic homogeneous half-plane in the presence of surface relief. The seismic waves are upward propagating, time-harmonic, horizontally polarized shear (SH) waves. More specifically, we examine: (a) the scattered wave fields along the free surface and inside the half-plane with the embedded tunnels; (b) the dynamic stress concentration factors that develop at the soil-liner interfaces; (c) the stresses and displacements that develop inside the tunnel liners. We use a sub-structuring technique that is based on the direct boundary element method to model each constituent part of the problem separately. Then, assembly of the full problem is accomplished through the imposition of compatibility and equilibrium conditions at all interfaces. Next, a detailed verification study is carried out based on comparisons against available analytical and/or numerical results for a series of test examples. Subsequently, detailed numerical simulations are conducted and the results of these parametric studies reveal the influence of the following key parameters on the soil-tunnel system response: (a) the shape of the free-surface relief; (b) the depth of placement of the tunnels and their separation distance; (c) the SH-wavelength to tunnel diameter ratio; (d) the elastic properties of the tunnel lining rings and (e) the dynamic interaction effects between the free-surface relief and the tunnels.     

Performance‐based optimum seismic design of reinforced concrete structures

A fully automated design methodology based on nonlinear response history analysis is proposed for the optimum seismic design of reinforced concrete (RC) structures. The conventional trial‐and‐error process is replaced by a structural optimization algorithm that serves as a search engine capable of locating the most efficient design in terms of cost and performance. Two variations of the proposed design methodology are introduced. The first approach treats the optimum design problem in a deterministic manner, while in the second variation the optimum design is sought in the framework of a reliability‐based optimization problem. The reliability‐based approach seems to be a more rational procedure since more meaningful design criteria that correlate better with the performance‐based design concept can be adopted. Thus, the practice of using the mean annual frequency of a limit‐state being exceeded to assess the candidate designs is compared with the use of deterministic criteria. Both formulations take into consideration the structural response for a number of limit‐states, from serviceability to collapse prevention. The proposed design procedure is specifically tailored to the design of RC structures, where a preliminary design step of generating tables of concrete sections is introduced. In order to handle the large size of the tables, the concept of multi‐database cascade optimization is implemented. The final design has to comply with the provisions of European design codes. The proposed methodology allows for a significant reduction of the direct construction cost combined with improved control of the seismic performance under earthquake loading. Copyright © 2008 John Wiley & Sons, Ltd.     

Building design based on energy dissipation: a critical assessment

The basic objective of this study is the assessment of the European seismic design codes and in particular of EC2 and EC8 with respect to the recommended behaviour factor q. The assessment is performed on two reinforced concrete multi-storey buildings, having symmetrical and non-symmetrical plan view respectively, which were optimally designed under four different values of the behaviour factor. In the mathematical formulation of the optimization problem the initial construction cost is considered as the objective function to be minimized while the cross sections and steel reinforcement of the beams and the columns constitute the design variables. The provisions of Eurocodes 2 and 8 are imposed as constraints to the optimization problem. Life-cycle cost analysis, in conjunction with structural optimization, is believed to be a reliable procedure for assessing the performance of structures during their life time. The two most important findings that can be deduced are summarized as follows: (1) The proposed Eurocode behaviour factor does not lead to a more economical design with respect to the total life-cycle cost compared to other values of q (q = 1, 2). (2) The differences of the total life-cycle cost values may be substantially greater than those observed for the initial construction cost for four different q (q = 1, 2, 3, 4).     

Are Solar Active Regions with Major Flares More Fractal, Multifractal, or Turbulent Than Others?

Multiple recent investigations of solar magnetic-field measurements have raised claims that the scale-free (fractal) or multiscale (multifractal) parameters inferred from the studied magnetograms may help assess the eruptive potential of solar active regions, or may even help predict major flaring activity stemming from these regions. We investigate these claims here, by testing three widely used scale-free and multiscale parameters, namely, the fractal dimension, the multifractal structure function and its inertial-range exponent, and the turbulent power spectrum and its power-law index, on a comprehensive data set of 370 timeseries of active-region magnetograms (17?733 magnetograms in total) observed by SOHO’s Michelson Doppler Imager (MDI) over the entire Solar Cycle 23. We find that both flaring and non-flaring active regions exhibit significant fractality, multifractality, and non-Kolmogorov turbulence but none of the three tested parameters manages to distinguish active regions with major flares from flare-quiet ones. We also find that the multiscale parameters, but not the scale-free fractal dimension, depend sensitively on the spatial resolution and perhaps the observational characteristics of the studied magnetograms. Extending previous works, we attribute the flare-forecasting inability of fractal and multifractal parameters to i)?a?widespread multiscale complexity caused by a possible underlying self-organization in turbulent solar magnetic structures, flaring and non-flaring alike, and ii)?a?lack of correlation between the fractal properties of the photosphere and overlying layers, where solar eruptions occur. However useful for understanding solar magnetism, therefore, scale-free and multiscale measures may not be optimal tools for active-region characterization in terms of eruptive ability or, ultimately, for major solar-flare prediction.     

Evaluation of a global ensemble flood prediction system in Peru

ABSTRACT

Flood early warning systems play a more substantial role in risk mitigation than ever before. Hydrological forecasts, which are an essential part of these systems, are used to trigger action against floods around the world. This research presents an evaluation framework, where the skills of the Global Flood Awareness System (GloFAS) are assessed in Peru for the years 2009–2015. Simulated GloFAS discharges are compared against observed ones for 10 river gauges. Forecasts skills are assessed from two perspectives: (i) by calculating verification scores at every river section against simulated discharges and (ii) by comparing the flood signals against reported events. On average, river sections with higher discharges and larger upstream areas perform better. Raw forecasts provide correct flood signals for 82% of the reported floods, but exhibit low verification scores. Post-processing of raw forecasts improves most verification scores, but reduces the percentage of the correctly forecasted reported events to 65%.     

The COVID-19 pandemic and global environmental change: Emerging research needs

The outbreak of COVID-19 raised numerous questions on the interactions between the occurrence of new infections, the environment, climate and health. The European Union requested the H2020 HERA project which aims at setting priorities in research on environment, climate and health, to identify relevant research needs regarding Covid-19. The emergence and spread of SARS-CoV-2 appears to be related to urbanization, habitat destruction, live animal trade, intensive livestock farming and global travel. The contribution of climate and air pollution requires additional studies. Importantly, the severity of COVID-19 depends on the interactions between the viral infection, ageing and chronic diseases such as metabolic, respiratory and cardiovascular diseases and obesity which are themselves influenced by environmental stressors. The mechanisms of these interactions deserve additional scrutiny. Both the pandemic and the social response to the disease have elicited an array of behavioural and societal changes that may remain long after the pandemic and that may have long term health effects including on mental health. Recovery plans are currently being discussed or implemented and the environmental and health impacts of those plans are not clearly foreseen. Clearly, COVID-19 will have a long-lasting impact on the environmental health field and will open new research perspectives and policy needs.     

The influence of anisotropy on compaction bands: The case of coaxiality between stress and fabric anisotropy tensors

Compaction bands are localized failure patterns that appear in highly porous rock material under the effect of relatively high confining pressure. Being affected mainly by volumetric compression, these bands appear to be almost perpendicular to the most compressive principal stress of a stress state at the so-called “cap” of the yield surface (YS). In this study, we focus on the mechanism that leads to the onset of compaction bands by using a viscoplasticity model able to describe the post-localization response of these materials. The proposed constitutive framework is based on the overstress theory of Perzyna (1966) and the anisotropic clay plasticity model of Dafalias (1986), which provides not only the necessary “cap” of the YS, but introduces a rotational hardening (RH) mechanism, thus, accounting for the effect of fabric anisotropy. Following the analysis of Veveakis and Regenauer-Lieb (2015), we identify the compaction bands as “static” cnoidal wave formations in the medium that occur at a post-yield regime, and we study the effect of rotational and isotropic hardening on their onset. Moreover, we determine a theoretical range of confining pressures in triaxial compression tests for the compaction bands to develop. Under the assumption of coaxiality between stress and anisotropy tensors, the results show that the isotropic hardening promotes compaction localization, whereas the RH has a slightly negative effect on the onset of compaction localization.