Efficient block preconditioners for the coupled equations of pressure and deformation in highly discontinuous media

Large‐scale simulations of flow in deformable porous media require efficient iterative methods for solving the involved systems of linear algebraic equations. Construction of efficient iterative methods is particularly challenging in problems with large jumps in material properties, which is often the case in geological applications, such as basin evolution at regional scales. The success of iterative methods for this type of problems depends strongly on finding effective preconditioners. This paper investigates how the block‐structured matrix system arising from single‐phase flow in elastic porous media should be preconditioned, in particular for highly discontinuous permeability and significant jumps in elastic properties. The most promising preconditioner combines algebraic multigrid with a Schur complement‐based exact block decomposition. The paper compares numerous block preconditioners with the aim of providing guidelines on how to formulate efficient preconditioners. Copyright © 2010 John Wiley & Sons, Ltd.     

Simulating stochastically excited oscillations: The mode lifetime of ξ Hya

The discovery of solar-like oscillations in the giant star ξ Hya (G7 III) was reported by Frandsen et al. (2002). Their frequency analysis was very limited due to alias problems in the data set (caused by single-site observations). The extent to which the aliasing affected their analysis was unclear due to the unknown damping time of the stellar oscillation modes. In this paper we describe a simulator created to generate time series of stochastically excited oscillations, which takes as input an arbitrary window function and includes both white and non-white noise. We also outline a new method to compare a large number of simulated time series with an observed time series to determine the damping time, amplitude, and limited information on the degree of the stochastically excited modes. For ξ Hya we find the most likely amplitude to be ∼ 2 m s⁻¹, in good agreement with theory (Houdek and Gough, 2002), and the most likely damping time to be ∼ 2 days, which is much shorter than the theoretical value of 15–20 days calculated by Houdek and Gough (2002).     

The influence of central black holes on gravitational lenses

Recent observations indicate that many if not all galaxies host massive central black holes. In this paper we explore the influence of black holes on the lensing properties. We model the lens as an isothermal ellipsoid with a finite core radius plus a central black hole. We show that the presence of the black hole substantially changes the critical curves and caustics. If the black hole mass is above a critical value, then it will completely suppress the central images for all source positions. Realistic central black holes are likely to have masses below this critical value. Even in such subcritical cases, the black hole can suppress the central image when the source is inside a zone of influence, which depends on the core radius and black hole mass. In the subcritical cases, an additional image may be created by the black hole in some regions, which for some radio lenses may be detectable with high-resolution and large dynamic range VLBI maps. The presence of central black holes should also be taken into account when one constrains the core radius from the lack of central images in gravitational lenses.     

Radio-silent isolated neutron stars as a new astronomical reality

Summary. As of today, seven X-ray sources have been tentatively identified as radio-quiet, isolated neutron stars. The family appears to be a rapidly growing one, although not all the objects have been identified with the same degree of certainty. The most convincing example of radio quiet pulsar is certainly Geminga, the neutron star nature of which, proposed in 1983 on the basis of its similarity with the Vela pulsar, has been firmly established with the discovery of its X and pulsation. Four more neutron star candidates, originally found in the Einstein data, have been confirmed by ROSAT, which has added to the list two more entries. All this is not the result of an unbiased search. The seven sources were not selected at random: four are inside supernova remnants, an obvious place to search for isolated neutron stars, while the remaining three were singled out because of some peculiarity. Intense -ray emission in the case of Geminga, very high X-ray counting rate for RXJ185635-3754, or being the brightest unidentified source in the Einstein medium sensitivity survey, MS 0317-6647. In spite of the limited number of objects and of the observational biases, these seven radio quiet neutron star candidates add valuable pieces of information to the observational panorama of known pulsars. Their properties, inferred from the X-ray emission, offer a coherent picture, pointing towards thermally emitting, cooling neutron stars. Received: April 1, 1996     

Quaternary volcanic ash zones on the Iceland Plateau,southern Norwegian Sea

Four Quaternary volcanic ash zones in the southern Norwegian Sea have been investigated in core P57-7 from the Iceland Plateau. Both the geochemical composition and morphological variation of each ash layer have been studied. The four volcanic ash zones appear in the light oxygen isotope stages 1, 5, 7 and 11. The ash zones are composed of transparent platy grains, light brown transparent grains, brown blocky and black blocky grains and white/transparent pumice, and each zone shows a distinct stratigraphic evolution. The geochemical results show a mixture of basaltic and silicic grains in each ash zone, and that each zone contains grains from more than one eruption. The geochemical investigations strongly suggest that all the ashes are derived from Iceland. The youngest ash zone includes two layers, which based on their geochemical composition and stratigraphic position are correlated with the ¹⁴C-dated Vedde Ash (10 600 yr BP) and Saksunarvatn Ash (9 100 yr BP). Possible sources on Iceland for these layers are discussed.     

