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.

Demonic Possession and Exorcism in Early Modern France

This is a highly original study of possession by demons and their exorcism, which was rife in early modern times, focussing on the place where they were most prevalent, France. Catholics at the time believed that the devil was present everywhere, in the rise of the heretics, in the activities of witches, and even in the bodies of the pious young women. The rite of exorcism was intended to heal the possessed and show the power of the church - but it generated as many problems as it resolved. Possessed nuns endured frequently violent exorcisms, exorcists were suspected of conjuring devils, and possession itself came to be seen as a form of holiness, elevating several women to the status of living saints. Sarah Ferber offers a challenging study of one of the most intriguing phenomena of early modern Europe; looking also at the present day, it argues that early modern conflicts over the devil still carry an unexpected force and significance for western Christianity (from publisher).

Hysteresis and phase transitions in a single partially internally wetted catalyst pellet: Thermogravimetric studies

The hydrogenation of cyclohexene over palladium supported in a microporous gamma-alumina pellet is studied thermogravimetrically with a view to measuring the extent of partial internal wetting associated with the different steady state branches. As many as three steady state branches having significantly different degrees of internal wetting and reaction rates, with transitions between them, are confirmed from observations of catalyst weight change. It is seen that with reduction in catalyst activity the middle branch, obtained by condensation from a vapor filled pellet, is much more prominent without showing an evaporative transition for the range of hydrogen partial pressures used here. The catalyst activity is therefore an important parameter affecting the structure of the steady state branches. Hysteresis effects are found to occur, and the thermogravimetric results also confirm the importance of history in determining the catalyst state. The measured degree of wetting is in accordance with that estimated from a mathematical model incorporating capillary condensation effects in addition to reaction-diffusion phenomena. The same model also satisfactorily interprets the reaction rate variations and transitions seen in the present work.

A comparison of lipopolysaccharide-induced hepatitis in ethanol-fed Wistar and Lewis rats

Elevated concentrations of plasma proinflammatory cytokines have been detected in patients with alcoholic hepatitis (AH) and in a model of lipopolysaccharide-induced hepatitis in ethanol-fed Wistar rats. These cytokines have been implicated in the pathogenesis of the liver damage. Considering the likely involvement of the immune system in AH, and the frequent use of Lewis rats in autoimmune disease models, Lewis rats were examined in the model to determine whether they would more closely mimic the immune status of a chronic alcoholic and be a preferable strain for use in future experiments. Lipopolysaccharide-induced hepatic tumor necrosis factor-cu, interleukin-1 alpha, interleukin-1 beta, and interleukin-6 mRNA expression was examined in both rat strains. The overall pattern of histological (panlobular piecemeal necrosis) and biochemical liver damage (plasma ALT levels), and cytokine expression was similar in both strains. Thus, it would appear that, despite the known susceptibility of Lewis rats to autoimmune phenomena, they do not respond to the experimental regime significantly better than Wistar rats. This study confirms that unknown mediators are contributing to the liver damage seen in this model and possibly in AH.

A relative gradient theory for layered materials

It is possible to remedy certain difficulties with the description of short wave length phenomena and interfacial slip in standard models of a laminated material by considering the bending stiffness of the layers. If the couple or moment stresses are assumed to be proportional to the relative deformation gradient, then the bending effect disappears for vanishing interface slip, and the model correctly reduces to an isotropic standard continuum. In earlier Cosserat-type models this was not the case. Laminated materials of the kind considered here occur naturally as layered rock, or at a different scale, in synthetic layered materials and composites. Similarities to the situation in regular dislocation structures with couple stresses, also make these ideas relevant to single slip in crystalline materials. Application of the theory to a one-dimensional model for layered beams demonstrates agreement with exact results at the extremes of zero and infinite interface stiffness. Moreover, comparison with finite element calculations confirm the accuracy of the prediction for intermediate interfacial stiffness.

Exploring sequential data: Commentary on Bowers, Jentsch, Salas, and Braun (1998)

In an investigation intended to determine training needs of night crews, Bowers et al. (1998, this issue) report two studies showing that the patterning of communication is a better discriminator of good and poor crews than is the content of communication. Bowers et al. characterize their studies as intended to generate hypotheses for training needs and draw connections with Exploratory Sequential Data Analysis (ESDA). Although applauding the intentions of Bowers ct al., we point out some concerns with their characterization and implementation of ESDA. Our principal concern is that the Bowers et al. exploration of the data does not convincingly lead them back to a better fundamental understanding of the original phenomena they are investigating.

