Comparison of Squeezing Prediction Methods: A Case Study on Nowsoud Tunnel

According to the International Society of Rock Mechanics, squeezing is a time dependent large deformation occurring during tunnel construction around the tunnel associated with creep caused by exceeding a limiting shear stress (Barla in ISRM News J 2:44–49, 1995). This research is conducted using a case study on the Nowsoud Tunnel, Iran. Being 14 km in length and 4.5 m in diameter, the tunnel is located in the western part of Iran near the Iraq border. Nowsoud tunnel, which was excavated using a double shield TBM, exhibited severe squeezing (with 8919 m) in its critical zone which resulted in excavation termination. In this research, the best approach for predicting squeezing among the recommended methods for reducing the damages caused by squeezing on TBM was determined. In this regard, approaches commonly used to predict squeezing are empirical, semi-empirical, and theoretical–analytical methods. Besides, these methods, numerical modeling is used to estimate convergence generated along the tunnel pathways, which is ultimately used to categorize squeezing. This paper compares squeezing prediction methods in 68 section of Nowsoud Tunnel. These 68 sections indicate that the empirical methods propose a general estimation/overview of squeezing. Among the semi-analytical approaches, the one proposed by Hoek and Marinos (Rock engineering in difficult rock conditions—soft rocks and karst, Taylor & Francis Group, London, pp 49–60, 2000) are compatible with the occurrence of squeezing in the critical zone. However, the degree of predicted squeezing is less than the real degree of squeezing in this zone. Based on the result of Aydan approach, 75 % of the tunnel sections are under squeezing condition. Theoretical–analytical approaches underestimate the possibility of squeezing in the critical zone. Barla?s approach (1995) demonstrated the possibility of squeezing in the critical zone with low intensity. The numerical computations in this paper were performed using Plaxis (version 8.5), a two-dimensional numerical program based on the finite element method. Plaxis results, classified by Hoek and Marinos (2000) method, show that 8800 m of the tunnel length is under the non-squeezing condition. According to all prediction methods, the squeezing zones depend on faulted zones, argillaceous limestone and shale formations such as J₁Kh, J₄Kh, J₅Kh, and Kgr. These formations were identified with a high quantity of shale and argillaceous limestone. Bedding of these geological formations is thin and their geotechnical properties are weaker than those of limestone formations. On the other hand, non-squeezing zones depend on limestone formations such as J₂Kh, J₃Kh, J₆Kh, K^abg, and Kbg. Moreover, all approaches predicted squeezing potential for the critical zone where TBM is jammed.     

Using uncertainty and sensitivity analysis for finding the best rainfall-runoff model in mountainous watersheds(Case study: the Navrood watershed in Iran)

In watersheds that have not sufficient meteorological and hydrometric data for simulating rainfall-runoff events, using geomorphologic and geomorphoclimatic characteristics of watershed is a conventional method for the simulation. A number of rainfall-runoff models utilize these characteristics such as Nash-IUH, Clark-IUH, Geomorphologic Instantaneous Unit Hydrograph(GIUH), Geomorphoclimatic Instantaneous Unit Hydrograph(GcIUH), GIUH-based Nash(GIUH-Nash) and GcIUH-based Clark(GcIUH-Clark). But all these models are not appropriate for mountainous watersheds. Therefore, the objective of this study is to select the best of them for the simulation. The procedure of this study is: a) selecting appropriate rainfall-runoff events for calibration and validation of six hybrid models, b) distinguishing the best model based on different performance criteria(Percentage Error in Volume(PEV); Percentage Error in Peak(PEP); Percentage Error in Time to Peak(PETP); Root Mean Square Error(RMSE) and Nash-Sutcliffe model efficiency coefficient(ENS)), c) Sensitivity analysis for determination of the most effective parameter at each model, d) Uncertainty determination of different parameters in each model and confirmation of the obtained results by application of the performance criteria. For application of this procedure, the Navrood watershed in the north of Iran as a mountainous watershed has been considered. The results showed that the ClarkIUH and GcIUH-Clark are suitable models for simulation of flood hydrographs, while other models cannot simulate flood hydrographs appropriately. The sensitivity analysis shows that the most sensitive parameters are the infiltration constant rate and time of concentration in the Clark-IUH model. Also, the most sensitive parameters include the infiltration constant rate and storage coefficient in the GcIUHClark model. The Clark-IUH and GcIUH-Clark models are more sensitive to their parameters. The Latin Hypercube Sampling(LHS) on Monte Carlo(MC) simulation method was used for evaluation of uncertainty of data in rainfall-runoff models. In this method 500 sets of data values are produced and then the peak discharge of flood hydrographs for each produced data set is simulated with rainfall-runoff models. The uncertainty of data changes the value of simulated peak discharge of flood hydrograph. The uncertainty analysis shows that the observed peak discharges of different rainfall-runoff events are within the range of values of simulated by the six hybrid rainfall-runoff models and IUH that inputs of these models were the produced data sets. The range of the produced peak discharge of flood hydrographs by the Clark-IUH and GcIUH-Clark models is wider than those of other models.     

