Investigation into the effect of fault properties on wave transmission

The estimation of wave transmission across the fractured rock masses is of great importance for rock engineers to assess the stability of rock slopes in open pit mines. Presence of fault, as a major discontinuity, in the jointed rock mass can significantly impact on the peak particle velocity and transmission of blast waves, particularly where a fault contains a thick infilling with weak mechanical properties. This paper aims to study the effect of fault properties on transmission of blasting waves using the distinct element method. First, a validation study was carried out on the wave transmission across a single joint and different rock mediums through undertaking a comparative study against analytical models. Then, the transmission of blast wave across a fault with thick infilling in the Golgohar iron mine, Iran, was numerically studied, and the results were compared with the field measurements. The blast wave was numerically simulated using a hybrid finite element and finite difference code which then the outcome was used as the input for the distinct element method analysis. The measured uplift of hanging wall, as a result of wave transmission across the fault, in the numerical model agrees well with the recorded field measurement. Finally, the validated numerical model was used to study the effect of fault properties on wave transmission. It was found that the fault inclination angle is the most effective parameter on the peak particle velocity and uplift. Copyright © 2017 John Wiley & Sons, Ltd.     

Constitutive and numerical framework for modeling joints and faults in rock

A three‐dimensional constitutive model for joints is described that incorporates nonlinear elasticity based on volumetric elastic strain, and plasticity for both compaction and shear with emphasis on compaction. The formulation is general in the sense that alternative specific functional forms and evolution equations can be easily incorporated. A corresponding numerical structure based on finite elements is provided so that a joint width can vary from a fraction of an element size to a width that occupies several elements. The latter case is particularly appropriate for modeling a fault, which is considered simply to be a joint with large width. For small joint widths, the requisite equilibrium and kinematic requirements within an element are satisfied numerically. The result is that if the constitutive equation for either the joint or the rock is changed, the numerical framework remains unchanged. A unique aspect of the general formulation is the capability to handle either pre‐existing gaps or the formation of gaps. Representative stress–strain plots are given to illustrate both the features of the model and the effects of changes in values of material parameters. Copyright © 2015 John Wiley & Sons, Ltd.     

Precessional forcing of tropical vegetation carbon storage

Since the Mid Pleistocene Revolution, which occurred about one million years ago, global temperatures have fluctuated with a quasi‐periodicity of ca. 100 ka. The pattern of past change in the extent of woodlands, and therefore by inference vegetation carbon storage, has been demonstrated to have a strong positive link with this global temperature change at high and mid latitudes. However, understanding of climate systems and ecosystem function indicates that the pattern of woodland change at low latitudes may follow a fundamentally different pattern. We present output from the intermediate complexity model GENIE‐1, comprising a single transient simulation over the last 800 ka and a 174‐member ensemble of 130 ka transient simulations over the last glacial cycle. These simulations suggest that while vegetation carbon storage in mid–high northern latitudes robustly follows the characteristic ca. 100 ka cycle, this signal is not a robust feature of tropical vegetation, which is subject to stronger direct forcing by the precessional (21 ka) orbital cycle (albeit with a highly uncertain response). We conclude that the correlation of palaeoenvironmental records from low latitudes with global temperature change must be done with caution. Copyright © 2011 John Wiley & Sons, Ltd.     

Three-dimensional active earth pressure from cohesive backfills with tensile strength cutoff

Most of previous analyses on the active earth pressure were performed in two-dimensional cases using the Mohr-Coulomb (M-C) failure function to describe the soil strength. However, all failures of retained slopes indicate a somewhat three-dimensional (3D) feature, and the M-C function is found to overestimate the tensile strength of cohesive soil. In this work, a kinematic limit analysis–based approach is developed for computing the 3D active earth pressure resulting from cohesive backfills. The concept of tensile strength cutoff is adopted to implement the reduction or elimination of tensile strength from the strength envelope. An extended 3D horn failure mechanism that is associated with the modified strength envelope is developed to characterize the collapse of retained slopes. The resultant of active earth pressure is evaluated from the work rate balance equation and expressed as an unfactored coefficient. The obtained results indicate that less support provided by the wall is required when allowing the existence of soil cohesion and 3D effects and that eliminating the tensile strength can observably increase the active earth pressure, especially for the backfill with a great level of cohesion.     

