Influence of seaward boundary condition on contaminant transport in unconfined coastal aquifers

Contaminant transport in coastal aquifers is complicated partly due to the conditions at the seaward boundary including seawater intrusion and tidal variations of sea level. Their inclusion in modelling this system will be computationally expensive. Therefore, it will be instructive to investigate the consequence of simplifying the seaward boundary condition by neglecting the seawater density and tidal variations in numerical predictions of contaminant transport in this zone. This paper presents a comparison of numerical predictions for a simplified seaward boundary condition with experimental results for a corresponding realistic one including a saltwater interface and tidal variations. Different densities for contaminants are considered. The comparison suggests that the neglect of the seawater intrusion and tidal variations does not affect noticeably the overall migration rate of the plume before it reaches the saltwater interface. However, numerical prediction shows that a more dense contaminant travels further seaward and part of the solute mass exits under the sea if the seawater density is not included. This is not consistent with the experimental result, which shows that the contaminant travels upwards to the shoreline along the saltwater interface. Neglect of seawater density, therefore, will result in an underestimation of the exit rate of solute mass around the coastline and fictitious migration paths under the seabed. For a less dense contaminant, neglect of seawater density has little effect on numerical prediction of migration paths. (C) 2001 Elsevier Science B.V. All rights reserved.

Associations involving delays (particularly long delays) between certain weather parameters and geomagnetic activity

Four sunspot-minimum periods (1963-1966, 1971-1977, 1983-1987 and 1992-1997) have been examined for the results which are presented. Using several different weather parameters, tropospheric gravity waves, enhanced cold fronts and two rainfall data sets in Eastern Australia, associations at reasonably high levels of significance have been found with enhanced geomagnetic activity (EGA). Statistically this EGA involved either short delays of several days or long delays of about 20 days. The geomagnetic parameters used were (a) the AE index (b) the hourly H component for a number of stations and (c) the daily K-P-sum value. The K-P-sum analyses have shown that the EGA associated with the delays form part of four or five cycles of recurrent geomagnetic activity for 27-day periodicities. Furthermore statistically two recurrent cycles are found to exist concurrently, one apparently related to the short delays and the other to the long delays. Periodicities of 13.5 days are created because the two sets are displaced from each other by approximately this interval. A brief reference is made to the 13.5 periodicity known to exist for geomagnetic activity and the evidence in the literature for active regions on the sun to be displaced by 180 degrees of solar longitude.

Delineation of sulphide ore-zones by borehole radar tomography at Hellyer Mine, Australia

Velocity and absorption tomograms are the two most common forms of presentation of radar tomographic data. However, mining personnel, geophysicists included, are often unfamiliar with radar velocity and absorption. In this paper, general formulae are introduced, relating velocity and attenuation coefficient to conductivity and dielectric constant. The formulae are valid for lossy media as well as high-resistivity materials. The transformation of velocity and absorption to conductivity and dielectric constant is illustrated via application of the formulae to radar tomograms from the Hellyer zinc-lead-silver mine, Tasmania, Australia. The resulting conductivity and dielectric constant tomograms constructed at Hellyer demonstrated the potential of radar tomography to delineate sulphide ore zones. (C) 2001 Elsevier Science B.V. All rights reserved.

Spontaneous-combustion propensity of New Zealand coals under adiabatic conditions

Pulverised New Zealand coal samples have been tested from an initial temperature of 40 degreesC and reacted adiabatisally in an oven with oxygen to provide a full temperature history of auto-oxidation up to the self-sustained process of combustion. This procedure produces a self-heating rate index, R-70, calculated as the ratio of the time taken to reach 70 degreesC (degreesC/h). The R-70 index is a measure of the coal's propensity to spontaneous combustion. R-70 values for New Zealand coals are much higher than any previously published results. They show a rank dependence, whereby subbituminous coals have the highest propensity to spontaneous combustion (14.91-17.23 degreesC/h). A lignite sample has an R-70 value of 7.76 degreesC/h, and high-volatile bituminous B coals have R-70 values of 0.31-2.23 degreesC/h. Samples stored for 2 years show the same rank trend. The nature of this trend is most likely a function of the internal surface area of the coal that governs the available sites for oxidation. Calculating the Suggate rank; for any New Zealand coal can be used to rare its propensity to spontaneous combustion. Resin bodies in the subbituminous coal show no propensity to spontaneous combustion. (C) 2001 Elsevier Science B.V. All rights reserved.

