Bio-physical interactions and feedbacks in a global climate model

This PhD thesis addresses the topic of large-scale interactions between climate and marine biogeochemistry. To this end, centennial simulations are performed under present and projected future climate conditions with a coupled ocean-atmosphere model containing a complex marine biogeochemistry model. The role of marine biogeochemistry in the climate system is first investigated. Phytoplankton solar radiation absorption in the upper ocean enhances sea surface temperatures and upper ocean stratification. The associated increase in ocean latent heat losses raises atmospheric temperatures and water vapor. Atmospheric circulation is modified at tropical and extratropical latitudes with impacts on precipitation, incoming solar radiation, and ocean circulation which cause upper-ocean heat content to decrease at tropical latitudes and to increase at middle latitudes. Marine biogeochemistry is tightly related to physical climate variability, which may vary in response to internal natural dynamics or to external forcing such as anthropogenic carbon emissions. Wind changes associated with the North Atlantic Oscillation (NAO), the dominant mode of climate variability in the North Atlantic, affect ocean properties by means of momentum, heat, and freshwater fluxes. Changes in upper ocean temperature and mixing impact the spatial structure and seasonality of North Atlantic phytoplankton through light and nutrient limitations. These changes affect the capability of the North Atlantic Ocean of absorbing atmospheric CO2 and of fixing it inside sinking particulate organic matter. Low-frequency NAO phases determine a delayed response of ocean circulation, temperature and salinity, which in turn affects stratification and marine biogeochemistry. In 20th and 21st century simulations natural wind fluctuations in the North Pacific, related to the two dominant modes of atmospheric variability, affect the spatial structure and the magnitude of the phytoplankton spring bloom through changes in upper-ocean temperature and mixing. The impacts of human-induced emissions in the 21st century are generally larger than natural climate fluctuations, with the phytoplankton spring bloom starting one month earlier than in the 20th century and with ~50% lower magnitude. This PhD thesis advances the knowledge of bio-physical interactions within the global climate, highlighting the intrinsic coupling between physical climate and biosphere, and providing a framework on which future studies of Earth System change can be built on.

Simulating short lived carbonaceous compounds with an atmospheric chemistry general circulation model

Ein neu entwickeltes globales Atmosphärenchemie- und Zirkulationsmodell (ECHAM5/MESSy1) wurde verwendet um die Chemie und den Transport von Ozonvorläufersubstanzen zu untersuchen, mit dem Schwerpunkt auf Nichtmethankohlenwasserstoffen. Zu diesem Zweck wurde das Modell durch den Vergleich der Ergebnisse mit Messungen verschiedenen Ursprungs umfangreich evaluiert. Die Analyse zeigt, daß das Modell die Verteilung von Ozon realistisch vorhersagt, und zwar sowohl die Menge als auch den Jahresgang. An der Tropopause gibt das Modell den Austausch zwischen Stratosphäre und Troposphäre ohne vorgeschriebene Flüsse oder Konzentrationen richtig wieder. Das Modell simuliert die Ozonvorläufersubstanzen mit verschiedener Qualität im Vergleich zu den Messungen. Obwohl die Alkane vom Modell gut wiedergeben werden, ergibt sich einige Abweichungen für die Alkene. Von den oxidierten Substanzen wird Formaldehyd (HCHO) richtig wiedergegeben, während die Korrelationen zwischen Beobachtungen und Modellergebnissen für Methanol (CH3OH) und Aceton (CH3COCH3) weitaus schlechter ausfallen. Um die Qualität des Modells im Bezug auf oxidierte Substanzen zu verbessern, wurden einige Sensitivitätsstudien durchgeführt. Diese Substanzen werden durch Emissionen/Deposition von/in den Ozean beeinflußt, und die Kenntnis über den Gasaustausch mit dem Ozean ist mit großen Unsicherheiten behaftet. Um die Ergebnisse des Modells ECHAM5/MESSy1 zu verbessern wurde das neue Submodell AIRSEA entwickelt und in die MESSy-Struktur integriert. Dieses Submodell berücksichtigt den Gasaustausch zwischen Ozean und Atmosphäre einschließlich der oxidierten Substanzen. AIRSEA, welches Informationen über die Flüssigphasenkonzentration des Gases im Oberflächenwasser des Ozeans benötigt wurde ausgiebig getestet. Die Anwendung des neuen Submodells verbessert geringfügig die Modellergebnisse für Aceton und Methanol, obwohl die Verwendung einer vorgeschriebenen Flüssigphasenkonzentration stark den Erfolg der Methode einschränkt, da Meßergebnisse nicht in ausreichendem Maße zu Verfügung stehen. Diese Arbeit vermittelt neue Einsichten über organische Substanzen. Sie stellt die Wichtigkeit der Kopplung zwischen Ozean und Atmosphäre für die Budgets vieler Gase heraus.

