EOT11a - Global Empirical Ocean Tide model from multi-mission satellite altimetry, with links to model results

EOT11a is a global (E)mpirical (O)cean (T)ide model derived in 2011 by residual analysis of multi-mission satellite (a)ltimeter data. EOT11a includes amplitudes and phases of the main astronomical tides M2, S2, N2, K2, 2N2, O1, K1, P2, and Q1, the non-linear constituent M4, the long period tides Mm and Mf, and the radiational tide S1. Ocean tides as well as loading tides are provided. EOT11a was computed by means of residual tidal analysis of multi-mission altimeter data from TOPEX/Poseidon, ERS-2, ENVISAT, and Jason-1/2, as far as acquired between September 1992 and April 2010. The resolution of 7.5'x7.5' is identical with FES2004 which was used as reference model for the residual tide analysis. The development of EOT11a was funded by the Deutsche Forschungsgemeinschaft (DFG) under grant BO1228/6-2.

Marine Isotope Stage (MIS) 11 results of model data with Community Climate System Model version 3 including the Dynamic Global Vegetation Model (CCSM3-DGVM) in NetCDF format at global and regional scale

The climate of Marine Isotope Stage (MIS) 11, the interglacial roughly 400,000 years ago, is investigated for four time slices, 416, 410, 400, and 394 ka. The overall picture is that MIS 11 was a relatively warm interglacial in comparison to preindustrial, with Northern Hemisphere (NH) summer temperatures early in MIS 11 (416-410 ka) warmer than preindustrial, though winters were cooler. Later in MIS 11, especially around 400 ka, conditions were cooler in the NH summer, mainly in the high latitudes. Climate changes simulated by the models were mainly driven by insolation changes, with the exception of two local feedbacks that amplify climate changes. Here, the NH high latitudes, where reductions in sea ice cover lead to a winter warming early in MIS 11, as well as the tropics, where monsoon changes lead to stronger climate variations than one would expect on the basis of latitudinal mean insolation change alone, are especially prominent. The results support a northward expansion of trees at the expense of grasses in the high northern latitudes early during MIS 11, especially in northern Asia and North America.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Global distributions of diatoms abundance, biovolume and biomass - Gridded data product (NetCDF) - Contribution to the MAREDAT World Ocean Atlas of Plankton Functional Types

This study is a first effort to compile the largest possible body of data available from different plankton databases as well as from individual published or unpublished datasets regarding diatom distribution in the world ocean. The data obtained originate from time series studies as well as spatial studies. This effort is supported by the Marine Ecosystem Data (MAREDAT) project, which aims at building consistent data sets for the main PFTs (Plankton Functional Types) in order to help validate biogeochemical ocean models by using converted C biomass from abundance data. Diatom abundance data were obtained from various research programs with the associated geolocation and date of collection, as well as with a taxonomic information ranging from group down to species. Minimum, maximum and average cell size information were mined from the literature for each taxonomic entry, and all abundance data were subsequently converted to biovolume and C biomass using the same methodology.

Modelled global ocean temperature from 6 ka to present

We compare the ocean temperature evolution of the Holocene as simulated by climate models and reconstructed from marine temperature proxies. This site provides informations about the Holocene temperature trends as simulated by the models. We use transient simulations from a coupled atmosphere-ocean general circulation model, as well as an ensemble of time slice simulations from the Paleoclimate Modelling Intercomparison Project. The general pattern of sea surface temperature (SST) in the models shows a high latitude cooling and a low latitude warming. The proxy dataset comprises a global compilation of marine alkenone- and Mg/Ca-derived SST estimates. Independently of the choice of the climate model, we observe significant mismatches between modelled and estimated SST amplitudes in the trends for the last 6000 years. Alkenone-based SST records show a similar pattern as the simulated annual mean SSTs, but the simulated SST trends underestimate the alkenone-based SST trends by a factor of two to five. For Mg/Ca, no significant relationship between model simulations and proxy reconstructions can be detected. We tested if such discrepancies can be caused by too simplistic interpretations of the proxy data. We tested different seasons and depths in the model to compare the proxy data trends, and can reconcile only part of the mismatches on a regional scale. We therefore considered the additional environmental factor changes in the planktonic organisms' habitat depth and a time-shift in the recording season to diagnose whether invoking those environmental factors can help reconciling the proxy records and the model simulations. We find that invoking shifts in the living season and habitat depth can remove some of the model-data discrepancies in SST trends. Regardless whether such adjustments in the environmental parameters during the Holocene are realistic, they indicate that when modeled temperature trends are set up to allow drastic shifts in the ecological behavior of planktonic organisms, they do not capture the full range of reconstructed SST trends. Our findings indicate that climate model and reconstructed temperature trends are to a large degree only qualitatively comparable, thus providing a challenge for the interpretation of proxy data as well as the models' sensitivity to orbital forcing.

