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.

Power of the joint segregation analysis method for testing mixed major-gene and polygene inheritance models of quantitative traits

Understanding the genetic architecture of quantitative traits can greatly assist the design of strategies for their manipulation in plant-breeding programs. For a number of traits, genetic variation can be the result of segregation of a few major genes and many polygenes (minor genes). The joint segregation analysis (JSA) is a maximum-likelihood approach for fitting segregation models through the simultaneous use of phenotypic information from multiple generations. Our objective in this paper was to use computer simulation to quantify the power of the JSA method for testing the mixed-inheritance model for quantitative traits when it was applied to the six basic generations: both parents (P-1 and P-2), F-1, F-2, and both backcross generations (B-1 and B-2) derived from crossing the F-1 to each parent. A total of 1968 genetic model-experiment scenarios were considered in the simulation study to quantify the power of the method. Factors that interacted to influence the power of the JSA method to correctly detect genetic models were: (1) whether there were one or two major genes in combination with polygenes, (2) the heritability of the major genes and polygenes, (3) the level of dispersion of the major genes and polygenes between the two parents, and (4) the number of individuals examined in each generation (population size). The greatest levels of power were observed for the genetic models defined with simple inheritance; e.g., the power was greater than 90% for the one major gene model, regardless of the population size and major-gene heritability. Lower levels of power were observed for the genetic models with complex inheritance (major genes and polygenes), low heritability, small population sizes and a large dispersion of favourable genes among the two parents; e.g., the power was less than 5% for the two major-gene model with a heritability value of 0.3 and population sizes of 100 individuals. The JSA methodology was then applied to a previously studied sorghum data-set to investigate the genetic control of the putative drought resistance-trait osmotic adjustment in three crosses. The previous study concluded that there were two major genes segregating for osmotic adjustment in the three crosses. Application of the JSA method resulted in a change in the proposed genetic model. The presence of the two major genes was confirmed with the addition of an unspecified number of polygenes.

A new method for identification of protein (sub)families in a set of proteins based on hydropathy distribution in proteins

Structural similarity among proteins is reflected in the distribution of hydropathicity along the amino acids in the protein sequence. Similarities in the hydropathy distributions are obvious for homologous proteins within a protein family. They also were observed for proteins with related structures, even when sequence similarities were undetectable. Here we present a novel method that employs the hydropathy distribution in proteins for identification of (sub)families in a set of (homologous) proteins. We represent proteins as points in a generalized hydropathy space, represented by vectors of specifically defined features. The features are derived from hydropathy of the individual amino acids. Projection of this space onto principal axes reveals groups of proteins with related hydropathy distributions. The groups identified correspond well to families of structurally and functionally related proteins. We found that this method accurately identifies protein families in a set of proteins, or subfamilies in a set of homologous proteins. Our results show that protein families can be identified by the analysis of hydropathy distribution, without the need for sequence alignment. (C) 2005 Wiley-Liss, Inc.

A fast cross-entropy method for estimating buffer overflows in queuing networks

In this paper, we propose a fast adaptive importance sampling method for the efficient simulation of buffer overflow probabilities in queueing networks. The method comprises three stages. First, we estimate the minimum cross-entropy tilting parameter for a small buffer level; next, we use this as a starting value for the estimation of the optimal tilting parameter for the actual (large) buffer level. Finally, the tilting parameter just found is used to estimate the overflow probability of interest. We study various properties of the method in more detail for the M/M/1 queue and conjecture that similar properties also hold for quite general queueing networks. Numerical results support this conjecture and demonstrate the high efficiency of the proposed algorithm.

Quantum trajectory simulations of the fluorescence intensity from a two-level atom driven by a multichromatic field

