Towards realistic simulations of lava dome growth using the level set method

The level set method has been implemented in a computational volcanology context. New techniques are presented to solve the advection equation and the reinitialisation equation. These techniques are based upon an algorithm developed in the finite difference context, but are modified to take advantage of the robustness of the finite element method. The resulting algorithm is tested on a well documented Rayleigh–Taylor instability benchmark [19], and on an axisymmetric problem where the analytical solution is known. Finally, the algorithm is applied to a basic study of lava dome growth.

Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.

Lattice-dome design using a knowledge-based system approach

In the design of lattice domes, design engineers need expertise in areas such as configuration processing, nonlinear analysis, and optimization. These are extensive numerical, iterative, and lime-consuming processes that are prone to error without an integrated design tool. This article presents the application of a knowledge-based system in solving lattice-dome design problems. An operational prototype knowledge-based system, LADOME, has been developed by employing the combined knowledge representation approach, which uses rules, procedural methods, and an object-oriented blackboard concept. The system's objective is to assist engineers in lattice-dome design by integrating all design tasks into a single computer-aided environment with implementation of the knowledge-based system approach. For system verification, results from design examples are presented.

Costs increase as ritualized fighting progresses within and between phases in the sierra dome spider, Neriene litigiosa

Magnitudes and patterns of energy expenditure in animal contests are seldom measured, but can be critical for predicting contest dynamics and understanding the evolution of ritualized fighting behaviour. In the sierra dome spider, males compete for sexual access to females and their webs. They show three distinct phases of fighting behaviour, escalating from ritualized noncontact display (phase 1) to cooperative wrestling (phase 2), and finally to unritualized, potentially fatal fighting (phase 3). Using CO2 respirometry, we estimated energetic costs of male-male combat in terms of mean and maximum metabolic rates and the rate of increase in energy expenditure. We also investigated the energetic consequences of age and body mass, and compared fighting metabolism to metabolism during courtship. All three phases involved mean energy expenditures well above resting metabolic rate (3.5 X, 7.4 X and 11.5 X). Both mean and maximum energy expenditure became substantially greater as fights escalated through successive phases. The rates of increase in energy use during phases 2 and 3 were much higher than in phase 1. In addition, age and body mass affected contest energetics. These results are consistent with a basic prediction of evolutionarily stable strategy contest models, that sequences of agonistic behaviours should be organized into phases of escalating energetic costs. Finally, higher energetic costs of escalated fighting compared to courtship provide a rationale for first-male sperm precedence in this spider species. (C) 2004 The Association for the Study of Animal Behaviour. Published by Elsevier Ltd. All rights reserved.

Factors affecting the design and construction of Lamella suspen-dome systems

The suspen-dome system is a new structural form that has become popular in the construction of long-span roof structures. These domes are very slender and lightweight, their configuration is complicated, and hence sequential consideration in the structural design is needed. This paper focuses on these considerations, which include the method for designing cable prestress force, a simplified analysis method, and the estimation of buckling capacity. Buckling is one of the most important problems for dome structures. This paper presents the findings of an intensive buckling study of the Lamella suspen-dome system that takes geometric imperfection, asymmetric loading, rise-to-span ratio, and connection rigidity into consideration. Finally, suggested design and construction guidelines are given in the conclusion of this paper. (c) 2005 Elsevier Ltd. All rights reserved.

Microchemistry and microstructures of hydrothermally altered shock-metamorphosed basement gneiss, Woodleigh impact structure, Southern Carnarvon Basin, Western Australia