Global glacier changes

Changes in glaciers and ice caps provide some of the clearest evidence of climate change, and as such they constitute key variables for early detection strategies in global climate-related observations. These changes have impacts on global sea level fluctuations, the regional to local natural hazard situation, as well as on societies dependent on glacier meltwater. Internationally coordinated collection and publication of standardised information about ongoing glacier changes was initiated back in 1894. The compiled data sets on the global distribution and changes in glaciers and ice caps provide the backbone of the numerous scientific publications on the latest findings about surface ice on land. Since the very beginning, the compiled data has been published by the World Glacier Monitoring Service and its predecessor organisations. However, the corresponding data tables, formats and meta-data are mainly of use to specialists.

The early Eocene global warming events: Implications from ODP Leg 208 Walvis Ridge depth transect

Five sections drilled in multiple holes over a depth transect of more than 2200 m at the Walvis Ridge (SE Atlantic) during Ocean Drilling Program (ODP) Leg 208 resulted in the first complete early Paleogene deep-sea record. Here we present high-resolution stratigraphic records spanning a ~4.3 million yearlong interval of the late Paleocene to early Eocene. This interval includes the Paleocene-Eocene thermal maximum (PETM) as well as the Eocene thermal maximum (ETM) 2 event. A detailed chronology was developed with nondestructive X-ray fluorescence (XRF) core scanning records and shipboard color data. These records were used to refine the shipboard-derived spliced composite depth for each site and with a record from ODP Site 1051 were then used to establish a continuous time series over this interval. Extensive spectral analysis reveals that the early Paleogene sedimentary cyclicity is dominated by precession modulated by the short (100 kyr) and long (405 kyr) eccentricity cycles. Counting of precession-related cycles at multiple sites results in revised estimates for the duration of magnetochrons C24r and C25n. Direct comparison between the amplitude modulation of the precession component derived from XRF data and recent models of Earth's orbital eccentricity suggests that the onset of the PETM and ETM2 are related to a 100-kyr eccentricity maximum. Both events are approximately a quarter of a period offset from a maximum in the 405-kyr eccentricity cycle, with the major difference that the PETM is lagging and ETM2 is leading a 405-kyr eccentricity maximum. Absolute age estimates for the PETM, ETM2, and the magnetochron boundaries that are consistent with recalibrated radiometric ages and recent models of Earth's orbital eccentricity cannot be precisely determined at present because of too large uncertainties in these methods. Nevertheless, we provide two possible tuning options, which demonstrate the potential for the development of a cyclostratigraphic framework based on the stable 405-kyr eccentricity cycle for the entire Paleogene.

Part of the global DOC versus AOU (dissolved organic carbon/apparent oxygen utilization) data compilation, Azores-I