Study of the inter-annual variability of particle vertical fluxes in two moorings in the Ross Sea (Antarctica)

The most ocean - atmosphere exchanges take place in polar environments due to the low temperatures which favor the absorption processes of atmospheric gases, in particular CO2. For this reason, the alterations of biogeochemical cycles in these areas can have a strong impact on the global climate. With the aim of contributing to the definition of the mechanisms regulating the biogeochemical fluxes we have analyzed the particles collected in the Ross Sea in different years (ROSSMIZE, BIOSESO 1 and 2, ROAVERRS and ABIOCLEAR projects) in two sites (mooring A and B). So it has been developed a more efficient method to prepare sediment trap samples for the analyses. We have also processed satellite data of sea ice, chlorophyll a and diatoms concentration. At both sites, in each year considered, there was a high seasonal and inter-annual variability of biogeochemical fluxes closely correlated with sea ice cover and primary productivity. The comparison between the samples collected at mooring A and B in 2008 highlighted the main differences between these two sites. Particle fluxes at Mooring A, located in a polynia area, are higher than mooring B ones and they happen about a month before. In the mooring B area it has been possible to correlate the particles fluxes to the ice concentration anomalies and with the atmospheric changes in response to El Niño Southern Oscillations. In 1996 and 1999, years subjected to La Niña, the concentrations of sea ice in this area have been less than in 1998, year subjected to El Niño. Inverse correlation was found for 2005 and 2008. In the mooring A area significant differences in mass and biogenic fluxes during 2005 and 2008 has been recorded. This allowed to underline the high variability of lateral advection processes and to connect them to the physical forcing.

Causes and consequences of past and projected Scandinavian summer temperatures, 500-2100 AD

Tree rings dominate millennium-long temperature reconstructions and many records originate from Scandinavia, an area for which the relative roles of external forcing and internal variation on climatic changes are, however, not yet fully understood. Here we compile 1,179 series of maximum latewood density measurements from 25 conifer sites in northern Scandinavia, establish a suite of 36 subset chronologies, and analyse their climate signal. A new reconstruction for the 1483–2006 period correlates at 0.80 with June–August temperatures back to 1860. Summer cooling during the early 17th century and peak warming in the 1930s translate into a decadal amplitude of 2.9°C, which agrees with existing Scandinavian tree-ring proxies. Climate model simulations reveal similar amounts of mid to low frequency variability, suggesting that internal ocean-atmosphere feedbacks likely influenced Scandinavian temperatures more than external forcing. Projected 21st century warming under the SRES A2 scenario would, however, exceed the reconstructed temperature envelope of the past 1,500 years.

A reconstruction of atmospheric carbon dioxide and its stable carbon isotopic composition from the penultimate glacial maximum to the last glacial inception

