Improved structural interpretation of legacy 3D seismic data from Karee platinum mine (South Africa) through the application of novel seismic attributes

Seismic detection of faults, dykes, potholes and iron-rich ultramafic pegmatitic bodies is of great importance to the platinum mining industry, as these structures affect safety and efficiency. The application of conventional seismic attributes (such as instantaneous amplitude, phase and frequency) in the hard-rock environment is more challenging than in soft-rock settings because the geology is often complex, reflections disrupted and the seismic energy strongly scattered. We have developed new seismic attributes that sharpen seismic reflections, enabling additional structural information to be extracted from hard-rock seismic data. The symmetry attribute is based on the invariance of an object with respect to transformations such as rotation and reflection; it is independent of the trace reflection amplitude, and hence a better indicator of the lateral continuity of thin and weak reflections. The reflection-continuity detector attribute is based on the Hilbert transform; it enhances the visibility of the peaks and troughs of the seismic traces, and hence the continuity of weak reflections. We demonstrate the effectiveness of these new seismic attributes by applying them to a legacy 3D seismic data set from the Bushveld Complex in South Africa. These seismic attributes show good detection of deep-seated thin (∼1.5 m thick) platinum ore bodies and their associated complex geological structures (faults, dykes, potholes and iron-rich ultramafic pegmatites). They provide a fast, cost-effective and efficient interpretation tool that, when coupled with horizon-based seismic attributes, can reveal structures not seen in conventional interpretations.     

DMSP and DMS dynamics and microzooplankton grazing in the Labrador Sea: application of the dilution technique

We adapted the dilution technique to study microzooplankton grazing of algal dimethylsulfoniopropionate (DMSP) vs. Chl a, and to estimate the impact of microzooplankton grazing on dimethyl sulfide (DMS) production in the Labrador Sea. Phytoplankton numbers were dominated by autotrophic nanoflagellates in the Labrador basin, but diatoms and colonial Phaeocystis pouchetii contributed significantly to phytomass at several high chlorophyll stations and on the Newfoundland and Greenland shelfs. Throughout the region, growth of algal Chl a and DMSP was generally high (0.2–1 d⁻¹), but grazing rates were lower and more variable, characteristic of the early spring bloom period. Production and consumption of Chl a vs. DMSP followed no clear pattern, and sometimes diverged greatly, likely because of their differing distributions among algal prey taxa and size class. In several experiments where Phaeocystis was abundant, we observed DMS production proportional to grazing rate, and we found clear evidence of DMS production by this haptophyte following physical stress such as sparging or filtration. It is possible that grazing-activated DMSP cleavage by Phaeocystis contributes to grazer deterrence: protozoa and copepods apparently avoided healthy colonies (as judged by relative growth and grazing rates of Chl a and DMSP), and grazing of Phaeocystis was significant only at one station where cells were in poor condition. Although we hoped to examine selective grazing on or against DMSP-containing algal prey, the dilution technique cannot differentiate selective ingestion and varying digestion rates of Chl a and DMSP. We also found that the dilution method alone was poorly suited for assessing the impact of grazing on dissolved sulfur pools, because of rapid microbial consumption and the artifactual release of DMSP and DMS during filtration. Measuring and understanding the many processes affecting organosulfur cycling by the microbial food web in natural populations remain a technical challenge that will likely require a combination of techniques to address.     

A guide to TAURUS-2 Fabry–Perot data reduction

The data reduction process for optical emission-line observations of galaxies using the TAURUS-2 Fabry–Perot interferometer mounted on the 3.9-m Anglo-Australian Telescope is described in detail. The initial steps (bias subtraction, flat-fielding, etc.) are the same as for calibration of CCD images, and the wavelength calibration is similar to that in optical spectroscopy. The final steps are specific to Fabry–Perot instruments, and include the fitting of several instrumental parameters and a phase correction to convert the raw ( x ,  y ,  z ) data cube into a useful position–velocity ( α ,  δ ,  v ) cube. Software has been written to assist with the latter steps of the data reduction. H α observations of NGC 1808, NGC 2442 and Circinus are used to demonstrate the reduction process.     

Integrated Global System Model for Climate Policy Assessment: Feedbacks and Sensitivity Studies

