Integrated Data Sets of the EU FP5 Research Project ORFOIS: Origin and fate of biogenic particle fluxes in the ocean and their interactions with atmospheric CO2 concentrations as well as the marine sediment

For a reliable simulation of the time and space dependent CO2 redistribution between ocean and atmosphere an appropriate time dependent simulation of particle dynamics processes is essential but has not been carried out so far. The major difficulties were the lack of suitable modules for particle dynamics and early diagenesis (in order to close the carbon and nutrient budget) in ocean general circulation models, and the lack of an understanding of biogeochemical processes, such as the partial dissolution of calcareous particles in oversaturated water. The main target of ORFOIS was to fill in this gap in our knowledge and prediction capability infrastructure. This goal has been achieved step by step. At first comprehensive data bases (already existing data) of observations of relevance for the three major types of biogenic particles, organic carbon (POC), calcium carbonate (CaCO3), and biogenic silica (BSi or opal), as well as for refractory particles of terrestrial origin were collated and made publicly available.

Modern and fossil corals from Palmyra and Christmas islands

The relationship between decadal to centennial changes in ocean circulation and climate is difficult to discern using the sparse and discontinuous instrumental record of climate and, as such, represents a large uncertainty in coupled ocean-atmosphere general circulation models. We present new modern and fossil coral radiocarbon (D14C) records from Palmyra (6°N, 162°W) and Christmas (2°N, 157°W) islands to constrain central tropical Pacific ocean circulation changes during the last millennium. Seasonally to annually resolved coral D14C measurements from the 10th, 12th-17th, and 20th centuries do not contain significant interannual to decadal-scale variations, despite large changes in coral d18O on these timescales. A centennial-scale increase in coral radiocarbon from the Medieval Climate Anomaly (~900-1200 AD) to the Little Ice Age (~1500-1800) can be largely explained by changes in the atmospheric D14C, as determined with a box model of Palmyra mixed layer D14C. However, large 12th century depletions in Palmyra coral D14C may reflect as much as a 100% increase in upwelling rates and/or a significant decrease in the D14C of higher-latitude source waters reaching the equatorial Pacific during this time. SEM photos reveal evidence for minor dissolution and addition of secondary aragonite in the fossil corals, but our results suggest that coral D14C is only compromised after moderate to severe diagenesis for these relatively young fossil corals.

Reconstructed glaciation limits in Eurasia (digitized from map Fig.1)

Based on field investigations in northern Russia and interpretation of offshore seismic data, we have made a preliminary reconstruction of the maximum ice-sheet extent in the Barents and Kara Sea region during the Early/Middle Weichselian and the Late Weichselian. Our investigations indicate that the Barents and Kara ice sheets attained their maximum Weichselian positions in northern Russia prior to 50 000 yr BP, whereas the northeastern flank of the Scandinavian Ice Sheet advanced to a maximum position shortly after 17000 calendar years ago. During the Late Weichselian (25 000-10000 yr BP), much of the Russian Arctic remained ice-free. According to our reconstruction, the extent of the ice sheets in the Barents and Kara Sea region during the Late Weichselian glacial maximum was less than half that of the maximum model which, up to now, has been widely used as a boundary condition for testing and refining General Circulation Models (GCMs). Preliminary numerical-modelling experiments predict Late Weichselian ice sheets which are larger than the ice extent implied for the Kara Sea region from dated geological evidence, suggesting very low precipitation.

Plankton abundances in water samples of cruise EBENE

The Etude du Broutage en Zone Equatoriale (EBENE) transect (8°S-8°N) explored the equatorial high-nutrient, low-chlorophyll (HNLC) zone and adjacent oligotrophic areas during a La Niña period (October-November 1996). During this time the passage of a tropical instability wave also influenced the region north of the equator. We present a brief summary of EBENE findings, with an emphasis on phytoplankton utilization by the assemblage of protistan and animal consumers. Despite significant variability over the diel cycle, phytoplankton biomass at the equator was relatively constant on a 24-hour timescale, denoting a dynamic balance between growth and losses. The magnitude of the daily cycle in phytoplankton biomass was well constrained by in situ observations of the diel variability in pigments and suspended particulates, by 14C uptake rates from in situ incubations, and from experimental determinations of specific growth and grazing rates. The general equilibrium of production and grazing processes is illustrated by applying biomass-specific grazing rates from the equatorial station to measured planktonic biomass along the EBENE transect and comparing them to measured 14C uptake. Most of the grazing turnover is supported by the production of Prochloroccus (31%) and picoeukaryotic algae (34%). Among the consumers, microzooplankton (<200 µm) account for 59-98% of the grazing losses. The coherence of the results obtained by independent methods suggests that the essential features of the system have been adequately represented by rate and standing stock assessments from the EBENE study.

