Impact of geomagnetic excursions on atmospheric chemistry and dynamics

Geomagnetic excursions, i.e. short periods in time with much weaker geomagnetic fields and substantial changes in the position of the geomagnetic pole, occurred repeatedly in the Earth's history, e.g. the Laschamp event about 41 kyr ago. Although the next such excursion is certain to come, little is known about the timing and possible consequences for the state of the atmosphere and the ecosystems. Here we use the global chemistry climate model SOCOL-MPIOM to simulate the effects of geomagnetic excursions on atmospheric ionization, chemistry and dynamics. Our simulations show significantly increased concentrations of nitrogen oxides (NOx) in the entire stratosphere, especially over Antarctica (+15%), due to enhanced ionization by galactic cosmic rays. Hydrogen oxides (HOx) are also produced in greater amounts (up to +40%) in the tropical and subtropical lower stratosphere, while their destruction by reactions with enhanced NOx prevails over the poles and in high altitudes (by −5%). Stratospheric ozone concentrations decrease globally above 20 km by 1–2% and at the northern hemispheric tropopause by up to 5% owing to the accelerated NOx-induced destruction. A 5% increase is found in the southern lower stratosphere and troposphere. In response to these changes in ozone and the concomitant changes in atmospheric heating rates, the Arctic vortex intensifies in boreal winter, while the Antarctic vortex weakens in austral winter and spring. Surface wind anomalies show significant intensification of the southern westerlies at their poleward edge during austral winter and a pronounced northward shift in spring. Major impacts on the global climate seem unlikely.

The coupled atmosphere–chemistry–ocean model SOCOL-MPIOM

The newly developed atmosphere–ocean-chemistry-climate model SOCOL-MPIOM is presented by demonstrating the influence of the interactive chemistry module on the climate state and the variability. Therefore, we compare pre-industrial control simulations with (CHEM) and without (NOCHEM) interactive chemistry. In general, the influence of the chemistry on the mean state and the variability is small and mainly restricted to the stratosphere and mesosphere. The largest differences are found for the atmospheric dynamics in the polar regions, with slightly stronger northern and southern winter polar vortices in CHEM. The strengthening of the vortex is related to larger stratospheric temperature gradients, which are attributed to a parametrization of the absorption of ozone and oxygen in the Lyman-alpha, Schumann–Runge, Hartley, and Higgins bands. This effect is parametrized in the version with interactive chemistry only. A second reason for the temperature differences between CHEM and NOCHEM is related to diurnal variations in the ozone concentrations in the higher atmosphere, which are missing in NOCHEM. Furthermore, stratospheric water vapour concentrations differ substantially between the two experiments, but their effect on the temperatures is small. In both setups, the simulated intensity and variability of the northern polar vortex is inside the range of present day observations. Sudden stratospheric warming events are well reproduced in terms of their frequency, but the distribution amongst the winter months is too uniform. Additionally, the performance of SOCOL-MPIOM under changing external forcings is assessed for the period 1600–2000 using an ensemble of simulations driven by a spectral solar forcing reconstruction. The amplitude of the reconstruction is large in comparison to other state-of-the-art reconstructions, providing an upper limit for the importance of the solar signal. In the pre-industrial period (1600–1850) the simulated surface temperature trends are in reasonable agreement with temperature reconstructions, although the multi-decadal variability is more pronounced. This enhanced variability can be attributed to the variability in the solar forcing. The simulated temperature reductions during the Maunder Minimum are in the lowest probability range of the proxy records. During the Dalton Minimum, when also volcanic forcing is an important driver of temperature variations, the agreement is better. In the industrial period from 1850 onward SOCOL-MPIOM overestimates the temperature increase in comparison to observational data sets. Sensitivity simulations show that this overestimation can be attributed to the increasing trend in the solar forcing reconstruction that is used in this study and an additional warming induced by the simulated ozone changes.

The diurnal variation in stratospheric ozone from the MACC reanalysis, the ERA-Interim reanalysis, WACCM and Earth observations: characteristics and intercomparison

