Modeling Seismic Attributes of Pn Waves using the Spectral-Element Method

To investigate the nature of Pn propagation, we have implemented the spectral-element method (SEM) for vertically and laterally varying media with and without attenuation. As a practical measure, essential features of the Pn waves are distilled into seismic attributes including arrival times, amplitudes and pulse frequencies. To validate the SEM simulations, we first compare the SEM results with reflectivity calculations of Braile and Smith (Geophys. J.R. Astr. Soc. 40, 145–176, 1975) and then to the asymptotic results of ?erveny and Ravindra (Theory of Seismic Headwaves, University of Toronto Press, pp. 235–250, 1971). Models with random, laterally varying Moho structures are then simulated, where the amplitude and pulse frequency characteristics are found to be stable to small Moho interface perturbations. SEM calculations for models with different upper-mantle velocity gradients are next performed where it is found that interference effects can strongly influence the Pn amplitudes and pulses frequencies. For larger-scale, laterally varying structures, SEM models similar to that found along the Hi-CLIMB array in Tibet are then performed. It is observed that large-scale structures, along with small-scale structures, upper-mantle velocity gradients and attenuation, can all significantly affect the Pn attributes. Ambiguities between upper-mantle velocity gradients and attenuation are also found when using Pn amplitudes and pulse frequency attributes. These ambiguities may be resolved, to some degree, by using the curvature of the travel times at longer regional distance, however, this would also be complicated by lateral variability.     

Speciation and Determination of Inorganic Mercury and Methylmercury by Headspace Single Drop Microextraction and Electrothermal Atomic Absorption Spectrometry in Water and Fish

In this study, headspace single drop microextraction (HS‐SDME) method in combination with electrothermal atomic absorption spectrometry (ETAAS) method was developed and validated for the speciation and determination of inorganic mercury (iHg) and methylmercury (MeHg). MeHg and iHg species were reduced to volatile methylmercury hydride (CH₃HgH) and elemental mercury, respectively, in the presence of NaBH₄ and trapped onto a drop of acceptor phase in the tip of a microsyringe. Thiourea and ammonium pyrrolydinedithiocarbamate (APDC) were tested as the acceptor phase. The experimental parameters of the method such as microextraction time, temperature, NaBH₄ concentration, acceptor phase concentration, and pH of the medium were investigated to obtain distinctive conditions for mercury species. Possible interference effects have also been investigated. In order to validation of the method, analytical figures of merits such as accuracy, precision, limit of detection (LOD), limit of quantitation (LOQ), and linear working range have been evaluated. Accuracy of the method has been verified by analyzing certified reference materials (BCR 453 Tuna fish) and spiked samples. The proposed method was applied for the speciation and determination of mercury species in water and fish samples. Mercury species (MeHg and iHg) have been determined in the real samples with a relative error less than 10%.     

Spatial distribution and temporal variation of reference evapotranspiration in arid and semi‐arid regions of Iran

Analysis of spatial and temporal variations of reference evapotranspiration (ET_o) is important in arid and semi‐arid regions where water resources are limited. The main aim of this study was to analyse the spatial distribution and the annual, seasonal and monthly trends of the Penman–Monteith ET_o for 21 stations in the arid and semi‐arid regions of Iran. Three statistical tests the Mann‐Kendall, Sen's slope estimator and linear regression were used for the analysis. The analysis revealed that ET_o increased from January to July and deceased from July to December at almost all stations. Additionally, higher annual ET_o values were found in the southeast of the study region and lower values in the northwest of the region. Although the results showed both positive and negative trends in annual ET_o series, ET_o generally increased, significantly so in six (～30%) of the stations. Analysis of the impacts of meteorological variables on the temporal trends of ET_o indicated that the increasing trend of ET_o was most likely due to a significant increase in minimum air temperature, while decreasing trend of ET_o was mainly caused by a significant decrease in wind speed. At the sites where increasing ET_o trends were statistically significant, the rate of increase varied from (+)8·36 mm/year at Mashhad station to (+)31·68 mm/year at Iranshahr station. On average, an increasing trend of (+)4·42 mm/year was obtained for the whole study area during the last four decades. Seasonal and monthly ET_o have also tended to increase at the majority of the stations. The greatest numbers of significant trends were observed in winter on the seasonal time‐scale and in September on the monthly time‐scale. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling of groundwater flow in heterogeneous porous media by finite element method

The HySuf‐FEM code (Hydrodynamic of Subsurface Flow by Finite Element Method) is proposed in this article in order to estimate the spatial variability of the transmissivity values of the Berrechid aquifer (Morocco). The calibration of the model is based on the hydraulic head, hydraulic conductivity and recharge. Three numerical tests are used to validate the model and verify its convergence. The first test case consists in using the steady analytical solution of the Poisson equation. In the second, the model has been compared with the hydrogeological system which is characterized by an unconfined monolayer (isotropic layer) and computed by using PMWIN‐MODFLOW software. The third test case is based on the comparison between the results of HySuf‐FEM and the multiple cell balance method in the aquifer system with natural boundaries case. Good agreement between the Hydrodynamic of Subsurface Flow, the numerical tests and the spatial distribution of the thickening of the hydrogeological system is deduced from the analysis and the interpretations of hydrogeological wells. Copyright © 2011 John Wiley & Sons, Ltd.     

