The expression of the end-Cenomanian oceanic anoxic event (OAE 2) in the Helvetic Alps

Upper Cenomanian pelagic sediments from the northern Alpine Helvetic fold-and-thrust belt (northern Tethyan margin) coeval with Oceanic Anoxic Event (OAE) 2 are characterized by the temporal persistence of micrite sedimentation and lack of organic carbon-rich layers. We studied an expanded section in the Chrummflueschlucht (east of Euthal, Switzerland), which encompasses the OAE 2 time interval. In order to identify the paleoceanographic and paleoenvironmental conditions during OAE 2 in this part of the northern Tethyan margin, and more specifically to trace eventual changes in nutrient levels and oxic conditions, we investigated the biostratigraphy (planktonic foraminifera), the bulk-rock mineralogy, and measured stable carbon- and oxygen-isotopes, total phosphorus (P) and redox-sensitive trace-element (RSTE) contents. We were able to determine – with some remaining uncertainties – the different planktonic foraminiferal biozones characteristic of the Cenomanian–Turonian boundary interval (Rotalipora cushmani, Whiteinella archaeocretacea and Helvetoglobotruncana helvetica zones). In the lower part of the section (R. cushmani total range zone), the bulk-rock δ13C record shows a long-term increase. Within sediments attributed to the W. archaeocretacea partial range zone, δ13C values reach a maximum of 3.3‰ (peak “a”). In the following the values decrease and increase again to arrive at a plateau with high δ13C values of around 3.1‰, which ends with a peak of 3.3‰ (peak “c”). At the top of the section, in sediments belonging to the H. helvetica total range zone, δ13C values decrease to post-OAE values of around 2.2‰. The last occurrence of R. cushmani is observed just above the positive δ13C shift characterizing OAE 2. P contents display small variations along the section with a long-term decreasing trend towards the top. Before the OAE 2 interval, P values show higher values and relatively good covariation with detrital input, indicating higher nutrient input before OAE 2. In sediments corresponding to the onset of the δ13C positive excursion, P content is marked by a sharp peak probably linked to a slowdown in sedimentation rates and/or the presence of a small hiatus, as is shown by the presence of glauconite and phosphatic grains. In the interval corresponding to OAE 2, P values remain low and increase slightly at the end of the positive shift in the δ13C record (in the H. helvetica total range zone). The average contents of RSTE (U, V, As, Co, Mo and Mn) remain low throughout the section and appreciable RSTE enrichments have not been observed for the sedimentary interval corresponding to OAE 2. No correlation is observed with stratigraphic trends in RSTE contents in organic-rich deeper-water sections. The presence of double-keeled planktonic foraminifera species during most of the Cenomanian/Turonian boundary event is another evidence of relatively well-oxygenated conditions in this part of the northern Tethyan outer shelf. Our results show that the Chrummflueschlucht section corresponds to one of the most complete section for the Cenomanian–Turonian boundary interval known from the Helvetic realm even if a small hiatus may be present at the onset of the δ13C record (peak “a”). The evolution of P contents suggests an increase in input of this nutritive element at the onset of OAE2. However, the trends in RSTE contents and the planktonic foraminifera assemblages show that the Helvetic realm has not been affected by strongly depleted oxygen conditions during OAE 2.

The Valanginian δ13C excursion may not be an expression of a global oceanic anoxic event

Marine and terrestrial sediments of the Valanginian age display a distinct positive δ13C excursion, which has recently been interpreted as the expression of an oceanic anoxic episode (OAE) of global importance. Here we evaluate the extent of anaerobic conditions in marine bottom waters and explore the mechanisms involved in changing carbon storage on a global scale during this time interval. We asses redox-sensitive trace-element distributions (RSTE; uranium, vanadium, cobalt, arsenic and molybdenum) and the quality and quantity of preserved organic matter (OM) in representative sections along a shelf-basin transect in the western Tethys, in the Polish Basin and on Shatsky Rise. OM-rich layers corresponding in time to the δ13C shift are generally rare in the Tethyan sections and if present, the layers are not thicker than several centimetres and total organic carbon (TOC) contents do not surpass 1.68 wt..%. Palynological observations and geochemical properties of the preserved OM suggest a mixed marine and terrestrial origin and deposition in an oxic environment. In the Polish Basin, OM-rich layers show evidence for an important continental component. RSTE exhibit no major enrichments during the δ13C excursion in all studied Tethyan sections. RSTE enrichments are, however, observed in the pre-δ13C excursion OM-rich “Barrande” levels of the Vocontian Trough. In addition, all Tethyan sections record higher Mn contents during the δ13C shift, indicating rather well-oxygenated bottom waters in the western Tethys and the presence of anoxic basins elsewhere, such as the restricted basins of the North Atlantic and Weddell Sea. We propose that the Valanginian δ13C shift is the consequence of a combination of increased OM storage in marginal seas and on continents (as a sink of 12C-enriched organic carbon), coupled with the demise of shallow-water carbonate platforms (diminishing the storage capacity of 13C-enriched carbonate carbon). As such the Valanginian provides a more faithful natural analogue to present-day environmental change than most other Mesozoic OAEs, which are characterized by the development of ocean-wide dysaerobic to anaerobic conditions.

