Cracking Processes in Rock-Like Material Containing a Single Flaw Under Uniaxial Compression: A Numerical Study Based on Parallel Bonded-Particle Model Approach

Cracking processes have been extensively studied in brittle rock and rock-like materials. Due to the experimental limitations and the complexity of rock texture, details of the cracking processes could not always be observed and assessed comprehensively. To contribute to this field of research, a numerical approach based on the particle element model was used in present study. It would give us insights into what is happening to crack initiation, propagation and coalescence. Parallel bond model, a type of bonded-particle model, was used to numerically simulate the cracking process in rock-like material containing a single flaw under uniaxial vertical compression. The single flaw’s inclinations varied from 0° to 75° measured from the horizontal. As the uniaxial compression load was increased, multiple new microcracks initiated in the rock, which later propagated and eventually coalesced into longer macrocracks. The inclination of the pre-existing flaw was found to have a strong influence on the crack initiation and propagation patterns. The simulations replicated most of the phenomena observed in the physical experiments, such as the type, the initiation location and the initiate angle of the first cracks, as well as the development of hair-line cracks, which later evolved to macrocracks. Analyses of the parallel bond forces and displacement fields revealed some important mechanisms of the cracking processes. The first cracks typically initiated from the tensile stress concentration regions, in which the tensile stress was partially released after their initiation. The tensile stress concentration regions subsequently shifted outwards close to the propagating tips of the first cracks. The initiation and propagation of the first cracks would not significantly influence the compressive stress singularity at the flaw tips, which was the driving force of the initiation of secondary cracks. The initiation of microcracking zone consisting almost exclusively of micro-tensile cracks, and that of microcracking zone consisting of micro-tensile cracks and mixed micro-tensile and shear cracks, were found to be correlated with two distinct types of displacement fields, namely type I (DF_I) and type II (DF_II), respectively.     

Major human Hepatitis A virus genotype in Hong Kong marine waters and detection by real-time PCR

Marine waters from seven sites around Hong Kong with varying levels of sewage pollution were analyzed for Hepatitis A virus (HAV) by PCR cloning and DNA sequencing of the highly variable VP1/2A junction of the HAV genome. Phylogenetic analysis of 10 PCR clones from each of the HAV-positive marine sites indicated that human HAV genotype IB is the most widely distributed type in Hong Kong waters. A sensitive and quantitative TaqMan-based PCR method targeting the 5′-noncoding region (5′-NCR) of HAV was used to quantify HAV particles in marine water samples along with the total Escherichiacoli counts being enumerated on TBX medium for comparison. Our results showed that no correlation of any significance between HAV and E. coli counts was observed which underscores the inadequacy in using E. coli as a sanitary standard to predict the levels of HAV in marine waters.     

Training Programme on Major Tectonics and Lithounits in the Indus and Shyok Suture Zones of Ladakh Himalaya — M. Ram Mohan, (NGRI-CSIR,Hyderabad. Email: rammohan@ngri.res.in) and Sita Bora (Kumaun University,Nainital; Email: sitabora@yahoo.com)

News and Notes

Training Programme on Major Tectonics and Lithounits in the Indus and Shyok Suture Zones of Ladakh Himalaya — M. Ram Mohan, (NGRI-CSIR, Hyderabad. Email: rammohan@ngri.res.in) and Sita Bora (Kumaun University, Nainital; Email: sitabora@yahoo.com)     

Age and geochemistry of the Newania dolomite carbonatites, India: implications for the source of primary carbonatite magma

