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.

Simulated temperature changes as a result of Antarctic ice sheet growth and CO2 reduction

During the Middle Miocene climate transition about 14 million years ago, the Antarctic ice sheet expanded to near-modern volume. Surprisingly, this ice sheet growth was accompanied by a warming in the surface waters of the Southern Ocean, whereas a slight deep-water temperature increase was delayed by more than 200 thousand years. Here we use a coupled atmosphere-ocean model to assess the relative effects of changes in atmospheric CO2 concentration and ice sheet growth on regional and global temperatures. In the simulations, changes in the wind field associated with the growth of the ice sheet induce changes in ocean circulation, deep-water formation and sea-ice cover that result in sea surface warming and deep-water cooling in large swaths of the Atlantic and Indian ocean sectors of the Southern Ocean. We interpret these changes as the dominant ocean surface response to a 100-thousand-year phase of massive ice growth in Antarctica. A rise in global annual mean temperatures is also seen in response to increased Antarctic ice surface elevation. In contrast, the longer-term surface and deep-water temperature trends are dominated by changes in atmospheric CO2 concentration. We therefore conclude that the climatic and oceanographic impacts of the Miocene expansion of the Antarctic ice sheet are governed by a complex interplay between wind field, ocean circulation and the sea-ice system.

Density of alkyl esters and its mixtures

Biodiesel density is a key parameter in biodiesel simulations and process development. In this work we selected, evaluated and improved two density models, one theoretical (Rackett-Soave) and one empirical (Lapuerta's method) for methanol based biodiesels (FAME) and ethanol based biodiesel (FAEE). For this purpose, biodiesel was produced from vegetable oils (sunflower, rapeseed, soybean, olive, safflower and other two commercial mixtures of vegetable oils) and animal fats (edible and crude pork fat and beef tallow) using both methanol and ethanol for the transesterification reactions, and blended to get 21 FAME and 21 FAEE, reporting their density and detailed composition. Bibliographic data have also been used. The Rackett-Soave method has been improved by the use of a new acentric factor correlation, whereas the parameters of the empirical one are improved by considering a bigger density data bank. Results show that the evaluated models could be used to estimate the biodiesel density with a good grade of accuracy but the performed modifications improve the accuracy of the models: ARD (%) for FAME; 0.33, and FAEE; 0.26, both calculated with the modification of Rackett-Soave method and ARD (%) for FAME; 0.40 calculated with the modification of the Lapuerta's method).

Transient Holocene experiments under orbital forcing using the comprehensive global climate model CCSM3 (Community Climate System Model 3)

The Southern Hemisphere Westerly Winds (SWW) have been suggested to exert a critical influence on global climate through wind-driven upwelling of deep water in the Southern Ocean and the potentially resulting atmospheric CO2 variations. The investigation of the temporal and spatial evolution of the SWW along with forcings and feedbacks remains a significant challenge in climate research. In this study, the evolution of the SWW under orbital forcing from the early Holocene (9 kyr BP) to pre-industrial modern times is examined with transient experiments using the comprehensive coupled global climate model CCSM3. Analyses of the model results suggest that the annual and seasonal mean SWW were subject to an overall strengthening and poleward shifting trend during the course of the early-to-late Holocene under the influence of orbital forcing, except for the austral spring season, where the SWW exhibited an opposite trend of shifting towards the equator.

Sedimentary record of the western Amundsen Sea Embayment

The Amundsen Sea Embayment (ASE) drains approximately 35% of the West Antarctic Ice Sheet (WAIS) and is one of the most rapidly changing parts of the cryosphere. In order to predict future ice-sheet behaviour, modellers require long-term records of ice-sheet melting to constrain and build confidence in their simulations. Here, we present detailed marine geological and radiocarbon data along three palaeo-ice stream tributary troughs in the western ASE to establish vital information on the timing of deglaciation of the WAIS since the Last Glacial Maximum (LGM). We have undertaken multi-proxy analyses of the cores (core description, shear strength, x-radiographs, magnetic susceptibility, wet bulk density, total organic carbon/nitrogen, carbonate content and clay mineral analyses) in order to: (1) characterise the sedimentological facies and depositional environments; and (2) identify the horizon(s) in each core that would yield the most reliable age for deglaciation. In accordance with previous studies we identify three key facies, which offer the most reliable stratigraphies for dating deglaciation by recording the transition from a grounded ice sheet to open marine environments. These facies are: i) subglacial, ii) proximal grounding-line, and iii) seasonal open-marine. In addition, we incorporate ages from other facies (e.g., glaciomarine diamictons deposited at some distance from the grounding line, such as glaciogenic debris flows and iceberg rafted diamictons and turbates) into our deglacial model. In total, we have dated 78 samples (mainly the acid insoluble organic (AIO) fraction, but also calcareous foraminifers), which include 63 downcore and 15 surface samples. Through careful sample selection prior to dating, we have established a robust deglacial chronology for this sector of the WAIS. Our data show that deglaciation of the western ASE was probably underway as early as 22,351 calibrated years before present (cal 44 yr BP), reaching the mid-shelf by 13,837 cal yr BP and the inner shelf to within c.10-12 km of the present ice shelf front between 12,618 and 10,072 cal yr BP. The deglacial steps in the western ASE broadly coincide with the rapid rises in sea-level associated with global meltwater pulses 1a and 1b, although given the potential dating uncertainty, additional, more precise ages are required before these findings can be fully substantiated. Finally, we show that the rate of ice-sheet retreat increased across the deep (up to1,600 m) basins of the inner shelf, highlighting the importance of reverse slope and pinning points in accelerated phases of deglaciation.