Radiocarbon ages and stable isotope measurements on sediment cores of the subarctic Pacific and its marginal seas

Under modern conditions only North Pacific Intermediate Water is formed in the northwest Pacific Ocean. This situation might have changed in the past. Recent studies with general circulation models indicate a switch to deep-water formation in the northwest Pacific during Heinrich Stadial 1 (17.5-15.0 ka) of the last glacial termination. Reconstructions of past ventilation changes based on paleoceanographic proxy records are still insufficient to test whether a deglacial mode of deep-water formation in the North Pacific Ocean existed. Here we present deglacial ventilation records based on radiocarbon-derived ventilation ages in combination with epibenthic stable carbon isotopes from the northwest Pacific including the Okhotsk Sea and Bering Sea, the two potential source regions for past North Pacific ventilation changes. Evidence for most rigorous ventilation of the intermediate-depth North Pacific occurred during Heinrich Stadial 1 and the Younger Dryas, simultaneous to significant reductions in Atlantic Meridional Overturning Circulation. Concurrent changes in d13C and ventilation ages point to the Okhotsk Sea as driver of millennial-scale changes in North Pacific Intermediate Water ventilation during the last deglaciation. Our records additionally indicate that changes in the d13C intermediate-water (700-1750 m water depth) signature and radiocarbon-derived ventilation ages are in antiphase to those of the deep North Pacific Ocean (>2100 m water depth) during the last glacial termination. Thus, intermediate- and deep-water masses of the northwest Pacific have a differing ventilation history during the last deglaciation.

Planktonic foraminifera distribution during the mid-Piacenzian warm period in the North Atlantic

The mid-Piacenzian (MP) warm period (3.264-3.025 Ma) has been identified as the most recent time in geologic history during which mean global surface temperatures were considerably warmer than today for a sustained period. This interval has therefore been proposed as a potential (albeit imperfect) analog for future climate change and as such, has received much scientific attention over the past two decades. Central to this research effort is the Pliocene Research, Interpretation, and Synoptic Mapping (PRISM) project, an iterative paleoenvironmental reconstruction of the MP focused on increasing our understanding of warm-period climate forcings, dynamics, and feedbacks by providing three-dimensional data sets for general circulation models. A mainstay of the PRISM project has been the development of a global sea surface temperature (SST) data set based primarily upon quantitative analyses of planktic foraminifer assemblages, supplemented with geochemical SST estimates wherever possible. In order to improve spatial coverage of the PRISM faunal data set in the low and mid-latitude North Atlantic, this study provides a description of the MP planktic foraminifer assemblage from five Ocean Drilling Program sites (951, 958, 1006, 1062, and 1063) in the subtropical gyre, a region critical to Atlantic Ocean circulation and tropical heat advection. Assemblages from each core provide evidence for a temperature- and circulation-driven 5-10° northward displacement of MP faunal provinces, as well as regional shifts in planktic foraminifer populations linked to species ecology and interactions. General biogeographic trends also indicate that, relative to modern conditions, gyre circulation was stronger (particularly the Gulf Stream, North Atlantic Current, and North Equatorial Current) and meridionally broader. A comparison of mid-Piacenzian and modern North Atlantic planktic foraminifer assemblages suggests that low latitude western boundary currents were less than 1 °C warmer while eastern boundary currents were ~1-2 °C warmer, supporting the hypothesis of enhanced northward heat advection along western boundary currents and warming of high latitude Northeast Atlantic source regions for the Canary Current. These findings are consistent with a model of reduced meridional SST gradients, with little-to-no low latitude warming, and more vigorous ocean circulation. Results therefore support the theory that enhanced meridional overturn circulation and associated northward heat advection made an important contribution, in conjunction with elevated atmospheric CO2 concentrations, to the 2-3 °C global surface temperature increase (relative to today) and strong polar amplification of SST warmth during the MP warm period.

Radiogenic isotope trace composition of the Labrador Sea

