A hybrid discrete–finite element model for numerical simulation of geomaterials

A hybrid discrete–finite element model is introduced for simulation of mechanical behavior of geomaterials. The soil or rock is modeled as a system of discrete balls that interact through normal and shear springs. The balls can be bonded at the contact points to withstand the applied deviatoric stresses. The important feature of this model is that the confining walls that can be imagined for example as the surrounding membrane or the mold in a physical test are modeled by deformable finite elements. This allows simulation of laboratory test features more realistically compared to the situations where the surrounding walls are rigid. The relationships between micro- and macro-properties are investigated in this paper as well. These relationships and the corresponding curves are helpful tools in calibration of the numerical model for the macroscopic elastic properties.     

Rainfall over Oman and its teleconnection with El Niño Southern Oscillation

The Sultanate of Oman is located in the south-eastern part of the Arabian Peninsula and covers the larger part of the southern coasts of the Arabian Peninsula in both arid and semi-arid environments except for the southern part which is swept by the monsoon affecting the Arabian Sea during the period from June to September. The summer rainfall over Oman shows year-to-year variability, and this is caused by oceanic and atmospheric influences. In the present study, we tried to explore the influence of El Niño on the rainfall over Oman using different data sets. The empirical orthogonal function (EOF) technique employed to the zonal wind at 850 hPa for the 30-year period shows that the second and third modes of EOF are showing high variability over the Oman regions. The corresponding PCs were subjected to FFT analysis, and it showed a peak about 5–6 years. In addition to this, the zonal wind over the Oman regions is correlated with the global zonal wind and found a significant correlation (1 % significant level). It has already been proved that the wind and rainfall during summer monsoon is in phase. Moreover, the spectral analysis of rainfall at Masirah station and the Niño3.4 index show the similar mode of variability indicating a direct relationship. The correlation between rainfall and the Niño3.4 index is also showing a positive significant value, and therefore, it can be concluded that the El Niño in the Pacific favours rainfall over the Oman region.     

Investigation of moment-rotation relation in different joint types and evaluation of their effects on segmental tunnel lining

This paper presents the results of a numerical study that has been performed to investigate the different joint types that affect the moment-rotation relation and ultimate bending moment capacity of a joint. A 3D finite element method was adapted to establish elaborate numerical models of segments. To evaluate the possible differences in moment-rotation behavior between different joint types, 10 different joint types were simulated. Additionally, the effect of different joint types on behavior of a lining ring was investigated. The validity of the peridynamic simulation was tested by comparing results obtained in this paper against the results obtained in a study performed by Hordijk and Gijsbers. Observations of the results demonstrate that in a flat joint, the expanding of joint height increased the rotation stiffness of the joint in the linear branch, and accordingly, the ultimate bending moment of the joint increased. In the ring model, it was observed that the expanding of joint height led to the decreased of ring displacement and stress concentration in the joint. Whenever there were full surface contacts (contact two segment in total cross section) in the joints, the rotation stiffness of the joints in the linear branch became equal, and, as a result, the displacement in the ring model was the same in all joint types. In addition, it was observed that using a convex joint in the ring model increases the displacement of the ring. The ultimate bending moment of bolted joints was higher than that of joints without bolts, especially in the case of a negative moment.     

Hydrologic behavior and flood probability for selected arid basins in Makkah area, western Saudi Arabia

In arid regions, flash floods often occur as a consequence of excessive rainfall. Occasionally causing major loss of property and life, floods are large events of relatively short duration. Makkah area in western Saudi Arabia is characterized by high rainfall intensity that leads to flash floods. This study quantifies the hydrological characteristics and flood probability of some major wadis in western Saudi Arabia, including Na’man, Fatimah, and Usfan. Flood responses in these wadis vary due to the nature and rainfall distribution within these wadis. Rainfall frequency analysis was performed using selected annual maximums of 24-h rainfall from eight stations located in the area. Two of the most applied methods of statistical distribution, Gumbel’s extreme value distribution and log Pearson type III distribution, were applied to maximum daily rainfall data over 26 to 40 years. The Gumbel’s model was found to be the best fitting model for identifying and predicting future rainfall occurrence. Rainfall estimations from different return periods were identified. Probable maximum floods of the major wadis studied were also estimated for different return periods, which were extrapolated from the probable maximum precipitation.     

