An individual based model spatial distribution, production and biomass for Calanus finmarchicus in the Norwegian Sea

The copepod Calanus finmarchicus is the dominant species of the meso-zooplankton in the Norwegian Sea, and constitutes an important link between the phytoplankton and the higher trophic levels in the Norwegian Sea food chain. An individualbased model for C. finmarchicus, based on super-individuals and evolving traits for behaviour, stages, etc., is two-way coupled to the NORWegian ECOlogical Model system (NORWECOM). One year of modelled C. finmarchicus spatial distribution, production and biomass are found to represent observations reasonably well. High C. finmarchicus abundance is found along the Norwegian shelf-break in the early summer, while the overwintering population is found along the slope and in the deeper Norwegian Sea basins. The timing of the spring bloom is generally later than in the observations. Annual Norwegian Sea production is found to be 29 million tonnes of carbon and a production to biomass (P/B) ratio of 4.3 emerges. Sensitivity tests show that the modelling system is robust to initial values of behavioural traits and with regards to the number of super-individuals simulated given that this is above about 50,000 individuals. Experiments with the model system indicate that it provides a valuable tool for studies of ecosystem responses to causative forces such as prey density or overwintering population size. For example, introducing C. finmarchicus food limitations reduces the stock dramatically, but on the other hand, a reduced stock may rebuild in one year under normal conditions. The NetCDF file contains model grid coordinates and bottom topography.

A global seasonal surface ocean climatology of phytoplankton types based on CHEMTAX analysis of HPLC pigments

Much advancement has been made in recent years in field data assimilation, remote sensing and ecosystem modeling, yet our global view of phytoplankton biogeography beyond chlorophyll biomass is still a cursory taxonomic picture with vast areas of the open ocean requiring field validations. High performance liquid chromatography (HPLC) pigment data combined with inverse methods offer an advantage over many other phytoplankton quantification measures by way of providing an immediate perspective of the whole phytoplankton community in a sample as a function of chlorophyll biomass. Historically, such chemotaxonomic analysis has been conducted mainly at local spatial and temporal scales in the ocean. Here, we apply a widely tested inverse approach, CHEMTAX, to a global climatology of pigment observations from HPLC. This study marks the first systematic and objective global application of CHEMTAX, yielding a seasonal climatology comprised of ~1500 1°x1° global grid points of the major phytoplankton pigment types in the ocean characterizing cyanobacteria, haptophytes, chlorophytes, cryptophytes, dinoflagellates, and diatoms, with results validated against prior regional studies where possible. Key findings from this new global view of specific phytoplankton abundances from pigments are a) the large global proportion of marine haptophytes (comprising 32 ± 5% of total chlorophyll), whose biogeochemical functional roles are relatively unknown, and b) the contrasting spatial scales of complexity in global community structure that can be explained in part by regional oceanographic conditions. These publicly accessible results will guide future parameterizations of marine ecosystem models exploring the link between phytoplankton community structure and marine biogeochemical cycles.

Analysis of an ancient global warming event, the Paleocene-Eocene thermal maximum, under examination from ODP Holes 113-689B and 113-690B

Spatiotemporal patterns of carbonate dissolution provide a critical constraint on carbon input during an ancient (~55.5 Ma) global warming event known as the Paleocene-Eocene thermal maximum (PETM), yet the magnitude of lysocline shoaling in the Southern Ocean is poorly constrained due to limited spatial coverage in the circum-Antarctic region. This shortcoming is partially addressed by comparing patterns of carbonate sedimentation at the Site 690 PETM reference section to those herein reconstructed for nearby Site 689. Biochemostratigraphic correlation of the two records reveals that the first ~36 ka of the carbon isotope excursion (CIE) signaling PETM conditions is captured by the Site 689 section, while the remainder of the CIE interval and nearly all of the CIE recovery are missing due to a coring gap. A relatively expanded stratigraphy and higher carbonate content at mid-bathyal Site 689 indicate that dissolution was less severe than at Site 690. Thus, the bathymetric transect delimited by these two PETM records indicates that the lysocline shoaled above Site 689 (~1,100 m) while the calcite compensation depth remained below Site 690 (~1,900 m) in the Weddell Sea region. The ensuing recovery of carbonate sedimentation conforms to a bathymetric trend best explained by gradual lysocline deepening as negative feedback mechanisms neutralized ocean acidification. Further, biochemostratigraphic evidence indicates the tail end of the CIE recovery interval at both sites has been truncated by a hiatus most likely related to vigorous production and advection of intermediate waters.

Acoustic multi-frequency indicator of four major groups on a spatial grid in the Bay of Biscay

