New equipment and processing for magnetotelluric remote reference observations

Robust estimates of magnetotelluric and geomagnetic response functions are determined using the coherency and expected uniformity of the magnetic source field as quality criteria. The method is applied on data sets of three simultaneously recording sites. For the data acquisition we used a new generation of geophysical equipment (S.P.A.M. MkIII), which comprises novel concepts of parallel computing and networked, digital data transmission. The data-processing results show that the amount of noise on the horizontal components of the magnetic field varies considerably in time, between sites and over the frequency range. The removal of such contaminated data beforehand is essential for most data-processing schemes, as the magnetic channels are usually assumed to be free of noise. The standard remote reference method is aimed at reducing bias in response function estimates. However, this does not necessarily improve their precision as our results clearly show. With our method, on the other hand, we can filter out source field irregularities, thereby providing suitable working conditions for the robust algorithm, and eventually obtain considerably improved results. Contrary to previous concepts, we suggest rejecting as much data as feasible in order to concentrate on the remaining parts of high-quality observations.     

Optical and infrared diagnostics of SDSS galaxies in the SWIRE survey

We present the rest-frame optical and infrared colours of a complete sample of  1114 z < 0.3  galaxies from the Spitzer Wide-Area Infrared Extragalactic (SWIRE) Legacy Survey and the Sloan Digital Sky Survey (SDSS). We discuss the optical and infrared colours of our sample and analyse in detail the contribution of dusty star-forming galaxies and active galactic nuclei (AGN) to optically selected red sequence galaxies.
We propose that the optical  ( g − r )  colour and infrared  log( L ₂₄/ L _3.6)  colour of galaxies in our sample are determined primarily by a bulge-to-disc ratio. The  ( g − r )  colour is found to be sensitive to the bulge-to-disc ratio for disc-dominated galaxies, whereas the  log( L ₂₄/ L _3.6)  colour is more sensitive for bulge-dominated systems.
We identify ∼18 per cent (195 sources) of our sample as having red optical colours and infrared excess. Typically, the infrared luminosities of these galaxies are found to be at the high end of star-forming galaxies with blue optical colours. Using emission-line diagnostic diagrams, 78 are found to have an AGN contribution and 117 are identified as star-forming systems. The red  ( g − r )  colour of the star-forming galaxies could be explained by extinction. However, their high optical luminosities cannot. We conclude that they have a significant bulge component.
The number densities of optically red star-forming galaxies are found to correspond to ∼13 per cent of the total number density of our sample. In addition, these systems contribute ∼13 per cent of the total optical luminosity density, and 28 per cent of the total infrared luminosity density of our SWIRE/SDSS sample. These objects may reduce the need for 'dry mergers'.     

Multiscale Gravitational Field Recovery from GPS-Satellite-to-Satellite Tracking

The purpose of GPS-satellite-to-satellite tracking (GPS-SST) is to determine the gravitational potential at the earth's surface from measured ranges (geometrical distances) between a low-flying satellite and the high-flying satellites of the Global Positioning System (GPS). In this paper, GPS-satellite-to-satellite tracking is reformulated as the problem of determining the gravitational potential of the earth from given gradients at satellite altitude. The uniqueness and stability of the solution are investigated. The essential tool is to split the gradient field into a normal part (i.e. the first-order radial derivative) and a tangential part (i.e. the surface gradient). Uniqueness is proved for polar, circular orbits corresponding to both types of data (first radial derivative and/or surface gradient). In both cases gravity recovery based on satellite-to-satellite tracking turns out to be an exponentially ill-posed problem. Regularization in terms of spherical wavelets is proposed as an appropriate solution method, based on the knowledge of the singular system. Finally, the extension of this method is generalized to a nonspherical earth and a non-spherical orbital surface, based on combined terrestrial and satellite data.     

Construction of the granitoid crust of an island arc. Part II: a quantitative petrogenetic model

Results of simple model calculations that integrate cumulate compositions from the Kohistan arc terrain are presented in order to develop a consistent petrogenetic model to explain the Kohistan island arc granitoids. The model allows a quantitative approximation of the possible relative roles of fractional crystallization and assimilation to explain the silica-rich upper crust composition of oceanic arcs. Depending in detail on the parental magma composition hydrous moderate-to-high pressure fractional crystallization in the lower crust/upper mantle is an adequate upper continental crust forming mechanism in terms of volume and compositions. Accordingly, assimilation and partial melting in the lower crust is not per se a necessary process to explain island arc granitoids. However, deriving few percent of melts using low degree of dehydration melting is a crucial process to produce volumetrically important amounts of upper continental crust from silica-poorer parental magmas. Even though the model can explain the silica-rich upper crustal composition of the Kohistan, the fractionation model does not predict the accepted composition of the bulk continental crust. This finding supports the idea that additional crustal refining mechanism (e.g., delamination of lower crustal rocks) and/or non-cogenetic magmatic process were critical to create the bulk continental crust composition.     

The electrical conductivity of Posidonia black shales — from magnetotelluric exploration to rock samples

