Assessment of relationships between precipitation and satellite derived vegetation condition in South Australia

The normalised difference vegetation index (NDVI) has evolved as a primary tool for monitoring continental-scale vegetation changes and interpreting the impact of short to long-term climatic events on the biosphere. The objective of this research was to assess the nature of relationships between precipitation and vegetation condition, as measured by the satellite-derived NDVI within South Australia. The correlation, timing and magnitude of the NDVI response to precipitation were examined for different vegetation formations within the State (forest, scrubland, shrubland, woodland and grassland). Results from this study indicate that there are strong relationships between precipitation and NDVI both spatially and temporally within South Australia. Differences in the timing of the NDVI response to precipitation were evident among the five vegetation formations. The most significant relationship between rainfall and NDVI was within the forest formation. Negative correlations between NDVI and precipitation events indicated that vegetation green-up is a result of seasonal patterns in precipitation. Spatial patterns in the average NDVI over the study period closely resembled the boundaries of the five classified vegetation formations within South Australia. Spatial variability within the NDVI data set over the study period differed greatly between and within the vegetation formations examined depending on the location within the state. ACRONYMS AVHRR Advanced Very High Resolution Radiometer ENVSAEnvironments of South Australia EOS Terra-Earth Observing System EVIEnhanced Vegetation Index MODIS Moderate Resolution Imaging Spectro-radiometer MVC Maximum Value Composite NDVINormalised Difference Vegetation Index NIRNear Infra-Red NOAANational Oceanic and Atmospheric Administration SPOT Systeme Pour l’Observation de la Terre. [ABSTRACT FROM AUTHOR]

The programming language python in earth system simulations

-scale vary from a planetary scale and million years for convection problems to 100km and 10 years for fault systems simulations. Various techniques are in use to deal with the time dependency (e.g. Crank-Nicholson), with the non-linearity (e.g. Newton-Raphson) and weakly coupled equations (e.g. non-linear Gauss-Seidel). Besides these high-level solution algorithms discretization methods (e.g. finite element method (FEM), boundary element method (BEM)) are used to deal with spatial derivatives. Typically, large-scale, three dimensional meshes are required to resolve geometrical complexity (e.g. in the case of fault systems) or features in the solution (e.g. in mantel convection simulations). The modelling environment escript allows the rapid implementation of new physics as required for the development of simulation codes in earth sciences. Its main object is to provide a programming language, where the user can define new models and rapidly develop high-level solution algorithms. The current implementation is linked with the finite element package finley as a PDE solver. However, the design is open and other discretization technologies such as finite differences and boundary element methods could be included. escript is implemented as an extension of the interactive programming environment python (see www.python.org). Key concepts introduced are Data objects, which are holding values on nodes or elements of the finite element mesh, and linearPDE objects, which are defining linear partial differential equations to be solved by the underlying discretization technology. In this paper we will show the basic concepts of escript and will show how escript is used to implement a simulation code for interacting fault systems. We will show some results of large-scale, parallel simulations on an SGI Altix system. Acknowledgements: Project work is supported by Australian Commonwealth Government through the Australian Computational Earth Systems Simulator Major National Research Facility, Queensland State Government Smart State Research Facility Fund, The University of Queensland and SGI.

Bajo de la Alumbrera copper-gold deposit: Stable isotope evidence for a porphyry-related hydrothermal system dominated by magmatic aqueous fluids

Alteration zones at the gold-rich Bajo de la Alumbrera porphyry copper deposit in northwestern Argentina are centered on several porphyritic intrusions. They are zoned from a central copper-iron sulfide and gold-mineralized potassic (biotite-K-feldspar +/- quartz) core outward to propylitic (chlorite-illite-epidote-calcite) assemblages. A mineralized intermediate argillic alteration assemblage (chlorite-illite +/- pyrite) has overprinted the potassic alteration zone across the top and sides of the deposit and is itself zoned outward into phyllic (quartzinuscovite-illite +/- pyrite) alteration. This study contributes new data to previously reported delta(18)O and delta D compositions of fluids responsible for the alteration at Bajo de la Alumbrera, and the data are used to infer likely ore-forming processes. Measured and calculated delta(18)O and delta D values of fluids (+8.3 to +10.2 and -33 to -81 parts per thousand, respectively) confirm a primary magmatic origin for the earliest potassic alteration phase. Lower temperature potassic alteration formed from magmatic fluids with lower delta D values (down to -123 parts per thousand). These depleted compositions are distinct from meteoric water and consistent with degassing and volatile exsolution of magmatic fluids derived from an underlying magma. Variability in the calculated composition of fluid associated with potassic alteration is explained in terms of phase separation (or boiling). if copper-iron sulfide deposition occurred during cooling (as proposed elsewhere), this cooling was largely a result of phase separation. Magmatic water was directly involved in the formation of overprinting intermediate argillic alteration assemblages at Bajo de la Alumbrera. Calculated delta(18)O and delta D values of fluids associated with this alteration range from +4.8 to +8.1 and -31 to -71 per mil, respectively Compositions determined for fluids associated with phyllic alteration (-0.8 to +10.2 and -31 to -119 parts per thousand) overlap with the values determined for the intermediate argillic alteration. We infer that phyllic alteration assemblages developed during two stages; the first was a high-temperature (400 degrees-300 degrees C) stage with D-depleted water (delta D = -66 to -119 parts per thousand). This compositional range may have resulted from magma degassing and/or the injection of new magmatic water into a compositionally evolved hydrothermal system. The isotopic variations also can be explained by increased fluid-rock interaction. The second stage of phyllic alteration occurred at a lower temperature (similar to 200 degrees C), and variations in the modeled isotopic compositions imply mixing of magmatic and meteoric waters. Ore deposition that occurred late in the evolution of the hydrothermal system was probably associated with further cooling of the magmatic fluid, in part caused by fluid-rock interaction and phase separation. Changing pH and/or oxygen fuoracity may have caused additional ore deposition. The ingress of meteoric water appears to postdate the bulk of mineralization and occurred as the system at Bajo de la Alumbrera waned.

