Comment on Pisarenko et al., “Characterization of the Tail of the Distribution of Earthquake Magnitudes by Combining the GEV and GPD Descriptions of Extreme Value Theory”

In this short note, I comment on the research of Pisarenko et al. (Pure Appl. Geophys 171:1599–1624, 2014) regarding the extreme value theory and statistics in the case of earthquake magnitudes. The link between the generalized extreme value distribution (GEVD) as an asymptotic model for the block maxima of a random variable and the generalized Pareto distribution (GPD) as a model for the peaks over threshold (POT) of the same random variable is presented more clearly. Inappropriately, Pisarenkoet al. (Pure Appl. Geophys 171:1599–1624, 2014) have neglected to note that the approximations by GEVD and GPD work only asymptotically in most cases. This is particularly the case with truncated exponential distribution (TED), a popular distribution model for earthquake magnitudes. I explain why the classical models and methods of the extreme value theory and statistics do not work well for truncated exponential distributions. Consequently, these classical methods should be used for the estimation of the upper bound magnitude and corresponding parameters. Furthermore, I comment on various issues of statistical inference in Pisarenkoet al. and propose alternatives. I argue why GPD and GEVD would work for various types of stochastic earthquake processes in time, and not only for the homogeneous (stationary) Poisson process as assumed by Pisarenko et al. (Pure Appl. Geophys 171:1599–1624, 2014). The crucial point of earthquake magnitudes is the poor convergence of their tail distribution to the GPD, and not the earthquake process over time.     

An analytical method for the evaluation of the in-plan irregularity of non-regularly asymmetric buildings

This paper describes a new method for the evaluation of the static eccentricity $e_{s}$ and the ratio $\Omega _{\uptheta } $ of uncoupled torsional to lateral frequencies in real multi-storey buildings. The above-mentioned parameters greatly affect the lateral-to-torsional coupling of the response of asymmetric systems and thus are of paramount importance in the assessment of the in-plan irregularity of buildings. The proposed method, which is a generalization of that suggested by Calderoni et al. (Earthq Spectra 18(2):219–231, 2002), allows the calculation of the static eccentricity $e_{s}$ and the ratio $\Omega _{\uptheta } $ from the structural response to arbitrary distributions of forces and torsional couples. The effectiveness of the method is validated on some regularly and non-regularly asymmetric buildings characterised by different in-plan irregularity. The analyses demonstrate that the results of the method are rigorous in the case of regularly asymmetric systems and only slightly depend upon the heightwise distribution of the forces in the case of non-regularly asymmetric systems. Finally, the values of the static eccentricity $e_{s}$ and the ratio $\Omega _{\uptheta } $ resulting from the proposed method are compared to those obtained by means of the procedure suggested by Makarios and Anastassiadis in (Struct Des Tall Spec Build 7(1):33–55, 1998a; Struct Des Tall Spec Build 7(1):57–71, 1998b) .     

Seismic hazard assessment in the megacity of Blida (Algeria) and its surrounding regions using parametric-historic procedure

The region of Blida is characterized by a relatively high seismic activity, pointed especially during the past two centuries. Indeed, it experienced a significant number of destructive earthquakes such as the earthquakes of March 2, 1825 and January 2, 1867, with intensity of X and IX, respectively. This study aims to investigate potential seismic hazard in Blida city and its surrounding regions. For this purpose, a typical seismic catalog was compiled using historical macroseismic events that occurred over a period of a few hundred years, and the recent instrumental seismicity dating back to 1900. The parametric-historic procedure introduced by Kijko and Graham (1998, 1999) was applied to assess seismic hazard in the study region. It is adapted to deal with incomplete catalogs and does not use any subjective delineation of active seismic zones. Because of the lack of recorded strong motion data, three ground prediction models have been considered, as they seem the most adapted to the seismicity of the study region. Results are presented as peak ground acceleration (PGA) seismic hazard maps, showing expected peak accelerations with 10% probability of exceedance in 50-year period. As the most significant result, hot spot regions with high PGA values are mapped. For example, a PGA of 0.44 g has been found in a small geographical area centered on Blida city.     