The quantum trajectories method is illustrated for the resonance fluorescence of a two-level atom driven by a multichromatic field. We discuss the method for the time evolution of the fluorescence intensity in the presence of bichromatic and trichromatic driving fields. We consider the special case wherein one multichromatic field component is strong and resonant with the atomic transition whereas the other components are much weaker and arbitrarily detuned from the atomic resonance. We find that the phase-dependent modulations of the Rabi oscillations, recently observed experimentally [Q. Wu, D. J. Gauthier, and T. W. Mossberg, Phys. Rev. A 49, R1519 (1994)] for the special case when the weaker component of a bichromatic driving field is detuned from the atomic resonance by the strong-field Rabi frequency, appear also for detunings close to the subharmonics of the Rabi frequency. Furthermore, we show that for the atom initially prepared in one of the dressed states of the strong field component the modulations are not sensitive to the phase. We extend the calculations to the case of a trichromatic driving field and find that apart from the modulations of the amplitude there is a modulation of the frequency of the Rabi oscillations. Moreover, the time evolution of the fluorescence intensity depends on the phase regardless of the initial conditions and a phase-dependent suppression of the Rabi oscillations can be observed when the sideband fields are tuned to the subharmonics of the strong-field Rabi frequency. [S1050-2947(98)03501-X].

Dental set-vice provision in private general practice for an insured population in New South Wales, Australia

Objective To determine patterns of dental set-vices provided to a cohort of the insured population 18 years and over, in private general practice in New South Wales, Australia. Basic research design A cohort study using the person-years method and Poisson regression for analysis. Setting Data were derived from claims records submitted by members of a health insurance fund (Government Employees Health Fund-GEHF) for rebates during the study period 1 January 1992-31 December 1995. Participants There were 133,467 members aged 18 years and over from New South Wales. Main outcome measures To determine, by age group, for those members who used private general practice and made a claim (referred to as 'patients') the annual number of visits, total number of services received per year and number of services received at a visit, Results The mean number of visits per patient was 2.4 per year with patients under 45 years making fewer visits than the 45-54 age group reference category. Mean number of services utilised per patient-year was 5.9, with services provided increasing from 3.5 for the 18-24-year-old group, reaching a plateau of approximately 6.2 for those aged 45 years or more. The number of services received per visit was 2.4 and there were no differences by age. Service mix was dominated by restorative (35%), diagnostic (27%), and preventive services (18%); with age specific variations observed. Conclusions Age was found to be an important determinant in the use of dental services, independent of membership duration and gender.

A method for setting the size of plant conservation target areas

Realistic time frames in which management decisions are made often preclude the completion of the detailed analyses necessary for conservation planning. Under these circumstances, efficient alternatives may assist in approximating the results of more thorough studies that require extensive resources and time. We outline a set of concepts and formulas that may be used in lieu of detailed population viability analyses and habitat modeling exercises to estimate the protected areas required to provide desirable conservation outcomes for a suite of threatened plant species. We used expert judgment of parameters and assessment of a population size that results in a specified quasiextinction risk based on simple dynamic models The area required to support a population of this size is adjusted to take into account deterministic and stochastic human influences, including small-scale disturbance deterministic trends such as habitat loss, and changes in population density through processes such as predation and competition. We set targets for different disturbance regimes and geographic regions. We applied our methods to Banksia cuneata, Boronia keysii, and Parsonsia dorrigoensis, resulting in target areas for conservation of 1102, 733, and 1084 ha, respectively. These results provide guidance on target areas and priorities for conservation strategies.

The design of planar gradient coils. Part I: A winding path correction method

In this work a new approach for designing planar gradient coils is outlined for the use in an existing MRI apparatus. A technique that allows for gradient field corrections inside the diameter-sensitive volume is deliberated. These corrections are brought about by making changes to the wire paths that constitute the coil windings, and hence, is called the path correction method. The existing well-known target held method is used to gauge the performance of a typical gradient coil. The gradient coil design methodology is demonstrated for planar openable gradient coils that can be inserted into an existing MRI apparatus. The path corrected gradient coil is compared to the coil obtained using the target field method. It is shown that using a wire path correction with optimized variables, winding patterns that can deliver high magnetic gradient field strengths and large imaging regions can be obtained.