Temporal variation in habitat use by nekton in a subtropical estuarine system

Utilisation by fish of different estuarine habitats is known to vary at many different temporal scales, however no study to date has examined how utilisation varies at all the relevant times scales simultaneously. Here, we compare the utilisation by fish of sandy, intertidal foreshore habitats in a subtropical estuary at four temporal scales: between major spawning periods (spring/ summer and winter), among months within spawning periods, between the full and new moon each month, and between night and day within those lunar phases. Comparisons of assemblage composition, abundance of individuals and of fish in seven different,ecological guilds' were used to identify the temporal scales at which fish varied their use of unvegetated sandy habitats in the lower Noosa Estuary, Queensland, Australia. Fish assemblages were sampled with a seine net at three different regions. The most numerically dominant species caught were southern herring (Herklotsichthys castelnaui: Clupeidae), sand whiting (Sillago ciliata: Sillaginidae), weeping toadfish (Torquigener pleurogramma: Tetraodomidae), and silver biddy (Gerres subfasciatus: Gerreidae). Considerable variation at a range of temporal scales from short term (day versus night) to longer term (spawning periods) was detected for all but one of the variables examined. The clearest patterns were observed for diurnal effects, where generally abundance was greater at night than during the day. There were also strong lunar effects, although there were no consistent patterns between full moon and new moon periods. Significant differences among months within spawning periods were more common than differences between the actual spawning periods. The results clearly indicate that utilisation of sandy, unvegetated estuarine habitats is very dynamic and highly variable in space and time. (c) 2006 Elsevier B.V. All rights reserved.

Aquatic geochemistry of the rare earth elements and yttrium in the Pioneer River catchment, Australia

The rare earth elements are strong provenance indicators in geological materials, yet the potential for tracing provinciality in surface freshwater samples has not been adequately tested. Rare earth element and yttrium concentrations were measured at 33 locations in the Pioneer River catchment, Mackay, central Queensland, Australia. The rare earth element patterns were compared on the basis of geological, topographical and land-use features in order to investigate the provenancing potential of these elements in a small freshwater system. The rare earth element patterns of streams draining single lithological units with minor land modification show strongly coherent normalised behaviour, with a loss of coherence in agricultural locations. Evidence is reported for an anthropogenic Gd anomaly that may provide a useful hydrological tracer in this region since the introduction of magnetic resonance imaging in 2003. Several samples display a superchondritic Y/Ho mass ratio (up to 44), which is not explainable within the constraints imposed by local geology. Instead, it is suggested that the additional Y is derived from a marine source, specifically marine phosphorites, which are a typical source of fertiliser phosphorus. The data indicate that, under some circumstances, scaled and normalised freshwater rare earth patterns behave conservatively.

Construction of an Intraplate Fault System Model of South Australia, and Simulation Tool for the iSERVO Institute Seed Project

To foster ongoing international cooperation beyond ACES (APEC Cooperation for Earthquake Simulation) on the simulation of solid earth phenomena, agreement was reached to work towards establishment of a frontier international research institute for simulating the solid earth: iSERVO = International Solid Earth Research Virtual Observatory institute (http://www.iservo.edu.au). This paper outlines a key Australian contribution towards the iSERVO institute seed project, this is the construction of: (1) a typical intraplate fault system model using practical fault system data of South Australia (i.e., SA interacting fault model), which includes data management and editing, geometrical modeling and mesh generation; and (2) a finite-element based software tool, which is built on our long-term and ongoing effort to develop the R-minimum strategy based finite-element computational algorithm and software tool for modelling three-dimensional nonlinear frictional contact behavior between multiple deformable bodies with the arbitrarily-shaped contact element strategy. A numerical simulation of the SA fault system is carried out using this software tool to demonstrate its capability and our efforts towards seeding the iSERVO Institute.