In this study we compare the diurnal variation in stratospheric ozone derived from free-running simulations of the Whole Atmosphere Community Climate Model (WACCM) and from reanalysis data of the atmospheric service MACC (Monitoring Atmospheric Composition and Climate) which both use a similar stratospheric chemistry module. We find good agreement between WACCM and the MACC reanalysis for the diurnal ozone variation in the high-latitude summer stratosphere based on photochemistry. In addition, we consult the ozone data product of the ERA-Interim reanalysis. The ERA-Interim reanalysis ozone system with its long-term ozone parametrization can not capture these diurnal variations in the upper stratosphere that are due to photochemistry. The good dynamics representations, however, reflects well dynamically induced ozone variations in the lower stratosphere. For the high-latitude winter stratosphere we describe a novel feature of diurnal variation in ozone where changes of up to 46.6% (3.3 ppmv) occur in monthly mean data. For this effect good agreement between the ERA-Interim reanalysis and the MACC reanalysis suggest quite similar diurnal advection processes of ozone. The free-running WACCM model seriously underestimates the role of diurnal advection processes at the polar vortex at the two tested resolutions. The intercomparison of the MACC reanalysis and the ERA-Interim reanalysis demonstrates how global reanalyses can benefit from a chemical representation held by a chemical transport model. The MACC reanalysis provides an unprecedented description of the dynamics and photochemistry of the diurnal variation of stratospheric ozone which is of high interest for ozone trend analysis and research on atmospheric tides. We confirm the diurnal variation in ozone at 5 hPa by observations of the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) experiment and selected sites of the Network for Detection of Atmospheric Composition Change (NDACC). The latter give valuable insight even to diurnal variation of ozone in the polar winter stratosphere.

Signatures of the 2-day wave and sudden stratospheric warmings in Arctic water vapour observed by ground-based microwave radiometry

The ground-based microwave radiometer MIAWARA-C recorded the upper stratospheric and lower mesospheric water vapour distribution continuously from June 2011 to March 2013 above the Arctic station of Sodankylä, Finland (67.4° N, 26.6° E) without major interruptions and offers water vapour profiles with temporal resolution of 1 h for average conditions. The water vapour time series of MIAWARA-C shows strong periodic variations in both summer and winter related to the quasi-2-day wave. Above 0.1 hPa the amplitudes are strongest in summer. The stratospheric wintertime 2-day wave is pronounced for both winters on altitudes below 0.1 hPa and reaches a maximum amplitude of 0.8 ppmv in November 2011. Over the measurement period, the instrument monitored the changes in water vapour linked to two sudden stratospheric warmings in early 2012 and 2013. Based on the water vapour measurements, the descent rate in the vortex after the warmings is 364 m d−1 for 2012 and 315 m d−1 for 2013.

Trajectory mapping of middle atmospheric water vapor by a mini network of NDACC instruments

The important task to observe the global coverage of middle atmospheric trace gases like water vapor or ozone usually is accomplished by satellites. Climate and atmospheric studies rely upon the knowledge of trace gas distributions throughout the stratosphere and mesosphere. Many of these gases are currently measured from satellites, but it is not clear whether this capability will be maintained in the future. This could lead to a significant knowledge gap of the state of the atmosphere. We explore the possibilities of mapping middle atmospheric water vapor in the Northern Hemisphere by using Lagrangian trajectory calculations and water vapor profile data from a small network of five ground-based microwave radiometers. Four of them are operated within the frame of NDACC (Network for the Detection of Atmospheric Composition Change). Keeping in mind that the instruments are based on different hardware and calibration setups, a height-dependent bias of the retrieved water vapor profiles has to be expected among the microwave radiometers. In order to correct and harmonize the different data sets, the Microwave Limb Sounder (MLS) on the Aura satellite is used to serve as a kind of traveling standard. A domain-averaging TM (trajectory mapping) method is applied which simplifies the subsequent validation of the quality of the trajectory-mapped water vapor distribution towards direct satellite observations. Trajectories are calculated forwards and backwards in time for up to 10 days using 6 hourly meteorological wind analysis fields. Overall, a total of four case studies of trajectory mapping in different meteorological regimes are discussed. One of the case studies takes place during a major sudden stratospheric warming (SSW) accompanied by the polar vortex breakdown; a second takes place after the reformation of stable circulation system. TM cases close to the fall equinox and June solstice event from the year 2012 complete the study, showing the high potential of a network of ground-based remote sensing instruments to synthesize hemispheric maps of water vapor.

Trajectory mapping of middle atmospheric water vapor by a mini network of NDACC instruments