The Newania carbonatite complex of India is one of the few dolomite-dominated carbonatites of the world. Intruding into Archean basement gneisses, the rocks of the complex have undergone limited diversification and are not associated with any alkaline silicate rock. Although the magmatic nature of the complex was generally accepted, its age of emplacement had remained equivocal because of the disturbed nature of radioisotope systems. Many questions about the nature of its mantle source and mode of origin had remained unanswered because of lack of geochemical and isotopic data. Here, we present results of our effort to date the complex using ¹⁴⁷Sm–¹⁴³Nd, ²⁰⁷Pb–²⁰⁶Pb and ⁴⁰Ar–³⁹Ar dating techniques. We also present mineral chemistry, major and trace element geochemistry and Sr–Nd isotopic ratio data for these carbonatites. Our age data reveal that the complex was emplaced at ~1,473 Ma and parts of it were affected by a thermal event at ~904 Ma. The older ²⁰⁷Pb–²⁰⁶Pb ages reported here (~2.4 Ga) and by one earlier study (~2.3 Ga; Schleicher et al. Chem Geol 140:261–273, 1997) are deemed to be a result of heterogeneous incorporation of crustal Pb during the post-emplacement thermal event. The thermal event had little effect on many magmatic signatures of these rocks, such as its dolomite–magnesite–ankerite–Cr-rich magnetite–magnesio-arfvedsonite–pyrochlore assemblage, mantle like δ¹³C and δ¹⁸O and typical carbonatitic trace element patterns. Newania carbonatites show fractional crystallization trend from high-Mg to high-Fe through high-Ca compositions. The least fractionated dolomite carbonatites of the complex possess very high Mg# (≥80) and have similar major element oxide contents as that of primary carbonatite melts experimentally produced from peridotitic sources. In addition, lower rare earth element (and higher Sr) contents than a typical calcio-carbonatite and mantle like Nb/Ta ratios indicate that the primary magma for the complex was a magnesio-carbonatite melt and that it was derived from a carbonate bearing mantle. The Sr–Nd isotopic data suggest that the primary magma originated from a metasomatized lithospheric mantle. Trace element modelling confirms such an inference and suggests that the source was a phlogopite bearing mantle, located within the garnet stability zone.     

Dynamic Study on Fracture Problems in Viscoelastic Sedimentary Rocks Using the Numerical Manifold Method

The viscoelastic deformation behavior of a sedimentary rock under different loading rates is numerically modeled and investigated by the numerical manifold method (NMM). By incorporating a modified 3-element viscoelastic constitutive mode in the NMM, crack initiation and propagation criteria, and crack identification and evolution techniques, the effects of the loading rates on the cracking behavior of a sedimentary rock, such as crack open displacement, crack sliding displacement, crack initiation, crack propagation and final failure mode, are successfully modeled. The numerical results reveal that under a high loading rate (>1,000 MPa/s), due to the viscoelastic property of the sedimentary rock, not only the structural behavior deviates from that of elastic model, but also different cracking processes and final failure modes are obtained.     

Thermal metamorphism of the Arunachal Himalaya, India: Raman thermometry and thermochronological constraints on the tectono-thermal evolution

Determination of the peak thermal condition is vital in order to understand tectono-thermal evolution of the Himalayan belt. The Lesser Himalayan Sequence (LHS) in the Western Arunachal Pradesh, being rich in carbonaceous material (CM), facilitates the determination of peak metamorphic temperature based on Raman spectroscopy of carbonaceous material (RSCM). In this study, we have used RSCM method of Beyssac et al. (J Metamorph Geol 20:859–871, 2002a) and Rahl et al. (Earth Planet Sci Lett 240:339–354, 2005) to estimate the thermal history of LHS and Siwalik foreland from the western Arunachal Pradesh. The study indicates that the temperature of 700–800 °C in the Greater Himalayan Sequence (GHS) decreases to 650–700 °C in the main central thrust zone (MCTZ) and decreases further to <200 °C in the Mio-Pliocene sequence of Siwaliks. The work demonstrates greater reliability of Rahl et al.’s (Earth Planet Sci Lett 240:339–354, 2005) RSCM method for temperatures >600 and <340 °C. We show that the higher and lower zones of Bomdila Gneiss (BG) experienced temperature of ~600 °C and exhumed at different stages along the Bomdila Thrust (BT) and Upper Main Boundary Thrust (U.MBT). Pyrolysis analysis of the CM together with the Fission Track ages from upper Siwaliks corroborates the RSCM thermometry estimate of ~240 °C. The results indicate that the Permian sequence north of Lower MBT was deposited at greater depths (>12 km) than the upper Siwalik sediments to its south at depths <8 km before they were exhumed. The ⁴⁰Ar/³⁹Ar ages suggest that the upper zones of Se La evolved ~13–15 Ma. The middle zone exhumed at ~11 Ma and lower zone close to ~8 Ma indicating erosional unroofing of the MCT sheet. The footwall of MCTZ cooled between 6 and 8 Ma. Analyses of P–T path imply that LHS between MCT and U.MBT zone falls within the kyanite stability field with near isobaric condition. At higher structural level, the temperatures increase gradually with P–T conditions in the sillimanite stability field. The near isothermal (700–800 °C) condition in the GHS, isobaric condition in the MCTZ together with T–t path evidence of GHS that experienced relatively longer duration of near peak temperatures and rapid cooling towards MCTZ, compares the evolution of GHS and inverted metamorphic gradient closely to channel flow predictions.     

