Seagrass time-series analysis products, in Moreton Bay, Australia, with links to GeoTIFFs

The spatial and temporal dynamics of seagrasses have been well studied at the leaf to patch scales, however, the link to large spatial extent landscape and population dynamics is still unresolved in seagrass ecology. Traditional remote sensing approaches have lacked the temporal resolution and consistency to appropriately address this issue. This study uses two high temporal resolution time-series of thematic seagrass cover maps to examine the spatial and temporal dynamics of seagrass at both an inter- and intra-annual time scales, one of the first globally to do so at this scale. Previous work by the authors developed an object-based approach to map seagrass cover level distribution from a long term archive of Landsat TM and ETM+ images on the Eastern Banks (~200 km**2), Moreton Bay, Australia. In this work a range of trend and time-series analysis methods are demonstrated for a time-series of 23 annual maps from 1988 to 2010 and a time-series of 16 monthly maps during 2008-2010. Significant new insight was presented regarding the inter- and intra-annual dynamics of seagrass persistence over time, seagrass cover level variability, seagrass cover level trajectory, and change in area of seagrass and cover levels over time. Overall we found that there was no significant decline in total seagrass area on the Eastern Banks, but there was a significant decline in seagrass cover level condition. A case study of two smaller communities within the Eastern Banks that experienced a decline in both overall seagrass area and condition are examined in detail, highlighting possible differences in environmental and process drivers. We demonstrate how trend and time-series analysis enabled seagrass distribution to be appropriately assessed in context of its spatial and temporal history and provides the ability to not only quantify change, but also describe the type of change. We also demonstrate the potential use of time-series analysis products to investigate seagrass growth and decline as well as the processes that drive it. This study demonstrates clear benefits over traditional seagrass mapping and monitoring approaches, and provides a proof of concept for the use of trend and time-series analysis of remotely sensed seagrass products to benefit current endeavours in seagrass ecology.

Grain size analysis on CRP sediment cores from the Ross Sea, Antarctica

The upper 1200 m of pre-Pliocene sediment recovered by Cape Roberts Project (CRP) drilling off the Victoria Land coast of Antarctica between 1997-1999 has been subdivided into 54 unconformity-bound stratigraphic sequences, spanning the period c. 32 to 17 Ma. The sequences are recognised on the basis of the cyclical vertical stacking of their constituent lithofacies, which are enclosed by erosion surfaces produced during the grounding of the advancing ice margin onto the sea floor. Each sequence represents deposition in a range of offshore shelf to coastal glacimarine sedimentary environments during oscillations in the ice margin across the Western Ross Sea shelf, and coeval fluctuations in water depth. This paper applies spectral analysis techniques to depth- and time-series of sediment grain size (500 samples) for intervals of the core with adequate chronological data. Time series analysis of 0.5-l.0m-spaced grainsize data spanning sequences 9-11 (CRP-2/2A) and sequences 1-7 (CRP-3) suggests that the length of individual sequences correspond to Milankovitch frequencies, probably 41 k.y., but possibly as low as 100 k.y. Higher frequency periodic components at 23 k.y. (orbital precession) and 15-10 k.y. (sub-orbital) are recognised at the intrasequence-scale, and may represent climatic cycles akin to the ice rafting episodes described in the North Atlantic Ocean during the Quaternary. The cyclicity recorded by glacimarine sequences in CRP core provides direct evidence from the periphery of Antarctica for orbital oscillations in the size of the Oligocene-Early Miocene East Antarctic Ice Sheet.

Oxygen isotopes, IRD and SST reconstruction for the interval MIS31-19 of IODP Site 306-U1314

Early and Mid-Pleistocene climate, ocean hydrography and ice sheet dynamics have been reconstructed using a high-resolution data set (planktonic and benthic d18O time series, faunal-based sea surface temperature (SST) reconstructions and ice-rafted debris (IRD)) record from a high-deposition-rate sedimentary succession recovered at the Gardar Drift formation in the subpolar North Atlantic (Integrated Ocean Drilling Program Leg 306, Site U1314). Our sedimentary record spans from late in Marine Isotope Stage (MIS) 31 to MIS 19 (1069-779 ka). Different trends of the benthic and planktonic oxygen isotopes, SST and IRD records before and after MIS 25 (~940 ka) evidence the large increase in Northern Hemisphere ice-volume, linked to the cyclicity change from the 41-kyr to the 100-kyr that occurred during the Mid-Pleistocene Transition (MPT). Beside longer glacial-interglacial (G-IG) variability, millennial-scale fluctuations were a pervasive feature across our study. Negative excursions in the benthic d18O time series observed at the times of IRD events may be related to glacio-eustatic changes due to ice sheets retreats and/or to changes in deep hydrography. Time series analysis on surface water proxies (IRD, SST and planktonic d18O) of the interval between MIS 31 to MIS 26 shows that the timing of these millennial-scale climate changes are related to half-precessional (10 kyr) components of the insolation forcing, which are interpreted as cross-equatorial heat transport toward high latitudes during both equinox insolation maxima at the equator.

