Mineralogical and chemical evolution of the Ernest Henry Fe oxide–Cu–Au ore system, Cloncurry district, northwest Queensland, Australia

The Ernest Henry Cu–Au deposit was formed within a zoned, post-peak metamorphic hydrothermal system that overprinted metamorphosed dacite, andesite and diorite (ca 1740–1660 Ma). The Ernest Henry hydrothermal system was formed by two cycles of sodic and potassic alteration where biotite–magnetite alteration produced in the first cycle formed ca 1514±24 Ma, whereas paragenetically later Na–Ca veining formed ca 1529 +11/−8 Ma. These new U–Pb_titanite age dates support textural evidence for incursion of hydrothermal fluids after the metamorphic peak, and overlap with earlier estimates for the timing of Cu–Au mineralization (ca 1540–1500 Ma). A distal to proximal potassic alteration zone correlates with a large (up to 1.5 km) K–Fe–Mn–Ba enriched alteration zone that overprints earlier sodic alteration. Mass balance analysis indicates that K–Fe–Mn–Ba alteration—largely produced during pre-ore biotite- and magnetite-rich alteration—is associated with K–Rb–Cl–Ba–Fe–Mn and As enrichment and Na, Ca and Sr depletion. The aforementioned chemical exchange almost precisely counterbalances the mass changes associated with regional Na–Ca alteration. This initial transition from sodic to potassic alteration may have been formed during the evolution of a single fluid that evolved via alkali exchange during progressive fluid-rock interaction. Cu–Au ore, dominated by co-precipitated magnetite, minor specular hematite, and chalcopyrite as breccia matrix, forms a pipe-like body at the core of a proximal alteration zone dominated by K-feldspar alteration. Both the core and K-feldspar alteration overprint Na–Ca alteration and biotite–magnetite (K–Fe) alteration. Ore was associated with the concentration of a diverse range of elements (e.g. Cu, Au, Fe, Mo, U, Sb, W, Sn, Bi, Ag, F, REE, K, S, As, Co, Ba and Ca). Mineralization also involved the deposition of significant barite, K(–Ba)–feldspar, calcite, fluorite and complexly zoned pyrite. The complexly zoned pyrite and variable K–(Ba)–feldspar versus barite associations are interpreted to indicate fluctuating sulphur and/or barium supply. Together with the alteration zonation geochemistry and overprinting criteria, these data are interpreted to indicate that Cu–Au mineralization occurred as a result of fluid mixing during dilation and brecciation, in the location of the most intense initial potassic alteration. A link between early alteration (Na–Ca and K–Fe) and the later K-feldspathization and the Cu–Au ore is possible. However, the ore-related enrichments in particular elements (especially Ba, Mn, As, Mo, Ag, U, Sb and Bi) are so extreme compared with earlier alteration that another fluid, possibly magmatic in origin, contributed the diverse element suite geochemically independently of the earlier stages. Structural focussing of successive stages produced the distinctive alteration zoning, providing a basis both for exploration for similar deposits, and for an understanding of ore genesis.     

High-T, low-P metamorphism in the Palaeoproterozoic Halls Creek Orogen, northern Australia: the middle crustal response to a mantle-related transient thermal pulse

