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.

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.

Rare earth element fractionation in authigenic illite–smectite from Late Permian clastic rocks, Bowen Basin, Australia: implications for physico-chemical environments of fluids during illitization

REE analyses were performed on authigenic illitic clay. minerals from Late Permian mudrocks, sandstones and bentonites from the Bowen Basin (Australia). The mixed-layer illite-smectite exhibit REE patterns with an obvious fractionation of the HREE from the LREE and MREE, which is an apparent function of degree of illitization reaction. The highly illitic (R greater than or equal to 3) illite-smectite from the northern Bowen Basin show a depletion of LREE relative to the less illitic (R=0 and 1) clays. In contrast, an enrichment of HREE for the illite-rich clays relative to less. illitic clays is evident for the southern Bowen Basin samples. The North American Shale Composite-normalized (La/Lu)(sn) ratios show negative correlations with the illite content in illite-smectite and positive correlations with the delta(18)O values of the clays for both the northern and southern Bowen Basin samples. These correlations indicate that the increasing depletion of LREE in hydrothermal fluids is a function of increasing water/rock ratios in the northern Bowen Basin. Good negative correlations between (La/Lu)(sn) ratios and illite content in illite-smectite from the southern Bowen Basin suggest the involvement of fluids with higher alkalinity and higher pH in low water/ rock ratio conditions. Increasing HREE enrichment with delta(18)O decrease indicates the effect of increasing temperature at low water/rock ratios in the southern Bowen Basin. Results of the present study confirm the conclusions of some earlier studies suggesting that REE in illitic clay minerals are mobile and fractionated during illitization and that this fact should be considered in studies of sedimentary processes and in identifying provenance. Moreover, our results show that REE systematic of illitic clay minerals can be applied as an useful technique to gain information about physico-chemical conditions during thermal and fluid flow events in certain sedimentary basins. (C) 2003 Elsevier Science B.V. All rights reserved.

Numerical investigation of seawater intrusion at Gooburrum, Bundaberg, Queensland, Australia

Seawater intrusion in coastal agricultural areas due to groundwater abstraction is a major environmental problem along the northeastern coast of Australia. Management options are being explored using numerical modelling, however, questions remain concerning the appropriate level of sophistication in models, choice of seaward boundary conditions, and how to accommodate heterogeneity and data uncertainty. The choice of seaward boundary condition is important since it affects the amount of salt transported into the aquifers and forms the focus of the present study. The impact of this boundary condition is illustrated for the seawater-intrusion problem in the Gooburrum aquifers, which occur within Tertiary sedimentary strata. A two-dimensional variable-density groundwater and solute-transport model was constructed using the computer code 2DFEMFAT (Cheng et al. 1998). The code was tested against an experiment for a steady-state freshwater-saltwater interface and against the Elder (Elder 1967) free-convection problem. Numerical simulations show that the imposition of the commonly-used equivalent hydrostatic freshwater heads, combined with a constant salt concentration at the seaward boundary, results in overestimated seawater intrusion in the lower Gooburrum aquifer. Since the imposition of this boundary condition allows water flow across the boundary, which subsequently takes salt into the aquifer, a careful check is essential to estimate whether too much mass of salt is introduced.

Geochemical changes recorded in Lynch's Crater, Northeastern Australia, over the past 50 ka

There are many geochemical reconstructions of environmental change in the mid and high latitudes but relatively few in the tropical latitudes, despite their considerable potential for reconstructing environmental processes that cannot be identified using more traditional proxies. Here we present one reconstruction of environmental change for the tropics. This reconstruction covers the past 50 ka using a suite of geochemical data from the high-resolution sequence of Lynch's Crater in northeast Queensland, Australia, a region highly sensitive to El Nino-Southern Oscillation (ENSO) activity. The 23 major oxides and trace elements measured Could be summarised by extracting three axes using principal components analysis (accounting for 72% of the variability). The data indicate that the greatest variability in the geochemical data accounted for erosional activity within the catchment that was associated with past changes in the frequency of ENSO activity (though this was less sensitive during wetter periods, probably as a result of buffering by high vegetation cover). The remaining variability was largely explained by elements that form complexes with organic compounds (e.g., humic acids) and those that are important nutrients for specific vegetation types (and therefore a measure of vegetation distribution). For more detailed reconstructions, further work is required to disentangle the complex controls of clements within sedimentary sequences. (c) 2005 Elsevier B.V. All rights reserved.

