Geochemical and isotopic constraints on the age and origin of the Nidar Ophiolitic Complex, Ladakh, India: Implications for the Neo-Tethyan subduction along the Indus suture zone

The Nidar ophiolite complex is exposed within the Indus suture zone in eastern Ladakh, India. The suture zone is considered to represent remnant Neo-Tethyan Ocean that closed via subduction as the Indian plate moved northward with respect to the Asian plate. The two plates ultimately collided during the Middle Eocene. The Nidar ophiolite complex comprises a sequence of ultra-mafic rocks at the base, gabbroic rocks in the middle and volcano-sedimentary assemblage on the top. Earlier studies considered the Nidar ophiolite complex to represent an oceanic floor sequence based on lithological assemblage. However, present study, based on new mineral and whole rock geochemical and isotopic data (on bulk rocks and mineral separates) indicate their generation and emplacement in an intra-oceanic subduction environment. The plutonic and volcanic rocks have nearly flat to slightly depleted rare earth element (REE) patterns. The gabbroic rocks, in particular, show strong positive Sr and Eu anomalies in their REE and spidergram patterns, probably indicating plagioclase accumulation. Depletion in high field strength elements (HFSE) in the spidergram patterns may be related to stabilization of phases retaining the HFSE in the subducting slab and / or fractional crystallization of titano-magnetite phases. The high radiogenic Nd- and low radiogenic Sr-isotopic ratios for these rocks exclude any influence of continental material in their genesis, implying an intra-oceanic environment.

Nine point mineral–whole rock Sm–Nd isochron corresponds to an age of 140 ± 32 Ma with an initial ¹⁴³Nd/¹⁴⁴Nd of 0.513835 ± 0.000053 (E_Nd t = + 7.4). This age is consistent with the precise Early Cretaceous age of Hauterivian (132 ± 2 to 127 ± 1.6 Ma) to Aptian (121 ± 1.4 to 112 ±1.1 Ma) for the overlying volcano-sedimentary (radiolarian bearing chert) sequences based on well-preserved radiolarian fossils (Kojima, S., Ahmad, T., Tanaka, T., Bagati, T.N., Mishra, M., Kumar, R. Islam, R., Khanna, P.P., 2001. Early Cretaceous radiolarians from the Indus suture zone, Ladakh, northern India. In: News of Osaka Micropaleontologists (NOM), Spec. Vol., 12, 257–270.) and cooling ages of 110–130 Ma based on ³⁹Ar/⁴⁰Ar for Nidar–Spontang ophiolitic rocks (Mahéo, G., Berttrand, H., Guillot, S., Villa, I. M., Keller, F., Capiez, P., 2004. The South Ladakh Ophiolites (NW Himalaya, India): an intra-oceanic tholeiitic arc origin with implications for the closure of the Neo-Tethys. Chem. Geol., 203, 273–303.). As these gabbroic and volcanic rocks are interpreted to be arc related, the new Sm–Nd age data may indicate that intra-ocean subduction in the Neo-Tethyan ocean may have started much before  140 ± 32 Ma as this date is interpreted as the age of crystallization of the arc magma. Present and published age data on the arc magmatic rocks from the Indus suture zone may collectively indicate episodic magmatism with increasing maturity of the arc from more basic (during ~ 140 ± 32 Ma) when the arc was immature through intermediate (andesitic/granodioritic) at ~ 100 Ma to more felsic (rhyolitic/dioritic) magmatism at ~ 50–45 Ma, when the Indian and the Asian plates collided.     

Assessment of agricultural drought in rainfed cereal production areas of northern China