Interannual variability of the air–sea CO2 flux in the north Indian Ocean

A simple biogeochemical model coupled to an offline ocean tracer transport model driven by reanalysis ocean data is used to simulate the seasonal and interannual CO $_2$ flux variability in the northern Indian Ocean. The maximum of seasonal and interannual CO $_2$ emission variances in the northern Indian Ocean are located in the coastal Arabian Sea (AS) and Southern Peninsular India (SP) with a basin-wide seasonal amplitude and standard deviation of 0.044 $\pm $ 0.04 Pg C year $^{-1}$ . The area integrated CO $_2$ emissions from these two regions in the model are significantly correlated (above a 95 % level) with the observations of Takahashi et al. (Deep-Sea Res-II, 56:554–577, 2009). The interannual anomalies of CO $_2$ emission from the AS and SP are found as 40 and 30 % of their respective seasonal amplitudes. Both the Arabian Sea (AS) and Southern Peninsular India (SP) interannual CO $_2$ emission anomalies show a 3–4-year variability. The correlations of AS and SP CO $_2$ emission anomalies with the Indian Ocean Dipole/Zonal Mode (IODZM) and Southern Oscillation (SO) indices from 1980 to 1999 are 0.35, 0.21 and 0.32, 0.01 respectively. A 5-year window moving correlation analysis shows that the relationship of AS and SP CO $_2$ emission to the SO and IODZM are complementary to each other. During the years when the correlation of air–sea CO $_2$ emission with the IODZM is stronger, the corresponding correlation with the SO is weaker or opposite. The total change in pCO $_2$ is broken down into changes induced by the individual components such as dissolved inorganic carbon (DIC), sea surface temperature (SST), alkalinity, and salinity and found that (1) the effect of SST in the AS CO $_2$ emission increases (decreases) when the correlation of CO $_2$ emission with the IODZM is positive (negative), and (2) the SP CO $_2$ emission is strongly controlled by the circulation-driven DIC changes; however, this relation is found to be weaker when the SO correlates negatively with the SP CO $_2$ emission.     

An interactive biplot implementation in R for modeling genotype-by-environment interaction

Classical and GGE biplot methods are graphical procedures that allow multivariate data to be analyzed. In particular, the GGE biplot displays the genotype main effect (G) and the genotype by environment interaction (GE) in two-way data. The GGE biplot originates from data graphical analysis of multi-environment trials (MET). Thus, agronomists, crop scientists and geneticists are potential users of this method. However, it can also be used to visualize and analyze other types of data. In this paper, we propose a new interactive computational implementation in R language to perform the main functions of the classical and GGE biplot methods, so it is also useful for MET data visual analysis. This implementation is organized in an R package named GGEBiplotGUI . This package is the only interactive, noncommercial and open source software that currently exists, offering a free alternative to available commercial software. In addition, it can be used without to practically have knowledge of the R programming language. Here, we present and discuss the capabilities and features of the GGEBiplotGUI package and illustrate them by using real data. The GGEBiplotGUI package graphically addresses the questions that a researcher likely asks. This R package is not only a tool for visual data analysis of multi-environment trials, useful for plant breeders and geneticists, in order to study yields from genotypes and interactions between genotype and environment, but also data from other areas can be analyzed by the GGEBiplotGUI package.     

An alternative depth-integrated formulation for granular avalanches over temporally varying topography with small curvature

An alternative formulation is proposed for deriving depth-integrated equations for gravity-driven granular avalanches over a non-trivial topography with small curvature. The coordinate system of Bouchut and Westdickenberg (2004 Bouchut, F and Westdickenberg, M. 2004. Gravity driven shallow water models for arbitrary topography. Commun. Math. Sci., 2: 359–389. [Crossref] , [Google Scholar]) is combined with the unified coordinate (UC) method, so that it can evolve in accordance with the entrainment–deposition processes at the basal surface. The resultant mass and momentum equations are formulated as a conservation system of the Cartesian components of the conservative physical variables. The motion of the flows is driven by the basal topography-induced pressure, pressure gradient, and resisted by the basal friction. The best benefit of this formulation is that it greatly simplifies the computation of the varying coordinate orientations. The features and advantages of this formulation are illustrated by the sliding-mass examples where we simulate the motion of a finite mass of granular material sliding down an inclined chute, running through a transition zone, and being deposited onto a horizontal plane.     

