Regional ocean model for the Mediterranean, links to model files

A regional ocean general circulation model of the Mediterranean is used to study the climate of the Last Glacial Maximum. The atmospheric forcing for these simulations has been derived from simulations with an atmospheric general circulation model, which in turn was forced with surface conditions from a coarse resolution earth system model. The model is successful in reproducing the general patterns of reconstructed sea surface temperature anomalies with the strongest cooling in summer in the northwestern Mediterranean and weak cooling in the Levantine, although the model underestimates the extent of the summer cooling in the western Mediterranean. However, there is a strong vertical gradient associated with this pattern of summer cooling, which makes the comparison with reconstructions complicated. The exchange with the Atlantic is decreased to roughly one half of its present value, which can be explained by the shallower Strait of Gibraltar as a consequence of lower global sea level. This reduced exchange causes a strong increase of salinity in the Mediterranean in spite of reduced net evaporation.

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.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Three time-slice simulations of the Eemian (eem), the mid-Holocene (hol), and a pre-industrial (pri) control run with the coupled atmosphere-ocean-model ECHAM5/MPI-OM

Orbital forcing does not only exert direct insolation effects, but also alters climate indirectly through feedback mechanisms that modify atmosphere and ocean dynamics and meridional heat and moisture transfers. We investigate the regional effects of these changes by detailed analysis of atmosphere and ocean circulation and heat transports in a coupled atmosphere-ocean-sea ice-biosphere general circulation model (ECHAM5/JSBACH/MPI-OM). We perform long term quasi equilibrium simulations under pre-industrial, mid-Holocene (6000 years before present - yBP), and Eemian (125 000 yBP) orbital boundary conditions. Compared to pre-industrial climate, Eemian and Holocene temperatures show generally warmer conditions at higher and cooler conditions at lower latitudes. Changes in sea-ice cover, ocean heat transports, and atmospheric circulation patterns lead to pronounced regional heterogeneity. Over Europe, the warming is most pronounced over the north-eastern part in accordance with recent reconstructions for the Holocene. We attribute this warming to enhanced ocean circulation in the Nordic Seas and enhanced ocean-atmosphere heat flux over the Barents Shelf in conduction with retreat of sea ice and intensified winter storm tracks over northern Europe.

Age model and stable isotope record of DSDP holes 48-401 and 86-577

Early Paleogene warm climates may have been linked to different modes and sources of deepwater formation. Warm polar temperatures of the Paleocene and Eocene may have resulted from either increased atmospheric trace gases or increased heat transport through deep and intermediate waters. The possibility of increasing ocean heat transport through the production of warm saline deep waters (WSDW) in the Tethyan region has generated considerable interest. In addition, General Circulation Model results indicate that deepwater source regions may be highly sensitive to changing basin configurations. To decipher deepwater changes, we examined detailed benthic foraminiferal faunal and isotopic records of the late Paleocene through the early Eocene (~60 to 50 Ma) from two critical regions: the North Atlantic (Bay of Biscay Site 401) and the Pacific (Shatsky Rise Site 577). These records are compared with published data from the Southern Ocean (Maud Rise Site 690, Islas Orcadas Rise Site 702). During the late Paleocene, similar benthic foraminiferal delta18O values were recorded at all four sites. This indicates uniform deepwater temperatures, consistent with a single source of deep water. The highest delta13C values were recorded in the Southern Ocean and were 0.5 per mil more positive than those of the Pacific. We infer that the Southern Ocean was proximal to a source of nutrient-depleted deep water during the late Paleocene. Upper Paleocene Reflector Ab was cut on the western Bermuda Rise by cyclonically circulating bottom water, also suggesting a vigorous source of bottom water in the Southern Ocean. A dramatic negative excursion in both carbon and oxygen isotopes occurred in the latest Paleocene in the Southern Ocean. This is a short-term (<100 kyr), globally synchronous event which also is apparent in both the Atlantic and Pacific records as a carbon isotopic excursion of approximately 1 per mil. Faunal analyses from the North Atlantic and Pacific sites indicate that the largest benthic foraminiferal faunal turnover of the Cenozoic was synchronous with the isotopic excursion, lending support to the hypothesis that the extinctions were caused by a change in deepwater circulation. We speculate that the Southern Ocean deepwater source was reduced or eliminated at the time of the excursion. During the early Eocene, Southern Ocean delta13C values remained enriched relative to the North Atlantic and Pacific. However, the Southern Ocean was also enriched in delta18O relative to these basins. We interpret that these patterns indicate that although the Southern Ocean was proximal to a source of cool, nutrient-depleted water, the intermediate to upper deep water sites of the North Atlantic and Pacific were ventilated by a different source that probably originated in low latitudes, i.e., WSDW.

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 changes in ocean physics and carbon cycle for LGM experiments with shifting positions of the Southern Hemisphere westerly winds using the MITgcm, links to model results in NetCDF format

