A global compolation of late Holocene planktonic foraminiferal d18O

We review the different sources of uncertainty affecting the oxygen isotopic composition of planktonic foraminifera and present a global planktonic foraminifera oxygen isotope data set that has been assembled within the MARGO project for the Late Holocene time slice. The data set consists of over 2100 data from recent sediment with thorough age control, that have been checked for internal consistency. We further examine how the oxygen isotopic composition of fossil foraminifera is related to hydrological conditions, based on published results on living foraminifera from plankton tows and cultures. Oxygen isotopic values (delta18O) of MARGO recent fossil foraminifera are 0.2-0.8 per mil higher than those of living foraminifera. Our results show that this discrepancy is related to the stratification of the upper water mass and generally increases at low latitudes. Therefore, as stratification of surface waters and seasonality depends on climatic conditions, the relationship between temperature and delta18O established on fossil foraminifera from recent sediment must be used with caution in paleoceanographic studies. Before models predicting seasonal flux, abundance and delta18O composition of a foraminiferal population in the sediment are available, we recommend studying relative changes in isotopic composition of fossil planktonic foraminifera. These changes primarily record variations in temperature and oxygen isotopic composition of sea water, although part of the changes might reflect modifications of planktonic foraminifera seasonality or depth habitat

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.

Major and trace element content of the 17 lithostratigraphic units recovered during DSDP Leg 49 (Table 1)

IPOD Leg 49 recovered basalts from 9 holes at 7 sites along 3 transects across the Mid-Atlantic Ridge: 63°N (Reykjanes), 45°N and 36°N (FAMOUS area). This has provided further information on the nature of mantle heterogeneity in the North Atlantic by enabling studies to be made of the variation of basalt composition with depth and with time near critical areas (Iceland and the Azores) where deep mantle plumes are thought to exist. Over 150 samples have been analysed for up to 40 major and trace elements and the results used to place constraints on the petrogenesis of the erupted basalts and hence on the geochemical nature of their source regions. It is apparent that few of the recovered basalts have the geochemical characteristics of typical "depleted" midocean ridge basalts (MORB). An unusually wide range of basalt compositions may be erupted at a single site: the range of rare earth patterns within the short section cored at Site 413, for instance, encompasses the total variation of REE patterns previously reported from the FAMOUS area. Nevertheless it is possible to account for most of the compositional variation at a single site by partial melting processes (including dynamic melting) and fractional crystallization. Partial melting mechanisms seem to be the dominant processes relating basalt compositions, particularly at 36°N and 45°N, suggesting that long-lived sub-axial magma chambers may not be a consistent feature of the slow-spreading Mid-Atlantic Ridge. Comparisons of basalts erupted at the same ridge segment for periods of the order of 35 m.y. (now lying along the same mantle flow line) do show some significant inter-site differences in Rb/Sr, Ce/Yb, 87Sr/86Sr, etc., which cannot be accounted for by fractionation mechanisms and which must reflect heterogeneities in the mantle source. However when hygromagmatophile (HYG) trace element levels and ratios are considered, it is the constancy or consistency of these HYG ratios which is the more remarkable, implying that the mantle source feeding a particular ridge segment was uniform with respect to these elements for periods of the order of 35 m.y. and probably since the opening of the Atlantic. Yet these HYG element ratios at 63°N are very different from those at 45°N and 36°N and significantly different from the values at 22°N and in "MORB". The observed variations are difficult to reconcile with current concepts of mantle plumes and binary mixing models. The mantle is certainly heterogeneous, but there is not simply an "enriched" and a "depleted" source, but rather a range of sources heterogeneous on different scales for different elements - to an extent and volume depending on previous depletion/enrichment events. HYG element ratios offer the best method of defining compositionally different mantle segments since they are little modified by the fractionation processes associated with basalt generation.

Magnetic stratigraphy of ODP Leg 133 sites

This study is a synthesis of paleomagnetic and mineral magnetic results for Sites 819 through 823 of Ocean Drilling Program (ODP) Leg 133, which lie on a transect from the outer edge of the Great Barrier Reef (GBR) down the continental slope to the bottom of the Queensland Trough. Because of viscous remagnetization and pervasive overprinting, few reversal boundaries can be identified in these extremely high-resolution Quaternary sequences. Some of the magnetic instability, and the differences in the quality of the paleomagnetic signal among sites, can be explained in terms of the dissolution of primary iron oxides in the high near-surface geochemical gradients. Well-defined changes in magnetic properties, notably susceptibility, reflect responses to glacio-eustatic sea-level fluctuations and changes in slope sedimentation processes resulting from formation of the GBR. Susceptibility can be used to correlate between adjacent holes at a given site to an accuracy of about 20 cm. Among-site correlation of susceptibility is also possible for certain parts of the sequences and permits (tentative) extension of the reversal chronology. The reversal boundaries that can be identified are generally compatible with the calcareous nannofossil biostratigraphy and demonstrate a high level of biostratigraphic consistency among sites. A revised chronology based on an optimum match with the susceptibility stratigraphy is presented. Throughout most of the sequences there is a strong inverse correlation both between magnetic susceptibility and calcium carbonate content, and between susceptibility and d18O. In the upper, post-GBR, sections a more complicated type of magnetic response occurs during glacial maxima and subsequent transgressions, resulting in a positive correlation between susceptibility and d18O. Prior to and during formation of the outer-reef barrier, the sediments have relatively uniform magnetic properties showing multidomain behavior and displaying cyclic variations in susceptibility related to sea-level change. The susceptibility oscillations are controlled more by carbonate dilution than by variation in terrigenous influx. Establishment of the outer reef between 1.01 and 0.76 Ma restricted the supply of sediment to the slope, causing a four-fold reduction in sedimentation rates and a transition from prograding to aggrading seismic geometries (see other chapters in this volume). The Brunhes/Matuyama boundary and the end of the transition period mark a change to lower and more subdued susceptibility oscillations with higher carbonate contents. The major change in magnetic properties comes at about 0.4 Ma in the aggrading sequence, which contains prominent sharp susceptibility peaks associated with glacial cycles, with distinctive single-domain magnetite and mixed single-domain/superparamagnetic characteristics. Bacterial magnetite has been found in the sediments, particularly where there are high susceptibility peaks, but its importance has not yet been assessed. A possible explanation for the characteristic pattern of magnetic properties in the post-GBR glacial cycles can be found in terms of fluvio-deltaic processes and inter-reefal lagoonal reservoirs that develop when the shelf becomes exposed at low sea-level.