Theoretical model simulations for variations in thermal reaction norms

Thermal reaction norms for growth rates of six Emiliania huxleyi isolates originating from the central Atlantic (Azores, Portugal) and five isolates from the coastal North Atlantic (Bergen, Norway) were assessed. We used the template mode of variation model to decompose variations in growth rates into modes of biological interest: vertical shift, horizontal shift, and generalist-specialist variation. In line with the actual habitat conditions, isolates from Bergen (Bergen population) grew well at lower temperatures, and isolates from the Azores (Azores population) performed better at higher temperatures. The optimum growth temperature of the Azores population was significantly higher than that of the Bergen population. Neutral genetic differentiation was found between populations by microsatellite analysis. These findings indicate that E. huxleyi populations are adapted to local temperature regimes. Next to between-population variation, we also found variation within populations. Genotype-by-environment interactions resulted in the most pronounced phenotypic differences when isolates were exposed to temperatures outside the range they naturally encounter. Variation in thermal reaction norms between and within populations emphasizes the importance of using more than one isolate when studying the consequences of global change on marine phytoplankton. Phenotypic plasticity and standing genetic variation will be important in determining the potential of natural E. huxleyi populations to cope with global climate change.

Bacteria growth, enzyme activities and TEP during KOSMOS 2011 mesocosm experiment in the Raunefjord

Marine bacteria are the main consumers of freshly produced organic matter. Many enzymatic processes involved in the bacterial digestion of organic compounds were shown to be pH sensitive in previous studies. Due to the continuous rise in atmospheric CO2 concentration, seawater pH is presently decreasing at a rate unprecedented during the last 300 million years but the consequences for microbial physiology, organic matter cycling and marine biogeochemistry are still unresolved. We studied the effects of elevated seawater pCO2 on a natural plankton community during a large-scale mesocosm study in a Norwegian fjord. Nine Kiel Off-Shore Mesocosms for Future Ocean Simulations (KOSMOS) were adjusted to different pCO2 levels ranging initially from ca. 280 to 3000 µatm and sampled every second day for 34 days. The first phytoplankton bloom developed around day 5. On day 14, inorganic nutrients were added to the enclosed, nutrient-poor waters to stimulate a second phytoplankton bloom, which occurred around day 20. Our results indicate that marine bacteria benefit directly and indirectly from decreasing seawater pH. During the first phytoplankton bloom, 5-10% more transparent exopolymer particles were formed in the high pCO2 mesocosms. Simultaneously, the efficiency of the protein-degrading enzyme leucine aminopeptidase increased with decreasing pH resulting in up to three times higher values in the highest pCO2/lowest pH mesocosm compared to the controls. In general, total and cell-specific aminopeptidase activities were elevated under low pH conditions. The combination of enhanced enzymatic hydrolysis of organic matter and increased availability of gel particles as substrate supported up to 28% higher bacterial abundance in the high pCO2 treatments. We conclude that ocean acidification has the potential to stimulate the bacterial community and facilitate the microbial recycling of freshly produced organic matter, thus strengthening the role of the microbial loop in the surface ocean.

Transient simulations for present and last interglacials using a comprehensive coupled climate model CCSM3 (Community Climate System Model 3.0), links to model result files in NetCDF format

Transient simulations are widely used in studying the past climate as they provide better comparison with any exisiting proxy data. However, multi-millennial transient simulations using coupled climate models are usually computationally very expensive. As a result several acceleration techniques are implemented when using numerical simulations to recreate past climate. In this study, we compare the results from transient simulations of the present and the last interglacial with and without acceleration of the orbital forcing, using the comprehensive coupled climate model CCSM3 (Community Climate System Model 3). Our study shows that in low-latitude regions, the simulation of long-term variations in interglacial surface climate is not significantly affected by the use of the acceleration technique (with an acceleration factor of 10) and hence, large-scale model-data comparison of surface variables is not hampered. However, in high-latitude regions where the surface climate has a direct connection to the deep ocean, e.g. in the Southern Ocean or the Nordic Seas, acceleration-induced biases in sea-surface temperature evolution may occur with potential influence on the dynamics of the overlying atmosphere. The data provided here are from both accelerated and non-accelerated runs as decadal mean values.