The reconstruction of the stable carbon isotope evolution in atmospheric CO2 (δ13Catm), as archived in Antarctic ice cores, bears the potential to disentangle the contributions of the different carbon cycle fluxes causing past CO2 variations. Here we present a new record of δ13Catm before, during and after the Marine Isotope Stage 5.5 (155 000 to 105 000 yr BP). The dataset is archived on the data repository PANGEA® (www.pangea.de) under 10.1594/PANGAEA.817041. The record was derived with a well established sublimation method using ice from the EPICA Dome C (EDC) and the Talos Dome ice cores in East Antarctica. We find a 0.4‰ shift to heavier values between the mean δ13Catm level in the Penultimate (~ 140 000 yr BP) and Last Glacial Maximum (~ 22 000 yr BP), which can be explained by either (i) changes in the isotopic composition or (ii) intensity of the carbon input fluxes to the combined ocean/atmosphere carbon reservoir or (iii) by long-term peat buildup. Our isotopic data suggest that the carbon cycle evolution along Termination II and the subsequent interglacial was controlled by essentially the same processes as during the last 24 000 yr, but with different phasing and magnitudes. Furthermore, a 5000 yr lag in the CO2 decline relative to EDC temperatures is confirmed during the glacial inception at the end of MIS5.5 (120 000 yr BP). Based on our isotopic data this lag can be explained by terrestrial carbon release and carbonate compensation.

A Model Study of the Seasonal Circulation in the Gulf of Maine

The Princeton Ocean Model is used to study the circulation in the Gulf of Maine and its seasonal transition in response to wind, surface heat flux, river discharge, and the M-2 tide. The model has an orthogonal-curvature linear grid in the horizontal with variable spacing from 3 km nearshore to 7 km offshore and 19 levels in the vertical. It is initialized and forced at the open boundary with model results from the East Coast Forecast System. The first experiment is forced by monthly climatological wind and heat flux from the Comprehensive Ocean Atmosphere Data Set; discharges from the Saint John, Penobscot, Kennebec, and Merrimack Rivers are added in the second experiment; the semidiurnal lunar tide (M-2) is included as part of the open boundary forcing in the third experiment. It is found that the surface heat flux plays an important role in regulating the annual cycle of the circulation in the Gulf of Maine. The spinup of the cyclonic circulation between April and June is likely caused by the differential heating between the interior gulf and the exterior shelf/slope region. From June to December the cyclonic circulation continues to strengthen, but gradually shrinks in size. When winter cooling erodes the stratification, the cyclonic circulation penetrates deeper into the water column. The circulation quickly spins down from December to February as most of the energy is consumed by bottom friction. While inclusion of river discharge changes details of the circulation pattern, the annual evolution of the circulation is largely unaffected. On the other hand, inclusion of the tide results in not only the anticyclonic circulation on Georges Bank but also modifications to the seasonal circulation.

Freshening of the Alaska Coastal Current Recorded By Coralline Algal Ba/Ca Ratios

Arctic Ocean freshening can exert a controlling influence on global climate, triggering strong feedbacks on ocean-atmospheric processes and affecting the global cycling of the world's oceans. Glacier-fed ocean currents such as the Alaska Coastal Current are important sources of freshwater for the Bering Sea shelf, and may also influence the Arctic Ocean freshwater budget. Instrumental data indicate a multiyear freshening episode of the Alaska Coastal Current in the early 21st century. It is uncertain whether this freshening is part of natural multidecadal climate variability or a unique feature of anthropogenically induced warming. In order to answer this, a better understanding of past variations in the Alaska Coastal Current is needed. However, continuous long-term high-resolution observations of the Alaska Coastal Current have only been available for the last 2 decades. In this study, specimens of the long-lived crustose coralline alga Clathromorphum nereostratum were collected within the pathway of the Alaska Coastal Current and utilized as archives of past temperature and salinity. Results indicate that coralline algal Mg/Ca ratios provide a 60 year record of sea surface temperatures and track changes of the Pacific Decadal Oscillation, a pattern of decadal-to-multidecadal ocean-atmosphere climate variability centered over the North Pacific. Algal Ba/Ca ratios (used as indicators of coastal freshwater runoff) are inversely correlated to instrumentally measured Alaska Coastal Current salinity and record the period of freshening from 2001 to 2006. Similar multiyear freshening events are not evident in the earlier portion of the 60 year Ba/Ca record. This suggests that the 21st century freshening of the Alaska Coastal Current is a unique feature related to increasing glacial melt and precipitation on mainland Alaska.

Sr/Ca ratios of coccoliths from the Paleocene-Eocene Thermal Maximum