Relative contributions of sublimation and volcanoes to Io's atmosphere inferred from its plasma interaction during solar eclipse

We present a model that describes Io's delayed electrodynamic response to a temporal change in Io's atmosphere. Our model incorporates the relevant physical processes involved in Io's atmosphere-ionosphere-magnetosphere electrodynamic interaction to predict the far-ultraviolet (FUV) radiation as Io enters Jupiter's shadow and re-emerges into sunlight. The predicted FUV brightnesses are highly nonlinear as the strength of the electrodynamic interaction depends on the ratios of ionospheric conductances to the torus Alfvén conductance, but the former are functions of electrodynamics and the atmospheric density, which decays rapidly upon entering eclipse. Key factors governing the time evolution are the column density due to sublimation and the column density due to volcanoes, which maintain the background atmosphere during eclipse. The plasma interaction does not react instantaneously, but lags to a temporarily changing atmosphere. We find three qualitatively different scenarios with two of them including a post-eclipse brightening. The brightness ratio of in-sunlight/in-eclipse coupled with the existence of a sub-jovian equatorial spot constrains the volcanic column density to several times 10¹⁸ m⁻², based on the currently available observations. Thus in sunlight, the sublimation driven part of Io's atmosphere dominates the volcanically driven contribution by roughly a factor of 10 or more.     

Ensemble standardization constraints on the influence of the tree growth trends in dendroclimatology

Climate Dynamics - Tree growth trends can affect the interpretation of the response of tree-ring proxies (especially tree-ring width) to climate in the low-frequency band, which in turn may limit...     

40Ar‐39Ar step‐heating of impact glasses from the Nördlinger Ries impact crater—Implications on excess argon in impact melts and tektites

Seven impact melts from various places in the Nördlinger Ries were dated by ⁴⁰Ar‐³⁹Ar step‐heating. The aim of these measurements was to increase the age data base for Ries impact glasses directly from the Ries crater, because there is only one Ar‐Ar step‐heating spectrum available in the literature. Almost all samples display saddle‐shaped age spectra, indicating the presence of excess argon in most Ries glass samples, most probably inherited argon from incompletely degassed melt and possibly also excess argon incorporated during cooling from adjacent phases. In contrast, moldavites usually contain no inherited argon, probably due to their different formation process implying solidification during ballistic transport. The plateau age of the only flat spectrum is 14.60 ± 0.16 (0.20) Ma (2σ), while the total age of this sample is 14.86 ± 0.20 (0.22) Ma (isochron age: 14.72 ± 0.18 [0.22] Ma [2σ]), proofing the chronological relationship of the Ries impact and moldavites. The total ages of the other samples range between 15.77 ± 0.52 and 20.4 ± 1.0 Ma (2σ), implying approximately 2–40% excess ⁴⁰Ar (compared to the nominal age of the Ries crater) in respective samples. Thus, the age of 14.60 ± 0.16 (0.20) (2σ) (14.75 ± 0.16 [0.20 Ma] [2σ], calculated using the most recent suggestions for the K decay constants) can be considered as reliable and is within uncertainties indistinguishable from the most recent compilation for the age of the moldavite tektites.     

Optimization routine for identification of model parameters in soil plasticity

The paper presents an optimization routine especially developed for the identification of model parameters in soil plasticity on the basis of different soil tests. Main focus is put on the mathematical aspects and the experience from application of this optimization routine. Mathematically, for the optimization, an objective function and a search strategy are needed. Some alternative expressions for the objective function are formulated. They capture the overall soil behaviour and can be used in a simultaneous optimization against several laboratory tests. Two different search strategies, Rosenbrock's method and the Simplex method, both belonging to the category of direct search methods, are utilized in the routine. Direct search methods have generally proved to be reliable and their relative simplicity make them quite easy to program into workable codes. The Rosenbrock and simplex methods are modified to make the search strategies as efficient and user‐friendly as possible for the type of optimization problem addressed here. Since these search strategies are of a heuristic nature, which makes it difficult (or even impossible) to analyse their performance in a theoretical way, representative optimization examples against both simulated experimental results as well as performed triaxial tests are presented to show the efficiency of the optimization routine. From these examples, it has been concluded that the optimization routine is able to locate a minimum with a good accuracy, fast enough to be a very useful tool for identification of model parameters in soil plasticity. Copyright © 2001 John Wiley & Sons, Ltd.