REXEL: computer aided record selection for code-based seismic structural analysis

In code-based seismic design and assessment it is often allowed the use of real records as an input for nonlinear dynamic analysis. On the other hand, international seismic guidelines, concerning this issue, have been found hardly applicable by practitioners. This is related to both the difficulty in rationally relating the ground motions to the hazard at the site and the required selection criteria, which do not favor the use of real records, but rather various types of spectrum matching signals. To overcome some of these obstacles a software tool for code-based real records selection was developed. REXEL, freely available at the website of the Italian network of earthquake engineering university labs (), allows to search for suites of waveforms, currently from the European Strong-motion Database, compatible to reference spectra being either user-defined or automatically generated according to Eurocode 8 and the recently released new Italian seismic code. The selection reflects the provisions of the considered codes and others found to be important by recent research on the topic. In the paper, record selection criteria are briefly reviewed first, and then the algorithms implemented in the software are discussed. Finally, via some examples, it is shown how REXEL can effectively be a contribution to code-based real records selection for seismic structural analysis.     

Lava flow hazard at Nyiragongo Volcano,DRC

Mt. Nyiragongo is one of the most dangerous volcanoes in the world for the risk associated with the propagation of lava flows. In 2002 several vents opened along a huge system of fractures, pouring out lava which reached and destroyed a considerable part of Goma, a town of about 500,000 inhabitants on the shore of Lake Kivu. In a companion paper (Favalli et al. in Bull Volcanol, this issue, 2008) we employed numerical simulations of probable lava flow paths to evaluate the lava flow hazard on the flanks of the volcano, including the neighbouring towns of Goma (DRC) and Gisenyi (Rwanda). In this paper we use numerical simulations to investigate the possibility of significantly reducing the lava flow hazard in the city through the construction of protective barriers. These barriers are added to the DEM of the area as additional morphological elements, and their effect is evaluated by repeating numerical simulations with and without the presence of barriers. A parametric study on barrier location, size, shape and orientation led to the identification of barriers which maximize protection while minimizing their impact. This study shows that the highest hazard area corresponding to eastern Goma, which was largely destroyed by lava flows in 2002, cannot be effectively protected from future lava flows towards Lake Kivu and should be abandoned. On the contrary, the rest of the town can be sheltered from lava flows by means of two barriers that deviate or contain the lava within the East Goma sector. A proposal for the future development of the town is formulated, whereby “new” Goma is completely safe from the arrival of lava flows originating from vents outside its boundaries. The proposal minimizes the risk of further destruction in town due to future lava flows.     

Effectiveness of a market-based instrument for the allocation of water in a tropical river environment

The seasonal nature of Australia’s tropical rivers means that connected groundwater aquifers are an important source of both consumptive and non-consumptive water, particularly during the dry season. The management of these common pool groundwater resources is one of the predominant water issues facing northern Australia. A national program of water reform stipulates the expansion of water trading as a key instrument for water allocation. The effectiveness of new institutional arrangements such as water markets will be determined mostly by how well they coordinate with local environmental requirements, local institutions and local norms. This paper describes a novel application of combined field work, institutional analysis, experimental economics and agent-based modeling to the analysis of a potential water market in the Katherine region of the Northern Territory, Australia. The effectiveness of different versions of the policy instrument is assessed in light of local conditions. Instruments that enable personal relationships and local institutions and norms to play a role in water management are found to be more effective in terms of both farming income and environmental impact.     

Fractal generation of surface area of porous media

Many natural porous geological rock formations, as well as engineered porous structures, have fractal properties, i.e., they are self-similar over several length scales. While there have been many experimental and theoretical studies on how to quantify a fractal porous medium and on how to determine its fractal dimension, the numerical generation of a fractal pore structure with predefined statistical and scaling properties is somewhat scarcer. In the present paper a new numerical method for generating a three-dimensional porous medium with any desired probability density function (PDF) and autocorrelation function (ACF) is presented. The well-known Turning Bands Method (TBM) is modified to generate three-dimensional synthetic isotropic and anisotropic porous media with a Gaussian PDF and exponential-decay ACF. Porous media with other PDF's and ACF's are constructed with a nonlinear, iterative PDF and ACF transformation, whereby the arbitrary PDF is converted to an equivalent Gaussian PDF which is then simulated with the classical TBM. Employing a new method for the estimation of the surface area for a given porosity, the fractal dimensions of the surface area of the synthetic porous media generated in this way are then measured by classical fractal perimeter/area relationships. Different 3D porous media are simulated by varying the porosity and the correlation structure of the random field. The performance of the simulations is evaluated by checking the ensemble statistics, the mean, variance and ACF of the simulated random field. For a porous medium with Gaussian PDF, an average fractal dimension of approximately 2.76 is obtained which is in the range of values of actually measured fractal dimensions of molecular surfaces. For a porous medium with a non-Gaussian quadratic PDF the calculated fractal dimension appears to be consistently higher and averages 2.82. The results also show that the fractal dimension is neither strongly dependent of the porosity nor of the degree of anisotropy assumed.     

