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.

Incorporating level set methods in Geographical Information Systems (GIS) for land-surface process modelling

Land-surface processes include a broad class of models that operate at a landscape scale. Current modelling approaches tend to be specialised towards one type of process, yet it is the interaction of processes that is increasing seen as important to obtain a more integrated approach to land management. This paper presents a technique and a tool that may be applied generically to landscape processes. The technique tracks moving interfaces across landscapes for processes such as water flow, biochemical diffusion, and plant dispersal. Its theoretical development applies a Lagrangian approach to motion over a Eulerian grid space by tracking quantities across a landscape as an evolving front. An algorithm for this technique, called level set method, is implemented in a geographical information system (GIS). It fits with a field data model in GIS and is implemented as operators in map algebra. The paper describes an implementation of the level set methods in a map algebra programming language, called MapScript, and gives example program scripts for applications in ecology and hydrology.

Human Fatalities from Cyanobacteria: Chemical and Biological Evidence for Cyanotoxins

An outbreak of acute liver failure occurred at a dialysis center in Caruaru, Brazil (8 degrees 17 'S, 35 degrees 58 'W), 134 km from Recife, the state capital of Pernambuco. At the clinic, 116 (89%) of 131 patients experienced visual disturbances, nausea, and vomiting after routine hemodialysis treatment on 13-20 February 1996. Subsequently, 100 patients developed acute liver failure, and of these 76 died. As of December 1996, 52 of the deaths could be attributed to a common syndrome now called Caruaru syndrome. Examination of phytoplankton from the dialysis clinic's water source, analyses of the clinic's water treatment system, plus serum and liver tissue of clinic patients led to the identification of two groups of cyanobacterial toxins, the hepatotoxic cyclic peptide microcystins and the hepatotoxic alkaloid cylindrospermopsin. Comparison of victims' symptoms and pathology using animal studies of these two cyanotoxins leads us to conclude that the major contributing factor to death of the dialyses patients was intravenous exposure to microcystins, specifically microcystin-YR, -LR, and -AR. From liver concentrations and exposure volumes, it was estimated that 19.5 mug/L microcystin was in the water used for dialysis treatments. This is 19.5 times the level set as a guideline for safe drinking water supplies by the World. Health Organization.

Brief note on the variation of constants formula for fuzzy differential equations

This note gives a theory of state transition matrices for linear systems of fuzzy differential equations. This is used to give a fuzzy version of the classical variation of constants formula. A simple example of a time-independent control system is used to illustrate the methods. While similar problems to the crisp case arise for time-dependent systems, in time-independent cases the calculations are elementary solutions of eigenvalue-eigenvector problems. In particular, for nonnegative or nonpositive matrices, the problems at each level set, can easily be solved in MATLAB to give the level sets of the fuzzy solution. (C) 2002 Elsevier Science B.V. All rights reserved.

Kidney dysfunction and hypertension: Role for cadmium, P450 and heme oxygenases?

Cadmium (Cd) is a metal toxin of continuing worldwide concern. Daily intake of Cd, albeit in small quantities, is associated with a number of adverse health effects which are attributable to distinct pathological changes in a variety of tissues and organs. In the present review, we focus on its renal tubular effects in people who have been exposed environmentally to Cd at levels below the provisional tolerable intake level set for the toxin. We highlight the data linking such low-level Cd intake with tubular injury, altered abundance of cytochromes P450 (CYPs) in the kidney and an expression of a hypertensive phenotype. We provide updated knowledge on renal and vascular effects of the eicosanoids 20-hydroxyeicosatetraenoic acid (20-HETE) and eicosatrienoic acids (EETs), which are biologically active metabolites from arachidonate metabolism mediated by certain CYPs in the kidney. We note the ability of Cd to elicit oxidative stress and to alter metal homeostasis notably of zinc which may lead to augmentation of the defense mechanisms involving induction of the antioxidant enzyme heme oxygenase-1 (HO-1) and the metal binding protein metallothionein (MT) in the kidney. We hypothesize that renal Cd accumulation triggers the host responses mediated by HO-I and MT in an attempt to protect the kidney against injurious oxidative stress and to resist a rise in blood pressure levels. This hypothesis predicts that individuals with less active HO-1 (caused by the HO-1 genetic polymorphisms) are more likely to have renal injury and express a hypertensive phenotype following chronic ingestion of low-level Cd, compared with those having more active HO-1. Future analytical and molecular epidemiologic research should pave the way to the utility of induction of heme oxygenases together with dietary antioxidants in reducing the risk of kidney injury and hypertension in susceptible people.