Impact of oceanic processes on the carbon cycle during the last termination

During the last termination (from ~18 000 years ago to ~9000 years ago), the climate significantly warmed and the ice sheets melted. Simultaneously, atmospheric CO2 increased from ~190 ppm to ~260 ppm. Although this CO2 rise plays an important role in the deglacial warming, the reasons for its evolution are difficult to explain. Only box models have been used to run transient simulations of this carbon cycle transition, but by forcing the model with data constrained scenarios of the evolution of temperature, sea level, sea ice, NADW formation, Southern Ocean vertical mixing and biological carbon pump. More complex models (including GCMs) have investigated some of these mechanisms but they have only been used to try and explain LGM versus present day steady-state climates. In this study we use a coupled climate-carbon model of intermediate complexity to explore the role of three oceanic processes in transient simulations: the sinking of brines, stratification-dependent diffusion and iron fertilization. Carbonate compensation is accounted for in these simulations. We show that neither iron fertilization nor the sinking of brines alone can account for the evolution of CO2, and that only the combination of the sinking of brines and interactive diffusion can simultaneously simulate the increase in deep Southern Ocean δ13C. The scenario that agrees best with the data takes into account all mechanisms and favours a rapid cessation of the sinking of brines around 18 000 years ago, when the Antarctic ice sheet extent was at its maximum. In this scenario, we make the hypothesis that sea ice formation was then shifted to the open ocean where the salty water is quickly mixed with fresher water, which prevents deep sinking of salty water and therefore breaks down the deep stratification and releases carbon from the abyss. Based on this scenario, it is possible to simulate both the amplitude and timing of the long-term CO2 increase during the last termination in agreement with ice core data. The atmospheric δ13C appears to be highly sensitive to changes in the terrestrial biosphere, underlining the need to better constrain the vegetation evolution during the termination.

The vertical structure of oceanic Rossby waves: a comparison of high-resolution model data to theoretical vertical structures

Tests of the new Rossby wave theories that have been developed over the past decade to account for discrepancies between theoretical wave speeds and those observed by satellite altimeters have focused primarily on the surface signature of such waves. It appears, however, that the surface signature of the waves acts only as a rather weak constraint, and that information on the vertical structure of the waves is required to better discriminate between competing theories. Due to the lack of 3-D observations, this paper uses high-resolution model data to construct realistic vertical structures of Rossby waves and compares these to structures predicted by theory. The meridional velocity of a section at 24° S in the Atlantic Ocean is pre-processed using the Radon transform to select the dominant westward signal. Normalized profiles are then constructed using three complementary methods based respectively on: (1) averaging vertical profiles of velocity, (2) diagnosing the amplitude of the Radon transform of the westward propagating signal at different depths, and (3) EOF analysis. These profiles are compared to profiles calculated using four different Rossby wave theories: standard linear theory (SLT), SLT plus mean flow, SLT plus topographic effects, and theory including mean flow and topographic effects. Our results support the classical theoretical assumption that westward propagating signals have a well-defined vertical modal structure associated with a phase speed independent of depth, in contrast with the conclusions of a recent study using the same model but for different locations in the North Atlantic. The model structures are in general surface intensified, with a sign reversal at depth in some regions, notably occurring at shallower depths in the East Atlantic. SLT provides a good fit to the model structures in the top 300 m, but grossly overestimates the sign reversal at depth. The addition of mean flow slightly improves the latter issue, but is too surface intensified. SLT plus topography rectifies the overestimation of the sign reversal, but overestimates the amplitude of the structure for much of the layer above the sign reversal. Combining the effects of mean flow and topography provided the best fit for the mean model profiles, although small errors at the surface and mid-depths are carried over from the individual effects of mean flow and topography respectively. Across the section the best fitting theory varies between SLT plus topography and topography with mean flow, with, in general, SLT plus topography performing better in the east where the sign reversal is less pronounced. None of the theories could accurately reproduce the deeper sign reversals in the west. All theories performed badly at the boundaries. The generalization of this method to other latitudes, oceans, models and baroclinic modes would provide greater insight into the variability in the ocean, while better observational data would allow verification of the model findings.

