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.     

Population displacement after cyclone and its consequences: empirical evidence from coastal Bangladesh

This study investigates the disaster-induced population displacement scenario at individual household level in Bangladesh. ‘Population displacement’ is seen here as an alternative adaptation option to natural hazard for the survivors after cyclone Aila. The changes both in origin and destination community due to population displacement are described here on the basis of social ‘inclusion’ and ‘exclusion’ concept. The field survey was conducted during March–July 2010, and a sample of 280 respondents from 12 villages in southwest coastal Bangladesh was interviewed. Findings show that at the end of emergency aid, male members of the family started moving towards nearer cities to find an income. Based on the income and asset distribution at the community level, this study developed a societal cluster of displacement and demonstrates the societal changes because of cyclone-induced population displacement.     

Provenance of the Bosnian Flysch

Sandwiched between the Adriatic Carbonate Platform and the Dinaride Ophiolite Zone, the Bosnian Flysch forms a c. 3000 m thick, intensely folded stack of Upper Jurassic to Cretaceous mixed carbonate and siliciclastic sediments in the Dinarides. New petrographic, heavy mineral, zircon U/Pb and fission-track data as well as biostratigraphic evidence allow us to reconstruct the palaeogeology of the source areas of the Bosnian Flysch basin in late Mesozoic times. Middle Jurassic intraoceanic subduction of the Neotethys was shortly followed by exhumation of the overriding oceanic plate. Trench sedimentation was controlled by a dual sediment supply from the sub-ophiolitic high-grade metamorphic soles and from the distal continental margin of the Adriatic plate. Following obduction onto Adria, from the Jurassic–Cretaceous transition onwards a vast clastic wedge (Vranduk Formation) was developed in front of the leading edge, fed by continental basement units of Adria that experienced Early Cretaceous synsedimentary cooling, by the overlying ophiolitic thrust sheets and by redeposited elements of coeval Urgonian facies reefs grown on the thrust wedge complex. Following mid-Cretaceous deformation and thermal overprint of the Vranduk Formation, the depozone migrated further towards SW and received increasing amounts of redeposited carbonate detritus released from the Adriatic Carbonate Platform margin (Ugar Formation). Subordinate siliciclastic source components indicate changing source rocks on the upper plate, with ophiolites becoming subordinate. The zone of the continental basement previously affected by the Late Jurassic–Early Cretaceous thermal imprint has been removed; instead, the basement mostly supplied detritus with a wide range of pre-Jurassic cooling ages. However, a c. 80 Ma, largely synsedimentary cooling event is also recorded by the Ugar Formation, that contrasts the predominantly Early Cretaceous cooling of the Adriatic basement and suggests, at least locally, a fast exhumation.     

Polymetamorphic relations in iron ores from the Iron Quadrangle,Brazil: The correlation of oxygen isotope variations with deformation history

This study deals with the oxygen isotope composition of hematite-rich ore bodies in the Iron Quadrangle, Brazil. The area studied can be divided into two different regions: a western (W) zone of greenschist facies assemblages and an eastern (E) zone of amphibolite facies with transitions into granulite facies assemblages.The ¹⁸O values of 136 quartz-iron oxide pairs have been determined and temperatures of formation have been calculated. The ¹⁸O values of quartz vary between +6 and 20 except one value near +23, whereas the iron oxides fall between –4 and +10, with nearly 80% of the iron oxide values between –0.5 and 4.0. The regional distribution of the ¹⁸O values is as follows: in the W-region 85% of the quartz are >12, whereas in the E-region only 46% fall in this range, In contrast to quartz the iron oxides do not show any regional differences.The variation of oxygen isotope fractionations between quartz and iron oxides is obviously related to the complex deformation history of the iron ores. Samples with a primary schistosity (S₁) only represent peak metamorphic conditions. In the E-region the (S₁) high temperatures >700° C seem to correspond to orogenic events in the Archaen basement 2,700 m.y. ago. In the W-region S₁-temperatures between 460° and 560° C seem to represent peak metamorphic conditions of the Proterozoic Minas metamorphism 2,000 m.y. ago. Iron ores which have been overprinted by later deformation events are selectively reset to lower isotopic temperatures. The more closely spaced the schistosity planes the larger the extent of a temperature lowering.The genetic processes associated with these hematite-rich ore bodies appear to be sedimentary-metamorphic rather than metasomatic processes. Furthermore, there is no evidence for secondary leaching by weathering solutions.     

