In ecology, a common form of statistical analysis relates a biological variable to variables that delineate the physical environment, typically by fitting a regression model or one of its extensions. Unfortunately, the biological data and the physical data are frequently obtained from eparate sources of data. In such cases there is no guarantee that the biological and physical data are co-located and the regression model cannot be used. A common and pragmatic solution is to predict the physical variables at the locations of the biological variables and then to use the predictions as if they were observations.We show that this procedure can cause potentially misleading inferences and we use generalized linear models as an example. We propose a Berkson error model which overcomes the limitations. The differences between using predicted covariates and the Berkson error model are illustrated by using data from the marine environment, and a simulation study based on these data.

The characterisation of benthic habitats based on their abiotic (physical and chemical) attributes remains poorly defined in the marine environment, but is becoming increasingly central in the development of marine management plans in Australia and elsewhere in the world. The current study tested this link between physical and biological datasets for the southern Gulf of Carpentaria, Australia. The results presented were based on a range of physical factors, including the sediment composition (grain size and carbonate content), sediment mobility, water depth and organic carbon flux, and their relationship to the distribution and diversity of benthic macrofauna was tested. The results reveal the importance of process-based indices, such as sediment mobility, in addition to other environmental factors in defining the distribution of the benthic macrofauna. The distribution of the benthic macrofauna changes gradationally across the south-eastern Gulf, associated with changes in the per cent mud and gravel, the seabed exposure and the water depth. Patterns of diversity also reveal the importance of physical processes such as sediment mobility in defining benthic habitats. The species' environment relationships observed at the small scale of the current study are consistent with broader associations observed for other organisms within the Gulf.

An integrated analysis of both airborne and field short-wave infrared hyperspectral measurements was used in conjunction with conventional field mapping techniques to map hydrothermal alteration in the central portion of the Mount Painter Inlier in the Flinders Ranges, South Australia. The airborne hyperspectral data show the spatial distribution of spectrally distinct minerals occurring as primary minerals and as weathering and alteration products. Field spectral measurements, taken with a portable infrared mineral analyzer spectrometer and supported by thin-section analyses, were used to verify the mineral maps and enhance the level of information obtainable from the airborne data. Hydrothermal alteration zones were identified and mapped separately from the background weathering signals. A main zone of alteration, coinciding with the Paralana Fault zone, was recognized, and found to contain kaolinite, muscovite, biotite, and K-feldspar. A small spectral variation associated with a ring-like feature around Mount Painter was tentatively determined to be halloysite and interpreted to represent a separate hydrothermal fluid and fluid source, and probably a separate system. The older parts of the alteration system are tentatively dated as Permo-Carboniferous. The remote sensing of alteration at Mount Painter confirms that hyperspectral imaging techniques can produce accurate mineralogical maps with significant details that can be used to identify and map hydrothermal activity. Application of hyperspectral surveys such as that conducted at Mount Painter would be likely to provide similar detail about putative hydrothermal deposits on Mars.

Effective management of the global ocean requires an inventory of its features and marine living and non-living resources. To help meet this need, a new global seafloor geomorphic features map (GSFM) has been created based on the analysis and interpretation of the Shuttle Radar Topography Mission (SRTM) 30 arc-second (~1 km) bathymetry grid. The new digital GSFM includes 131,190 separate polygons in 29 geomorphic feature categories. We present the first comprehensive identification, enumeration, inventory and quantitative analysis of the ocean's seafloor geomorphic features. The GSFM allows a more accurate assessment of features (proxies for benthic habitats, ecosystems and resources). GIS analysis of the GSFM illustrates that more than 50% of the area of 11 feature categories are located beyond the area of national jurisdiction, and less than 10% of 21 feature categories are protected in marine reserves globally, including shelf valleys, submarine canyons, mid-ocean spreading ridges and rift valleys.

