Uncoupling of coupled flows in soil—a finite element method

Coupled flow of water, chemicals, heat and electrical potential in soil are of significance in a variety of circumstances. The problem is characterized by the coupling between different flows, i.e. a flow of one type driven by gradients of other types, and by the dual nature of certain flows, i.e. combined convection and conduction. Effective numerical solutions to the problem are challenged due to the coupling and the dual nature. In this paper, we first present a general expression that can be used to represent various types of coupled flows in soil. A finite element method is then proposed to solve the generalized coupled flows of convection-conduction pattern. The unknown vector is first decomposed into two parts, a convective part forming a hyperbolic system and a conductive part forming a parabolic system. At each time step, the hyperbolic system is solved analytically to give an initial solution. To solve the multi-dimensional hyperbolic system, we assume that a common eigenspace exists for the coefficient matrices, so that the system can be uncoupled by transforming the unknown vector to the common eigenspace. The uncoupled system is solved by the method of characteristics. Using the solution of the hyperbolic system as the initial condition, we then solve the parabolic system by a Galerkin finite element method for space discretization and a finite difference scheme for time stepping. The proposed technique can be used for solving multi-dimensional, transient, coupled or simultaneous flows of convection-conduction type. Application to a flow example shows that the technique indeed exhibits optimality in convergence and in stability.     

Alkaline fluid circulation in ultramafic rocks and formation of nucleotide constituents: a hypothesis

Seawater is constantly circulating through oceanic basement as a low-temperature hydrothermal fluid (<150°C). In cases when ultramafic rocks are exposed to the fluids, for instance during the initial phase of subduction, ferromagnesian minerals are altered in contact with the water, leading to high pH and formation of secondary magnesium hydroxide, among other – brucite, that may scavenge borate and phosphate from seawater. The high pH may promote abiotic formation of pentoses, particularly ribose. Pentoses are stabilized by borate, since cyclic pentoses form a less reactive complex with borate. Analyses have shown that borate occupies the 2' and 3' positions of ribose, thus leaving the 5' position available for reactions like phosphorylation. The purine coding elements (adenine, in particular) of RNA may be formed in the same general hydrothermal environments of the seafloor.     

The significance of rock structure,lithology and pre‐glacial deep weathering for the shape of intermediate‐scale glacial erosional landforms

The main landforms within the glacially scoured Precambrian rocks of the Swedish west coast are closely connected to the principal structural pattern and have lately been explained as mainly a result of etch processes, probably during the Mesozoic and with a possible second period of etching during the Neogene. To explore the effect of multiple glacial erosion on the rock surfaces, an island with two different lithologies and with striae from different directions was selected for a detailed study, focusing on the shape of roches moutonnées. Air‐photo interpretation of bedrock lineaments and roches moutonnées combined with detailed field mapping and striae measurements are used to interpret the structural and lithological control on the resulting shape. The study reveals a significant difference in shape between roches moutonnées in augen‐granite and orthogneiss. Low elongated and streamlined roches moutonnées occur in the gneiss area, striated by a Late Weichselian ice flow from the NE. This ice flow is subparallel with both the local dominant trend of topographically well‐expressed joints and the schistosity of the gneiss. Frequently, there are no signs of quarrying on the lee‐sides of the gneiss roches moutonnées and hence they resemble the shape of whalebacks, or ruwares, as typically associated with the exposed basal weathering surface found in tropical areas. The granite roches moutonnées were formed by an older ice flow from the ESE, which closely followed the etched WNW–ESE joint system of the granite. Late Weichselian ice flow from the NE caused only minor changes of the landforms. On the contrary, marks of the early ESE ice flow are poorly preserved in the gneiss area, where it probably never had any large effect as the flow was perpendicular to both schistosity and structures and, accordingly, also to the pre‐glacial relief. The study demonstrates that coincidence between ice flow direction and pre‐glacially etched structures is most likely to determine the effects of glacial erosion. Copyright © 2002 John Wiley & Sons, Ltd.     

The Late Cretaceous Klepa basalts in Macedonia (FYROM)—Constraints on the final stage of Tethys closure in the Balkans