Hydrothermally altered shock-metamorphosed gneisses consisting of relic igneous biotite-K-feldspor-Na-rich alkali feldspar - plagioclase - quartz assemblages ( accessory garnet, corundum, titanite, monazite, zircon), and showing extensive replacement by montmorillonite, illite, sericite, and to a lesser extent chlorite, calcite, epidote, zoisite and pyrite, occur in the basement core uplift of the Woodleigh impact structure, Western Australia. The rocks display extensive hydrothermal clay alteration, complicating identification of pre-hydrothermal and pre-impact textures and compositions. Analysis of quartz-hosted planar deformation features (PDFs) indicates a majority of indexed sets parallel to omega{10 (1) over bar3}, a lesser abundance of sets parallel to pi{10 (1) over bar2}, and some sets parallel to the basal plane (0001) and r,z {10 (1) over bar1}, consistent with pressures about or over 20 GPa. Feldspar-hosted FDFs form reticulate vein networks displaying checkerboard-like to irregular and serrated patterns attributable to preferential replacement of shock-damaged PDFs and/or perthitic twin lamella by clay minerals. The gneisses are pervaded by clay-dominated intergranular and intragranular veins of cryptocrystalline material that display marked departures from bulk-rock chemistry and from mineral compositions. XRD analysis identifies the cryptocrystalline components as illite - montmorillonite, illite and chlorite, while laser Raman analysis identifies high-fluorescence sub-micrometre clay assemblage, feldspar, quartz and minor mica. SEM/EDS-probe and laser-ICPMS analysis indicate low-K high-Mg clay mineral compositions consistent with montmorillonite. Quartz PDF-hosted cryptocrystalline laminae display distinct enrichments in Al, Mg, Ca and K. Altered intergranular veins and feldspar-hosted cryptocrystalline components show consistent enrichment in the relatively refractory elements (Al, Cc, Mg, Fe) and depletion in relatively volatile elements (Si, K, Na). The clay alteration retards determination whether clay-dominated vein networks represent altered shock-induced pseudotachylite veins, diaplectic zones and/or shock-damaged twin lamella, and/or result from purely mineralogical and chemical differentiation affected by hydrothermal fluids, Overall enrichment of the shocked gneiss and of the cryptocrystalline components in Mg and trace ferromagnesian elements (Ni, Cc, Cr) may be attributed alternatively to introduction of siderophile element-rich fluid from the projectile, or/and contamination of hydrothermal fluids by MgO from dolomites surrounding the basement uplift. High Ni/Co and Ni/Cr and anomalous DGE (platinum group elements) may support the former model.

A combined spectrophotometric and morphometric study of the lunar mare dome fields near Cauchy, Arago, Hortensius, and Milichius

In this Study we examine the spectral and morphometric properties of the four important lunar mare dome fields near Cauchy, Arago, Hortensius. and Milichius. We utilize Clementine UV vis mulfispectral data to examine the soil composition of the mare domes while employing telescopic CCD imagery to compute digital elevation maps in order to determine their morphometric properties, especially flank slope, height, and edifice Volume. After reviewing previous attempts to determine topographic data for lunar domes, we propose an image-based 3D reconstruction approach which is based on a combination of photoclinometry and shape from shading. Accordingly, we devise a classification scheme for lunar Marc domes which is based on a principal component analysis of the determined spectral and morphometric features. For the effusive mare domes of the examined fields we establish four Classes, two of which are further divided into two subclasses, respectively, where each class represents distinct combinations of spectral and morphometric dome properties. As a general trend, shallow and steep domes formed out of low-TiO2 basalts are observed in the Hortensius and Milichius dome fields, while the domes near Cauchy and Arago that consist of high-TiO2 basalts are all very shallow. The intrusive domes of our data set cover a wide continuous range of spectral and morphometric quantities, generally characterized by larger diameters and shallower flank slopes than effusive domes. A comparison to effusive and intrusive mare domes in other lunar regions, highland domes, and lunar cones has shown that the examined four mare dome fields display Such a richness in spectral properties and 3D dome shape that the established representation remains valid in a more global context. Furthermore, we estimate the physical parameters of dome formation for the examined domes based on a rheologic model. Each class of effusive domes defined in terms of spectral and morphometric properties is characterized by its specific range of values for lava viscosity, effusion rate, and duration of the effusion process. For our data set we report lava viscosities between about 10(2) and 10(8) Pas, effusion rates between 25 and 600 m(3) s(-1), and durations of the effusion process between three weeks and 18 years. Lava viscosity decreases with increasing R-415/R-750 spectral ratio and thus TiO2 content; however, the correlation is not strong, implying an important influence of further parameters like effusion temperature on lava viscosity.

Adenocarcinoma of colon differentiating as dome epithelium of gut-associated lymphoid tissue

