Simultaneous estimation of shape factor and depth of subsurface cavities from residual gravity anomalies using feed-forward back-propagation neural networks

We develop a new method of using feed-forward back-propagation (FFBP) neural networks to simultaneously estimate shape factor and depth of gravity anomalies. The advantages compared to neural network methods are the following: no pre-assumptions are made on source shape, the FFBP neural network estimates both depth and shape factor of source bodies and, once trained, works well for any new data in the training space, without repeating the initial calculations.     

Local site effect of a clay site in Shiraz based on seismic hazard of Shiraz Plain

Propagation of seismic waves through soil layers would drastically change the frequency content and amplitude-based features of ground motions at the surface. These alterations are known as seismic site effects. Computation of site effects of high-populated areas such as large cities is of great importance (e.g., it is used in development of seismic microzonation of a region). Shiraz is one of the most populous cities of Iran and is located in a high seismic hazardous region. A representative clay site in this city is selected to assess local site effects. The time series and random vibration theory procedure in the frequency domain are implemented to analyze the aforementioned site. Furthermore, the nonlinear dynamic soil behavior is simulated by the equivalent linear method and the nonlinear method via DEEPSOIL program. Three types of soil column uncertainties such as shear wave velocity, modulus reduction, and damping ratio of soil layers as well as depth of underlying rock half-space (D _bed) are considered herein. The mean amplification and standard deviation of natural logarithm of amplification factors are computed for a variety of analysis types. The results of the current study show that the computed mean and standard deviation of amplification factor in ln units by considering only V _S uncertainty are in good agreement with the corresponding ones by considering V _S and modulus reduction and damping ratio variabilities simultaneously for the studied site. Furthermore, it seems that the effect of bedrock depth in definition of spectral shapes of the Iranian seismic building code should be taken into account.     

Mantle flow patterns at the Neyriz Paleo-spreading center, Iran

The mantle peridotites of Neyriz record two successive episodes of plastic deformations; the first one related to the igneous accretion of the lithosphere and the second one developed during the first stage of the emplacement of the peridotites. These two events have been distinguished on the basis of microstructural criteria. The diapiric pattern, particularly relevant to the mantle process beneath spreading ridges, features vertical flow lines and elliptic flow plane trajectories in a pipe and extends along the ridge axis about 5 km. These structures rotate to horizontal and diverge in every direction in a narrow transition zone, a few hundred meters thick, below the Moho discontinuity. Such a diapiric pattern has been recognized in a few places along the Neyriz paleo-ridge. A large amount of magma passed through these mantle diapirs that were probably the main zones feeding the overlying magma chamber. The most common pattern features very regular structures over several kilometers along the strike of the paleo-ridge: the flow plane dips away from the ridge axis, and the flow line is parallel to the spreading direction. This flow pattern is frozen during the gradual accretion of the lithospheric mantle away from the ridge in a steady-state spreading regime. A shear-sense inversion at just below the Moho is commonly observed, pointing to forced asthenospheric flow. The reconstructed orientation of the Neyriz paleo-spreading center is 105°, compatible with the geometry and orientation of harzburgite foliations and lineations and sheeted dikes.     

Flood frequency analysis based on simulated peak discharges

Flood frequency approaches vary from statistical methods, directly applied on the observed annual maximum flood series, to adopting rainfall–runoff simulation models that transform design rainfalls to flood discharges. Reliance on statistical flood frequency analysis depends on several factors such as the selected probability distribution function, estimation of the function parameters, possible outliers, and length of the observed flood series. Through adopting the simulation approach in this paper, watershed-average rainfalls of various occurrence probabilities were transformed into the corresponding peak discharges using a calibrated hydrological model. A Monte Carlo scheme was employed to consider the uncertainties involved in rainfall spatial patterns and antecedent soil moisture condition (AMC). For any given rainfall depth, realizations of rainfall spatial distribution and AMC conditions were entered as inputs to the model. Then, floods of different return periods were simulated by transforming rainfall to runoff. The approach was applied to Tangrah watershed in northeastern Iran. It was deduced that the spatial rainfall distribution and the AMCs exerted a varying influence on the peak discharge of different return periods. Comparing the results of the simulation approach with those of the statistical frequency analysis revealed that, for a given return period, flood quantiles based on the observed series were greater than the corresponding simulated discharges. It is also worthy to note that existence of outliers and the selection of the statistical distribution function has a major effect in increasing the differences between the results of the two approaches.     

