Simulated change in atmospheric CO2 and D14C around 14.6 kyr BP, at the onset of the northern hemispheric warming into the Bølling/Allerød

One of the most abrupt and yet unexplained past rises in atmospheric CO2 (10 p.p.m.v. in two centuries) occurred in quasi-synchrony with abrupt northern hemispheric warming into the Bølling/Allerød, 14,600 years ago. Here we use a U/Th-dated record of atmospheric D14C from Tahiti corals to provide an independent and precise age control for this CO2 rise. We also use model simulations to show that the release of old (nearly 14C-free) carbon can explain these changes in CO2 and D14C. The D14C record provides an independent constraint on the amount of carbon released (125 Pg C). We suggest, in line with observations of atmospheric CH4 and terrigenous biomarkers, that thawing permafrost in high northern latitudes could have been the source of carbon, possibly with contribution from flooding of the Siberian continental shelf during meltwater pulse 1A. Our findings highlight the potential of the permafrost carbon reservoir to modulate abrupt climate changes via greenhouse-gas feedbacks.

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.

Simulated climatically disturbed emergence of agricultures in Western Eurasia 8500-3000 BC

Was the spread of agropastoralism from the Fertile Crescent throughout Europe influenced by rapid climatic shifts? We here generate idealized climate events using palaeoclimate records. In a mathematical model of regional sociocultural development, these events disturb the subsistence base of simulated forager and farmer societies. We evaluate the regional simulated transition timings and durations against a published large set of radiocarbon dates for western Eurasia; the model is able to realistically hindcast much of the inhomogeneous space-time evolution of regional Neolithic transitions. Our study shows that the inclusion of climate events improves the simulation of typical lags between cultural complexes, but that the overall difference to a model without climate events is not significant. Climate events may not have been as important for early sociocultural dynamics as endogenous factors.

Simulated transition to agropastoralism in the Indus valley 7500-3000 BC

The Indus Valley Civilization (IVC) was one of the first great civilizations in prehistory. This bronze age civilization flourished from the end of the fourth millennium BC. It disintegrated during the second millennium BC; despite much research effort, this decline is not well understood. Less research has been devoted to the emergence of the IVC, which shows continuous cultural precursors since at least the seventh millennium BC. To understand the decline, we believe it is necessary to investigate the rise of the IVC, i.e., the establishment of agriculture and livestock, dense populations and technological developments 7000-3000 BC. Although much archaeologically typed information is available, our capability to investigate the system is hindered by poorly resolved chronology, and by a lack of field work in the intermediate areas between the Indus valley and Mesopotamia. We thus employ a complementary numerical simulation to develop a consistent picture of technology, agropastoralism and population developments in the IVC domain. Results from this Global Land Use and technological Evolution Simulator show that there is (1) fair agreement between the simulated timing of the agricultural transition and radiocarbon dates from early agricultural sites, but the transition is simulated first in India then Pakistan; (2) an independent agropas- toralism developing on the Indian subcontinent; and (3) a positive relationship between archeological artifact richness and simulated population density which remains to be quantified.

Simulated agricultural subsistence, timing of transition, and percentage of immigrant farmers

Farming and herding were introduced to Europe from the Near East and Anatolia; there are, however, considerable arguments about the mechanisms of this transition. Were it the people who moved and either outplaced, or admixed with, the indigenous hunter-gatherer groups? Or was it material and information that moved---the Neolithic Package---consisting of domesticated plants and animals and the knowledge of their use? The latter process is commonly referred to as cultural diffusion and the former as demic diffusion. Despite continuous and partly combined efforts by archaeologists, anthropologists, linguists, palaeontologists and geneticists, a final resolution of the debate has not yet been reached. In the present contribution we interpret results from the Global Land Use and technological Evolution Simulator (GLUES). GLUES is a mathematical model for regional sociocultural development, embedded in the geoenvironmental context, during the Holocene. We demonstrate that the model is able to realistically hindcast the expansion speed and the inhomogeneous space-time evolution of the transition to agropastoralism in western Eurasia. In contrast to models that do not resolve endogenous sociocultural dynamics, our model describes and explains how and why the Neolithic advanced in stages. We uncouple the mechanisms of migration and information exchange and also of migration and the spread of agropastoralism. We find that: (1) An indigenous form of agropastoralism could well have arisen in certain Mediterranean landscapes, but not in Northern and Central Europe, where it depended on imported technology and material. (2) Both demic diffusion by migration and cultural diffusion by trade may explain the western European transition equally well. (3) Migrating farmers apparently contribute less than local adopters to the establishment of agropastoralism. Our study thus underlines the importance of adoption of introduced technologies and economies by resident foragers.

