Non-stationary extreme value analysis in a changing climate

This paper introduces a framework for estimating stationary and non-stationary return levels, return periods, and risks of climatic extremes using Bayesian inference. This framework is implemented in the Non-stationary Extreme Value Analysis (NEVA) software package, explicitly designed to facilitate analysis of extremes in the geosciences. In a Bayesian approach, NEVA estimates the extreme value parameters with a Differential Evolution Markov Chain (DE-MC) approach for global optimization over the parameter space. NEVA includes posterior probability intervals (uncertainty bounds) of estimated return levels through Bayesian inference, with its inherent advantages in uncertainty quantification. The software presents the results of non-stationary extreme value analysis using various exceedance probability methods. We evaluate both stationary and non-stationary components of the package for a case study consisting of annual temperature maxima for a gridded global temperature dataset. The results show that NEVA can reliably describe extremes and their return levels.     

Global trends and patterns of drought from space

This paper analyzes changes in areas under droughts over the past three decades and alters our understanding of how amplitude and frequency of droughts differ in the Southern Hemisphere (SH) and Northern Hemisphere (NH). Unlike most previous global-scale studies that have been based on climate models, this study is based on satellite gauge-adjusted precipitation observations. Here, we show that droughts in terms of both amplitude and frequency are more variable over land in the SH than in the NH. The results reveal no significant trend in the areas under drought over land in the past three decades. However, after investigating land in the NH and the SH separately, the results exhibit a significant positive trend in the area under drought over land in the SH, while no significant trend is observed over land in the NH. We investigate the spatial patterns of the wetness and dryness over the past three decades, and we show that several regions, such as the southwestern United States, Texas, parts of the Amazon, the Horn of Africa, northern India, and parts of the Mediterranean region, exhibit a significant drying trend. The global trend maps indicate that central Africa, parts of southwest Asia (e.g., Thailand, Taiwan), Central America, northern Australia, and parts of eastern Europe show a wetting trend during the same time span. The results of this satellite-based study disagree with several model-based studies which indicate that droughts have been increasing over land. On the other hand, our findings concur with some of the observation-based studies.     

Efficiency analysis of open trench for impact pile driving through a single-variable method

Abstract

The ground vibrations during pile driving operation have a drastic potential to undermine the surrounding structures both in land and reclaimed land. Particularly, reclaimed land necessitates ample application of pile driving due to the weak land condition. To prevent the structural damage, attenuation of the ground vibrations to an allowable level through active isolation of circular open trench is the scope of this study. In this research, finite element simulations of continuous impact pile driving process from the ground surface was executed with particular attention to the pile-soil interaction, and thereby, the efficiency of open trench application in attenuation of the unsafe distance of different structures was surveyed using the vibration sensitivity degree. Regarding the crucial parameters of an open trench (depth, width, and location), it was concluded that a sufficient high depth can attenuate the unsafe distance up to 68%, the trench width variations are less effective, and an average pile-trench distance is the most efficient option. The excavation volume was also concluded as another crucial parameter in open trench design which takes all three parameters into account. The trench depth equal to the pile’s maximum critical depth of vibration was inferred for an optimum design.     

Geochemistry and zircon U-Pb geochronology of Aligoodarz granitoid complex, Sanandaj-Sirjan Zone, Iran

The Aligoodarz granitoid complex (AGC) is located in the Sanandaj-Sirjan Zone (SSZ), western Iran and consists of quartz-diorites, granodiorites and subordinate granites. Whole rock major and trace element data mostly define linear trends on Harker diagrams suggesting a cogenetic origin of the different rock types. (⁸⁷Sr/⁸⁶Sr)_i and εNd_t ratios are in the ranges 0.7074-0.7110 and −3.56 to −5.50, respectively. The trace elements and Sr-Nd isotopic composition suggest that the granitoids from the AGC are similar to crustal derived I-type granitoids of continental arcs. The whole rock suite was produced by assimilation and fractional crystallization starting from a melt with intermediate composition likely possessing a mantle component. In situ zircon U-Pb data on the granites with LA-ICP-MS yield a crystallization age of ∼165 Ma. Inherited grains spanning in age from ∼180 Ma up to 2027 Ma were also found and confirm that assimilation of country rock has occurred.Chemical and chronological data on the AGC were compared with those available for other granitoid complexes of the central SSZ (e.g., Dehno, Boroujerd and Alvand). The comparison reveals that in spite of the different origins that have been proposed, all these granitoid complexes are likely genetically related. They share many chemical features and are derived from crustal melts with minor differences. Alvand granites have the most peculiar compositions most likely related to the presence of abundant pelitic component. All these intrusions are coeval and reveal the presence of an extensive magmatic activity in the central sector of the SSZ during middle Jurassic.     

Forward induced seismic hazard assessment: application to a synthetic seismicity catalogue from hydraulic stimulation modelling