The first data set contains the mean and cofficient of variation (standard deviation divided by mean) of a multi-frequency indicator I derived from ER60 acoustic information collected at five frequencies (18, 38, 70, 120, and 200 kHz) in the Bay of Biscay in May of the years 2006, 2008, 2009 and 2010 (Pelgas surveys). The multi-frequency indicator was first calculated per voxel (20 m long × 5 m deep sampling unit) and then averaged on a spatial grid (approx. 20 nm × 20 nm) for five 5-m depth layers in the surface waters (10-15m, 15-20m, 20-25m, 25-30m below sea surface); there are missing values in particular in the shallowest layer. The second data set provides for each grid cell and depth layer the proportion of voxels for which the multi-frequency indicator I was indicative of a certain group of organisms. For this the following interpretation was used: I < 0.39 swim bladder fish or large gas bubbles, I = 0.39-0.58 small resonant bubbles present in gas bearing organisms such as larval fish and phytoplankton, I = 0.7-0.8 fluidlike zooplankton such as copepods and euphausiids, and I > 0.8 mackerel. These proportions can be interpreted as a relative abundance index for each of the four organism groups.

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.

Biomass and stable isotope composition of plankton in five size fractions along a zonal transect (24°N) from the Canary Islands to Florida in January-March 2011

The spatial variability of biomass and stable isotopes in plankton size fractions in the upper 200 m was studied in a high spatial resolution transect along 24°N from Canary Islands to Florida (January - March 2011) during Leg 8 of the Malaspina-2010 expedition (http://www.expedicionmalaspina.es) to determine nitrogen and carbon sources. Plankton samples were collected by vertical tows of a microplankton net (40 mm mesh size) and a mesoplankton net (200 mm mesh size) through the upper 200 m of the water column. Sampling was between 10:00 and 16:00 h GMT. Plankton was separated into five size fractions (40 - 200, 200 - 500, 500 - 1000, 1000 - 2000 and > 2000 mm) by gentle filtration of the samples by a graded series of nylon sieves (2000, 1000, 500, 200 and 40 mm). Large gelatinous organisms were removed before filtration. Aliquots for each size fraction were collected on pre-weighed glass-fibre filters, dried (60°C, 48 h) and stored in a desiccator before determination of biomass (dry weight), carbon and nitrogen content and natural abundance of stable carbon and nitrogen isotopes ashore. Vertical advection of waters predominated in lateral zones while the central Atlantic (30-70°W) was characterized by a strong stratification and oligotrophic surface waters. Plankton biomass was low in the central zone and high in both eastern and western sides, with most of the variability due to either large (>2000 µm) and small plankton (<500 µm). Carbon isotopes reflected mainly the advection the deep water in lateral zones. Stable nitrogen isotopes showed a nearly symmetrical spatial distribution in all fractions, with the lowest values (delta15N <1per mill) in the central zone, and were inversely correlated to carbon stable isotopes (delta13C) and to the abundance of the nitrogen-fixer Trichodesmium. Diazotrophy was estimated to account for >50% of organic nitrogen in the central zone, and even >30% in eastern and western zones. The impact of diazotrophy increased with the size of the organisms, supporting the wide participation of all trophic levels in the processing of recently fixed nitrogen. These results indicate that atmospheric sources of carbon and nitrogen prevail over deep water sources in the subtropical North Atlantic and that the zone influenced by diazotrophy is much larger than reported in previous studies.

U-Pb, trace element, and Hf isotope analysis for Macquarie Island detrital zircon

Oceanic zircon trace element and Hf-isotope geochemistry offers a means to assess the magmatic evolution of a dying spreading ridge and provides an independent evaluation of the reliability of oceanic zircon as an indicator of mantle melting conditions. The Macquarie Island ophiolite in the Southern Ocean provides a unique testing ground for this approach due to its formation within a mid-ocean ridge that gradually changed into a transform plate boundary. Detrital zircon recovered from the island records this change through a progressive enrichment in incompatible trace elements. Oligocene age (33-27 Ma) paleo-detrital zircon in ophiolitic sandstones and breccias interbedded with pillow basalt have trace element compositions akin to a MORB crustal source, whereas Late Miocene age (8.5 Ma) modern-detrital zircon collected from gabbroic colluvium on the island have highly enriched compositions unlike typical oceanic zircon. This compositional disparity between age populations is not complimented by analytically equivalent eHf data that primarily ranges from 14 to 13 for sandstone and modern-detrital populations. A wider compositional range for the sandstone population reflects a multiple pluton source provenance and is augmented by a single cobble clast with eHf equivalent to the maximum observed composition in the sandstone (~17). Similar sandstone and colluvium Hf-isotope signatures indicate inheritance from a similar mantle reservoir that was enriched from the depleted MORB mantle average. The continuity in Hf-isotope signature relative to trace element enrichment in Macquarie Island zircon populations, suggests the latter formed by reduced partial melting linked to spreading-segment shortening and transform lengthening along the dying spreading ridge.

Growth rates of late Miocene corals from Crete (Greece)

Modern scleractinian corals are classical components of marine shallow warm water ecosystems. Their occurrence and diversity patterns in the geological record have been widely used to infer past climates and environmental conditions. Coral skeletal composition data reflecting the nature of the coral environment are often affected by diagenetic alteration. Ghost structures of annual growth rhythms are, however, often well preserved in the transformed skeleton. We show that these relicts represent a valuable source of information on growth conditions of fossil corals. Annual growth bands were measured in massive hemispherical Porites of late Miocene age from the island of Crete (Greece) that were found in patch reefs and level bottom associations of attached mixed clastic environments as well as isolated carbonate environments. The Miocene corals grew slowly, about 2-4 mm/yr, compatible with present-day Porites from high-latitude reefs. Slow annual growth of the Miocene corals is in good agreement with the position of Crete at the margin of the Miocene reef belt. Within a given time slice, extension rates were lowest in level bottom environments and highest in attached inshore reef systems. Because sea surface temperatures (SSTs) can be expected to be uniform within a time slice, spatial variations in extension rates must reflect local variations in light levels (low in the level bottom communities) and nutrients (high in the attached reef systems). During the late Miocene (Tortonian-early Messinian), maximum linear extension rates remained remarkably constant within seven chronostratigraphic units, and if the relationship of SSTs and annual growth rates observed for modern massive Indo-Pacific Porites spp. applies to the Neogene, minimum (winter) SSTs were 20°-21°C. Although our paleoclimatic record has a low resolution, it fits the trends revealed by global data sets. In the near future we expect this new and easy to use Porites thermometer to add important new information to our understanding of Neogene climate.

Global LAI uncertainty and relative uncertainty datasets derived from a Triple Collocation Error Model (TCEM)

Uncertainty information for global leaf area index (LAI) products is important for global modeling studies but usually difficult to systematically obtain at a global scale. Here, we present a new method that cross-validates existing global LAI products and produces consistent uncertainty information. The method is based on a triple collocation error model (TCEM) that assumes errors among LAI products are not correlated. Global monthly absolute and relative uncertainties, in 0.05° spatial resolutions, were generated for MODIS, CYCLOPES, and GLOBCARBON LAI products, with reasonable agreement in terms of spatial patterns and biome types. CYCLOPES shows the lowest absolute and relative uncertainties, followed by GLOBCARBON and MODIS. Grasses, crops, shrubs, and savannas usually have lower uncertainties than forests in association with the relatively larger forest LAI. With their densely vegetated canopies, tropical regions exhibit the highest absolute uncertainties but the lowest relative uncertainties, the latter of which tend to increase with higher latitudes. The estimated uncertainties of CYCLOPES generally meet the quality requirements (± 0.5) proposed by the Global Climate Observing System (GCOS), whereas for MODIS and GLOBCARBON only non-forest biome types have met the requirement. Nevertheless, none of the products seems to be within a relative uncertainty requirements of 20%. Further independent validation and comparative studies are expected to provide a fair assessment of uncertainties derived from TCEM. Overall, the proposed TCEM is straightforward and could be automated for the systematic processing of real time remote sensing observations to provide theoretical uncertainty information for a wider range of land products.

Tie points and goechemical analysis from different Holes of ODP Site 165-1002

The fine-grained sediments of the Cariaco Basin, Venezuela, of the last 130 ky, whose deposition history is well characterized, were analyzed geochemically in order to test the validity of sediment bulk geochemistry as an indicator of detrital provenance. Several binary and ternary diagrams as well as the chemical index of alteration (CIA) were tested for their capacity to discriminate the poorly contrasted detrital sources to the Cariaco Basin, and to describe the temporal evolution of the contributions of these different sources. Most of the diagrams tested did not allow a good discrimination of sources or, when sources were well discriminated, did not allow an interpretation of the temporal variations consistent with the known history. A relatively good discrimination of sources and a consistent interpretation of temporal variations were however obtained using Hf vs. Th and La/Yb vs. Gd/Yb binary diagrams, as well as Ti-Zr-Th, Ti-Zr-La, and Lu-Hf-Th ternary diagrams. Compared to the previous studies of the detrital content of the Cariaco Basin sediments, the geochemical approach permitted the recognition of a sediment contribution eroded from the Unare platform and Gulf of Cariaco during rapid sea level oscillations, and the contribution of Saharan eolian particles during the Younger Dryas-Preboreal and MIS6-5 transition. The choice of plotted elements was determined after considering carrier minerals, so that different elements may be informative in different sedimentary contexts. Overall, mineral sorting during transport appears as a major limit to quantitative estimation of the different contributions. In particular mineral sorting leads to the selective enrichment of elements associated with clays (Al, Rb, Th and LREE) in sediments deposited in the basin. Unless the geochemical effect of mineral sorting can be measured, it appears that quantitative provenance analysis should be performed on fractions of similar grain size instead of bulk sediment.