Basin evolution of the Vlaming Sub-basin and the deep-water Mentelle Basin, both located offshore on the southwest Australian continental margin, were investigated using 2D and 3D petroleum system modelling. Compositional kinetics, determined on the main source sequences, were used to predict timing of hydrocarbon generation and migration as well as GOR evolution and phase behaviour in our 2D and 3D basin models. The main phase of petroleum generation in the Vlaming Sub-basin occurred at 150 Ma and ceased during following inversion and erosion episodes. Only areas which observed later burial have generated additional hydrocarbons during the Tertiary and up to present day. The modelling results indicate the likely generation and trapping of light oils for the Jurassic intervals for a variety of structural traps. It is these areas which are of greatest interest from an exploration point of view. The 2D numerical simulations in the Mentelle Basin indicate the presence of active hydrocarbon generating kitchen areas. Burial histories and generalized petroleum evolutionary histories are investigated.

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.

Integration of disparate sets of geophysical data, such as Bouguer gravity, magnetotelluric and seismic travel time data for a robust interpretation of architectural settings of the subsurface is carried out. At the outset, the layered 2D model space is appropriately gridded. A spline node layer boundary parameterisation with a sigmoid basis function is used to relate local 1D layered model parameterisation to the 2D model space. Joint 1D inversion of seismic travel time and magnetotelluric data is carried out at the spline nodes using empirical relationships between seismic velocity and resistivity. The two objective functions corresponding to each of the input data types are combined through a weighting factor, the appropriate value of which is determined using the L-curve technique. The particle swarm optimization scheme is used as a robust optimiser for the layer depths and property values. The inverted velocity model is transformed to a density model using a second empirical rock property relationship. A 2D inversion of Bouguer gravity data is then carried out producing adjustments to the depths for the layer boundaries. This completes the initialisation phase of the procedure. A second iterative phase during which only the depths to the layer boundaries are modified is carried out to build a coherent model which is consistent with all three kinds of data. This involves re-inverting jointly the seismic travel time and magnetotelluric data, where parameters corresponding to the rock properties are kept unaltered, but the depth to the layer interfaces are updated. The method is trialled with a synthetic situation and is implemented to interpret the architectural settings of newly discovered Millungera basin of North Queensland, Australia.

To perform a realistic 3D inversion of gravity data covering a significant proportion of the surface of the Earth, it is necessary to take into account the curvature of the Earth. We have developed an algorithm for inverting gravity data in spherical coordinates and have demonstrated this using data covering the continental mass of Australia and surrounding ocean areas. The density structures evident in the crust and uppermost mantle of the resultant 3D inversion model are in broad agreement with knowledge of the geological features for the region and with variations in the depth to the Moho that are present in the AusMoho model.

The Collaborative Research Centre for Greenhouse Gas Technologies (CO2CRC) Program 3.2 Risk Assessment is working toward a risk assessment procedure that integrates risk across the complete CCS system and can be used to meet the needs of a range of stakeholders. Any particular CCS project will hold the interest of multiple stakeholders who will have varied interests in the type of information and in the level of detail they require. It is unlikely that any single risk assessment tool will be able to provide the full range of outputs required to meet the needs of regulators, the general public and project managers; however, in many cases the data and structure behind the outputs will be the same. In using a suite of tools, a well designed procedure will optimize the interaction between the scientists, engineers and other experts contributing to the assessment and will allow for the required information to be presented in a manner appropriate for each stakeholder. Discussions of risk in CCS, even amongst the risk assessment community, often become confused because of the differing emphases on what the risks of interest are. A key question that must be addressed is: 'What questions is the risk analysis trying to answer?' Ultimately, this comes down to the stakeholders, whose interests can be broken into four target questions: - Which part of the capture-transport-storage CCS system? - Which timeline? (project planning, project lifespan, post closure, 1,000 years, etc) - Which risk aspect? (technical, regulatory, economic, public acceptance, or heath safety and environment) - Which risk metric? (Dollars, CO2 lost, dollars/tonne CO2 avoided, etc.) Once the responses to these questions are understood a procedure and suite of tools can be selected that adequately addresses the questions. The key components of the CO2CRC procedure we describe here are: etc

The Fitzoy Estuary is one of several macrotidal estuaries in tropical northern Australia that face ecological change due to agricultural activities in their catchments. The biochemical functioning of such macrotidal estuaries is not well understood in Australia, and there is a pressing need to identify sediment, nutrient and agrochemical pathways, sinks and accumulation rates in these extremely dynamic environments. This is particularly the case in coastal northern Queensland because the impact of water quality changes in rivers resulting from vegetation clearing, changes in land-use and modern agricultural practices are the single greatest threat to the Great Barrier Reef Marine Park. This report includes: 1 Aims and Research questions 2 Study Area 3 Climate and Hydrology 4 Geology 5 Vegetation and land use 6 Methods 7 Sampling strategy 8 Water column observations and samples 9 Bottom sediment properties 10 Core and bottle incubations 11 Data analysis 12 Results 13 Discussion 14 The roll of Keppel Bay in accumulating and redirecting sediment and nutrients from the catchment 15 Sediment biogeochemistry 16 Links between primary production, biogeochemistry and sediment dynamics: A preliminary zonation for Keppel Bay 17 Conclusions