Sediment-hosted Pb-Zn (SH Pb-Zn) deposits can be divided into two sub-types: 'clastic-dominated lead-zinc' (CD Pb-Zn) ores hosted in shale, sandstone, siltstone, mixed clastic or as carbonate replacement within a clastic dominated sedimentary sequence and Mississippi Valley-type (MVT Pb-Zn) ores that occurs in platform carbonate sequences, typically in a passive margin tectonic setting. The emergence of CD and MVT deposits in the rock record between 2.02 Ga, the age of the earliest known deposit of these ores, and 1.85-1.58 Ga, a major period of CD Pb-Zn mineralization in Australia and India, corresponds to a time after the 'Great Oxygenation Event' (GOE) ca 2.4 to 1.8 Ga. Contributing to the blooming of CD deposits at ca 1.85-1.58 Ga was the following: a) enhanced oxidation of sulfides in the Earth's crust that provided sulfate and lead and zinc to the hydrosphere; b) development of major redox and compositional gradients in the oceans; c) first formation of significant sulfate-bearing evaporites; d) formation of red beds and oxidized aquifers: e) evolution of sulfate-reducing bacteria; and f) the formation of large and long-lived basins on stable cratons. A significant limitation imposed on interpreting the secular distribution of the deposits is that presently, there is no way to quantitatively evaluate the removal of deposits from the rock record through tectonic recycling. Considering that most of the sedimentary rock record has been recycled, probably most SH Pb-Zn deposits have also been destroyed by subduction and erosion or modified by metamorphism and tectonism so that they are no longer recognizable. Thus, the uneven secular distribution of SH Pb-Zn deposits reflects the genesis of these deposits, linked to Earth's evolving tectonic and geochemical systems, as well a record severely censored by an unknown amount of recycling of the sedimentary rock record.

The middle to lower Jurassic sequence in Australia's Surat Basin has been identified as a potential reservoir system for geological CO2 storage. The sequence comprises three major formations with distinctly different mineral compositions, and generally low salinity formation water (TDS<3000 mg/L). Differing geochemical responses between the formations are expected during geological CO2 storage. However, given the prevailing use of saline reservoirs in CCS projects elsewhere, limited data are available on CO2-water-rock dynamics during CO2 storage in such low-salinity formations. Here, a combined batch experiment and numerical modelling approach is used to characterise reaction pathways and to identify geochemical tracers of CO2 migration in the low-salinity Jurassic sandstone units. Reservoir system mineralogy was characterized for 66 core samples from stratigraphic well GSQ Chinchilla 4, and six representative samples were reacted with synthetic formation water and high-purity CO2 for up to 27 days at a range of pressures. Low formation water salinity, temperature, and mineralization yield high solubility trapping capacity (1.18 mol/L at 45°C, 100 bar), while the paucity of divalent cations in groundwater and the silicate reservoir matrix results in very low mineral trapping capacity under storage conditions. Formation water alkalinity buffers pH at elevated CO2 pressures and exerts control on mineral dissolution rates. Non-radiogenic, regional groundwater-like 87Sr/86Sr values (0.7048-0.7066) indicate carbonate and authigenic clay dissolution as the primary reaction pathways regulating solution composition, with limited dissolution of the clastic matrix during the incubations. Several geochemical tracers are mobilised in concentrations greater than found in regional groundwater, most notably cobalt, concentrations of which are significantly elevated regardless of CO2 pressure or sample mineralogy.

The surface sedimentary record from six sediment cores collected from beneath the Amery Ice Shelf, East Antarctica, provides a unique view of the sedimentary and oceanographic processes in this sub-ice shelf setting. The composition and age of the surface sediments indicate spatial variations in ice shelf cavity-ocean interaction, which are consistent with patterns of ocean inflow and outflow modelled and observed beneath the ice shelf. Sediments within 100 km of the ice shelf front (site AM01b) show the greatest open ocean influence with a young surface age and the highest total diatom abundance, compared to older ages and lower diatom abundances at sites deeper in the cavity (AM03 to AM06). The variable marine influence between sites determines the nature of benthic communities, with seabed imagery indicating the existence of sessile suspension feeders in areas of strong marine inflow (site AM01b), while grazers, deposit feeders and a few suspension feeders occur at sites more distal from the shelf calving front where the food supply is lower (sites AM03 and AM04). Understanding the sedimentary and oceanographic processes within the sub-ice shelf environment allows better constraint of interpretations of down core sediment records, an improved understanding of the nature of biological communities in sub-ice shelf environments, and a baseline for determining the sensitivity of the system to any future changes in ocean dynamics.