We explore the impact of a latitudinal shift in the westerly wind belt over the Southern Ocean on the Atlantic meridional overturning circulation (AMOC) and on the carbon cycle for Last Glacial Maximum background conditions using a state-of-the-art ocean general circulation model. We find that a southward (northward) shift in the westerly winds leads to an intensification (weakening) of no more than 10% of the AMOC. This response of the ocean physics to shifting winds agrees with other studies starting from preindustrial background climate, but the responsible processes are different. In our setup changes in AMOC seemed to be more pulled by upwelling in the south than pushed by downwelling in the north, opposite to what previous studies with different background climate are suggesting. The net effects of the changes in ocean circulation lead to a rise in atmospheric pCO2 of less than 10 atm for both northward and southward shift in the winds. For northward shifted winds the zone of upwelling of carbon- and nutrient-rich waters in the Southern Ocean is expanded, leading to more CO2 outgassing to the atmosphere but also to an enhanced biological pump in the subpolar region. For southward shifted winds the upwelling region contracts around Antarctica, leading to less nutrient export northward and thus a weakening of the biological pump. These model results do not support the idea that shifts in the westerly wind belt play a dominant role in coupling atmospheric CO2 rise and Antarctic temperature during deglaciation suggested by the ice core data.

On the effects of circulation, sediment resuspension and biological incorporation by diatoms in an ocean model of aluminium, link to model data

The distribution of dissolved aluminium in the West Atlantic Ocean shows a mirror image with that of dissolved silicic acid, hinting at intricate interactions between the ocean cycling of Al and Si. The marine biogeochemistry of Al is of interest because of its potential impact on diatom opal remineralisation, hence Si availability. Furthermore, the dissolved Al concentration at the surface ocean has been used as a tracer for dust input, dust being the most important source of the bio-essential trace element iron to the ocean. Previously, the dissolved concentration of Al was simulated reasonably well with only a dust source, and scavenging by adsorption on settling biogenic debris as the only removal process. Here we explore the impacts of (i) a sediment source of Al in the Northern Hemisphere (especially north of ~ 40° N), (ii) the imposed velocity field, and (iii) biological incorporation of Al on the modelled Al distribution in the ocean. The sediment source clearly improves the model results, and using a different velocity field shows the importance of advection on the simulated Al distribution. Biological incorporation appears to be a potentially important removal process. However, conclusive independent data to constrain the Al / Si incorporation ratio by growing diatoms are missing. Therefore, this study does not provide a definitive answer to the question of the relative importance of Al removal by incorporation compared to removal by adsorptive scavenging.

Model simulation output data (CCSM3) related to uplift of Africa and Himalaya/Tibet Plateau and the influence on tropical African vegetation

Hominid evolution in the late Miocene has long been hypothesized to be linked to the retreat of the tropical rainforest in Africa. One cause for the climatic and vegetation change often considered was uplift of Africa, but also uplift of the Himalaya and the Tibetan Plateau was suggested to have impacted rainfall distribution over Africa. Recent proxy data suggest that in East Africa open grassland habitats were available to the common ancestors of hominins and apes long before their divergence and do not find evidence for a closed rainforest in the late Miocene. We used the coupled global general circulation model CCSM3 including an interactively coupled dynamic vegetation module to investigate the impact of topography on African hydro-climate and vegetation. We performed sensitivity experiments altering elevations of the Himalaya and the Tibetan Plateau as well as of East and Southern Africa. The simulations confirm the dominant impact of African topography for climate and vegetation development of the African tropics. Only a weak influence of prescribed Asian uplift on African climate could be detected. The model simulations show that rainforest coverage of Central Africa is strongly determined by the presence of elevated African topography. In East Africa, despite wetter conditions with lowered African topography, the conditions were not favorable enough to maintain a closed rainforest. A discussion of the results with respect to other model studies indicates a minor importance of vegetation-atmosphere or ocean-atmosphere feedbacks and a large dependence of the simulated vegetation response on the land surface/vegetation model.

Glacial-interglacial bioproductivity of the Mexican Margin