Geological reconstructions and general circulation models suggest that the onset of both Northern Hemisphere glaciation, 2.7 Myr ago, and convection of Labrador Sea Water (LSW) were caused by the closure of the Panama Gateway ~4.5 Myr ago. Time series data that have been obtained from studies of ferromanganese crusts from the northwestern Atlantic suggest that radiogenic isotopes of intermediate ocean residence time (Pb and Nd) can serve as suitable tracers to reconstruct these events. However, it has been unclear until now as to whether the changes that have been observed in isotope composition at this time are the result of increased thermohaline circulation or due to the effects of increased glacial weathering. In this paper we adopt a box model approach to demonstrate that the shifts in radiogenic isotope compositions are unlikely to be due to changes in convection in LSW but can be explained in terms of increases of erosion levels due to the glaciation of Greenland and Canada. Furthermore, we provide experimental evidence for the incongruent release of a labile fraction of strongly radiogenic Pb and nonradiogenic Nd from continental detritus eroding into the Labrador Sea. This can be attributed to the glacial weathering of old continents and accounts for the paradox that one of the areas of the world most deficient in radiogenic Pb should provide such a rich supply of radiogenic Pb to the oceans. An important general conclusion is that the compositions of radiogenic isotopes in seawater are not always a reflection of their continental sources. Perhaps more importantly, the transition from chemical weathering to mechanical erosion is likely to result in significant variations in radiogenic tracers in seawater.

(Table 1) 87Sr/86Sr ratios and ages of DSDP Site 144A samples

Glassy Turonian foraminifera preserved in clay-rich sediments from the western tropical Atlantic yield the warmest equivalent d18O sea-surface temperatures (SSTs) yet reported for the entire Cretaceous-Cenozoic. We estimate Turonian SSTs that were at least as warm as (conservative mean ~30 °C) to significantly warmer (warm mean ~33 °C) than those in the region today. However, if independent evidence for high middle Cretaceous pCO2 is reliable and resulted in greater isotopic fractionation between seawater and calcite because of lower sea-surface pH, our conservative and warm SST estimates would be even higher (32 and 36°C, respectively). Our new tropical SSTs help reconcile geologic data with the predictions of general circulation models that incorporate high Cretaceous pCO2 and lend support to the hypothesis of a Cretaceous greenhouse. Our data also strengthen the case for a Turonian age for the Cretaceous thermal maximum and highlight a 20-40 m.y. mismatch between peak Cretaceous-Cenozoic global warmth and peak inferred tectonic CO2 production. We infer that this mismatch is either an artifact of a hidden Turonian pulse in global ocean-crust cycling or real evidence of the influence of some other factor on atmospheric CO2 and/or SSTs. A hidden pulse in crust cycling would explain the timing of peak Cretaceous-Cenozoic sea level (also Turonian), but other factors are needed to explain high-frequency (~10-100 k.y.) instability in middle Cretaceous SSTs reported elsewhere.

A 300 kyr record of upwelling off Oman during the Late Quaternary: evidence of the Asian southwest monsoon

In the northwest Arabian Sea upwelling occurs each summer, driven by the strong SW monsoon winds. Upwelling results in high biological productivity and a distinctive assemblage of plankton species in the surface waters off Oman that are preserved in the sediments along the Oman continental margin, creating a geologic record of monsoon-driven upwelling. Sediments recovered from the Oman continental margin during Ocean Drilling Program leg 117 provide an opportunity to examine how upwelling has varied during the late Quaternary, spanning a longer interval than piston cores recovered prior to the ODP cruise. Variations in foraminifer shell accumulation and in the relative abundance of Globigerina bulloides indicate dominant cycles of variation at 1/100 kyr, the dominant frequency of glacial-interglacial variations, and at 1/23 kyr, the frequency of precessionally driven cycles in seasonal insolation. The strongest monsoon winds (indicated by increased upwelling) occurred during interglacial times when perihelion was aligned with the summer solstice, an orbital change that increased the insolation received during summer in the northern hemisphere. During glacial times upwelling was reduced, and although the precessional cycles were still present their amplitude was smaller. At both frequencies the upwelling cycles are in phase with minimum ice volume, evidence that glacial-interglacial climate changes also include changes to the climate system that influence the low-latitude monsoon. We attribute the decrease in the monsoon winds observed during glacial times to changes in bare land albedo over Asia and/or to changes in the areal extent and seasonal cycle in Asian snow cover that decrease the summer land-sea temperature contrast. Other mechanisms may also be involved. These new upwelling time series differ substantially from previous results, however the previous work relied on cores located farther offshore where upwelling is less intense and other physical mechanisms become important. Our results support the observations derived from atmospheric general circulation models of the atmosphere that indicate that both glacial boundary conditions, and the strength of summer insolation are important variables contributing to cycles in the monsoon winds during the late Quaternary.

On the Southern Ocean CO2 uptake and the role of the biological carbon pump in the 21st century, links to supplementary material