Alternative policies to address global climate change are being debated in many nations and within the United Nations Framework Convention on Climate Change. To help provide objective and comprehensive analyses in support of this process, we have developed a model of the global climate system consisting of coupled sub-models of economic growth and associated emissions, natural fluxes, atmospheric chemistry, climate, and natural terrestrial ecosystems. The framework of this Integrated Global System Model is described and the results of sample runs and a sensitivity analysis are presented. This multi-component model addresses most of the major anthropogenic and natural processes involved in climate change and also is computationally efficient. As such, it can be used effectively to study parametric and structural uncertainty and to analyze the costs and impacts of many policy alternatives. Initial runs of the model have helped to define and quantify a number of feedbacks among the sub-models, and to elucidate the geographical variations in several variables that are relevant to climate science and policy. The effect of changes in climate and atmospheric carbon dioxide levels on the uptake of carbon and emissions of methane and nitrous oxide by land ecosystems is one potentially important feedback which has been identified. The sensitivity analysis has enabled preliminary assessment of the effects of uncertainty in the economic, atmospheric chemistry, and climate sub-models as they influence critical model results such as predictions of temperature, sea level, rainfall, and ecosystem productivity. We conclude that uncertainty regarding economic growth, technological change, deep oceanic circulation, aerosol radiative forcing, and cloud processes are important influences on these outputs.     

Seismic anisotropy in granite at the Underground Research Laboratory, Manitoba

The Shear-Wave Experiment at Atomic Energy of Canada Limited's Underground Research Laboratory was probably the first controlled-source shear-wave survey in a mine environment. Taking place in conjunction with the excavation of the Mine-by test tunnel at 420 m depth, the shear-wave experiment was designed to measure the in situ anisotropy of the rockmass and to use shear waves to observe excavation effects using the greatest variety of raypath directions of any in situ shear-wave survey to date. Inversion of the shear-wave polarizations shows that the anisotropy of the in situ rockmass is consistent with hexagonal symmetry with an approximate fabric orientation of strike 023° and dip 35°. The in situ anisotropy is probably due to microcracks with orientations governed by the in situ stress field and to mineral alignment within the weak gneissic layering. However, there is no unique interpretation as to the cause of the in situ anisotropy as the fabric orientation agrees approximately with both the orientation expected from extensive-dilatancy anisotropy and that of the gneissic layering. Eight raypaths with shear waves propagating wholly or almost wholly through granodiorite, rather than granite, do not show the expected shear-wave splitting and indicate a lower in situ anisotropy, which may be due to the finer grain size and/or the absence of gneissic layering within the granodiorite. These results suggest that shear waves may be used to determine crack and mineral orientations and for remote monitoring of a rockmass. This has potential applications in mining and waste monitoring.     

The heating of the temperature minimum region in solar flares — A reassessment

As a sequel to the work by Machado et al. (1978), we discuss and evaluate the suggestions made by these authors on how to possibly reconcile the observed temperature enhancements at temperature-minimum levels in solar flares with some form of theoretical heating mechanism. After establishing the H^– LTE assumption used by Machado et al., we then consider EUV irradiation, and joule heating by steady currents, as heating mechanisms. We find that, unless there are strong inhomogeneities associated with either mechanism, neither can reasonably be reconciled with observations. It is concluded that detailed, high resolution (both spatial and temporal) measurements are necessary to further our understanding of the flare process at temperature-minimum levels.On leave from: Department of Astronomy, The University, Glasgow G12 8QQ, Scotland, U. K.     

Studies in the application of recurrence relations to special perturbation methods

The procedure of numerical integration of the elliptic three dimensional restricted threebody problem by the use of recurrence relations to evaluate successively higher derivatives of the relative position and velocity vectors of the bodies and of the variational matrix is investigated. A set of recurrence relations is developed which involves the introduction of fewer auxiliary variables than in previous papers of this series, while the recurrence relations themselves are of a simpler form than those in other treatments involving the same number of such auxiliary variables. A technique for automatic adjustment of the integration step-length at each point in the orbit, such that the local truncation error remains close to, but always less than, some specified amount, is incorporated. This technique involves the comparison of pre-integration values with those obtained after consecutive forward and reverse integration steps, and has decided advantages over step-adjustment methods currently in use.Both these modifications to previous techniques are shown, by presentation of sample computational results, to represent considerable savings in machine time for a given calculation and desired accuracy; these savings are generally around a factor of two and become greater as the desired accuracy in the computations increases.     

