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.

Dynamik und Klimatologie von diabatischen Rossby-Wellen

Der Wintersturm Lothar zog am 26. Dezember 1999 über Europa und richtete in Frankreich, in Deutschland, in der Schweiz und in Österreich ungewöhnlich hohe Schäden an. Lothar entstand aus einer diabatischen Rossby Welle (DRW) und erreichte erst wenige Stunden vor dem europäischen Kontinent Orkanstärke. DRWs weisen ein interessantes atmosphärisches Strömungsmuster auf. Sie bestehen aus einer positiven PV-Anomalie in der unteren Troposphäre, die sich in einer Region mit starkem meridionalen Temperaturgradient befindet. Die positive PV-Anomalie löst eine zyklonale Strömung aus, dadurch wird östlich der PV-Anomalie warme Luft aus dem Süden herantransportiert. Während des Aufstieg der warmen Luft finden diabatische Prozesse statt, die zur Bildung einer neuen positiven PV-Anomalie in der unteren Troposphäre (PVA) führen. DRWs entstehen unabhängig von PV-Anomalien an der Tropopause. Falls sie jedoch mit ihnen in Wechselwirkung treten, kann - wie im Falle von Lothar - eine explosive Zyklogenese daraus resultieren. Im ersten Teil wird die Dynamik einer DRW am Beispiel des Wintersturms Lothar untersucht. Es wird insbesondere auf das Potential einer DRW zur explosiven Zyklogenese eingegangen. Im zweiten Teil wird das Aufretreten von DRWs in ECMWF-Vorhersagen untersucht. Es werden Unterschiede zwischen DRWs und anderen PV-Anomalien in der unteren Troposphäre hervorgehoben. Die Dynamik von DRWs wird mit Hilfe eines ECMWF-"Ensemble Prediction System" (EPS) des Wintersturms Lothar untersucht. Die 50 Modellläufe des EPS starten am 24. Dezember 1999 um 12 UTC und reichen bis zum 26. Dezember 1999 um 12 UTC. Nur 16 der 50 Modellläufe sagen einen ähnlich starken Sturm wie Lothar vorher. 10 Modellläufen sagen am 26. Dezember keine Zyklone mehr vorher. Die Ausprägung der baroklinen Zone, in der sich die DRW befindet, ist ausschlaggebend für die Intensität der DRW. Weitere wichtige Parameter sind der Feuchtegehalt der unteren Troposphäre und der latente Wärmefluss über dem Ozean. Diejenigen DRWs, die sich zu am 25. Dezember um 12 UTC näher als 400 km am Tropopausenjet befinden, entwickeln sich zu einer starken Zyklone. Alle anderen lösen sich auf oder bleiben schwache Zyklonen. Es ist schwierig, diabatische Prozesse in Wettervorhersagemodellen abzubilden, dementsprechend treten Schwierigkeiten bei der Vorhersage von PVAs auf. In den operationellen ECMWF-Vorhersagen von Juni 2004 bis Mai 2005 werden mit Hilfe eines Tracking- Algorithmus PVAs im Nordpazifik und Nordatlantik bestimmt und in fünf Kategorien eingeteilt. Die fünf Kategorien unterscheiden sich in ihrer Häufigkeit, ihrer Zugbahn und ihrer Gestalt. Im Nordpazifik entstehen doppelt so viele PVAs wie im Nordatlantik. Durchschnittlich werden im Winter weniger PVAs gefunden als im Sommer. Die Baroklinität und die Geschwindigkeit des Tropopausenjets ist in der Nähe von DRWs besonders hoch. Verglichen mit anderen PVAs weisen DRWs eine ähnliche Verteilung des reduzierten Bodendrucks auf. DRWs können in etwa gleich gut vorhergesagt werden wie andere PVAs.

Temporal variation of genetic diversity and differentiation of albacore, Thunnus alalunga, in the Mediterranean Sea