Static and dynamic analysis on slope stability using a DFN-DEM approach on the right abutment of the Karun 4 dam

Uncertainty in input fracture geometric parameters during analysis of the stability of jointed rock slopes is inevitable and therefore the stochastic discrete fracture network (DFN) — distinct element method (DEM) is an efficient modeling tool. In this research, potentially unstable conditions are detected in the right abutment of the Karun 4 dam and downstream of the dam body as a case study. Two extreme states with small and relatively large block sizes are selected and a series of numerical DEM models are generated using a number of validated DFN models. Stability of the rock slope is assessed in both static and dynamic loading states. Based on the design basis earthquake (DBE) and maximum credible earthquake (MCE) expected in the dam site, histories of seismic waves are applied to analyze the stability of the slope in dynamic earthquake conditions. The results indicate that a MCE is likely to trigger sliding of rock blocks on the rock slope major joint. Furthermore, the dynamic analysis also shows a local block failure by the DBE, which can consequently lead to slope instability over the long term. According to the seismic behavior of the two models, larger blocks are prone to greater instability and are less safe against earthquakes.

     

Paleoenvironmental reconstruction of the upper Cretaceous succession (Abtalkh Formation) of the Kopeh-Dagh Basin,northeastern Iran based on foraminiferal and palynological analyses

Arabian Journal of Geosciences - In order to reconstruct sea surface water productivity and sea floor oxygenation during late Campanian-Maastrichtian, planktonic and benthic foraminiferal...     

Mathematical modeling of competition for ammonium among Bacteria, Archaea and cyanobacteria within cyanobacterial mats: Can ammonia-oxidizers force nitrogen fixation?

Molecular analysis of cyanobacterial mat communities indicated that cyanobacteria, ammonia-oxidizing Archaea (AOA), and ammonia-oxidizing bacteria (AOB) coexist in those systems, competing for ammonium; this situation would imply competitive exclusion. We attempted to model how ammonia utilization niche partitioning occurs, and how ammonium levels can influence the interaction between those groups in a one-dimensional diffusionlimited system using Michaelis-Menten kinetics to describe ammonium consumption by each of those three groups. In our model, AOAs were able to dominate ammonium uptake by the community under most circumstances, except for unrealistically high (millimolar) levels of ammonium, where AOBs gained advantage. Cyanobacteria were unable to effectively compete for ammonium with either AOBs or AOAs throughout the mat at all ammonium concentrations and cell counts, suggesting that the presence of AOAs or AOBs forces cyanobacteria into nitrogen fixation mode. Such interaction can make cyanobacterial mats a net nitrogen source, as well as provide a carbon-independent energy transfer pathway from primary producers to the rest of the ecosystem.     

Laboratory Study of Topographic Effects on the Near-surface Tornado Flow Field

