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.     

Wave data assimilation using a hybrid approach in the Persian Gulf

The main goal of this study is to develop an efficient approach for the assimilation of the hindcasted wave parameters in the Persian Gulf. Hence, the third generation SWAN model was employed for wave modeling forced by the 6-h ECMWF wind data with a resolution of 0.5°. In situ wave measurements at two stations were utilized to evaluate the assimilation approaches. It was found that since the model errors are not the same for wave height and period, adaptation of model parameter does not result in simultaneous and comprehensive improvement of them. Therefore, an approach based on the error prediction and updating of output variables was employed to modify wave height and period. In this approach, artificial neural networks (ANNs) were used to estimate the deviations between the simulated and measured wave parameters. The results showed that updating of output variables leads to significant improvement in a wide range of the predicted wave characteristics. It was revealed that the best input parameters for error prediction networks are mean wind speed, mean wind direction, wind duration, and the wave parameters. In addition, combination of the ANN estimated error with numerically modeled wave parameters leads to further improvement in the predicted wave parameters in contrast to direct estimation of the parameters by ANN.     

Agrobiodiversity of Muda Rice Agroecosystem:A case Study in Largest Granary Area of Malaysia

A survey was carried out at the largest rice cultivation area in Peninsular Malaysia,the Muda rice agroecosystem.The main objective of this study was to document the overall biodiversity associated with this unique agroecosystem by using a combination of sampling techniques in order to record different groups of fauna and flora.The total number of biota recorded and identified from the rice field ecosystem during the study period consisted of 46 species of zooplankton,81 species of aquatic insects,5 species of rodents,7 species of bats,87 species of birds,11 species of fishes and 58 species of weeds.A long-term study should be carried out as more species are expected to be recorded when more of the Muda rice agroecosystem area has been sampled to obtain sufficient information on the Muda rice agrobiodiversity.     

Forecasting flood-prone areas using Shannon’s entropy model

With regard to the lack of quality information and data in watersheds, it is of high importance to present a new method for evaluating flood potential. Shannon’s entropy model is a new model in evaluating dangers and it has not yet been used to evaluate flood potential. Therefore, being a new model in determining flood potential, it requires evaluation and investigation in different regions and this study is going to deal with this issue. For to this purpose, 70 flooding areas were recognized and their distribution map was provided by ArcGIS10.2 software in the study area. Information layers of altitude, slope angle, slope aspect, plan curvature, drainage density, distance from the river, topographic wetness index (TWI), lithology, soil type, and land use were recognized as factors affecting flooding and the mentioned maps were provided and digitized by GIS environment. Then, flood susceptibility forecasting map was provided and model accuracy evaluation was conducted using ROC curve and 30% flooding areas express good precision of the model (73.5%) for the study area.     

Stability Analysis and Design of Cantilever Retaining Walls with Regard to Possible Failure Mechanisms: An Upper Bound Limit Analysis Approach

This study presents a new algorithm for design of cantilever retaining walls based on the proposed failure mechanisms and considers the effects of wall geometric parameters using an upper-bound limit analysis approach. All previous work on this subject has only focused on the optimum design of the retaining walls assuming constant forces, irrespective of the total stability and critical conditions of failure mechanisms. In the present study, the upper-bound limit analysis method was used to determine the shape of the critical failure mechanisms for a retaining wall simultaneously with its optimal dimensions. The safety factors against overturning, sliding, and bearing capacity failure were assessed by the limit analysis approach. The current results show good agreement with the results obtained using the limit equilibrium methods and finite element analyses. The results obtained based on the proposed failure mechanism show that the geometry and dimensions of the wall affect its stability safety factors, the shape of the critical failure mechanisms and the active pressure on the wall; therefore, the process of determining the shape of the critical failure mechanisms, checking the stability of the wall and the procedure of finding its optimal dimensions should be performed simultaneously.     

Estimation of enterococci input from bathers and animals on a recreational beach using camera images

Enterococci, are used nationwide as a water quality indicator of marine recreational beaches. Prior research has demonstrated that enterococci inputs to the study beach site (located in Miami, FL) are dominated by non-point sources (including humans and animals). We have estimated their respective source functions by developing a counting methodology for individuals to better understand their non-point source load impacts. The method utilizes camera images of the beach taken at regular time intervals to determine the number of people and animal visitors. The developed method translates raw image counts for week days and weekend days into daily and monthly visitation rates. Enterococci source functions were computed from the observed number of unique individuals for average days of each month of the year, and from average load contributions for humans and for animals. Results indicate that dogs represent the larger source of enterococci relative to humans and birds.     

Probabilistic data integration for characterization of spatial distribution of residual LNAPL

The spatial distribution of residual light non-aqueous phase liquid (LNAPL) is an important factor in reactive solute transport modeling studies. There is great uncertainty associated with both the areal limits of LNAPL source zones and smaller scale variability within the areal limits. A statistical approach is proposed to construct a probabilistic model for the spatial distribution of residual NAPL and it is applied to a site characterized by ultra-violet-induced-cone-penetration testing (CPT–UVIF). The uncertainty in areal limits is explicitly addressed by a novel distance function (DF) approach. In modeling the small-scale variability within the areal limits, the CPT–UVIF data are used as primary source of information, while soil texture and distance to water table are treated as secondary data. Two widely used geostatistical techniques are applied for the data integration, namely sequential indicator simulation with locally varying means (SIS–LVM) and Bayesian updating (BU). A close match between the calibrated uncertainty band (UB) and the target probabilities shows the performance of the proposed DF technique in characterization of uncertainty in the areal limits. A cross-validation study also shows that the integration of the secondary data sources substantially improves the prediction of contaminated and uncontaminated locations and that the SIS–LVM algorithm gives a more accurate prediction of residual NAPL contamination. The proposed DF approach is useful in modeling the areal limits of the non-stationary continuous or categorical random variables, and in providing a prior probability map for source zone sizes to be used in Monte Carlo simulations of contaminant transport or Monte Carlo type inverse modeling studies.     

Mixed estimation methods for Halphen distributions with applications in extreme hydrologic events

The Halphen family of distributions is a flexible and complete system to fit sets of observations independent and identically distributed. Recently, it is shown that this family of distributions represents a potential alternative to the generalized extreme value distributions to model extreme hydrological events. The existence of jointly sufficient statistics for parameter estimation leads to optimality of the method of maximum likelihood (ML). Nevertheless, the ML method requires numerical approximations leading to less accurate values. However, estimators by the method of moments (MM) are explicit and their computation is fast. Even though MM method leads to good results, it is not optimal. In order to combine the advantages of the ML (optimality) and MM (efficiency and fast computations), two new mixed methods were proposed in this paper. One of the two methods is direct and the other is iterative, denoted respectively direct mixed method (MMD) and iterative mixed method (MMI). An overall comparison of the four estimation methods (MM, ML, MMD and MMI) was performed using Monte Carlo simulations regarding the three Halphen distributions. Generally, the MMI method can be considered for the three Halphen distributions since it is recommended for a majority of cases encountered in hydrology. The principal idea of the mixed methods MMD and MMI could be generalized for other distributions with complicated density functions.