Sedimentology,mineralogy, and geochemistry of the Late Quaternary Meyghan Playa sediments,NE Arak,Iran: palaeoclimate implications

Playas are shallow ephemeral lakes that form in arid and semi-arid regions. Iran has a large number of playas such as Meyghan Playa, which is located in the northeast of Arak city that borders the central Iran and Sanandaj-Sirjan zones. This study aims to investigate the mineralogical, sedimentological, and geochemical characteristics of the playa sediments. In order to determine the palaeoenvironment, we carried out X-ray diffraction (XRD), X-ray fluorescence (XRF), and scanning electron microscopy (SEM) studies. Meyghan Playa sediments consist of very fine-grained sediments and contain both evaporite and clastic minerals. The evaporite minerals include calcite, gypsum, halite, glauberite, and thenardite, whereas clastic minerals are quartz and clay. The calcite abundance decreases from the margin to the central portion of the playa but gypsum and halite abundances show an increasing trend from the margin to the center. This observation is consistent with the general zonation of other playas. Variations of calcite and gypsum concentration profiles present increasing and decreasing trends with depth, which could be ascribed to the changes in climatic factors. These factors include brine chemical modifications owing to changes in evaporation and precipitation rates and variations in relative abundance of anions-cations or in the rate of clastic and evaporite minerals due to variations in the freshwater influx (climatic changes) with time. A decrease in calcite and increase in sulfate minerals (especially gypsum) with depth is probably due to the higher water level and rainfall, a more humid climate, and salinity variations.     

A new insight into onset of inertial flow in porous media using network modeling with converging/diverging pores

The network modeling approach is applied to provide a new insight into the onset of non-Darcy flow through porous media. The analytical solutions of one-dimensional Navier-Stokes equation in sinusoidal and conical converging/diverging throats are used to calculate the pressure drop/flow rate responses in the capillaries of the network. The analysis of flow in a single pore revealed that there are two different regions for the flow coefficient ratio as a function of the aspect ratio. It is found that the critical Reynolds number strongly depends on the pore geometrical properties including throat length, average aspect ratio, and average coordination number of the porous media, and an estimation of such properties is required to achieve more reliable predictions. New criteria for the onset of non-Darcy flow are also proposed to overcome the lack of geometrical data. Although the average aspect ratio is the main parameter which controls the inertia effects, the effect of tortuosity on the onset of non-Darcy flow increases when the coordination number of media decreases. In addition, the higher non-Darcy coefficient does not essentially accelerate the onset of inertial flow. The results of this work can help to better understand how the onset of inertial flow may be controlled/changed by the pore architecture of porous media.     

A numerical study on the impact of thermal alterations in porous media during hot fluid injection process employing a modified Boussinesq model

The Oberbeck-Boussinesq (OB) approximation is widely employed as a simplifying assumption for density-dependent flow problems. It reduces the governing differential equations to simpler forms, which can be handled analytically or numerically. In this study, a modified OB model is formulated to account for the variation of rock permeability and porosity with temperature during the hot fluid injection process in an oil-saturated porous medium under the assumption of local thermal equilibrium (LTE). The mathematical model is solved numerically using a fully implicit control volume finite difference discretization with the successive over relaxation (SOR) method to handle the non-linearity. Subsequently, the numerical model is validated with the analytical solution of the simplified problem successfully. Through detailed sensitivity analyses, the simulation results reveal the hot fluid injection rate as the most important operational parameter to be optimized for a successful thermal flood. The numerical runs show that that for single-phase core-flood simulation, the effect of temperature on the rock absolute permeability and porosity can be neglected without introducing any significant errors in the estimated recovery and temperature profile.     

Analysis of subdiffusion in disordered and fractured media using a Grünwald-Letnikov fractional calculus model

The increasing applications of fractional calculus in simulating the anomalous transport behavior in disordered and fractured heterogeneous porous media has grown rapidly over the past decade. In the present study, a temporal fractional flux relationship is employed as a constitutive equation to relate the volumetric flow rate to the gradient of the pore pressure. The novelty of this paper entails interpreting the time fractional derivative operator in the flux relationship by the Grünwald-Letnikov (G-L) definition as opposed to the Caputo interpretation which has been widely considered. Subsequently, a numerical scheme based on the block-centered finite-difference discretization is formulated to handle the resulting non-linear fractional diffusion model. In addition, a linear stability analysis is successfully performed to establish the stability criterion of the developed numerical scheme. An expression for the modified incremental material balance index was derived to assess the effectiveness of the numerical discretization process. Finally, numerical experiments were performed to provide qualitative insights into the nature of pressure evolution in a hydrocarbon reservoir under the influence subdiffusion. In summary, the results establish that subdiffusion regime results in the development of higher pressure drop in the reservoir. This paper will provide a strong foundation for researchers interested in investigating anomalous diffusion phenomena in porous media.     

Spatial variability of climatic hazards in Bangladesh

Natural Hazards - Climate change is evident with the extreme climatic indices changing all over the world. Bangladesh is one of the most vulnerable countries to climate change. The patterns of...     