This study is on albacore (Thunnus alalunga, Bonnaterre 1788), an epi- and mesopelagic oceanic tuna species cosmopolitan in the tropical and temperate waters of all oceans including the Mediterranean Sea, extending in a broad band between 40°N and 40°S. What it’s known about albacore population structure is based on different studies that used fisheries data, RFLP, mtDNA control region and nuDNA markers, blood lectins analysis, individual tags and microsatellite. At the moment, for T. alalunga six management units are recognized: the North Pacific, South Pacific, Indian, North Atlantic, South Atlantic and Mediterranean stocks. In this study I have done a temporal and spatial comparison of genetic variability between different Mediterranean populations of Thunnus alalunga matching an historical dataset ca. from 1920s composed of 43 individuals divided in 3 populations (NADR, SPAIN and CMED) with a modern dataset composed of 254 individuals and 7 populations (BAL, CYP, LIG, TYR, TUR, ADR, ALB). The investigation was possible using a panel of 94 nuclear SNPs, built specifically for the target species at the University of Basque Country UPV/EHU. First analysis done was the Hardy-Weinberg, then the number of clusters (K) was determined using STRUCTURE and to assess the genetic variability, allele frequencies, the average number of alleles per locus, expected (He) and observed (Ho) heterozygosis, and the index of polymorphism (P) was used the software Genetix. Historical and modern samples gives different results, showing a clear loss of genetic diversity over time leading to a single cluster in modern albacore instead of the two found in historical samples. What this study reveals is very important for conservation concerns, and additional research endeavours are needed.

Historical weather extremes in the “Twentieth Century Reanalysis”

Meteorological or climatological extremes are rare and hence studying them requires long meteorological data sets. Moreover, for addressing the underlying atmospheric processes, detailed three-dimensional data are desired. Until recently the two requirements were incompatible as long meteorological series were only available for a few locations, whereas detailed 3-dimensional data sets such as reanalyses were limited to the past few decades. In 2011, the “Twentieth Century Reanalysis” (20CR) was released, a 6-hourly global atmospheric data set covering the past 140 years, thus combining the two properties. The collection of short papers in this volume contains case studies of individual extreme events in the 20CR data set. In this overview paper we introduce the first six cases and summarise some common findings. All of the events are represented in 20CR in a physically consistent way, allowing further meteorological interpretations and process studies. Also, for most of the events, the magnitudes are underestimated in the ensemble mean. Possible causes are addressed. For interpreting extrema it may be necessary to address individual ensemble members. Also, the density of observations underlying 20CR should be considered. Finally, we point to problems in wind speeds over the Arctic and the northern North Pacific in 20CR prior to the 1950s.

Last interglacial temperature evolution — a model inter-comparison

There is a growing number of proxy-based reconstructions detailing the climatic changes that occurred during the last interglacial period (LIG). This period is of special interest, because large parts of the globe were characterized by a warmer-than-present-day climate, making this period an interesting test bed for climate models in light of projected global warming. However, mainly because synchronizing the different palaeoclimatic records is difficult, there is no consensus on a global picture of LIG temperature changes. Here we present the first model inter-comparison of transient simulations covering the LIG period. By comparing the different simulations, we aim at investigating the common signal in the LIG temperature evolution, investigating the main driving forces behind it and at listing the climate feedbacks which cause the most apparent inter-model differences. The model inter-comparison shows a robust Northern Hemisphere July temperature evolution characterized by a maximum between 130–125 ka BP with temperatures 0.3 to 5.3 K above present day. A Southern Hemisphere July temperature maximum, −1.3 to 2.5 K at around 128 ka BP, is only found when changes in the greenhouse gas concentrations are included. The robustness of simulated January temperatures is large in the Southern Hemisphere and the mid-latitudes of the Northern Hemisphere. For these regions maximum January temperature anomalies of respectively −1 to 1.2 K and −0.8 to 2.1 K are simulated for the period after 121 ka BP. In both hemispheres these temperature maxima are in line with the maximum in local summer insolation. In a number of specific regions, a common temperature evolution is not found amongst the models. We show that this is related to feedbacks within the climate system which largely determine the simulated LIG temperature evolution in these regions. Firstly, in the Arctic region, changes in the summer sea-ice cover control the evolution of LIG winter temperatures. Secondly, for the Atlantic region, the Southern Ocean and the North Pacific, possible changes in the characteristics of the Atlantic meridional overturning circulation are crucial. Thirdly, the presence of remnant continental ice from the preceding glacial has shown to be important when determining the timing of maximum LIG warmth in the Northern Hemisphere. Finally, the results reveal that changes in the monsoon regime exert a strong control on the evolution of LIG temperatures over parts of Africa and India. By listing these inter-model differences, we provide a starting point for future proxy-data studies and the sensitivity experiments needed to constrain the climate simulations and to further enhance our understanding of the temperature evolution of the LIG period.