The important task to observe the global coverage of middle atmospheric trace gases like water vapor or ozone usually is accomplished by satellites. Climate and atmospheric studies rely upon the knowledge of trace gas distributions throughout the stratosphere and mesosphere. Many of these gases are currently measured from satellites, but it is not clear whether this capability will be maintained in the future. This could lead to a significant knowledge gap of the state of the atmosphere. We explore the possibilities of mapping middle atmospheric water vapor in the Northern Hemisphere by using Lagrangian trajectory calculations and water vapor profile data from a small network of five ground-based microwave radiometers. Four of them are operated within the frame of NDACC (Network for the Detection of Atmospheric Composition Change). Keeping in mind that the instruments are based on different hardware and calibration setups, a height-dependent bias of the retrieved water vapor profiles has to be expected among the microwave radiometers. In order to correct and harmonize the different data sets, the Microwave Limb Sounder (MLS) on the Aura satellite is used to serve as a kind of traveling standard. A domain-averaging TM (trajectory mapping) method is applied which simplifies the subsequent validation of the quality of the trajectory-mapped water vapor distribution towards direct satellite observations. Trajectories are calculated forwards and backwards in time for up to 10 days using 6 hourly meteorological wind analysis fields. Overall, a total of four case studies of trajectory mapping in different meteorological regimes are discussed. One of the case studies takes place during a major sudden stratospheric warming (SSW) accompanied by the polar vortex breakdown; a second takes place after the reformation of stable circulation system. TM cases close to the fall equinox and June solstice event from the year 2012 complete the study, showing the high potential of a network of ground-based remote sensing instruments to synthesize hemispheric maps of water vapor.

Investigation of hemodynamics in an in vitro system simulating left ventricular support through the right subclavian artery using 4-dimensional flow magnetic resonance imaging.

OBJECTIVES Left ventricular assist devices are an important treatment option for patients with heart failure alter the hemodynamics in the heart and great vessels. Because in vivo magnetic resonance studies of patients with ventricular assist devices are not possible, in vitro models represent an important tool to investigate flow alterations caused by these systems. By using an in vitro magnetic resonance-compatible model that mimics physiologic conditions as close as possible, this work investigated the flow characteristics using 4-dimensional flow-sensitive magnetic resonance imaging of a left ventricular assist device with outflow via the right subclavian artery as commonly used in cardiothoracic surgery in the recent past. METHODS An in vitro model was developed consisting of an aorta with its supra-aortic branches connected to a left ventricular assist device simulating the pulsatile flow of the native failing heart. A second left ventricular assist device supplied the aorta with continuous flow via the right subclavian artery. Four-dimensional flow-sensitive magnetic resonance imaging was performed for different flow rates of the left ventricular assist device simulating the native heart and the left ventricular assist device providing the continuous flow. Flow characteristics were qualitatively and quantitatively evaluated in the entire vessel system. RESULTS Flow characteristics inside the aorta and its upper branching vessels revealed that the right subclavian artery and the right carotid artery were solely supported by the continuous-flow left ventricular assist device for all flow rates. The flow rates in the brain-supplying arteries are only marginally affected by different operating conditions. The qualitative analysis revealed only minor effects on the flow characteristics, such as weakly pronounced vortex flow caused by the retrograde flow via the brachiocephalic artery. CONCLUSIONS The results indicate that, despite the massive alterations in natural hemodynamics due to the retrograde flow via the right subclavian and brachiocephalic arteries, there are no drastic consequences on the flow in the brain-feeding arteries and the flow characteristics in the ascending and descending aortas. It may be beneficial to adjust the operating condition of the left ventricular assist device to the residual function of the failing heart.

Meteorological measurements from the Fourcade and Polar Club Glacier, Warszawa Icefield, King George Island, West Antarctica

