Methodology and use of tensor invariants for satellite gravity gradiometry

Although its use is widespread in several other scientific disciplines, the theory of tensor invariants is only marginally adopted in gravity field modeling. We aim to close this gap by developing and applying the invariants approach for geopotential recovery. Gravitational tensor invariants are deduced from products of second-order derivatives of the gravitational potential. The benefit of the method presented arises from its independence of the gradiometer instrument’s orientation in space. Thus, we refrain from the classical methods for satellite gravity gradiometry analysis, i.e., in terms of individual gravity gradients, in favor of the alternative invariants approach. The invariants approach requires a tailored processing strategy. Firstly, the non-linear functionals with regard to the potential series expansion in spherical harmonics necessitates the linearization and iterative solution of the resulting least-squares problem. From the computational point of view, efficient linearization by means of perturbation theory has been adopted. It only requires the computation of reference gravity gradients. Secondly, the deduced pseudo-observations are composed of all the gravitational tensor elements, all of which require a comparable level of accuracy. Additionally, implementation of the invariants method for large data sets is a challenging task. We show the fundamentals of tensor invariants theory adapted to satellite gradiometry. With regard to the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite gradiometry mission, we demonstrate that the iterative parameter estimation process converges within only two iterations. Additionally, for the GOCE configuration, we show the invariants approach to be insensitive to the synthesis of unobserved gravity gradients.     

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.     

The enigma of helium and anomalous xenon

Isotopically anomalous xenon in chondrites is closely associated with low-Z noble gases, but there is no helium (or neon) in the noble gas component with normal xenon. The correlation of elemental and isotopic heterogeneities in meteoritic noble gases places stringent limits on the origin of isotopically anomalous elements in meteorites and on the formation of the solar system.     

Fe57 Mossbauer spectroscopy of the Leoville meteorite

The Mossbauer spectra has been measured for the Leoville meteorite. Experimental results identify it as a type III carbonaceous meteorite. The spectrum shows the principle iron components to be olivine, magnetite, metallic iron, troilite and phyllosilicate material.     

Proximity of galactic civilizations

Approximate expressions are derived for the number of civilizations within a few tens of light years of each other since intelligent life first evolved in the Galaxy. The number is proportional to the square of the usual selectivity factors and to the first power of the longevity. Arguments are presented for expecting intelligent life in certain multiple star systems, and the number of coexistent civilizations in such systems is estimated.     

Impact of Flood Disasters on Taiwan in the Last Quarter Century

The increasing natural disasters, especially floods during the last quarter century, are raising the economic losses in Taiwan. The most severe hazard in Taiwan is flooding induced by typhoons and storms in summer and autumn. By comparing the rivers around the world, the ones in Taiwan have the steepest slopes, the largest discharge per unit drainage area, and the shortest time of concentrations. Rapid urbanization without proper land uses managements usually worsen the flood problems. Consequently, flood hazards mitigation has become the most essential task for Taiwan to deal with. Although the government keeps improving flood defense structures, the flood damage grows continuously. In this article, possible flood mitigation strategies are identified for coping with complex environmental and social decisions with flood risk involved.     

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.     

Rapid orbit determination of LEO satellites using IGS clock and ephemeris products

Different types of GPS clock and orbit data provided by the International GPS Service (IGS) have been used to assess the accuracy of rapid orbit determination for satellites in low Earth orbit (LEO) using spaceborne GPS measurements. To avoid the need for reference measurements from ground-based reference receivers, the analysis is based on an undifferenced processing of GPS code and carrier-phase measurements. Special attention is therefore given to the quality of GPS clock data that directly affects the resulting orbit determination accuracy. Interpolation of clock data from the available 15 min grid points is identified as a limiting factor in the use of IGS ultra-rapid ephemerides. Despite this restriction, a 10-cm orbit determination accuracy can be obtained with these products data as demonstrated for the GRACE-B spacecraft during selected data arcs between 2002 and 2004. This performance may be compared with a 5-cm orbit determination accuracy achievable with IGS rapid and final products using 5 min clock samples. For improved accuracy, high-rate (30 s) clock solutions are recommended that are presently only available from individual IGS centers. Likewise, a reduced latency and more frequent updates of IGS ultra-rapid ephemerides are desirable to meet the requirements of upcoming satellite missions for near real-time and precise orbit determination.     

Evaluation of the Reduction in Uncertainty Obtained by Conditioning a 3D Stochastic Channel to Multiwell Pressure Data

A stochastic channel embedded in a background facies is conditioned to data observed at wells. The background facies is a fixed rectangular box. The model parameters consist of geometric parameters that describe the shape, size, and location of the channel, and permeability and porosity in the channel and nonchannel facies. We extend methodology previously developed to condition a stochastic channel to well-test pressure data, and well observations of the channel thickness and the depth of the top of the channel. The main objective of this work is to characterize the reduction in uncertainty in channel model parameters and predicted reservoir performance that can be achieved by conditioning to well-test pressure data at one or more wells. Multiple conditional realizations of the geometric parameters and rock properties are generated to evaluate the uncertainty in model parameters. The ensemble of predictions of reservoir performance generated from the suite of realizations provides a Monte Carlo estimate of the uncertainty in future performance predictions. In addition, we provide some insight on how prior variances, data measurement errors, and sensitivity coefficients interact to determine the reduction in model parameters obtained by conditioning to pressure data and examine the value of active and observation well data in resolving model parameters.     

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.