Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

A semi-empirical approach to optimise the quantity of drying aids required to spray dry sugar-rich foods

A semi-empirical linear equation has been developed to optimise the amount of maltodextrin additive (DE 6) required to successfully spray dry a sugar-rich product on the basis of its composition. Based on spray drying experiments, drying index values for individual sugars (sucrose, glucose, frutose) and citric acid were determined, and us;ng these index values an equation for model mixtures of these components was established. This equation has been tested with two sugar-rich natural products, pineapple juice and honey. The relationship was found to be valid for these products.

Determination of the hemispheric character (i.e. southern or northern latitudinal position) of equatorial ionosonde stations as deduced from long-term recordings of the minimum height of the F-region (h ' F)

The virtual (or minimum) height of the F-region (h'F), recorded over a number of solar cycles for I I equatorial and mid-latitude ionosonde stations, was used to deduce the hemispheric (i.e. southern or northern hemisphere) character of equatorial stations. The semi-annual median monthly height (h'F) variations consist of two components: major local summer maximum and winter sub-maximum (about 5 percent of the summer maximum). This hemispheric pattern was most consistently observed for equatorial stations (within 5degrees of the geomagnetic equator) in a period centred on the local midnight (21-03 LT) but was also present, to a lesser extent, at mid-latitude stations and at other time intervals. It is evident that the physical parameter h'F defines the hemispheric character of an equatorial station which has different (sometimes opposite) geographic and geomagnetic latitudes. There is a sharp transition in the latitudinal character of the stations on both sides of the equator leading to hypothesis that the equal maxima in h'F in December and June solstices are observed at a near-equator position labelled as ionosonde deduced equator (IDE). Such a signature was observed for an American equatorial (both geographic and geomagnetic) station Talara (Peru) which is an experimental support of the hypothesis. The IDE can be another useful parameter characterising the equatorial ionosphere. This finding reveals a new application of the standard ionosonde data in defining the geophysical character of equatorial stations, being an important contribution to space climatology. (C) 2002 Elsevier Science Ltd. All rights reserved.

Diversity of pigmentation in cultured human melanocytes is due to differences in the type as well as quantity of melanin

Cultured human melanocytes differ tremendously in visual pigmentation, and recapitulate the pigmentary phenotype of the donor's skin. This diversity arises from variation in type as well as quantity of melanin produced. Here, we measured contents of eumelanin (EM) and pheomelanin (PM) in 60 primary human melanocyte cultures (51 neonatal and nine adults), and correlated some of these values with the respective activity and protein levels of tyrosinase, and the melanocortin-1 receptor (MC1R) genotype. Melanocytes were classified into four phenotypes (L, L+, D, D+) as depicted by visual pigmentation using light microscopy, and by the pigmentary phenotype of the donor's skin. There were large differences in total melanin (TM) and EM, which increased progressively for L, L+, D and D+ melanocytes. TM content, the sum of EM and PM, showed a good correlation with TM measured spectrophotometrically, and with the activity and protein levels of tyrosinase. Log EM/PM ratio did not correlate with MC1R genotype. We conclude that: (i) EM consistently correlates with the visual phenotype; (ii) lighter melanocytes tend to be more pheomelanic in composition than darker melanocytes; (iii) in adult melanocyte cultures, EM correlates with the ethnic background of the donors (African-American > Indian > Caucasian); and (iv) MC1R loss-of-function mutations do not necessarily alter the phenotype of cultured melanocytes.

Timing and chemistry of fluid-flow events in the Lawn Hill Platform, Northern Australia

Mudrocks and carbonates of the Isa superbasin in the Lawn Hill platform in northern Australia host major base metal sulfide mineralization, including the giant strata-bound Century Zn-Pb deposit. Mineral paragenesis, stable isotope, and K-Ar dating studies demonstrate that long-lived structures such as the Termite Range fault acted as hot fluid conduits several times during the Paleoproterozoic and Mesoproterozoic in response to major tectonic events. Illite and chlorite crystallinity studies suggest the southern part of the platform has experienced higher temperatures (up to 300 degrees C) than similar stratigraphic horizons in the north. The irregular downhole variation of illite crystallinity values provides further information oil the thermal regime in the basin and shows that clay formation was controlled not only by temperature increase with depth but also by high water/rock ratios along relatively permeable zones. K-Ar dating of illite, in combination with other data, may indicate three major thermal events in the central and northern Lawn Hill platform Lit 1500, 1440 to 1400, and 1250 to 1150 Ma. This study did not detect the earlier Century base metal mineralizing event at 1575 Ma. 1500 Ma ages are recorded only in the south and correspond to the age of the Late Isan orogeny and deposition of the Lower Roper superbasin. They may reflect exhumation of a provenance region. The 1440 to 1300 Ma ages are related to fault reactivation and a thermal pulse at similar to 1440 to 1400 Ma possibly accompanied by fluid flow, with subsequent enhanced cooling possibly due to thermal relaxation or further crustal exhumation. The youngest thermal and/or fluid-flow event at 1250 to 1150 Ma is recorded mainly to the cast of the Tern-lite Range fault and may be related to the assembly of the Rodinian supercontinent. Fluids in equilibrium with illite that formed over a range of temperatures, at different times in different parts of the platform. have relatively uniform oxygen isotope compositions and more variable hydrogen isotope compositions (delta O-18 = 3.5-9.7 parts per thousand V-SMOW; delta D = -94 to -36 parts per thousand V-SMOW). The extent of the 180 enrichment and the variably depleted hydrogen isotope compositions suggest the illite interacted with deep-basin hypersaline brines that were composed of evaporated seawater and/or highly evolved meteoric water. Siderite is the most abundant iron-rich gangue phase in the Century Zn-Pb deposit, which is surrounded by all extensive ferroan carbonate alteration halo. Modeling suggests that the ore siderite formed at temperatures of 120 degrees to 150 degrees C, whereas siderite and ankerite in the alteration halo formed at temperatures of 150 degrees to 180 degrees C. The calculated isotopic compositions of the fluids are consistent with O-18-rich basinal brines and mixed inorganic and organic carbon Sources (6180 = 3-10 parts per thousand V-SMOW, delta C-13 = -7 to -3 parts per thousand V-PDB). in the northeast Lawn Hill platform carbonate-rich rocks preserve marine to early diagenetic carbon and oxygen isotope compositions, whereas ferroan carbonate cements in siltstones and shales in the Desert Creek borehole are O-18 and C-13 depleted relative to the sedimentary carbonates. The good agreement between temperature estimates from illite crystallinity and organic reflectance (160 degrees-270 degrees C) and inverse correlation with carbonate delta O-18 values indicates that organic maturation and carbonate precipitation in the northeast Lawn Hill platform resulted from interaction with the 1250 to 1150 Ma fluids. The calculated isotopic compositions of the fluid are consistent with evolved basinal brine (delta O-18 = 5.1-9.4 parts per thousand V-SMOW; delta C-13 = -13.2 to -3.7 parts per thousand V-PDB) that contained a variable organic carbon component from the oxidation and/or hydrolysis of organic matter in the host sequence. The occurrence of extensive O-18- and C-13-depleted ankerite and siderite alteration in Desert Creek is related to the high temperature of the 1250 to 1150 Ma fluid-flow event in the northeast Lawn Hill platform, in contrast to the lower temperature fluids associated with the earlier Century Zn-Pb deposit in the central Lawn Hill platform.

Different mineralization styles in a volcanic-hosted ore deposit: the fluid and isotopic signatures of the Mt Morgan Au-Cu deposit, Australia

Quantitative laser ablation (LA)-ICP-MS analyses of fluid inclusions, trace element chemistry of sulfides, stable isotope (S), and Pb isotopes have been used to discriminate the formation of two contrasting mineralization styles and to evaluate the origin of the Cu and Au at Mt Morgan. The Mt Morgan Au-Cu deposit is hosted by Devonian felsic volcanic rocks that have been intruded by multiple phases of the Mt Morgan Tonalite, a low-K, low-Al2O3 tonalite-trondhjemite-dacite (TTD) complex. An early, barren massive sulfide mineralization with stringer veins is conforming to VHMS sub-seafloor replacement processes, whereas the high-grade Au-Cu. ore is associated with a later quartz-chalcopyrite-pyrite stock work mineralization that is related to intrusive phases of the Tonalite complex. LA-ICP-MS fluid inclusion analyses reveal high As (avg. 8850 ppm) and Sb (avg. 140 ppm) for the Au-Cu mineralization and 5 to 10 times higher Cu concentration than in the fluids associated with the massive pyrite mineralization. Overall, the hydrothermal system of Mt Morgan is characterized by low average fluid salinities in both mineralization styles (45-80% seawater salinity) and temperatures of 210 to 270 degreesC estimated from fluid inclusions. Laser Raman Spectroscopic analysis indicates a consistent and uniform array Of CO2-bearing fluids. Comparison with active submarine hydrothermal vents shows an enrichment of the Mt Morgan fluids in base metals. Therefore, a seawater-dominated fluid is assumed for the barren massive sulfide mineralization, whereas magmatic volatile contributions are implied for the intrusive related mineralization. Condensation of magmatic vapor into a seawater-dominated environment explains the CO2 occurrence, the low salinities, and the enriched base and precious metal fluid composition that is associated with the Au-Cu. mineralization. The sulfur isotope signature of pyrite and chalcopyrite is composed of fractionated Devonian seawater and oxidized magmatic fluids or remobilized sulfur from existing sulfides. Pb isotopes indicate that Au and Cu. originated from the Mt Morgan intrusions and a particular volcanic strata that shows elevated Cu background. (C) 2002 Elsevier Science B.V. All rights reserved.

Fluid-mobile trace element constraints on the role of slab melting and implications for Archaean crustal growth models

We use published and new trace element data to identify element ratios which discriminate between arc magmas from the supra-subduction zone mantle wedge and those formed by direct melting of subducted crust (i.e. adakites). The clearest distinction is obtained with those element ratios which are strongly fractionated during refertilisation of the depleted mantle wedge, ultimately reflecting slab dehydration. Hence, adakites have significantly lower Pb/Nd and B/Be but higher Nb/Ta than typical arc magmas and continental crust as a whole. Although Li and Be are also overenriched in continental crust, behaviour of Li/Yb and Be/Nd is more complex and these ratios do not provide unique signatures of slab melting. Archaean tonalite-trondhjemite-granodiorites (TTGs) strongly resemble ordinary mantle wedge-derived arc magmas in terms of fluid-mobile trace element content, implying that they-did not form by slab melting but that they originated from mantle which was hydrated and enriched in elements lost from slabs during prograde dehydration. We suggest that Archaean TTGs formed by extensive fractional crystallisation from a mafic precursor. It is widely claimed that the time between the creation and subduction of oceanic lithosphere was significantly shorter in the Archaean (i.e. 20 Ma) than it is today. This difference was seen as an attractive explanation for the presumed preponderance of adakitic magmas during the first half of Earth's history. However, when we consider the effects of a higher potential mantle temperature on the thickness of oceanic crust, it follows that the mean age of oceanic lithosphere has remained virtually constant. Formation of adakites has therefore always depended on local plate geometry and not on potential mantle temperature.

A simple, rapid method for mapping bathymetry of small wetland basins

Many tools exist for determining and mapping the bathymetry and topography of aquatic systems, such as freshwater wetlands. However, these tools often require time-consuming survey work to produce accurate maps. In particular, the large quantity of data necessary may be prohibitive for projects where determining bathymetry is not a central focus, but instead a necessary step in achieving some other goal. We present a method to produce bathymetric surface maps with a minimum amount of effort using global positioning system receiver and laser transit survey data. We also demonstrate that this method is surprisingly accurate, given the small amount of data we use to generate the bathymetry maps.

Pendant drop formation of shear-thinning and yield stress fluids

A volume-of-fluid numerical method is used to predict the dynamics of shear-thinning liquid drop formation in air from a circular orifice. The validity of the numerical calculation is confirmed for a Newtonian liquid by comparison with experimental measurements. For particular values of Weber number and Froude number, predictions show a more rapid pinch-off, and a reduced number of secondary droplets, with increasing shear-thinning. Also a minimum in the limiting drop length occurs for the smallest Weber number as the zero-shear viscosity is varied. At the highest viscosity, the drop length is reduced due to shear-thinning, whereas at lower viscosities there is little effect of shear-thinning. The evolution of predicted drop shape, drop thickness and length, and the configuration at pinch-off are discussed for shear-thinning drops. The evolution of a drop of Bingham yield stress liquid is also considered as a limiting case. In contrast to the shear-thinning cases, it exhibits a plug flow prior to necking, an almost step-change approach to pinch-off of a torpedo shaped drop following the onset of necking, and a much smaller neck length; no secondary drops are formed. The results demonstrate the potential of the numerical model as a design tool in tailoring the fluid rheology for controlling drop formation behaviour. (c) 2006 Elsevier Inc. All rights reserved.

Dam break wave of thixotropic fluid

Thixotropy is the characteristic of a fluid to form a gelled structure over time when it is not subjected to shearing, and to liquefy when agitated. Thixotropic fluids are commonly used in the construction industry (e.g., liquid concrete and drilling fluids), and related applications include some forms of mud flows and debris flows. This paper describes a basic study of dam break wave with thixotropic fluid. Theoretical considerations were developed based upon a kinematic wave approximation of the Saint-Venant equations down a prismatic sloping channel. A very simple thixotropic model, which predicts the basic theological trends of such fluids, was used. It describes the instantaneous state of fluid structure by a single parameter. The analytical solution of the basic flow motion and theology equations predicts three basic flow regimes depending upon the fluid properties and flow conditions, including the initial degree of jamming of the fluid (related to its time of restructuration at rest). These findings were successfully compared with systematic bentonite suspension experiments. The present work is the first theoretical analysis combining the basic principles of unsteady flow motion with a thixotropic fluid model and systematic laboratory experiments.