Temperature-induced bleaching of corals begins with impairment of the CO2 fixation mechanism in zooxanthellae

The early effects of heat stress on the photosynthesis of symbiotic dinoflagellates (zooxanthellae) within the tissues of a reef-building coral were examined using pulse-amplitude-modulated (PAM) chlorophyll fluorescence and photorespirometry. Exposure of Stylophora pistillata to 33 and 34 degrees C for 4 h resulted in (1) the development of strong non-photochemical quenching (qN) of the chlorophyll fluorescence signal, (2) marked decreases in photosynthetic oxygen evolution, and (3) decreases in optimal quantum yield (F-v/F-m) of photosystern II (PSII), Quantum yield decreased to a greater extent on the illuminated surfaces of coral branches than on lower (shaded) surfaces, and also when high irradiance intensities were combined with elevated temperature (33 degrees C as opposed to 28 degrees C), qN collapsed in heat-stressed samples when quenching analysis was conducted in the absence of oxygen, Collectively, these observations are interpreted as the initiation of photoprotective dissipation of excess absorbed energy as heat (qN) and O-2-dependent electron flow through the Mehler-Ascorbate-Peroxidase cycle (MAP-cycle) following the point at which the rate of light-driven electron transport exceeds the capacity of the Calvin cycle. A model for coral bleaching is proposed whereby the primary site of heat damage in S, pistillata is carboxylation within the Calvin cycle, as has been observed during heat damage in higher plants, Damage to PSII and a reduction in F-v/F-m (i.e. photoinhibition) are secondary effects following the overwhelming of photoprotective mechanisms by light. This secondary factor increases the effect of the primary variable, temperature. Potential restrictions of electron flow in heat-stressed zooxanthellae are discussed with respect to Calvin cycle enzymes and the unusual status of the dinoflagellate Rubisco, Significant features of our model are that (1) damage to PSII is not the initial step in the sequence of heat stress in zooxanthellae, acid (2) light plays a key secondary role in the initiation of the bleaching phenomena.

Quantum oscillations in quasi-one-dimensional metals with spin-density-wave ground states

We consider the magnetoresistance oscillation phenomena in the Bechgaard salts (TMTSF)(2)X, where X = ClO4, PF6, and AsF6 in pulsed magnetic fields to 51 T. Of particular importance is the observation of a new magnetoresistance oscillation for X = ClO4 in its quenched state. In the absence of any Fermi-surface reconstruction due to anion order at low temperatures, all three materials exhibit nonmonotonic temperature dependence of the oscillation amplitude in the spin-density-wave (SDW) state. We discuss a model where, below a characteristic temperature T* within the SDW state, a magnetic breakdown gap opens. [S0163-1829(99)00904-2].

Accuracy of predicting a child's response to potential threat: A comparison of children and their mothers

Considerable research has indicated that children and their parents often demonstrate marked discrepancies in their reporting of anxiety-related phenomena. In such cases, the question arises as to whether children are capable of accurately reporting on their anxiety. In the present study, 50 children (aged 5 to 14 years) were asked to approach a large, German Shepherd dog. Prior to the task, both the mother and child independently predicted the closest point likely to be reached by the child and the degree of anxiety likely to be experienced. These predictions were then compared with the actual phenomena displayed by the child during the task. On the behavioural measure (closest step reached), both the child and mother demonstrated equivalent predictive accuracy. On the subjective measure (fear ratings) children were considerably more accurate than their mothers. The data were not influenced by gender, age, or clinical status. The results indicate the ability of children to accurately predict their anxious responses, and support the value of incorporating children's self-reports in the assessment of emotional disorders.

Chemical synthesis and folding pathways of large cyclic polypeptide: Studies of the cystine knot polypeptide kalata B1

Kalata B1 is a member of a new family of polypeptides, isolated from. plants, which have a cystine knot structure embedded within an amide-cyclized backbone. This family of molecules are the largest known cyclic peptides, and thus, the mechanism of synthesis and folding is of great interest. To provide information about both these phenomena, we have synthesized kalata B1 using two distinct strategies. In the first, oxidation of the cysteine residues of a linear precursor peptide to form the correct disulfide bonds results in folding of the three-dimensional structure and preorganization of the termini in close proximity for subsequent cyclization. The second approach involved cyclization prior to oxidation. In the first method, the correctly folded peptide was produced only in the presence of partially hydrophobic solvent conditions. These conditions are presumably required to stabilize the surface-exposed hydrophobic residues. However,; in the synthesis,involving cyclization prior to oxidation, the cyclic reduced peptide folded to a significant degree in the absence of hydrophobic solvents and even more efficiently in the presence of hydrophobic solvents. Cyclization clearly has a major effect on the folding pathway and facilitates formation of the correctly disulfide-bonded form in aqueous solution; In addition to facilitating folding to a compact stable structure cyclization has an important effect on biological activity as assessed by hemolytic activity.