The waning stage(s) of the Tethyan ocean(s) in the Balkans are not well understood. Controversy centres on the origin and life‐span of the Cretaceous Sava Zone, which is allegedly a remnant of the last oceanic domain in the Balkan Peninsula, defining the youngest suture between Eurasia‐ and Adria‐derived plates. In order to investigate to what extent Late‐Cretaceous volcanism within the Sava Zone is consistent with this model we present new age data together with trace‐element and Sr–Nd–Pb isotope data for the Klepa basaltic lavas from the central Balkan Peninsula. Our new geochemical data show marked differences between the Cretaceous Klepa basalts (Sava Zone) and the rocks of other volcanic sequences from the Jurassic ophiolites of the Balkans. The Klepa basalts mostly have Sr–Nd–Pb isotopic and trace‐element signatures that resemble enriched within‐plate basalts substantially different from Jurassic ophiolite basalts with MORB, BAB and IAV affinities. Trace‐element modelling of the Klepa rocks indicates 2%–20% polybaric melting of a relatively homogeneously metasomatised mantle source that ranges in composition from garnet lherzolite to ilmenite+apatite bearing spinel–amphibole lherzolite. Thus, the residual mineralogy is characteristic of a continental rather than oceanic lithospheric mantle source, suggesting an intracontinental within‐plate origin for the Klepa basalts. Two alternative geodynamic models are internally consistent with our new findings: (1) if the Sava Zone represents remnants of the youngest Neotethyan Ocean, magmatism along this zone would be situated within the forearc region and triggered by ridge subduction; (2) if the Sava Zone delimits a diffuse tectonic boundary between Adria and Europe which had already collided in the Late Jurassic, the Klepa basalts together with a number of other magmatic centres represent volcanism related to transtensional tectonics.     

Rutile R632 – A New Natural Reference Material for U‐Pb and Zr Determination

A new natural rutile reference material is presented, suitable for U‐Pb dating and Zr‐in‐rutile thermometry by microbeam methods. U‐Pb dating of rutile R632 using laser ablation ICP‐MS with both magnetic sector field and quadrupole instruments as well as isotope dilution‐thermal ionisation mass spectrometry yielded a concordia age of 496 ± 2 Ma. The high U content (> 300 μg g^?1) enabled measurement of high‐precision U‐Pb ages despite its young age. The sample was found to have a Zr content of 4294 ± 196 μg g^?1, which makes it an excellent complementary reference material for Zr‐in‐rutile thermometry. Individual rutile grains have homogeneous compositions of a number of other trace elements including V, Cr, Fe, Nb, Mo, Sn, Sb, Hf, Ta and W. This newly characterised material significantly expands the range of available rutile reference materials relevant for age and temperature determinations.     

Chemistry, reaction mechanisms and micro-structures during retrograde metamorphism of gedrite-biotite-plagioclase bearing rocks from Bergslagen, south-central Sweden

The Proterozoic rhyolitic volcanics constituting the foot-wall rocks in the Stollberg ore-field, Bergslagen, south-central Sweden, locally contain gedrite altered to chlorite and serpentine, biotite altered to chlorite and plagioclase altered to epidote.

The intergrowths between the host gedrite and the chlorite/serpentine inclusions are oriented with the a^* of gedrite parallel to the c^* of serpentine and chlorite. The biotite has been altered to chlorite by brucitization of both the K-interlayer and talc-like layer. In both cases the net change in volume during chloritization is small.

The assumption that Al is conserved during alteration of gedrite and biotite agrees very well with the micro-structures and orientation relations observed by transmission electron microscopy. Normalizing the chlorite to 1.00 mole, the overall chemical change that took place during the retrograde metamorphism of the Stollberg rocks can be written as: 0.84Ged+0.14Bio+0.65Mg+4.76H₂O+0.42H=1.00Chl+0.57Alb+0.72Fe+0.01Na+0.12K+0.01Ti+0.05Mn+0.95H₄SiO₄ The reaction results in ca 9% increase in volume for the solid phases. Thus, a slightly acidic Mg-rich fluid started the reaction and, upon leaving the system, the metasomatic fluid was enriched in Na, Fe, K, and Si.     

Nonlinear interaction between acoustic gravity waves

The resonant interaction between three acoustic gravity waves is considered. We improve on the results of previous authors and write the new coupling coefficients in a symmetric form. Particular attention is paid to the low-frequency limit.     

Effect of glacier loading/deloading on volcanism: postglacial volcanic production rate of the Dyngjufjöll area,central Iceland

Tephrochronological dating of postglacial volcanism in the Dyngjufjöll volcanic complex, a major spreading center in the Icelandic Rift Zone, indicates a high production rate in the millennia following deglaciation as compared to the present low productivity. The visible and implied evidence indicates that lava production in the period 10 000–4500 bp was at least 20 to 30 times higher than that in the period after 2900 bp but the results are biased towards lower values for lava volumes during the earlier age periods since multiple lava layers are buried beneath younger flows. The higher production rate during the earlier period coincides with the disappearance of glaciers of the last glaciation. Decreasing lithostatic pressure as the glacier melts and vigorous crustal movements caused by rapid isostatic rebound may trigger intense volcanism until a new pressure equilibrium has been established.