Diatom abundance in surface sediments of the northern North Pacific

In order to map the modern distribution of diatoms and to establish a reliable reference data set for paleoenvironmental reconstruction in the northern North Pacific, a new data set including the relative abundance of diatom species preserved in a total of 422 surface sediments was generated, which covers a broad range of environmental variables characteristic of the subarctic North Pacific, the Sea of Okhotsk and the Bering Sea between 30° and 70°N. The biogeographic distribution patterns as well as the preferences in sea surface temperature of 38 diatom species and species groups are documented. A Q-mode factor analysis yields a three-factor model representing assemblages associated with the Arctic, Subarctic and Subtropical water mass, indicating a close relationship between the diatom composition and the sea surface temperatures. The relative abundance pattern of 38 diatom species and species groups was statistically compared with nine environmental variables, i.e. the summer sea surface temperature and salinity, annual surface nutrient concentration (nitrate, phosphate, silicate), summer and winter mixed layer depth and summer and winter sea ice concentrations. Canonical Correspondence Analysis (CCA) indicates 32 species and species groups have strong correspondence with the pattern of summer sea surface temperature. In addition, the total diatom flux data compiled from ten sediment traps reveal that the seasonal signals preserved in the surface sediments are mostly from spring through autumn. This close relationship between diatom composition and the summer sea surface temperature will be useful in deriving a transfer function in the subarctic North Pacific for the quantitative paleoceanographic and paleoenvironmental studies. The relative abundance of the sea-ice indicator diatoms Fragilariopsis cylindrus and F. oceanica of >20% in the diatom composition is used to represent the winter sea ice edge in the Bering Sea. The northern boundary of the distribution of F. doliolus in the open ocean is suggested to be an indicator of the Subarctic Front, while the abundance of Chaetoceros resting spores may indicate iron input from nearby continents and shelves and induced productivity events in the study area.

Diatom abundance in surface sediments of the Southern Ocean

Based on the quantitative analysis of diatom assemblages preserved in 274 surface sediment samples recovered in the Pacific, Atlantic and western Indian sectors of the Southern Ocean we have defined a new reference database for quantitative estimation of late-middle Pleistocene Antarctic sea ice fields using the transfer function technique. The Detrended Canonical Analysis (DCA) of the diatom data set points to a unimodal distribution of the diatom assemblages. Canonical Correspondence Analysis (CCA) indicates that winter sea ice (WSI) but also summer sea surface temperature (SSST) represent the most prominent environmental variables that control the spatial species distribution. To test the applicability of transfer functions for sea ice reconstruction in terms of concentration and occurrence probability we applied four different methods, the Imbrie and Kipp Method (IKM), the Modern Analog Technique (MAT), Weighted Averaging (WA), and Weighted Averaging Partial Least Squares (WAPLS), using logarithm-transformed diatom data and satellite-derived (1981-2010) sea ice data as a reference. The best performance for IKM results was obtained using a subset of 172 samples with 28 diatom taxa/taxa groups, quadratic regression and a three-factor model (IKM-D172/28/3q) resulting in root mean square errors of prediction (RMSEP) of 7.27% and 11.4% for WSI and summer sea ice (SSI) concentration, respectively. MAT estimates were calculated with different numbers of analogs (4, 6) using a 274-sample/28-taxa reference data set (MAT-D274/28/4an, -6an) resulting in RMSEP's ranging from 5.52% (4an) to 5.91% (6an) for WSI as well as 8.93% (4an) to 9.05% (6an) for SSI. WA and WAPLS performed less well with the D274 data set, compared to MAT, achieving WSI concentration RMSEP's of 9.91% with WA and 11.29% with WAPLS, recommending the use of IKM and MAT. The application of IKM and MAT to surface sediment data revealed strong relations to the satellite-derived winter and summer sea ice field. Sea ice reconstructions performed on an Atlantic- and a Pacific Southern Ocean sediment core, both documenting sea ice variability over the past 150,000 years (MIS 1 - MIS 6), resulted in similar glacial/interglacial trends of IKM and MAT-based sea-ice estimates. On the average, however, IKM estimates display smaller WSI and slightly higher SSI concentration and probability at lower variability in comparison with MAT. This pattern is a result of different estimation techniques with integration of WSI and SSI signals in one single factor assemblage by applying IKM and selecting specific single samples, thus keeping close to the original diatom database and included variability, by MAT. In contrast to the estimation of WSI, reconstructions of past SSI variability remains weaker. Combined with diatom-based estimates, the abundance and flux pattern of biogenic opal represents an additional indication for the WSI and SSI extent.

