A new empirical estimator of coseismic landslide displacement for Zagros Mountain region (Iran)

Earthquake-induced landslides are responsible worldwide for significant socioeconomic losses and historically have a prominent position in the list of natural hazards affecting the Iran plateau. As a step toward the development of tools for the assessment and the management of this kind of hazard at regional scale, an empirical estimator of coseismic displacements along potential sliding surfaces was obtained through a regression analysis for the Zagros region, a mountainous Iranian region subjected to earthquake-induced landslides. This estimator, based on the Newmark’s model, allows to evaluate the expected permanent displacement (named “Newmark displacement”) induced by seismic shaking of defined energy on potential sliding surface characterized by a given critical acceleration. To produce regression models for Newmark displacement estimators, a data set was constructed for different critical acceleration values on the basis of 108 accelerometric recordings from 80 Iranian earthquakes with moment magnitudes between 3.6 and 7. The empirical estimator has a general form, proposed by Jibson (Eng Geol 91:209–218, 2007), relating Newmark displacement to Arias intensity (as parameter representing the energy of the seismic forces) and to critical acceleration (as parameter representing the dynamic shear resistance of the sliding mass). As an example of application, this relation was employed to provide a basic document for earthquake-induced landslide hazard assessment at regional scale, according to a method proposed by Del Gaudio et al. (Bull Seismol Soc Am 93:557–569, 2003), applied to the whole Iranian territory, including Zagros region. This method consists in evaluating the shear resistance required to slopes to limit the occurrence of seismically induced failures, on the basis of the Newmark’s model. The obtained results show that the exposure to landslide seismic induction is maximum in the Alborz Mountains region, where critical accelerations up to ~0.1 g are required to limit the probability of seismic triggering of coherent type landslides within 10% in 50 years.     

A Study on the Surface Wettability of Clastic Rocks with Potential Application for CO2 Storage Sites

Natural Resources Research - There have been many studies carried out in the past decades attempting to develop strategies for a safe injection of CO2 into storage sites without leakage and...     

Three-parameter Radon transform in layered transversely isotropic media

Radon transform is a powerful tool with many applications in different stages of seismic data processing, because of its capability to focus seismic events in the transform domain. Three-parameter Radon transform can optimally focus and separate different seismic events, if its basis functions accurately match the events. In anisotropic media, the conventional hyperbolic or shifted hyperbolic basis functions lose their accuracy and cannot preserve data fidelity, especially at large offsets. To address this issue, we propose an accurate traveltime approximation for transversely isotropic media with vertical symmetry axis, and derive two versions of Radon basis functions, time-variant and time-invariant. A time-variant basis function can be used in time domain Radon transform algorithms while a time-invariant version can be used in, generally more efficient, frequency domain algorithms. Comparing the time-variant and time-invariant Radon transform by the proposed basis functions, the time-invariant version can better focus different seismic events; it is also more accurate, especially in presence of vertical heterogeneity. However, the proposed time-invariant basis functions are suitable for a specific type of layered anisotropic media, known as factorized media. We test the proposed methods and illustrate successful applications of them for trace interpolation and coherent noise attenuation.     

Steady critical discharge rates from vertical and horizontal wells in fresh–saline aquifers with sharp interfaces

In this study, the critical (or maximum) discharge rates before saline water enters a well were determined for vertical and horizontal wells in a freshwater aquifer which is separated from a static saline aquifer by a sharp interface. Flow around the well was solved by integration of a point sink solution along the well axis, and both the critical discharge rate and critical interface rise were determined through a comparison of the heads and vertical gradients at the saline–fresh water interface. The rates were determined for vertical and horizontal wells with various lengths and depths for different aquifer salinities. Results were generalized by drawing dimensionless type curves. The results showed that the dimensionless total critical discharge rates are higher for the longer horizontal wells and longer vertical wells with a certain bottom depth, and they almost linearly decrease with well depth at rates of 0.7–0.9. For the dimensionless well length of 0.2, the dimensionless total discharge rate of a horizontal well is about 0.1 more than that of a vertical well with the same length and well-top depth. Also, the critical discharge rates per unit length of well are inversely proportional to well length and remarkably higher for shallower wells. Additionally, the critical pumping rate is proportional to the salinity difference of the aquifers. These results were confirmed by comparison to existing solutions for vertical wells with dimensionless lengths of 0.2, 0.5 and 0.6, and for critical interface rises in the range of 0.75–0.9.     

