Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

A 2D numerical model for simulating the physics of fault systems

Simulations provide a powerful means to help gain the understanding of crustal fault system physics required to progress towards the goal of earthquake forecasting. Cellular Automata are efficient enough to probe system dynamics but their simplifications render interpretations questionable. In contrast, sophisticated elasto-dynamic models yield more convincing results but are too computationally demanding to explore phase space. To help bridge this gap, we develop a simple 2D elastodynamic model of parallel fault systems. The model is discretised onto a triangular lattice and faults are specified as split nodes along horizontal rows in the lattice. A simple numerical approach is presented for calculating the forces at medium and split nodes such that general nonlinear frictional constitutive relations can be modeled along faults. Single and multi-fault simulation examples are presented using a nonlinear frictional relation that is slip and slip-rate dependent in order to illustrate the model.

Closely related Symbiodinium spp. differ in relative dominance in coral reef host communities across environmental, latitudinal and biogeographic gradients

The diversity and community structures of symbiotic dinoflagellates are described from reef invertebrates in southern and central provinces of the Great Barrier Reef (GBR), Australia, and Zamami Island, Okinawa, Japan. The symbiont assemblages from region to region were dominated by Clade C Symbiodinium spp. and consisted of numerous host-specific and/or rare types (specialists), and several types common to many hosts (generalists). Prevalence in the host community among certain host-generalist symbionts differed between inshore and offshore environments, across latitudinal (central versus southern GBR) gradients, and over wide geographic ranges (GBR versus Okinawa). One particular symbiont (C3h) from the GBR had a dramatic shift in dominance. Its prevalence ranged from being extremely rare, or absent on high-latitude reefs to dominating the scleractinian diversity on a mid-latitude inshore reef. These changes occurred among coral fauna whose larvae must acquire symbionts from environmental sources (horizontal symbiont acquisition). Such differences did not occur among 'vertical transmitters' such as Porites spp., Montipora spp. and pocilloporids (corals that directly transmit symbionts to their offspring) or among those hosts displaying 'horizontal acquisition', but that associate with specific symbionts. Most host-specialized types were found to be characteristic of a particular geographic region (i.e. Okinawa versus Central GBR versus Southern GBR). The mode of symbiont acquisition may play an important role in how symbiont composition may shift in west Pacific host communities in response to climate change. There is no indication that recent episodes of mass bleaching have provoked changes in host-symbiont combinations from the central GBR.

Monocular transparency is not a new form of unpaired stereopsis

Howard and Duke (2003) generated stereograms with a grey transparent square offset from a vertical bar in one eye, and a vertical bar with a gap in the other eye. They argued that these displays were 'without conventional disparity' and that the metrical depth experienced was a new form of unpaired stereopsis due to 'transparency rather than occlusion'. Another possibility is that the perceived depth in these displays was obtained from horizontal contours. To test this possibility, we generated three displays that contained similar horizontal contourterminations, but were inconsistent with transparency. Reliable depth was seen in all stimuli. We conclude that in our stimuli, and those of Howard and Duke, transparency is not responsible for the perception of depth, which appears to be based instead on disparate horizontal contour terminations. Our results also show that disparate contours of opposite contrast polarity can generate depth.

Observations of anti-winds in a deep alpine valley, Lake Tekapo, New Zealand

This paper presents observations of summertime anti-winds monitored under ideal conditions in the Lake Tekapo hydro-catchment situated in the central Southern Alps, New Zealand. Onset and cessation of anti-winds was observed to coincide with the change in phase of the surface limbs of thermally generated valley and mountain winds under settled anti-cyclonic conditions. Anti-winds were best developed in the early morning before surface heating and associated convective mixing of the valley atmosphere began to mask the boundaries between the surface based limb of the mountain-valley wind and the corresponding anti-wind. By mid-day, the anti-valley wind exceeded the height of the surrounding ridgeline and became embedded in the topographically channeled gradient wind. Observations presented here have both theoretical and applied implications with regard to the development of thermally generated wind systems in deep alpine valleys, and their role in the dispersion of air pollution.