Are Greenhouse Gas Signals of Northern Hemisphere winter extra-tropical cyclone activity dependent on the identification and tracking algorithm?

For Northern Hemisphere extra-tropical cyclone activity, the dependency of a potential anthropogenic climate change signal on the identification method applied is analysed. This study investigates the impact of the used algorithm on the changing signal, not the robustness of the climate change signal itself. Using one single transient AOGCM simulation as standard input for eleven state-of-the-art identification methods, the patterns of model simulated present day climatologies are found to be close to those computed from re-analysis, independent of the method applied. Although differences in the total number of cyclones identified exist, the climate change signals (IPCC SRES A1B) in the model run considered are largely similar between methods for all cyclones. Taking into account all tracks, decreasing numbers are found in the Mediterranean, the Arctic in the Barents and Greenland Seas, the mid-latitude Pacific and North America. Changing patterns are even more similar, if only the most severe systems are considered: the methods reveal a coherent statistically significant increase in frequency over the eastern North Atlantic and North Pacific. We found that the differences between the methods considered are largely due to the different role of weaker systems in the specific methods.

State of the Antarctic and Southern Ocean Climate System

This paper reviews developments in our understanding of the state of the Antarctic and Southern Ocean climate and its relation to the global climate system over the last few millennia. Climate over this and earlier periods has not been stable, as evidenced by the occurrence of abrupt changes in atmospheric circulation and temperature recorded in Antarctic ice core proxies for past climate. Two of the most prominent abrupt climate change events are characterized by intensification of the circumpolar westerlies (also known as the Southern Annular Mode) between similar to 6000 and 5000 years ago and since 1200-1000 years ago. Following the last of these is a period of major trans-Antarctic reorganization of atmospheric circulation and temperature between A. D. 1700 and 1850. The two earlier Antarctic abrupt climate change events appear linked to but predate by several centuries even more abrupt climate change in the North Atlantic, and the end of the more recent event is coincident with reorganization of atmospheric circulation in the North Pacific. Improved understanding of such events and of the associations between abrupt climate change events recorded in both hemispheres is critical to predicting the impact and timing of future abrupt climate change events potentially forced by anthropogenic changes in greenhouse gases and aerosols. Special attention is given to the climate of the past 200 years, which was recorded by a network of recently available shallow firn cores, and to that of the past 50 years, which was monitored by the continuous instrumental record. Significant regional climate changes have taken place in the Antarctic during the past 50 years. Atmospheric temperatures have increased markedly over the Antarctic Peninsula, linked to nearby ocean warming and intensification of the circumpolar westerlies. Glaciers are retreating on the peninsula, in Patagonia, on the sub-Antarctic islands, and in West Antarctica adjacent to the peninsula. The penetration of marine air masses has become more pronounced over parts of West Antarctica. Above the surface, the Antarctic troposphere has warmed during winter while the stratosphere has cooled year-round. The upper kilometer of the circumpolar Southern Ocean has warmed, Antarctic Bottom Water across a wide sector off East Antarctica has freshened, and the densest bottom water in the Weddell Sea has warmed. In contrast to these regional climate changes, over most of Antarctica, near-surface temperature and snowfall have not increased significantly during at least the past 50 years, and proxy data suggest that the atmospheric circulation over the interior has remained in a similar state for at least the past 200 years. Furthermore, the total sea ice cover around Antarctica has exhibited no significant overall change since reliable satellite monitoring began in the late 1970s, despite large but compensating regional changes. The inhomogeneity of Antarctic climate in space and time implies that recent Antarctic climate changes are due on the one hand to a combination of strong multidecadal variability and anthropogenic effects and, as demonstrated by the paleoclimate record, on the other hand to multidecadal to millennial scale and longer natural variability forced through changes in orbital insolation, greenhouse gases, solar variability, ice dynamics, and aerosols. Model projections suggest that over the 21st century the Antarctic interior will warm by 3.4 degrees +/- 1 degrees C, and sea ice extent will decrease by similar to 30%. Ice sheet models are not yet adequate enough to answer pressing questins about the effect of projected warming on mass balance and sea level. Considering the potentially major impacts of a warming climate on Antarctica, vigorous efforts are needed to better understand all aspects of the highly coupled Antarctic climate system as well as its influence on the Earth's climate and oceans.