Southern Ocean CO2 uptake and buffer factor from model run 2012-2100, with links to model results

The Southern Ocean is a key region for global carbon uptake and is characterised by a strong seasonality with the annual CO2 uptake being mediated by biological carbon draw-down in summer. Here, we show that the contribution of biology to CO2 uptake will become even more important until 2100. This is the case even if biological production remains unaltered and can be explained by the decreasing buffer capacity of the ocean as its carbon content increases. The same amount of biological carbon draw-down leads to a more than twice as large reduction in CO2 (aq) concentration and hence to a larger CO2 gradient between ocean and atmosphere that drives the gas-exchange. While the winter uptake south of 44°S changes little, the summer uptake increases largely and is responsible for the annual mean response. The combination of decreasing buffer capacity and strong seasonality of biological carbon draw-down introduces a strong and increasing seasonality in the anthropogenic carbon uptake.

Ocean temperature response to idealized Gleissberg and de Vries solar cycles in a comprehensive climate model

The ~90-year Gleissberg and ~200-year de Vries cycles have been identified as two distinctive quasi-periodic components of Holocene solar activity. Evidence exists for the impact of such multi-decadal to centennial-scale variability in total solar irradiance (TSI) on climate, but concerning the ocean, this evidence is mainly restricted to the surface response. Here we use a comprehensive global climate model to study the impact of idealized solar forcing, representing the Gleissberg and de Vries cycles, on global ocean potential temperature at different depth levels, after a recent proxy record indicates a signal of TSI anomalies in the northeastern Atlantic at mid-depth. Potential impacts of TSI anomalies on deeper oceanic levels are climatically relevant due to their possible effect on ocean circulation by altering water mass characteristics. Simulated solar anomalies are shown to penetrate the ocean down to at least deep-water levels. Despite the fact that the two forcing periods differ only by a factor of ~2, the spatial pattern of response is significantly distinctive between the experiments, suggesting different mechanisms for solar signal propagation. These are related to advection by North Atlantic Deep Water flow (200-year forcing), and barotropic adjustment in the South Atlantic in response to a latitudinal shift of the westerly wind belt (90-year forcing).

Microcharcoal concentration and elongation, and age model of sediment core MD96-2098

Although grassland and savanna occupy only a quarter of the world's vegetation, burning in these ecosystems accounts for roughly half the global carbon emissions from fire. However, the processes that govern changes in grassland burning are poorly understood, particularly on time scales beyond satellite records. We analyzed microcharcoal, sediments, and geochemistry in a high-resolution marine sediment core off Namibia to identify the processes that have controlled biomass burning in southern African grassland ecosystems under large, multimillennial-scale climate changes. Six fire cycles occurred during the past 170,000 y in southern Africa that correspond both in timing and magnitude to the precessional forcing of north-south shifts in the Intertropical Convergence Zone. Contrary to the conventional expectation that fire increases with higher temperatures and increased drought, we found that wetter and cooler climates cause increased burning in the study region, owing to a shift in rainfall amount and seasonality (and thus vegetation flammability). We also show that charcoal morphology (i.e., the particle's length-to-width ratio) can be used to reconstruct changes in fire activity as well as biome shifts over time. Our results provide essential context for understanding current and future grassland-fire dynamics and their associated carbon emissions.

Mineral dust deposition in interglacial and glacial climate conditions: model results