Comment on “Complete Eulerian-mean tracer equation for coarse resolution OGCMs” by M. S. Dubovikov and V. M. Canuto

Dubovikov and Canuto (Dubovikov, M.S. and Canuto, V.M., Complete Eulerian-mean tracer equation for coarse resolution OGCMs. Geophys. Astrophys. Fluid Dyn., 2006, 100, 197–214), after averaging the tracer conservation equation in density coordinates and transforming to height coordinates, concluded that present ocean models are missing important terms in their mean tracer equations. Here we point out some errors made by Dubovikov and Canuto (2006 Dubovikov, MS and Canuto, VM. 2006. Complete Eulerian-mean tracer equation for coarse resolution OGCMs. Geophys. & Astrophys. Fluid Dynam., 100: 197–214. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) in their isopycnal averaging procedure. We draw on the temporal-residual-mean (TRM) theory and show that when the isopycnal averaging and coordinate transformation are performed correctly, the tracer equations of present ocean circulation models are recovered. We therefore conclude that present ocean circulation models are not neglecting the leading order terms identified by Dubovikov and Canuto (2006 Dubovikov, MS and Canuto, VM. 2006. Complete Eulerian-mean tracer equation for coarse resolution OGCMs. Geophys. & Astrophys. Fluid Dynam., 100: 197–214. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]).     

New Seismic Noise Models Obtained Using Very Broadband Stations

It has been two decades since the last comprehensive standard model of ambient earth noise was published Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The PETERSON model was updated by analyzing the absolute quietest conditions for stations within the GSN (Berger et al. in J Geophys Res 109, 2005; Mcnamara and Buland in Bull Seism Soc Am 94:1517–1527, 2004; Ringler et al. in Seismol Res Lett 81(4) doi:10.1785/gssrl.81.4.605, 2010). Unfortunately, both the original model and the updated models did not include any deployed station in North Africa and Middle East, which reflects the noise levels within the desert environment of those regions. In this study, a survey was conducted to create a new seismic noise model from very broadband stations which recently deployed in North Africa. For this purpose, 1 year of continuous recording of seismic noise data of the Egyptian National Seismic Network (ENSN) was analyzed in order to create a new noise model. Seasonal and diurnal variations in noise spectra were recorded at each station. Moreover, we constructed a new noise model for each individual station. Finally, we obtained a new cumulative noise model for all the stations. We compared the new high-noise model (EHNM) and new low-noise model (ELNM) with both the high-noise model (NHNM) and low-noise model (NLNM) of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The obtained noise levels are considerably lower than low-noise model of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993) at ultra long period band (ULP band), but they are still below the high-noise model of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The results of this study could be considered as a first step to create permanent seismic noise models for North Africa and Middle East regions.     

Simple laminar dynamos: from two rolls to one

The study of laminar flows with simple streamlines in conducting fluids that might support kinematic dynamo action has a long history. Early successful examples include the Gailitis (1970 Gailitis, A. 1970. Self-excitation of a magnetic field by a pair of annular vortices. Mag. Gidrod, 6: 19–22. (English translation: Magnetohydrodynamics, 6:14–17) [Google Scholar]) dynamo, which has two meridional rolls. We show here that a single meridional roll, considered as a limit of a modified Gailitis-like dynamo, can also support kinematic dynamo action in a deep spherical shell. This flow has a claim to be the simplest such flow possible.     

Empirical ground-motion relations using moderate earthquakes recorded by Medellín–Aburrá Valley (Colombia) strong-motion networks

