Assessing desertification by using soil indices

Desertification generally refers to land degradation in arid, semiarid, and dry semi-humid climatic zones. It involves five principal processes: vegetation degradation, water erosion, wind erosion, salinization and waterlogging, and soil crusting and compaction. The aim of this study is assessing desertification using soil criteria. For this purpose, nine indices including sodium absorption ratio (SAR), soil gypsum percentage, soil texture, the content of HCO₃ ^?, the percentage of the organic matter, electrical conductivity (EC), pH, the content of the soil sodium, and chloride were used. The soil samples were taken in the north of Zayandeh-Rood River in Isfahan province of Iran, using soil data randomly sampled in a depth of 0–20 cm. After assessing the normality of the samples using Kolmogorov-Smirnov test, indices were imported into GIS environment and interpolated with IDW and normal and discrete kriging methods for delineating soil characteristics maps based on MEDALUS model. In this model, the data were firstly changed from 100 to 200. Thus 100 and 200 are estimated as the best and worst quality, respectively. Then the final map of soil criteria has been created by geometric mean of its indicators. The results showed that the maximum area is related to the medium class of desertification and is equal to 44,746 ha. The areas of severe and very severe classes of desertification are equal to 30,949 and 351 ha, respectively. The results also revealed that the indices of the organic matter and soil gypsum percentage are the most influential indices which affect desertification phenomenon.     

Probabilistic seismic hazard assessment in the northeastern part of Algeria

This work involves updating the evaluation of seismic hazard in Northeast Algeria by a probabilistic approach. This reassessment attempts to resolve inconsistencies between seismic zoning in regional building codes and is further motivated by the need to refine the input data that are used to evaluate seismic hazard scenarios. We adopted a seismotectonic model that accounts for differences in interpretations of regional seismicity. We then performed a probabilistic assessment of regional seismic hazard in Northeast Algeria. Based on a homogeneous earthquake catalog and geological and seismotectonic data gathered in the first part of the study, a seismotectonic zoning map was created and seven risk areas were identified. For each area, peak ground acceleration hazard maps were produced. Details of the calculations are provided, including hazard curves at periods of 0.1, 0.2, 0.33, 0.5, 1.0, and 2.0 s and uniform hazard spectra at urban locations in the area, including Sétif, Constantine, Kherrata, Bejaia, and Jijel.     

Paleoenvironment of geomorphic surfaces of an alluvial fan in the eastern Isfahan,Iran, in the light of micromorphology and clay mineralogy

Alluvial fans are one of the most important landforms in geomorphological and paloenvironmental studies. The objective of this study was the application of clay mineral assemblages and micromorphological properties of the studied paleosols in the geomorphic surfaces of an alluvial fan in the eastern Isfahan as proxies for paleoenvironmental and paleoclimatic changes. Micromorphology, X-ray diffraction, and scanning electron microscopy approaches were used to study the representative pedons. The results indicated that the illuviation process in the calcareous soils of the arid regions of the eastern Isfahan was probably in response to Quaternary moist conditions. There was no significant difference between clay coating properties of the studied relict and buried paleosols. Clay mineralogical study suggested that kaolinite and illite were inherited from the parent materials, while smectite and palygorskite were formed in the soil environment. Paleoargillic horizon was characterized by smectite and calcic (especially the calcrete) horizons were dominated by palygorskite. Palygorskite was accumulated by both neoformation and illuviation processes. High clay content, high intensity of smectite peak, and activity of the illuviation process in paleoargillic horizon demonstrated the seasonality of climate (rainfall) even in the moist periods of Quaternary in Central Iran. Clay mineralogical assemblages suggested a trend of increasing environmental aridity in the study area. This study, therefore, highlighted the role of clay mineralogical investigations in arid lands’ geomorphological and paleoenvironmental researches.     

Discussion of “Interfacing the geographic information system,remote sensing,and the soil conservation service-curve number method to estimate curve number and runoff volume in the ASIR region of Saudi Arabia” by Fawzi S. Mohammad,Jan Adamowski

The soil conservation service curve number (SCS-CN) method is one of the most commonly used methods to compute the direct runoff from a rainfall event. Since the method was established, numerous researches were undertaken to improve the method through accurate estimation of its parameter and especially the curve number (CN). However, the essence of the SCS method, as an event-based Hortonian mechanism method, remained unchanged. The main assumption of the method related to the rainfall input is that the rainfall is continuous in time and uniform over the watershed. Mohammad and Adamowski (2015) paper apparently used the SCS method to estimate the annual runoff using the annual rainfall as one cumulative rainfall input value, which is a violation of the event-based principle of the method and of the assumption of the continuity of the rainfall event.To re-estimate the average annual runoff more realistically for the Asir region, Saudi Arabia, daily rainfall data from 14 rainfall stations are used for calculating the resulting runoff depths, on a daily event-by-event rainfall basis, throughout the whole simulation period. The resulting runoff depths are added for each year, and the total cumulative annual runoff values for each year are averaged to get the average annual runoff. The runoff values based on the previously mentioned procedure are an upper limit of the actual average annual runoff as the underlying SCS equations discard evaporation and similar long-term losses. Nevertheless, the average runoff values obtained in the discussion paper are an order of magnitude (at least five to tenfold) lower than the ones of the original paper. An equation is proposed to obtain a more realistic estimate of the average annual runoff, to be used with the average annual rainfall as an input, if the annual value is the only available rainfall information.     