Welding and rheomorphism of phonolitic fallout deposits from the Las Canadas caldera, Tenerife, Canary Islands

The Las Canadas caldera is a nested collapse caldera formed by the successive migration and collapse of shallow magmatic chambers. Among the pyroclastic products of this caldera are phonolitic fallout deposits that crop out in the caldera wall and on the extracaldera slopes. These deposits exhibit an uninterrupted facies gradation from nonwelded to lava-like and record continuous volcanic deposition. Densely welded and lava-like facies result from the extreme attenuation and complete homogenization of juvenile clasts that destroy original clast outlines and any evidence of fallout deposition. Agglutination contributes significantly to the final degree of flattening observed in the welded facies. After deposition, rheomorphic flowage occurs. Emplacement temperatures for one of the welding sequences are calculated from magmatic temperatures and a model of tephra cooling during fallout. Results are 486 degreesC for the nonwelded facies and 740 degreesC for the moderately welded facies. For the same welding sequence, a cooling time between 25 and 54 days is estimated from published experimental and computational data as the possible duration of welding and rheomorphism. Following deposition and agglutination, the lava-like pyroclastic facies had the rheological properties of viscous lavas and flowed down the outer slopes away from the caldera. Some lava-like masses detached from proximal areas to more distal regions. During deposition, the eruptive style evolved from Plinian fallout to fountain-fed spatter deposition. This evolution was accompanied by a decrease in explosive power and a lower height of the eruptive column, which produce higher emplacement temperatures and more effective heat retention of pyroclasts.

New arsenite-oxidizing bacteria isolated from Australian gold mining environments - Phylogenetic relationships

Nine novel arsenite-oxidizing bacteria have been isolated from two different gold mine environments in Australia. Four of these organisms grow chemolithoautotrophically with oxygen as the terminal electron acceptor, arsenite as the electron donor, and carbon dioxide-bicarbonate as the sole carbon source. Five heterotrophic arsenite-oxidizing bacteria were also isolated, one of which was found to be both phylogenetically and physiologically identical to the previously described heterotrophic arsenite oxidizer misidentified as Alcaligenes faecalis. The results showed that this strain belongs to the genus Achromobacter. Phylogenetically, the arsenite-oxidizing bacteria fall within two separate subdivisions of the Proteobacteria. Interestingly, the chemolithoautotrophic arsenite oxidizers belong to the alpha-Proteobacteria, whereas the heterotrophic arsenite oxidizers belong to the beta-Proteobacteria.

The DMSO reductase family of microbial molybdenum enzymes; Molecular properties and role in the dissimilatory reduction of toxic elements

The dimethylsulfoxide (DMSO) reductase family of molybdenum enzymes is a large and diverse group that is found in bacteria and archaea. These enzymes are characterised by a bis(molybdopterin guanine dinucleotide)Mo form of the molybdenum cofactor, and they are particularly important in anaerobic respiration including the dissimilatory reduction of certain toxic oxoanions. The structural and phylogenetic relationship between the proteins of this family is discussed. High-resolution crystal structures of enzymes of the DMSO reductase family have revealed a high degree of similarity in tertiary structure. However, there is considerable variation in the structure of the molybdenum active site and it seems likely that these subtle but important differences lead to the great diversity of function seen in this family of enzymes. This diversity of catalytic capability is associated with several distinct pathways of electron transport.

Marine and estuarine reservoir effects in central Queensland, Australia: Determination of Delta R values

As a component of archaeological investigations on the central Queensland coast, a series of five marine shell specimens live-collected between A.D. 1904 and A.D. 1929 and 11 shell/ charcoal paired samples from archaeological contexts were radiocarbon dated to determine local DeltaR values. The object of the study was to assess the potential influence of localized variation in marine reservoir effect in accurately determining the age of marine and estuarine shell from archaeological deposits in the area. Results indicate that the routinely applied DeltaR value of -5 +/- 35 for northeast Australia is erroneously calculated. The determined values suggest a minor revision to Reimer and Reimer's (2000) recommended value for northeast Australia from DeltaR = +11 +/- 5 to + 12 +/- 7, and specifically for central Queensland to DeltaR = +10 +/- 7, for near-shore open marine environments. In contrast, data obtained from estuarine shell/charcoal pairs demonstrate a general lack of consistency, suggesting estuary-specific patterns of variation in terrestrial carbon input and exchange with the open ocean. Preliminary data indicate that in some estuaries, at some time periods, a DeltaR value of more than - 155 +/- 55 may be appropriate, In estuarine contexts in central Queensland, a localized estuary-specific correction factor is recommended to account for geographical and temporal variation in C-14 activity. (C) 2002 Wiley Periodicals.