The M _w 3.2-induced seismic event in 2006 due to fluid injection at the Basel geothermal site in Switzerland was the starting point for an ongoing discussion in Europe on the potential risk of hydraulic stimulation in general. In particular, further development of mitigation strategies of induced seismic events of economic concern became a hot topic in geosciences and geoengineering. Here, we present a workflow to assess the hazard of induced seismicity in terms of occurrence rate of induced seismic events. The workflow is called Forward Induced Seismic Hazard Assessment (FISHA) as it combines the results of forward hydromechanical-numerical models with methods of time-dependent probabilistic seismic hazard assessment. To exemplify FISHA, we use simulations of four different fluid injection types with various injection parameters, i.e. injection rate, duration and style of injection. The hydromechanical-numerical model applied in this study represents a geothermal reservoir with preexisting fractures where a routine of viscous fluid flow in porous media is implemented from which flow and pressure driven failures of rock matrix and preexisting fractures are simulated, and corresponding seismic moment magnitudes are computed. The resulting synthetic catalogues of induced seismicity, including event location, occurrence time and magnitude, are used to calibrate the magnitude completeness M _c and the parameters a and b of the frequency-magnitude relation. These are used to estimate the time-dependent occurrence rate of induced seismic events for each fluid injection scenario. In contrast to other mitigation strategies that rely on real-time data or already obtained catalogues, we can perform various synthetic experiments with the same initial conditions. Thus, the advantage of FISHA is that it can quantify hazard from numerical experiments and recommend a priori a stimulation type that lowers the occurrence rate of induced seismic events. The FISHA workflow is rather general and not limited to the hydromechanical-numerical model used in this study and can therefore be applied to other fluid injection models.     

Assessment of geomorphological bank evolution of the alluvial threshold rivers based on entropy concept parameters

The complex stream bank profiles in alluvial channels and rivers that are formed after reaching equilibrium has been a popular topic of research for many geomorphologists and river engineers. The entropy theory has recently been successfully applied to this problem. However, the existing methods restrict the further application of the entropy parameter to determine the cross-section slope of the river banks. To solve this limitation, we introduce a novel approach in the extraction of the equation based on the calculation of the entropy parameter (λ) and the transverse slope of the bank profile at threshold channel conditions. The effects of different hydraulic and geometric parameters are evaluated on a variation of the entropy parameter. Sensitivity analysis on the parameters affecting the entropy parameter shows that the most effective parameter on the λ-slope multiplier is the maximum slope of the bank profile and the dimensionless lateral distance of the river banks.     

An effective stress-based DSC model for predicting the coefficient of lateral soil pressure in unsaturated soils

Acta Geotechnica - In this study, an analytical model is developed to establish a framework for predicting the coefficient of lateral soil pressure in unsaturated soils. To this end, the disturbed...     

APG:A novel python-based ArcGIS toolbox to generate absence-datasets for geospatial studies

One important step in binary modeling of environmental problems is the generation of absence-datasets that are traditionally generated by random sampling and can undermine the quality of outputs.To solve this problem,this study develops the Absence Point Generation(APG)toolbox which is a Python-based ArcGIS toolbox for automated construction of absence-datasets for geospatial studies.The APG employs a frequency ratio analysis of four commonly used and important driving factors such as altitude,slope degree,topographic wetness index,and distance from rivers,and considers the presence locations buffer and density layers to define the low potential or susceptibility zones where absence-datasets are gener-ated.To test the APG toolbox,we applied two benchmark algorithms of random forest(RF)and boosted regression trees(BRT)in a case study to investigate groundwater potential using three absence datasets i.e.,the APG,random,and selection of absence samples(SAS)toolbox.The BRT-APG and RF-APG had the area under receiver operating curve(AUC)values of 0.947 and 0.942,while BRT and RF had weaker per-formances with the SAS and Random datasets.This effect resulted in AUC improvements for BRT and RF by 7.2,and 9.7％from the Random dataset,and AUC improvements for BRT and RF by 6.1,and 5.4％from the SAS dataset,respectively.The APG also impacted the importance of the input factors and the pattern of the groundwater potential maps,which proves the importance of absence points in environmental bin-ary issues.The proposed APG toolbox could be easily applied in other environmental hazards such as landslides,floods,and gully erosion,and land subsidence.     

Earthquake Stability Analysis of Rock Slopes: a Case Study

Stability analysis of Surabhi landslide in the Dehradun and Tehri districts of Uttaranchal located in Mussoorie, India, has been simulated numerically using the distinct element method focusing on the weak zones (fracture). This is an active landslide on the main road toward the town centre, which was triggered after rainfall in July–August 1998. Understanding the behaviour of this landslide will be helpful for planning and implementing mitigation measures. The first stage of the study includes the total area of the landslide. The area identified as the zone of detachment is considered the most vulnerable part of the landslide. Ingress of water and increased pore pressures result in reduced mobilized effective frictional resistance, causing the top layer of the zone of detachment to start moving. The corresponding total volume of rock mass that is potentially unstable is estimated to 11.58 million m³. The second stage of this study includes a 2D model focussing only on the zone of detachment. The result of the analyses including both static and dynamic loading indicates that most of the total displacement observed in the slide model is due to the zone of detachment. The discontinuum modelling in the present study gives reasonable agreement with actual observations and has improved understanding of the stability of the slide slope.