A method of analysing stream channel response to environmental change: gauge data for the Tully River

Gauging data are available from numerous streams throughout Australia, and these data provide a basis for historical analysis of geomorphic change in stream channels in response to both natural phenomena and human activities. We present a simple method for analysis of these data, and a briefcase study of an application to channel change in the Tully River, in the humid tropics of north Queensland. The analysis suggests that this channel has narrowed and deepened, rather than aggraded: channel aggradation was expected, given the intensification of land use in the catchment, upstream of the gauging station. Limitations of the method relate to the time periods over which stream gauging occurred; the spatial patterns of stream gauging sites; the quality and consistency of data collection; and the availability of concurrent land-use histories on which to base the interpretation of the channel changes.

Social identity and self-categorization processes in organizational contexts

Although aspects of social identity theory are familiar to organizational psychologists, its elaboration, through self-categorization theory, of how social categorization and prototype-based depersonalization actually produce social identity effects is less well known. We describe these processes, relate self-categorization theory to social identity theory, describe new theoretical developments in detail, and show how these developments can address a: range of organizational phenomena. We discuss cohesion and deviance, leadership, subgroup and sociodemographic structure, and mergers and acquisitions.

Three dimensional preparatory trunk motion precedes asymmetrical upper limb movement

Three-dimensional trunk motion. trunk muscle electromyography and intra-abdominal pressure were evaluated to investigate the preparatory control of the trunk associated with voluntary unilateral upper limb movement. The directions of angular motion produced by moments reactive to limb movement in each direction were predicted using a three-dimensional model of the body. Preparatory motion of the trunk occurred in three dimensions in the directions opposite to the reactive moments. Electromyographic recordings from the superficial trunk muscles were consistent with preparatory trunk motion. However, activation of transversus abdominis was inconsistent with control of direction-specific moments acting on the trunk. The results provide evidence that anticipatory postural adjustments result in movements and not simple rigidification of the trunk. (C) 2000 Elsevier Science B.V. All rights reserved.

Robust Mixture Modelling Using the t Distribution

Normal mixture models are being increasingly used to model the distributions of a wide variety of random phenomena and to cluster sets of continuous multivariate data. However, for a set of data containing a group or groups of observations with longer than normal tails or atypical observations, the use of normal components may unduly affect the fit of the mixture model. In this paper, we consider a more robust approach by modelling the data by a mixture of t distributions. The use of the ECM algorithm to fit this t mixture model is described and examples of its use are given in the context of clustering multivariate data in the presence of atypical observations in the form of background noise.

Large amplitude folding in finely layered viscoelastic rock structures

We analyze folding phenomena in finely layered viscoelastic rock. Fine is meant in the sense that the thickness of each layer is considerably smaller than characteristic structural dimensions. For this purpose we derive constitutive relations and apply a computational simulation scheme (a finite-element based particle advection scheme; see MORESI et al., 2001) suitable for problems involving very large deformations of layered viscous and viscoelastic rocks. An algorithm for the time integration of the governing equations as well as details of the finite-element implementation is also given. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered plate. The domain is compressed along the layer axis by prescribing velocities along the sides. First, for the viscous limit we consider the response to a series of harmonic perturbations of the director orientation. The Fourier spectra of the initial folding velocity are compared for different viscosity ratios. Turning to the nonlinear regime we analyze viscoelastic folding histories up to 40% shortening. The effect of layering manifests itself in that appreciable buckling instabilities are obtained at much lower viscosity ratios (1:10) as is required for the buckling of isotropic plates (1:500). The wavelength induced by the initial harmonic perturbation of the director orientation seems to be persistent. In the section of the parameter space considered here elasticity seems to delay or inhibit the occurrence of a second, larger wavelength. Finally, in a linear instability analysis we undertake a brief excursion into the potential role of couple stresses on the folding process. The linear instability analysis also provides insight into the expected modes of deformation at the onset of instability, and the different regimes of behavior one might expect to observe.