Sedimentary accumulation of biogenic components (organic carbon, opal, and biogenic barium) on the northwestern Mexican margin declined during every glacial interval of the past 140 kyr, indicating decreases in upwelling-induced productivity during cold periods. The glacial-interglacial contrasts in upwelling on this margin are attributed to reversals in land-ocean thermal contrast, the waxing and waning of the Laurentide Ice Sheet, and consequent responses of the western hemisphere wind fields. This scenario is consistent with three independent lines of evidence: terrestrial paleoclimatic data, general circulation model results, and our marine records. This pattern of glacial-interglacial variability in upwelling off NW Mexico is opposite to that observed in other low-latitude and midlatitude upwelling areas, such as the eastern equatorial Pacific. These results add to a growing pool of observations that the response of oceanic upwelling to glacial climatic forcing has been regionally variable.

Sea-surface temperature reconstruction of sediments from the North Pacific and North Atlantic

Holocene climate variability is investigated in the North Pacific and North Atlantic realms, using alkenone-derived sea-surface temperature (SST) records as well as a millennial scale simulation with a coupled atmosphere-ocean general circulation model (AOGCM). The alkenone SST data indicate a temperature increase over almost the entire North Pacific from 7 cal kyr BP to the present. A dipole pattern with a continuous cooling in the northeastern Atlantic and a warming in the eastern Mediterranean Sea and the northern Red Sea is detected in the North Atlantic realm. Similarly, SST variations are opposite in sign between the northeastern Pacific and the northeastern Atlantic. A 2300 year long AOGCM climate simulation reveals a similar SST seesaw between the northeastern Pacific and the northeastern Atlantic on centennial time scales. Our analysis of the alkenone SST data and the model results suggests fundamental inter-oceanic teleconnections during the Holocene.

Simulated Methyl iodide emission and concentration using the ocean biogeochemistry model MPIOM/HAMOCC forced by OMIP data

Production pathways of the prominent volatile organic halogen compound methyl iodide (CH3I) are not fully understood. Based on observations, production of CH3I via photochemical degradation of organic material or via phytoplankton production has been proposed. Additional insights could not be gained from correlations between observed biological and environmental variables or from biogeochemical modeling to identify unambiguously the source of methyl iodide. In this study, we aim to address this question of source mechanisms with a three-dimensional global ocean general circulation model including biogeochemistry (MPIOM-HAMOCC (MPIOM - Max Planck Institute Ocean Model HAMOCC - HAMburg Ocean Carbon Cycle model)) by carrying out a series of sensitivity experiments. The simulated fields are compared with a newly available global data set. Simulated distribution patterns and emissions of CH3I differ largely for the two different production pathways. The evaluation of our model results with observations shows that, on the global scale, observed surface concentrations of CH3I can be best explained by the photochemical production pathway. Our results further emphasize that correlations between CH3I and abiotic or biotic factors do not necessarily provide meaningful insights concerning the source of origin. Overall, we find a net global annual CH3I air-sea flux that ranges between 70 and 260 Gg/yr. On the global scale, the ocean acts as a net source of methyl iodide for the atmosphere, though in some regions in boreal winter, fluxes are of the opposite direction (from the atmosphere to the ocean).

Stable oxygen isotope data of Globorotalia inflata and radiocarbon dates for sediment cores GeoB6211-2, GeoB13862-1, and GeoB6308-3

The Southern Westerly Winds (SWW) exert a crucial influence over the world ocean and climate. Nevertheless, a comprehensive understanding of the Holocene temporal and spatial evolution of the SWW remains a significant challenge due to the sparsity of high-resolution marine archives and appropriate SWW proxies. Here, we present a north-south transect of high-resolution planktonic foraminiferal oxygen isotope records from the western South Atlantic. Our proxy records reveal Holocene migrations of the Brazil- Malvinas Confluence (BMC), a highly sensitive feature for changes in the position and strength of the northern portion of the SWW. Through the tight coupling of the BMC position to the large-scale wind field, the records allow a quantitative reconstruction of Holocene latitudinal displacements of the SWW across the South Atlantic. Our data reveal a gradual poleward movement of the SWW by about 1-1.5° from the early to the mid-Holocene. Afterwards variability in the SWW is dominated by millennial-scale displacements in the order of 1° in latitude with no recognizable longer-term trend. These findings are confronted with results from a state-of-the-art transient Holocene climate simulation using a comprehensive coupled atmosphere-ocean general circulation model. Proxy-inferred and modeled SWW shifts compare qualitatively, but the model underestimates both orbitally forced multi-millennial and internal millennial SWW variability by almost an order of magnitude. The underestimated natural variability implies a substantial uncertainty in model projections of future SWW shifts.