Environmental controls on the Emiliania huxleyi calcite mass

Although ocean acidification is expected to impact (bio)calcification by decreasing the seawater carbonate ion concentration, [CO3]2-, there exists evidence of non-uniform response of marine calcifying plankton to low seawater [CO3]2-. This raises questions on the role of environmental factors other than acidification and on the complex physiological responses behind calcification. Here we investigate the synergistic effect of multiple environmental parameters, including temperature, nutrient (nitrate and phosphate) availability, and seawater carbonate chemistry on the coccolith calcite mass of the cosmopolitan coccolithophore Emiliania huxleyi, the most abundant species in the world ocean. We use a suite of surface (late Holocene) sediment samples from the South Atlantic and southwestern Indian Ocean taken from depths lying well above the modern lysocline. The coccolith calcite mass in our results presents a latitudinal distribution pattern that mimics the main oceanographic features, thereby pointing to the potential importance of phosphorus and temperature in determining coccolith mass by affecting primary calcification and possibly driving the E. huxleyi morphotype distribution. This evidence does not necessarily argue against the potentially important role of the rapidly changing seawater carbonate chemistry in the future, when unabated fossil fuel burning will likely perturb ocean chemistry beyond a critical point. Rather our study highlights the importance of evaluating the combined effect of several environmental stressors on calcifying organisms to project their physiological response(s) in a high CO2 world and improve interpretation of paleorecords.

Alkenones and SST of surface sediments (Appendix)