The global aerosol/climate model ECHAM5-HAM is used in order to investigate the dust cycle for four interglacial and one glacial climate conditions. The 20-year time-slices are the pre-industrial control (CTRL), mid-Holocene (6000 years BP), last glacial inception (115000 years BP), Eemian (126000 years BP) and Last Glacial Maximum (LGM) (21000 years BP) time intervals. The study is focused on the Antarctic region. The model is able to reproduce the magnitude order of dust deposition globally for the pre-industial and LGM climates. Correlation coefficient of the natural logarithm of the observed and modeled values is 0.78 for the CTRL and 0.81 for the LGM. For the pre-industrial simulation the model overestimates observed values in Antarctica by a factor of about 2-3 due to overestimation of the Australian dust source and too high wet deposition in the Antarctica interior. In the LGM, the model underestimates dust deposition in eastern Antarctica by a factor of about 4-5 due to underestimation of the South American dust source. More records are needed to validate dust deposition for the past interglacial time-slices. The modeled results show that dust deposition in Antarctica in the past interglacial time-slices is higher than in the CTRL simulation. The largest increase of dust deposition in Antarctica is simulated for the LGM, showing about 10-fold increase compared to CTRL.

Model simulation output data (CCSM3) related to uplift of African topography and the intensification of upwelling in the Benguela region

The Benguela Current, located off the west coast of southern Africa, is tied to a highly productive upwelling system**1. Over the past 12 million years, the current has cooled, and upwelling has intensified**2, 3, 4. These changes have been variously linked to atmospheric and oceanic changes associated with the glaciation of Antarctica and global cooling**5, the closure of the Central American Seaway**1, 6 or the further restriction of the Indonesian Seaway**3. The upwelling intensification also occurred during a period of substantial uplift of the African continent**7, 8. Here we use a coupled ocean-atmosphere general circulation model to test the effect of African uplift on Benguela upwelling. In our simulations, uplift in the East African Rift system and in southern and southwestern Africa induces an intensification of coastal low-level winds, which leads to increased oceanic upwelling of cool subsurface waters. We compare the effect of African uplift with the simulated impact of the Central American Seaway closure9, Indonesian Throughflow restriction10 and Antarctic glaciation**11, and find that African uplift has at least an equally strong influence as each of the three other factors. We therefore conclude that African uplift was an important factor in driving the cooling and strengthening of the Benguela Current and coastal upwelling during the late Miocene and Pliocene epochs.

Age model of ODP Site 162-983

At ODP Site 983, relative geomagnetic paleointensity and planktic and benthic delta18O records have been acquired for the last 350 kyr. The mean sedimentation rate in this interval is 11.3 cm/kyr. Magnetic properties and hysteresis ratios indicate that pseudo-single domain magnetite is the remanence carrier. Volume susceptibility (kappa), anhysteretic (ARM) and isothermal (IRM) remanence values vary by a factor of 3-4, well within the criteria usually cited for paleointensity studies. Natural remanent magnetization (NRM) is normalized by ARM and IRM to acquire the paleointensity proxy. Arithmetic means of NRM/ARM and NRM/IRM, calculated for five demagnetization steps in the 25-45 mT range, constitute the relative paleointensity estimates. Some paleointensity lows (particularly those at ~40, ~120 and ~188 ka) are associated with directional excursions of the field, especially the event at ~188 ka (referred to here as the Iceland Basin Event) that constitutes a short-lived polarity reversal. For the last 200 kyr, the records can be correlated with other high-resolution paleointensity records such as those from the Labrador Sea, Mediterranean/Somali Basin and Sulu Sea, implying that the millennial scale features are globally synchronous. A labeling system for paleointensity features is proposed that ties prominent highs and lows to oxygen isotope stages.

Model of interactive visualization of diffusion processes from different strains, and source routes and spatial diffusion of capuchin monkeys over the past 6 million years