Geochemistry of basalts from ODP Leg 49, North Atlantic Ocean

The geochemistry of basalts recovered from seven sites in the North Atlantic is described with particular reference to minor elements. Three sites (407, 408, and 409) along the same mantle flow line, transverse to the Reykjanes Ridge at about 63°N, provide information on the composition of basalts erupted over a 34-m.y. interval between 2.3 and 36 m.y. ago. At Site 410, at 45°N, penetration into 10 m.y.-old crust west of the ridge axis permits comparisons with young basalts dredged from the median valley at 45°N. Three sites in the FAMOUS area at about 36°N provided material from very young (1 m.y.) basaltic crust (Site 411), and material to test the geochemical coherence of basalts of different ages (1.5 and 3.5 m.y.) on either side of a fracture zone (Sites 412 and 413). These sites complement earlier data from dredged and drilled sites (Leg 37) in the FAMOUS area. At Site 407, four geochemically distinct basalt units occur, with different normative and rare-earth element (REE) characteristics, and there is a clear correlation with magnetic stratigraphy. Yet there is a remarkable consistency in incompatible element ratios between these units, indicating derivation from an essentially similar mantle source. The basalts from the younger sites, 408 and 409, show a similar range of normative and REE variation, but incompatible element ratios are identical to those at Site 407, indicating that basalts at all three sites were produced from a mantle source which was geochemically relatively uniform. Rare-earth differences between the basalts can be interpreted in terms of variations in the degree and depth of partial melting causing HREE (+Y) retention in the source, although there may be some inter-site differences with respect to REE. A similar picture is presented at 45°N. Apparently a range of tholeiitic, transitional, and alkalic basalts were being erupted 10 m.y. ago, which have almost identical geochemical characteristics to those recently erupted in the median valley at 45°N. Incompatible element ratios are markedly different from those recorded at the Reykjanes Ridge. Basalts recovered from the FAMOUS sites are geochemically similar to previous samples recovered from the FAMOUS area, and their incompatible element ratios are similar, but not identical, to those at 45°N. However, total trace element levels are consistently lower than in 45°N basalts, which might imply smaller degrees of partial melting and/or greater depths of magma generation at 45°N, or higher trace element levels in the mantle source at 45°N. Few of the basalts recovered on Leg 49 have the geochemical characteristics of typical "MORB" (e.g., Nazca Plate, Leg 34). The data strongly support models invoking geochemical inhomogeneity in the source regions of basalts produced at the Mid-Atlantic Ridge. However, the data also introduce an additional time factor into such models and demonstrate the uniformity of the mantle source at a particular ridge sector (over periods in excess of 30 m.y.), while emphasizing the marked differences along the ridge. Mixing models invoking "depleted" and "enriched" mantle sources would seem to be inadequate to account for the observed variations.

Towards convection-resolving, global atmospheric simulations with the Model for Prediction Across Scales (MPAS) v3.1: an extreme scaling experiment, links to sourcecode of model and grides

The Model for Prediction Across Scales (MPAS) is a novel set of Earth system simulation components and consists of an atmospheric model, an ocean model and a land-ice model. Its distinct features are the use of unstructured Voronoi meshes and C-grid discretisation to address shortcomings of global models on regular grids and the use of limited area models nested in a forcing data set, with respect to parallel scalability, numerical accuracy and physical consistency. This concept allows one to include the feedback of regional land use information on weather and climate at local and global scales in a consistent way, which is impossible to achieve with traditional limited area modelling approaches. Here, we present an in-depth evaluation of MPAS with regards to technical aspects of performing model runs and scalability for three medium-size meshes on four different high-performance computing (HPC) sites with different architectures and compilers. We uncover model limitations and identify new aspects for the model optimisation that are introduced by the use of unstructured Voronoi meshes. We further demonstrate the model performance of MPAS in terms of its capability to reproduce the dynamics of the West African monsoon (WAM) and its associated precipitation in a pilot study. Constrained by available computational resources, we compare 11-month runs for two meshes with observations and a reference simulation from the Weather Research and Forecasting (WRF) model. We show that MPAS can reproduce the atmospheric dynamics on global and local scales in this experiment, but identify a precipitation excess for the West African region. Finally, we conduct extreme scaling tests on a global 3?km mesh with more than 65 million horizontal grid cells on up to half a million cores. We discuss necessary modifications of the model code to improve its parallel performance in general and specific to the HPC environment. We confirm good scaling (70?% parallel efficiency or better) of the MPAS model and provide numbers on the computational requirements for experiments with the 3?km mesh. In doing so, we show that global, convection-resolving atmospheric simulations with MPAS are within reach of current and next generations of high-end computing facilities.