Ice Core Paleovolcanic Records from the St. Elias Mountains, Yukon, Canada

We previously reported a record of regionally significant volcanic eruptions in the North Pacific using an ice core from Eclipse Icefield (St. Elias Mountains, Yukon, Canada). The acquisition of two new ice cores from Eclipse Icefield, along with the previously available Eclipse Icefield and Mount Logan Northwest Col ice cores, allows us to extend our record of North Pacific volcanism to 550 years before present using a suite of four ice cores spanning an elevation range of 3 - 5 km. Comparison of volcanic sulfate flux records demonstrates that the results are highly reproducible, especially for the largest eruptions such as Katmai ( A. D. 1912). Correlation of volcanic sulfate signals with historically documented eruptions indicates that at least one-third of the eruptions recorded in St. Elias ice cores are from Alaskan and Kamchatkan volcanoes. Although there are several moderately large ( volcanic explosivity index (VEI) >= 4) eruptions recorded in only one core from Eclipse Icefield, the use of multiple cores provides signals in at least one core from all known VEI >= 4 eruptions in Alaska and Kamchatka since A. D. 1829. Tephrochronological evidence from the Eclipse ice cores documents eruptions in Alaska (Westdahl, Redoubt, Trident, and Katmai), Kamchatka (Avachinsky, Kliuchevoskoi, and Ksudach), and Iceland (Hekla). Several unidentified tephra-bearing horizons, with available geochemical evidence suggesting Alaskan and Kamchatkan sources, were also found. We present a reconstruction of annual volcanic sulfate loading for the North Pacific troposphere based on our ice core data, and we provide a detailed assessment of the atmospheric and climatic effects of the Katmai eruption.

Kelp Forest Ecosystems: Biodiversity, Stability, Resilience and Future

Kelp forests are phyletically diverse, structurally complex and highly productive components of cold-water rocky marine coastlines. This paper reviews the conditions in which kelp forests develop globally and where, why and at what rate they become deforested. The ecology and long archaeological history of kelp forests are examined through case studies from southern California, the Aleutian Islands and the western North Atlantic, well-studied locations that represent the widest possible range in kelp forest biodiversity. Global distribution of kelp forests is physiologically constrained by light at high latitudes and by nutrients, warm temperatures and other macrophytes at low latitudes. Within mid-latitude belts (roughly 40-60degrees latitude in both hemispheres) well-developed kelp forests are most threatened by herbivory, usually from sea urchins. Overfishing and extirpation of highly valued vertebrate apex predators often triggered herbivore population increases, leading to widespread kelp deforestation. Such deforestations have the most profound and lasting impacts on species-depauperate systems, such as those in Alaska and the western North Atlantic. Globally urchin-induced deforestation has been increasing over the past 2-3 decades. Continued fishing down of coastal food webs has resulted in shifting harvesting targets from apex predators to their invertebrate prey, including kelp-grazing herbivores. The recent global expansion of sea urchin harvesting has led to the widespread extirpation of this herbivore, and kelp forests have returned in some locations but, for the first time, these forests are devoid of vertebrate apex predators. In the western North Atlantic, large predatory crabs have recently filled this void and they have become the new apex predator in this system. Similar shifts from fish- to crab-dominance may have occurred in coastal zones of the United Kingdom and Japan, where large predatory finfish were extirpated long ago. Three North American case studies of kelp forests were examined to determine their long history with humans and project the status of future kelp forests to the year 2025. Fishing impacts on kelp forest systems have been both profound and much longer in duration than previously thought. Archaeological data suggest that coastal peoples exploited kelp forest organisms for thousands of years, occasionally resulting in localized losses of apex predators, outbreaks of sea urchin populations and probably small-scale deforestation. Over the past two centuries, commercial exploitation for export led to the extirpation of sea urchin predators, such as the sea otter in the North Pacific and predatory fishes like the cod in the North Atlantic. The largescale removal of predators for export markets increased sea urchin abundances and promoted the decline of kelp forests over vast areas. Despite southern California having one of the longest known associations with coastal kelp forests, widespread deforestation is rare. It is possible that functional redundancies among predators and herbivores make this most diverse system most stable. Such biodiverse kelp forests may also resist invasion from non-native species. In the species-depauperate western North Atlantic, introduced algal competitors carpet the benthos and threaten future kelp dominance. There, other non-native herbivores and predators have become established and dominant components of this system. Climate changes have had measurable impacts on kelp forest ecosystems and efforts to control the emission of greenhouse gasses should be a global priority. However, overfishing appears to be the greatest manageable threat to kelp forest ecosystems over the 2025 time horizon. Management should focus on minimizing fishing impacts and restoring populations of functionally important species in these systems.

