Results of wave modelling for Moreton Bay, Southeast Queensland, and a reef lagoon off Lizard Island, Great Barrier Reef, Australia

The distribution, abundance, behaviour, and morphology of marine species is affected by spatial variability in the wave environment. Maps of wave metrics (e.g. significant wave height Hs, peak energy wave period Tp, and benthic wave orbital velocity URMS) are therefore useful for predictive ecological models of marine species and ecosystems. A number of techniques are available to generate maps of wave metrics, with varying levels of complexity in terms of input data requirements, operator knowledge, and computation time. Relatively simple "fetch-based" models are generated using geographic information system (GIS) layers of bathymetry and dominant wind speed and direction. More complex, but computationally expensive, "process-based" models are generated using numerical models such as the Simulating Waves Nearshore (SWAN) model. We generated maps of wave metrics based on both fetch-based and process-based models and asked whether predictive performance in models of benthic marine habitats differed. Predictive models of seagrass distribution for Moreton Bay, Southeast Queensland, and Lizard Island, Great Barrier Reef, Australia, were generated using maps based on each type of wave model. For Lizard Island, performance of the process-based wave maps was significantly better for describing the presence of seagrass, based on Hs, Tp, and URMS. Conversely, for the predictive model of seagrass in Moreton Bay, based on benthic light availability and Hs, there was no difference in performance using the maps of the different wave metrics. For predictive models where wave metrics are the dominant factor determining ecological processes it is recommended that process-based models be used. Our results suggest that for models where wave metrics provide secondarily useful information, either fetch- or process-based models may be equally useful.

(Table 1) Abundance of tunicate spicules along the Great Barrier Reef-Queensland Plateau transect

The abundance patterns of tunicate spicules are documented for the Pliocene-Pleistocene sediments at seven sites along the Great Barrier Reef-Queensland Plateau transect. The spatial distribution pattern indicates that tunicate spicules were limited to waters shallower than 900 m. The occurrences of tunicate spicules at Sites 822 and 823 that are deeper than 900 m are ascribed to downslope transport, and their distribution patterns can be used to monitor downslope transport processes. The first common occurrence of tunicate spicules at Sites 822 and 823 around 1.6 Ma may suggest the initiation of the central Great Barrier Reef at this time. The morphology of tunicate spicules varies greatly and appears to be gradational among different forms. Older tunicate assemblages are less diverse than those in younger sediments, presumably because of diagenesis. Tunicate spicules do not appear to be a promising biostratigraphic tool for the Pliocene-Pleistocene.

Sedimentation across the Great Barrier Reef

The continental margin off northeast Australia, comprising the Great Barrier Reef (GBR) platform and Queensland Trough, is the largest tropical mixed siliciclastic/carbonate depositional system in existence. We describe a suite of 35 piston cores and two Ocean Drilling Program (ODP) sites from a 130*240 km rectangular area of the Queensland Trough, the slope and basin setting east of the central GBR platform. Oxygen isotope records, physical property (magnetic susceptibility and greyscale) logs, analyses of bulk carbonate content and radiocarbon ages at these locations are used to construct a high resolution stratigraphy. This information is used to quantify mass accumulation rates (MARs) for siliciclastic and carbonate sediments accumulating in the Queensland Trough over the last 31,000 years. For the slope, highest MARs of siliciclastic sediment occur during transgression (1.0 Million Tonnes per year; MT/yr), and lowest MARs of siliciclastic (<0.1 MT/yr) and carbonate (0.2 MT/yr) sediment occur during sea level lowstand. Carbonate MARs are similar to siliciclastic MARs for transgression and highstand (1.1-1.4 MT/yr). In contrast, for the basin, MARs of siliciclastic (0-0.1 MT/yr) and carbonate sediment (0.2-0.4 MT/yr) are continuously low, and within a factor of two, for lowstand, transgression, and highstand. Generic models for carbonate margins predict that maximum and minimum carbonate MARs on the slope will occur during highstand and lowstand, respectively. Conversely, most models for siliciclastic margins suggest maximum and minimum siliciclastic MARs will occur during lowstand and transgression, respectively. Although carbonate MARs in the Queensland Trough are similar to those predicted for carbonate depositional systems, siliciclastic MARs are the opposite. Given uniform siliciclastic MARs in the basin through time, we conclude that terrigenous material is stored on the shelf during sea level lowstand, and released to the slope during transgression as wave driven currents transport shelf sediment offshore.

A 20 m spatial resolution seamless multisource Digital Elevation/Depth Model for Lizard Island, northern Great Barrier Reef

