Prediction of longitudinal dispersion coefficient in natural rivers using a cluster-based Bayesian network

The longitudinal dispersion coefficient is a key element in determining the distribution and transmission of pollution, especially when cross-sectional mixing is completed. However, the existing predictive techniques for this purpose exhibit great amounts of uncertainty. The main objective of this study is to present a more accurate model for predicting longitudinal dispersion coefficient in natural rivers and streams. Bayesian network (BN) approach was considered in the modeling procedure. Two forms of input variables including dimensional and dimensionless parameters were examined to find the best model structure. In order to increase the performance of the model, the clustering method as a preprocessing data technique was applied to categorize the data in separate groups with similar characteristics. An expansive data set consisting of 149 field measurements was used for training and testing steps of the developed models. Three performance evaluation criteria were adopted for comparison of the results of the different models. Comparison of the present results with the artificial neural network (ANN) model and also well-known existing equations showed the efficiency of the present model. The performance of dimensionless BN model 30% is more than dimensional ones in terms of the root mean square error. The accuracy criterion was increased from 70 to 83% by performing clustering analysis on the BN model. The BN-cluster model 43% is more accurate than ANN model in terms of the accuracy criterion. The results indicate that the BN-cluster model give 16% better results than the best available considered model in terms of the accuracy criterion. The developed model provides a suitable approach for predicting pollutant transport in natural rivers.     

Applying different scenarios for landslide spatial modeling using computational intelligence methods

Landslides every year impose extensive damages to human beings in various parts of the world; therefore, identifying prone areas to landslides for preventive measures is essential. The main purpose of this research is applying different scenarios for landslide susceptibility mapping by means of combination of bivariate statistical (frequency ratio) and computational intelligence methods (random forest and support vector machine) in landslide polygon and point formats. For this purpose, in the first step, a total of 294 landslide locations were determined from various sources such as aerial photographs, satellite images, and field surveys. Landslide inventory was randomly split into a testing dataset 70% (206 landslide locations) for training the different scenarios, and the remaining 30% (88 landslides locations) was used for validation purposes. To providing landslide susceptibility maps, 13 conditioning factors including altitude, slope angle, plan curvature, slope aspect, topographic wetness index, lithology, land use/land cover, distance from rivers, drainage density, distance from fault, distance from roads, convergence index, and annual rainfall are used. Tolerance and the variance inflation factor indices were used for considering multi-collinearity of conditioning factors. Results indicated that the smallest tolerance and highest variance inflation factor were 0.31 and 3.20, respectively. Subsequently, spatial relationship between classes of each landslide conditioning factor and landslides was obtained by frequency ratio (FR) model. Also, importance of the mentioned factors was obtained by random forest (RF) as a machine learning technique. The results showed that according to mean decrease accuracy, factors of altitude, aspect, drainage density, and distance from rivers had the greatest effect on the occurrence of landslide in the study area. Finally, the landslide susceptibility maps were produced by ten scenarios according to different ensembles. The receiver operating characteristics, including the area under the curve (AUC), were used to assess the accuracy of the models. Results of validation of scenarios showed that AUC was varying from 0.668 to 0.749. Also, FR and seed cell area index indicators show a high correlation between the susceptibility classes with the landslide pixels and field observations in all scenarios except scenarios 10_RF and 10_SVM. The results of this study can be used for landslides management and mitigation and development activities such as construction of settlements and infrastructure in the future.     

Effect of silica fume on compressive strength of oil-polluted concrete in different marine environments

In the present research, effect of silica fume as an additive and oil polluted sands as aggregates on compressive strength of concrete were investigated experimentally. The amount of oil in the designed mixtures was assumed to be constant and equal to 2% of the sand weight. Silica fume accounting for 10%, 15% and 20% of the weight is added to the designed mixture. After preparation and curing, concrete specimens were placed into the three different conditions: fresh, brackish and saltwater environments (submerged in fresh water, alternation of exposed in air & submerged in sea water and submerged in sea water). The result of compressive strength tests shows that the compressive strength of the specimens consisting of silica fume increases significantly in comparison with the control specimens in all three environments. The compressive strength of the concrete with 15% silica fume content was about 30% to 50% higher than that of control specimens in all tested environments under the condition of using polluted aggregates in the designed mixture.     

Modeling wave propagation by using normalized input-output minimization (NIOM)