We carried out a magnetotelluric field campaign in the South–East Lower Saxony Basin, Germany, with the main goal of testing this method for imaging regional Posidonia black shale sediments. Two‐dimensional inversion results of the magnetotelluric data show a series of conductive structures correlating with brine‐saturated sediments but also with deeper, anthracitic Westphalian/Namurian coals. None of these structures can be directly related with the Posidonia black shale, which appears to be generally resistive and therefore difficult to resolve with the magnetotelluric method. This assumption is supported by measurements of electrical resistivity on a set of Posidonia shale samples from the Hils syncline in the Lower Saxony basin. These rock samples were collected in shallow boreholes and show immature (0.53% Ro), oil (0.88% Ro), and gas (1.45% Ro) window thermal maturities. None of the black shale samples showed low electrical resistivity, particularly those with oil window maturity show resistivity exceeding 10⁴ Ωm. Moreover, we could not observe a direct correlation between maturity and electrical resistivity; the Harderode samples showed the highest resistivity, whereas the Haddessen samples showed the lowest. A similar trend has been seen for coals in different states of thermal maturation. Saturation of the samples with distilled and saline water solutions led to decreasing electrical resistivity. Moreover, a positive correlation of electrical resistivity with porosity is observed for the Wickensen and Harderode samples, which suggests that the electrical resistivity of the Posidonia black shale is mainly controlled by porosity.     

Forward induced seismic hazard assessment: application to a synthetic seismicity catalogue from hydraulic stimulation modelling

The M _w 3.2-induced seismic event in 2006 due to fluid injection at the Basel geothermal site in Switzerland was the starting point for an ongoing discussion in Europe on the potential risk of hydraulic stimulation in general. In particular, further development of mitigation strategies of induced seismic events of economic concern became a hot topic in geosciences and geoengineering. Here, we present a workflow to assess the hazard of induced seismicity in terms of occurrence rate of induced seismic events. The workflow is called Forward Induced Seismic Hazard Assessment (FISHA) as it combines the results of forward hydromechanical-numerical models with methods of time-dependent probabilistic seismic hazard assessment. To exemplify FISHA, we use simulations of four different fluid injection types with various injection parameters, i.e. injection rate, duration and style of injection. The hydromechanical-numerical model applied in this study represents a geothermal reservoir with preexisting fractures where a routine of viscous fluid flow in porous media is implemented from which flow and pressure driven failures of rock matrix and preexisting fractures are simulated, and corresponding seismic moment magnitudes are computed. The resulting synthetic catalogues of induced seismicity, including event location, occurrence time and magnitude, are used to calibrate the magnitude completeness M _c and the parameters a and b of the frequency-magnitude relation. These are used to estimate the time-dependent occurrence rate of induced seismic events for each fluid injection scenario. In contrast to other mitigation strategies that rely on real-time data or already obtained catalogues, we can perform various synthetic experiments with the same initial conditions. Thus, the advantage of FISHA is that it can quantify hazard from numerical experiments and recommend a priori a stimulation type that lowers the occurrence rate of induced seismic events. The FISHA workflow is rather general and not limited to the hydromechanical-numerical model used in this study and can therefore be applied to other fluid injection models.     

Final analysis of ELAIS 15-μm observations: method, reduction and catalogue

We present the final analysis of the European Large Area Infrared Space Observatory ( ISO ) Survey (ELAIS) 15-μm observations, carried out with the ISO Camera (ISOCAM) instrument on board the ISO .
The data-reduction method, known as the Lari Method, is based on a mathematical model of the behaviour of the detector and was specifically designed for the detection of faint sources in ISOCAM/ISO Photopolarimeter (ISOPHOT) data. The method is fully interactive and leads to very reliable and complete source lists.
The resulting catalogue includes 1923 sources detected with signal-to-noise ratio of  > 5  in the 0.5–100 mJy flux range and over an area of 10.85 deg² split into four fields, making it the largest non-serendipitous extragalactic source catalogue obtained to date from the ISO data.
This paper presents the concepts underlying the data-reduction method together with its latest enhancements. The data-reduction process, the production and basic properties of the resulting catalogue are discussed. The catalogue quality is assessed by means of detailed simulations, optical identifications and comparison with previous analyses.     

In-flight performance analysis of the CHAMP BlackJack GPS Receiver

JPL's BlackJack receiver currently represents the most widely used geodetic grade GPS receiver for space applications. Using data from the CHAMP science mission, the in-flight performance of the BlackJack receiver has been assessed and the impact of various software updates performed during the 2.5 years since launch is described. Key aspects of the study comprise the channel allocation, anomalous data points, and the noise level of the code and carrier data. In addition, it has been demonstrated that the code measurements collected onboard the CHAMP satellite are notably affected by multipath errors in the aft-looking hemisphere, which can be attributed to cross-talk between the occultation antenna string and the primary precise orbit determination antenna. For carrier smoothed 10 s normal points, the code noise itself varies between a minimum of 5 cm at high elevations and 0.5 m (C/A) to 1.0 m (P₁, P₂) at 10° elevation. Carrier-phase data exhibit representative errors of 0.2 to 2.5 mm. The results of the CHAMP GPS data analysis contribute to a better understanding and possible improvement of the BlackJack receiver and support the design of optimal data editing and weighting strategies in precise orbit determination applications.     

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.