How configurations of binocular disparity determine whether stereoscopic slant or stereoscopic occlusion is seen

A partially occluded contour and a slanted contour may generate identical binocular horizontal disparities. We investigated conditions promoting an occlusion resolution indicated by an illusory contour in depth along the aligned ends of horizontally disparate line sets. For a set of identical oblique lines with a constant width added to one eye's view, strength, depth, and stability of the illusory contour were poor, whereas for oblique lines of alternating orientations the illusory contours were strong, indicating a reliance on vertical size disparities rather than vertical positional disparities in generating perceived occlusion. For horizontal lines, occlusion was seen when the lines were of different lengths and absolute width disparity was invariant across the set. In all line configurations, when the additional length was on the wrong eye to be attributed to differential occlusion, lines appeared slanted consistent with their individual horizontal disparities. This rules out monocular illusory contours as the determining factor.

Spatial bias and temporal lag in tactile space perception

Constructing a veridical spatial map by touch poses at least two problems for a perceptual system. First, as the hand is moved through space, the locations of features may be displaced if there is an uncorrected lag between the moment the hand encounters a feature and the time that feature is encoded on a spatial map. Second, due to the sequential nature of the process, some form of memory, which itself may be subject to spatial distortions, is required for integration of spatial samples. We investigated these issues using a task involving active haptic exploration with a stylus swept back and forth in the horizontal plane at the wrist. Remembered locations of tactile targets were shifted towards the medial axis of the forearm, suggesting a central tendency in haptic spatial memory, while evidence for a displacement of perceived locations in the direction of sweep motion was consistent with processing delays.

A synthesis of the findings from the Quake Impact Study: A two year investigation of the psychosocial sequelae of the 1989 Newcastle earthquake

This paper summarises the major findings from the Quake Impact Study (QIS), a four-phase longitudinal project that was conducted in the aftermath of the 1989 Newcastle (Australia) earthquake. A total of 3,484 subjects participated in at least one component of the QIS, comprising a stratified sample of 3,007 drawn from community electoral rolls and 477 from specially targeted supplementary samples (the injured, the displaced, the owners of damaged businesses, and the helpers). Subjects' initial earthquake experiences were rated in terms of weighted indices of exposure to threat and disruption. Psychological morbidity was measured at each phase using the General Health Questionnaire (GHQ-12) and the Impact of Event Scale (IES). Selected findings and key conclusions are presented for each of six areas of investigation: service utilisation during the first 6 months post-disaster; patterns of earthquake experience and short-term (6-month) psychosocial outcome; earthquake exposure and medium term (2-year) psychosocial outcome; vulnerability factors and medium-term psychosocial outcome: specific community groups at increased risk (e.g., the elderly and immigrants from non-English-speaking backgrounds); the effects of stress debriefing for helpers. Threshold morbidity (i.e., likely caseness) rates are also presented for a broad range of subgroups. In addition to presenting an overview of the QIS, this paper synthesises the major findings and discusses their implications for future disaster management and research from a mental health perspective.

The Edinburgh-Cape Blue Object Survey - I. Description of the survey

The Edinburgh-Cape Blue Object Survey is a major survey to discover blue stellar objects brighter than B similar to 18 in the southern sky. It is planned to cover an area of sky of 10 000 deg(2) with \b\ > 30 degrees and delta < 0 degrees. The blue stellar objects are selected by automatic techniques from U and B pairs of UK Schmidt Telescope plates scanned with the COSMOS measuring machine. Follow-up photometry and spectroscopy are being obtained with the SAAO telescopes to classify objects brighter than B = 16.5. This paper describes the survey, the techniques used to extract the blue stellar objects, the photometric methods and accuracy, the spectroscopic classification, and the limits and completeness of the survey.

Transfer of an autonomously replicating vector between vegetatively incompatible biotypes of Colletotrichum gloeosporioides