Evidence of a spotted fever-like rickettsia and a potential new vector from northeastern Australia

A spotted fever-like rickettsia was identified in a Hemaphysalis tick by polymerase chain reaction (PCR) amplification and sequencing of the 16S rDNA, ompA, and ompB genes. A comparison of these nucleotide sequences with those of other spotted fever group (SFG) rickettsiae revealed that the Hemaphysalis tick rickettsia was distinct from other previously reported strains. Phylogenetic analysis based on both ompA and ompB also indicates that the strain’s closest relatives are the agents of Thai tick typhus (Rickettsia honei strain TT-118) and Flinders Island spotted fever (R. honei). This study represents the first report of an R. honei-like agent from a Hemaphysalis tick in Australia and of a spotted fever group rickettsia from Cape York Peninsula, Queensland.

Sedimentologic, petrographic, and sulfur isotope constraints on fine-grained pyrite formation at Mount Isa Mine and environs, Northwest Queensland, Australia

Fine-grained pyrite is the earliest generation of pyrite and the most abundant sulfide within the Urquhart Shale at Mount Isa, northwest Queensland. The pyrite is intimately interbanded with ore-grade Pb-Zn miner alization at the Mount Isa mine but is also abundant north and south of the mine at several stratigraphic horizons within the Urquhart Shale. Detailed sedimentologic, petrographic, and sulfur isotope studies of the Urquhart Shale, mostly north of the mine, reveal that the fine-grained pyrite (delta(34)S = -3.3 to +26.3 parts per thousand) formed by thermochemical sulfate reduction during diagenesis. The sulfate source was local sulfate evaporites, pseudo morphs of which are present throughout the Urquhart Shale (i.e., gypsum, anhydrite, and barite). Deep-burial diagenetic replacement of these evaporites resulted in sulfate-bearing ground waters which migrated parallel to bedding. Fine-grained pyrite formed where these fluids infiltrated and then interacted with carbon-rich laminated siltstones. Comparison of the sulfur isotope systematics of fine-grained pyrite and spatially associated base metal sulfides from the Mount Isa Pb-Zn and Cu orebodies indicates a common sulfur source of ultimately marine origin for all sulfide types. Different sulfur isotope ratio distributions for the various sulfides are the result of contrasting formation mechanisms and/or depositional conditions rather than differing sulfur sources. The sulfur isotope systematics of the base metal and associated iron sulfide generations are consistent with mineralization by reduced hydrothermal fluids, perhaps generated by bulk reduction of evaporite-sourced sulfate-bearing waters generated deeper in the Mount Isa Group, the sedimentary sequence which contains the Urquhart Shale. The available sulfur isotope data from the Mount Isa orebodies are consistent with either a chemically and thermally zoned, evolving Cu-Pb-Zn system, or discrete Cu and Pb-Zn mineralizing events linked by a common sulfur source.

The last glacial cycle from the humid tropics of northeastern Australia: comparison of a terrestrial and a marine record