Geochemistry of black shales from Breitbach Creek

Lower Cretaceous sediments are frequently characterized by a well expressed cyclicity. While the processes influencing environments above the carbonate compensation depth (CCD) are reasonably well understood, almost nothing is known about the deep ocean. Cretaceous sub-CCD sediments from the Tethys and Atlantic Oceans typically show rhythmic black/green shale successions. To gain insight into the nature of these black/green shale cycles, we performed detailed geochemical analyses (X-ray fluorescence, Rock-Eval and reactive iron analysis) on a 3 m long section of latest Aptian age. The major-element distribution of the analyzed shale sequence indicates a periodic change from a high-productivity and well-oxygenated green shale mode to a low-productivity oxygen-deficient black shale mode. It is proposed here that the preservation of organic matter was dependent on the strength of salinity-driven deepwater generation. Furthermore, the data show that the Corg content covaries with changes in the detrital composition. Therefore we hypothesize that Tethyan deepwater circulation was sensitive to changes in the monsoonal system. Time series analysis suggests that these changes are periodic in nature, although we are currently unable to prove that the dominant periodicity is related to the precession component of the Milankovitch frequencies.

Quaternary carbonate and stable isotope record of ODP Holes 133-817A and 133-818B

The Quaternary history of metastable CaCO3 input and preservation within Antarctic Intermediate Water (AAIW) was examined by studying sediments from ODP Holes 818B (745 mbsl) and 817A (1015 mbsl) drilled in the Townsville Trough on the southern slope of the Queensland Plateau. These sites lie within the core of modern AAIW, and near the aragonite saturation depth (~1000 m). Thus, they are well positioned to monitor chemical changes that may have occurred within this watermass during the past 1.6 m.y. The percent of fine aragonite content, percent of fine magnesian calcite content, and percent of whole pteropods (>355 µm) were used to separate the fine aragonite input signal from the CaCO3 preservation signal. Stable d18O and d13C isotopic ratios were determined for the planktonic foraminifer Globigerinoides sacculifer and, in Hole 818B, for the benthic foraminifer Cibicidoides spp. to establish the oxygen isotope stratigraphy and to study the relationship between intermediate and shallow water d13C of Sum CO2 and the relationship between benthic foraminiferal d13C and CaCO3 preservation within intermediate waters of the Townsville Trough. Data were converted from depth to age using oxygen isotope stratigraphy, nannostratigraphy, and foraminiferal biostratigraphy. Several long hiatuses and the absence of magnetostratigraphy did not permit time series analysis. The principal results of the CaCO3 preservation study include the following (1) a general increase in CaCO3 preservation between 0.9 and 1.6 Ma; (2) a CaCO3 dissolution maximum near 0.9 Ma, primarily expressed in the Hole 818B fine aragonite record; (3) an abrupt and permanent increase of fine aragonite content between 0.86 and 0.875 Ma in both Holes 818B and 817A probably reflecting a dramatic increase of fine carbonate sediment production on the Queensland Plateau; (4) an improvement in CaCO3 preservation near 0.87 Ma, which accompanied the increase of sediment input, indicated by the first appearance of whole pteropods in the deeper Hole 817A and a "spike" in the percent whole pteropods in Hole 818B; (5) a period of strong CaCO3 dissolution during the mid-Brunhes Chron from 0.36 to 0.41 Ma; and (6) a complex CaCO3 preservation pattern between 0.36 Ma and the present characterized by a general increase in CaCO3 preservation through time with good preservation during interglacial stages and poor preservation during glacial stages. The long-term aragonite preservation histories for Holes 818B and 817A appear to be similar in general shape, although different in detail, to CaCO3 preservation records from the deep Indian and central equatorial Pacific oceans as well as from intermediate water sites in the Bahamas and the Maldives. All of these areas have experienced CaCO3 dissolution at about 0.9 Ma and during the mid-Brunhes Chron. However, the late Quaternary (0 to 0.36 Ma) glacial to interglacial preservation pattern in Holes 818B and 817A is out of phase with CaCO3 preservation records for sediments deposited in Pacific deep and bottom waters. The sharp increase in bank production and export from the Queensland Plateau and the coincident improvement of CaCO3 preservation between 0.86 and 0.875 Ma may have been synchronous with the initiation of the Great Barrier Reef and roughly coincides with an increase in carbonate accumulation on the Bahama banks, in the western North Atlantic Ocean, and on Mururoa atoll, in the central South Pacific Ocean. The development of these reef systems during the middle Quaternary may be related to the transition in the frequency and amplitude of global sea level change from 41 k.y. low amplitude cycles prior to 0.9 Ma to 100 k.y. high amplitude cycles after 0.73 Ma. Carbon isotopic analyses show that benthic foraminiferal d13C values (Cibicidoides spp.) have been heavier than planktonic foraminiferal d13C values (G. sacculifer) throughout most of the last 0.54 m.y., which may indicate that 13C-enriched intermediate water (AAIW) occupied the Townsville Trough during much of the late Quaternary. Furthermore, both planktonic and benthic foraminiferal d13C values are often observed to be heaviest during interglacial to glacial transitions, and lightest during glacial to interglacial transitions. We suggest that this pattern is the result of changes in the preformed d13C of Sum CO2 of AAIW and may reflect changes in nutrient utilization by primary producers in Antarctic surface waters, changes in the d13C of upwelled Circumpolar Deep Water, or changes in the extent and/or temperature of equilibration between surface water and atmospheric CO2 within the Antarctic Polar Frontal Zone (the source area for AAIW). Finally, the poor correlation between percent of whole pteropods (aragonite preservation) and d13C of Cibicidoides spp. may be the result of a decoupling of d13C from CO2 due to the numerous and complex variables that combine to produce the preformed d13C of AAIW.

