Importance of oceanic resolution and mean state on the extra-tropical response to El Niño in a matrix of coupled models

The extra-tropical response to El Niño in configurations of a coupled model with increased horizontal resolution in the oceanic component is shown to be more realistic than in configurations with a low resolution oceanic component. This general conclusion is independent of the atmospheric resolution. Resolving small-scale processes in the ocean produces a more realistic oceanic mean state, with a reduced cold tongue bias, which in turn allows the atmospheric model component to be forced more realistically. A realistic atmospheric basic state is critical in order to represent Rossby wave propagation in response to El Niño, and hence the extra-tropical response to El Niño. Through the use of high and low resolution configurations of the forced atmospheric-only model component we show that, in isolation, atmospheric resolution does not significantly affect the simulation of the extra-tropical response to El Niño. It is demonstrated, through perturbations to the SST forcing of the atmospheric model component, that biases in the climatological SST field typical of coupled model configurations with low oceanic resolution can account for the erroneous atmospheric basic state seen in these coupled model configurations. These results highlight the importance of resolving small-scale oceanic processes in producing a realistic large-scale mean climate in coupled models, and suggest that it might may be possible to “squeeze out” valuable extra performance from coupled models through increases to oceanic resolution alone.

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.

On the connection between the 1984 Atlantic warm event and the 1982-1983 ENSO

The warm event which spread in the tropical Atlantic during Spring-Summer 1984 is assumed to be partially initiated by atmospheric disturbances, themselves related to the major 1982–1983 El-Niño which occurred 1 year earlier in the Pacific. This paper tests such an hypothesis. For that purpose, an atmospheric general circulation model (AGCM) is forced by different conditions of climatic and observed sea surface temperature and an Atlantic ocean general circulation model (OGCM) is subsequently forced by the outputs of the AGCM. It is firstly shown that both the AGCM and the OGCM correctly behave when globally observed SST are used: the strengthening of the trades over the tropical Atlantic during 1983 and their subsequent weakening at the beginning of 1984 are well captured by the AGCM, and so is the Spring 1984 deepening of the thermocline in the eastern equatorial Atlantic, simulated by the OGCM. As assumed, the SST anomalies located in the El-Niño Pacific area are partly responsible for wind signal anomaly in the tropical Atlantic. Though this remotely forced atmospheric signal has a small amplitude, it can generate, in the OGCM run, an anomalous sub-surface signal leading to a flattening of the thermocline in the equatorial Atlantic. This forced oceanic experiment cannot explain the amplitude and phase of the observed sub-surface oceanic anomaly: part of the Atlantic ocean response, due to local interaction between ocean and atmosphere, requires a coupled approach. Nevertheless this experiment showed that anomalous conditions in the Pacific during 82–83 created favorable conditions for anomaly development in the Atlantic.

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.

Interdecadal enhancement of the Walker circulation over the Tropical Pacific in the late 1990s

The Walker circulation is one of the major components of the large-scale tropical atmospheric circulation and variations in its strength are critical to equatorial Pacific Ocean circulation. It has been argued in the literature that during the 20th century the Walker circulation weakened, and that this weakening was attributable to anthropogenic climate change. By using updated observations, we show that there has been a rapid interdecadal enhancement of the Walker circulation since the late 1990s. Associated with this enhancement is enhanced precipitation in the tropical western Pacific, anomalous westerlies in the upper troposphere, descent in the central and eastern tropical Pacific, and anomalous surface easterlies in the western and central tropical Pacific. The characteristics of associated oceanic changes are a strengthened thermocline slope and an enhanced zonal SST gradient across the tropical Pacific. Many characteristics of these changes are similar to those associated with the mid-1970s climate shift with an opposite sign. We also show that the interdecadal variability of the Walker circulation in the tropical Pacific is inversely correlated to the interdecadal variability of the zonal circulation in the tropical Atlantic. An enhancement of the Walker circulation in the tropical Pacific is associated with a weakening zonal circulation in the tropical Atlantic and vise versa, implying an inter-Atlantic-Pacific connection of the zonal overturning circulation variation. Whether these recent changes will be sustained is not yet clear, but our research highlights the importance of understanding the interdecadal variability, as well as the long-term trends, that influence tropical circulation.

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.