An oceanic origin for the increase of atmospheric radiocarbon during the Younger Dryas

Variations in carbon-14 to carbon-12 ratio in the atmosphere (Δ14Catm) provide a powerful diagnostic for elucidating the timing and nature of geophysical and anthropological change. The (Atlantic) marine archive suggests a rapid Δ14Catm increase of 50‰ at the onset of the Younger Dryas (YD) cold reversal (12.9–11.7 kyr BP), which has not yet been satisfactorily explained in terms of magnitude or causal mechanism, as either a change in ocean ventilation or production rate. Using Earth-system model simulations and comparison of marine-based radiocarbon records from different ocean basins, we demonstrate that the YD Δ14Catm increase is smaller than suggested by the marine archive. This is due to changes in reservoir age, predominantly caused by reduced ocean ventilation.

Global monsoons in the mid-Holocene and oceanic feedback

The response of the six major summer monsoon systems (the North American monsoon, the northern Africa monsoon, the Asia monsoon, the northern Australasian monsoon, the South America monsoon and the southern Africa monsoon) to mid-Holocene orbital forcing has been investigated using a coupled ocean–atmosphere general circulation model (FOAM), with the focus on the distinct roles of the direct insolation forcing and oceanic feedback. The simulation result is also found to compare well with the NCAR CSM. The direct effects of the change in insolation produce an enhancement of the Northern Hemisphere monsoons and a reduction of the Southern Hemisphere monsoons. Ocean feedbacks produce a further enhancement of the northern Africa monsoon and the North American monsoon. However, ocean feedbacks appear to weaken the Asia monsoon, although the overall effect (direct insolation forcing plus ocean feedback) remains a strengthened monsoon. The impact of ocean feedbacks on the South American and southern African monsoons is relatively small, and therefore these regions, especially the South America, experienced a reduced monsoon regime compared to present. However, there is a strong ocean feedback on the northern Australian monsoon that negates the direct effects of orbital changes and results in a strengthening of austral summer monsoon precipitation in this region. A new synthesis is made for mid-Holocene paleoenvironmental records and is compared with the model simulations. Overall, model simulations produce changes in regional climates that are generally consistent with paleoenvironmental observations.

Impact of variable atmospheric and oceanic form drag on simulations of Arctic sea ice

Over Arctic sea ice, pressure ridges and floe andmelt pond edges all introduce discrete obstructions to the flow of air or water past the ice and are a source of form drag. In current climate models form drag is only accounted for by tuning the air–ice and ice–ocean drag coefficients, that is, by effectively altering the roughness length in a surface drag parameterization. The existing approach of the skin drag parameter tuning is poorly constrained by observations and fails to describe correctly the physics associated with the air–ice and ocean–ice drag. Here, the authors combine recent theoretical developments to deduce the total neutral form drag coefficients from properties of the ice cover such as ice concentration, vertical extent and area of the ridges, freeboard and floe draft, and the size of floes and melt ponds. The drag coefficients are incorporated into the Los Alamos Sea Ice Model (CICE) and show the influence of the new drag parameterization on the motion and state of the ice cover, with the most noticeable being a depletion of sea ice over the west boundary of the Arctic Ocean and over the Beaufort Sea. The new parameterization allows the drag coefficients to be coupled to the sea ice state and therefore to evolve spatially and temporally. It is found that the range of values predicted for the drag coefficients agree with the range of values measured in several regions of the Arctic. Finally, the implications of the new form drag formulation for the spinup or spindown of the Arctic Ocean are discussed.

The response of the sea ice edge to atmospheric and oceanic jet formation