Problematic Soft Soil Improvement with Both Polypropylene Fiber and Polyvinyl Acetate Resin

Several methods are used to improve mechanical properties of loose soils including rewetting, soil replacement, compaction control, chemical additives, moisture control, thermal methods, and more recently, discrete fibers. All the methods are applied to soft soil to increase load bearing capacity and to improve other properties such as prevention of erosion and dust generation. In the present study, a new method of soil improvement using both discrete polypropylene (PP) fibers and polyvinyl acetate (PVAc) is introduced. The method is applied to improve load bearing capacity of a problematic sandy soil in both dry and saturated states. Based on the results from CBR tests on various specimens, it has been revealed that the combination of PP fiber and PVAc resin with weight percentages of 0.1 and 0.6 %, respectively, had the optimum effect in increasing the CBR value in both saturated and dry soil specimens. It should be mentioned that this method has caused a great increase in the CBR value in the saturated soil.     

Morphological controls on the dynamics of carbonate landscapes under a mediterranean climate

The influence of the morphological setting on the denudation of carbonate landscapes and the respective contributions of mechanical and chemical weathering processes are still debated. We have addressed these questions by measuring ³⁶Cl concentrations in 40 samples from the Luberon mountain, SE France, to constrain the denudation of various landscape elements. We observe a clear contrast between the local denudation rates from the flat summit surface, clustered around 30 mm/ka, and the basin‐average denudation rates across the flanks, ranging from 100 to 200 mm/ka. This difference highlights the transient evolution of the range, whose topography is still adjusting to previous uplift events. Such a pattern also suggests that carbonate dissolution is not the only driver of denudation in this setting, which appears to be significantly controlled by slope‐dependent processes.     

Estimation of Tasuj aquifer response to main meteorological parameter variations under Shared Socioeconomic Pathways scenarios

Theoretical and Applied Climatology - In this study, statistical and soft-computing methods are compared in forecasting groundwater levels under Shared Socioeconomic Pathways (SSPs) SSP1-2.6,...     

Estimation of joint trace length probability distribution function in igneous,sedimentary, and metamorphic rocks

To predict the behavior of structures in and on jointed rock masses, it is necessary to characterize the geomechanical properties of joints and intact rock. Among geometry properties of joints, trace length has a vital importance, because it affects rock mass strength and controls the stability of the rock structures in jointed rock masses. Since joint length has a range of values, it is useful to have an understanding of the distribution of these values in order to predict how the extreme values may be compared to the values obtained from a small sample. For this purpose, three datasets of joint systems from nine exposures of igneous, metamorphic, and sedimentary rocks are studied. Joint trace length is one of the most difficult properties to measure accurately, but it may be possible to record other geometrical properties of exposed joints accurately; thereby, support vector machine (SVM) model is used to predict the joint trace length. SVM is a novel machine learning method, which is a powerful tool used to solve the problem characterized by small sample and non-linearity with a good generalization performance. Consequently, goodness-of-fit (GOF) tests were applied on these data. According to these GOF tests, the lognormal distribution was found to be the best probability distribution function for representing a joint trace length distribution.     

Plio-Quaternary stress states in NE Iran: Kopeh Dagh and Allah Dagh-Binalud mountain ranges