Probabilistic aftershock hazard assessment II: application of strong ground motion simulations

Probabilistic aftershock hazard assessment (PAHA) has been introduced by Wiemer (Geophys Res Lett 27:3405–3408, 2000). The method, in its original form, utilizes attenuation relations in evaluating peak ground velocity (PGV) exceedence probability. We substitute the attenuation relations together with their uncertainties by strong ground motion simulations for a set of scenarios. The main advantage of such an approach is that the simulations account for specific details of the aftershock source effects (faulting style, slip distribution, position of the nucleation point, etc.). Mean PGVs and their standard deviations are retrieved from the simulation results obtained by the new hybrid k-squared source model, and they are used for the PAHA analysis at a station under study. The model chosen for the testing purposes is inspired by the Izmit A25 aftershock (M _w  = 5.8) that occurred 26 days after the mainshock. The PAHA maps are compared with (1) those obtained by the use of attenuation relations and (2) the peak values of ten selected strong-motion recordings written by the aftershock at epicentral distances <50 km. We conclude that, although the overall hazard decay with increasing fault distance is similar, the PAHA maps obtained by the use of simulations exhibit remanent radiation pattern effect and prolongation in the strike direction due to the directivity effect pronounced for some of the scenarios. As regard the comparison with real data, we conclude that the PAHA maps agree with observed peak values due to appropriate attenuation model adopted in the analysis.     

Truncation of the distribution of ground-motion residuals

Recent studies to assess very long-term seismic hazard in the USA and in Europe have highlighted the importance of the upper tail of the ground-motion distribution at the very low annual frequencies of exceedance required by these projects. In particular, the use of an unbounded lognormal distribution to represent the aleatory variability of ground motions leads to very high and potentially unphysical estimates of the expected level of shaking. Current practice in seismic hazard analysis consists of truncating the ground-motion distribution at a fixed number (ε _max) of standard deviations (σ). However, there is a general lack of consensus regarding the truncation level to adopt. This paper investigates whether a physical basis for choosing ε _max can be found, by examining records with large positive residuals from the dataset used to derive one of the ground-motion models of the Next Generation Attenuation (NGA) project. In particular, interpretations of the selected records in terms of causative physical mechanisms are reviewed. This leads to the conclusion that even in well-documented cases, it is not possible to establish a robust correlation between specific physical mechanisms and large values of the residuals, and thus obtain direct physical constraints on ε _max. Alternative approaches based on absolute levels of ground motion and numerical simulations are discussed. However, the choice of ε _max is likely to remain a matter of judgment for the foreseeable future, in view of the large epistemic uncertainties associated with these alternatives. Additional issues arise from the coupling between ε _max and σ, which causes the truncation level in terms of absolute ground motion to be dependent on the predictive equation used. Furthermore, the absolute truncation level implied by ε _max will also be affected if σ is reduced significantly. These factors contribute to rendering a truncation scheme based on a single ε _max value impractical.     

综合RS与GIS方法提取青藏高原冲川变化信息：以布喀塔格峰为例

冰川变化的常规观测方法在高原的大部分地区无法实施。遥感与地理信息系统技术的发展为研究冲川变化提供了有效的手段,本文以位于昆仑山脉中段的布喀塔格山峰冰川为例,利用１９７３年至１９９４年的ＲＢＶ,ＭＳＳ和ＴＭ遥感资料为信息源,综合目视判读与统计分析方法提取冰川界限,形成冰川边界图,并规一化坐标系统,以ＧＩＳ为工具分析该冲川群的变化情况,计算冲川进退变化速率。     

重庆市经济发展水平的区域差异研究

建立了重庆市区域经济发展水平的指标体系,利用层次分析法确定指标权重,根据１９９６年统计数据对全市４３个区（市,县）进行综合评价,并依据综合评价价值划分４个经济发展梯度区,据此揭示了区域差异的空间特点,并讨论了区域差异的成因及缩小区域差异的对策。     

星状沙丘研究进展综述

本文阐述了近年来的星状沙丘方面的研究进展，介绍了国内外星状沙丘的形成发育模式，粒度特征和沉积构造特征方面的理论及研究进展。