High-frequency climatic oscillations recorded in a Holocene coral reef at Leizhou Peninsula, South China sea

A detailed study of the Goniopora reef profile at Dengloujiao, Xuwen County, Leizhou Peninsula, the northern coast of the South China Sea suggests that a series of high-frequency, large-amplitude and abrupt cold events occurred during the Holocene Hypsithermal, an unusual phenomenon termed Leizhou Events in this paper. This period (corresponding to C-14 age of 6.2 -6.7 kaBP or calendar age of 6.7-7.2 kaBP), when the climatic conditions were ideal for coral. reefs to develop, can be divided into at least nine stages. Each stage (or called a climate optimum), lasting about 20 to 50 a, was terminated by an abrupt cold nap and (or) a sea-level lowering event in winter, leading to widespread emergence and death of the Goniopora corals, and growth discontinuities on the coral surface. Such a cyclic process resulted in the creation of a > 4m thick Goniopora reef flat. During this period, the crust subsided periodically but the sea level was rising. The reef profile provides valuable archives for the study of decadal-scale mid-Holocene climatic oscillations in the tropical area of South China. Our results provide new evidence for high-frequency climate instability in the Holocene Hypsithermal, and challenge the traditional understanding of Holocene climate.

Ductile fabrics in the zone of active oblique convergence near the Alpine Fault, New Zealand: identifying the neotectonic overprint

The mid-crustal Alpine Schist in central Southern Alps, New Zealand has been exhumed during the past similar to3 m.y. on the hanging wall of the oblique-slip Alpine Fault. These rocks underwent ductile deformation during their passage through the similar to 150-km-wide Pacific-Australia plate boundary zone. Likely to be Cretaceous in age, peak metamorphism predates the largely Pliocene and younger oblique convergence that continues to uplift the Southern Alps today. Late Cenozoic ductile deformation constructively reinforced a pre-existing fabric that was well oriented to accommodate a dextral-transpressive overprint. Quartz microstructures below a recently exhumed brittle-ductile transition zone reflect a late Cenozoic increment of ductile strain that was distributed across deeper levels of the Pacific Plate. Deformation was transpressive, including a dextral-normal shear component that bends and rotates a delaminated panel of Pacific Plate crust onto the oblique footwall ramp of the Alpine Fault. Progressive ductile shear in mylonites at the base of the Pacific Plate overprints earlier fabrics in a dextral-reverse sense, a deformation that accompanies translation of the schists up the Alpine Fault. Ductile shear along that structure affects not only the 12-km-thick section of Alpine mylonites, but is distributed across several kilometres of overlying nonmylonitic rocks. (C) 2001 Elsevier Science Ltd. All rights reserved.

Kinematics of oblique collision and ramping inferred from microstructures and strain in middle crustal rocks, central Southern Alps, New Zealand

The hanging wall of the Alpine Fault near Franz Josef Glacier has been exhumed during the past similar to2-3 m.y. providing a sample of the ductilely deformed middle crust of a modem obliquely convergent orogen. Presently exposed rocks of the Pacific Plate are inferred to have undergone several phases of ductile deformation as they moved westward above a mid-crustal detachment. Initially they were transpressed across the outboard part of the orogen, resulting in oblate fabrics with a down-dip stretch. Later, they encountered the Alpine Fault, experiencing an oblique-slip backshearing on vertical planes. This escalator-like deformation tilted and thinned the incoming crust onto that crustal-scale oblique ramp. This style of hanging wall deformation may affect only the most rapidly uplifting, central part of the Southern Alps because of the low flexural rigidity of the crust in that region and its displacement over a relatively sharp ramp-angle at depth. A 3D transpressive flow affected mylonites locally near the fault, but their shear direction remained parallel to plate motion, ruling out ductile 'extrusion' as an important process in this orogen. Outside the mylonite zone, late Cenozoic shortening is inferred to be modest (30-40%), as measured from deformation of younger biotite grains. Oblique collision is dominated by translation on the Alpine Fault, and rocks migrate rapidly through the deforming zone, preventing the accumulation of large finite strains. Transpression may play a minor role in oblique collision. (C) 2001 Elsevier Science Ltd. All rights reserved.