High‐T, low‐P metamorphic rocks of the Palaeoproterozoic central Halls Creek Orogen in northern Australia are characterised by low radiogenic heat production, high upper crustal thermal gradients (locally exceeding 40 °C km^?1) sustained for over 30 Myr, and a large number of layered mafic‐ultramafic intrusions with mantle‐related geochemical signatures. In order to account for this combination of geological and thermal characteristics, we model the middle crustal response to a transient mantle‐related heat pulse resulting from a temporary reduction in the thickness of the mantle lithosphere. This mechanism has the potential to raise mid‐crustal temperatures by 150–400 °C within 10–20 Myr following initiation of the mantle temperature anomaly, via conductive dissipation through the crust. The magnitude and timing of maximum temperatures attained depend strongly on the proximity, duration and lateral extent of the thermal anomaly in the mantle lithosphere, and decrease sharply in response to anomalies that are seated deeper than 50–60 km, maintained for <5 Myr in duration and/or have half‐widths <100 km. Maximum temperatures are also intimately linked to the thermal properties of the model crust, primarily due to their influence on the steady‐state (background) thermal gradient. The amplitudes of temperature increases in the crust are principally a function of depth, and are broadly independent of crustal thermal parameters. Mid‐crustal felsic and mafic plutonism is a predictable consequence of perturbed thermal regimes in the mantle and the lowermost crust, and the advection of voluminous magmas has the potential to raise temperatures in the middle crust very quickly. Although pluton‐related thermal signatures significantly dissipate within <10 Myr (even for very large, high‐temperature intrusive bodies), the interaction of pluton‐ and mantle‐related thermal effects has the potential to maintain host rock temperatures in excess of 400–450 °C for up to 30 Myr in some parts of the mid‐crust. The numerical models presented here support the notion that transient mantle‐related heat sources have the capacity to contribute significantly to the thermal budget of metamorphism in high‐T, low‐P metamorphic belts, especially in those characterised by low surface heat flow, very high peak metamorphic geothermal gradients and abundant mafic intrusions.     

Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.     

Total nitrogen content of deep sea basalts

Various models have been suggested concerning the origin and evolution of the earth's atmosphere. An estimate of the nitrogen content of the mantle could further constrain atmospheric models. Total nitrogen content was determined by thermal neutron activation analysis via ¹⁴N(n,p)¹⁴C. The ¹⁴C was converted to carbon dioxide and counted in miniature low level proportional counters. The total nitrogen content of U.S.G.S. standards BCR-1 and G-2 as determined by different laboratories is variable, probably due to atmospheric adsorption by the finely ground samples. Total nitrogen content was determined in deep sea basalt glasses from three regions: East Pacific Rise ($\text{15 ± 4, 18 ± 4}$, and $\text{7 ±}\text{ppm}\text{2 N}$), Mid-Atlantic Rift (FAMOUS Region:$\text{22 ± 5, 18 ± 3}$, and $\text{10 ± 2}\text{ppm}\text{N}$) and the Juan de Fuca Ridge ($\text{17 ± 4}\text{ppm}\text{N}$). Matrix material from the same samples as the glasses was available from the East Pacific Rise ($\text{37 ± 6, 26 ± 4}$, and $\text{34 ± 6}\text{ppm}\text{N}$) and the Mid-Atlantic Rift ($\text{39 ± 4}\text{ppm}\text{N}$) which are about 50 to 100% greater than the associated glasses. The increased matrix abundance may be due to incorporation of chemically bound nitrogen from sea water rather than dissolved molecular nitrogen. The nitrogen content of the FAMOUS samples are inconsistent with the model of Langmuir et al. (1977) for petrogenesis based on trace element data. Factors which can affect the observed nitrogen content in the basalts and the interpretation in terms of the mantle nitrogen abundance are discussed (e.g. partial melting and degassing of the basalts). A lower limit of about 2 ppm N in the mantle can be estimated.     

Present and Future Observing Trends in Atmospheric Magnetoseismology

With modern imaging and spectral instruments observing in the visible, EUV, X-ray, and radio wavelengths, the detection of oscillations in the solar outer atmosphere has become a routine event. These oscillations are considered to be the signatures of a wave phenomenon and are generally interpreted in terms of magnetohydrodynamic (MHD) waves. With multiwavelength observations from ground- and space-based instruments, it has been possible to detect waves in a number of different wavelengths simultaneously and, consequently, to study their propagation properties. Observed MHD waves propagating from the lower solar atmosphere into the higher regions of the magnetized corona have the potential to provide excellent insight into the physical processes at work at the coupling point between these different regions of the Sun. High-resolution wave observations combined with forward MHD modeling can give an unprecedented insight into the connectivity of the magnetized solar atmosphere, which further provides us with a realistic chance to reconstruct the structure of the magnetic field in the solar atmosphere. This type of solar exploration has been termed atmospheric magnetoseismology. In this review we will summarize some new trends in the observational study of waves and oscillations, discussing their origin and their propagation through the atmosphere. In particular, we will focus on waves and oscillations in open magnetic structures (e.g., solar plumes) and closed magnetic structures (e.g., loops and prominences), where there have been a number of observational highlights in the past few years. Furthermore, we will address observations of waves in filament fibrils allied with a better characterization of their propagating and damping properties, the detection of prominence oscillations in UV lines, and the renewed interest in large-amplitude, quickly attenuated, prominence oscillations, caused by flare or explosive phenomena.     