The sea ice edge presents a region of many feedback processes between the atmosphere, ocean, and sea ice (Maslowski et al.). Here the authors focus on the impact of on-ice atmospheric and oceanic flows at the sea ice edge. Mesoscale jet formation due to the Coriolis effect is well understood over sharp changes in surface roughness such as coastlines (Hunt et al.). This sharp change in surface roughness is experienced by the atmosphere and ocean encountering a compacted sea ice edge. This paper presents a study of a dynamic sea ice edge responding to prescribed atmospheric and oceanic jet formation. An idealized analytical model of sea ice drift is developed and compared to a sea ice climate model [the Los Alamos Sea Ice Model (CICE)] run on an idealized domain. The response of the CICE model to jet formation is tested at various resolutions. It is found that the formation of atmospheric jets at the sea ice edge increases the wind speed parallel to the sea ice edge and results in the formation of a sea ice drift jet in agreement with an observed sea ice drift jet (Johannessen et al.). The increase in ice drift speed is dependent upon the angle between the ice edge and wind and results in up to a 40% increase in ice transport along the sea ice edge. The possibility of oceanic jet formation and the resultant effect upon the sea ice edge is less conclusive. Observations and climate model data of the polar oceans have been analyzed to show areas of likely atmospheric jet formation, with the Fram Strait being of particular interest.

Rapid changes in the redox conditions of the western Tethys Ocean during the early Aptian oceanic anoxic event

The early Aptian (125 to 121 Ma) records an episode of severe environmental change including a major perturbation of the carbon cycle, an oceanic anoxic event (OAE 1a, 122.5 Ma), a platform drowning episode and a biocalcification crisis. We propose to trace changes in the oxygenation state of the ocean during the early Aptian anoxic event using the redox-sensitive trace-element (RSTE) distribution, phosphorus accumulation rates (PARs) and organic-matter characterization in three different basins of the western Tethys. The following sections have been investigated: Gorgo a Cerbara (central Italy) in the Umbria Marche basin, Glaise (SE France) in the Vocontian basin and Cassis/La Bédoule (SE France) located in the Provencal basin. In the Gorgo a Cerbara section, RSTE distributions show a low background level along the main part of the section, contrasted by different maxima in concentrations within the Selli level. In the Glaise section, the Goguel level displays a weak increase in RSTE contents coeval with moderate TOC values. At Cassis/La Bédoule, no significant RSTE enrichments have been observed in sediments equivalent to the Selli level. These differences in the records of the geochemical proxies of the Selli level or its equivalent indicate the deposition under different redox conditions, probably related to the paleogeography. Our data indicate the development of anoxic–euxinic conditions in the deeper part of the Tethys during OAE 1a, whereas in the shallower environments, conditions were less reducing. Moreover, at Gorgo a Cerbara, the Selli level is characterized by rapid changes in the intensity of reducing conditions in the water column. Ocean eutrophication seems to be a major factor in the development and the persistence of anoxia as suggested by the PAR evolution. Higher PAR values at the onset of OAE 1a suggest an increase in nutrient input, whereas the return to lower values through the first part of the OAE 1a interval may be related to the weakened capacity to retain P in the sedimentary reservoir due to bottom-water oxygen depletion. This general pattern is contrasted by the data of Gorgo a Cerbara, where the sediments deposited during the OAE 1a interval show P-enrichments (mainly authigenic P). This is associated with maxima in TOC values and Corg:Ptot ratios, suggesting that a part of the remobilized P was trapped in the sediments and as such prevented from returning to the water column.

The expression of the Cenomanian-Turonian oceanic anoxic event in Tibet

The Gongzha section of Tibet, China is located at the northern margin of the Indian Plate (SE Tethys) and is characterized by hemipelagic grey marls and marly limestones, light grey limestones and silty limestones, but no organic-rich sediments. High-resolution biostratigraphy reveals an expanded Cenomanian–Turonian (CT) boundary interval and the δ13C record includes the main features of the classical positive carbon-isotope excursion that characterizes the CT oceanic anoxic event. The biotic response inferred from the foraminifera suggests that oxic to dysoxic conditions prevailed, except for a short interval marked by peak abundance of Heterohelix that indicates a significantly dysoxic environment during the δ13C “b” peak excursion. The overall decreasing trend in redox-sensitive trace elements (RSTE) during the maximum δ13C excursion confirms the absence of significant longer-lasting anoxia in the Gongzha section. Enrichments in RSTE are linked to phases of increased detrital input. Chemical weathering indices suggest that the upper Cenomanian sediments accumulated under an increasingly hot and humid climate that culminated near the CT boundary. In the early Turonian lower weathering indices suggest a warm, drier climatic regime with reduced continental runoff. Phosphorus mass-accumulation rates show a significant peak at the onset of the positive δ13C excursion, followed by a decrease up to the basal Turonian. This pattern is positively correlated with the long-term decrease in detrital index as also observed in numerous other CT boundary sections (e.g., Eastbourne, Pueblo, and Whadi El Ghaib, Sinaï). Long-term phosphorus accumulation in the Gongzha section is therefore associated with changes in detrital input. The overall decreased detrital input can be explained by the increasingly remote continental sources due to the major transgression at the end of Cenomanian, coupled with changes in continental weathering intensity linked to increasingly more arid climate conditions.

