Origin of fluids in iron oxide–copper–gold deposits: constraints from δ 37Cl, 87Sr/86Sri and Cl/Br

The origin of the hypersaline fluids (magmatic or basinal brine?), associated with iron oxide (Cu–U–Au–REE) deposits, is controversial. We report the first chlorine and strontium isotope data combined with Cl/Br ratios of fluid inclusions from selected iron oxide–copper–gold (IOCG) deposits (Candelaria, Raúl–Condestable, Sossego), a deposit considered to represent a magmatic end member of the IOCG class of deposit (Gameleira), and a magnetite–apatite deposit (El Romeral) from South America. Our data indicate mixing of a high δ ³⁷Cl magmatic fluid with near 0‰ δ ³⁷Cl basinal brines in the Candelaria, Raúl–Condestable, and Sossego IOCG deposits and leaching of a few weight percent of evaporites by magmatic-hydrothermal (?) fluids at Gameleira and El Romeral. The Sr isotopic composition of the inclusion fluids of Candelaria, Raúl–Condestable, and El Romeral confirms the presence of a non-magmatic fluid component in these deposits. The heavy chlorine isotope signatures of fluids from the IOCG deposits (Candelaria, Raúl–Condestable, Sossego), reflecting the magmatic-hydrothermal component of these fluids, contrast with the near 0‰ δ ³⁷Cl values of porphyry copper fluids known from the literature. The heavy chlorine isotope compositions of fluids of the investigated IOCG deposits may indicate a prevailing mantle Cl component in contrast to porphyry copper fluids, an argument also supported by Os isotopes, or could result from differential Cl isotope fractionation processes (e.g. phase separation) in fluids of IOCG and porphyry Cu deposits.     

Deccan volcanism linked to the Cretaceous-Tertiary boundary mass extinction: New evidence from ONGC wells in the Krishna-Godavari Basin

A scientific challenge is to assess the role of Deccan volcanism in the Cretaceous-Tertiary boundary (KTB) mass extinction. Here we report on the stratigraphy and biologic effects of Deccan volcanism in eleven deep wells from the Krishna-Godavari (K-G) Basin, Andhra Pradesh, India. In these wells, two phases of Deccan volcanism record the world’s largest and longest lava mega-flows interbedded in marine sediments in the K-G Basin about 1500 km from the main Deccan volcanic province. The main phase-2 eruptions (∼80% of total Deccan Traps) began in C29r and ended at or near the KTB, an interval that spans planktic foraminiferal zones CF1–CF2 and most of the nannofossil Micula prinsii zone, and is correlative with the rapid global warming and subsequent cooling near the end of the Maastrichtian. The mass extinction began in phase-2 preceding the first of four mega-flows. Planktic foraminifera suffered a 50% drop in species richness. Survivors suffered another 50% drop after the first mega-flow, leaving just 7 to 8 survivor species. No recovery occurred between the next three mega-flows and the mass extinction was complete with the last phase-2 mega-flow at the KTB. The mass extinction was likely the consequence of rapid and massive volcanic CO₂ and SO₂ gas emissions, leading to high continental weathering rates, global warming, cooling, acid rains, ocean acidification and a carbon crisis in the marine environment.     

Eccentric geomagnetic dipole drift

A recalculation of the rate of westward drift of the eccentric geomagnetic dipole since 1900 has confirmed its correlation with irregularities in the excess length-of-day (LOD) as originally reported by Vestine. The major changes in the LOD curve since 1900 are shown to be well correlated with the westward drift rate from a number of magnetic models covering the interval. A new Magsat model of secular variation for 1980.0 shows no increase in the drift rate which dropped by a factor of three in the interval 1965–1970.     

Recent progress in the operational forecasting of summer severe weather

Abstract

