Fracture mechanics analysis of multiple cracks in anisotropic media

This paper presents a single‐domain boundary element method (BEM) for linear elastic fracture mechanics analysis in the two‐dimensional anisotropic material. In this formulation, the displacement integral equation is collocated on the un‐cracked boundary only, and the traction integral equation is collocated on one side of the crack surface only. A special crack‐tip element was introduced to capture exactly the crack‐tip behavior. A computer program with the FORTRAN language has been developed to effectively calculate the stress intensity factors of an anisotropic material. This BEM program has been verified having a good accuracy with the previous researches. Furthermore, by analyzing the different anisotropic degree cracks in a finite plate, we found that the stress intensity factors of crack tips had apparent influence by the geometry forms of cracks and media with different anisotropic degrees. Copyright © 2010 John Wiley & Sons, Ltd.     

Chemical weathering and CO2 consumption in the glaciated Karuxung River catchment,Tibetan Plateau

Climate factors play critical roles in controlling chemical weathering, while chemically weathered surface material can regulate climate change. To estimate global chemical weathering fluxes and CO₂ balance, it is important to identify the characteristics and driving factors of chemical weathering and CO₂ consumption on the Tibetan Plateau, especially in glaciated catchments. The analysis of the hydro-geochemical data indicated that silicate weathering in this area was inhibited by low temperatures, while carbonate weathering was promoted by the abundant clastic rocks with fresh surfaces produced by glacial action. Carbonate weathering dominated the riverine solute generation (with a contribution of 58%, 51%, and 43% at the QiangYong Glacier (QYG), the WengGuo Hydrological Station (WGHS), and the lake estuary (LE), respectively). The oxidation of pyrite contributed to 35%, 42%, and 30% of the riverine solutes, while silicate weathering contributed to 5%, 6%, and 26% of the riverine solutes at the QYG, WGHS, and LE, respectively. The alluvial deposit of easily weathering fine silicate minerals, the higher air temperature, plant density, and soil thickness at the downstream LE in comparison to upstream and midstream may lead to longer contact time between pore water and mineral materials, thus enhancing the silicate weathering. Because of the involvement of sulfuric acid produced by the oxidation of pyrite, carbonate weathering in the upstream and midstream did not consume atmospheric CO₂, resulting in the high rate of carbonate weathering (73.9 and 75.6 t km⁻² yr⁻¹, respectively, in maximum) and potential net release of CO₂ (with an upper constraint of 35.6 and 35.2 t km⁻² yr⁻¹, respectively) at the QYG and WGHS. The above results indicate the potential of the glaciated area of the Tibetan Plateau with pyrite deposits being a substantial natural carbon source, which deserves further investigation.     

Petrological evidence for shock‐induced high‐P metamorphism in a gabbro