A case study of impact‐induced hydrothermal activity: The Haughton impact structure,Devon Island,Canadian High Arctic

Abstract— The well‐preserved state and excellent exposure at the 39 Ma Haughton impact structure, 23 km in diameter, allows a clearer picture to be made of the nature and distribution of hydrothermal deposits within mid‐size complex impact craters. A moderate‐ to low‐temperature hydrothermal system was generated at Haughton by the interaction of groundwaters with the hot impact melt breccias that filled the interior of the crater. Four distinct settings and styles of hydrothermal mineralization are recognized at Haughton: a) vugs and veins within the impact melt breccias, with an increase in intensity of alteration towards the base; b) cementation of brecciated lithologies in the interior of the central uplift; c) intense veining around the heavily faulted and fractured outer margin of the central uplift; and d) hydrothermal pipe structures or gossans and mineralization along fault surfaces around the faulted crater rim. Each setting is associated with a different suite of hydrothermal minerals that were deposited at different stages in the development of the hydrothermal system. Minor, early quartz precipitation in the impact melt breccias was followed by the deposition of calcite and marcasite within cavities and fractures, plus minor celestite, barite, and fluorite. This occurred at temperatures of at least 200 °C and down to ?100–120 °C. Hydrothermal circulation through the faulted crater rim with the deposition of calcite, quartz, marcasite, and pyrite, occurred at similar temperatures. Quartz mineralization within breccias of the interior of the central uplift occurred in two distinct episodes (?250 down to ?90 °C, and <60 °C). With continued cooling (<90 °C), calcite and quartz were precipitated in vugs and veins within the impact melt breccias. Calcite veining around the outer margin of the central uplift occurred at temperatures of ?150 °C down to <60 °C. Mobilization of hydrocarbons from the country rocks occurred during formation of the higher temperature calcite veins (>80 °C). Appreciation of the structural features of impact craters has proven to be key to understanding the distribution of hydrothermal deposits at Haughton.     

Bayesian calibrated significance levels applied to the spectral tilt and hemispherical asymmetry

Bayesian model selection provides a formal method of determining the level of support for new parameters in a model. However, if there is not a specific enough underlying physical motivation for the new parameters it can be hard to assign them meaningful priors, an essential ingredient of Bayesian model selection. Here we look at methods maximizing the prior so as to work out what is the maximum support the data could give for the new parameters. If the maximum support is not high enough then one can confidently conclude that the new parameters are unnecessary without needing to worry that some other prior may make them significant. We discuss a computationally efficient means of doing this which involves mapping p-values on to upper bounds of the Bayes factor (or odds) for the new parameters. A p-value of 0.05 (1.96σ) corresponds to odds less than or equal to 5:2, which is below the 'weak' support at best threshold. A p-value of 0.0003 (3.6σ) corresponds to odds of less than or equal to 150:1, which is the 'strong' support at best threshold. Applying this method we find that the odds on the scalar spectral index being different from one are 49:1 at best. We also find that the odds that there is primordial hemispherical asymmetry in the cosmic microwave background are 9:1 at best.     

Thermal plasma enhancements in the topside ionosphere and their relationship to the auroral electrojets

Data from the Topside sounder aboard ISIS II are used to calculate the latitudinal distribution of thermal plasma at 1400 km when the satellite moves along a line of constant geographic longitude. The distribution of thermal plasma is related to the latitudinal regime of the auroral electrojets as inferred from data from a meridian line of magnetometers. It is found that thermal plasma at high altitude tends to be found above the poleward portion of the auroral electrojet. This finding is explained in terms of the spectrum of precipitating electrons across the auroral oval. It is found that the thermal plasma distribution in the post-noon sector is distinctly different from that in the premidnight sector despite the fact that both quadrants feature clear latitudinally confined eastward current flow. This difference can be used to define whether or not the polar cleft penetrates into the local time sector traversed by the satellite. On the nightside, the peak in thermal plasma poleward of the ionospheric trough can be used to identify the magnetic field lines which map to the boundary between the tail lobe and the plasma sheet.