Calculation of the Inverse of the Covariance

In reservoir characterization, the covariance is often used to describe the spatial correlation and variation in rock properties or the uncertainty in rock properties. The inverse of the covariance, on the other hand, is seldom discussed in geostatistics. In this paper, I show that the inverse is required for simulation and estimation of Gaussian random fields, and that it can be identified with the differential operator in regularized inverse theory. Unfortunately, because the covariance matrix for parameters in reservoir models can be extremely large, calculation of the inverse can be a problem. In this paper, I discuss four methods of calculating the inverse of the covariance, two of which are analytical, and two of which are purely numerical. By taking advantage of the assumed stationarity of the covariance, none of the methods require inversion of the full covariance matrix.     

Characterization of Compass M-1 signals

An analysis of observations from China’s first medium earth orbit satellite Compass M-1 is presented, with main focus on the first orbit and clock solution for this satellite. The orbit is computed from laser ranging measurements. Based on this orbit solution, the apparent clock offset is estimated using measurements from two GNSS receivers, which allow Compass tracking. The analysis of the clock solutions reveals unexpectedly high dynamics in the pseudorange and carrier-phase observations. Furthermore, carrier-to-noise density ratio, pseudorange noise, and multipath are analyzed and compared to GPS and GIOVE. The results of the clock analysis motivate further research on the signals of the geostationary satellites of the Compass constellation.     

GNSS satellite transmit power and its impact on orbit determination

Antenna thrust is a small acceleration acting on Global Navigation Satellite System satellites caused by the transmission of radio navigation signals. Knowledge about the transmit power and the mass of the satellites is required for the computation of this effect. The actual transmit power can be obtained from measurements with a high-gain antenna and knowledge about the properties of the transmit and receive antennas as well as losses along the propagation path. Transmit power measurements for different types of GPS, GLONASS, Galileo, and BeiDou-2 satellites were taken with a 30-m dish antenna of the German Aerospace Center (DLR) located at its ground station in Weilheim. For GPS, total L-band transmit power levels of 50–240 W were obtained, 20–135 W for GLONASS, 95–265 W for Galileo, and 130–185 W for BeiDou-2. The transmit power differs usually only slightly for individual spacecraft within one satellite block. An exception are the GLONASS-M satellites where six subgroups with different transmit power levels could be identified. Considering the antenna thrust in precise orbit determination of GNSS satellites decreases the orbital radius by 1–27 mm depending on the transmit power, the satellite mass, and the orbital period.     

Evaluation of magmatic processes for the products of the Nevado del Ruiz Volcano, Colombia from geochemical and petrological data

Trace-element concentrations of the eruptive products from the Nevado del Ruiz volcano indicate that at least for the last 1 Ma, the major controlling factor in the evolution of the magma has been a simple fractional crystallization.The Zr/Hf ratios suggest that the source material is mantle-like in origin and whilst subduction-related contamination of the source material has taken place, there is no evidence that assimilation of continental crust has affected the series.No convincing geochemical evidence for magma mixing was found, and data on banded pumices of Pleistocene, Historic and the 1985 eruptions strongly support a crystal fractionation model to explain the compositional variation of their glasses. Therefore, if mixing has taken place it can only be a mixing of layers within a stratified or partially stratified magma chamber.Modeling of fractional crystallization using a stepwise program has been undertaken and least-squares approximations with small residuals are commensurate with analyzed samples. Crystal fractionation models with plagioclase dominating two-pyroxenes and iron oxides give good fits for all calculations. We conclude that simple crystal fractionation is the process most consistent with major- and trace-element variations in the Ruiz series. This closed-system model requires at least 77% crystal fractionation of a basaltic parent to generate the observed compositions of dacites.