Despite growing concerns about potential enhancement of global warming and slope failure by methane produced by gas hydrate dissociation, much uncertainty surrounds estimates of gas hydrate reservoir sizes, as well as methane fluxes and oxidation rates at the sea floor. For cold seep sediments of the eastern Mediterranean Sea, depth-dependent methane concentrations and rates of anaerobic oxidation of methane (AOM) are constrained by modeling the measured pore-water sulfate profile. The calculated dissolved methane distribution and flux are sensitive to the advective flow velocity, which is estimated from the depth distributions of conservative pore-water constituents (Na, B). Near-complete anaerobic oxidation of the upward methane flux is supported by the depth distributions of indicative biomarkers, and the carbon isotopic compositions of organic matter and dissolved inorganic carbon. Pore-water and solid-phase data are consistent with a narrow depth interval of AOM, 14-18 cm below the sediment-water interface. Based on an isotopic mass balance, the biomass of the microbial population carrying out oxidation of methane coupled to sulfate reduction at the given methane flux represents about 20% of the total organic carbon, which is a significant pool of in situ formed organic matter. Model results indicate that the asymptotic methane concentration is reached a few meters below the sediment surface. The predicted asymptotic concentration is close to the in situ saturation value with respect to gas hydrate, suggesting that the rate of shallow gas hydrate formation is controlled by the ascending methane flux. The proposed model approach can be used to predict the formation of gas hydrate, and to quantify methane fluxes plus transformation rates in surface sediments where fluid advection is an important transport mechanism.

The Paterson National Geoscience Agreement project is using a number of tools to better understand the time-space evolution of the northwest Paterson Orogen in Western Australia. One of these tools, 3D Geomodeller, is an emerging technology that constructs three-dimensional (3D) volumetric models based on a range of geological information. The Paterson project is using 3D Geomodeller to build geologically-constrained 3D models for the northwest Paterson Orogen. This report documents the model building capability and benefits of 3D Geomodeller and highlights some of the geological insights gained from the model building exercise. The principal benefit of 3D Geomodeller is that it provides geoscientists with a rapid tool for testing multiple working hypotheses. The Cottesloe Syncline district was selected as the focus for a trial of the 3D Geomodeller software. The 3D model was built by members of the Paterson Project, as well as model building specialists within Geoscience Australia. The resultant Cottesloe Syncline model including two dimensional sections, maps and images was exported from 3D GeoModeller and transformed into a Virtual Reality Modelling Language (VRML), enabling a wide audience to view the model using readily available software.

Geoscience Australia's Risk & Impact Analysis Group has developed a statistical model of wind hazard utilising the Generalised Pareto Distribution (GPD). The model calculates the return period of severe winds based on daily maximum wind gust observations. The model utilises an automated procedure to partition the data into the hazard constituents (thunderstorms, synoptic winds, tornadoes, etc) based on the World Meteorological Observation Codes 3-hourly coded observations. This observational data set records the archived and present weather at the station site. The model fits the GPD to the station data (daily maximum wind gust) by automating the selection of the appropriate threshold above which data is included in the extreme value distribution. This threshold <em>u</em> is selected as the maximum of all feasible return period values obtained by fitting the GPD. Published comparative findings, including same region results, demonstrate the model can produce similar results in a more efficient, fully computational way. Confidence intervals for return periods are calculated automatically to allow wind analysts to distinguish regions of greater reliability.

For many basins along the western Australian margin, knowledge of basement and crustal structure is limited, yet both play an important role in controlling basin evolution. To provide new insight into these fundamental features of a continental margin, we present the results of process-oriented gravity modelling along a NW-SE profile across the Browse Basin through the Brecknock field. Process-oriented gravity modelling is a method that considers the rifting, sedimentation and magmatism that led to the present-day gravity field. By backstripping the sediment load under different isostatic assumptions (i.e. range of flexural rigidities), the crustal structure associated with rifting can be inferred. Combining the gravity anomalies caused by rifting and sedimentation and comparing them to observed gravity provides insight into the presence of magmatic underplating, the location of the continent-ocean boundary and the thermal history of a margin. For an effective elastic thickness of 25 km, backstripping syn- and post-rift sediments (Jurassic and younger) along the Browse Basin profile suggests moderate Jurassic stretching (beta-1-2) and shows that rifting and sedimentation generally explain the observed free-air gravity signature. The gravity fit is reasonable for most of the Scott Plateau and Caswell Sub-basin, but over the Leveque Shelf and Wilson Spur, predicted gravity is less than observed and predicted Moho is also shallower than suggested by seismic refraction data. These misfits suggest the presence of magmatic underplating beneath the Leveque Shelf and outermost parts of the basin, an inference that has mixed support from refraction and crustal-scale seismic reflection data.