Real surface conductivity component as indicator for the hydraulic conductivity

Estimation of hydraulic conductivity from surface resistivity measurements is one of the most difficult and challenging hydrogeophysical targets. The promising side of this relation is the analogy between electric current flow and water flow, whereas the grand ambiguity is the non-dimensionality between both two quantities. Imaginary surface conductivity component is used recently to deduce the hydraulic conductivity via complex resistivity measurements. Since there are similar properties between imaginary (out-of-phase) and real (in-phase) surface conductivity components, the latter is used in this paper to predict the hydraulic conductivity. Two mathematical parameters were determined to express the electrical equivalent of hydraulic conductivity in sand and clay systems based on the mode of electrical double-layer formation in both systems. The reliability of the proposed method is tested through applying on two datasets representing sand and clay systems. The first dataset is a clean sand and gravel aquifer in the Keritis basin in Chania, Crete, Greece. The second is mostly clayey sand aquifer in Wadi El-Assuity, Egypt. Application of the present approach in these two cases resulted in promising nearly identical values with the measured hydraulic conductivity via pumping test or geometric hydraulic conductivity via grain size analysis.     

Spatial variability of soil erodibility and its correlation with soil properties in semi-arid mountainous watershed,Saudi Arabia

Soil erodibility values are best estimated from long-term direct measurements on runoff-plots; however, in lack of field tests, these values can be estimated using relationships based on physico-chemical soil properties. The study objective was to assess the erodibility and its correlation with soil properties. The average erodibility value was estimated 0.043 t ha h ha^?1 MJ^?1 mm^?1. The areas with heavy textured soil and low organic matter content had the lowest values of erodibility. The erodibility decreases as the sand content increases, whereas silt showed a positive correlation. The erodibility factors and its relation to soil properties were evaluated using multiple regression analysis. Results revealed that sand and organic matter content of soil combinedly explained 78% of variation. Altitudinal increases also seem to affect the soil texture. This study has demonstrated that soil properties and erodibility values can be used as assistance for soil conservation practices and modelling of landscape processes.     

Structures and microstructures related to steady-state mantle flow in the Neyriz ophiolite, Iran

Structural analysis in the well-exposed and well-preserved Neyriz ophiolite suggests that it is a relatively undisturbed piece of oceanic lithosphere. Detailed structural mapping of high-T deformation mantle flow revealed the presence of three elliptical shaped diapirs trending NW–SE. These diapirs are characterized by vertical mantle foliations associated with vertical plunging stretching lineations, which progressively incline toward parallelism with the gently NE-dipping Moho. The NW–SE direction of asthenospheric upwelling of diapirs is parallel with the orientations of the vertical sheeted dike complex. This suggests that the Neyriz ophiolite was created by two NW–SE palaeo-ridge axes. These palaeo-ridges are classified as fast-spreading ridges. These ridges are segmented by a dextral palaeo-transform fracture zone. This is consistent with fast-spreading ridges. Comparison between the Neyriz and Oman ophiolites reveals that they show similar characteristics. Most of the Oman palaeo-ridge systems are oriented NW–SE and NNW–SE. They also display similar sheeted dike complex orientations and crustal thickness variations. These two ophiolites originally were part of the Neo-Tethyan oceanic lithosphere and afterwards were separated by the Oman line during continental collision of the Iranian micro-continent and Afro-Arabian continent.     

Mixed model for interoccurrence times of earthquakes based on the expectation-maximization algorithm

In this paper, the goodness-of-fit test based on a convex combination of Akaike and Bayesian information criteria is used to explain the features of interoccurrence times of earthquakes. By analyzing the seismic catalog of Iran for different tectonic settings, we have found that the probability distributions of time intervals between successive earthquakes can be described by the generalized normal distribution. This indicates that the sequence of successive earthquakes is not a Poisson process. It is found that by decreasing the threshold magnitude, the interoccurrence time distribution changes from the generalized normal distribution to the gamma distribution in some seismotectonic regions. As a new insight, the probability distribution of time intervals between earthquakes is described as a mixture distribution via the expectation-maximization algorithm.     

Coastal boulder deposit as evidence of an ocean-wide prehistoric tsunami originated on the Atacama Desert coast (northern Chile)