Source routes and Spatial Diffusion of capuchin monkeys over the past 6 million years, rebuilt in the SPREAD 1.0.6 from the MCC tree. The map shows the 10 different regions to which distinctive samples were associated. The different transmission routes have been calculated from the average rate over time. Only rates with Bayes factor> 3 were considered as significantly different from zero. Significant diffusion pathways are highlighted with color varying from dark brown to red, being the dark brown less significant rates and deep red the most significant rates.

Coccolithophore taxa and PCA factor analysis of ODP Site 1233

We provide high-resolution sea surface temperature (SST) and paleoproductivity data focusing on Termination 1. We describe a new method for estimating SSTs based on multivariate statistical analyses performed on modern coccolithophore census data, and we present the first downcore reconstructions derived from coccolithophore assemblages at Ocean Drilling Project (ODP) Site 1233 located offshore Chile. We compare our coccolithophore SST record to alkenone-based SSTs as well as SST reconstructions based on dinoflagellates and radiolaria. All reconstructions generally show a remarkable concordance. As in the alkenone SST record, the Last Glacial Maximum (LGM, 19-23 kyr B.P.) is not clearly defined in our SST reconstruction. After the onset of deglaciation, three major warming steps are recorded: from 18.6 to 18 kyr B.P. (~2.6°C), from 15.7 to 15.3 kyr B.P. (~2.5°C), and from 13 to 11.4 kyr B.P. (~3.4°C). Consistent with the other records from Site 1233 and Antarctic ice core records, we observed a clear Holocene Climatic Optimum (HCO) from ~8-12 kyr B.P. Combining the SST reconstruction with coccolith absolute abundances and accumulation rates, we show that colder temperatures during the LGM are linked to higher coccolithophore productivity offshore Chile and warmer SSTs during the HCO to lower coccolithophore productivity, with indications of weak coastal upwelling. We interpret our data in terms of latitudinal displacements of the Southern Westerlies and the northern margin of the Antarctic Circumpolar Current system over the deglaciation and the Holocene.

Clay mineralogical assemblages, Sm-Nd concentrations, Sm and Nd isotopic ratios and oxygen isotope ratios, OIS number and the age model of ODP Site 105-646

We present 40 Sm-Nd isotope measurements of the clay-size (<2 µm) fractions of sediments from the Southern Greenland rise (ODP-646) that span the last 365 kyr. These data track changes in the relative supply of fine particles carried into the deep Labrador Sea by the Western Boundary Under Current (WBUC) back to the fourth glacial-interglacial cycles. Earlier studies revealed three general sources of particles to the core site: (i) Precambrian crustal material from Canada, Greenland, and/or Scandinavia (North American Shield - NAS), (ii) Palaeozoic or younger crustal material from East Greenland, NW Europe, and/or western Scandinavia (Young Crust - YC) and (iii) volcanic material from Iceland and the Mid-Atlantic Ridge (MAR). Clay-size fractions from glacial sediments have the lowest Nd isotopic ratios. Supplies of young crustal particles were similar during glacial oxygen isotope stages (OIS) 2, 6, and 10. In contrast the mean volcanic contributions decreased relative to old craton material from OIS 10 to OIS 6 and then from OIS 6 to OIS 2. The glacial OIS 8 interval displays a mean Sm/Nd ratio similar to those of interglacials OIS 1, 5, and 9. Compared with other interglacials, OIS 7 was marked by a higher YC contribution but a similar ~30% MAR supply. The overall NAS contribution dropped by a factor of 2 during each glacial/interglacial transition, with the MAR contribution broadly replacing it during interglacials. To decipher between higher supplies and/or dilution, particle fluxes from each end member were estimated. Glacial NAS fluxes were systematically higher than interglacial fluxes. During the time interval examined, fine particle supplies to the Labrador Sea were strongly controlled by proximal ice-margin erosion and thus echoed the glacial stage intensity. In contrast, the WBUC-carried MAR supplies from the eastern basins did not change significantly throughout the last 365 kyr, except for a marked increase in surface-sediments that suggests unique modern conditions. Distal WBUC-controlled inputs from the Northern and NE North Atlantic seem to have been less variable than proximal supplies linked with glacial erosion rate.