Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Wetland hydroperiod classification in the western prairies using multitemporal synthetic aperture radar

Wetlands represent one of the world's most biodiverse and threatened ecosystem types and were diminished globally by about two‐thirds in the 20th century. There is continuing decline in wetland quantity and function due to infilling and other human activities. In addition, with climate change, warmer temperatures and changes in precipitation and evapotranspiration are reducing wetland surface and groundwater supplies, further altering wetland hydrology and vegetation. There is a need to automate inventory and monitoring of wetlands, and as a study system, we investigated the Shepard Slough wetlands complex, which includes numerous wetlands in urban, suburban, and agricultural zones in the prairie pothole region of southern Alberta, Canada. Here, wetlands are generally confined to depressions in the undulating terrain, challenging wetlands inventory and monitoring. This study applied threshold and frequency analysis routines for high‐resolution, single‐polarization (HH) RADARSAT‐2, synthetic aperture radar mapping. This enabled a growing season surface water extent hyroperiod‐based wetland classification, which can support water and wetland resource monitoring. This 3‐year study demonstrated synthetic aperture radar‐derived multitemporal open‐water masks provided an effective index of wetland permanence class, with overall accuracies of 89% to 95% compared with optical validation data, and RMSE between 0.2 and 0.7 m between model and field validation data. This allowed for characterizing the distribution and dynamics of 4 marsh wetlands hydroperiod classes, temporary, seasonal, semipermanent, and permanent, and mapping of the sequential vegetation bands that included emergent, obligate wetland, facultative wetland, and upland plant communities. Hydroperiod variation and surface water extent were found to be influenced by short‐term rainfall events in both wet and dry years. Seasonal hydroperiods in wetlands were particularly variable if there was a decrease in the temporary or semipermanent hydroperiod classes. In years with extreme rain events, the temporary wetlands especially increased relative to longer lasting wetlands (84% in 2015 with significant rainfall events, compared with 42% otherwise).     

Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

Assessment of alternative adsorption models and global sensitivity analysis to characterize hexavalent chromium loss from soil to surface runoff

We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AIC_C] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component.     

Metamorphic phase relations in orthopyroxene-bearing granitoids: implication for high-pressure metamorphism and prograde melting in the continental crust

In this work, the factors controlling the formation and preservation of high-pressure mineral assemblages in the metamorphosed orthopyroxene-bearing metagranitoids of the Sandmata Complex, Aravalli-Delhi Mobile Belt (ADMB), northwestern India have been modelled. The rocks range in composition from farsundite through quartz mangerite to opdalite, and with varying K₂O, Ca/(Ca + Na)_rock and FeO_tot + MgO contents. A two stage metamorphic evolution has been recorded in these rocks.
An early hydration event stabilized biotite with or without epidote at the expense of magmatic orthopyroxene and plagioclase. Subsequent high-pressure granulite facies metamorphism (∼15 kbar, ∼800 °C) of these hydrated rocks produced two rock types with contrasting mineralogy and textures. In the non-migmatitic metagranitoids, spectacular garnet ± K-feldspar ± quartz corona was formed around reacting biotite, plagioclase, quartz and/or pyroxene. In contrast, biotite ± epidote melting produced migmatites, containing porphyroblastic garnet incongruent solids and leucosomes.
Applying NCKFMASHTO T–M (H₂O) and P–T pseudosection modelling techniques, it is demonstrated that the differential response of these magmatic rocks to high-pressure metamorphism is primarily controlled by the scale of initial hydration. Rocks, which were pervasively hydrated, produced garnetiferous migmatites, while for limited hydration, the same metamorphism formed sub-solidus garnet-bearing coronae. Based on the sequence of mineral assemblage evolution and the mineral compositional zoning features in the two metagranitoids, a clockwise metamorphic P–T path is constrained for the high-pressure metamorphic event. The finding has major implications in formulating geodynamic model of crustal amalgamation in the ADMB.     

An atlas of the optical spectrum of supernova 1987A in the large magellanic cloud

Photographic spectra of SN1987A in the LMC have been obtained from 1987 February 25 to 1988 June 30. Microdensitometer tracings of these have been reduced to intensity and corrections for instrumental response have been applied to the spectra. This paper presents these data in an atlas format, discusses the reduction procedures in detail, and presents radial velocity measurements of selected lines in the spectra     

利用自回归模型检测地球极移重力效应

同时采用4个台站的国际超导重力仪长期连续重力观测资料和国际地球自转服务中心提供的同步地球自转参数，研究了极移引起的地球重力场变化特征。利用自回归模型估计了各序列的功率谱密度和积谱密度，结果表明极移导致的重力效应的主要能量集中在Chandler摆动和周年项附近，叠积后实际重力观测与极移重力信号理论值之间的差异分别为0．4％和3．9％，说明超导重力仪可有效监测极移导致的重力变化。     

The evolution of binary fractions in globular clusters

We study the evolution of binary stars in globular clusters using a new Monte Carlo approach combining a population synthesis code ( startrack ) and a simple treatment of dynamical interactions in the dense cluster core using a new tool for computing three- and four-body interactions ( fewbody ). We find that the combination of stellar evolution and dynamical interactions (binary–single and binary–binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster such as 47 Tuc, assuming an initial binary fraction of 100 per cent, is only ∼ 5–10 per cent. We show that this is in good agreement with recent Hubble Space Telescope observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modelling of globular clusters, suggesting that 'realistic models' should incorporate much larger initial binary fractions than has usually been the case in the past.     

Atmospheric radiative transfer simulation for atmospheric correction of remote sensing data

IntroductionThe radiance leaving the earth-atmosphere sys-tem which can be sensed by a satellite borne ra-diometer is the sum of radiation emission fromtheearth surface and each atmospheric level that aretransmittedtothe top of the atmosphere.The radia-tion emissionfromthe earthsurface andthe radianceof each atmospheric level can be separated fromtheradiance at the top the atmospheric level a satellitemeasured.Thus,the earth surface parameters willbe retrieved from the surface radiance after a…     

Reliability assessment of granular filters in embankment dams

Empirical criteria have been used successfully to design filters of most embankment large dam projects throughout the world. However, these empirical rules are only applicable to a particular range of soils tested in laboratory and do not take into account the variability of the base material and filter particle sizes. In addition, it is widely accepted that the safety of fill dams is mainly dependent on the reliability of their filter performance. The work herein presented consists in a new general method for assessing the probability of fulfilling any empirical filter design criteria accounting for base and filter heterogeneity by means of first‐order reliability methods (FORM), so that reliability indexes and probabilities of fulfilling any particular criteria are obtained. This method will allow engineers to estimate the safety of existing filters in terms of probability of fulfilling their design criteria and might also be used as a decision tool on sampling needs and material size tolerances during construction. In addition, sensitivity analysis makes possible to analyse how reliabilities are influenced by different sources of input data. Finally, in case of a portfolio risk assessment, this method will allow engineers to compare the safety of several existing dams in order to prioritize safety investments and it is expected to be a very useful tool to evaluate probabilities of failure due to internal erosion. Copyright © 2006 John Wiley & Sons, Ltd.