Defining organominerals: Comment on ‘Defining biominerals and organominerals: Direct and indirect indicators of life’ by Perry et al. (2007, Sedimentary Geology, 201, 157–179)

The paper by Perry et al. (2007, Defining biominerals and organominerals: Direct and indirect indicators of life, Sedimentary Geology, 201, 157–179) proposes to introduce “the new term ‘organomineral’” to describe mineral products whose formation is induced by by-products of biological activity, dead and decaying organisms, or nonbiological organic compounds, to be distinguished from the biomineral components of living organisms. The substantive ‘organomineral’, however, is not new: it was first introduced in 1993, with basically the same definition and distinction from biominerals, at the 7th International Symposium on Biomineralization (Défarge and Trichet, 1995, From biominerals to ‘organominerals’: The example of the modern lacustrine calcareous stromatolites from Polynesian atolls, Bulletin de l'Institut Océanographique de Monaco, n° spéc. 14, vol. 2, pp. 265–271). Thereafter, more than twenty-five papers by various authors have been devoted to organominerals and organomineral formation (‘organomineralization’) processes. Only two of these papers are cited by Perry et al., and without any reference to the definitions, or even the terms ‘organomineral’ or ‘organomineralization’, which they included. Moreover, Perry et al. tend to enlarge the original concept of organomineral to encompass all minerals containing organic matter, whether these organic compounds are active or passive in the mineralization, which introduces ambiguities detrimental to a fine understanding of present and past geobiological processes. Finally, Perry et al. propose to consider organominerals as indirect biosignatures that could be used in the search for evidence of life in the geological record and extraterrestrial bodies. This latter proposition also is problematical, in that organominerals may be formed in association with prebiotic or abiotic organic matter.     

Experimental growth of åkermanite reaction rims between wollastonite and monticellite: evidence for volume diffusion control

Growth rates of monomineralic, polycrystalline åkermanite (Ca₂MgSi₂O₇) rims produced by solid-state reactions between monticellite (CaMgSiO₄) and wollastonite (CaSiO₃) single crystals were determined at 0.5 GPa dry argon pressure, 1,000–1,200°C and 5 min to 60 h, using an internally heated pressure vessel. Inert Pt-markers, initially placed at the monticellite–wollastonite interface, indicate symmetrical growth into both directions. This and mass balance considerations demonstrate that rim growth is controlled by transport of MgO. At 1,200°C and run durations between 5 min and 60 h, rim growth follows a parabolic rate law with rim widths ranging from 0.4 to 16.3 μm indicating diffusion-controlled rim growth. The effective bulk diffusion coefficient \( D_{\text{eff,MgO}}^{\text{Ak}} \) is calculated to 10^?15.8±0.1 m² s^?1. Between 1,000°C and 1,200°C, the effective bulk diffusion coefficient follows an Arrhenius law with E _a = 204 ± 18 kJ/mol and D ₀ = 10^?8.6±1.6 m² s^?1. Åkermanite grains display a palisade texture with elongation perpendicular to the reaction interface. At 1,200°C, average grain widths measured normal to elongation, increase with the square root of time and range from 0.4 to 5.4 μm leading to a successive decrease in the grain boundary area fraction, which, however, does not affect \( D_{\text{eff,MgO}}^{\text{Ak}} \) to a detectible extent. This implies that grain boundary diffusion only accounts for a minor fraction of the overall chemical mass transfer, and rim growth is essentially controlled by volume diffusion. This is corroborated by the agreement between our estimates of the effective MgO bulk diffusion coefficient and experimentally determined volume diffusion data for Mg and O in åkermanite from the literature. There is sharp contrast to the MgO–SiO₂ binary system, where grain boundary diffusion controls rim growth.     

Mixing and chemical interdiffusion of trachytic and latitic magma in a subvolcanic complex of the Tertiary Westerwald (Germany)

A special kind of magma mixing is extraordinarily well exposed in the Bittersberg subvolcanic complex in the Tertiary volcanic field of the German Westerwald: A trachytic melt has been penetrated by a latitic dyke which has been dispersed within the host magma as small spherical enclaves (globules). Whole rock analyses of the globules show a change in composition that cannot be explained by a simple mechanical mixing between the endmembers. The most evolved globules have a phonolitic composition. Microprobe measurements in the microlithic matrix of the host rock and the guest indicate a diffusive motion of the alkalis from the host into the globules. On the other hand, an opposite trend can be observed for Ca, Mg, Fe and Ti, which are impoverished in the globules. The trace elements and the middle rare earth elements (MREE) has also been involved in the diffusive exchange. The REE-pattern of the most evolved (phonolitic) globules shows a characteristic trough in the area of the MREE which is almost identical to the REE-pattern of many phonolites. The phonolites and the alkali-rich trachytes of the Westerwald show similar globular textures as the Bittersberg volcanics. Therefore, generation of these rocks involving diffusive element exchange during mixing processes in a magma reservoir situated on a deeper crustal level may be possible.     

A high precision U–Pb age of metamorphic rutile in coesite-bearing eclogite from the Dabie Mountains in central China: a new constraint on the cooling history

This paper first reports a high precision U–Pb age of 218±1.2 Ma for rutile in coesite-bearing eclogite from Jinheqiao in the Dabie Mounteins, east–central China. This work shows that the U–Pb mineral (rutile+omphacite) isochron age of 218±2.5 Ma and conventional rutile U–Pb concordia age of 218±1.2 Ma obtained by common Pb correction based on the Pb isotopic composition of omphacite in the same eclogite sample are consistent, proving that the omphacite with low U/Pb ratio (μ=2.8) can be used for common Pb correction in U–Pb dating of rutile. Oxygen isotope analysis of rutile aliquots gave the consistent δ¹⁸O values of −6.1±0.1%, demonstrating oxygen isotope homogenization in the rutile of different grains as inclusion in garnet and grain in matrix. Oxygen isotope thermometry yields temperatures of 695±35 and 460±15 °C for quartz–garnet and quartz–rutile pairs, respectively. These oxygen isotopic observations suggest that the diffusion of oxygen in rutile as inclusion in garnet is not controlled by garnet. According to field-based thermochronological studies of rutile, an estimate of the T_c of about 460 °C for U–Pb system in rutile under rapid cooling conditions (20 °C/Ma) was advised. Based on this U–Pb age as well as the reported chronological data with their corresponding metamorphic and/or closure temperature, an improved T–t path has been constructed. The T–t path confirms that the UHPM rocks in South Dabie experienced a rapid cooling following the peak metamorphism before 220 Ma and a long isothermal stage from 213 to 180 Ma around 425 °C.     

The impact of shale gas on the costs of climate policy

The use of shale gas is commonly considered as a low-cost option for meeting ambitious climate policy targets. This article explores global and country-specific effects of increasing global shale gas exploitation on the energy markets, on greenhouse gas emissions, and on mitigation costs. The global techno-economic partial equilibrium model POLES (Prospective Outlook on Long-term Energy Systems) is employed to compare policies which limit global warming to 2°C and baseline scenarios when the availability of shale gas is either high or low. According to the simulation results, a high availability of shale gas has rather small effects on the costs of meeting climate targets in the medium and long term. In the long term, a higher availability of shale gas increases baseline emissions of greenhouse gases for most countries and for the world, and leads to higher compliance costs for most, but not all, countries. Allowing for global trading of emission certificates does not alter these general results. In sum, these findings cast doubt on shale gas’s potential as a low-cost option for meeting ambitious global climate targets.

POLICY RELEVANCE

Many countries with a large shale gas resource base consider the expansion of local shale gas extraction as an option to reduce their GHG emissions. The findings in this article imply that a higher availability of shale gas in these countries might actually increase emissions and mitigation costs for these countries and also for the world. An increase in shale gas extraction may spur a switch from coal to gas electricity generation, thus lowering emissions. At the global level and for many countries, though, this effect is more than offset by a crowding out of renewable and nuclear energy carriers, and by lower energy prices, leading to higher emissions and higher mitigation costs in turn. These findings would warrant a re-evaluation of the climate strategy in most countries relying on the exploitation of shale gas to meet their climate targets.