A new method for modeling wave propagation is described here, its application is discussed, and its results are compared with those obtained using the conventional correlation and unit impulse response methods. The method uses spectral analysis by minimizing the mean square values of the system input and output when subjected to a constraint, and is effective in detecting arrival times of incident and reflected waves and in revealing their relative amplitudes as well. The method is applied to simple models and to the strong motion records obtained at the TTRL (Koto-ku, Minamisuna) and Chiba vertical arrays during three earthquakes in Japan. The travel times evaluated by the NIOM method agree with the results obtained by the geophysical measurements of S-wave velocity at the sites. The method is also effective in showing the amplification property of shallow layers at the TTRL and Chiba sites.     

Aeolian sediment fingerprinting using a Bayesian mixing model

Identifying sand provenance in depositional aeolian environments (e.g. dunefields) can elucidate sediment pathways and fluxes, and inform potential land management strategies where windblown sand and dust is a hazard to health and infrastructure. However, the complexity of these pathways typically makes this a challenging proposition, and uncertainties on the composition of mixed‐source sediments are often not reported. This study demonstrates that a quantitative fingerprinting method within the Bayesian Markov Chain Monte Carlo (MCMC) framework offers great potential for exploring the provenance and uncertainties associated with aeolian sands. Eight samples were taken from dunes of the small (~58 km²) Ashkzar erg, central Iran, and 49 from three distinct potential sediment sources in the surrounding area. These were analyzed for 61 tracers including 53 geochemical elements (trace, major and rare earth elements (REE)) and eight REE ratios. Kruskal–Wallis H‐tests and stepwise discriminant function analysis (DFA) allowed the identification of an optimum composite fingerprint based on six tracers (Rb, Sr, ⁸⁷Sr, (La/Yb)_n, Ga and δCe), and a Bayesian mixing model was applied to derive the source apportionment estimates within an uncertainty framework. There is substantial variation in the uncertainties in the fingerprinting results, with some samples yielding clear discrimination of components, and some with less clear fingerprints. Quaternary terraces and fans contribute the largest component to the dunes, but they are also the most extensive surrounding unit; clay flats and marls, however, contribute out of proportion to their small outcrop extent. The successful application of these methods to aeolian sediment deposits demonstrates their potential for providing quantitative estimates of aeolian sediment provenances in other mixed‐source arid settings, and may prove especially beneficial where sediment is derived from multiple sources, or where other methods of provenance (e.g. detrital zircon U–Pb dating) are not possible due to mineralogical constraints. Copyright © 2017 John Wiley & Sons, Ltd.     

Comparison of Bayesian estimates of peak ground acceleration (<Emphasis Type="Italic">A</Emphasis>max) with PSHA in Iran

Bayesian probability theory is an appropriate and useful method for estimating parameters in seismic hazard analysis. The analysis in Bayesian approaches is based on a posterior belief, also their special ability is to take into account the uncertainty of parameters in probabilistic relations and a priori knowledge. In this study, we benefited the Bayesian approach in order to estimate maximum values of peak ground acceleration (A_max) also quantiles of the relevant probabilistic distributions are figured out in a desired future interval time in Iran. The main assumptions are Poissonian character of the seismic events flow and properties of the Gutenberg-Richter distribution law. The map of maximum possible values of A_max and also map of 90% quantile of distribution of maximum values of A_max on a future interval time 100 years is presented. According to the results, the maximum value of the A_max is estimated for Bandar Abbas as 0.3g and the minimum one is attributed to Esfahan as 0.03g. Finally, the estimated values in Bayesian approach are compared with what was presented applying probabilistic seismic hazard (PSH) methods based on the conventional Cornel (1968) method. The distribution function of A_max for future time intervals of 100 and 475 years are calculated for confidence limit of probability level of 90%.     

Geostatistical seismic inversion for non‐stationary patterns using direct sequential simulation and co‐simulation