Observational constraints on atmospheric and oceanic cross-equatorial heat transports: revisiting the precipitation asymmetry problem in climate models

Satellite based top-of-atmosphere (TOA) and surface radiation budget observations are combined with mass corrected vertically integrated atmospheric energy divergence and tendency from reanalysis to infer the regional distribution of the TOA, atmospheric and surface energy budget terms over the globe. Hemispheric contrasts in the energy budget terms are used to determine the radiative and combined sensible and latent heat contributions to the cross-equatorial heat transports in the atmosphere (AHT_EQ) and ocean (OHT_EQ). The contrast in net atmospheric radiation implies an AHT_EQ from the northern hemisphere (NH) to the southern hemisphere (SH) (0.75 PW), while the hemispheric difference in sensible and latent heat implies an AHT_EQ in the opposite direction (0.51 PW), resulting in a net NH to SH AHT_EQ (0.24 PW). At the surface, the hemispheric contrast in the radiative component (0.95 PW) dominates, implying a 0.44 PW SH to NH OHT_EQ. Coupled model intercomparison project phase 5 (CMIP5) models with excessive net downward surface radiation and surface-to-atmosphere sensible and latent heat transport in the SH relative to the NH exhibit anomalous northward AHT_EQ and overestimate SH tropical precipitation. The hemispheric bias in net surface radiative flux is due to too much longwave surface radiative cooling in the NH tropics in both clear and all-sky conditions and excessive shortwave surface radiation in the SH subtropics and extratropics due to an underestimation in reflection by clouds.

Atmospheric and oceanic contributions to irreducible forecast uncertainty of Arctic surface climate

Uncertainty of Arctic seasonal to interannual predictions arising from model errors and initial state uncertainty has been widely discussed in the literature, whereas the irreducible forecast uncertainty (IFU) arising from the chaoticity of the climate system has received less attention. However, IFU provides important insights into the mechanisms through which predictability is lost, and hence can inform prioritization of model development and observations deployment. Here, we characterize how internal oceanic and surface atmospheric heat fluxes contribute to IFU of Arctic sea ice and upper ocean heat content in an Earth system model by analyzing a set of idealized ensemble prediction experiments. We find that atmospheric and oceanic heat flux are often equally important for driving unpredictable Arctic-wide changes in sea ice and surface water temperatures, and hence contribute equally to IFU. Atmospheric surface heat flux tends to dominate Arctic-wide changes for lead times of up to a year, whereas oceanic heat flux tends to dominate regionally and on interannual time scales. There is in general a strong negative covariance between surface heat flux and ocean vertical heat flux at depth, and anomalies of lateral ocean heat transport are wind-driven, which suggests that the unpredictable oceanic heat flux variability is mainly forced by the atmosphere. These results are qualitatively robust across different initial states, but substantial variations in the amplitude of IFU exist. We conclude that both atmospheric variability and the initial state of the upper ocean are key ingredients for predictions of Arctic surface climate on seasonal to interannual time scales.

The impact of oceanic heat transport on the atmospheric circulation

A general circulation model of intermediate complexity with an idealized Earth-like aquaplanet setup is used to study the impact of changes in the oceanic heat transport on the global atmospheric circulation. Focus is on the atmospheric mean meridional circulation and global thermodynamic properties. The atmosphere counterbalances to a large extent the imposed changes in the oceanic heat transport, but, nonetheless, significant modifications to the atmospheric general circulation are found. Increasing the strength of the oceanic heat transport up to 2.5 PW leads to an increase in the global mean near-surface temperature and to a decrease in its equator-to-pole gradient. For stronger transports, the gradient is reduced further, but the global mean remains approximately constant. This is linked to a cooling and a reversal of the temperature gradient in the tropics. Additionally, a stronger oceanic heat transport leads to a decline in the intensity and a poleward shift of the maxima of both the Hadley and Ferrel cells. Changes in zonal mean diabatic heating and friction impact the properties of the Hadley cell, while the behavior of the Ferrel cell is mostly controlled by friction. The efficiency of the climate machine, the intensity of the Lorenz energy cycle and the material entropy production of the system decline with increased oceanic heat transport. This suggests that the climate system becomes less efficient and turns into a state of reduced entropy production as the enhanced oceanic transport performs a stronger large-scale mixing between geophysical fluids with different temperatures, thus reducing the available energy in the climate system and bringing it closer to a state of thermal equilibrium.