A detailed pollen record from the Ocean Drilling Program Site 820 core, located on the upper part of the continental slope off the coast of northeast Queensland, was constructed to compare with the existing pollen record from Lynch's Crater on the adjacent Atherton Tableland and allow the production of a regional picture of vegetation and environmental change through the last glacial cycle. Some broad similarities in patterns of vegetation change are revealed, despite the differences between sites and their pollen catchments, which can be related largely to global climate and sea-level changes. The original estimated time scale of the Lynch's Crater record is largely confirmed from comparison with the more thoroughly dated ODP record. Conversely, the Lynch's Crater pollen record has assisted in dating problematic parts of the ODP record. In contrast to Lynch's Crater, which reveals a sharp and sustained reduction in drier araucarian forest around 38,000 yrs BP, considered to have been the result of burning by Aboriginal people, the ODP record indicates, most likely, a stepwise reduction, dating from 140,000 yrs BP or beyond. The earliest reduction shows lack of a clear connection between Araucaria decline and increased burning and suggests that people may not have been involved at this stage. However, a further decline in araucarian forest, possibly around 45,000 yrs BP, which has a more substantial environmental impact and is not related to a time of major climate change, is likely, at least partially, the result of human burning. The suggestion, from the ODP core oxygen isotope record, of a regional sea-surface temperature increase of around 4 degrees C between about 400,000 and 250,000 yrs BP, may have had some influence on the overall decline in Araucaria and its replacement by sclerophyll vegetation. (C) 2000 Elsevier Science B.V. All rights reserved.

Yarrol terrane of the northern New England Fold Belt: forearc or backarc?

The Upper Devonian to Lower Carboniferous volcanosedimentary rocks of the Yarrol terrane of the northern New England Fold Belt have previously been ascribed to a forearc basin setting. New data presented here, however, suggest that the Yarrol terrane developed as a backarc basin during the Middle to early Late Devonian. Based on field studies, we recognise four regionally applicable strati graphic units: (i) a basal, ?Middle to Upper Devonian submarine mafic volcanic suite (Monal volcanic facies association); (ii) the lower Frasnian Lochenbar beds that locally unconformably overlie the Monal volcanic facies association: (iii) the Three Moon Conglomerate (Upper Devonian - Lower Carboniferous): and (iv) the Lower Carboniferous Rockhampton Group characterised by the presence of oolitic limestone. Stratigraphic and compositional differences suggest the Monal volcanic facies association post-dates Middle Devonian silicic-dominated magmatism that was coeval with gold-copper mineralisation at Mt Morgan. The Lochenbar beds, Three Moon Conglomerate and Rockhampton Group represent a near-continuous sedimentary record of volcanism that changed in composition and style from mafic effusive (Late Devonian) to silicic explosive volcanism (Early Carboniferous). Palaeocurrent data from the Three Moon Conglomerate and Rockhampton Group indicate dispersal of sediment to the west and northwest, and are inconsistent with derivation from a volcanic-are source situated to the west (Connors-Auburn Arch). Geochemical data show that the Monal volcanic facies association ranges from tholeiitic subalkaline basalts to calc-alkaline basaltic andesite. Trace and rare-earth element abundances are distinctly MORE-like (e.g, light rare earth element depletion), with only moderate enrichment of the large-ion lithophile elements in some units, and negative Nb anomalies, suggesting a subduction-related signature. Basalts of the Monal volcanic facies association are best described as transitional between calc-alkali basalts and N-MORB. The elevated high field strength element contents (e.g. Zr, Y, Ti) are higher than modern island-are basalts, but comparable to basalts that floor modern backarc basins. This geochemical study, coupled with stratigraphic relationships, suggest that the eruption of backarc basin basalts followed widespread Middle Devonian, extension-related silicic magmatism (e.g. Retreat Batholith, Mt Morgan), and floored the Yarrol terrane. The Monal volcanic facies association thus shows similarities in its tectonic environment to the Lower Permian successions (e.g. Rookwood Volcanics) of the northern New England Fold Belt. These mafic volcanic sequences are interpreted to record two backarc basin-forming periods (Middle - Late Devonian and Late Carboniferous - Early Permian) during the Late Palaeozoic history of the New England Orogen. Silicic-dominated explosive volcanism, occurring extensively across the northern New England Fold Belt in the Early Carboniferous (Varrol terrane, Campwyn Volcanics, Drummond and Burdekin Basins), reflects another period of crustal melting and extension, most likely related to the opening of the Drummond Basin.

Anthropogenic and Natural Organohalogen Compounds in Blubber of Dolphins and Dugongs (Dugong dugon) from Northeastern Australia