NE Iran, including the Kopeh Dagh and Allah Dagh-Binalud deformation domains, comprises the northeastern boundary of the Arabia–Eurasia collision zone. This study focuses on the evolution of the Plio-Quaternary tectonic regimes of northeast Iran. We present evidence for drastic temporal changes in the stress state by inversion of both geologically and seismically determined fault slip vectors. The inversions of fault kinematics data reveal distinct temporal changes in states of stress during the Plio-Quaternary (since ～ 5 Ma). The paleostress state is characterized by a regional transpressional tectonic regime with a mean N140 ± 10°E trending horizontal maximum stress axis (σ₁). The youngest (modern) state of stress shows two distinct strike-slip and compressional tectonic regimes with a regional mean of N030 ± 15°E trending horizontal σ₁. The change from the paleostress to modern stress states has occurred through an intermediate stress field characterized by a mean regional N trending σ₁. The inversion analysis of earthquake focal mechanisms reveals a homogeneous, transpressional tectonic regime with a regional N023 ± 5°E trending σ₁. The modern stress state, deduced from the youngest fault kinematics data, is in close agreement with the present-day stress state given by the inversions of earthquake focal mechanisms. According to our data and the deduced results, in northeast Iran, the Arabia–Eurasia convergence is taken up by strike-slip faulting along NE trending left-lateral and NNW trending right-lateral faults, as well as reverse to oblique-slip reverse faulting along NW trending faults. Such a structural assemblage is involved in a mechanically compatible and homogeneous modern stress field. This implies that no strain and/or stress partitioning or systematic block rotations have occurred in the Kopeh Dagh and Allah Dagh-Binalud deformation domains. The Plio-Quaternary stress changes documented in this paper call into question the extrapolation of the present-day seismic and GPS-derived deformation rates over geological time intervals encompassing tens of millions of years.     

Design of a sensitive air quality monitoring network using an integrated optimization approach

A new method is developed to design a multi-objective and multi-pollutant sensitive air quality monitoring network (AQMN) for an industrial district. A dispersion model is employed to estimate the ground level concentration of the air pollutants emitted from different emission sources. The primary objective of AQMN is providing the maximum information about the pollutant with respect to (1) maximum coverage area, (2) maximum detection of violations over ambient air standards and (3) sensitivity of monitoring stations to emission sources. Ant Colony Optimization algorithm (ACO) and Genetic Algorithm (GA) are adopted as the optimization tools to identify the optimal configuration of the monitoring network. The comparison between the results of ACO and GA shows that the performance of both algorithms is acceptable in finding the optimal configuration of AQMN. The application of the method to a network of existing refinery stacks indicates that three stations are suitable to cover the study area. The sensitivity of the three optimal station locations to emission sources is investigated and a database including the sensitivity of stations to each source is created.     

Analysis of Spatial Autocorrelation Patterns of Heavy and Super-Heavy Rainfall in Iran

Rainfall is a highly variable climatic element, and rainfall-related changes occur in spatial and temporal dimensions within a regional climate. The purpose of this study is to investigate the spatial autocorrelation changes of Iran's heavy and super-heavy rainfall over the past 40 years. For this purpose, the daily rainfall data of 664 meteorological stations between 1971 and 2011 are used. To analyze the changes in rainfall within a decade, geostatistical techniques like spatial autocorrelation analysis of hot spots, based on the Getis-Ord Gistatistic, are employed. Furthermore, programming features in MATLAB, Surfer, and GIS are used. The results indicate that the Caspian coast, the northwest and west of the western foothills of the Zagros Mountains of Iran, the inner regions of Iran, and southern parts of Southeast and Northeast Iran,have the highest likelihood of heavy and super-heavy rainfall. The spatial pattern of heavy rainfall shows that, despite its oscillation in different periods, the maximum positive spatial autocorrelation pattern of heavy rainfall includes areas of the west, northwest and west coast of the Caspian Sea. On the other hand, a negative spatial autocorrelation pattern of heavy rainfall is observed in central Iran and parts of the east, particularly in Zabul. Finally, it is found that patterns of super-heavy rainfall are similar to those of heavy rainfall.