Local dispersion of the Eucalyptus leaf beetle Chrysophtharta bimaculata (Coleoptera: Chrysomelidae), and implications for forest protection

1. Chrysophtharta bimaculata is a native chrysomelid species that can cause chronic defoliation of plantation and regrowth Eucalyptus forests in Tasmania, Australia. Knowledge of the dispersion pattern of C. bimaculata was needed in order to assess the efficiency of an integrated pest management (IPM) programme currently used for its control. 2. Using data from yellow flight traps, local populations of C. bimaculata adults were monitored over a season at spatial scales relevant to commercial forestry: within a 50-ha operational management unit (a forestry 'coupe') and between coupes. In addition, oviposition was monitored over a season at a subset of the between-coupe sites. 3. Dispersion indices (Taylor's Power Law and Iwao's Mean Crowding regression method) demonstrated that C. bimaculata adults were spatially aggregated within and between coupes, although the number of egg-batches laid at the between-coupe scale was uniform. Spatial autocorrelation analysis showed that trap-catches at the within-coupe level were similar (positively autocorrelated) to a radius distance of approximately 110 m, and then dissimilar (negatively autocorrelated) at approximately 250 m. At the between-coupe scale, no repeatable spatial autocorrelation patterns were observed. 4. For any individual site, rapid changes in beetle density were observed to be associated with loosely aggregated flights of beetles into and out of that site. Peak adult catches (> the weekly mean plus standard deviation trap-catch) for a site occurred for a period of 2.0 +/- 0.22 weeks at a time (n = 37), with normally only one or two peaks per site per season. Peak oviposition events for a site occurred on average 1.4 +/- 0.11 times per season and lasted 1.5 +/- 0.12 weeks. 5. Analysis of an extensive data set (n = 417) demonstrated that adult abundance at a site was positively correlated with egg density, but negatively correlated with tree damage (caused by conspecifics) and the presence of conspecific larvae. There was no relationship between adult abundance and a visual estimate of the amount of young foliage on trees. 6. Adults of C. bimaculata are show n to occur in relatively small, mobile aggregations. This means that pest surveys must be both regular (less than 2 weeks apart) and intensive (with sampling points no more than 150 m apart) if beetle populations are to be monitored with confidence. Further refinement of the current IPM strategy must recognize the problems posed by this temporal and spatial patchiness, particularly with regard to the use of biological insecticides, such as Bacillus thuringiensis, for which only a very short operational window exists.

Changes in conceptions of learning for indigenous Australian university students

Background. Conceptions of learning have been investigated for students in higher. education in different countries. Some studies found that students' conceptions change and develop over time while others have found no changes. Investigating conceptions of learning for Australian Aboriginal and Torres Strait Islander university students is a relatively new area of research. Aims. This study set out to investigate conceptions of learning for Aboriginal and Torres Strait Islander university students during the first two years of their undergraduate degree courses in three Australian universities. Conceptions for each year were compared. Knowing, more about learning as conceived by this cultural group may facilitate more productive higher educational experiences. Sample. The sample comprised 17 students studying various degrees; Il were male and 6 were female. Ages ranged from 18 to 48 years; mean age was 26 years. Method. This was a phenomenographic, longitudinal study. Individual semistructured interviews were conducted each year to ascertain students' conceptions of learning. Conceptions for second year were derived independently of those From first year. A comparative analysis then took place to determine ally changes. Results. These students held conceptions of learning that were similar to those of other university students; however there were some intrinsic differences. On a group level, conceptions changed somewhat over the two years as did core conceptions reported by some individual students. Some students also exhibited a greater awareness of learning during their second year that resulted in three dimensions of changed awareness. Conclusions. We believe the changed conceptions and awareness resulted from learning at university where there is some need to understand and explain phenomena in relation to theory. This brought about new understandings which allowed students to see their own learning in a relational sense.