Two-level atom in a squeezed vacuum with finite bandwidth

The resonance fluorescence of a two-level atom driven by a coherent laser field and damped by a finite bandwidth squeezed vacuum is analysed. We extend the Yeoman and Barnett technique to a non-zero detuning of the driving field from the atomic resonance and discuss the role of squeezing bandwidth and the detuning in the level shifts, widths and intensities of the spectral lines. The approach is valid for arbitrary values of the Rabi frequency and detuning but for the squeezing bandwidths larger than the natural linewidth in order to satisfy the Markoff approximation. The narrowing of the spectral lines is interpreted in terms of the quadrature-noise spectrum. We find that, depending on the Rabi frequency, detuning and the squeezing phase, different factors contribute to the line narrowing. For a strong resonant driving field there is no squeezing in the emitted field and the fluorescence spectrum exactly reveals the noise spectrum. In this case the narrowing of the spectral lines arises from the noise reduction in the input squeezed vacuum. For a weak or detuned driving field the fluorescence exhibits a large squeezing and, as a consequence, the spectral lines have narrowed linewidths. Moreover, the fluorescence spectrum can be asymmetric about the central frequency despite the symmetrical distribution of the noise. The asymmetry arises from the absorption of photons by the squeezed vacuum which reduces the spontaneous emission. For an appropriate choice of the detuning some of the spectral lines can vanish despite that there is no population trapping. Again this process can be interpreted as arising from the absorption of photons by the squeezed vacuum. When the absorption is large it may compensate the spontaneous emission resulting in the vanishing of the fluorescence lines.

Field quantization, photons and non-Hermitean modes

Field quantization in unstable optical systems is treated by expanding the vector potential in terms of non-Hermitean (Fox-Li) modes. We define non-Hermitean modes and their adjoints in both the cavity and external regions and make use of the important bi-orthogonality relationships that exist within each mode set. We employ a standard canonical quantization procedure involving the introduction of generalized coordinates and momenta for the electromagnetic (EM) field. Three-dimensional systems are treated, making use of the paraxial and monochromaticity approximations for the cavity non-Hermitean modes. We show that the quantum EM field is equivalent to a set of quantum harmonic oscillators (QHOs), associated with either the cavity or the external region non-Hermitean modes, and thus confirming the validity of the photon model in unstable optical systems. Unlike in the conventional (Hermitean mode) case, the annihilation and creation operators we define for each QHO are not Hermitean adjoints. It is shown that the quantum Hamiltonian for the EM field is the sum of non-commuting cavity and external region contributions, each of which can be expressed as a sum of independent QHO Hamiltonians for each non-Hermitean mode, except that the external field Hamiltonian also includes a coupling term responsible for external non-Hermitean mode photon exchange processes. The non-commutativity of certain cavity and external region annihilation and creation operators is associated with cavity energy gain and loss processes, and may be described in terms of surface integrals involving cavity and external region non-Hermitean mode functions on the cavity-external region boundary. Using the essential states approach and the rotating wave approximation, our results are applied to the spontaneous decay of a two-level atom inside an unstable cavity. We find that atomic transitions leading to cavity non-Hermitean mode photon absorption are associated with a different coupling constant to that for transitions leading to photon emission, a feature consequent on the use of non-Hermitean mode functions. We show that under certain conditions the spontaneous decay rate is enhanced by the Petermann factor.

Development of a mungbean (Vigna radiata) RFLP linkage map and its comparison with lablab (Lablab purpureus) revelas a high level of colinearity between the two genomes