We tested attenuation relations obtained for different regions of the world to verify their suitability to predict strong-motion data recorded by Medellín and Aburrá Valley Accelerographic Networks. We used as comparison criteria, the average of the difference between the observed and the predicted data as a function of epicenter distance and its standard deviation. We also used the approach developed by Sherbaum et al. (Bull Seism Soc Am 94:2164–2185, 2004) that provides a method to evaluate the overall goodness-of-fit of ground-motion prediction equations. The predictive models selected use a generic focal depth. We found that this parameter has an important influence in the ground-motion predictions and must be taken into account as an independent variable. We also found important to characterize the local soil amplification to improve the attenuation relations. We found empirical relations for peak horizontal acceleration PGA and velocity PGV based on the Kamiyama and Yanagisawa (Soils Found 26:16–32, 1986) approach. $$\begin{aligned} \log _{10} (PGA)=0.5886M_L -1.0902\log _{10}(R)-0.0035H+C_{st}\pm 0.\text{29} \end{aligned}$$ $$\begin{aligned} \log _{10} (PGV)=0.7255M_L -1.8812\log _{10}(R)-0.0016H+C_{st}\pm 0.36 \end{aligned}$$ where PGA is measured in cm/s $^{2}$ and PGV in cm/s, $M_{L}$ is local magnitude in the range 2.8–6.5, $R$ is epicentral distance up to 290 km, $H$ is focal depth in km and $C_{st}$ is a coefficient that accounts for the site response due to soil conditions of each recording station. The introduction of focal depth and local site conditions as independent variables, minimize the residuals and the dispersion of the predicted data. We conclude that $H$ and $C_{st}$ are sensitive parameters, having a strong influence on the strong-motion predictions. Using the same functional form, we also propose an empirical relation for the root mean square acceleration a $_\mathrm{rms}$ : $$\begin{aligned} \log _{10} \left( {a_{rms} } \right)=0.4797M_L -1.1665\log _{10} (R)-0.00201H+C_{st}\pm 0.40 \end{aligned}$$ where a $_\mathrm{rms}$ is measured in cm/s $^{2}$ , from the S-wave arrival and using a window length equal to the rupture duration. The other variables are the same as those for PGA and PGV. The site correction coefficients $C_{st}$ found for PGA, PGV and a $_\mathrm{rms}$ show a similar trend indicating a good correlation with the soil conditions of the recording sites.     

Convection in rotating annular channels heated from below. Part 2. Transitions from steady flow to turbulence

We report the results of fully three-dimensional numerical simulations of nonlinear convection in a Boussinesq fluid in an annular channel rotating about a vertical axis with lateral no-slip or stress-free sidewalls, stress-free top and bottom, uniformly heated from below, a problem first studied by Davies-Jones and Gilman (1971 Davies-Jones, RP and Gilman, PA. 1971. Convection in a rotating annulus uniformly heated from below.. J. Fluid Mech., 46: 65–81.  [Google Scholar]) and Gilman (1973 Gilman, PA. 1973. Convection in a rotating annulus uniformly heated from below. Part 2. Nonlinear results. J. Fluid Mech., 57: 381–400.  [Google Scholar]). A substantial range of the Rayleigh number R (R_c≤R≤O(100 R_c)), where R_c denotes the critical value at the onset of convection) is considered. It is found that the wall-localized convection mode, unaffected by the velocity boundary condition imposed on the sidewalls, is nonlinearly robust. Both directions of travelling waves, one propagating against the sense of rotation near the outer sidewall and the other propagating in the same sense as the rotation in the vicinity of the inner sidewall, are always present in the nonlinear solutions. In contrast to nonlinear convection in a rotating Bénard layer, neither convection rolls nor the Küpper–Lortz instability can exist in a rotating annular channel because of the effect of the sidewalls. It is the nonlinear interaction between the wall-localized modes and the internal mode that plays an essential role in determining the nonlinear properties of convection in a rotating annular channel. Our studies reveal systematically the various nonlinear phenomena, from steady travelling waves trapped in the vicinities of the sidewalls to convective turbulence exhibiting columnar structure.     

Publications received by the Editor