Quantitative Analysis and Efficiency Assessment of Floodwater Harvesting System in Arid Region: Case of Touzouz ephemeral stream—Mzab Valley

Water Resources - The Mzab Valley is characterized by an arid climate with limited amounts of rainfall. In the context of water scarcity, the Mzab water management system recovers floodwater for...     

A new integration scheme for application to seismic hybrid simulation

Hybrid simulation is a powerful test method for evaluating the seismic performance of structural systems. This method makes it feasible that only critical components of a structure be experimentally tested. This paper presents a newly proposed integration algorithm for seismic hybrid simulation which is aimed to extend its capabilities to a wide range of systems where existing methods encounter some limitations. In the proposed method, which is termed the variable time step (VTS) integration method, an implicit scheme is employed for hybrid simulation by eliminating the iterative phase on experimental element, the phase which is necessary in regular implicit applications. In order to study the effectiveness of the VTS method, a series of numerical investigations are conducted which show the successfulness of the VTS method in obtaining accurate, stable and converged responses. Then, in a comparative approach, the improved accuracy of the VTS method over commonly used integration methods is demonstrated. The stability of the VTS method is also studied and the results show that it provides conditional stability; however, its stability limit is well beyond the accuracy limit. The effect of time delay on the VTS method results is also investigated and it is shown that the VTS method is quite successful in handling this experimental error.     

Comparative study of different wavelets for developing parsimonious Volterra model for rainfall-runoff simulation

Although the Volterra models are non-parsimonious ones, they are being used because they can mimic dynamics of complex systems. However, applying and identification of the Volterra models using data may result in overfitting problem and uncertainty. In this investigation we evaluate capability of different wavelet forms for decomposing and compressing the Volterra kernels in order to overcome this problem by reducing the number of the model coefficients to be estimated and generating smooth kernels. A simulation study on a rainfall?runoff process over the Cache River watershed showed that the method performance is successful due to multi-resolution capacity of the wavelet analysis and high capability of the Volterra model. The results also revealed that db2 and sym2 wavelets have the same high potential in improving the linear Volterra model performance. However, QS wavelet was more successful in yielding smooth kernels. Moreover, the probability of overfitting while identifying the nonlinear Volterra model may be less than the linear model.     

Integration of multiple soft data sets in MPS thru multinomial logistic regression: a case study of gas hydrates

A new approach is described to allow conditioning to both hard data (HD) and soft data for a patch- and distance-based multiple-point geostatistical simulation. The multinomial logistic regression is used to quantify the link between HD and soft data. The soft data is converted by the logistic regression classifier into as many probability fields as there are categories. The local category proportions are used and compared to the average category probabilities within the patch. The conditioning to HD is obtained using alternative training images and by imposing large relative weights to HD. The conditioning to soft data is obtained by measuring the probability–proportion patch distance. Both 2D and 3D cases are considered. Synthetic cases show that a stationary TI can generate non-stationary realizations reproducing the HD, keeping the texture indicated by the TI and following the trends identified in probability maps obtained from soft data. A real case study, the Mallik methane-hydrate field, shows perfect reproduction of HD while keeping a good reproduction of the TI texture and probability trends.     

Probabilistic analysis of soil‐structure interaction effects on the seismic performance of structures

This paper revisits the phenomenon of dynamic soil‐structure interaction (SSI) with a probabilistic approach. For this purpose, a twofold objective is pursued. First, the effect of SSI on inelastic response of the structure is studied considering the prevailing uncertainties. Second, the consequence of practicing SSI provisions of the current seismic design codes on the structural performance is investigated in a probabilistic framework. The soil‐structure system is modeled by the sub‐structure method. The uncertainty in the properties of the soil and the structure is described by random variables that are input to this model. Monte Carlo sampling analysis is employed to compute the probability distribution of the ductility demand of the structure, which is selected as the metrics for the structural performance. In each sample, a randomly generated soil‐structure system is subjected to a randomly selected and scaled ground motion. To comprehensively model the uncertainty in the ground motion, a suite of 3269 records is employed. An extensive parametric study is conducted to cover a wide range of soil‐structure systems. The results reveal the probability that SSI increases the ductility demand of structures designed based on the conventional fixed‐based assumption but built on flexible soil in reality. The results also show it is highly probable that practicing SSI provisions of modern seismic codes increase the ductility demand of the structure. Copyright © 2016 John Wiley & Sons, Ltd.