We use a suite of eight ocean biogeochemical/ecological general circulation models from the MAREMIP and CMIP5 archives to explore the relative roles of changes in winds (positive trend of Southern Annular Mode, SAM) and in warming- and freshening-driven trends of upper ocean stratification in altering export production and CO2 uptake in the Southern Ocean at the end of the 21st century. The investigated models simulate a broad range of responses to climate change, with no agreement ona dominance of either the SAM or the warming signal south of 44° S. In the southernmost zone, i.e., south of 58° S, they concur on an increase of biological export production, while between 44 and 58° S the models lack consensus on the sign of change in export. Yet, in both regions, the models show an enhanced CO2 uptake during spring and summer. This is due to a larger CO 2 (aq) drawdown by the same amount of summer export production at a higher Revelle factor at the end of the 21st century. This strongly increases the importance of the biological carbon pump in the entire Southern Ocean. In the temperate zone, between 30 and 44° S all models show a predominance of the warming signal and a nutrient-driven reduction of export production. As a consequence, the share of the regions south of 44° S to the total uptake of the Southern Ocean south of 30° S is projected to increase at the end of the 21st century from 47 to 66% with a commensurable decrease to the north. Despite this major reorganization of the meridional distribution of the major regions of uptake, the total uptake increases largely in line with the rising atmospheric CO2. Simulations with the MITgcm-REcoM2 model show that this is mostly driven by the strong increase of atmospheric CO2, with the climate-driven changes of natural CO2 exchange offsetting that trend only to a limited degree (~10%) and with negligible impact of climate effects on anthropogenic CO2 uptake when integrated over a full annual cycle south of 30° S.

Last Interglacial synthesis of high-latitude temperature: temperature anomalies and associated errors for 4 time slices

The Last Interglacial (LIG, 129-116 thousand of years BP, ka) represents a test bed for climate model feedbacks in warmer-than-present high latitude regions. However, mainly because aligning different palaeoclimatic archives and from different parts of the world is not trivial, a spatio-temporal picture of LIG temperature changes is difficult to obtain. Here, we have selected 47 polar ice core and sub-polar marine sediment records and developed a strategy to align them onto the recent AICC2012 ice core chronology. We provide the first compilation of high-latitude temperature changes across the LIG associated with a coherent temporal framework built between ice core and marine sediment records. Our new data synthesis highlights non-synchronous maximum temperature changes between the two hemispheres with the Southern Ocean and Antarctica records showing an early warming compared to North Atlantic records. We also observe warmer than present-day conditions that occur for a longer time period in southern high latitudes than in northern high latitudes. Finally, the amplitude of temperature changes at high northern latitudes is larger compared to high southern latitude temperature changes recorded at the onset and the demise of the LIG. We have also compiled four data-based time slices with temperature anomalies (compared to present-day conditions) at 115 ka, 120 ka, 125 ka and 130 ka and quantitatively estimated temperature uncertainties that include relative dating errors. This provides an improved benchmark for performing more robust model-data comparison. The surface temperature simulated by two General Circulation Models (CCSM3 and HadCM3) for 130 ka and 125 ka is compared to the corresponding time slice data synthesis. This comparison shows that the models predict warmer than present conditions earlier than documented in the North Atlantic, while neither model is able to produce the reconstructed early Southern Ocean and Antarctic warming. Our results highlight the importance of producing a sequence of time slices rather than one single time slice averaging the LIG climate conditions.

Modeled monthly local topographic lapse rates, Pacific Northwest, USA

An emerging approach to downscaling the projections from General Circulation Models (GCMs) to scales relevant for basin hydrology is to use output of GCMs to force higher-resolution Regional Climate Models (RCMs). With spatial resolution often in the tens of kilometers, however, even RCM output will likely fail to resolve local topography that may be climatically significant in high-relief basins. Here we develop and apply an approach for downscaling RCM output using local topographic lapse rates (empirically-estimated spatially and seasonally variable changes in climate variables with elevation). We calculate monthly local topographic lapse rates from the 800-m Parameter-elevation Regressions on Independent Slopes Model (PRISM) dataset, which is based on regressions of observed climate against topographic variables. We then use these lapse rates to elevationally correct two sources of regional climate-model output: (1) the North American Regional Reanalysis (NARR), a retrospective dataset produced from a regional forecasting model constrained by observations, and (2) a range of baseline climate scenarios from the North American Regional Climate Change Assessment Program (NARCCAP), which is produced by a series of RCMs driven by GCMs. By running a calibrated and validated hydrologic model, the Soil and Water Assessment Tool (SWAT), using observed station data and elevationally-adjusted NARR and NARCCAP output, we are able to estimate the sensitivity of hydrologic modeling to the source of the input climate data. Topographic correction of regional climate-model data is a promising method for modeling the hydrology of mountainous basins for which no weather station datasets are available or for simulating hydrology under past or future climates.