Previous research has indicated that biotypes A and B of Colletotrichum gloeosporioides that infect Stylosanthes spp. in Australia are asexual and vegetatively incompatible. Selectable marker genes conferring resistance either to hygromycin or phleomycin were introduced into isolates of these biotypes. Vectors conferring resistance to hygromycin and carrying telomeric sequences from Fusarium oxysporum replicated autonomously in C. gloeosporioides and gave frequencies of transformation 100-times higher than vectors that integrated into the genome. Monoconidial colonies resistant to both antibiotics were recovered when hygromycin-resistant biotype-A transformants carrying an autonomously replicating vector were paired in culture with a phleomycin-resistant biotype-B transformant carrying integrative vector sequences. Molecular analysis of double antibiotic-resistant progeny indicated that they contained the autonomous vector in a biotype-B genetic background, Results indicate that transfer of the autonomous vector had occurred from biotype A to biotype B, demonstrating the potential for transfer of genetic information between these biotypes.

A simple solution of the laminar dam break wave

The sudden release of a mass of fluid in a channel generates a highly unsteady flow motion, called dam break wave. While industrial fluids exhibit sometimes non-Newtonian behaviours, the viscous fluid flow assumption remains a useful approximation for simplified analyses. In this study, new solutions of laminar dam break wave are proposed for a semi-infinite reservoir based upon the method of characteristics. The solutions yield simple explicit expressions of the wave front location, wave front celerity and instantaneous free-surface profiles that compare favourably with experimental observations. Both horizontal and sloping channel configurations are treated. The simplicity of the equations may allow future extension to more complicated fluid flows.

Infiltration effects on sediment mobility under waves

Sediment mobility measurements with a horizontal sand bed under non-breaking waves are reported. Conditions include no seepage and steady downward seepage corresponding to head gradients up to 2.5. The results indicate that infiltration tends to inhibit sediment mobility for a horizontal bcd of 0.2 mm quartz sand exposed to moderated wave induced bed shear stresses. The effect is weak for the parameter range of the present study. The two opposing effects of shear stress increase due to boundary layer thinning and the stabilizing downward drag are successfully accounted for through the modified Shields parameter of Nielsen [Nielsen, P., 1997. Coastal groundwater dynamics. Proc. Coastal Dynamics '97, Plymouth, ASCE, Dp, 546-555] using coefficients derived from independent studies. That is, from the shear stress experiments of Conley [Conley, D.C., 1993. Ventilated oscillatory boundary layers. PhD Thesis, University of California, San Diego, 74 pp.] and the slope stability experiments of Martin and Aral [Martin, C.S. and M.M. Aral, 1971. Seepage force on interfacial bed particles. J. Hydraulics Div., proc. ASCE, Vol. 97, No. Hy7, pp. 1081-1100]. (C) 2001 Elsevier Science B.V. All rights reserved.

The application of electron microscopy and microanalysis in conjunction with organic petrology to further the understanding of organic-mineral association: examples from Mount Isa and McArthur Basins, Australia

Organic petrology supported by electron microscopical and micro-analytical techniques was applied to organic matter in Proterozoic sediments to better understand hydrothermal processes responsible for ore-grade mineralisation. It was shown that organic maturation was not only closely linked to the geological history of the sediments, but also highlighted heat transfer by convection as differentiated from conduction solely through sediment burial and step-wise subsidence. Water-rock ratios effect organic maturation in hydrothermal systems, and erratic reflectance profiles are indicators of convective heat transfer. Identification and characterisation of organic materials in terms of source rock and migrated hydrocarbons was shown to be a powerful tool in reconstructing the thermal history of sediments, identifying hydrothermal episodes, fluid pathways and heat source in the northern Australian Proterozoic basins. Higher reflectance of organic matter towards the central parts of the Mount Isa Basin and some of the most northerly parts point to proximity to higher heat flow at times, in contrast to relatively low temperatures prevailing in the western parts of the basin, next to the Murphy Inlier. A close correlation shown between peak organic reflectance values and super-sequence boundaries farther highlighted the valuable information to be gained from organic petrology, by allowing the separation of processes responsible for metal dissolution and transportation from those inducing precipitation. (C) 2001 Elsevier Science B.V All rights reserved.