Occurrence and stratigraphy of Coccolithus pelagicus in ODP Hole 120-747A

A record based on counts of the relative abundance of the dominant calcareous nannofossil taxa Coccolithus pelagicus and Reticulofenestra spp. in sediments recovered from Ocean Drilling Program Hole 747A (Kerguelen Plateau, Southern Indian Ocean) is established in this paper. This record (17 m.y. long) virtually spans the entire Miocene. Broad, steplike variations in the abundance of C. pelagicus range between 0% and 96%. Based on these variations, five stratigraphic units characterized by high abundance in C. pelagicus are delineated. We suggest that these variations are caused by water-mass movements (such as the north/south shifting of a front). This pronounced signal is compared with paleoceanographic events revealed by isotopic (d18O and d13C) studies. The five defined units are tentatively correlated to well-known global isotopic events. In particular, Units A and D correlate respectively with the Oligocene/Miocene boundary glaciation and the middle Miocene cooling event. Time-series analysis indicates the presence of the three main periodic components of the eccentricity of the Earth's orbit. A 200-k.y. cycle is also present. The stratigraphic and paleoceanographic significance of this record is discussed.

Mg/Ca ratios in planktonic foraminifera

The primary Mg/Ca ratio of foraminiferal shells is a potentially valuable paleoproxy for sea surface temperature (SST) reconstructions. However, the reliable extraction of this ratio from sedimentary calcite assumes that we can overcome artifacts related to foraminiferal ecology and partial dissolution, as well as contamination by secondary calcite and clay. The standard batch method for Mg/Ca analysis involves cracking, sonicating, and rinsing the tests to remove clay, followed by chemical cleaning, and finally acid-digestion and single-point measurement. This laborious procedure often results in substantial loss of sample (typically 30-60%). We find that even the earliest steps of this procedure can fractionate Mg from Ca, thus biasing the result toward a more variable and often anomalously low Mg/Ca ratio. Moreover, the more rigorous the cleaning, the more calcite is lost, and the more likely it becomes that any residual clay that has not been removed by physical cleaning will increase the ratio. These potentially significant sources of error can be overcome with a flow-through (FT) sequential leaching method that makes time- and labor-intensive pretreatments unnecessary. When combined with time-resolved analysis (FT-TRA) flow-through, performed with a gradually increasing and highly regulated acid strength, produces continuous records of Mg, Sr, Al, and Ca concentrations in the leachate sorted by dissolution susceptibility of the reacting material. Flow-through separates secondary calcite from less susceptible biogenic calcite and clay, and further resolves the biogenic component into primary and more resistant fractions. FT-TRA reliably separates secondary calcite (which is not representative of original life habitats) from the more resistant biogenic calcite (the desired signal) and clay (a contaminant of high Mg/Ca, which also contains Al), and further resolves the biogenic component into primary and more resistant fractions that may reflect habitat or other changes during ontogeny. We find that the most susceptible fraction of biogenic calcite in surface dwelling foraminifera gives the most accurate value for SST and therefore best represents primary calcite. Sequential dissolution curves can be used to correct the primary Mg/Ca ratio for clay, if necessary. However, the temporal separation of calcite from clay in FT-TRA is so complete that this correction is typically <=2%, even in clay-rich sediments. Unlike hands-on batch methods, that are difficult to reproduce exactly, flow-through lends itself to automation, providing precise replication of treatment for every sample. Our automated flow-through system can process 22 samples, two system blanks, and 48 mixed standards in <12 hours of unattended operation. FT-TRA thus represents a faster, cheaper, and better way to determine Mg/Ca ratios in foraminiferal calcite.