The South Shetland Islands are located at the northern tip of the AP which is among the fastest warming regions on Earth. The islands are especially vulnerable to climate change due to their exposure to transient low-pressure systems and their maritime climate. Surface air temperature increases (2.5K in 50 years) are concurrent with retreating glacier fronts, an increase in melt areas, ice surface lowering and rapid break-up and disintegration of ice shelves. We have compiled a unique meteorological data set for the King George Island (KGI)/Isla 25 de Mayo, the largest of the South Shetland Islands. It comprises high-temporal resolution and spatially distributed observations of surface air temperature, wind directions and wind velocities, as well as glacier ice temperatures in profile with a fully equipped automatic weather station on the Warszawa Icefield, from November 2010 and ongoing. In combination with two long-term synoptic datasets (40 and 10 years, respectively) and NCEP/NCAR reanalysis data, we have looked at changes in the climatological drivers of the glacial melt processes, and the sensitivity of the inland ice cap with regard to winter melting periods and pressure anomalies. The analysis has revealed, a positive trend of 5K over four decades in minimum surface air temperatures for winter months, clearly exceeding the published annual mean statistics, associated to a decrease in mean monthly winter sea level pressure. This concurs with a positive trend in the Southern Annular Mode (SAM) index, which gives a measure for the strength and extension of the Antarctic vortex. We connect this trend with a higher frequency of low-pressure systems hitting the South Shetland Islands during austral winter, bringing warm and moist air masses from lower latitudes. Due to its exposure to the impact of transient synoptic weather systems, the ice cap of KGI is especially vulnerable to changes during winter glacial mass accumulation period. A revision of seasonal changes in adiabatic air temperature lapse rates and their dependency on exposure and elevation has shown a clear decoupling of atmospheric surface layers between coastal areas and the higher-elevation ice cap, showing the higher sensitivity to free atmospheric flow and synoptic changes. Observed surface air temperature lapse rates show a high variability during winter months (standard deviations up to ±1.0K/100 m), and a distinct spatial variability reflecting the impact of synoptic weather patterns. The observed advective conditions bringing warm, moist air with high temperatures and rain, lead to melt conditions on the ice cap, fixating surface air temperatures to the melting point. This paper assesses the impact of large-scale atmospheric circulation variability and climatic changes on the atmospheric surface layer and glacier mass accumulation of the upper ice cap during winter season for the Warszawa Icefield on KGI.

Biological, physical and chemical characteristics of water and sediment samples from the Volga River delta, northern Caspian

Based on the data of synchronous observations of hydrophysical and biogeochemical parameters in the near-mouth and shallow-water areas of the northern Caspian in 2000-2001, the scale of spatiotemporal variability in the following characteristics of the water-bottom system was estimated (1) flow velocity and direction within vortex structures formed by the combined effect of wind, discharge current, and the presence of higher aquatic plants; (2) dependence of the spatial distribution of the content and composition of suspended particulate matter on the hydrodynamic regime of waters and development of phytoplankton; (3) variations in the grain-size, petrographic, mineralogical, and chemical compositions of the upper layer of bottom sediments at several sites in the northern Caspian related to the particular local combination of dominant natural processes; and (4) limits of variability in the group composition of humus compounds in bottom sediments. The acquired data are helpful in estimating the geochemical consequences of a sea level rise and during the planning of preventive environmental protection measures in view of future oil and gas recovery in this region.

Maestrichtian through Neogene benthic foraminifers of Maud Rise

Upper abyssal to lower bathyal benthic foraminifers from ODP Sites 689 (present water depth 2080 m) and 690 (present water depth 2941 m) on Maud Rise (eastern Weddell Sea, Antarctica) are reliable indicators of Maestrichtian through Neogene changes in the deep-water characteristics at high southern latitudes. Benthic foraminiferal faunas were divided into eight assemblages, with periods of faunal change at the early/late Maestrichtian boundary (69 Ma), at the early/late Paleocene boundary (62 Ma), in the latest Paleocene (57.5 Ma), in the middle early Eocene to late early Eocene (55-52 Ma), in the middle middle Eocene (46 Ma), in the late Eocene (38.5 Ma), and in the middle-late Miocene (14.9-11.5 Ma). These periods of faunal change may have occurred worldwide at the same time, although specific first and last appearances of deep-sea benthic foraminifers are commonly diachronous. There were minor faunal changes at the Cretaceous/Tertiary boundary (less than 14?7o of the species had last appearances at Site 689, less than 9% at Site 690). The most abrupt benthic foraminiferal faunal event occurred in the latest Paleocene, when the diversity dropped by 50% (more than 35% of species had last appearances) over a period of less than 25,000 years; after the extinction the diversity remained low for about 350,000 years. The highest diversities of the post-Paleocene occurred during the middle Eocene; from that time on the diversity decreased steadily at both sites. Data on faunal composition (percentage of infaunal versus epifaunal species) suggest that the waters bathing Maud Rise were well ventilated during the Maestrichtian through early Paleocene as well as during the latest Eocene through Recent. The waters appeared to be less well ventilated during the late Paleocene as well as the late middle through early late Eocene, with the least degree of ventilation during the latest Paleocene through early Eocene. The globally recognized extinction of deep-sea benthic foraminifers in the latest Paleocene may have been caused by a change in formational processes of the deep to intermediate waters of the oceans: from formation of deep waters by sinking at high latitudes to formation of deep to intermediate water of the oceans by evaporation at low latitudes. Benthic foraminiferal data (supported by carbon and oxygen isotopic data) suggest that there was a short period of intense formation of warm, salty deep water at the end of the Paleocene (with a duration of about 0.35 m.y.), and that less intense, even shorter episodes might have occurred during the late Paleocene and early Eocene. The faunal record from the Maud Rise sites agrees with published faunal and isotopic records, suggesting cooling of deep to intermediate waters in the middle through late Eocene.