During the Paleocene-Eocene Thermal Maximum (PETM), rapid release of isotopically light C to the ocean-atmosphere system elevated the greenhouse effect and warmed temperatures by 5-7 °C for 105 yr. The response of the planktic ecosystems and productivity to the dramatic climate changes of the PETM may represent a significant feedback to the carbon cycle changes, but has been difficult to document. We examine Sr/Ca ratios in calcareous nannofossils in sediments spanning the PETM in three open ocean sites as a new approach to examine productivity and ecological shifts in calcifying plankton. The large heterogeneity in Sr/Ca among different nannofossil genera indicates that nannofossil Sr/Ca reflects primary productivity-driven geochemical signals and not diagenetic overprinting. Elevated Sr/Ca ratios in several genera and constant ratios in other genera suggest increased overall productivity in the Atlantic sector of the Southern Ocean during the PETM. Dominant nannofossil genera in tropical Atlantic and Pacific sites show Sr/Ca variations during the PETM which are comparable to background variability prior to the PETM. Despite acidification of the ocean there was not a productivity crisis among calcifying phytoplankton. We use the Pandora ocean box model to explore possible mechanisms for PETM productivity change. If independent proxy evidence for more stratified conditions in the Southern Ocean during the PETM is robust, then maintenance of stable or increased productivity there likely reflects increased nutrient inventories of the ocean. Increased nutrient inventories could have resulted from climatically enhanced weathering and would have important implications for burial rates of organic carbon and stabilization of climate and the carbon cycle.

Modern and fossil corals from Palmyra and Christmas islands

The relationship between decadal to centennial changes in ocean circulation and climate is difficult to discern using the sparse and discontinuous instrumental record of climate and, as such, represents a large uncertainty in coupled ocean-atmosphere general circulation models. We present new modern and fossil coral radiocarbon (D14C) records from Palmyra (6°N, 162°W) and Christmas (2°N, 157°W) islands to constrain central tropical Pacific ocean circulation changes during the last millennium. Seasonally to annually resolved coral D14C measurements from the 10th, 12th-17th, and 20th centuries do not contain significant interannual to decadal-scale variations, despite large changes in coral d18O on these timescales. A centennial-scale increase in coral radiocarbon from the Medieval Climate Anomaly (~900-1200 AD) to the Little Ice Age (~1500-1800) can be largely explained by changes in the atmospheric D14C, as determined with a box model of Palmyra mixed layer D14C. However, large 12th century depletions in Palmyra coral D14C may reflect as much as a 100% increase in upwelling rates and/or a significant decrease in the D14C of higher-latitude source waters reaching the equatorial Pacific during this time. SEM photos reveal evidence for minor dissolution and addition of secondary aragonite in the fossil corals, but our results suggest that coral D14C is only compromised after moderate to severe diagenesis for these relatively young fossil corals.

Paleoceanographic data from cores M23414 and PS1243

Variations in the poleward-directed Atlantic heat transfer was investigated over the past 135 ka with special emphasis on the last and present interglacial climate development (Eemian and Holocene). Both interglacials exhibited very similar climatic oscillations during each preceding glacial terminations (deglacial TI and TII). Like TI, also TII has pronounced cold-warm-cold changes akin to events such as H1, Bølling/Allerød, and the Younger Dryas. But unlike TI, the cold events in TII were associated with intermittent southerly invasions of an Atlantic faunal component which underscores quite a different water mass evolution in the Nordic Seas. Within the Eemian interglaciation proper, peak warming intervals were antiphased between the Nordic Seas and North Atlantic. Moreover, inferred temperatures for the Nordic Seas were generally colder in the Eemian than in the Holocene, and vice versa for the North Atlantic. A reduced intensity of Atlantic Ocean heat transfer to the Arctic therefore characterized the Eemian, requiring a reassessment of the actual role of the ocean-atmosphere system behind interglacial, but also, glacial climate changes.

Deuterium analysis of ice core GISP2