Temperate, transitional and subtropical waters of the remote Azores Front region east of Azores (24-40°N, 22-32°W) were sampled during three cruises conducted under increasing stratification conditions (April 1999, May 1997 and August 1998). Despite the temporal increase of surface temperature (by 5 °C) and stratification (by 2.1 1/min**2), as well as the thermocline shoaling (by ~15 m), dissolved organic carbon (DOC) and nitrogen (DON) in the surface layer were not significantly different for the early spring, late spring and summer periods, with average concentrations of 69±2 µM-C and 5.2±0.4 µM-N, respectively. The surface excess of semi-labile DOC, compared with the baseline DOC concentration in the deep ocean (47±2 µM-C), represents 33% of the bulk DOC concentration and as much as 85% of the TOC (=POC+DOC) excess. When compared with the winter baseline (56±2 µM-C), the seasonal surface DOC excess is 20% of the bulk DOC concentration and 87% of the seasonal TOC excess. These results confirm the major role played by DOC in the carbon cycle of surface waters of the Azores Front region. The total amount of bioreactive DOC transported from the temperate to the subtropical North Atlantic by the Ekman flux between March and December represents only ~15% of the average annual primary production, and ~15% and ~30% of the measured sinking POC flux+vertical DOC eddy diffusion during early spring and summer, respectively. Vertical eddy diffusion is 35% and 2% of the spring and summer sinking POC flux, respectively. On the other hand, DOC only contributes 13% to the local oxidation of organic matter in subsurface waters (between the pycnocline and 500 m) of the study region.

Global compilation of last glacial maximum oxygen isotopic ratios for planktonic and benthic foraminifera

We use the fully coupled atmosphere-ocean three-dimensional model of intermediate complexity iLOVECLIM to simulate the climate and oxygen stable isotopic signal during the Last Glacial Maximum (LGM, 21 000 yr). By using a model that is able to explicitly simulate the sensor (d18O), results can be directly compared with data from climatic archives in the different realms. Our results indicate that iLOVECLIM reproduces well the main feature of the LGM climate in the atmospheric and oceanic components. The annual mean d18O in precipitation shows more depleted values in the northern and southern high latitudes during the LGM. The model reproduces very well the spatial gradient observed in ice core records over the Greenland ice-sheet. We observe a general pattern toward more enriched values for continental calcite d18O in the model at the LGM, in agreement with speleothem data. This can be explained by both a general atmospheric cooling in the tropical and subtropical regions and a reduction in precipitation as confirmed by reconstruction derived from pollens and plant macrofossils. Data-model comparison for sea surface temperature indicates that iLOVECLIM is capable to satisfyingly simulate the change in oceanic surface conditions between the LGM and present. Our data-model comparison for calcite d18O allows investigating the large discrepancies with respect to glacial temperatures recorded by different microfossil proxies in the North Atlantic region. The results argue for a trong mean annual cooling between the LGM and present (>6°C), supporting the foraminifera transfer function reconstruction but in disagreement with alkenones and dinocyst reconstructions. The data-model comparison also reveals that large positive calcite d18O anomaly in the Southern Ocean may be explained by an important cooling, although the driver of this pattern is unclear. We deduce a large positive d18Osw anomaly for the north Indian Ocean that contrasts with a large negative d18Osw anomaly in the China Sea between the LGM and present. This pattern may be linked to changes in the hydrological cycle over these regions. Our simulation of the deep ocean suggests that changes in d18Osw between the LGM and present are not spatially homogenous. This is supported by reconstructions derived from pore fluids in deep-sea sediments. The model underestimates the deep ocean cooling thus biasing the comparison with benthic calcite d18O data. Nonetheless, our data-model comparison support a heterogeneous cooling of few degrees (2-4°C) in the LGM Ocean.

Gridded climate data from 5 Global Climate Models (GCM) of the Last Glacial Maximum (LGM) downscaled to 30 arc seconds for Europe, with links to NetCDF files

Studies on the impact of historical, current and future global change require very high-resolution climate data (less or equal 1km) as a basis for modelled responses, meaning that data from digital climate models generally require substantial rescaling. Another shortcoming of available datasets on past climate is that the effects of sea level rise and fall are not considered. Without such information, the study of glacial refugia or early Holocene plant and animal migration are incomplete if not impossible. Sea level at the last glacial maximum (LGM) was approximately 125m lower, creating substantial additional terrestrial area for which no current baseline data exist. Here, we introduce the development of a novel, gridded climate dataset for LGM that is both very high resolution (1km) and extends to the LGM sea and land mask. We developed two methods to extend current terrestrial precipitation and temperature data to areas between the current and LGM coastlines. The absolute interpolation error is less than 1°C and 0.5 °C for 98.9% and 87.8% of all pixels for the first two 1 arc degree distance zones. We use the change factor method with these newly assembled baseline data to downscale five global circulation models of LGM climate to a resolution of 1km for Europe. As additional variables we calculate 19 'bioclimatic' variables, which are often used in climate change impact studies on biological diversity. The new LGM climate maps are well suited for analysing refugia and migration during Holocene warming following the LGM.