A 2-D crustal velocity model has been derived from a 1997 364 km north-south wide-angle seismic profile that passed from Ordovician volcanic and volcaniclastic rocks (Molong Volcanic Belt of the Macquarie Arc) in the north, across the Lachlan Transverse Zone into Ordovician turbidites and Early Devonian intrusive granitoids in the south. The Lachlan Transverse Zone is a proposed west-northwest to east-southeast structural feature in the Eastern Lachlan Orogen and is considered to be a possible early lithospheric feature controlling structural evolution in eastern Australia; its true nature, however, is still contentious. The velocity model highlights significant north to south lateral variations in subsurface crustal architecture in the upper and middle crust. In particular, a higher P-wave velocity (6.24-6.32 km/s) layer identified as metamorphosed arc rocks (sensu lato) in the upper crust under the arc at 5-15 km depth is juxtaposed against Ordovician craton-derived turbidites by an inferred south-dipping fault that marks the southern boundary of the Lachlan Transverse Zone. Near-surface P-wave velocities in the Lachlan Transverse Zone are markedly less than those along other parts of the profile and some of these may be attributed to mid-Miocene volcanic centres. In the middle and lower crust there are poorly defined velocity features that we infer to be related to the Lachlan Transverse Zone. The Moho depth increases from 37 km in the north to 47 km in the south, above an underlying upper mantle with a P-wave velocity of 8.19 km/s. Comparison with velocity layers in the Proterozoic Broken Hill Block supports the inferred presence of Cambrian oceanic mafic volcanics (or an accreted mafic volcanic terrane) as substrate to this part of the Eastern Lachlan Orogen. Overall, the seismic data indicate significant differences in crustal architecture between the northern and southern parts of the profile. The crustal-scale P-wave velocity differences are attributed to the different early crustal evolution processes north and south of the Lachlan Transverse Zone.

Benthic habitats on the continental shelf are strongly influenced by exposure to the effects of surface ocean waves, and tidal, wind and density driven ocean currents. These processes combine to induce a combined flow bed shear stress upon the seabed which can mobilise sediments or directly influence organisms disturbing the benthic environment. Output from a suite of numerical models predicting these oceanic processes have been utilised to compute the combined flow bed shear stresses over the entire Australian continental shelf for an 8-year period (March 1997- February 2005 inclusive). To quantify the relative influence of extreme or catastrophic combined flow bed shear stress events and more frequent events of smaller magnitude, three methods of classifying the oceanographic levels of exposure are presented: 1. A spectral regionalisation method, 2. A method based on the shape of the probability distribution function, and 3. A method which assesses the balance between the amount of work a stress does on the seabed, and the frequency with which it occurs. Significant relationships occur between the three regionalisation maps indicating seabed exposure to oceanographic processes and physical sediment properties (mean grain size and bulk carbonate content), and water depth, particularly when distinction is made between regions dominated by high-frequency (diurnal or semi-diurnal) events and low-frequency (synoptic or annual) events. It is concluded that both magnitude and frequency of combined-flow bed shear stresses must be considered when characterising the benthic environment. The regionalisation outputs of the Australian continental shelf presented in this study are expected to be of benefit to quantifying exposure of seabed habitats on the continental shelf to oceanographic processes in future habitat classification schemes for marine planning and policy procedures.

In most circumstances the conventional radioelement ratio method is sufficient for the enhancement of the differences between radioelement concentrations across map areas. However, there are areas where the range of radioelement concentration values are such that the ratio image is dominated by one or other of the radioelements. In this paper we demonstrate that, in some areas, the use of these ratios can be enhanced through suitable normalisation of the radioelement data prior to the calculation of ratios.