Response of A. aurita ephyrae to three stressors: high-temperature, high-CO2 and low-oxygen

Global change is affecting marine ecosystems through a combination of different stressors such as warming, ocean acidification and oxygen depletion. Very little is known about the interactions among these factors, especially with respect to gelatinous zooplankton. Therefore, in this study we investigated the direct effects of pH, temperature and oxygen availability on the moon jellyfish Aurelia aurita, concentrating on the ephyral life stage. Starved one-day-old ephyrae were exposed to a range of pCO2 (400-4000 ppm) and three different dissolved oxygen levels (from saturated to hypoxic conditions), in two different temperatures (5 and 15 °C) for 7 days. Carbon content and swimming activity were analysed at the end of the incubation period, and mortality noted. General linearized models were fitted through the data, with the best fitting models including two- and three-way interactions between pCO2, temperature and oxygen concentration. The combined effect of the stressors was small but significant, with the clearest negative effect on growth caused by the combination of all three stressors present (high temperature, high CO2, low oxygen). We conclude that A. aurita ephyrae are robust and that they are not likely to suffer from these environmental stressors in a near future.

Bomb-tritium input function calculated from the details of the nuclear atmospheric bomb tests released by the UNSCEAR [2000]

Improving the representation of the hydrological cycle in Atmospheric General Circulation Models (AGCMs) is one of the main challenges in modeling the Earth's climate system. One way to evaluate model performance is to simulate the transport of water isotopes. Among those available, tritium (HTO) is an extremely valuable tracer, because its content in the different reservoirs involved in the water cycle (stratosphere, troposphere, ocean) varies by order of magnitude. Previous work incorporated natural tritium into LMDZ-iso, a version of the LMDZ general circulation model enhanced by water isotope diagnostics. Here for the first time, the anthropogenic tritium injected by each of the atmospheric nuclear-bomb tests between 1945 and 1980 has been first estimated and further implemented in the model; it creates an opportunity to evaluate certain aspects of LDMZ over several decades by following the bomb-tritium transient signal through the hydrological cycle. Simulations of tritium in water vapor and precipitation for the period 1950-2008, with both natural and anthropogenic components, are presented in this study. LMDZ-iso satisfactorily reproduces the general shape of the temporal evolution of tritium. However, LMDZ-iso simulates too high a bomb-tritium peak followed by too strong a decrease of tritium in precipitation. The too diffusive vertical advection in AGCMs crucially affects the residence time of tritium in the stratosphere. This insight into model performance demonstrates that the implementation of tritium in an AGCM provides a new and valuable test of the modeled atmospheric transport, complementing water stable isotope modeling.

Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo 1.0.5-beta)

Sub-ice shelf circulation and freezing/melting rates in ocean general circulation models depend critically on an accurate and consistent representation of cavity geometry. Existing global or pan-Antarctic data sets have turned out to contain various inconsistencies and inaccuracies. The goal of this work is to compile independent regional fields into a global data set. We use the S-2004 global 1-minute bathymetry as the backbone and add an improved version of the BEDMAP topography for an area that roughly coincides with the Antarctic continental shelf. Locations of the merging line have been carefully adjusted in order to get the best out of each data set. High-resolution gridded data for upper and lower ice surface topography and cavity geometry of the Amery, Fimbul, Filchner-Ronne, Larsen C and George VI Ice Shelves, and for Pine Island Glacier have been carefully merged into the ambient ice and ocean topographies. Multibeam survey data for bathymetry in the former Larsen B cavity and the southeastern Bellingshausen Sea have been obtained from the data centers of Alfred Wegener Institute (AWI), British Antarctic Survey (BAS) and Lamont-Doherty Earth Observatory (LDEO), gridded, and again carefully merged into the existing bathymetry map. The global 1-minute dataset (RTopo-1 Version 1.0.5) has been split into two netCDF files. The first contains digital maps for global bedrock topography, ice bottom topography, and surface elevation. The second contains the auxiliary maps for data sources and the surface type mask. A regional subset that covers all variables for the region south of 50 deg S is also available in netCDF format. Datasets for the locations of grounding and coast lines are provided in ASCII format.