Benthic nutrient fluxes along the Laurentian Channel: Impacts on the N budget of the St. Lawrence marine system

Water column concentrations and benthic fluxes of dissolved inorganic nitrogen (DIN) and oxygen (DO) were measured in the Gulf of St. Lawrence and the Upper and Lower St. Lawrence Estuary (USLE and LSLE, respectively) to assess the nitrogen (N) budget in the St. Lawrence (SL) system, as well as to elucidate the impact of bottom water hypoxia on fixed-N removal in the LSLE. A severe nitrate deficit, with respect to ambient phosphate concentrations (N*∼−10 μmol L⁻¹), was observed within and in the vicinity of the hypoxic bottom water of the LSLE. Given that DO concentrations in the water column have remained above 50 μmol L⁻¹, nitrate reduction in suboxic sediments, rather than in the water column, is most likely responsible for the removal of fixed N from the SL system. Net nitrate fluxes into the sediments, derived from pore water nitrate concentration gradients, ranged from 190 μmol m⁻² d⁻¹ in the hypoxic western LSLE to 100 μmol m⁻² d⁻¹ in the Gulf. The average total benthic nitrate reduction rate for the Laurentian Channel (LC) is on the order of 690 μmol m⁻² d⁻¹, with coupled nitrification-nitrate reduction accounting for more than 70%. Using average nitrate reduction rates derived from the observed water column nitrate deficit, the annual fixed-N elimination within the three main channels of the Gulf of St. Lawrence and LSLE was estimated at 411 × 10⁶ t N, yielding an almost balanced N budget for the SL marine system.     

Investigation on the Effect of Geometrical and Geotechnical Parameters on Elongated Offshore Piles Using Fuzzy Inference Systems

Among numerous offshore structures used in oil extraction, jacket platforms are still the most favorable ones in shallow waters. In such structures, log piles are used to pin the substructure of the platform to the seabed. The pile’s geometrical and geotechnical properties are considered as the main parameters in designing these structures. In this study, ANSYS was used as the FE modeling software to study the geometrical and geotechnical properties of the offshore piles and their effects on supporting jacket platforms. For this purpose, the FE analysis has been done to provide the preliminary data for the fuzzy-logic post-process. The resulting data were implemented to create Fuzzy Inference System (FIS) classifications. The resultant data of the sensitivity analysis suggested that the orientation degree is the main factor in the pile’s geometrical behavior because piles which had the optimal operational degree of about 5° are more sustained. Finally, the results showed that the related fuzzified data supported the FE model and provided an insight for extended offshore pile designs.     

Thermal stability of a layer with permeable boundaries and variable gravitational field

The aim of this paper is to investigate the thermal stability of a fluid layer with permeable boundaries and a variable gravitational field. It is observed that the principle of exchange of stabilities is valid when the layer is heated form below and the complex growth rate of an arbitrary oscillatory mode exists outside of a circle whose radius depends upon the permeability parameter, Prandtl number and wavelength of the mode. In the case of a layer heated from below, gravity increasing upward has a destabilizing effect whereas the permeability parameter has a stabilizing effect.     

Weakly compressible SPH simulation of cnoidal waves with strong plunging breakers

Ocean Dynamics - Hydrodynamics of highly nonlinear cnoidal waves and their subsequent strong plunging breakers are among the least understood and most significant issues in coastal engineering. In...     

Makran ophiolitic basalts (SE Iran) record Late Cretaceous Neotethys plume-ridge interaction

ABSTRACT