Simulated Bromoform emission and concentration using the ocean biogeochemistry model MPIOM/HAMOCC forced by 6-hourly NCEP data

Bromoform (CHBr3) is one important precursor of atmospheric reactive bromine species that are involved in ozone depletion in the troposphere and stratosphere. In the open ocean bromoform production is linked to phytoplankton that contains the enzyme bromoperoxidase. Coastal sources of bromoform are higher than open ocean sources. However, open ocean emissions are important because the transfer of tracers into higher altitude in the air, i.e. into the ozone layer, strongly depends on the location of emissions. For example, emissions in the tropics are more rapidly transported into the upper atmosphere than emissions from higher latitudes. Global spatio-temporal features of bromoform emissions are poorly constrained. Here, a global three-dimensional ocean biogeochemistry model (MPIOM-HAMOCC) is used to simulate bromoform cycling in the ocean and emissions into the atmosphere using recently published data of global atmospheric concentrations (Ziska et al., 2013) as upper boundary conditions. Our simulated surface concentrations of CHBr3 match the observations well. Simulated global annual emissions based on monthly mean model output are lower than previous estimates, including the estimate by Ziska et al. (2013), because the gas exchange reverses when less bromoform is produced in non-blooming seasons. This is the case for higher latitudes, i.e. the polar regions and northern North Atlantic. Further model experiments show that future model studies may need to distinguish different bromoform-producing phytoplankton species and reveal that the transport of CHBr3 from the coast considerably alters open ocean bromoform concentrations, in particular in the northern sub-polar and polar regions.

Experiment: Organic matter exudation by Emiliania huxleyi under simulated future ocean conditions

Emiliania huxleyi (strain B 92/11) was exposed to different nutrient supply, CO2 and temperature conditions in phosphorus controlled chemostats to investigate effects on organic carbon exudation and partitioning between the pools of particulate organic carbon (POC) and dissolved organic carbon (DOC). 14C incubation measurements for primary production (PP) and extracellular release (ER) were performed. Chemical analysis included the amount and composition of high molecular weight (>1 kDa) dissolved combined carbohydrates (HMW-dCCHO), particulate combined carbohydrates (pCCHO) and the carbon content of transparent exopolymer particles (TEP-C). Applied CO2 and temperature conditions were 300, 550 and 900 µatm pCO2 at 14 °C, and additionally 900 µatm pCO2 at 18 °C simulating a greenhouse ocean scenario. Enhanced nutrient stress by reducing the dilution rate (D) from D = 0.3 /d to D = 0.1 /d (D = µ) induced the strongest response in E. huxleyi. At µ = 0.3 /d, PP was significantly higher at elevated CO2 and temperature and DO14C production correlated to PO14C production in all treatments, resulting in similar percentages of extracellular release (PER; (DO14C production/PP) × 100) averaging 3.74 ± 0.94%. At µ = 0.1 /d, PO14C production decreased significantly, while exudation of DO14C increased. Thus, indicating a stronger partitioning from the particulate to the dissolved pool. Maximum PER of 16.3 ± 2.3% were observed at µ = 0.1 /d at elevated CO2 and temperature. While cell densities remained constant within each treatment and throughout the experiment, concentrations of HMW-dCCHO, pCCHO and TEP were generally higher under enhanced nutrient stress. At µ= 0.3 /d, pCCHO concentration increased significantly with elevated CO2 and temperature. At µ = 0.1 /d, the contribution (mol % C) of HMW-dCCHO to DOC was lower at elevated CO2 and temperature while pCCHO and TEP concentrations were higher. This was most pronounced under greenhouse conditions. Our findings suggest a stronger transformation of primary produced DOC into POC by coagulation of exudates under nutrient limitation. Our results further imply that elevated CO2 and temperature will increase exudation by E. huxleyi and may affect organic carbon partitioning in the ocean due to an enhanced transfer of HMW-dCCHO to TEP by aggregation processes.

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.

