Drought characterization using a new copula-based trivariate approach

Meteorological drought is a natural climatic phenomenon that occurs over various time scales and may cause significant economic, environmental and social damages. Three drought characteristics, namely duration, average severity and peak intensity, are important variables in water resources planning and decision making. This study presents a new method for construction of three-dimensional copulas to describe the joint distribution function of meteorological drought characteristics. Using the inference function for margins, the parameters for six types of copulas were tested to select the best-fitted copulas. According to the values of the log-likelihood function, Galambos, Frank and Clayton were the selected copula models to describe the dependence structure for pairs of duration–severity, severity–peak and duration–peak, respectively. Trivariate cumulative probability, conditional probability and drought return period were also investigated based on the derived copula-based joint distributions. The proposed model was evaluated over the observed data of a Qazvin synoptic station, and the results were compared with the empirical probabilities. For measuring the model accuracy, R ², root mean square error (RMSE) and the Nash–Sutcliffe efficiency (NSE) criteria were used. Results indicated that R ², RMSE and NSE were equal to 0.91, 0.098 and 0.668, respectively, which demonstrate sufficient accuracy of the proposed model. Drought probabilistic characteristics can provide useful information for water resource planning and management.     

An investigation of near surface wind speed trends in arid and semiarid regions of Iran

Wind plays an important role on the ecosystems and hydrological cycles besides other meteorological parameters such as temperature, precipitation, sunshine, and relative humidity. It strongly affects evapotranspiration, especially in arid and semiarid regions where there are serious problems in regard to water resource management. Evaluating the wind speed trend can provide good information for future agricultural planning. This study was conducted in order to investigate the wind speed trends over 24 synoptic meteorological stations located in arid and semiarid regions of Iran from 1975 to 2005. Near-surface wind speed was trended by nonparametric Mann–Kendall test spatially and temporally in three time scales (annual, seasonal, and monthly). Then, Sen’s slope estimator was used to determine the amount of the changes; furthermore, 10-year moving average low-pass filter was applied to show general trends. Finally, the smoothed time series derived from the mentioned filter were classified in three clusters for each time series and then mapped to show their spatial distribution pattern. Results showed insignificant and significant, increasing and decreasing trends during the surveyed time. Wind speeds in less than 50 % of stations changed statistically in all time scales, and in most cases, the frequency of the upward trends was more than that of downward ones. The spatial distribution of significant wind speed showed that the increase mostly occurred in eastern part. Clustering gave us the turning point around 1990. Clearly, when clusters were mapped, they indicated the same pattern as the Z value maps derived from Mann–Kendall test which meant that the outputs of the mentioned method confirmed the other one. As the wind speed trends in different stations likely to follow the previous evapotranspiration (ET₀) trend results in Iran, it confirms that wind speed was an effective parameter on ET₀, even though other parameters should be considered too.     

Comparing grain size distribution of sediment and original soil under raindrop detachment and raindrop-induced and flow transport mechanism

The grain size distribution (GSD) of sediment in comparison with the original soil GSD is discussed under different slopes (5, 15 and 25%) and rainfall intensities (30, 60 and 90 mm h^–1 with respective duration of 30, 15 and 10 min) but identical runoff (15 mm). The sediment quantification was carried out by raindrop-induced flow transport (RIFT) or/and transport by flow (FT) using a rainfall simulator and a 6 × 1 m² erosion plot and a silt loam. The results show a high degree of enrichment for size classes of 2–4 and 4–8 μm and a high degree of depletion for size classes of >63 μm under different slopes and rainfall intensities. In addition, the results show that the experimental enrichment ratio (ER) for particle size <16 μm under different slopes and rainfall intensities was greater than 1, while the ER for particle size >32 μm was less than 1.     

Sedimentary facies and sequence stratigraphy of the Asmari Formation at Chaman-Bolbol, Zagros Basin, Iran

