Evolution of melt composition during intrusion of basalts into a silicic magma chamber

The article describes heat exchange between basaltic and rhyolite melts accompanied by fractional crystallization of phases in a basaltic melt. A numerical model has been developed for the homogenization mechanism of magma composition during intrusion of basaltic magma batches into felsic magma chambers. The results of numerical modeling demonstrate that the time needed for cooling the basalts and their fractionation to rhyolite melts is much shorter than the time required for chemical interaction based on diffusive mechanisms.     

Mineralogy and composition of Kuroko deposits from northeastern Honshu and their possible modern analogues from the Izu-Ogasawara (Bonin) Arc south of Japan: Implications for mode of formation

The Kuroko deposits of NE Honshu are a key type deposit for the study of volcanogenic massive sulfide deposits. However, these deposits have not been studied in detail since the early 1980's and knowledge of their mode of formation is now dated. In this study, we present the analysis of 12 samples of the Kuroko deposits, 12 samples of submarine hydrothermal minerals from the Sunrise deposit and 6 samples from Suiyo Seamount, both of which are located on the Izu-Ogasawara (Bonin) Arc, for 27 elements. For the Kuroko deposit, Cd>Sb>Ag>Pb>Hg>As>Zn>Cu are highly enriched, Au>Te>Bi>Ba>Mo are moderately enriched, In>Tl are somewhat enriched and Fe is not significantly enriched relative to the average continental crust. Within each of these deposits, a similar pattern of element associations is apparent: Zn–Pb with As, Sb, Cd, Ag, Hg, Tl and Au; Fe–Cu–Ba with As, Sb, Ag, Tl, Mo, Te and Au; Si–Ba with Ag and Au; CaSO₄. The enrichment of the chalcophilic elements in these deposits is consistent with hydrothermal leaching of these elements from the host rocks which are dominantly rhyolite–dacite in the case of the Kuroko deposits, rhyolite in the case of the Sunrise deposit and dacite–rhyolite in the case of the Suiyo Seamount deposit. However, this pattern of element enrichment is also similar to that observed in fumarolic gas condensates from andesitic volcanoes. This suggests that there may be a significant magmatic contribution to the composition of the hydrothermal fluids responsible for the formation of the Kuroko deposits, although it is not yet possible to quantify the relative contributions of these two sources of elements.The compositional data show that Sunrise and Suiyo Seamount deposits are much closer compositionally to the Kuroko deposits from NE Honshu than are the submarine hydrothermal deposits from the JADE site in the Okinawa Trough which contain, on average, significantly higher concentrations of Pb, Zn, Sb, As and Ag than each of these deposits. In spite of the greater similarity in tectonic setting of the Hokuroku Basin in which the Kuroko deposits formed to the Okinawa Trough (intracontinental rifted back-arc basin) compared to Myojin Knoll and Suiyo Seamount (active arc volcanoes), it appears that submarine hydrothermal deposits from Myojin Knoll and Suiyo Seamount are closer analogues of the Kuroko deposit than are those from the Okinawa Trough. The present data are consistent with the magmatic hydrothermal model for the formation of Kuroko-type deposits as formulated by Urabe and Marumo [Urabe, T., Marumo, K., 1991. A new model for Kuroko-type deposits of Japan. Episodes 14, 246–251].     

Ultrahigh pressure macro diamonds from Copeton (New South Wales, Australia), based on Raman spectroscopy of inclusions

Mining of Cenozoic alluvial deposits at Copeton and Bingara (Eastern Australia) has produced two million macrodiamonds (0.25 ct median size). Raman spectroscopy is used to identify included minerals within uncut Copeton diamonds, with sealed chamber remnant pressures of 31.7 to 35.6 kbar for coesite, 13.6 and 22.7 kbar for clinopyroxene, and 7.6 kbar for grossular garnet. Assuming elastic behaviour, these values generate inclusion entrapment PT loci which intersect, restricting diamond formation conditions: from 250 °C, 43 kbar to 800 °C, 52 kbar. Larger than error (± 100 °C and ± 4 kbar), this range shows a systematic variation in inclusion composition with diamond zoning and N properties. Published research shows 1) Copeton and Bingara diamonds are unique, and 2) modern alluvium in the Bingara district carries mantle-formed garnet, captured by post-tectonic alkali basalt from an extensive diamondiferous ultrahigh pressure (UHP) terrane that stalled at depth because it is dominated by mafic eclogite. The combined Raman and geological results indicate two sets of subduction UHP diamond formation conditions/protolith are required, firstly cooler oceanic slab and secondly including higher temperature continental crust. The Copeton and Bingara stones are UHP macrodiamonds, and Carboniferous ⁴⁰Ar/³⁹Ar age dates on clinopyroxene inclusions should be interpreted as ages of crystallisation, representing the termination of subduction. The characteristic features of ruptured inclusions and etched percussion marks on Copeton and Bingara diamond indicate volcanic delivery to the earth's surface. Alluvial deposits elsewhere in Eastern Australia may carry similar diamond along with diamond of different origin.     