We have analysed alkenones in 149 surface sediments from the eastern South Atlantic in order to establish a sediment-based calibration of the U37K' paleotemperature index. Our study covers the major tropical to subpolar production systems and sea-surface temperatures (SST's) between 0° and 27°C. In order to define the most suitable calibration for this region, the U37K' values were correlated to seasonal, annual, and production-weighted annual mean atlas temperatures and compared to previously published culture and core-top calibrations. The best linear correlation between U37K' and SST was obtained using annual mean SST from 0 to 10 m water depth (U37K' = 0.033 T + 0.069, r**2 = 0.981). Data scattering increased significantly using temperatures of waters deeper than 20 m, suggesting that U37K' reflects mixed-layer SST and that alkenone production at thermocline depths was not high enough to significantly bias the mixed-layer signal. Regressions based on both production-weighted and on actual annual mean atlas SST were virtually identical, indicating that regional variations in the seasonality of primary production have no discernible effect on the U37K' vs. SST relationship. Comparison with published core-top calibrations from other oceanic regions revealed a high degree of accordance. We, therefore, established a global core-top calibration using U37K' data from 370 sites between 60°S and 60°N in the Atlantic, Indian, and Pacific Oceans and annual mean atlas SST (0-29°C) from 0 m water depth. The resulting relationship (U37K' = 0.033 T + 0.044, r**2 = 958) is identical within error limits to the widely used E. huxleyi calibrations of and attesting their general applicability. The observation that core-top calibrations extending over various biogeographical coccolithophorid zones are strongly linear and in better accordance than culture calibrations suggests that U37K' is less species-dependent than is indicated by culture experiments. The results also suggest that variations in growth rate of algae and nutrient availability do not significantly affect the sedimentary record of U37K' in open ocean environments.

Global compilation of benthic data sets I

Approaches to quantify the organic carbon accumulation on a global scale generally do not consider the small-scale variability of sedimentary and oceanographic boundary conditions along continental margins. In this study, we present a new approach to regionalize the total organic carbon (TOC) content in surface sediments (<5 cm sediment depth). It is based on a compilation of more than 5500 single measurements from various sources. Global TOC distribution was determined by the application of a combined qualitative and quantitative-geostatistical method. Overall, 33 benthic TOC-based provinces were defined and used to process the global distribution pattern of the TOC content in surface sediments in a 1°x1° grid resolution. Regional dependencies of data points within each single province are expressed by modeled semi-variograms. Measured and estimated TOC values show good correlation, emphasizing the reasonable applicability of the method. The accumulation of organic carbon in marine surface sediments is a key parameter in the control of mineralization processes and the material exchange between the sediment and the ocean water. Our approach will help to improve global budgets of nutrient and carbon cycles.

Compilation of planktonic foraminiferal census counts from the Southern Hemisphere Oceans

We present an improved database of planktonic foraminiferal census counts from the Southern Hemisphere Oceans (SHO) from 15°S to 64°S. The SHO database combines 3 existing databases. Using this SHO database, we investigated dissolution biases that might affect faunal census counts. We suggest a depth/[DCO3]2- threshold of ~3800 m/[DCO3]2- = ~-10 to -5 µmol/kg for the Pacific and Indian Oceans, and ~4000 m/[DCO3]2- = ~0 to 10 µmol/kg for the Atlantic Ocean, under which core-top assemblages can be affected by dissolution and are less reliable for paleo-sea surface temperature (SST) reconstructions. We removed all core-tops beyond these thresholds from the SHO database. This database has 598 core-tops and is able to reconstruct past SST variations from 2° to 25.5°C, with a root mean square error of 1.00°C, for annual temperatures. To inspect dissolution affects SST reconstruction quality, we tested the data base with two "leave-one-out" tests, with and without the deep core-tops. We used this database to reconstruct Summer SST (SSST) over the last 20 ka, using the Modern Analog Technique method, on the Southeast Pacific core MD07-3100. This was compared to the SSST reconstructed using the 3 databases used to compile the SHO database. Thus showing that the reconstruction using the SHO database is more reliable, as its dissimilarity values are the lowest. The most important aspect here is the importance of a bias-free, geographic-rich, database. We leave this dataset open-ended to future additions; the new core-tops must be carefully selected, with their chronological frameworks, and evidence of dissolution assessed.

Glacial North Atlantic: Sea surface conditions reconstructed by GLAMAP 2000

The response of the tropical ocean to global climate change and the extent of sea ice in the glacial nordic seas belong to the great controversies in paleoclimatology. Our new reconstruction of peak glacial sea surface temperatures (SSTs) in the Atlantic is based on census counts of planktic foraminifera, using the Maximum Similarity Technique Version 28 (SIMMAX-28) modern analog technique with 947 modern analog samples and 119 well-dated sediment cores. Our study compares two slightly different scenarios of the Last Glacial Maximum (LGM), the Environmental Processes of the Ice Age: Land, Oceans, Glaciers (EPILOG), and Glacial Atlantic Ocean Mapping (GLAMAP 2000) time slices. The comparison shows that the maximum LGM cooling in the Southern Hemisphere slightly preceeded that in the north. In both time slices sea ice was restricted to the north western margin of the nordic seas during glacial northern summer, while the central and eastern parts were ice-free. During northern glacial winter, sea ice advanced to the south of Iceland and Faeroe. In the central northern North Atlantic an anticyclonic gyre formed between 45° and 60°N, with a cool water mass centered west of Ireland, where glacial cooling reached a maximum of >12°C. In the subtropical ocean gyres the new reconstruction supports the glacial-to-interglacial stability of SST as shown by CLIMAP Project Members (CLIMAP) [1981]. The zonal belt of minimum SST seasonality between 2° and 6°N suggests that the LGM caloric equator occupied the same latitude as today. In contrast to the CLIMAP reconstruction, the glacial cooling of the tropical east Atlantic upwelling belt reached up to 6°-8°C during Northern Hemisphere summer. Differences between these SIMMAX-based and published U37[k]- and Mg/Ca-based equatorial SST records are ascribed to strong SST seasonalities and SST signals that were produced by different planktic species groups during different seasons.

Distribution of planktic foraminifera in surface sediments of the Atlantic Ocean

We present a data set of 738 planktonic foraminiferal species counts from sediment surface samples of the eastern North Atlantic and the South Atlantic between 87°N and 40°S, 35°E and 60°W including published Climate: Long-Range Investigation, Mapping, and Prediction (CLIMAP) data. These species counts are linked to Levitus's [1982] modern water temperature data for the four caloric seasons, four depth ranges (0, 30, 50, and 75 m), and the combined means of those depth ranges. The relation between planktonic foraminiferal assemblages and sea surface temperature (SST) data is estimated using the newly developed SIMMAX technique, which is an acronym for a modern analog technique (MAT) with a similarity index, based on (1) the scalar product of the normalized faunal percentages and (2) a weighting procedure of the modern analog's SSTs according to the inverse geographical distances of the most similar samples. Compared to the classical CLIMAP transfer technique and conventional MAT techniques, SIMMAX provides a more confident reconstruction of paleo-SSTs (correlation coefficient is 0.994 for the caloric winter and 0.993 for caloric summer). The standard deviation of the residuals is 0.90°C for caloric winter and 0.96°C for caloric summer at 0-m water depth. The SST estimates reach optimum stability (standard deviation of the residuals is 0.88°C) at the average 0- to 75-m water depth. Our extensive database provides SST estimates over a range of -1.4 to 27.2°C for caloric winter and 0.4 to 28.6°C for caloric summer, allowing SST estimates which are especially valuable for the high-latitude Atlantic during glacial times.

Compilation of 220 sediment core d13C data from the glacial Atlantic Ocean

We compare a compilation of 220 sediment core d13C data from the glacial Atlantic Ocean with three-dimensional ocean circulation simulations including a marine carbon cycle model. The carbon cycle model employs circulation fields which were derived from previous climate simulations. All sediment data have been thoroughly quality controlled, focusing on epibenthic foraminiferal species (such as Cibicidoides wuellerstorfi or Planulina ariminensis) to improve the comparability of model and sediment core carbon isotopes. The model captures the general d13C pattern indicated by present-day water column data and Late Holocene sediment cores but underestimates intermediate and deep water values in the South Atlantic. The best agreement with glacial reconstructions is obtained for a model scenario with an altered freshwater balance in the Southern Ocean that mimics enhanced northward sea ice export and melting away from the zone of sea ice production. This results in a shoaled and weakened North Atlantic Deep Water flow and intensified Antarctic Bottom Water export, hence confirming previous reconstructions from paleoproxy records. Moreover, the modeled abyssal ocean is very cold and very saline, which is in line with other proxy data evidence.

Global data compilation of benthic data sets II

In this study we present a global distribution pattern and budget of the minimum flux of particulate organic carbon to the sea floor (J POC alpha). The estimations are based on regionally specific correlations between the diffusive oxygen flux across the sediment-water interface, the total organic carbon content in surface sediments, and the oxygen concentration in bottom waters. For this, we modified the principal equation of Cai and Reimers [1995] as a basic monod reaction rate, applied within 11 regions where in situ measurements of diffusive oxygen uptake exist. By application of the resulting transfer functions to other regions with similar sedimentary conditions and areal interpolation, we calculated a minimum global budget of particulate organic carbon that actually reaches the sea floor of ~0.5 GtC yr**-1 (>1000 m water depth (wd)), whereas approximately 0.002-0.12 GtC yr**-1 is buried in the sediments (0.01-0.4% of surface primary production). Despite the fact that our global budget is in good agreement with previous studies, we found conspicuous differences among the distribution patterns of primary production, calculations based on particle trap collections of the POC flux, and J POC alpha of this study. These deviations, especially located at the southeastern and southwestern Atlantic Ocean, the Greenland and Norwegian Sea and the entire equatorial Pacific Ocean, strongly indicate a considerable influence of lateral particle transport on the vertical link between surface waters and underlying sediments. This observation is supported by sediment trap data. Furthermore, local differences in the availability and quality of the organic matter as well as different transport mechanisms through the water column are discussed.

LGM, planktic foraminifera (EPILOG version)

[1] We used planktic foraminiferal assemblages in 70 sediment cores from the tropical and subtropical South Atlantic Ocean (10°N-37°S) to estimate annual mean sea surface temperatures (SSTs) and seasonality for the Last Glacial Maximum with a modified version of the Imbrie-Kipp transfer function method (IKTF) that takes into account the abundance of rare but temperature sensitive species. In contrast to CLIMAP Project Members [1981], the reconstructed SSTs indicate cooler glacial SSTs in the entire tropical/subtropical South Atlantic with strongest cooling in the upwelling region off Namibia (7-10°C) and smallest cooling (1-2°C) in the western subtropical gyre. In the western Atlantic, our data support recent temperature estimates from other proxies. In the upwelling regions in the eastern Atlantic, our data conflict with SST reconstructions from alkenones, which may be due to an environmental preference of the alkenone-producing algae or to an underestimation of foraminiferal SSTs due to anomalous high abundances of N. pachyderma (sinistral).

Dissolution index of Globigerina bulloides in recent and Last Glacial Maximum sediments

The modern Atlantic Ocean, dominated by the interactions of North Atlantic Deep Water (NADW) and Antarctic Bottom Water (AABW), plays a key role in redistributing heat from the Southern to the Northern Hemisphere. In order to reconstruct the evolution of the relative importance of these two water masses, the NADW/AABW transition, reflected by the calcite lysocline, was investigated by the Globigerina bulloides dissolution index (BDX?). The depth level of the Late Glacial Maximum (LGM) calcite lysocline was elevated by several hundred metres, indicating a more corrosive water mass present at modern NADW level. Overall, the small range of BDX? data and the gradual decrease in preservation below the calcite lysocline point to a less stratified Atlantic Ocean during the LGM. Similar preservation patterns in the West and East Atlantic demonstrate that the modern west-east asymmetry did not exist due to an expansion of southern deep waters compensating for the decrease in NADW formation.