Why asymptotic limit of the acid dissolution capacity can lead to a sharp dissolution front in chemical dissolution of porous rocks?

The use of the asymptotic limit can greatly simplify the theoretical analysis of chemical dissolution front instabilities in fluid‐saturated rocks and therefore make it possible to obtain mathematical solutions, which often play a crucial role in understanding the propagation behavior of chemical dissolution fronts in chemical dissolution systems. However, there has been a debate in recent years that the asymptotic limit of the acid dissolution capacity (i.e., the acid dissolution capacity number approaching zero) alone cannot lead to a sharp dissolution front of the Stefan type in the acidization dissolution system, in which the dissolvable minerals of carbonate rocks are chemically dissolved by the injected acid flow. The acid dissolution capacity number is commonly defined as the ratio of the volume of the carbonate rock dissolved by an acid to that of the acid. In this paper, we use four different proof methods, including (i) direct use of the fundamental concepts; (ii) use of the mathematical governing equations of an acidization dissolution system; (iii) use of the different time scaling approach; and (iv) use of a moving coordinate system approach, to demonstrate that the asymptotic limit of the acid dissolution capacity can indeed lead to sharp dissolution fronts of the Stefan type in acidization dissolution systems on a much larger time scale (than the dissolution time scale). Our new finding is that on the reaction time scale, the condition of the conventional time derivative of porosity approaching zero alone can ensure that the acidization dissolution front has a sharp shape of the Stefan type. Copyright © 2017 John Wiley & Sons, Ltd.     

Inverse modelling using PS‐InSAR data for improved land subsidence simulation in Las Vegas Valley,Nevada

Las Vegas Valley has had a long history of groundwater development and subsequent surface deformation. InSAR interferograms have revealed detailed and complex spatial patterns of subsidence in the Las Vegas Valley area that do not coincide with major pumping regions. This research represents the first effort to use high spatial and temporal resolution subsidence observations from InSAR and hydraulic head data to inversely calibrate transmissivities (T), elastic and inelastic skeletal storage coefficients (S_ke and S_kv) of the developed‐zone aquifer and conductance (CR) of the basin‐fill faults for the entire Las Vegas basin. The results indicate that the subsidence observations from InSAR are extremely beneficial for accurately quantifying hydraulic parameters, and the model calibration results are far more accurate than when using only groundwater levels as observations, and just a limited number of subsidence observations. The discrepancy between distributions of pumping and greatest levels of subsidence is found to be attributed to spatial variations in clay thickness. The Eglington fault separates thicker interbeds to the northwest from thinner interbeds to the southeast and the fault may act as a groundwater‐flow barrier and/or subsidence boundary, although the influence of the groundwater barrier to this area is found to be insignificant. Copyright © 2016 John Wiley & Sons, Ltd.     

An analytical study on three‐dimensional versus two‐dimensional water table‐induced flow patterns in a Tóthian basin

Although it has been increasingly acknowledged that groundwater flow pattern is complicated in the three‐dimensional (3‐D) domain, two‐dimensional (2‐D) water table‐induced flow models are still widely used to delineate basin‐scale groundwater circulation. However, the validity of 2‐D cross‐sectional flow field induced by water table has been seldom examined. Here, we derive the analytical solution of 3‐D water table‐induced hydraulic head in a Tóthian basin and then examine the validity of 2‐D cross‐sectional models by comparing the flow fields of selected cross sections calculated by the 2‐D cross‐sectional model with those by the 3‐D model, which represents the “true” cases. For cross sections in the recharge or discharge area of the 3‐D basin, even if head difference is not significant, the 2‐D cross‐sectional models result in flow patterns absolutely different from the true ones. For the cross section following the principal direction of groundwater flow, although 2‐D cross‐sectional models would overestimate the penetrating depth of local flow systems and underestimate the recharge/discharge flux, the flow pattern from the cross‐sectional model is similar to the true one and could be close enough to the true one by adjusting the decay exponent and anisotropy ratio of permeability. Consequently, to determine whether a 2‐D cross‐sectional model is applicable, a comparison of hydraulic head difference between 2‐D and 3‐D solutions is not enough. Instead, the similarity of flow pattern should be considered to determine whether a cross‐sectional model is applicable. This study improves understanding of groundwater flow induced by more natural water table undulations in the 3‐D domain and the limitations of 2‐D models accounting for cross‐sectional water table undulation only.     