Abstract

Temperature profiles (temperature as function of depth) can be used to derive vertical flow velocities or recharge rates. In many cases, solutions to the one-dimensional (1-D) heat transport equation are used, considering steady-state boundary conditions. Factors which can influence the derivation of the mean vertical flow velocity are studied here. Therefore, an explicit finite-difference approximation to the 1-D heat transport equation coupled with an inverse scheme was used to interpret temperature profiles. Measurement error (larger than 0.05°C) can result in important deviation of the derived mean flow velocity. Variation of vertical flow velocity as a function of time leads to asymmetric temperature envelopes. Yearly variation in vertical flow velocities, or temperature variations of the recharge water, also results in asymmetric temperature envelopes. Interpretation of these asymmetric envelopes leads to important differences between derived and actual mean vertical flow velocities.

Citation Vandenbohede, A. &; Lebbe, L. (2010) Vandenbohede, A. and Lebbe, L. 2010. Parameter estimation based on vertical heat transport in the surficial zone. Hydrogeol. J., 18 , 931–943,doi:10.1007/s10040-009–0557–5 [Google Scholar] Recharge assessment by means of vertical temperature profiles: analysis of possible influences. Hydrol. Sci. J. 55(5), 792–804.     

Integrated approach to choosing suitable areas for the realization of productive wells in rural areas of sub-Saharan Africa

Abstract

Comprehensive geological, hydrogeological and hydrogeochemical investigations were carried out in the south of Hodh El Chargui (southeast Mauritania). Obtaining a hydrogeological conceptual model is crucial for groundwater resources development and management. This is especially true in developing countries and in the rural areas of sub-Saharan Africa. The chosen areas are represented by lithologies referring to a long geological time frame, dating from the Upper Neoproterozoic to the Quaternary age. We developed a methodology in order to identify sites that were suitable for the realization of productive, protected and correct wells to supply safe water to the rural community. A multicriteria approach to studying hydrogeology was used in the project area. In order to identify some main areas in which to carry out pilot interventions, criteria relating to water accessibility and availability, and to hydrogeological and water quality, were considered. Moreover, during the project, it was possible to transfer know-how and hand over responsibilities to the local population and bodies.

Citation Ghiglieri, G. & Carletti, A. (2010 Ghiglieri, G., Balia, R., Oggiano, G. and Pittalis, D. 2010. Prospecting for safe (low fluoride) groundwater in the Eastern African Rift: the Arumeru District (Northern Tanzania). Hydrol. Earth System Sci., 14: 1081–1091.  [Google Scholar]) Integrated approach to choosing suitable areas for the realization of productive wells in rural areas of sub-Saharan Africa. Hydrol. Sci. J. 55(8), 1357–1370.     

Intensity predictive attenuation models calibrated in Mw for metropolitan France