Upper mantle anisotropy of southeast Arabia passive margin [Gulf of Aden northern conjugate margin], Oman

In this study, we used data recorded by two consecutive passive broadband deployments on the Gulf of Aden northern margin, Dhofar region, Sultanate of Oman. The objective of these deployments is to map the young eastern Gulf of Aden passive continental margin crust and upper mantle structure and rheology. In this study, we use shear-wave splitting analysis to map lateral variations of upper mantle anisotropy beneath the study area. In this study, we found splitting magnitudes to vary between 0.33 and 1.0 s delay times, averaging about 0.6 s for a total of 17 stations from both deployment periods. Results show distinct abrupt lateral anisotropy variation along the study area. Three anisotropy zones are identified: a western zone dominated by NW–SE anisotropy orientations, an eastern zone dominated with NE–SW anisotropy orientations, and central zone with mixed anisotropy orientations similar to the east and west zones. We interpret these shorter wavelength anisotropy zones to possibly represent fossil lithospheric mantle anisotropy. We postulate that the central anisotropy zone may be representing a Proterozoic suture zone that separates two terranes to the east and west of it. The anisotropy zones west and east were being used indicative of different terranes with different upper mantle anisotropy signatures.     

Extension at the coast of the Makran subduction zone (Iran)

In the Makran subduction zone, earthquake focal mechanisms and geodetic data indicate that the deforming prism currently experiences N–S compression. However, palaeostress inversions performed on normal faults observed along the coast reveal local stress components consistent with N‐S extension. Previously proposed mechanisms such as gravitational collapse are not favoured by N–S compression and surface uplift. We propose that the observed kinematics result from transient stress reversals following large earthquakes. During the interseismic period (now), the region experiences N–S compression. However, following a large reverse rupture on the subduction interface, stresses in the inner wedge relax, enabling a brief period of extensional faulting before a compressive stress state is re‐established. This mechanism, also observed in other subduction zones, requires low overall stresses in the upper plate and that the margin ruptures in large megathrust earthquakes that result in nearly complete stress drops.     

Bed sediment entrainment by streamflow: State of the science

In fluvial sedimentology, bed sediment entrainment by streamflow has a decisive role in controlling several fluvial processes. Owing to its huge practical importance, the subject has been painstakingly explored for over a century. However, a detailed understanding of the mechanism of the bed sediment entrainment phenomenon achieved hitherto is far from complete. The central theme of bed sediment entrainment is occupied by the sediment entrainment threshold, which varies enormously in its qualitative definition, identification and quantification encompassing a broad range of spatiotemporal scales. This article presents the state of the science of the entrainment of non‐cohesive bed sediments under a steady‐unidirectional streamflow. It begins with the diverse definitions and representations of the entrainment threshold criterion from both qualitative and quantitative perspectives, scrutinising its suitability and ambit of applicability. Then, the effects of energetic factors that drive the entrainment threshold criterion are critically appraised. The indispensable mechanisms of bed sediment entrainment, including the theoretical background and modelling strategies, the role of turbulent bursting phenomenon and the phenomenological perspective into the origin of the scaling laws of sediment entrainment, are explained. Throughout the article, special emphasis is given to the strengths and weaknesses of the current state of the science. In addition, a deliberate attempt is made to invoke the thought‐provoking ideas on the multifarious features of bed sediment entrainment. Finally, the innovative perspectives on the bed sediment entrainment are provided and the concluding remarks are made, elucidating the major challenges and suggesting the prospective ways to resolve them as a future scope of research.     

Spectral wheat yield prediction modeling using SPOT satellite imagery and leaf area index