Organic carbon deposition and phosphorus accumulation during Oceanic Anoxic Event 2 in Tarfaya, Morocco

With a multi-proxy approach, an attempt was made to constrain productivity and bottom-water redox conditions and their effects on the phosphorus accumulation rate at the Mohammed Plage section on the Tarfaya coast, Morocco, during the Cenomanian-Turonian Anoxic Event (OAE 2). A distinct δ¹³C_org isotope excursion of +2.5‰ occurs close to the top of the section. The unusually abrupt shift of the isotope excursion and disappearance of several planktonic foraminiferal species (e.g. Rotalipora cushmani and Rotalipora greenhornensis) in this level suggests a hiatus of between 40–60 kyrs at the excursion onset. Nevertheless, it was possible to determine both the long-term environmental history as well as the processes that took place immediately prior to and during OAE 2. TOC% values increase gradually from the base of the section to the top (from 2.5% to 10%). This is interpreted as the consequence of a long-term eustatic sea-level rise and subsidence causing the encroachment of less oxic waters into the Tarfaya Basin. Similarly a reduction in the mineralogically constructed ‘detrital index’ can be explained by the decrease in the continental flux of terrigenous material due to a relative sea-level rise. A speciation of phosphorus in the upper part of the section, which spans the start and mid-stages of OAE 2, shows overall higher abundances of P_reactive mass accumulation rates before the isotope excursion onset and lower values during the plateau. Due to the probable short hiatus, the onset of the decrease in phosphorus content relative to the isotope excursion is uncertain, although the excursion plateau already contains lower concentrations. The C_org/P_total and V/Al ratios suggest that this reduction was mostly likely caused by a decrease in the available bottom oxygen content (probably as a result of higher productivity) and a corresponding fall in the phosphorus retention ability of the sediment. Productivity appears to have remained high during the isotope plateau possibly due to a combination of ocean-surface fertilisation via increased aridity (increased K/Al and Ti/Al ratios) and/or higher dissolved inorganic phosphorus content in the water column as a result of the decrease in sediment P retention. The evidence for decreased P-burial has been observed in many other palaeoenvironments during OAE 2. Tarfaya's unique upwelling paleosituation provides strong evidence that the nutrient recycling was a global phenomenon and therefore a critical factor in starting and sustaining OAE 2.     

Magnetotelluric investigations for imaging electrical structure of Garhwal Himalayan corridor,Uttarakhand, India

Magnetotelluric investigations have been carried out in the Garhwal Himalayan corridor to delineate the electrical structure of the crust along a profile extending from Indo-Gangetic Plain to Higher Himalayan region in Uttarakhand, India. The profile passing through major Himalayan thrusts: Himalayan Frontal Thrust (HFF), Main Boundary Thrust (MBT) and Main Central Thrust (MCT), is nearly perpendicular to the regional geological strike. Data processing and impedance analysis indicate that out of 44 stations MT data recorded, only 27 stations data show in general, the validity of 2D assumption. The average geoelectric strike, N70°W, was estimated for the profile using tensor decomposition. 2D smooth geoelectrical model has been presented, which provides the electrical image of the shallow and deeper crustal structure. The major features of the model are (i) a low resistivity (<50Ωm), shallow feature interpreted as sediments of Siwalik and Indo-Gangetic Plain, (ii) highly resistive (> 1000Ωm) zone below the sediments at a depth of 6 km, interpreted as the top surface of the Indian plate, (iii) a low resistivity (< 10Ωm) below the depth of 6 km near MCT zone coincides with the intense micro-seismic activity in the region. The zone is interpreted as the partial melting or fluid phase at mid crustal depth. Sensitivity test indicates that the major features of the geoelectrical model are relevant and desired by the MT data.     

Characteristics of spectral aerosol optical depths over India during ICARB

Spectral aerosol optical depth (AOD) measurements, carried out regularly from a network of observatories spread over the Indian mainland and adjoining islands in the Bay of Bengal and Arabian Sea, are used to examine the spatio-temporal and spectral variations during the period of ICARB (March to May 2006). The AODs and the derived Ångström parameters showed considerable variations across India during the above period. While at the southern peninsular stations the AODs decreased towards May after a peak in April, in the north Indian regions they increased continuously from March to May. The Ångström coefficients suggested enhanced coarse mode loading in the north Indian regions, compared to southern India. Nevertheless, as months progressed from March to May, the dominance of coarse mode aerosols increased in the columnar aerosol size spectrum over the entire Indian mainland, maintaining the regional distinctiveness. Compared to the above, the island stations showed considerably low AODs, so too the northeastern station Dibrugarh, indicating the prevalence of cleaner environment. Long-range transport of aerosols from tshe adjoining regions leads to remarkable changes in the magnitude of the AODs and their wavelength dependencies during March to May. HYSPLIT back-trajectory analysis shows that enhanced long-range transport of aerosols, particularly from the west Asia and northwest coastal India, contributed significantly to the enhancement of AOD and in the flattening of the spectra over entire regions; if it is the peninsular regions and the island Minicoy are more impacted in April, the north Indian regions including the Indo Gangetic Plain get affected the most during May, with the AODs soaring as high as 1.0 at 500 nm. Over the islands, the Ångström exponent (α) remained significantly lower (～1) over the Arabian Sea compared to Bay of Bengal (BoB) (～1.4) as revealed by the data respectively from Minicoy and Port Blair. Occurrences of higher values of α, showing dominance of accumulation mode aerosols, over BoB are associated well with the advection, above the boundary layer, of fine particles from the east Asian region during March and April. The change in the airmass to marine in May results in a rapid decrease in α over the BoB.     