The Oligocene–Miocene Asmari Formation of the Zagros Basin is a thick sequence of shallow water carbonate. In the study area, it is subdivided into 14 microfacies that are distinguished on the basis of their depositional textures, petrographic analysis and fauna. Based on the paleoecology and lithology, four distinct depositional settings can be recognized: tidal flat, lagoon, barrier, and open marine. The Asmari Formation represents sedimentation on a carbonate ramp. In the inner ramp, the most abundant lithofacies are medium grained wackestone–packstone with imperforated foraminifera. The middle ramp is represented by packstone–grainstone to floatstone with a diverse assemblage of larger foraminifera with perforate wall, red algae, bryozoa, and echinoids. The outer ramp is dominated by argillaceous wackestone characterized by planktonic foraminifera and large and flat nummulitidae and lepidocyclinidae. Three third-order depositional sequences are recognized from deepening and shallowing trends in the depositional facies, changes in cycle stacking patterns, and sequence boundary features.     

Upper-mantle structure of the Baltic Shield below the Swedish National Seismological Network (SNSN) resolved by teleseismic tomography

Upper-mantle structure under the Baltic Shield is studied using non-linear high resolution teleseismic P -phase tomography. Observed relative arrival-time residuals from 52 teleseismic earthquakes recorded by the Swedish National Seismological Network (SNSN) are inverted to delineate the structure of the upper mantle. The network consists of 47 (currently working) three-component broad-band stations located in an area about 450 km wide and 1450 km long. In order to reduce complications due to possible significant three-dimensionality of Earth structure, events chosen for this study lay close to in-line with the long-axis of the array  (±30°)  . Results indicate P -wave velocity perturbations of ±3 per cent down to at least 470 km below the network. The size of the array allows inversion for structures even at greater depths, and lateral variations of velocity at depths of up to 680 km appear to be resolved. Below the central part of the array (60°–64° N), where ray coverage is best, the data reveals a large region of relatively low velocity at depths of over about 300 km. At depths less than about 250–300 km, the models include a number of features, including an apparent slab-like structure dipping gently towards the north.     

Analysis of pervasive precipitation in similar gradient areas of Iran

Precipitation is a major climatic element with high spatial variations. Temporal and spatial variations may differ in large and small scales. It is, therefore, of utmost importance to study areas with similar gradients in terms of precipitation patterns in order to shed light on the complexities of precipitation variations. In the present study, attempts were made to identify areas with similar gradients experiencing the same precipitation pattern over a 50-year period (1964–2013). To this end, data were collected from synoptic stations in Iran in two phases (i.e., 1434 stations in the first phase and 673 stations in the second one). Alexanderson’s technique was adopted to examine sudden changes in precipitation patterns. The results showed that five regions with similar gradients could be identified in terms of precipitation patterns: negative and high variations, negative and moderate variations, positive and high variations, positive and moderate variations, and little or no variations. The distribution of such regions indicated that the regions with positive trends experienced more annual variations and had further spatial distribution. Furthermore, the findings revealed that the regions with negative precipitation patterns experienced more sudden changes in comparison with those with positive precipitation patterns. Additionally, more variations were observed in the precipitation patterns in recent years.     

Jurassic igneous rocks of the central Sanandaj–Sirjan zone (Iran) mark a propagating continental rift,not a magmatic arc

Jurassic igneous bodies of the Sanandaj–Sirjan zone (SaSZ) in SW Iran are generally considered as a magmatic arc but critical evaluation of modern geochronology, geochemistry and radiogenic isotopes challenges this conclusion. There is no evidence for sustained igneous activity along the ~1,200 km long SaSZ, as expected for a convergent plate margin; instead activity was brief at most sites and propagated NW at ~20 mm/a. Jurassic igneous rocks define a bimodal suite of gabbro‐diorite and granite. Chemical and isotopic compositions of mafic rocks indicate subcontinental lithospheric mantle sources that mostly lacked subduction‐related modifications. The arc‐like features of S‐type granites reflect massive involvement of Cadomian crust and younger sediments to generate felsic melts in response to mafic intrusions. We conclude that Jurassic SaSZ igneous activity occurred in a continental rift, not an arc. SaSZ igneous rocks do not indicate that subduction along the SW margin of Eurasia began in Jurassic time.     