In the framework of the SIGMA project, a study was launched to develop a parametric earthquake catalog for the historical period, covering the metropolitan territory and calibrated in Mw. A set of candidate calibration events was selected corresponding to earthquakes felt over a part of the French metropolitan territory, which are fairly well documented both in terms of macroseismic intensity distributions (SisFrance BRGM-EDF-IRSN) and magnitude estimates. The detailed analysis of the macroseismic data led us to retain only 30 events out of 65 with Mw ranging from 3.6 to 5.8. In order to supplement the dataset with data from larger magnitude events, Italian earthquakes were also considered (11 events posterior to 1900 with Mw?≥?6.0 out of 15 in total), using both the DBMI11 macroseismic database (Locati et al. in Seismol Resour Lett 85(3):727–734, 2014) and the parametric information from the CPTI11 (Rovida et al. in CPTI11, la versione 2011 del Catalogo Parametrico dei Terremoti Italiani Istituto Nazionale di Geofisica et Vulcanologia, Milano, Bologna, 2011.  https://doi.org/10.6092/ingv.it-cpti11). To avoid introducing bias related to the differences in terms of intensity scales (MSK vs. MCS), only intensities smaller than or equal to VII were considered (Traversa et al. in On the use of cross-border macroseismic data to improve the estimation of past earthquakes seismological parameters, 2014). Mw and depth metadata were defined according to the Si-Hex catalogue (Cara et al. in Bull Soc Géol Fr 186:3–19, 2015.  https://doi.org/10.2113/qssqfbull.186.1.3), published information, and to the specific worked conducted within SIGMA related to early instrumental recordings (Benjumea et al. in Study of instrumented earthquakes that occurred during the first part of the 20th century (1905–1962), 2015). For the depth estimates, we also performed a macroseismic analysis to evaluate the range of plausible estimates and check the consistency of the solutions. Uncertainties on the metadata related to the calibration earthquakes were evaluated using the range of available alternative estimates. The intensity attenuation models were developed using a one-step maximum likelihood scheme. Several mathematical formulations and sub-datasets were considered to evaluate the robustness of the results (similarly to Baumont and Scotti in Accounting for data and modeling uncertainties in empirical macroseismic predictive equations (EMPEs). Towards “European” EMPEs based on SISFRANCE, DBMI, ECOS macroseismic database, 2008). In particular, as the region of interest may be characterized by significant laterally varying attenuation properties (Bakun and Scotti in Geophys J Int 164:596–610, 2006; Gasperini in Bull Seismol Soc Am 91:826–841, 2001), we introduced regional attenuation terms to account for this variability. Two zonation schemes were tested, one at the national scale (France/Italy), another at the regional scale based on the studies of Mayor et al. (Bull Earthq Eng, 2017.  https://doi.org/10.1007/s10518-017-0124-8) for France and Gasperini (2001) for Italy. Between and within event residuals were analyzed in detail to identify the best models, that is, the ones associated with the best misfit and most limited residual trends with intensity and distance. This analysis led us to select four sets of models for which no significant trend in the between- and within-event residuals is detected. These models are considered to be valid over a wide range of Mw covering?~?3.5–7.0.     

On the selection of GMPEs for Vrancea subcrustal seismic source

The Vrancea subcrustal seismic source is characterized by large magnitude ( $M_{W} \ge 7$ ) intermediate-depth earthquakes that occur two or three times during a century on average. In this study several procedures are used to grade four candidate ground motion prediction equations proposed for Vrancea source in the SHARE project. In the work of Delavaud et al. (J Seismol 16(3):451–473, 2012) four ground motion prediction models developed for subduction zones (Zhao et al. in Bull Seism Soc Am 96(3):898–913, 2006; Atkinson and Boore in Bull Seism Soc Am 93(4):1703–1729, 2003; Youngs et al. in Seism Res Lett 68(1):58–73, 1997; Lin and Lee in Bull Seism Soc Am 98(1):220–240, 2008) are suggested as suitable for Vrancea subcrustal seismic source. The paper presents the appropriateness analysis of the four suggested ground motion prediction equations done using a dataset of 109 triaxial accelerograms recorded during seven Vrancea seismic events with moment magnitude $M_{W}$ between 5.4 and 7.4, occurred in the past 35 years. The strong ground motions were recorded in Romania, as well as in Bulgaria, Republic of Moldova and Serbia. Based on the ground motion dataset several goodness-of-fit measures are used in order to quantify how well the selected models match with the recorded data. The compatibility of the four ground motion prediction models with respect to magnitude scaling and distance scaling implied by strong ground motion dataset is investigated as well. The analyses show that the Youngs et al. (Seism Res Lett 68(1):58–73, 1997) and Zhao et al. (Bull Seism Soc Am 96(3):898–913, 2006) ground motion prediction models have a better fit with the data and can be candidate models for Probabilistic Seismic Hazard Assessment.     

Streamflow dynamics in a Mediterranean temporary river

Abstract

Mediterranean rivers are characterized by the irregularity of flow, harsh hydrological fluctuations and a profound transformation as the result of human activity. In this study, we investigate the streamflow response of a Mediterranean temporary river in which different groundwater, agriculture and urban contributions play an important role. Streamflow was measured at three nested gauging stations installed along Na Borges River, a lowland agricultural catchment (319 km²) on the island of Mallorca. Based on two hydrological years (2004/05 and 2005/06), potential evapotranspiration and surface water–groundwater interaction generated a succession of four different hydrological periods playing an important role in baseflow dynamics. The runoff coefficients were very low (<5%). At the event scale, groundwater also controlled runoff response, being very different according to hydrogeology, antecedent conditions and human impacts. During dry seasons, wastewater and karstic spring discharges maintain an influent regime into some streams. As a result, intense rainstorms in late summer generated water volumes over the impervious urban surfaces involved, with the result that quickflow was significant because the hydrological pathways were active. Citation Estrany, J., Garcia, C. & Alberich, R. (2010 Estrany, J., Garcia, C. and Batalla, R. J. 2010. Hydrological response of a small Mediterranean agricultural catchment. J. Hydrol., 380(1-2): 180–190. doi:10.1016/j.jhydrol.2009.10.035[Crossref] , [Google Scholar]) Streamflow dynamics in a Mediterranean temporary river. Hydrol. Sci. J. 55(5), 717–736.      

Proposal of an empirical site classification method based on target simulated horizontal over vertical spectral ratio

Nowadays, most of the site classifications schemes are based on the predominant period of the site as determined from the average horizontal to vertical spectral ratios of seismic motion or microtremor. However, the difficulty lies in the identification of the predominant period in particular if the observed average response spectral ratio does not present a clear peak but rather a broadband amplification or multiple peaks. In this work, based on the Eurocode-8 (2004) site classification, and assuming bounded random fields for both shear and compression waves-velocities, damping coefficient, natural period and depth of soil profile, one propose a new site-classification approach, based on “target” simulated average \( H/V \) spectral ratios, defined for each soil class. Taking advantage of the relationship of Kawase et al. (Bull Seismol Soc Am 101:2001–2014, 2011), which link the \( H/V \) spectral ratio to the horizontal (\( HTF \)) over the vertical (\( VTF \)) transfer functions, statistics of \( H/V \) spectral ratio via deterministic visco-elastic seismic analysis using the wave propagation theory are computed for the 4 soil classes. The obtained results show that \( H/V \) and \( HTF \) have amplitudes and shapes remarkably different among the four soil classes and exhibit fundamental peaks in the period ranges remarkably similar. Moreover, the “target” simulated average \( H/V \) spectral ratios for the 4 soil classes are in good agreement with the experimental ones obtained by Zhao et al. (Bull Seismol Soc Am 96:914–925, 2006) from the abundant and reliable Japanese strong motions database Kik-net, Ghasemi et al. (Soil Dyn Earthq Eng 29:121–132, 2009) from the Iranian strong motion data, and Di Alessandro et al. (Bull Sesismol Soc Am 106:2, 2011.  https://doi.org/10.1785/0120110084) from the Italian strong motion data. In addition to the 4 EC-8 standard soil classes (A, B, C and D), the superposition of the 4 target \( H/V \) reveals 3 new boundary site classes; AB, BC and CD, for overlapping \( V_{s,30} \) ranges when the predominant peak is not clearly consistent with any of the 4 proposed classes. Finally, one proposes a site classification index based on the ratio between the cross-correlation and the mean quadratic error between the in situ \( H/V \) spectral ratio and the “target” one. In order to test the reliability of the proposed approach, data from 139 sites were used, 132 collected from the Kik-net network database from Japan and 7 from Algeria. The site classification success rate per site class are around 93, 82, 89 and 100% for rock, hard soil, medium soil and soft soil, respectively. Zhao et al. (2006) found an average success for the 4 classes of soil close to 60%, similar to what one found in the present study (63%) without considering the new soil classes, but much smaller if one considers them (86%). In the absence of \( V_{s,30} \) data, the proposed approach can be an alternative to site classification.     

Effects of Location Errors in Pattern Informatics

The effect of location errors in the performance of seismicity-based forecasting methods was studied here using one particular binary forecast technique, the Pattern Informatics (PI) technique (Rundle et al., Proc Nat Acad Sci USA 99, 2514–2521, 2002; Tiampo et al., Pure Appl Geophys 159, 2429–2467, 2002). The Southern Californian dataset was used to generate a series of perturbed catalogs by adding different levels of noise to epicenter locations. The PI technique was applied to these perturbed datasets to perform retrospective forecasts that were evaluated by means of skill scores, commonly used in atmospheric sciences. These results were then compared to the effectiveness obtained from the original dataset. Isolated instances of decline of the PI performance were observed due to the nature of the skill scores themselves, but no clear trend of degradation was identified. Dependence on the total number of events in a catalog also was studied, with no systematic degradation in the performance of the PI for catalogs with events in the cases studied. These results suggest that the stability of the PI method is due to the invariance of the clustering patterns identified by the TM metric (Thirumalai and Mountain, Phys Rev A 39, 3563–3573, 1989) when applied to seismicity.     

Catchment-scale vulnerability assessment of groundwater pollution from diffuse sources using the DRASTIC method: a case study

Abstract

The catchment-scale groundwater vulnerability assessment that delineates zones representing different levels of groundwater susceptibility to contaminants from diffuse agricultural sources has become an important element in groundwater pollution prevention for the implementation of the EU Water Framework Directive (WFD). This paper evaluates the DRASTIC method using an ArcGIS platform for assessing groundwater vulnerability in the Upper Bann catchment, Northern Ireland. Groundwater vulnerability maps of both general pollutants and pesticides in the study area were generated by using data on the factors depth to water, net recharge, aquifer media, soil media, topography, impact of vadose zone, and hydraulic conductivity, as defined in DRASTIC. The mountain areas in the study area have “high” (in 4.5% of the study area) or “moderate” (in 25.5%) vulnerability for general pollutants due to high rainfall, net recharge and soil permeability. However, by considering the diffuse agricultural sources, the mountain areas are actually at low groundwater pollution risk. The results of overlaying the maps of land use and the groundwater vulnerability are closer to the reality. This study shows that the DRASTIC method is helpful for guiding the prevention practices of groundwater pollution at the catchment scale in the UK.

Citation Yang, Y. S. & Wang, L. (2010 Yang, Y. S. and Wang, L. 2010. A review of modelling tools for implementation of the EU Water Framework Directive in handling diffuse water pollution. Water Resour. Manage., 24: 1819–1843.  [Google Scholar]) Catchment scale vulnerability assessment of groundwater pollution from diffuse sources using the DRASTIC method: a case study. Hydrol. Sci. J. 55(7), 1206–1216.     

\tau_p^{{\rm max}} magnitude estimation, the case of the April 6, 2009 L’Aquila earthquake

Rapid magnitude estimate procedures represent a crucial part of proposed earthquake early warning systems. Most of these estimates are focused on the first part of the P-wave train, the earlier and less destructive part of the ground motion that follows an earthquake. Allen and Kanamori (Science 300:786–789, 2003) proposed to use the predominant period of the P-wave to determine the magnitude of a large earthquake at local distance and Olivieri et al. (Bull Seismol Soc Am 185:74–81, 2008) calibrated a specific relation for the Italian region. The Mw 6.3 earthquake hit Central Italy on April 6, 2009 and the largest aftershocks provide a useful dataset to validate the proposed relation and discuss the risks connected to the extrapolation of magnitude relations with a poor dataset of large earthquake waveforms. A large discrepancy between local magnitude (ML) estimated by means of $\tau_p^{{\rm max}}$ evaluation and standard ML (6.8 ± 1.5 vs. 5.9 ± 0.4) suggests using caution when ML vs. $\tau_p^{{\rm max}}$ calibrations do not include a relevant dataset of large earthquakes. Effects from large residuals could be mitigated or removed introducing selection rules on τ _p function, by regionalizing the ML vs. $\tau_p^{{\rm max}}$ function in the presence of significant tectonic or geological heterogeneity, and using probabilistic and evolutionary methods.