As wheat represents the main staple food and strategic crop in Egypt and worldwide and since remote sensing satellite imagery is the tool to obtain synoptic, multi-temporal, dynamic, and time-efficient information about any target on the Earth, the main objective of the current study is to use remote sensing satellite imagery to generate remotely sensed empirical preharvest wheat yield prediction models. The main input parameters of these models are spectral data either in the form of spectral reflectance data released from Satellite Pour l’Observation de la Terre (SPOT) 4 satellite imagery or in the form of spectral vegetation indices. The other input factor is leaf area index (LAI) that was measured by LAI Plant Canopy Analyzer. The four spectral bands of SPOT4 imagery are green, red, near-infrared, and middle infrared; the five vegetation indices that are forms of ratios between red and near-infrared bands are normalized difference vegetation index, ratio vegetation index, soil-adjusted vegetation index, difference vegetation index, and infrared percentage vegetation index. Another vegetation index is green vegetation index that is calculated through a ratio between green band and near-infrared band. Each of the above-mentioned factors was used as an input factor against wheat yield to generate wheat yield prediction models. All generated models are site-specific limited to the area and the environment and could be applicable under similar conditions in Egypt. The study was carried out in Sakha experimental station by using the dataset from two wheat season 2007/2008 and 2009/2010. The total wheat area was 1.3 ha cultivated by Sakha 93 cultivar. Modeling and validation process were carried out for each season independently. Modeled yield was tested against reported yield through two common statistical tests; the standard error of estimate between modeled yield and reported yield, and the correlation coefficient for a direct regression analysis between modeled and reported yield with each generated model. Generally, as shown from the correlation coefficient of the generated models, green and middle infrared bands did not show good accuracy to predict wheat yield, while the other spectral bands (red and near-infrared) bands showed high accuracy and sufficiency to predict yield. This was proven through the correlation coefficient of the generated models and through the generated models with the wheat crops for the two seasons. Accordingly, the green vegetation index that is generally calculated from green and near-infrared bands showed relatively lower accuracy than the rest of the vegetation index models that are calculated from red and near-infrared bands. LAI showed high accuracy to predict yield as shown from the statistical analysis. The models are applicable after 90 days from sowing stage and applicable in similar regions with the same conditions.     

A new look on the Jurassic formations of the western part of Iraq

The Jurassic succession of western Iraq includes the Ubaid, Hussainiyat, Amij, Muhaiwir, and Najmah formations. Each formation is composed of basal clastic unit overlain by upper carbonate unit. Extensive and huge erosional unconformity occurred at the Triassic–Jurassic boundary due to marked shifting of structural (E–W) strike of the Triassic (Rhaetic) Zor Hauran Formation to (NE–SW) Jurassic formations. Sea level falling (lowstand system tracts) would result in the progradation of the land on the expense of the sea forming the recognizable progradation of the fluvial and deltaic deposits of the lower clastic units of the Jurassic formations, whereas sea level rising (highstand system tract), i.e., sea prograding, causes deposition of the carbonate units of the Jurassic formations system. This progradation resulted to various carbonate environments of deposition ranging from subtidal, intertidal, to supratidal. The main target on most of the exploration blocks in the western part of Iraq focused on the lower Paleozoic successions, whereas prospects in Triassic, Jurassic, and lower Cretaceous targets are less extensive but may have significant potential on certain blocks in both stratigraphic and structural traps. The western part of Iraq was subjected to intermittent pulses of uplifting (sea regression) and subsidence (sea transgression) to form the Jurassic basin system in the area. The Jurassic formations lack the characteristics of petroleum systems. In contrast, in the central and northern parts of Iraq, the Jurassic formations (Najmah and Gotnia formations) were deposited in subsiding basins in which the reservoir and sealed evaporitic rocks existed. In turn, in the western desert of Iraq, the Jurassic formations lack these petroleum system characteristics. Hence, it can be proposed that the petroleum–nonpetroleum inflection could be proposed in the east of area km 160.     

Risk assessment of polycyclic aromatic hydrocarbons in the West Port semi-enclosed basin (Malaysia)

Sediment samples collected from the West Port, the west coastal waters of Malaysia, were analyzed by standard methods to determine the degree of hydrocarbon contamination and identify the sources of polyaromatic hydrocarbons (PAHs). Concentrations of PAHs in the port sediments ranged from 100.3 to 3,446.9 μg/kg dw. The highest concentrations were observed in stations close to the coastline, locations affected by intensive shipping activities and industrial input. These were dominated by high-molecular-weight PAHs (4–6 rings). Source identification showed that PAHs originated mostly pyrogenically, from the combustion of fossil fuels, grass, wood, and coal or from petroleum combustion. Regarding ecological risk estimation, only station 7 was moderately polluted, the rest of the stations suffered rare or slight adverse biological effects with PAH exposure in surface sediment, suggesting that PAHs are not considered as contaminants of concern in the West Port.