Cenozoic oxygen isotopic values for benthic foraminifera from DSDP and ODP low and high latitude marine sediment cores

The climate during the Cenozoic era changed in several steps from ice-free poles and warm conditions to ice-covered poles and cold conditions. Since the 1950s, a body of information on ice volume and temperature changes has been built up predominantly on the basis of measurements of the oxygen isotopic composition of shells of benthic foraminifera collected from marine sediment cores. The statistical methodology of time series analysis has also evolved, allowing more information to be extracted from these records. Here we provide a comprehensive view of Cenozoic climate evolution by means of a coherent and systematic application of time series analytical tools to each record from a compilation spanning the interval from 4 to 61 Myr ago. We quantitatively describe several prominent features of the oxygen isotope record, taking into account the various sources of uncertainty (including measurement, proxy noise, and dating errors). The estimated transition times and amplitudes allow us to assess causal climatological-tectonic influences on the following known features of the Cenozoic oxygen isotopic record: Paleocene-Eocene Thermal Maximum, Eocene-Oligocene Transition, Oligocene-Miocene Boundary, and the Middle Miocene Climate Optimum. We further describe and causally interpret the following features: Paleocene-Eocene warming trend, the two-step, long-term Eocene cooling, and the changes within the most recent interval (Miocene-Pliocene). We review the scope and methods of constructing Cenozoic stacks of benthic oxygen isotope records and present two new latitudinal stacks, which capture besides global ice volume also bottom water temperatures at low (less than 30°) and high latitudes. This review concludes with an identification of future directions for data collection, statistical method development, and climate modeling.

Stable isotope record of benthic foraminifera in late Pleistocene sediments

Variations in the contribution of North Atlantic Deep Water (NADW), relative to North Pacific Deep Water (NPDW), to the Southern Ocean, are assessed by comparing delta13C records from the mid-depth North Atlantic, deep Southern Ocean, and deep equatorial Pacific Ocean. In general, the relative contribution of NADW was greater during interglaciations than glaciations of the past 550,000 years. An increase in the NADW flux to the Southern Ocean since the last glaciation was proposed to have resulted in higher atmospheric CO2 in the Holocene (Broecker and Peng, 1989, doi:10.1029/GB003i003p00215). Glacial-interglacial variations in the proportion of NADW in the Southern Ocean may have also influenced atmospheric CO2 levels over the past 550,000 years. The greatest relative flux of NADW to the Southern Ocean occurred during interglacial stage 11. Faunal data suggest that the North Atlantic polar front and southern Indian Ocean subtropical convergence zone were located farthest poleward during stage 11. Warmth in these locations and a strong southward flux of NADW during stage 11 may be causally linked by the NADW formation process/warm water return route (Gordon, 1986, doi:10.1029/JC091iC04p05037). Time series analysis indicates that delta13C variations in the deep Southern Ocean occur at the same frequencies as the Earth's orbital variations and are coherent and in phase with delta18O. At most, 50% of the glacial-interglacial delta13C amplitude in the Southern Ocean is due changes in the contribution of NADW. The remainder is probably due to mean ocean delta13C changes.