Global and Regional Sea level budget components from GRACE and radar altimetry (2002-2014)

This dataset contains the result of a joint least squares inversion of GRACE and altimetry data. The results are evaluated in terms of sea level change for the global ocean as well as dedicated areas. In addition, some auxiliary data is provided to enable reproducibility of the results in Rietbroek et al. 2016, and a google Earth kmz file is provided which visualizes the trends derived from the inversion results.

Global distributions of coccolithophores abundance and biomass - Gridded data product (NetCDF) - Contribution to the MAREDAT World Ocean Atlas of Plankton Functional Types

Coccolithophores are calcifying marine phytoplankton of the class Prymnesiophyceae. They are considered to play an import role in the global carbon cycle through the production and export of organic carbon and calcite. We have compiled observations of global coccolithophore abundance from several existing databases as well as individual contributions of published and unpublished datasets. We estimate carbon biomass using standardised conversion methods and provide estimates of uncertainty associated with these values. The database contains 58 384 individual observations at various taxonomic levels. This corresponds to 12 391 observations of total coccolithophore abundance and biomass. The data span a time period of 1929-2008, with observations from all ocean basins and all seasons, and at depths ranging from the surface to 500 m. Highest biomass values are reported in the North Atlantic, with a maximum of 501.7 ?gCl-1. Lower values are reported for the Pacific (maximum of 79.4 ?gCl-1) and Indian Ocean (up to 178.3 ?gCl-1). Coccolithophores are reported across all latitudes in the Northern Hemisphere, from the Equator to 89degN, although biomass values fall below 3 ?gCl-1 north of 70degN. In the Southern Hemisphere, biomass values fall rapidly south of 50degS, with only a single non-zero observation south of 60degS. Biomass values show a clear seasonal cycle in the Northern Hemisphere, reaching a maximum in the summer months (June-July). In the Southern Hemisphere the seasonal cycle is less evident, possibly due to a greater proportion of low-latitude data.

Global distributions of Phaeocystis sp. abundance and biomass - Gridded data product (NetCDF) - Contribution to the MAREDAT World Ocean Atlas of Plankton Functional Types

The planktonic haptophyte Phaeocystis has been suggested to play a fundamental role in the global biogeochemical cycling of carbon and sulphur, but little is known about its global biomass distribution. We have collected global microscopy data of the genus Phaeocystis and converted abundance data to carbon biomass using species-specific carbon conversion factors. Microscopic counts of single-celled and colonial Phaeocystis were obtained both through the mining of online databases and by accepting direct submissions (both published and unpublished) from Phaeocystis specialists. We recorded abundance data from a total of 1595 depth-resolved stations sampled between 1955-2009. The quality-controlled dataset includes 5057 counts of individual Phaeocystis cells resolved to species level and information regarding life-stages from 3526 samples. 83% of stations were located in the Northern Hemisphere while 17% were located in the Southern Hemisphere. Most data were located in the latitude range of 50-70° N. While the seasonal distribution of Northern Hemisphere data was well-balanced, Southern Hemisphere data was biased towards summer months. Mean species- and form-specific cell diameters were determined from previously published studies. Cell diameters were used to calculate the cellular biovolume of Phaeocystis cells, assuming spherical geometry. Cell biomass was calculated using a carbon conversion factor for Prymnesiophytes (Menden-Deuer and Lessard, 2000). For colonies, the number of cells per colony was derived from the colony volume. Cell numbers were then converted to carbon concentrations. An estimation of colonial mucus carbon was included a posteriori, assuming a mean colony size for each species. Carbon content per cell ranged from 9 pg (single-celled Phaeocystis antarctica) to 29 pg (colonial Phaeocystis globosa). Non-zero Phaeocystis cell biomasses (without mucus carbon) range from 2.9 - 10?5 µg l-1 to 5.4 - 103 µg l-1, with a mean of 45.7 µg l-1 and a median of 3.0 µg l-1. Highest biomasses occur in the Southern Ocean below 70° S (up to 783.9 µg l-1), and in the North Atlantic around 50° N (up to 5.4 - 103 µg l-1).