Extended finite element framework for fault rupture dynamics including bulk plasticity

We present an explicit extended finite element framework for fault rupture dynamics accommodating bulk plasticity near the fault. The technique is more robust than the standard split‐node method because it can accommodate a fault propagating freely through the interior of finite elements. To fully exploit the explicit algorithmic framework, we perform mass lumping on the enriched finite elements that preserve the kinetic energy of the rigid body and enrichment modes. We show that with this technique, the extended FE solution reproduces the standard split‐node solution, but with the added advantage that it can also accommodate randomly propagating faults. We use different elastoplastic constitutive models appropriate for geomaterials, including the Mohr–Coulomb, Drucker–Prager, modified Cam‐Clay, and a conical plasticity model with a compression cap, to capture off‐fault bulk plasticity. More specifically, the cap model adds robustness to the framework because it can accommodate various modes of deformation, including compaction, dilatation, and shearing. Copyright © 2013 John Wiley & Sons, Ltd.     

Paleogene of the Huanghekou Sag in the Bohai Bay Basin,NE China: deposition–erosion response to a slope break system of rift lacustrine basins

The sequence architecture and depositional systems of the Paleogene lacustrine rift succession in the Huanghekou Sag, Bohai Bay Basin, NE China were investigated based on seismic profiles, combined with well log and core data. Four second‐order or composite sequences and seven third‐order sequences were identified. The depositional systems identified in the basin include: fan delta, braid delta, meander fluvial delta, lacustrine and sublacustrine fan. Identification of the slope break was conducted combining the interpretation of faults of each sequence and the identification of syndepositional faults, based on the subdivision of sequence stratigraphy and analysis of depositional systems. Multiple geomorphologic units were recognized in the Paleogene of the Huanghekou Sag including faults, flexures, depositional slope break belts, ditch‐valleys and sub‐uplifts in the central sag. Using genetic division principles and taking into consideration tectonic features of the Paleogene of the Huanghekou Sag, the study area was divided into the Northern Steep Slope/Fault Slope Break System, the Southern Gentle Slope Break System and T10 Tectonic Slope Break System/T10 Tectonic Belt. Responses of slope break systems to deposition–erosion are shown as: (1) basin marginal slope break is the boundary of the eroded area and provenance area; (2) ditch‐valley formed by different kinds of slope break belts is a good transport bypass for source materials; (3) shape of the slope break belt of the slope break system controls sediments types; (4) the ditch‐valley and sub‐sag of a slope break system is an unloading area for sediments; and (5) due to their different origins, association characteristics and developing patterns, the Paleogene slope break belt systems in the Huanghekou Sag show different controls on depositional systems. The Northern Fault Slope Break system controls the deposition of a fan delta‐lacustrine‐subaqueous fan, the Southern Gentle Slope Break system controls the deposition of a fluvial–deltaic–shallow lacustrine and sublacustrine fan, and the T10 Tectonic Slope Break System controls the deposition of shallow lacustrine beach bar sandbodies. The existence of a slope break system is a necessary but not a sufficient condition for studying sandbody development. The formation of effective sandbodies along the slope break depends on the reasonable coupling of effective provenance, necessary association patterns of slope break belt, adequate unloading space and creation of definite accommodation space. Copyright © 2013 John Wiley & Sons, Ltd.