Hydrographers have traditionally referred to the nearshore area as the "white ribbon" area due to the challenges associated with the collection of elevation data in this highly dynamic transitional zone between terrestrial and marine environments. Accordingly, available information in this zone is typically characterised by a range of datasets from disparate sources. In this paper we propose a framework to 'fill' the white ribbon area of a coral reef system by integrating multiple elevation and bathymetric datasets acquired by a suite of remote-sensing technologies into a seamless digital elevation model (DEM). A range of datasets are integrated, including field-collected GPS elevation points, terrestrial and bathymetric LiDAR, single and multibeam bathymetry, nautical chart depths and empirically derived bathymetry estimations from optical remote sensing imagery. The proposed framework ranks data reliability internally, thereby avoiding the requirements to quantify absolute error and results in a high resolution, seamless product. Nested within this approach is an effective spatially explicit technique for improving the accuracy of bathymetry estimates derived empirically from optical satellite imagery through modelling the spatial structure of residuals. The approach was applied to data collected on and around Lizard Island in northern Australia. Collectively, the framework holds promise for filling the white ribbon zone in coastal areas characterised by similar data availability scenarios. The seamless DEM is referenced to the horizontal coordinate system MGA Zone 55 - GDA 1994, mean sea level (MSL) vertical datum and has a spatial resolution of 20 m.

Benthic and substrate cover data derived from photo-transect surveys in Lizard Island, Great Barrier Reef conducted on October 3-7, 2012

Underwater georeferenced photo-transect surveys were conducted on October 3-7, 2012 at various sections of the reef and lagoon at Lizard Island, Great Barrier Reef. For this survey a snorkeler swam while taking photos of the benthos at a set distance from the benthos using a standard digital camera and towing a GPS in a surface float which logged the track every five seconds. A Canon G12 digital camera was placed in a Canon underwater housing and photos were taken at 1 m height above the benthos. Horizontal distance between photos was estimated by three fin kicks of the survey snorkeler, which corresponded to a surface distance of approximately 2.0 - 4.0 m. The GPS was placed in a dry bag and logged the position at the surface while being towed by the photographer (Roelfsema, 2009). A total of 1,265 benthic photos were taken. Approximation of coordinates of each benthic photo was conducted based on the photo timestamp and GPS coordinate time stamp, using GPS Photo Link Software (www.geospatialexperts.com). Coordinates of each photo were interpolated by finding the GPS coordinates that were logged at a set time before and after the photo was captured. Benthic or substrate cover data was derived from each photo by randomly placing 24 points over each image using the Coral Point Count for Microsoft Excel program (Kohler and Gill, 2006). Each point was then assigned to 1 of 79 cover types, which represented the benthic feature beneath it. Benthic cover composition summary of each photo scores was generated automatically using CPCE program. The resulting benthic cover data of each photo was linked to GPS coordinates, saved as an ArcMap point shapefile, and projected to Universal Transverse Mercator WGS84 Zone 55 South.

Seawater carbonate chemistry and calcification during tropical reef studies at Younge reef (Great Barrier Reef), 1996

Community metabolism was investigated using a Lagrangian flow respirometry technique on 2 reef flats at Moorea (French Polynesia) during austral winter and Yonge Reef (Great Barrier Reef) during austral summer. The data were used to estimate related air-sea CO2 disequilibrium. A sine function did not satisfactorily model the diel light curves and overestimated the metabolic parameters. The ranges of community gross primary production and respiration (Pg and R; 9 to 15 g C m-2 d-1) were within the range previously reported for reef flats, and community net calcification (G; 19 to 25 g CaCO3 m-2 d-1) was higher than the 'standard' range. The molar ratio of organic to inorganic carbon uptake was 6:1 for both sites. The reef flat at Moorea displayed a higher rate of organic production and a lower rate of calcification compared to previous measurements carried out during austral summer. The approximate uncertainty of the daily metabolic parameters was estimated using a procedure based on a Monte Carlo simulation. The standard errors of Pg,R and Pg/R expressed as a percentage of the mean are lower than 3% but are comparatively larger for E, the excess production (6 to 78%). The daily air-sea CO2 flux (FCO2) was positive throughout the field experiments, indicating that the reef flats at Moorea and Yonge Reef released CO2 to the atmosphere at the time of measurement. FCO2 decreased as a function of increasing daily irradiance.

Great Barrier Reef coral element/Ca and stable isotope data and U-Th ages from IODP Expedition 325