Model results in NetCDF format of North African vegetation-precipitation in early and mid-Holocene climate simulated by CCSM3-DGVM//

The present study analyses the sign, strength, and working mechanism of the vegetation-precipitation feedback over North Africa in middle (6 ka BP) and early Holocene (9 ka BP) simulations using the comprehensive coupled climate-vegetation model CCSM3-DGVM (Community Climate System Model version 3 and a dynamic global vegetation model). The coupled model simulates enhanced summer rainfall and a northward migration of the West African monsoon trough along with an expansion of the vegetation cover for the early and middle Holocene compared to the pre-industrial period. It is shown that dynamic vegetation enhances the orbitally triggered summer precipitation anomaly by approximately 20% in the Sahara-Sahel region (10-25° N, 20° W-30° E) in both the early and mid-Holocene experiments compared to their fixed-vegetation counterparts. The primary vegetation-rainfall feedback identified here operates through surface latent heat flux anomalies by canopy evaporation and transpiration and their effect on the mid-tropospheric African easterly jet, whereas the effects of vegetation changes on surface albedo and local water recycling play a negligible role. Even though CCSM3-DGVM simulates a positive vegetation-precipitation feedback in the North African region, this feedback is not strong enough to produce multiple equilibrium climate-ecosystem states on a regional scale.

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).

Effect of simulated ocean acidification on the acute toxicity of Cu and Cd to Tigriopus japonicus

Heavy metals pollution in marine environments has caused great damage to marine biological and ecological systems. Heavy metals accumulate in marine creatures, after which they are delivered to higher trophic levels of marine organisms through the marine food chain, which causes serious harm to marine biological systems and human health. Additionally, excess carbon dioxide in the atmosphere has caused ocean acidification. Indeed, about one third of the CO2 released into the atmosphere by anthropogenic activities since the beginning of the industrial revolution has been absorbed by the world's oceans, which play a key role in moderating climate change. Modeling has shown that, if current trends in CO2 emissions continue, the average pH of the ocean will reach 7.8 by the end of this century, corresponding to 0.5 units below the pre-industrial level, or a three-fold increase in H+ concentration. The ocean pH has not been at this level for several millions of years. Additionally, these changes are occurring at speeds 100 times greater than ever previously observed. As a result, several marine species, communities and ecosystems might not have time to acclimate or adapt to these fast changes in ocean chemistry. In addition, decreasing ocean pH has the potential to seriously affect the growth, development and reproduction reproductive processes of marine organisms, as well as threaten normal development of the marine ecosystem. Copepods are an important part of the meiofauna that play an important role in the marine ecosystem. Pollution of the marine environment can influence their growth and development, as well as the ecological processes they are involved in. Accordingly, there is important scientific value to investigation of the response of copepods to ocean acidification and heavy metals pollution. In the present study, we evaluated the effects of simulated future ocean acidification and the toxicological interaction between ocean acidity and heavy metals of Cu and Cd on T. japonicus. To accomplish this, harpacticoids were exposed to Cu and Cd concentration gradient seawater that had been equilibrated with CO2 and air to reach pH 8.0, 7.7, 7.3 and 6.5 for 96 h. Survival was not significantly suppressed under single sea water acidification, and the final survival rates were greater than 93% in both the experimental groups and the controls. The toxicity of Cu to T. japonicus was significantly affected by sea water acidification, with the 96h LC50 decreasing by nearly threefold from 1.98 to 0.64 mg/L with decreasing pH. The 96 h LC50 of Cd decreased with decreasing pH, but there was no significant difference in mortality among pH treatments. The results of the present study demonstrated that the predicted future ocean acidification has the potential to negatively affect survival of T. japonicus by exacerbating the toxicity of Cu. The calculated safe concentrations of Cu were 11.9 (pH 7.7) and 10.5 (pH 7.3) µg/L, which were below the class I value and very close to the class II level of the China National Quality Standard for Sea Water. Overall, these results indicate that the Chinese coastal sea will face a