A range of organohalogen compounds (10 polychlorinated biphenyl [PCB] congeners, DDT and metabolites, chlordane-related compounds, the potential natural organochlorine compound Q1, toxaphene, hexachlorobenzene, hexachlorocyclohexanes, dieldrin, and several yet unidentified brominated compounds) were detected in the blubber of four bottlenose dolphins (Tursiops truncatus), one common dolphin (Delphinus delphis), and seven dugongs (Dugong dugon), as well as in adipose tissue of a green turtle (Chelonia mydas) and a python (Morelia spilota) from northeast Queensland (Australia). The green turtle and dugongs accumulated lower organohalogen levels than the dolphins. Lower levels in dugongs were expected because this species is exclusively herbivorous. Highest PCB and DDT levels recorded in dugongs were 209 and 173 mug/kg lipids, respectively. Levels of the nonanthropogenic heptachlorinated compound Q1 (highest level in dugongs was 160 mug/kg lipids) were estimated using the ECD response factor of trans-nonachlor. Highest organohalogen levels were found in blubber of dolphins for sumDDT (575-52,500 mug/kg) and PCBs (600-25,500 mug/kg lipids). Furthermore, Q1 was a major organohalogen detected in all samples analyzed, ranging from 450 -9,100 mug/kg lipids. The highest concentration of Q1 determined in this study represents the highest concentration reported to date in an environmental sample. Levels of chlordane-related compounds were also high (280-7,700 mug/kg, mainly derived from trans-nonachlor), but concentrations of hexachlorobenzene, hexachlorocyclohexanes, dieldrin, and toxaphene were relatively low and contributed little to the overall organohalogen contamination. Furthermore, a series of three major (BC-1, BC-2, and BC-3) and six minor (BC-4 through BC-9) unknown brominated compounds were observable by extracting m/z 79 and m/z 81 from the GC/ECNI-MS full scan run. Structural proposals were made for the two major recalcitrant compounds (referred to as BC-1 and BC-2). BC-2 appears to be a tetrabromo-methoxy-diphenylether (512 u) and BC-1 has 14 u (corresponding with an additional CH2 group) more relative to BC-1. In general the organohalogen pattern observed in blubber of dolphins was different compared to similar samples from other locations in the world, which is apparent from the fact that the four major abundant signals in the GC/ECD chromatogram. of D. delphis originated from the four unknown compounds Q1, BC-1, BC-2, and BC-3.

Stratigraphy, facies architecture and tectonic implications of the Upper Devonian to Lower Carboniferous Campwyn Volcanics of the northern New England Fold Belt

Upper Devonian to Lower Carboniferous strata of the Campwyn Volcanics of east central Queensland preserve a substantial sequence of first-cycle volcaniclastic sedimentary and coeval volcanic rocks that record prolonged volcanic activity along the northern New England Fold Belt. The style and scale of volcanism varied with time, producing an Upper Devonian sequence of mafic volcano-sedimentary rocks overlain by a rhyolitic ignimbrite-dominated sequence that passes upward into a Lower Carboniferous limestone-bearing sedimentary sequence. We define two facies associations for the Campwyn Volcanics. A lower facies association is dominated by mafic volcanic-derived sedimentary breccias with subordinate primary mafic volcanic rocks comprising predominantly hyaloclastite and peperite. Sedimentary breccias record episodic and high energy, subaqueous depositional events with clastic material sourced from a mafic lava-dominated terrain. Some breccias contain a high proportion of attenuated dense, glassy mafic juvenile clasts, suggesting a syn-eruptive origin. The lower facies association coarsens upwards from a lithic sand-dominated sequence through a thick interval of pebble- to boulder-grade polymict volcaniclastic breccias, culminating in facies that demonstrate subaerial exposure. The silicic upper facies association marks a significant change in eruptive style, magma composition and the nature of eruptive sources, as well as the widespread development of subaerial depositional conditions. Crystal-rich, high-grade, low- to high-silica rhyolite ignimbrites dominate the base of this facies association. Biostratigraphic age controls indicate that the ignimbrite-bearing sequences are Famennian to lower-mid Tournaisian in age. The ignimbrites represent extra-caldera facies with individual units up to 40 m thick and mostly lacking coarse lithic breccias. Thick deposits of pyroclastic material interbedded with fine-grained siliceous sandstone and mudstone (locally radiolarian-bearing) were deposited from pyroclastic flows that crossed palaeoshorelines or represent syn-eruptive, resedimented pyroclastic material. Some block-bearing lithic-pumice-crystal breccias may also reflect more proximal subaqueous silicic explosive eruptions. Crystal-lithic sandstones interbedded with, and overlying the ignimbrites, contain abundant detrital volcanic quartz and feldspar derived from the pyroclastic deposits. Limestone is common in the upper part of the upper facies association, and several beds are oolitic (cf. Rockhampton Group of the Yarrol terrane). Overall, the upper facies association fines upward and is transgressive, recording a return to shallow-marine conditions. Palaeocurrent data from all stratigraphic levels in the Campwyn Volcanics indicate that the regional sediment-dispersal direction was to the northwest, and opposed to the generally accepted notion of easterly sediment dispersal from a volcanic arc source. The silicic upper facies association correlates in age and lithology to Early Carboniferous silicic volcanism in the Drummond (Cycle 1) and Burdekin Basins, Connors Arch, and in the Yarrol terranes of eastern Queensland. The widespread development of silicic volcanism in the Early Carboniferous indicates that silicic (rift-related) magmatism was not restricted to the Drummond Basin, but was part of a more substantial silicic igneous province.