Microlites (minute spherulitic, dendritic, skeletal, acicular and poikilitic crystals) diagnostic of crystallization in quenched melt or glass in fault rocks have been used to infer fossil earthquakes. High‐P microlites and crystallites are described here in a variably eclogitized gabbro, the wallrock to the coesite‐bearing eclogite breccia at Yangkou in the Chinese Su‐Lu high‐P metamorphic belt. The studied hand specimens are free of discernible shear deformation, although microfractures are not uncommon under the microscope. In the least eclogitized gabbro, the metagabbro, stellate growths of high‐P minerals on the relict igneous minerals are common. Dendritic garnet crystals (<1?5 μm) grew around rutile and/or phengite replacing ilmenite and biotite, respectively. Skeletal garnet also rims broken flakes of igneous biotite and mechanically twinned augite. Radial intergrowths of omphacite and quartz developed around relict igneous orthopyroxene and are rimmed by skeletal or poikilitic garnet where a Ti‐bearing mineral relict is present. Acicular epidote, kyanite and phengite crystallites are randomly distributed in a matrix of Na‐rich plagioclase, forming the pseudomorphs after igneous plagioclase. In the more eclogitized gabbro, the coronitic eclogite located closer to the eclogite breccia, all the igneous minerals broke down into high‐P assemblages. Thick coronas of poikilitic garnet grew between the pseudomorphs after igneous plagioclase and ferromagnesian minerals. The igneous plagioclase is replaced by omphacite crystallites, with minor amounts of phengite and kyanite. Thermodynamic modelling of the plagioclase pseudomorphs shows an increase in P–T in the wallrock from the metagabbro to the coronitic eclogite, and the P–T variation is unrelated to H₂O content. The fluid‐poor pressure overstepping scenario is unsupported both by phase diagram modelling and by whole‐rock chemical data, which show that the various types of eclogitized gabbro are all fairly dry. A large pressure difference of >2 GPa between the metagabbro and the coesite‐bearing eclogites ~20 m apart cannot be explained by the subduction hypothesis because this would require a depth difference of >60 km. The microlites and crystallites are evidence for dynamic crystallization due to rapid cooling because constitutional supercooling was unlikely for the plagioclase pseudomorphs. The lack of annealing of the broken biotite and augite overgrown by strain free skeletal garnet is consistent with a transient high‐P–T event at a low ambient temperature (<300 °C), probably in the crust. Therefore, the eclogitization of the wallrock to the eclogite breccia was also coseismic, as proposed earlier for the eclogite facies fault rocks. The outcrop‐scale P–T variation and the transient nature of the high‐P–T event are inconsistent with the other existing tectonic models for high‐P metamorphism. The fact that the less refractory but denser biotite is largely preserved while the more refractory but less dense plagioclase broke down completely into high‐P microlite assemblages in the metagabbro indicates a significant rise in pressure rather than temperature. Given that the metamorphic temperatures are far below the melting temperatures of most of the gabbroic minerals under fluid‐absent conditions, stress‐induced amorphization appears to be the more likely mechanism of the coseismic high‐P metamorphism.     

Using raw regional climate model outputs for quantifying climate change impacts on hydrology

General circulation model outputs are rarely used directly for quantifying climate change impacts on hydrology, due to their coarse resolution and inherent bias. Bias correction methods are usually applied to correct the statistical deviations of climate model outputs from the observed data. However, the use of bias correction methods for impact studies is often disputable, due to the lack of physical basis and the bias nonstationarity of climate model outputs. With the improvement in model resolution and reliability, it is now possible to investigate the direct use of regional climate model (RCM) outputs for impact studies. This study proposes an approach to use RCM simulations directly for quantifying the hydrological impacts of climate change over North America. With this method, a hydrological model (HSAMI) is specifically calibrated using the RCM simulations at the recent past period. The change in hydrological regimes for a future period (2041–2065) over the reference (1971–1995), simulated using bias‐corrected and nonbias‐corrected simulations, is compared using mean flow, spring high flow, and summer–autumn low flow as indicators. Three RCMs driven by three different general circulation models are used to investigate the uncertainty of hydrological simulations associated with the choice of a bias‐corrected or nonbias‐corrected RCM simulation. The results indicate that the uncertainty envelope is generally watershed and indicator dependent. It is difficult to draw a firm conclusion about whether one method is better than the other. In other words, the bias correction method could bring further uncertainty to future hydrological simulations, in addition to uncertainty related to the choice of a bias correction method. This implies that the nonbias‐corrected results should be provided to end users along with the bias‐corrected ones, along with a detailed explanation of the bias correction procedure. This information would be especially helpful to assist end users in making the most informed decisions.     

A tree‐ring width based drought reconstruction for southeastern China: links to Pacific Ocean climate variability