To study topographic effects on the near-surface tornado flow field, the Iowa State University tornado simulator was used to simulate a translating tornado passing over three different two-dimensional topographies: a ridge, an escarpment and a valley. The effect of the translation speed on maximum horizontal wind speeds is observed for translation speeds of 0.15 and 0.50 \(\hbox {m}\,\hbox {s}^{-1}\), with the lower value resulting in a larger maximum horizontal wind speed. The tornado translation over the three topographies with respect to flat terrain is assessed for changes in: (a) the maximum horizontal wind speeds in terms of the flow-amplification factor; (b) the maximum aerodynamic drag in terms of the tornado speed-up ratio; (c) the maximum duration of exposure at any location to high wind speeds of a specific range in terms of the exposure amplification factor. Results show that both the maximum wind amplification factor of 14%, as well as the maximum speed-up ratio of 14%, occur on the ridge. For all topographies, the increase in aerodynamic drag is observed to be maximized for low-rise buildings, which illustrates the importance of the vertical profiles of the horizontal wind speed near the ground. The maximum exposure amplification factors, estimated for the range of wind speeds corresponding to the EF2 (50–60 \(\hbox {m}\,\hbox {s}^{-1}\)) and EF3 (61–75 \(\hbox {m}\,\hbox {s}^{-1})\) scales, are 86 and 110% for the ridge, 4 and 60% for the escarpment and ? 6 and 47% for the valley, respectively.     

Experimental studies to constrain parental magma of Malangtoli volcanics from Singhbhum craton of the eastern Indian shield

Malangtoli volcanics of the Singhbhum craton of the eastern Indian shield is one of the important Proterozoic lava suites. Experimental studies on 1 atmosphere pressure constrain the parental magma type and temperature range of crystallization of the parent magma (deduced to be in the range of 1500°C to 1200°C). The experimental studies show that at 1500°C, plagioclase is the first phase to crystallize, followed by few opaques which join along with plagioclase at 1450°C. At subsequent lower temperature (1400°C-1300°C), plagioclase and opaque continue to crystallize. At 1250°C plagioclase and opaque still persist while pyroxene appears first and liquid (glass) still remains. Appearance of opaque minerals (magnetite and illmenite) at both ~1400°C and ~1300°C indicate oscillation of oxygen fugacity in the parent magma, petrographically documented by coarser phenocrysts as well as finer or peripheral tiny grains. Use of tectonic discrimination diagrams (based on discrimination factors F₁-F₂ and FeO^t/MgO vs. TiO₂) shows an island arc tholeiitic affinity for Malangtoli volcanic, suggests that the role of proto-plate convergence in Singhbhum architecture played an important role to build up Malangtoli volcanics during Proterozoic.     

Experimental estimation of velocities and anisotropy of a series of Swedish crystalline rocks and ores

To provide a guide for future deep (<1.5 km) seismic mineral exploration and to better understand the nature of reflections imaged by surface reflection seismic data in two mining camps and a carbonatite complex of Sweden, more than 50 rock and ore samples were collected and measured for their seismic velocities. The samples are geographically from the northern and central parts of Sweden, ranging from metallic ore deposits, meta‐volcanic and meta‐intrusive rocks to deformed and metamorphosed rocks. First, ultrasonic measurements of P‐ and S‐wave velocities at both atmospheric and elevated pressures, using 0.5 MHz P‐ and S‐wave transducers were conducted. The ultrasonic measurements suggest that most of the measured velocities show positive correlation with the density of the samples with an exception of a massive sulphide ore sample that shows significant low P‐ and S‐wave velocities. The low P‐ and S‐wave velocities are attributed to the mineral texture of the sample and partly lower pyrite content in comparison with a similar type sample obtained from Norway, which shows significantly higher P‐ and S‐wave velocities. Later, an iron ore sample from the central part of Sweden was measured using a low‐frequency (0.1–50 Hz) apparatus to provide comparison with the ultrasonic velocity measurements. The low‐frequency measurements indicate that the iron ore sample has minimal dispersion and attenuation. The iron ore sample shows the highest acoustic impedance among our samples suggesting that these deposits are favourable targets for seismic methods. This is further demonstrated by a real seismic section acquired over an iron ore mine in the central part of Sweden. Finally, a laser‐interferometer device was used to analyse elastic anisotropy of five rock samples taken from a major deformation zone in order to provide insights into the nature of reflections observed from the deformation zone. Up to 10% velocity‐anisotropy is estimated and demonstrated to be present for the samples taken from the deformation zone using the laser‐interferometery measurements. However, the origin of the reflections from the major deformation zone is attributed to a combination of anisotropy and amphibolite lenses within the deformation zone.