Intra‐ and interannual export of nitrogen and phosphorus in the subtropical Richmond River catchment,Australia

Nitrogen and phosphorus concentrations were measured and exports were calculated in the subtropical Richmond River catchment between July 1994 and June 1996. A stratified sampling approach was adopted owing to the extreme and rapid changes in riverine discharge, which varied by up to 10 000 times over the study period. Nutrient concentrations were lowest during baseflow. During storm discharge, dissolve inorganic and organic and particulate nitrogen and phosphorus concentrations increased two‐ to five fold, and followed hysteresis patterns that were attributed to the integration and/or depletion of catchment nutrient sources during an event. Dissolved organic nitrogen and particulate phosphorus were the dominant nutrient forms. Land use and antecedent conditions had a large influence on nutrient concentrations and exports. During the 1995–96 year (slightly above the mean annual discharge) 96% of nitrogen and 98% of phosphorus loads were transported in less than 6% of the time. When averaged across the catchment, monthly riverine nutrient loads varied by up to 1061‐fold during the study and exports were approximately four‐fold greater during the second year relative to the first. The subtropical Richmond River catchment has greater intra‐ and potential interannual variability in nutrient loads and exports when compared with temperate catchments from other parts of the world. It is suggested that in tropical and subtropical Australian catchments with large intra‐ and interannual variation in discharge, the need for parameterizing the antecedent conditions, such as the degree of nutrient storage, may make it difficult to model spatial and temporal (short time‐scale) nutrient exports. Copyright © 2000 John Wiley & Sons, Ltd.     

Planar and nonplanar quantum dust ion-acoustic Gardner double layers in multi-ion dusty plasma

The longitudinal fast solitary waves induced by weakly relativistic positron showers of astrophysical origin are studied in a plasma system contaminated with some massive impurities in presence of superthermal effects. The superthermal effects are due to the high energy electrons. The impurities are dust corpuscles with positive and negative charges. It is noticed that increase in the kappa parameter of electrons and relativistic streaming factor of weakly relativistic positron shower, negative dust concentration invoke an enhancement in the strength of solitary wave. On the other hand increase in the shower’s temperature as well as positive dust concentration diminish the solitary hump strength. It is worth to mention that only hump type compressive fast solitary waves are predicted by our model, for the given set of plasma parameters, because the convective coefficient of the nonlinear governing equation for solitary wave remains positive in considered regime of interaction for plasma and positron shower. Our calculations in linear regime predict both the fast and slow positron shower induced longitudinal, electrostatic perturbations. Our results may be of importance in understanding the nonlinear propagation of waves in doped astrophysical superthermal plasmas with relativistic positron showers.     

Site response and source spectra of S waves in the Zagros region, Iran

S wave amplitude spectra from shallow earthquakes with magnitudes ranging between 4.2 and 6.2 in the Zagros region of Iran that occurred between 1998 and 2008 are used to examine source parameters and site response of S waves. A generalized inversion scheme has been used to separate the source, propagation path, and local site effects from S wave spectra. For removing the trade-off between source and site terms and propagation effects (including geometric and anelastic attenuation), the spectral amplitudes of the records used were corrected for attenuation and geometrical spreading function using a path model proposed by Zafarani and Soghrat (Bull Seism Soc Am 102:2031–2045, 2012) for the region. We assume a Brune’s point source model to retrieve source parameters like corner frequency, moment magnitude, and high-frequency fall off coefficient, for each event. When the source spectra are interpreted in terms of Brune’s model, the average stress drops obtained are about 7.1 and 5.9 MPa (71 and 59 bars), respectively for the eastern and western Zagros regions. Stress drops range from 1.4 to 35.0 MPa (14 to 350 bars), with no clear dependence on magnitude. The results in terms of stress drop and S wave seismic energy indicate that the Zagros events are more similar to interplate earthquakes of western North America than to intraplate events of eastern North America. The method also provides us with site responses for all 40 stations individually and is an interesting alternative to other methods, such as the H/V method. A new empirical relationship between body-wave magnitudes and moment magnitude has been proposed for the Iranian plateau using derived seismic moment from the inversion.     

The cyclical reversible transitions of the universe’s evolution

In this work we propose cyclical reversible transitions as the scenario in which the universe evolves, through a series consisting of reversible expansion, temporary stability, and contraction. Our model is based on the comparison between local and global time-dependent densities {ρ ₀(τ ₀),ρ(τ)} instead of the critical density ρ _c, local and global time-dependent Hubble parameters {H ₀(τ ₀),H(τ)}, and the variations {Δρ(τ),ΔH(τ)} due to cosmological chaotic fluctuations, which are generally ignored in certain oscillating models. By taking into account all these factors, a rate equation in the form of (H ₀/H)² ∝(ρ ₀/ρ) has been established, and from it we derive some others, to provide a mechanism that is responsible for the cyclical reversible transitions. Also, the problems of singularities, black hole overproduction, and the second law of thermodynamics arising in oscillating universe models are conceptually resolved.