A genetic linkage map of mungbean (Vigna radiata, 2n = 2x = 22) consisting of 255 RFLP loci was developed using a recombinant inbred population of 80 individuals. The population was derived from an intersubspecific cross between the cultivated mungbean variety 'Berken' and a wild mungbean genotype 'ACC 41' (V radiata subsp. sublobata). The total length of the map, which comprised 13 linkage groups, spanned 737.9 cM with an average distance between markers of 3.0 cM and a maximum distance between linked markers of 15.4 cM. The mungbean map was compared to a previously published map of lablab (Lablab purpureus, 2n = 2x = 24) using a common set of 65 RFLP probes. In contrast to some other comparative mapping studies among members of the Fabaceae, where a high level of chromosomal rearrangement has been observed, marker order between mungbean and lablab was found to be highly conserved. However, the two genomes have apparently accumulated a large number of duplications/deletions after they diverged.

Exploration of the cell-cycle genes found within the RIKEN FANTOM2 data set

The cell cycle is one of the most fundamental processes within a cell. Phase-dependent expression and cell-cycle checkpoints require a high level of control. A large number of genes with varying functions and modes of action are responsible for this biology. In a targeted exploration of the FANTOM2-Variable Protein Set, a number of mouse homologs to known cell-cycle regulators as well as novel members of cell-cycle families were identified. Focusing on two prototype cell-cycle families, the cyclins and the NIMA-related kinases (NEKs), we believe we have identified all of the mouse members of these families, 24 cyclins and 10 NEKs, and mapped them to ENSEMBL transcripts. To attempt to globally identify all potential cell cycle-related genes within mouse, the MGI (Mouse Genome Database) assignments for the RIKEN Representative Set (RPS) and the results from two homology-based queries were merged. We identified 1415 genes with possible cell-cycle roles, and 1758 potential paralogs. We comment on the genes identified in this screen and evaluate the merits of each approach.

Higher-level language deficits resulting from left primary cerebellar lesions

Background: Contemporary neuropsychological studies suggest that cerebellar lesions may impact upon higher-level cognitive functioning via mechanisms of crossed cerebello-cerebral diaschisis. Accordingly, right cerebellar lesions have been previously associated with linguistic impairments such as reduced word fluency and agrammatic output. Recently, however, neuroimaging investigations have also identified ipsilateral cerebral hypoperfusion as a consequence of cerebellar lesions, implicating a potential role for the left cerebellum in the mediation of language processes. Aims: The purpose of this research was to investigate the effects of left cerebellar lesions of vascular origin, on general as well as high-level language skills. Methods & Procedures: Linguistic profiles were compiled for five individuals with left primary cerebellar lesions utilising a comprehensive language test battery. Individual scores relevant to each subtest were compared to a group of non-neurologically impaired controls. The criterion for anomalous performance was established as greater than or equal to 1.5 SD below the mean of the control group. Outcomes & Results: The findings of this research suggest that higher-level language deficits may result from left primary cerebellar lesions. All participants demonstrated deficits on measures of word fluency, sentence construction within a set context, producing word definitions, and producing multiple definitions for the same word. Deficits were also noted for several participants on measures of understanding figurative language, forming word associations, identifying and correcting semantic absurdities, and producing synonyms and antonyms. Conclusions: The results presented challenge the notion of a lateralised linguistic cerebellum, supporting a potential role for the left as well as right cerebellar hemispheres in the regulation of language processes, presumably via cerebellar-basal ganglia/thalamo-cortical pathways.

Adverse Health Effects of Chronic Exposure to Low-Level Cadmium in Foodstuffs and Cigarette Smoke