Aerosol mass loading over the marine environment of Arabian Sea during ICARB: Sea-salt and non-sea-salt components

Mass loading and chemical composition of atmospheric aerosols over the Arabian Sea during the pre-monsoon months of April and May have been studied as a part of the Integrated Campaign for Aerosols, gases and Radiation Budget (ICARB). These investigations show large spatial variabilities in total aerosol mass loading as well as that of individual chemical species. The mass loading is found to vary between 3.5 and 69.2 μg m^?3, with higher loadings near the eastern and northern parts of Arabian Sea, which decreases steadily to reach its minimum value in the mid Arabian Sea. The decrease in mass loading from the coast of India towards west is estimated to have a linear gradient of 1.53 μg m^?3/° longitude and an e^?1 scale distance of ～2300 km. SO ₄ ^2? , Cl^? and Na⁺ are found to be the major ionic species present. Apart from these, other dominating watersoluble components of aerosols are NO ₃ ^? (17%) and Ca²⁺ (6%). Over the marine environment of Arabian Sea, the non-sea-salt component dominates accounting to ～76% of the total aerosol mass. The spatial variations of the various ions are examined in the light of prevailing meteorological conditions and airmass back trajectories.     

Urban effects of Chennai on sea breeze induced convection and precipitation

Doppler radar derived wind speed and direction profiles showed a well developed sea breeze circulation over the Chennai, India region on 28 June, 2003. Rainfall totals in excess of 100 mm resulted from convection along the sea breeze front. Inland propagation of the sea breeze front was observed in radar reflectivity imagery. High-resolution MM5 simulations were used to investigate the influence of Chennai urban land use on sea breeze initiated convection and precipitation. A comparison of observed and simulated 10m wind speed and direction over Chennai showed that the model was able to simulate the timing and strength of the sea breeze. Urban effects are shown to increase the near surface air temperature over Chennai by 3.0K during the early morning hours. The larger surface temperature gradient along the coast due to urban effects increased onshore flow by 4.0m s⁻¹. Model sensitivity study revealed that precipitation totals were enhanced by 25mm over a large region 150 km west of Chennai due to urban effects. Deficiency in model physics related to night-time forecasts are addressed.     

Seismic hazard assessment of Chennai city considering local site effects

Chennai city suffered moderate tremors during the 2001 Bhuj and Pondicherry earthquakes and the 2004 Sumatra earthquake. After the Bhuj earthquake, Indian Standard IS: 1893 was revised and Chennai city was upgraded from zone II to zone III which leads to a substantial increase of the design ground motion parameters. Therefore, a comprehensive study is carried out to assess the seismic hazard of Chennai city based on a deterministic approach. The seismicity and seismotectonic details within a 100 km radius of the study area have been considered. The one-dimensional ground response analysis was carried out for 38 representative sites by the equivalent linear method using the SHAKE91 program to estimate the ground motion parameters considering the local site effects. The shear wave velocity profile was inferred from the corrected blow counts and it was verified with the Multichannel Analysis of Surface Wave (MASW) test performed for a representative site. The seismic hazard is represented in terms of characteristic site period and Spectral Acceleration Ratio (SAR) contours for the entire city. It is found that structures with low natural period undergo significant amplification mostly in the central and southern parts of Chennai city due to the presence of deep soil sites with clayey or sandy deposits and the remaining parts undergo marginal amplification.     

The WRF model performance for the simulation of heavy precipitating events over Ahmedabad during August 2006

The summer monsoon season of the year 2006 was highlighted by an unprecedented number of monsoon lows over the central and the western parts of India, particularly giving widespread rainfall over Gujarat and Rajasthan. Ahmedabad had received 540.2mm of rainfall in the month of August 2006 against the climatological mean of 219.8mm. The two spells of very heavy rainfall of 108.4mm and 97.7mm were recorded on 8 and 12 August 2006 respectively. Due to meteorological complexities involved in replicating the rainfall occurrences over a region, the Weather Research and Forecast (WRF-ARW version) modeling system with two different cumulus schemes in a nested configuration is chosen for simulating these events. The spatial distributions of large-scale circulation and moisture fields have been simulated reasonably well in this model, though there are some spatial biases in the simulated rainfall pattern. The rainfall amount over Ahmedabad has been underestimated by both the cumulus parameterization schemes. The quantitative validation of the simulated rainfall is done by calculating the categorical skill scores like frequency bias, threat scores (TS) and equitable threat scores (ETS). In this case the KF scheme has outperformed the GD scheme for the low precipitation threshold.