A Late Holocene cliff-top deposit of large boulders well above the limits of modern storm waves is described from the southern coast of the Atacama Desert (northern Chile). The largest moved boulder weighs >40 t and field data point to a flood height >18·5 m above high tide level and an inland penetration greater than 284 m from the cliff edge. The minimum flow velocity needed for particle entrainment was estimated as 10·1 ms⁻¹ and the most likely processes of sediment deposition for different boulders were deduced. The boulder distribution, sorting and orientation of imbricated debris, together with the significant wave height of extreme storms reported and the occurrence of interplate earthquakes in the study area indicate that the deposit records a single event, interpreted here as a tsunami wave train rather than exceptional storm waves. The boulder field was dated to between the 13th and the 16th Centuries ce and possibly correlates with the 1420 Oei orphan tsunami, that affected the eastern coast of Japan. A magnitude of 8·8 to 9·4 has been estimated for the earthquake, which may be one of the larger events of a super-cycle of earthquakes in the southern Atacama Desert. These cycle-ending earthquakes involve large rupture areas (lengths in excess of 600 km) and highly destructive ocean-wide tsunamigenic events.     

Anthropogenic Decline of Ancient,Sustainable Water Systems: Qanats

Qanat is an ancient underground structure to abstract groundwater without the need for external energy. A recognized world heritage, Qanat has enabled civilization in arid and semi-arid regions that lack perennial surface water resources. These important structures, however, have faced significant challenges in recent decades due to increasing anthropogenic pressures. This study uses remote sensing to investigate land-use changes and the loss of 15,983 Qanat shafts in the Mashhad plain, northeast of Iran, during the past six decades. This entails obtaining a rare aerial imagery from 1961, as well as recent satellite imagery, over a region with the highest density of Qanats in Iran, the birthplace of Qanat. Results showed that only 5.59% of the Qanat shafts in 1961 remained intact in 2021. The most prominent Qanat-impacting land-use changes were agriculture and urban areas, that accounted for 42.93 and 31.81% Qanat shaft destruction in the study area, respectively. This study also showed that groundwater table decline, demographic changes, and reduction in the appeal of working in the Qanat maintenance and construction industry among the new generation are existential threats to Qanats, and may result in the demise of these ancient structures in the future. Findings of this study can be used for urban planning in arid and semi-arid areas with the aim of protecting these historic water structures.     

Pattern recognition of major asperities using local recurrence time in Alborz Mountains, Northern Iran

In this study, seismic data recorded during the period 01/01/1996 to 09/01/2009 has been used to evaluate the seismic hazard potential along the Alborz region, Northern Iran. The technique of mapping local recurrence time, T _L, is used to map major asperities, which are considered as the areas with maximum hazard. We calculated T _L from a and b values which are in turn derived from the frequency–magnitude relation constants within a radius of 30 km about every corner point of a 10-km spacing grid. Since b value is inversely related to applied stress, the areas with lowest b values and/or shortest T _L are interpreted to locate the asperities or the areas of maximum seismic hazard. To test this method, we computed T _L map using seismic catalogues before and after the 2004 Baladeh earthquake of M _w 6.2. The local recurrence time map before the earthquake shows anomalously short T _L in the epicentral region of the Baladeh earthquake a decade before its occurrence. The T _L map after the earthquake indicates that this large event has redistributed the applied stress in the Alborz region. The microseismicity of the region after the Baladeh earthquake, however, suggests that there are two anomalies in T _L map positioned in Alborz. The places where these anomalies are observed can be considered as the areas with maximum seismic hazard for future large earthquake in the Alborz region.     

Computing derivatives of a gravity potential by using automatic differentiation

A new method, based on automatic differentiation technique, has been proposed in this paper to compute the derivatives of the gravity potential. Using this method we can obtain derivatives up to any order. Instead of explicit expressions of the derivatives we use an iterative scheme to simultaneously compute the value of all the desired derivatives. The algorithm here presented can be easily parallelized by using OpenMP with the consequent improvement in CPU-time efficiency.     

Modeling of the radial and transverse components of the albedo function

In this work, we shall confine ourselves to solve analytically the integrals called , for the two components of albedo radiation pressure on an Earth’s satellite. When the Earth’s albedo is variable, as far as we know, this case has not been dealt with analytically. We shall solve these two integrals when the satellite’ horizon is illuminated and when the sun lies on the equator. This approach will facilitate the evaluation of the mentioned two equations. We will also compare our results with previous works obtained numerically.