Predictive value of Keetch-Byram Drought Index for cereal yields in a semi-arid environment

The study examines the potential of urban roofs to reduce the urban heat island (UHI) effect by changing their reflectivity and implementing vegetation (green roofs) using the example of the City of Vienna. The urban modelling simulations are performed based on high-resolution orography and land use data, climatological observations, surface albedo values from satellite imagery and registry of the green roof potential in Vienna. The modelling results show that a moderate increase in reflectivity of roofs (up to 0.45) reduces the mean summer temperatures in the densely built-up environment by approximately 0.25 °C. Applying high reflectivity materials (roof albedo up to 0.7) leads to average cooling in densely built-up area of approximately 0.5 °C. The green roofs yield a heat load reduction in similar order of magnitude as the high reflectivity materials. However, only 45 % of roof area in Vienna is suitable for greening and the green roof potential mostly applies to industrial areas in city outskirts and is therefore not sufficient for substantial reduction of the UHI effect, particularly in the city centre which has the highest heat load. The strongest cooling effect can be achieved by combining the green roofs with high reflectivity materials. In this case, using 50 or 100 % of the green roof potential and applying high reflectivity materials on the remaining surfaces have a similar cooling effect.     

A GIS‐based multi‐criteria seismic vulnerability assessment using the integration of granular computing rule extraction and artificial neural networks

This study proposes multi‐criteria group decision‐making to address seismic physical vulnerability assessment. Granular computing rule extraction is combined with a feed forward artificial neural network to form a classifier capable of training a neural network on the basis of the rules provided by granular computing. It provides a transparent structure despite the traditional multi‐layer neural networks. It also allows the classifier to be applied on a set of rules for each incoming pattern. Drawbacks of original granular computing (GrC) are covered, where some input patterns remained unclassified. The study was applied to classify seismic vulnerability of the statistical units of the city of Tehran, Iran. Slope, seismic intensity, height and age of the buildings were effective parameters. Experts ranked 150 randomly selected sample statistical units with respect to their degree of seismic physical vulnerability. Inconsistency of the experts' judgments was investigated using the induced ordered weighted averaging (IOWA) operator. Fifty‐five classification rules were extracted on which a neural network was based. An overall accuracy of 88%, κ = 0.85 and R² = 0.89 was achieved. A comparison with previously implemented methodologies proved the proposed method to be the most accurate solution to the seismic physical vulnerability of Tehran.     

A Geometric Computational Model for Calculation of Longwall Face Effect on Gate Roadways

In this paper a geometric computational model (GCM) has been developed for calculating the effect of longwall face on the extension of excavation-damaged zone (EDZ) above the gate roadways (main and tail gates), considering the advance longwall mining method. In this model, the stability of gate roadways are investigated based on loading effects due to EDZ and caving zone (CZ) above the longwall face, which can extend the EDZ size. The structure of GCM depends on four important factors: (1) geomechanical properties of hanging wall, (2) dip and thickness of coal seam, (3) CZ characteristics, and (4) pillar width. The investigations demonstrated that the extension of EDZ is a function of pillar width. Considering the effect of pillar width, new mathematical relationships were presented to calculate the face influence coefficient and characteristics of extended EDZ. Furthermore, taking GCM into account, a computational algorithm for stability analysis of gate roadways was suggested. Validation was carried out through instrumentation and monitoring results of a longwall face at Parvade-2 coal mine in Tabas, Iran, demonstrating good agreement between the new model and measured results. Finally, a sensitivity analysis was carried out on the effect of pillar width, bearing capacity of support system and coal seam dip.