Benthic foraminifera in Tertiary sediments of DSDP Leg 90 holes

Eocene through Pliocene benthic foraminifers were examined from seven sites located at middle and lower bathyal depths on the Lord Howe Rise in the Tasman Sea, from another site at lower bathyal depths in the Coral Sea, and from a site in the intermediate-depth, hemipelagic province of the Chatham Rise, east of southern New Zealand. Age-related, depth-related, and bioprovincial faunal variations are documented in this chapter. One new species, Rectuvigerina tasmana, is named. The paleoecologic indications of several key groups, including the miliolids, uvigerinids, nuttallitids, and cibicidids, are combined with sedimentologic and stable isotopic tracers to interpret paleoceanographic changes in the Tasman Sea. Because the total stratigraphic ranges of many bathyal benthic foraminifers are not yet known, most endpoints in the Tasman Sea are considered ecologically controlled events. The disappearances of Uvigerina rippensis and Cibicidoidesparki and the first appearances of U. pigmaea, Sphaeroidina bulloides, and Rotaliatina sulcigera at the Eocene/Oligocene boundary can be considered evolutionary events, as also can the first appearance of Cibicides wuellerstorfi in Zone NN5. Species which are restricted to the lower bathyal zone except during discrete pulses, most of which are related to the development of glacial conditions, include Melonis pompilioides, M. sphaeroides, Pullenia quinqueloba, Nuttallides umbonifera, and U. hispido-costata. Middle bathyal indigenes include U. spinulosa, U. gemmaeformis, Ehrenbergina marwicki, R. sulcigera, and all rectuvigerinids except Rectuvigerina spinea. Although the miliolids first occurred at lower bathyal depths, they were more common in the middle bathyal zone. Although the Neogene hispido-costate uvigerinids first developed at lower bathyal depths and at higher middle latitude sites, in the later Neogene this group migrated to shallower depths and became predominant also in the middle bathyal zone. Despite the relatively similar sedimentologic settings at the six middle bathyal Tasman sites, there was extensive intrageneric and intraspecific geographic variation. Mililiolids, strongly ornamented brizalinids, bolivinitids, Bulimina aculeata, Osangularia culter, and strongly porous morphotypes were more common at higher latitudes. Osangularia bengalensis, striate brizalinids such as Brizalina subaenariensis, Gaudryina solida, osangularids in general, and finely porous morphotypes were more common in the subtropics. There was strong covariance between faunas at lower middle latitude, lower bathyal Site 591, and higher middle latitude, middle bathyal Site 593. The following oceanographic history of the Tasman Sea is proposed; using the stable isotopic record as evidence for glacials and examining the ecologic correlations between (1) miliolids and carbonate saturation, (2) nuttallitids and undersaturated, cooled, or "new" water masses, (3) uvigerinids with high organic carbon in the sediment and high rates of sediment accumulation, and (4) cibicidids and terrestrial organic carbon. The glacial located near the Eocene/Oligocene boundary is characterized by the penetration of cooler, more corrosive waters at intermediate depths in high southern latitudes. This may have caused overturn, upwelling pulses, in other Tasman areas. The development of Neogenelike conditions began in the late Oligocene (Zone NP24/NP25) with the evolution of several common Neogene species. A large number of Paleogene benthics disappeared gradually through the course of the early Miocene, which was not well preserved at any Tasman site. Corrosive conditions shallowed into the middle bathyal zone in several pulses during the early Miocene. The development of glacial conditions in the middle Miocene was accompanied by major changes throughout the Tasman Sea. Sediment accumulation rates increased and high-productivity faunas and corrosive conditions developed at all but the lowest-latitude Site 588. This increase in productivity and accumulation rate is attributed to the eutrophication of Antarctic water masses feeding Tasman current systems, as well as to invigorated circulation in general. It overlaps with the beginning of the Pacific High-productivity Episode (10-5 Ma). During the latest Miocene glacial episode, corrosive conditions developed at lower bathyal depths, while cooler water and lower nutrient levels shallowed to middle bathyal depths. Lower input of terrestrial organic carbon may be related to the lower nutrient levels of this time and to the termination of the Pacific High-productivity Episode. The moderate glacial episode during the mid-Pliocene (Zone NN15/NN16, ~3.2 Ma) corresponds to a decline in sediment accumulation rates and a reorganization of faunas unlike that of all other times. New genera proliferate and indices for cool, noncorrosive conditions and high organic carbon expand throughout the middle bathyal zone coeval with the sedimentation rate decreases. By the latest Pliocene (about 2.5 Ma), however, during another glacial episode, faunal patterns typical of this and later glacials develop throughout the Tasman Sea. Benthic foraminiferal patterns suggest increased input of terrestrial organic matter to Tasman Sea sediments during this episode and during later glacials.