Regional modeling of dust storm of February 8, 2015 in the southwest of Iran

According to World Meteorological Organization report in 2015, the southwest of Iran has become one of the dust sources in the region. And the objective of this research is to study the dust storms originating in this region. For this purpose, based on the weather data of 14 stations, the dust storms of the region were investigated, and the dust storm of February 7, 2015, was selected due to its very high concentration of dust particles (66 times the normal values). For the analysis of source areas and storm paths, the FNL data was used. The regional models of Navy Aerosol Analysis and Prediction System (NAAPS), Barcelona Supercomputing Centre-Dust Regional Atmospheric Model 8b (DREAM 8b), Non-hydrostatic Multiscale Model Barcelona Supercomputing Center (NMMB/BSC), and hybrid single-particle Lagrangian integrated trajectory (HYSPLIT) were used to study and analyze the selected storm. The results showed that the dust event in February 8, 2016, has been the result of the polar front jet stream (PJF) caused by western immigrant system that had been over the Sahara in Africa, the deserts of Iraq, Syria, Saudi Arabia, and finally southwest of Iran in making the extreme dust event. According to the Moderate Resolution Imaging Spectroradiometer data and point models of NAAPS, optical depth was very high. The DREAM 8b and NMMB/BSC models confirmed the impact of the local factors and closeness to the dust source regions. The backward tracking of the model with the HYSPLIT model showed three tracks transporting the dust particles to the region. This software also showed that the dust particles occupied an atmospheric tunnel of 1.5 km in diameter.     

Applying an optimized proxy-based workflow for fast history matching

History matching is still one of the main challenging parts of reservoir study especially in giant brown oil fields with lots of wells. In these cases, history matching with conventional manual technique needs many runs and takes months to get a match. In this work, an innovative approach was suggested for fast history matching in a real brown field. The workflow was employed based on an optimized proxy model for history matching of a field consisting of 14 active wells with multiple responses (which are production rate and pressure data) in the south part of Iran. The main important features of the proposed algorithm were defining a proxy model which is response surface method in which 21 model parameters were incorporated based on cubic centered face method. The proxy model was then optimized by one of the most famous algorithms which is genetic algorithm. Proxy model was successfully performed using 256 samples leading into p- value of 0.531 and R ² of 0.91 dataset. As a result, the proposed workflow and algorithm showed good and acceptable results for history matching of studied real model.     

Assessment of proposed approaches for bathymetry calculations using multispectral satellite images in shallow coastal/lake areas: a comparison of five models

Bathymetric information for shallow coastal/lake areas is essential for hydrological engineering applications such as sedimentary processes and coastal studies. Remotely sensed imagery is considered a time-effective, low-cost, and wide-coverage solution for bathymetric measurements. This study assesses the performance of three proposed empirical models for bathymetry calculations in three different areas: Alexandria port, Egypt, as an example of a low-turbidity deep water area with silt-sand bottom cover and a depth range of 10.5 m; the Lake Nubia entrance zone, Sudan, which is a highly turbid, unstable, clay bottom area with water depths to 6 m; and Shiraho, Ishigaki Island, Japan, a coral reef area with varied depths ranging up to 14 m. The proposed models are the ensemble regression tree-fitting algorithm using bagging (BAG), ensemble regression tree-fitting algorithm of least squares boosting (LSB), and support vector regression algorithm (SVR). Data from Landsat 8 and Spot 6 satellite images were used to assess the performance of the proposed models. The three models were used to obtain bathymetric maps using the reflectance of green, red, blue/red, and green/red band ratios. The results were compared with corresponding results yielded by two conventional empirical methods, the neural network (NN) and the Lyzenga generalised linear model (GLM). Compared with echosounder data, BAG, LSB, and SVR results demonstrate higher accuracy ranges from 0.04 to 0.35 m more than Lyzenga GLM. The BAG algorithm, producing the most accurate results, proved to be the preferable algorithm for bathymetry calculation.     

An alternative direct method towards mean dynamic topography computations

An alternative “direct method” to “mean dynamic topography” (MDT) computations using satellite altimetry-derived “mean sea surface” (MSS) and “global geopotential model” (GGM), without direct application of the geoid, is devised. The developed approach, which is based on derivation of an equipotential surface of the gravity field of the Earth that fits to global MSS in least squares sense, is formulated via a constrained optimization problem. The validity of our method is numerically tested by computing a global MDT model based on DNSC08 MSS model and EGM2008 GGM as input data.