We present a drought reconstruction for southeastern China based on a tree‐ring width chronology of Cryptomeria fortunei developed from two sampling sites in central Fujian. A reconstruction of July–February drought variability, spanning AD 1855–2011, was developed by calibrating total tree‐ring width data with the self‐calibrating Palmer drought severity index (scPDSI). The reconstruction was verified against an independent data set, and accounts for 36% of the actual scPDSI variance during the period 1955–2011. Relatively dry intervals were reconstructed between AD 1859–1880, 1899–1911, 1927–1933, 1946–1959, 1964–1970 and 1987–1997. Relatively wet conditions prevailed during 1855–1858, 1881–1898, 1912–1926, 1934–1945, 1960–1963, 1971–1986 and 1998–2011. Comparisons between our scPDSI reconstruction and a moisture‐sensitive tree‐ring width record from Vietnam revealed consistencies between the two data sets, suggesting similar drought regimes. Spectral peaks of 2.2–6.4 years may be indicative of El Niño‐Southern Oscillation (ENSO) activity, as also suggested by the significant correlations with sea surface temperatures (SSTs) in the eastern equatorial and southeastern Pacific Ocean and an extreme event analysis. The analysis of links between our scPDSI reconstruction and the large‐scale regional climatic variation shows that there is a relationship between regional drought variation and East Asian summer monsoon (EASM) intensity.     

Stable axisymmetric SPH formulation with no axis singularity

The axisymmetric formulation of the governing equations for geomechanics in the framework of smoothed particle hydrodynamics (SPH) is presented in this study. Two forms of SPH discretization for the motion equations, which are labeled as form I and form II, are proposed, and the methods to compute the hoop stress and strain terms including hoop strain rate and the acceleration introduced by the hoop stress are compared. To avoid possible singularity problem near the axis of symmetry, a perfectly smooth contact along with ghost particles are applied to prevent the real particles from overly approaching the axis of symmetry to remove this potential singularity. In addition, the Mohr–Coulomb constitutive model is implemented into the SPH formulation in describing soil behavior. Four numerical tests are carried out to validate and compare the accuracy and stability of the proposed algorithms, and their results are compared with analytical solutions and results from FEM analysis. The performance in these comparisons suggests that SPH II with hoop terms computed through direct hoop method is more stable than the others, and the adoption of contact for the symmetric axis is efficient in eliminating the singularity problem. Copyright © 2015 John Wiley & Sons, Ltd.     

The influence of coupled physical swelling and chemical reactions on deformable geomaterials

Coupled thermo‐hydro‐mechanical‐chemical modelling has attracted attention in past decades due to many contemporary geotechnical engineering applications (e.g., waste disposal, carbon capture and storage). However, molecular‐scale interactions within geomaterials (e.g., swelling and dissolution/precipitation) have a significant influence on the mechanical behaviour, yet are rarely incorporated into existing Thermal‐Hydro‐Mechanical‐Chemical (THMC) frameworks. This paper presents a new coupled hydro‐mechanical‐chemical constitutive model to bridge molecular‐scale interactions with macro‐physical deformation by combining the swelling and dissolution/precipitation through an extension of the new mixture‐coupling theory. Entropy analysis of the geomaterial system provides dissipation energy, and Helmholtz free energy gives the relationship between solids and fluids. Numerical simulation is used to compare with the selected recognized models, which demonstrates that the swelling and dissolution/precipitation processes may have a significant influence on the mechanical deformation of the geomaterials.     

Effects of temporal scales and space mismatches on the TRMM 3B42 v7 precipitation product in a remote mountainous area

The accurate measurement of precipitation is essential to understanding regional hydrological processes and hydrological cycling. Quantification of precipitation over remote regions such as the Tibetan Plateau is highly unreliable because of the scarcity of rain gauges. The objective of this study is to evaluate the performance of the satellite precipitation product of tropical rainfall measuring mission (TRMM) 3B42 v7 at daily, weekly, monthly, and seasonal scales. Comparison between TRMM grid precipitation and point‐based rain gauge precipitation was conducted using nearest neighbour and bilinear weighted interpolation methods. The results showed that the TRMM product could not capture daily precipitation well due to some rainfall events being missed at short time scales but provided reasonably good precipitation data at weekly, monthly, and seasonal scales. TRMM tended to underestimate the precipitation of small rainfall events (less than 1 mm/day), while it overestimated the precipitation of large rainfall events (greater than 20 mm/day). Consequently, TRMM showed better performance in the summer monsoon season than in the winter season. Through comparison, it was also found that the bilinear weighted interpolation method performs better than the nearest neighbour method in TRMM precipitation extraction.