Benthic foraminifera and stable isotope composition in Paleocene-Eocene sediments

In the late Paleocene to early Eocene, deep sea benthic foraminifera suffered their only global extinction of the last 75 million years and diversity decreased worldwide by 30-50% in a few thousand years. At Maud Rise (Weddell Sea, Antarctica; Sites 689 and 690, palaeodepths 1100 m and 1900 m) and Walvis Ridge (Southeastern Atlantic, Sites 525 and 527, palaeodepths 1600 m and 3400 m) post-extinction faunas were low-diversity and high-dominance, but the dominant species differed by geographical location. At Maud Rise, post-extinction faunas were dominated by small, biserial and triserial species, while the large, thick-walled, long-lived deep sea species Nuttallides truempyi was absent. At Walvis Ridge, by contrast, they were dominated by long-lived species such as N. truempyi, with common to abundant small abyssaminid species. The faunal dominance patterns at the two locations thus suggest different post-extinction seafloor environments: increased flux of organic matter and possibly decreased oxygen levels at Maud Rise, decreased flux at Walvis Ridge. The species-richness remained very low for about 50 000 years, then gradually increased. The extinction was synchronous with a large, negative, short-term excursion of carbon and oxygen isotopes in planktonic and benthic foraminifera and bulk carbonate. The isotope excursions reached peak negative values in a few thousand years and values returned to pre-excursion levels in about 50 000 years. The carbon isotope excursion was about -2 per mil for benthic foraminifera at Walvis Ridge and Maud Rise, and about -4 per mil for planktonic foraminifera at Maud Rise. At the latter sites vertical gradients thus decreased, possibly at least partially as a result of upwelling. The oxygen isotope excursion was about -1.5 per mil for benthic foraminifera at Walvis Ridge and Maud Rise, -1 per mil for planktonic foraminifera at Maud Rise. The rapid oxygen isotope excursion at a time when polar ice-sheets were absent or insignificant can be explained by an increase in temperature by 4-6°C of high latitude surface waters and deep waters world wide. The deep ocean temperature increase could have been caused by warming of surface waters at high latitudes and continued formation of the deep waters at these locations, or by a switch from dominant formation of deep waters at high latitudes to formation at lower latitudes. Benthic foraminiferal post-extinction biogeographical patterns favour the latter explanation. The short-term carbon isotope excursion occurred in deep and surface waters, and in soil concretions and mammal teeth in the continental record. It is associated with increased CaC03-dissolution over a wide depth range in the oceans, suggesting that a rapid transfer of isotopically light carbon from lithosphere or biosphere into the ocean-atmosphere system may have been involved. The rapidity of the initiation of the excursion (a few thousand years) and its short duration (50 000 years) suggest that such a transfer was probably not caused by changes in the ratio of organic carbon to carbonate deposition or erosion. Transfer of carbon from the terrestrial biosphere was probably not the cause, because it would require a much larger biosphere destruction than at the end of the Cretaceous, in conflict with the fossil record. It is difficult to explain the large shift by rapid emission into the atmosphere of volcanogenic CO2, although huge subaerial plateau basalt eruptions occurred at the time in the northern Atlantic. Probably a complex combination of processes and feedback was involved, including volcanogenic emission of CO2, changing circulation patterns, changing productivity in the oceans and possibly on land, and changes in the relative size of the oceanic and atmospheric carbon reservoirs.