Effects of ocean acidification and global warming on reef bioerosion-lessons from a clionaid sponge

Coral reefs are under threat, exerted by a number of interacting effects inherent to the present climate change, including ocean acidification and global warming. Bioerosion drives reef degradation by recycling carbonate skeletal material and is an important but understudied factor in this context. Twelve different combinations of pCO2 and temperature were applied to elucidate the consequences of ocean acidification and global warming on the physiological response and bioerosion rates of the zooxanthellate sponge Cliona orientalis-one of the most abundant and effective bioeroders on the Great Barrier Reef, Australia. Our results confirm a significant amplification of the sponges' bioerosion capacity with increasing pCO2, which is expressed by more carbonate being chemically dissolved by etching. The health of the sponges and their photosymbionts was not affected by changes in pCO2, in contrast to temperature, which had significant negative impacts at higher levels. However, we could not conclusively explain the relationship between temperature and bioerosion rates, which were slightly reduced at both colder as well as warmer temperatures than ambient. The present findings on the effects of ocean acidification on chemical bioerosion, however, will have significant implications for predicting future reef carbonate budgets, as sponges often contribute the lion's share of internal bioerosion on coral reefs.

(Table 1) CPI17-23, hopanoic acid, and n-C35 + lycopane/n-C31 ratios as well as 17beta(H), 21beta(H) homohopane, long-chain alkanes, and alkenones concentrations in IODP Exp302

During Integrated Ocean Drilling Program Expedition 302 (Arctic Coring Expedition (ACEX)) a more than 200 m thick sequence of Paleogene organic carbon (OC)-rich (black shale type) sediments was drilled. Here we present new biomarker data determined in ACEX sediment samples to decipher processes controlling OC accumulation and their paleoenvironmental significance during periods of Paleogene global warmth and proposed increased freshwater discharge in the early Cenozoic. Specific source-related biomarkers including n-alkanes, fatty acids, isoprenoids, carotenoids, hopanes/hopenes, hopanoic acids, aromatic terpenoids, and long-chain alkenones show a high variability of components, derived from marine and terrestrial origin. The distribution of hopanoic acid isomers is dominated by compounds with the biological 17beta(H), 21beta(H) configuration indicating a low level of maturity. On the basis of the biomarker data the terrestrial OC supply was significantly enriched during the late Paleocene and part of the earliest Eocene, whereas increased aquatic contributions and euxinic conditions of variable intensity were determined for the Paleocene-Eocene thermal maximum and Eocene thermal maximum 2 events as well as the middle Eocene time interval. Furthermore, samples from the middle Eocene are characterized by the occurrence of long-chain alkenones, high proportions of lycopane, and high ratios (>0.6) of (n-C35 + lycopane)/n-C31. The occurrence of C37-alkenenones, which were first determined toward the end of the Azolla freshwater event, indicates that the OC becomes more marine in origin during the middle Eocene. Preliminary UK'37- based sea surface temperature (SST) values display a longterm temperature decrease of about 15C during the time interval 49-44.5 Ma (25° to 10°C), coinciding with the global benthic d18O cooling trend after the early Eocene climatic optimum. At about 46 Ma, parallel with onset of ice-rafted debris, SST (interpreted as summer temperatures) decreased to values <15°C. For the late early Miocene a SST of 11°-15°C was determined. Most of the middle Eocene ACEX sediments are characterized by a smooth short-chain n-alkane distribution, which may point to natural oil-type hydrocarbons from leakage of petroleum reservoirs or erosion of related source rocks and redeposition.