An Efficient System for Creating Synthetic InSAR Images from Simulations

In this work we visualize tsunami and earthquake simulation results with graphics hardware acceleration. The rapid improvement in the computational power of graphics hardware and its programmability has made general computation on Graphics Processing Units (GPUs) very compelling. We generate Synthetic InSAR images using GPUs. Interference phenomena have formed the underlying theory for Interferometric Synthetic Aperture Radar (InSAR) in unveiling dynamical Earth movements. In our approach light path differences are defined by the surface values to be visualized. These path differences then modulate the lighting intensity to generate the interference patterns. We can interactively visualize surface deformation patterns by leveraging the computational power of GPUs. Our visualization method is applied to simulations of rupture fault displacements during the tsunamogenic earthquake events, which are vital to understanding the subsequent wave propagation. We also integrate the visualization results into Google Earth virtual globe to provide the geological context of the visualized regions.     

Differential developmental toxicity of naphthoic acid isomers in medaka (Oryzias latipes) embryos

Polycyclic aromatic hydrocarbons (PAHs) are widespread persistent pollutants that readily undergo biotic and abiotic conversion to numerous transformation products in rivers, lakes and estuarine sediments. Here we characterize the developmental toxicity of four PAH transformation products each structural isomers of hydroxynaphthoic acid: 1H2NA, 2H1NA, 2H3NA, and 6H2NA. Medaka fish (Oryzias latipes) embryos and eleutheroembryos were used to determine toxicity. A 96-well micro-plate format was used to establish a robust, statistically significant platform for assessment of early life stages. Individual naphthoic acid isomers demonstrated a rank order of toxicity with 1H2NA > 2H1NA > 2H3NA > 6H2NA being more toxic. Abnormalities of circulatory system were most pronounced including pericardial edema and tube heart. To determine if HNA isomers were AhR ligands, spatial-temporal expression and activity of CYP1A was measured via in vivo EROD assessments. qPCR measurement of CYP1A induction proved different between isomers dosed at respective concentrations affecting 50% of exposed individuals (EC50s). In vitro, all ANH isomers transactivated mouse AhR using a medaka CYP1A promoter specific reporter assay. Circulatory abnormalities followed P450 induction and response was consistent with PAH toxicity. A 96-well micro-plates proved suitable as exposure chambers and provided statistically sound evaluations as well as efficient toxicity screens. Our results demonstrate the use of medaka embryos for toxicity analysis thereby achieving REACH objectives for the reduction of adult animal testing in toxicity evaluations.     

The role of magnetic buoyancy in a Babcock-Leighton type solar dynamo

We study the effects of incorporating magnetic buoyancy in a model of the solar dynamo&amp;#x2014;which draws inspiration from the Babcock-Leighton idea of surface processes generating the poloidal field. We present our main results here.