A broad-scale analysis of links between coastal fisheries production and mangrove extent: A case-study for northeastern Australia

The paradigm that mangroves are critical for sustaining production in coastal fisheries is widely accepted, but empirical evidence has been tenuous. This study showed that links between mangrove extent and coastal fisheries production could be detected for some species at a broad regional scale (1000s of kilometres) on the east coast of Queensland, Australia. The relationships between catch-per-unit-effort for different commercially caught species in four fisheries (trawl, line, net and pot fisheries) and mangrove characteristics, estimated from Landsat images were examined using multiple regression analyses. The species were categorised into three groups based on information on their life history characteristics, namely mangrove-related species (banana prawns Penaeus merguiensis, mud crabs Scylla serrata and barramundi Lates calcarifer), estuarine species (tiger prawns Penaeus esculentus and Penaeus semisulcatus, blue swimmer crabs Portunus pelagicus and blue threadfin Eleutheronema tetradactylum) and offshore species (coral trout Plectropomus spp.). For the mangrove-related species, mangrove characteristics such as area and perimeter accounted for most of the variation in the model; for the non-mangrove estuarine species, latitude was the dominant parameter but some mangrove characteristics (e.g. mangrove perimeter) also made significant contributions to the models. In contrast, for the offshore species, latitude was the dominant variable, with no contribution from mangrove characteristics. This study also identified that finer scale spatial data for the fisheries, to enable catch information to be attributed to a particular catchment, would help to improve our understanding of relationships between mangroves and fisheries production. (C) 2005 Elsevier B.V. All rights reserved.

Pollen transport and deposition in riverine and marine environments within the humid tropics of northeastern Australia

Mechanisms of pollen transport in the humid tropics region of northeastern Australia were investigated to support the interpretation of a long Quaternary pollen record from ODP Site 820 located on the adjacent continental slope. Pollen analysis of surface sediment samples from the channels of two major river catchments demonstrated internal consistency in pollen spectra and little fluvial pollen sorting in relation to sediment variation. Differences in modem pollen spectra between catchments reflect existing variation in vegetation cover that, in turn, reflects climatic differences between catchments. Recent pollen spectra from top samples of the ODP core have sufficient in common with the riverine samples to suggest that the rivers are contributing a major pollen component to the offshore sediments, but these have been size sorted by marine action. Recent pollen samples from core tops taken from the Grafton Passage on the continental shelf that was thought to be the major passage for pollen transport to ODP Site 820 show significant differences to both riverine and ODP samples and suggest that pollen are dispersed across the continental shelf and through the outer Great Barrier Reef system in a more complex way than anticipated. (c) 2004 Elsevier B.V. All rights reserved.