The Greenland Ice Sheet Project 2 (GISP2) core can enhance our understanding of the relationship between parameters measured in the ice in central Greenland and variability in the ocean, atmosphere, and cryosphere of the North Atlantic Ocean and adjacent land masses. Seasonal (summer, winter) to annual responses of dD and deuterium excess isotopic signals in the GISP2 core to the seesaw in winter temperatures between West Greenland and northern Europe from A.D. 1840 to 1970 are investigated. This seesaw represents extreme modes of the North Atlantic Oscillation, which also influences sea surface temperatures (SSTs), atmospheric pressures, geostrophic wind strength, and sea ice extents beyond the winter season. Temperature excursions inferred from the dD record during seesaw/extreme NAO mode years move in the same direction as the West Greenland side of the seesaw. Symmetry with the West Greenland side of the seesaw suggests a possible mechanism for damping in the ice core record of the lowest decadal temperatures experienced in Europe from A.D. 1500 to 1700. Seasonal and annual deuterium excess excursions during seesaw years show negative correlation with dD. This suggests an isotopic response to a SST/ land temperature seesaw. The isotopic record from GISP2 may therefore give information on both ice sheet and sea surface temperature variability. Cross-plots of dD and d show a tendency for data to be grouped according to the prevailing mode of the seesaw, but do not provide unambiguous identification of individual seesaw years. A combination of ice core and tree ring data sets may allow more confident identification of GA and GB (extreme NAO mode) years prior to 1840.

Age models and stable isotope ratios from the Mid-Miocene

Mid-Miocene pelagic sedimentary sections can be correlated using intermediate and high resolution oxygen and carbon isotopic records of benthic foraminifera. Precision of a few tens of thousands of years is readily achievable at sites with high sedimentation rates, for example, Deep Sea Drilling Project sites 289 and 574. The mid-Miocene carbon isotope records are characterized by an interval of high d13C values between 17 and 13.5 Ma (the Monterey Excursion of Vincent and Berger 1985) upon which are superimposed a series of periodic or quasi-periodic fluctuations in d13C values. These fluctuations have a period of approximately 440 kyr, suggestive of the 413 kyr cycle predicted by Milankovitch theory. Vincent and Berger proposed that the Monterey Excursion was the result of increased organic carbon burial in continental margins sediments. The increased d13C values (called 13C maxima) superimposed on the generally high mid-Miocene signal coincide with increases in d18O values suggesting that periods of cooling and/or ice buildup were associated with exceptionally rapid burial of organic carbon and lowered atmospheric CO2 levels. It is likely that during the Monterey Excursion the ocean/atmosphere system became progressively more sensitive to small changes in insolation, ultimately leading to major cooling of deep water and expansion of continental ice. We have assigned an absolute chronology, based on biostratigraphic and magneto-biostratigraphic datum levels, to the isotope stratigraphy and have used that chronology to correlate unconformities, seismic reflectors, carbonate minima, and dissolution intervals. Intervals of sediment containing 13C maxima are usually better preserved than the overlying and underlying sediments, indicating that the d13C values of TCO2 in deep water and the corrosiveness of seawater are inversely correlated. This again suggests that the 13C maxima were associated with rapid burial of organic carbon and reduced levels of atmospheric CO2. The absolute chronology we have assigned to the isotopic record indicates that the major mid-Miocene deepwater cooling/ice volume expansion took 2 m.y. and was not abrupt as had been reported previously. The cooling appears abrupt at many sites because the interval is characterized by a number of dissolution intervals. The cooling was not monotonic, and the 2 m.y. interval included an episode of especially rapid cooling as well as a brief return to warmer conditions before the final phase of the cooling period. The increase in d18O values of benthic foraminifera between 14.9 and 12.9 Ma was greatest at deeper water sites and at sites closest to Antarctica. The data suggest that the d18O value of seawater increased by no more than about 1.1 per mil during this interval and that the remainder of the change in benthic d18O values resulted from cooling in Antarctic regions of deepwater formation. Equatorial planktonic foraminifera from sites 237 and 289 exhibit a series of 0.4 per mil steplike increases in d13C values. Only one of these increases in planktonic d13C is correlated with any of the features in the mid-Miocene benthic carbon isotope record.