Confined-unconfined flow in a horizontal confined aquifer

Mixed confined and unconfined groundwater flow occurs in a bounded initially dry aquifer when the hydraulic head at the side boundary suddenly rises above the elevation of the aquifer's top boundary. The flow problem as modelled by the Boussinesq equation is non-trivial because of the involvement of two moving boundaries. The transformed equation (based on a similarity transformation) can, however, be dealt with more easily. Here, we present an approximate analytical solution for this flow problem. The approximate solution is compared with an 'exact' numerical solution and found to be a very accurate description for describing the mixed confined and unconfined flow in the confined aquifer. (C) 2002 Elsevier Science B.V. All rights reserved.

Conditional simulation of random fields by successive residuals

This paper presents a new approach to the LU decomposition method for the simulation of stationary and ergodic random fields. The approach overcomes the size limitations of LU and is suitable for any size simulation. The proposed approach can facilitate fast updating of generated realizations with new data, when appropriate, without repeating the full simulation process. Based on a novel column partitioning of the L matrix, expressed in terms of successive conditional covariance matrices, the approach presented here demonstrates that LU simulation is equivalent to the successive solution of kriging residual estimates plus random terms. Consequently, it can be used for the LU decomposition of matrices of any size. The simulation approach is termed conditional simulation by successive residuals as at each step, a small set (group) of random variables is simulated with a LU decomposition of a matrix of updated conditional covariance of residuals. The simulated group is then used to estimate residuals without the need to solve large systems of equations.

Testing landscape metrics as indicators of habitat loss and fragmentation in continuous eucalypt forests (Queensland, Australia)

Landscape metrics are widely applied in landscape ecology to quantify landscape structure. However, many are poorly tested and require rigorous validation if they are to serve as reliable indicators of habitat loss and fragmentation, such as Montreal Process Indicator 1.1e. We apply a landscape ecology theory, supported by exploratory and confirmatory statistical techniques, to empirically test landscape metrics for reporting Montreal Process Indicator 1.1e in continuous dry eucalypt forests of sub-tropical Queensland, Australia. Target biota examined included: the Yellow-bellied Glider (Petaurus australis); the diversity of nectar and sap feeding glider species including P. australis, the Sugar Glider P. breviceps, the Squirrel Glider P. norfolcensis, and the Feathertail Glider Acrobates pygmaeus; six diurnal forest birds species; total diurnal bird species diversity; and the density of nectar-feeding diurnal bird species. Two scales of influence were considered: the stand-scale (2 ha), and a series of radial landscape extents (500 m - 2 km; 78 - 1250 ha) surrounding each fauna transect. For all biota, stand-scale structural and compositional attributes were found to be more influential than landscape metrics. For the Yellow-bellied Glider, the proportion of trace habitats with a residual element of old spotted-gum/ironbark eucalypt trees was a significant landscape metric at the 2 km landscape extent. This is a measure of habitat loss rather than habitat fragmentation. For the diversity of nectar and sap feeding glider species, the proportion of trace habitats with a high coefficient of variation in patch size at the 750 m extent was a significant landscape metric. None of the landscape metrics tested was important for diurnal forest birds. We conclude that no single landscape metric adequately captures the response of the region's forest biota per se. This poses a major challenge to regional reporting of Montreal Process Indicator 1.1e, fragmentation of forest types.

Late Devonian acanthodians from Colombia

Acanthodian remains occur in micaceous siltstone lenses (presumed to have been deposited during a marine incursion) in the Cuche Formation (?Frasnian) of northeast Colombia. The acanthodians are represented by patches of scales from climatiidid Nostolepis sp. cf. N. gatijensis and a fin spine and scales from a new diplacanthid. Type material of N. gaujensis is from the Frasnian Sventoji regional stage in the Baltic, and Nostolepis sp. cf. N. gaujensis has been recorded in the Frasnian of Iran, as well as from Colombia. The new diplacanthid taxon shows affinity to Baltic and Antarctic diplacanthids. The fauna thus shows possible links to both Gondwanan and Euramerican acanthodian assemblages. (C) 2003 Elsevier Science Ltd. All rights reserved.