Summer severe weather (SSW) can strike suddenly and unexpectedly with disastrous consequences for human activity. Considerable progress has been made in the past ten years in the operational forecasting of SSW. Traditionally, SSW was defined to consist of tornadoes, strong winds, hail, lightning and heavy rain. Hazardous types of strong winds have recently been expanded to include microbursts, macrobursts and surfacing rear inflow jet damage behind mesoscale convective systems. Doppler radar was used to relate surface damage to the appropriate atmospheric phenomena, first diagnostically and then prognostically. This improvement in classification has fedback to and improved the forecast process. Concurrent progress has been made in the use of synoptic observations. The concept of helical wind profiles and improved knowledge of the role of dry mid‐level air has improved the forecasting of tornadoes and strong gusty winds. Moisture flux convergence, derived from surface measurements, shows great promise in identifying areas of storm initiation. Satellite imagery has been used to identify dynamical atmospheric boundaries. Numerical modelling of the interaction of environmental wind profiles and individual thunderstorms has greatly contributed to the understanding of SSW. Studies of spatial and temporal patterns of lightning, both specific cases and climatology, contribute to the forecasting of severe storms. Polarization radar results have shown progress in separating the signals of hail from those of rain and in the improved measurement of heavy rainfalls. Radar observation of clear air boundaries and their interactions show potential for the forecasting of thunderstorm initiation. Though not traditionally considered part of SSW, hurricanes that evolve into extra‐tropical storms share many of the same hazardous features. The progress in computing, communications and display technologies has also made substantial contributions to operational forecasting and to the dissemination of weather warnings.     

Clustering and visualization of earthquake data in a grid environment

We present a web client-server service WEB-IS, which we have developed for remote analysis and visualization of seismic data consisting of both small magnitude events and large earthquakes. We show that the problem-solving environment (PSE) intended for prediction of large magnitude earthquakes can be based on this WEB-IS idea. The clustering schemes, feature generation, feature extraction techniques and rendering algorithms form a computational framework of this environment. On the other hand, easy and fast access both to the seismic data distributed among distant computing resources and to computational and visualization resources can be realized in a GRID framework. We discuss the usefulness of NaradaBrokering (iNtegrated Asynchronous Real-time Adaptive Distributed Architecture) as a middleware, allowing for flexibility and high throughput for remote visualization of geophysical data. The WEB-IS functionality was tested both on synthetic and the actual earthquake catalogs. We consider the application of similar methodology for tsunami alerts.     

A comparison between cokriging and ordinary kriging: Case study with a polymetallic deposit

Cokriging is applied to the estimation of mineral resources in a polymetallic deposit. Several major steps, which should be taken in using cokriging, are highlighted as necessary practical considerations. The case study is related to an ultramafic copper-nickel deposit. Six elements, Cu, Ni, Au, Ag, Pt, and Pd, occurring in the deposit, are partitioned into three subgroups and the elements within each group are simultaneously estimated on the basis of over 4000 drill assays. A comparison was made between ordinary kriging and cokriging methods through cross-validation. The results show that cokriging has significantly improved the estimates of resources by reducing the overall estimation error by over 15% and the variance of error by over 20%.     

Grid cell distortion and MODFLOW's integrated finite-difference numerical solution

The ground water flow model MODFLOW inherently implements a nongeneralized integrated finite-difference (IFD) numerical scheme. The IFD numerical scheme allows for construction of finite-difference model grids with curvilinear (piecewise linear) rows. The resulting grid comprises model cells in the shape of trapezoids and is distorted in comparison to a traditional MODFLOW finite-difference grid. A version of MODFLOW-88 (herein referred to as MODFLOW IFD) with the code adapted to make the one-dimensional DELR and DELC arrays two dimensional, so that equivalent conductance between distorted grid cells can be calculated, is described. MODFLOW IFD is used to inspect the sensitivity of the numerical head and velocity solutions to the level of distortion in trapezoidal grid cells within a converging radial flow domain. A test problem designed for the analysis implements a grid oriented such that flow is parallel to columns with converging widths. The sensitivity analysis demonstrates MODFLOW IFD's capacity to numerically derive a head solution and resulting intercell volumetric flow when the internal calculation of equivalent conductance accounts for the distortion of the grid cells. The sensitivity of the velocity solution to grid cell distortion indicates criteria for distorted grid design. In the radial flow test problem described, the numerical head solution is not sensitive to grid cell distortion. The accuracy of the velocity solution is sensitive to cell distortion with error <1% if the angle between the nonparallel sides of trapezoidal cells is <12.5 degrees. The error of the velocity solution is related to the degree to which the spatial discretization of a curve is approximated with piecewise linear segments. Curvilinear finite-difference grid construction adds versatility to spatial discretization of the flow domain. MODFLOW-88's inherent IFD numerical scheme and the test problem results imply that more recent versions of MODFLOW 2000, with minor modifications, have the potential to make use of a curvilinear grid.