The Makran complex in southeast Iran provides a spectacular subduction-related accretionary complex to understand the mechanism of oceanic accretion and the evolution of subduction zones. In this paper, we present new major and trace element data as well as isotopic compositions of mafic volcanic blocks from the Makran ophiolitic mélange complex (OMC). Our aim is to assess the genesis of these rocks and discuss their implications on the evolution of Neotethys Ocean. These volcanic blocks are composed mainly of basalts with minor trachytes. The Makran lavas are occasionally interlayered with tuff layers. Zircons from these tuffs give U-Pb ages of 95 Ma, which is well in accordance with the reposted microfossil data for the interlayered pelagic limestones with pillow lavas. Makran basalts can be geochemically subdivided into four groups; normal to transitional MORB, enriched-MORB, Plume-type MORB and alkaline (-OIB-like) basalts. The OIB-like pillow lavas are represented by high values of Th/Tb (6.3–7.4) which are higher than other basalts (group 1 = 0.3–0.8; groups 2 = 0.7–1.6; group 3 = 1.58–1.36).¹⁴³Nd/¹⁴⁴Nd_(t) ratios for basalts ranges from 0.51247 to 0.51292, whereas ⁸⁷Sr/⁸⁶Sr_(t) isotopic composition of the OMC lavas varies from 0.704433 to 0.709466. The Pb isotopic composition of the lavas are quite high, ranging from 15.49–15.66 for ²⁰⁷Pb/²⁰⁴Pb_(t), 18.09–19.12 for ²⁰⁶Pb/²⁰⁴Pb_(t) and 37.80–39.23 for ²⁰⁸Pb/²⁰⁴Pb_(t). The chemistry of these rocks suggests that they were formed most likely in an oceanic setting with clear plume-ridge interaction. These rocks can form from partial melting of a highly heterogeneous mantle source, which is extensively metasomatized with deep mantle OIB-type components. We suggest these rocks have been generated in an oceanic ridge with plume-ridge interaction, similar to the Iceland-Reykjanes Ridge, before being fragmented and accreted into the Makran accretionary complex.     

Woody species diversity and forest structure from lowland to montane forest in Hyrcanian forest ecoregion

Alborz Mountains host Caspian Hyrcanian forest ecoregion along the northern slopes and forest steppe ecoregion in highlands. Hyrcanian forest covers the southeastern part of Caucasus biodiversity hotspot and is of great biogeographic importance. Altitudinal pattern and correlation between woody species biodiversity (DIV), forest structure ((stem density (DEN), mean basal area (MBA) and mean height class (MHC)) and disturbance (DIS) were explored along 2,400 m altitudinal gradient in Hyrcanian relict forest, Central Alborz Mountains. Vegetation changes from lowland forest (LoF) to mid- altitude forest (MiF) and montane forest (MoF) in this area. The altitudinal gradient was divided into twelve 200 m elevational belts. Point centered quarter method (PCQM) with 96 sampling points and 83 vegetation samples by plot method (PM) were used to record field data. Shannon-Wiener index and Pearson coefficient were used for diversity and correlation analysis. The results showed that DEN decreased linearly, MBA and MHC showed relatively hump shaped and DIS showed a reverse hump shaped pattern of change along altitudinal gradient. Woody species diversity decreased non-steadily from LoF to MoF. Transitional vegetations of Carpinus-Fagus and Fagus-Quercus represented higher diversity of woody taxa compared to adjacent homogenous communities. Significant correlation was observed between altitude and all parameters: DEN with MBA, DIS and DIV; MBA with DIS; MHC with DIS along with DIV; and DIS with DIV at the study area scale. Surprisingly,correlation between studied parameters differed within each vegetation type. Altitude probably acts as a proxy for human and environmental driving forces in this area. Stability of warm and wet condition, season length, soil depth along with forest accessibility probably influences the altitudinal pattern of the studied parameters. Disturbance affects forest structure and consequently diversity; especially in lowlands. The obtained results recommend using both forest biodiversity and mensuration data in management process of forest ecosystems.