Tropical south-western Pacific temperatures are of vital importance to the Great Barrier Reef (GBR), but the role of sea surface temperatures (SSTs) in the growth of the GBR since the Last Glacial Maximum remains largely unknown. Here we present records of Sr/Ca and d18O for Last Glacial Maximum and deglacial corals that show a considerably steeper meridional SST gradient than the present day in the central GBR. We find a 1-2 °C larger temperature decrease between 17° and 20°S about 20,000 to 13,000 years ago. The result is best explained by the northward expansion of cooler subtropical waters due to a weakening of the South Pacific gyre and East Australian Current. Our findings indicate that the GBR experienced substantial meridional temperature change during the last deglaciation, and serve to explain anomalous deglacial drying of northeastern Australia. Overall, the GBR developed through significant SST change and may be more resilient than previously thought.

Seagrass and associated benthic community data derived from field surveys at Low Isles, Great Barrier Reef, conducted July-August, 1997

The distribution of seagrass and associated benthic communities on the reef and lagoon of Low Isles, Great Barrier Reef, was mapped between the 29 July and 29 August 1997. For this survey, observers walked or free-dived at survey points positioned approximately 50 m apart along a series of transects. Visual estimates of above-ground seagrass biomass and % cover of each benthos and substrate type were recorded at each survey point. A differential handheld global positioning system (GPS) was used to locate each survey point (accuracy ±3m). A total of 349 benthic survey points were examined. To assist with mapping meadow/habitat type boundaries, an additional 177 field points were assessed and a georeferenced 1:12,000 aerial photograph (26th August 1997) was used as a secondary source of information. Bathymetric data (elevation below Mean Sea Level) measured at each point assessed and from Ellison (1997) supplemented information used to determine boundaries, particularly in the subtidal lagoon. 127.8 ±29.6 hectares was mapped. Seagrass and associated benthic community data was derived by haphazardly placing 3 quadrats (0.25m**2) at each survey point. Seagrass above ground biomass (standing crop, grams dry weight (g DW m**-2)) was determined within each quadrat using a non-destructive visual estimates of biomass technique and the seagrass species present identified. In addition, the cover of all benthos was measured within each of the 3 quadrats using a systematic 5 point method. For each quadrat, frequency of occurrence for each benthic category was converted to a percentage of the total number of points (5 per quadrat). Data are presented as the average of the 3 quadrats at each point. Polygons of discrete seagrass meadow/habitat type boundaries were created using the on-screen digitising functions of ArcGIS (ESRI Inc.), differentiated on the basis of colour, texture, and the geomorphic and geographical context. The resulting seagrass and benthic cover data of each survey point and for each seagrass meadow/habitat type was linked to GPS coordinates, saved as an ArcMap point and polygon shapefile, respectively, and projected to Universal Transverse Mercator WGS84 Zone 55 South.

(Table 1) Bulk carbonate content of ODP Hole 194-1198A sediments, east of the Great Barrier Reef

Upper Quaternary sediment sequences east of the Great Barrier Reef are characterized by alternating siliciclastic- and carbonate-rich horizons caused by changes in the input of various sedimentary components and reflected in cores by variations in bulk carbonate content. A total of 153 measurements of bulk carbonate content were determined using the carbonate-bomb technique for late Pleistocene sediments between 0 and 23.69 meters below sea floor (mbsf) in Ocean Drilling Program Hole 1198A. Average sample resolution was 15 cm and multiple analyses were performed on each sample. Bulk carbonate content ranges from a maximum of 94 wt% at 13.63 mbsf to a minimum of 73 wt% at 14.54 mbsf. Five cyclic trends are observed that may relate to five major glacial events during the last 500 k.y. of the Quaternary.

Stable oxygen isotope record of Globigerinoides ruber of ODP Site 133-820 off the Great Barriere Reef

Oxygen-isotope ratio measurements are presented for the planktonic species Globigerinoides ruber collected from shallow-water, upper-slope sediments from Holes 820A and 820B in 280 m of water, on the seaward edge of the Great Barrier Reef. Correlation of the Site 820 isotope curve with deep-sea reference curves of the Pacific Ocean (Core V28-238, Hole 677A, Hole 607A) permits the definition of isotope stages 1 to 19 in the top 145 m of Holes 820A and 820B. However, paleontological data indicate that stages 4 and 7 might be missing and that two hiatuses occur at a depth of 8.05 to 12.1 and 34.55 to 35.8 mbsf. Using deep-sea Hole 677A as a reference for ice-volume variations, we determine the difference in isotopic signature between it and Site 820. We propose that this difference is a regional signal representing a progressive 4°C increase in surface-water temperature at Site 820. The proposed temperature change was initiated at about 400 k.y. and corresponds to a change from high-to-low frequency variations in Pleistocene isotope signals. We postulate that these changes may have catalyzed the growth of the Great Barrier Reef. The shift also coincides with changes in seismic character and some physical and chemical sediment characteristics.