On the Co-Occurrence of Warm Conveyor Belt Outflows and PV Streamers*

The co-occurrence of warm conveyor belts (WCBs), strongly ascending moist airstreams in extratropical cyclones, and stratospheric potential vorticity (PV) streamers, indicators for breaking Rossby waves on the tropopause, is investigated for a 21-yr period in the Northern Hemisphere using Interim European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-Interim) data. WCB outflows and PV streamers are respectively identified as two- and three-dimensional objects and tracked during their life cycle. PV streamers are more frequent than WCB outflows and nearly 15% of all PV streamers co-occur with WCBs during their life cycle, whereas about 60% of all WCB outflows co-occur with PV streamers. Co-occurrences are most frequent over the North Atlantic and North Pacific in spring and winter. WCB outflows are often located upstream of the PV streamers and form earlier, indicating the importance of diabatic processes for downstream Rossby wave breaking. Less frequently, PV streamers occur first, leading to the formation of new WCBs.

Impact of oceanic circulation changes on atmospheric δ¹³ CO₂

δ¹³ CO₂ measured in Antarctic ice cores provides constraints on oceanic and terrestrial carbon cycle processes linked with millennial-scale changes in atmospheric CO₂. However, the interpretation of δ¹³ CO₂ is not straight-forward. Using carbon isotope-enabled versions of the LOVECLIM and Bern3D models, we perform a set of sensitivity experiments in which the formation rates of North Atlantic Deep Water (NADW), North Pacific Deep Water (NPDW), Antarctic Bottom Water (AABW), and Antarctic Intermediate Water (AAIW) are varied. We study the impact of these circulation changes on atmospheric δ¹³ CO₂ as well as on the oceanic δ¹³ CO₂ distribution. In general, we find that the formation rates of AABW, NADW, NPDW, and AAIW are negatively correlated with changes in δ¹³ CO₂: namely, strong oceanic ventilation decreases atmospheric δ¹³ CO₂. However, since large-scale oceanic circulation reorganizations also impact nutrient utilization and the Earth’s climate, the relationship between atmospheric δ¹³ CO₂ levels and ocean ventilation rate is not unequivocal. In both models atmospheric δ¹³ CO₂ is very sensitive to changes in AABW formation rates: increased AABW formation enhances the transport of low δ¹³ CO₂ waters to the surface and decreases atmospheric δ¹³ CO₂. By contrast, the impact of NADW changes on atmospheric δ¹³ CO₂ is less robust and might be model dependent. This results from complex interplay between global climate, carbon cycle, and the formation rate of NADW, a water body characterized by relatively high δ¹³ CO₂.

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.