Cadmium is a cumulative nephrotoxicant that is absorbed into the body from dietary sources and cigarette smoking. The levels of Cd in organs such as liver and kidney cortex increase with age because of the lack of an active biochemical process for its elimination coupled with renal reabsorption. Recent research has provided evidence linking Cd-related kidney dysfunction and decreases in bone mineral density in nonoccupationally exposed populations who showed no signs of nutritional deficiency. This challenges the previous view that the concurrent kidney and bone damage seen in Japanese itai-itai disease patients was the result of Cd toxicity in combination with nutritional deficiencies, notably, of zinc and calcium. Further, such Cd-linked bone and kidney toxicities were observed in people whose dietary Cd intakes were well within the provisional tolerable weekly intake (PTWI) set by the Joint Food and Agriculture Organization/World Health Organization Expert Committee on Food Additives of 1 mug/kg body weight/day or 70 mug/day. This evidence points to the much-needed revision of the current PTWI for Cd. Also, evidence for the carcinogenic risk of chronic Cd exposure is accumulating and Cd effects on reproductive outcomes have begun to emerge.

Phobic anxiety in 11 nations: part II. Hofstede's dimensions of national cultures predict national-level variations

Hofstede's dimensions of national cultures termed Masculinity-Femininity (MAS) and Uncertainty Avoidance (UAI) (Hofstede, 2001) are proposed to be of relevance for understanding national-level differences in self-assessed fears. The potential predictive role of national MAS was based on the classical work of Fodor (Fodor, 1974). Following Fodor, it was predicted that masculine (or tough) societies in which clearer differentiations are made between gender roles (high MAS) would report higher national levels of fears than feminine (or soft/modest) societies in which such differentiations are made to a clearly lesser extent (low MAS). In addition, it was anticipated that nervous-stressful-emotionally-expressive nations (high UAI) would report higher national levels of fears than calm-happy and low-emotional countries (low UAI), and that countries high on both MAS and UAI would report the highest national levels of fears. A data set comprising 11 countries (N > 5000) served as the basis for analyses. As anticipated, (a) high MAS predicted higher national levels of Agoraphobic fears and of Bodily Injury-Illness-Death fears; (b) higher scores on both UAI and MAS predicted higher national scores on Bodily Injury-Illness-Death fears, fears of Sexual and Aggressive Scenes, and Harmless Animals fears; (c) higher UAI predicted higher national levels of Harmless Animals, Bodily Injury-Illness-Death, and Agoraphobic fears. (C) 2003 Elsevier Ltd. All rights reserved.

Theory of strongly correlated electron systems. I. Intersite coulomb interaction and the approximation of renormalized fermions in total energy calculations

The diagrammatic strong-coupling perturbation theory (SCPT) for correlated electron systems is developed for intersite Coulomb interaction and for a nonorthogonal basis set. The construction is based on iterations of exact closed equations for many - electron Green functions (GFs) for Hubbard operators in terms of functional derivatives with respect to external sources. The graphs, which do not contain the contributions from the fluctuations of the local population numbers of the ion states, play a special role: a one-to-one correspondence is found between the subset of such graphs for the many - electron GFs and the complete set of Feynman graphs of weak-coupling perturbation theory (WCPT) for single-electron GFs. This fact is used for formulation of the approximation of renormalized Fermions (ARF) in which the many-electron quasi-particles behave analogously to normal Fermions. Then, by analyzing: (a) Sham's equation, which connects the self-energy and the exchange- correlation potential in density functional theory (DFT); and (b) the Galitskii and Migdal expressions for the total energy, written within WCPT and within ARF SCPT, a way we suggest a method to improve the description of the systems with correlated electrons within the local density approximation (LDA) to DFT. The formulation, in terms of renormalized Fermions LIDA (RF LDA), is obtained by introducing the spectral weights of the many electron GFs into the definitions of the charge density, the overlap matrices, effective mixing and hopping matrix elements, into existing electronic structure codes, whereas the weights themselves have to be found from an additional set of equations. Compared with LDA+U and self-interaction correction (SIC) methods, RF LDA has the advantage of taking into account the transfer of spectral weights, and, when formulated in terms of GFs, also allows for consideration of excitations and nonzero temperature. Going beyond the ARF SCPT, as well as RF LIDA, and taking into account the fluctuations of ion population numbers would require writing completely new codes for ab initio calculations. The application of RF LDA for ab initio band structure calculations for rare earth metals is presented in part 11 of this study (this issue). (c) 2005 Wiley Periodicals, Inc.

Upper level, C House, Coorparoo, Brisbane

View along North-West elevation upper level.