Agricultural drought assessment is an important tool for water management in water-scarce regions such as Inner Mongolia and northeastern China. Conventional methods have difficulty of clarifying long-term influences of drought on regional agricultural production. To accurately evaluate regional agricultural drought, we assessed the performance of drought indices by constructing a new assessment framework with three components: crop model calibration and validation, drought index calculation, and index assessment (standard period setting, mean value and agreement assessments). The Environmental Policy Integrated Climate (EPIC) model simulated well of county-level wheat and maize yields in the nine investigated counties. We calculated a major crop-specific index yield reduction caused by water stress (WSYR) in the EPIC crop model, by relating potential and rainfed yields. Using 26 agricultural drought cases, we compared WSYR with two meteorological drought indices: precipitation (P) and aridity index (AI). The results showed that WSYR had greater agreement (85 %) than either the precipitation (65 %) or aridity index (68 %). The temporal trend of the indices over the period 1962–2010 was tested using three approaches. The result via WSYR revealed a significant increase in the trend of agricultural drought in drought-prone counties, which could not be shown by the precipitation and aridity indices. Total number of dry year via WSYR from 1990s to 2000s increases more sharply than via P or AI. As shown by WSYR, the number of dry years in northeastern China and Inner Mongolia is generally increasing, particularly after the 2000s, in the western part of the study area. The study reveals the usefulness of the framework for drought index assessment and indicates the potential of WSYR and possible drought cases for drought classification.     

Storm runoff processes in a small forested drainage basin

A brief review was made of storm runoff processes and the mechanisms of its generation in relation to subsurface water behaviors in a small forested drainage basin located in the western suburbs of Tokyo, Japan. The results of field investigations showed that the main source of storm runoff was groundwater flow and that the rapid and large amounts of groundwater discharge during a storm event could not be explained solely by the traditional concept of Darcian matrix flow. Several mechanisms such as pipe flow, air pressure effect, and capillary barrier effect were recognized that would induce a rapid response of groundwater to storm events depending on differences in local hydrologic conditions. All of these mechanisms were chiefly attributed to inhomogeneities of the soil deposits. The importance of dynamic behaviors of subsurface water during a storm event was emphasized in considering the mechanism of storm runoff generation.     

Dynamic Modelling of Fault Slip Induced by Stress Waves due to Stope Production Blasts

Seismic events can take place due to the interaction of stress waves induced by stope production blasts with faults located in close proximity to stopes. The occurrence of such seismic events needs to be controlled to ensure the safety of the mine operators and the underground mine workings. This paper presents the results of a dynamic numerical modelling study of fault slip induced by stress waves resulting from stope production blasts. First, the calibration of a numerical model having a single blast hole is performed using a charge weight scaling law to determine blast pressure and damping coefficient of the rockmass. Subsequently, a numerical model of a typical Canadian metal mine encompassing a fault parallel to a tabular ore deposit is constructed, and the simulation of stope extraction sequence is carried out with static analyses until the fault exhibits slip burst conditions. At that point, the dynamic analysis begins by applying the calibrated blast pressure to the stope wall in the form of velocities generated by the blast holes. It is shown from the results obtained from the dynamic analysis that the stress waves reflected on the fault create a drop of normal stresses acting on the fault, which produces a reduction in shear stresses while resulting in fault slip. The influence of blast sequences on the behaviour of the fault is also examined assuming several types of blast sequences. Comparison of the blast sequence simulation results indicates that performing simultaneous blasts symmetrically induces the same level of seismic events as separate blasts, although seismic energy is more rapidly released when blasts are performed symmetrically. On the other hand when nine blast holes are blasted simultaneously, a large seismic event is induced, compared to the other two blasts. It is concluded that the separate blasts might be employed under the adopted geological conditions. The developed methodology and procedure to arrive at an ideal blast sequence can be applied to other mines where faults are found in the vicinity of stopes.     

Preparation of fine particles of astrophysical interest by gas evaporation techniques

Fine particles of various chemical substances—carbon, iron, iron oxide and silica—which are expected to exist in interstellar space are prepared in argon gas. The size, shape, and crystal structure of the powder particles are studied by electron microscopy and electron diffraction, and may have some bearing on models of comets and of the solar nebula.The largest size we have obtained is about 0.1 μm and the smallest about 40Å. Generally the size becomes smaller for lower temperature and also for lower pressure. This tendency is discussed in connection with nucleation theories. Except for iron and iron oxide the powder particles have no crystal structure.     

Horizontal electric fields in the middle latitude