Geostatistical seismic inversion methods are routinely used in reservoir characterisation studies because of their potential to infer the spatial distribution of the petro‐elastic properties of interest (e.g., density, elastic, and acoustic impedance) along with the associated spatial uncertainty. Within the geostatistical seismic inversion framework, the retrieved inverse elastic models are conditioned by a global probability distribution function and a global spatial continuity model as estimated from the available well‐log data for the entire inversion grid. However, the spatial distribution of the real subsurface elastic properties is complex, heterogeneous, and, in many cases, non‐stationary since they directly depend on the subsurface geology, i.e., the spatial distribution of the facies of interest. In these complex geological settings, the application of a single distribution function and a spatial continuity model is not enough to properly model the natural variability of the elastic properties of interest. In this study, we propose a three‐dimensional geostatistical inversion technique that is able to incorporate the reservoir's heterogeneities. This method uses a traditional geostatistical seismic inversion conditioned by local multi‐distribution functions and spatial continuity models under non‐stationary conditions. The procedure of the proposed methodology is based on a zonation criterion along the vertical direction of the reservoir grid. Each zone can be defined by conventional seismic interpretation, with the identification of the main seismic units and significant variations of seismic amplitudes. The proposed method was applied to a highly non‐stationary synthetic seismic dataset with different levels of noise. The results of this work clearly show the advantages of the proposed method against conventional geostatistical seismic inversion procedures. It is important to highlight the impact of this technique in terms of higher convergence between real and inverted reflection seismic data and the more realistic approximation towards the real subsurface geology comparing with traditional techniques.     

Seismic source zoning and maximum credible earthquake prognosis of the Greater Kashmir Territory,NW Himalaya

We present the seismic source zoning of the tectonically active Greater Kashmir territory of the Northwestern Himalaya and seismicity analysis (Gutenberg-Richter parameters) and maximum credible earthquake (m _max) estimation of each zone. The earthquake catalogue used in the analysis is an extensive one compiled from various sources which spans from 1907 to 2012. Five seismogenic zones were delineated, viz. Hazara-Kashmir Syntaxis, Karakorum Seismic Zone, Kohistan Seismic Zone, Nanga Parbat Syntaxis, and SE-Kashmir Seismic Zone. Then, the seismicity analysis and maximum credible earthquake estimation were carried out for each zone. The low b value (<1.0) indicates a higher stress regime in all the zones except Nanga Parbat Syntaxis Seismic Zone and SE-Kashmir Seismic Zone. The m _max was estimated following three different methodologies, the fault parameter approach, convergence rates using geodetic measurements, and the probabilistic approach using the earthquake catalogue and is estimated to be M _w 7.7, M _w 8.5, and M _w 8.1, respectively. The maximum credible earthquake (m _max) estimated for each zone shows that Hazara Kashmir Syntaxis Seismic Zone has the highest m _max of M _w 8.1 (±0.36), which is espoused by the historical 1555 Kashmir earthquake of M _w 7.6 as well as the recent 8 October 2005 Kashmir earthquake of M _w 7.6. The variation in the estimated m _max by the above discussed methodologies is obvious, as the definition and interpretation of the m _max change with the method. Interestingly, historical archives (～900 years) do not speak of a great earthquake in this region, which is attributed to the complex and unique tectonic and geologic setup of the Kashmir Himalaya. The convergence is this part of the Himalaya is distributed not only along the main boundary faults but also along the various active out-of-sequence faults as compared to the Central Himalaya, where it is mainly adjusted along the main boundary fault.     

Development of fragility curves by incorporating new spectral shape indicators and a weighted damage index: case study of steel braced frames in the city of Mashhad,Iran

In this study, strong ground motion record(SGMR) selection based on Eta(η) as a spectral shape indicator has been investigated as applied to steel braced frame structures. A probabilistic seismic hazard disaggregation analysis for the definition of the target Epsilon(ε) and the target Eta(η) values at different hazard levels is presented, taking into account appropriately selected SGMR's. Fragility curves are developed for different limit states corresponding to three representative models of typical steel braced frames having significant irregularities in plan, by means of a weighted damage index. The results show that spectral shape indicators have an important effect on the predicted median structural capacities, and also that the parameter η is a more robust predictor of damage than searching for records with appropriate ε values.     

Effect of seismic waves on the hydro-mechanical properties of fractured rock masses

The transmission of seismic waves in a particular region may influence the hydraulic properties of a rock mass, including permeability, which is one of the most important. To determine the effect of a seismic wave on the hydraulic behavior of a fractured rock mass, systematic numerical modeling was conducted. A number of discrete fracture network(DFN) models with a size of 20 m × 20 m were used as geometrical bases, and a discrete element method(DEM) was employed as a numerical simulation tool. Three different boundary conditions without(Type Ⅰ) and with static(Type Ⅱ) and dynamic(Type Ⅲ) loading were performed on the models, and then their permeability was calculated. The results showed that permeability in Type Ⅲ models was respectively 62.7% and 44.2% higher than in Type I and Type Ⅱ models. This study indicates that seismic waves can affect deep earth, and, according to the results, seismic waves increase the permeability and change the flow rate patterns in a fractured rock mass.