Aims: An early adenocarcinoma of the ascending colon was confined to a mass of gut-associated lymphoid tissue (GALT). The first description of an adenocarcinoma of colon differentiating as dome epithelium is presented. Methods and results: A plaque-like carcinoma was identified opposite the ileocaecal valve in an asymptomatic 56-year-old man with a family history of colorectal cancer. Malignant epithelium was confined to a mass of GALT filling but limited to the submucosa, Characterization of the neoplasm was undertaken by means of mucin histochemistry, immunohistochemistry, electron microscopy and assessment of DNA microsatellite instability status. The malignant epithelium comprised well differentiated columnar cells with a microvillous brush border and expressing MUC1, but no goblet cells or expression of MUC2. The demonstration of focal clusters of intraepithelial B-lymphocytes supported the presence of functioning M-cells within the malignant neoplasm. The cancer was DNA microsatellite stable despite the finding of tumour infiltrating lymphocytes. Conclusions: There is evidence for the origin of colorectal neoplasia from dome epithelium in both experimental models and microreconstruction studies of early adenomas in nonpolypotic human colorectal mucose, It is suggested that the lymphocyte-rich subset of colorectal cancer that expresses MUC1 but not MUC2 may be differentiating as dome epithelium of gut-associated lymphoid tissue.

On the overabundance of light rare earth elements in terrestrial zircons and its implication for Earth's earliest magmatic differentiation

We present whole-rock and zircon rare earth element (REE) data from two early Archaean gneisses (3.81 Ga and 3.64 Ga) from the Itsaq gneiss complex, south-west Greenland. Both gneisses represent extremely rare examples of unaltered, fresh and relatively undeformed igneous rocks of such antiquity. Cathodoluminescence imaging of their zircons indicates a single crystallisation episode with no evidence for either later metamorphic and/or anatectic reworking or inheritance of earlier grains. Uniform, single-population U/Pb age data confirm the structural simplicity of these zircons. One sample, a 3.64 Ga granodioritic gneiss from the Gothabsfjord, yields a chondrite-normalised REE pattern with a positive slope from La to Lu as well as substantial positive Ce and slight negative Eu anomalies, features generally considered to be typical of igneous zircon. In contrast, the second sample, a 3.81 Ga tonalite from south of the Isua Greenstone Belt, has variable but generally much higher light REE abundances, with similar middle to heavy REE. Calculation of zircon/melt distribution coefficients (D-REE(zircon/melt)) from each sample yields markedly different values for the trivalent REE (i.e. Ce and Eu omitted) and simple application of one set of D-REE(zircon/melt) to model the melt composition for the other sample yields concentrations that are in error by up to two orders of magnitude for the light REE (La-Nd). The observed light REE overabundance in the 3.81 Ga tonalite is a commonly observed feature in terrestrial zircons for which a number of explanations ranging from lattice strain to disequilibrium crystallisation have been proposed and are further investigated herein. Regardless of the cause of light REE overabundance, our study shows that simple application of zircon/melt distribution coefficients is not an unambiguous method for ascertaining original melt composition. In this context, recent studies that use REE data to claim that > 4.3 Ga Hadean detrital zircons originally crystallised from an evolved magma, in turn suggesting the operation of geological processes in the early Earth analogous to those of the present day (e.g. subduction and melting of hydrated oceanic crust), must be regarded with caution. Indeed, comparison of terrestrial Hadean and > 3.9 Ga lunar highland zircons shows remarkable similarities in the light REE, even though subduction processes that have been used to explain the terrestrial zircons have never operated on the Moon. (C) 2002 Elsevier Science B.V. All rights reserved.

A refined solution to Earth's hidden niobium: implications for evolution of continental crust and mode of core formation