Three dimensional electric fields were measured at the altitude of about 27 km in the stratosphere over the Pacific Ocean about 200–400 km away from the Sanriku coast of Honsyu Island (L = 1·4) on 16–17 October 1973, which was magnetically disturbed. The average horizontal electric field thus measured is about 10 mV/m, and the electric field vectors made clockwise semidiurnal rotations rather than diurnal. Daily variation of this electric field was compared with data at L = 2·7–3·5 published by Mozer (1973) and was found to be very similar. This suggests that these electric fields are of common origin in the plasmasphere. From their mean daily variation it is estimated that the plasmaspheric convection is decreased in the night side and is increased in the day side by 200–300 m/sec, and there is an outward flow in the first half of the afternoon and an inward flow in the plasma bulge region of about 500 m/sec.     

Rare-earth elements in matrix,inclusions, and chondrules of the Allende meteorite

Rare-earth elements in a whole-rock sample and in major components of the Allende meteorite were investigated; for a few samples, abundances of Ba, Sr, Ca, and Al were also determined. Of the materials investigated in the present work, CaAl-rich inclusions G and O seem to be of the greatest significance. In spite of the minor difference in mineralogy between them, the apparent chondrite-normalized RE pattern is much different between these two inclusions. (Yb and Eu in inclusion G appear exceptionally irregular). This observation is inferred to reflect a rather subtle difference in condition of condensation. It is also worthwhile to note that, while two portions (pink and white) of the inclusion G show similar aspects in the abundances of lithophile trace elements investigated, they show a remarkable difference at the same time. The white portion (G_w) of inclusion G can be considered to be a mixture of chondritic material and highly fractionated material like the faintly pink portion (G_p) picked from the same inclusion. This would suggest the possibility that the G_p-like material was produced from chondritic dust.The “matrix” separated from Allende was found to be fractionated with respect to the RE abundances relative to representative chondrite. It has also a very high value for the Ba abundance.     

Sedimentary records of multidecadal-scale variability of diatom productivity in the Bungo Channel, Japan, associated with the Pacific Decadal Oscillation

In order to examine the responses of primary productivity in the southern coastal sea of Japan to the Pacific Decadal Oscillation (PDO) in the 20th century, sedimentary records of diatom productivity (diatom valve fluxes) were reconstructed using core samples from the Bungo Channel (BC) in southwest Japan. The record of the Thalassionema spp. flux—the best index of fall primary productivity in the BC—indicated a multidecadal-scale duration with a low flux (1943–1982) and those with a high flux (1913–1943 and 1982–2001); apparent shifts were recognized in 1943 and 1982. The shift in 1982 was also recognized in the flux records of other early summer to fall predominant genera in the BC and, previously, in the biogenic silica records from a broad region of the southeast BC. This indicates that in our records, this shift reflects a general trend in the primary production in the southeast BC. A comparison among the Thalassionema spp. flux records, meteorological data from an observatory adjacent to the core site, and the PDO index showed that the flux records were more similar to the PDO index than the other meteorological records, which suggests that the multidecadal-scale variability of the BC primary productivity may be associated with some marine-derived forcing. The bottom intrusions of nutrient-rich water that upwelled from the shelf slope into the BC, the axis movement or the transport of the Kuroshio Current off the BC, and a basin-scale wind stress in the North Pacific might play an important role in this forcing and mediate between the BC primary productivity and the PDO.     

Stationary Rossby waves and shocks on the Sverdrup coordinate

Gyre-scale frontal structures, e.g., the Subtropical Front, have been well documented in the oceans. Although the generation mechanism of such a front remains unclear, it may be ascribed to the steepening of nonlinear planetary waves as discussed by Dewar (1992). To discuss the stationary wave characteristics and their shock formation in a two-layer wind-driven gyre, the present paper introduces a new coordinate, referred to as theSverdrup coordinate here, in which the Sverdrup function is used instead of the longitudinal coordinate, and especially, investigates the possibility of the frontgenesis caused by the Rossby waves emanating from the western boundary region. On theSverdrup coordinate, since the advection by the Sverdrup flow is in the direction normal to that of the Rossby wave propagation, the system becomes much simpler than that on the (x,y)-coordinate, and the solution to this coordinate does not explicitly depend on the distribution of the Ekman pumping, i.e., we can treat cases, in which the Ekman pumping is a function of bothx andy, in a similar way to the case with zonally uniform Ekman pumping.