New high-precision niobium (Nb) and tantalum (Ta) concentration data are presented for early Archaean metabasalts, metabasaltic komatiites and their erosion products (mafic metapelites) from SW Greenland and the Acasta gneiss complex, Canada. Individual datasets consistently show sub-chondritic Nb/Ta ratios averaging 15.1+/-11.6. This finding is discussed with regard to two competing models for the solution of the Nb-deficit that characterises the accessible Earth. Firstly, we test whether Nb could have sequestered into the core due to its slightly siderophile (or chalcophile) character under very reducing conditions, as recently proposed from experimental evidence. We demonstrate that troilite inclusions of the Canyon Diablo iron meteorite have Nb and V concentrations in excess of typical chondrites but that the metal phase of the Grant, Toluca and Canyon Diablo iron meteorites do not have significant concentrations of these lithophile elements. We find that if the entire accessible Earth Nb-deficit were explained by Nb in the core, only ca. 17% of the mantle could be depleted and that by 3.7 Ga, continental crust would have already achieved ca. 50% of its present mass. Nb/Ta systematics of late Archaean metabasalts compiled from the literature would further require that by 2.5 Ga, 90% of the present mass of continental crust was already in existence. As an alternative to this explanation, we propose that the average Nb/Ta ratio (15.1+/-11.6) of Earth's oldest mafic rocks is a valid approximation for bulk silicate Earth. This would require that ca. 13% of the terrestrial Nb resided in the Ta-free core. Since the partitioning of Nb between silicate and metal melts depends largely on oxygen fugacity and pressure, this finding could mean that metal/silicate segregation did not occur at the base of a deep magma ocean or that the early mantle was slightly less reducing than generally assumed. A bulk silicate Earth Nb/Ta ratio of 15.1 allows for depletion of up to 40% of the total mantle. This could indicate that in addition to the upper mantle, a portion of the lower mantle is depleted also, or if only the upper mantle were depleted, an additional hidden high Nb/Ta reservoir must exist. Comparison of Nb/Ta systematics between early and late Archaean metabasalts supports the latter idea and indicates deeply subducted high Nb/Ta eclogite slabs could reside in the mantle transition zone or the lower mantle. Accumulation of such slabs appears to have commenced between 2.5 and 2.0 Ga. Regardless of these complexities of terrestrial Nb/Ta systematics, it is shown that the depleted mantle Nb/Th ratio is a very robust proxy for the amount of extracted continental crust, because the temporal evolution of this ratio is dominated by Th-loss to the continents and not Nb-retention in the mantle. We present a new parameterisation of the continental crust volume versus age curve that specifically explores the possibility of lithophile element loss to the core and storage of eclogite slabs in the transition zone. (C) 2003 Elsevier Science B.V. All rights reserved.

Inheritance of early Archaean Pb-isotope variability from long-lived Hadean protocrust

Comparison of initial Pb-isotope signatures of several early Archaean (3.65-3.82 Ga) lithologies (orthogneisses and metasediments) and minerals (feldspar and galena) documents the existence of substantial isotopic heterogeneity in the early Archaean, particularly in the Pb-207/Pb-204 ratio. The magnitude of isotopic variability at 3.82-3.65 Ga requires source separation between 4.3 and 4.1 Ga, depending on the extent of U/Pb fractionation possible in the early Earth. The isotopic heterogeneity could reflect the coexistence of enriched and depleted mantle domains or the separation of a terrestrial protocrust with a U-238/Pb-204 (mu) that was ca. 20-30% higher than coeval mantle. We prefer this latter explanation because the high-p signature is most evident in metasediments (that formed at the Earth's surface). This interpretation is strengthened by the fact that no straightforward mantle model can be constructed for these high-mu lithologies without violating bulk silicate Earth constraints. The Pb-isotope evidence for a long-lived protocrust complements similar Hf-isotope data from the Earth's oldest zircons, which also require an origin from an enriched (low Lu/Hf) environment. A model is developed in which greater than or equal to3.8-Ga tonalite and monzodiorite gneiss precursors (for one of which we provide zircon U-Pb data) are not mantle-derived but formed by remelting or differentiation of ancient (ca. 4.3 Ga) basaltic crust which had evolved with a higher U/Pb ratio than coeval mantle in the absence of the subduction process. With the initiation of terrestrial subduction at, we propose, ca. 3.75 Ga, most of the greater than or equal to3.8-Ga basaltic shell (and its differentiation products) was recycled into the mantle, because of the lack of a stabilising mantle lithosphere. We argue that the key event for preservation of all greater than or equal to3.8-Ga terrestrial crust was the intrusion of voluminous granitoids immediately after establishment of global subduction because of complementary creation of a lithospheric keel. Furthermore, we argue that preservation of !3.8-Ga material (in situ rocks and zircons) globally is restricted to cratons with a high U/Pb source character (North Atlantic, Slave, Zimbabwe, Yilgarn, and Wyoming), and that the Pb-isotope systematics of these provinces are ultimately explained by reworking of material that was derived from ca. 4.3 Ga (i.e. Hadean) basaltic crust.