The impact of westerly wind bursts and ocean initial state on the development, and diversity of El Niño events

Westerly wind bursts (WWBs) that occur in the western tropical Pacific are believed to play an important role in the development of El Niño events. Here, following the study of Lengaigne et al. (Clim Dyn 23(6):601–620, 2004), we conduct numerical simulations in which we reexamine the response of the climate system to an observed wind burst added to a coupled general circulation model. Two sets of twin ensemble experiments are conducted (each set has control and perturbed experiments). In the first set, the initial ocean heat content of the system is higher than the model climatology (recharged), while in the second set it is nearly normal (neutral). For the recharged state, in the absence of WWBs, a moderate El Niño with a maximum warming in the central Pacific (CP) develops in about a year. In contrast, for the neutral state, there develops a weak La Niña. However, when the WWB is imposed, the situation dramatically changes: the recharged state slides into an El Niño with a maximum warming in the eastern Pacific, while the neutral set produces a weak CP El Niño instead of previous La Niña conditions. The different response of the system to the exact same perturbations is controlled by the initial state of the ocean and the subsequent ocean–atmosphere interactions involving the interplay between the eastward shift of the warm pool and the warming of the eastern equatorial Pacific. Consequently, the observed diversity of El Niño, including the occurrence of extreme events, may depend on stochastic atmospheric processes, modulating El Niño properties within a broad continuum.

Decadal predictability of the Atlantic meridional overturning circulation and climate in the IPSL-CM5A-LR model

This study explores the decadal potential predictability of the Atlantic Meridional Overturning Circulation (AMOC) as represented in the IPSL-CM5A-LR model, along with the predictability of associated oceanic and atmospheric fields. Using a 1000-year control run, we analyze the prognostic potential predictability (PPP) of the AMOC through ensembles of simulations with perturbed initial conditions. Based on a measure of the ensemble spread, the modelled AMOC has an average predictive skill of 8 years, with some degree of dependence on the AMOC initial state. Diagnostic potential predictability of surface temperature and precipitation is also identified in the control run and compared to the PPP. Both approaches clearly bring out the same regions exhibiting the highest predictive skill. Generally, surface temperature has the highest skill up to 2 decades in the far North Atlantic ocean. There are also weak signals over a few oceanic areas in the tropics and subtropics. Predictability over land is restricted to the coastal areas bordering oceanic predictable regions. Potential predictability at interannual and longer timescales is largely absent for precipitation in spite of weak signals identified mainly in the Nordic Seas. Regions of weak signals show some dependence on AMOC initial state. All the identified regions are closely linked to decadal AMOC fluctuations suggesting that the potential predictability of climate arises from the mechanisms controlling these fluctuations. Evidence for dependence on AMOC initial state also suggests that studying skills from case studies may prove more useful to understand predictability mechanisms than computing average skill from numerous start dates.

Intercomparison and validation of the mixed layer depth fields of global ocean syntheses

Intercomparison and evaluation of the global ocean surface mixed layer depth (MLD) fields estimated from a suite of major ocean syntheses are conducted. Compared with the reference MLDs calculated from individual profiles, MLDs calculated from monthly mean and gridded profiles show negative biases of 10–20 m in early spring related to the re-stratification process of relatively deep mixed layers. Vertical resolution of profiles also influences the MLD estimation. MLDs are underestimated by approximately 5–7 (14–16) m with the vertical resolution of 25 (50) m when the criterion of potential density exceeding the 10-m value by 0.03 kg m−3 is used for the MLD estimation. Using the larger criterion (0.125 kg m−3) generally reduces the underestimations. In addition, positive biases greater than 100 m are found in wintertime subpolar regions when MLD criteria based on temperature are used. Biases of the reanalyses are due to both model errors and errors related to differences between the assimilation methods. The result shows that these errors are partially cancelled out through the ensemble averaging. Moreover, the bias in the ensemble mean field of the reanalyses is smaller than in the observation-only analyses. This is largely attributed to comparably higher resolutions of the reanalyses. The robust reproduction of both the seasonal cycle and interannual variability by the ensemble mean of the reanalyses indicates a great potential of the ensemble mean MLD field for investigating and monitoring upper ocean processes.

Southern Ocean isopycnal mixing and ventilation changes driven by winds

Observed and predicted changes in the strength of the westerly winds blowing over the Southern Ocean have motivated a number of studies of the response of the Antarctic Circumpolar Current and Southern Ocean Meridional Overturning Circulation (MOC) to wind perturbations and led to the discovery of the``eddy-compensation" regime, wherein the MOC becomes insensitive to wind changes. In addition to the MOC, tracer transport also depends on mixing processes. Here we show, in a high-resolution process model, that isopycnal mixing by mesoscale eddies is strongly dependent on the wind strength. This dependence can be explained by mixing-length theory and is driven by increases in eddy kinetic energy; the mixing length does not change strongly in our simulation. Simulation of a passive ventilation tracer (analogous to CFCs or anthropogenic CO$_2$) demonstrates that variations in tracer uptake across experiments are dominated by changes in isopycnal mixing, rather than changes in the MOC. We argue that, to properly understand tracer uptake under different wind-forcing scenarios, the sensitivity of isopycnal mixing to winds must be accounted for.

Interactions between wind-blown snow redistribution and melt ponds in a coupled ocean–sea ice model

Introducing a parameterization of the interactions between wind-driven snow depth changes and melt pond evolution allows us to improve large scale models. In this paper we have implemented an explicit melt pond scheme and, for the first time, a wind dependant snow redistribution model and new snow thermophysics into a coupled ocean–sea ice model. The comparison of long-term mean statistics of melt pond fractions against observations demonstrates realistic melt pond cover on average over Arctic sea ice, but a clear underestimation of the pond coverage on the multi-year ice (MYI) of the western Arctic Ocean. The latter shortcoming originates from the concealing effect of persistent snow on forming ponds, impeding their growth. Analyzing a second simulation with intensified snow drift enables the identification of two distinct modes of sensitivity in the melt pond formation process. First, the larger proportion of wind-transported snow that is lost in leads directly curtails the late spring snow volume on sea ice and facilitates the early development of melt ponds on MYI. In contrast, a combination of higher air temperatures and thinner snow prior to the onset of melting sometimes make the snow cover switch to a regime where it melts entirely and rapidly. In the latter situation, seemingly more frequent on first-year ice (FYI), a smaller snow volume directly relates to a reduced melt pond cover. Notwithstanding, changes in snow and water accumulation on seasonal sea ice is naturally limited, which lessens the impacts of wind-blown snow redistribution on FYI, as compared to those on MYI. At the basin scale, the overall increased melt pond cover results in decreased ice volume via the ice-albedo feedback in summer, which is experienced almost exclusively by MYI.

Impact of instantaneous sea ice removal in a coupled general circulation model

The impact of extreme sea ice initial conditions on modelled climate is analysed for a fully coupled atmosphere ocean sea ice general circulation model, the Hadley Centre climate model HadCM3. A control run is chosen as reference experiment with greenhouse gas concentration fixed at preindustrial conditions. Sensitivity experiments show an almost complete recovery from total removal or strong increase of sea ice after four years. Thus, uncertainties in initial sea ice conditions seem to be unimportant for climate modelling on decadal or longer time scales. When the initial conditions of the ocean mixed layer were adjusted to ice-free conditions, a few substantial differences remained for more than 15 model years. But these differences are clearly smaller than the uncertainty of the HadCM3 run and all the other 19 IPCC fourth assessment report climate model preindustrial runs. It is an important task to improve climate models in simulating the past sea ice variability to enable them to make reliable projections for the 21st century.

The BALTIMOS (BALTEX Integrated Model System) fi eld experiments: A comprehensive atmospheric boundary layer data set for model validation over the open and ice- covered Baltic Sea

A comprehensive atmospheric boundary layer (ABL) data set was collected in eight fi eld experiments (two during each season) over open water and sea ice in the Baltic Sea during 1998–2001 with the primary objective to validate the coupled atmospheric- ice-ocean-land surface model BALTIMOS (BALTEX Integrated Model System). Measurements were taken by aircraft, ships and surface stations and cover the mean and turbulent structure of the ABL including turbulent fl uxes, radiation fl uxes, and cloud conditions. Measurement examples of the spatial variability of the ABL over the ice edge zone and of the stable ABL over open water demonstrate the wide range of ABL conditions collected and the strength of the data set which can also be used to validate other regional models.

Reconstructing paleosalinity from δ18O: coupled model simulations of the Last Glacial Maximum, Last Interglacial and Late Holocene

Reconstructions of salinity are used to diagnose changes in the hydrological cycle and ocean circulation. A widely used method of determining past salinity uses oxygen isotope (δOw) residuals after the extraction of the global ice volume and temperature components. This method relies on a constant relationship between δOw and salinity throughout time. Here we use the isotope-enabled fully coupled General Circulation Model (GCM) HadCM3 to test the application of spatially and time-independent relationships in the reconstruction of past ocean salinity. Simulations of the Late Holocene (LH), Last Glacial Maximum (LGM), and Last Interglacial (LIG) climates are performed and benchmarked against existing compilations of stable oxygen isotopes in carbonates (δOc), which primarily reflect δOw and temperature. We find that HadCM3 produces an accurate representation of the surface ocean δOc distribution for the LH and LGM. Our simulations show considerable variability in spatial and temporal δOw-salinity relationships. Spatial gradients are generally shallower but within ∼50% of the actual simulated LH to LGM and LH to LIG temporal gradients and temporal gradients calculated from multi-decadal variability are generally shallower than both spatial and actual simulated gradients. The largest sources of uncertainty in salinity reconstructions are found to be caused by changes in regional freshwater budgets, ocean circulation, and sea ice regimes. These can cause errors in salinity estimates exceeding 4 psu. Our results suggest that paleosalinity reconstructions in the South Atlantic, Indian and Tropical Pacific Oceans should be most robust, since these regions exhibit relatively constant δOw-salinity relationships across spatial and temporal scales. Largest uncertainties will affect North Atlantic and high latitude paleosalinity reconstructions. Finally, the results show that it is difficult to generate reliable salinity estimates for regions of dynamic oceanography, such as the North Atlantic, without additional constraints.

Glacial lake drainage in Patagonia (13-8 kyr) and response of the adjacent Pacific Ocean

Large freshwater lakes formed in North America and Europe during deglaciation following the Last Glacial Maximum. Rapid drainage of these lakes into the Oceans resulted in abrupt perturbations in climate, including the Younger Dryas and 8.2 kyr cooling events. In the mid-latitudes of the Southern Hemisphere major glacial lakes also formed and drained during deglaciation but little is known about the magnitude, organization and timing of these drainage events and their e ect on regional climate. We use 16 new single-grain optically stimulated luminescence (OSL) dates to de ne three stages of rapid glacial lake drainage in the Lago General Carrera/Lago Buenos Aires and Lago Cohrane/ Pueyrredón basins of Patagonia and provide the rst assessment of the e ects of lake drainage on the Paci c Ocean. Lake drainage occurred between 13 and 8 kyr ago and was initially gradual eastward into the Atlantic, then subsequently reorganized westward into the Paci c as new drainage routes opened up during Patagonian Ice Sheet deglaciation. Coupled ocean-atmosphere model experiments using HadCM3 with an imposed freshwater surface “hosing” to simulate glacial lake drainage suggest that a negative salinity anomaly was advected south around Cape Horn, resulting in brief but signi cant impacts on coastal ocean vertical mixing and regional climate.

Evaluation of moisture sources for the Central European summer flood of May/June 2013 based on regional climate model simulations

Heavy precipitation affected Central Europe in May/June 2013, triggering damaging floods both on the Danube and the Elbe rivers. Based on a modelling approach with COSMO-CLM, moisture fluxes, backward trajectories, cyclone tracks and precipitation fields are evaluated for the relevant time period 30 May–2 June 2013. We identify potential moisture sources and quantify their contribution to the flood event focusing on the Danube basin through sensitivity experiments: Control simulations are performed with undisturbed ERA-Interim boundary conditions, while multiple sensitivity experiments are driven with modified evaporation characteristics over selected marine and land areas. Two relevant cyclones are identified both in reanalysis and in our simulations, which moved counter-clockwise in a retrograde path from Southeastern Europe over Eastern Europe towards the northern slopes of the Alps. The control simulations represent the synoptic evolution of the event reasonably well. The evolution of the precipitation event in the control simulations shows some differences in terms of its spatial and temporal characteristics compared to observations. The main precipitation event can be separated into two phases concerning the moisture sources. Our modelling results provide evidence that the two main sources contributing to the event were the continental evapotranspiration (moisture recycling; both phases) and the North Atlantic Ocean (first phase only). The Mediterranean Sea played only a minor role as a moisture source. This study confirms the importance of continental moisture recycling for heavy precipitation events over Central Europe during the summer half year.

Seasonality of submesoscale flows in the ocean surface boundary layer

A signature of submesoscale flows in the upper ocean is skewness in the distribution of relative vorticity. Expected to result for high Rossby-number flows, such skewness has implications for mixing, dissipation and stratification within the upper ocean. An array of moorings deployed in the Northeast Atlantic for one year as part of the OSMOSIS experiment reveals that relative vorticity is positively skewed during winter even though the scale of the Rossby number is less than 0.5. Furthermore, this skewness is reduced to zero during spring and autumn. There is also evidence of modest seasonal variations in the gradient Rossby number. The proposed mechanism by which relative vorticity is skewed is that the ratio of lateral to vertical buoyancy gradients, as summarized by the inverse gradient Richardson number, restricts its range during winter but less so at other times of the year. These results support recent observations and model simulations suggesting the upper ocean is host to a seasonal cycle in submesoscale turbulence.

Simulating decadal variability in the North Atlantic Ocean

Observations and climate models suggest significant decadal variability within the North Atlantic subpolar gyre (NA SPG), though observations are sparse and models disagree on the details of this variability. Therefore, it is important to understand 1) the mechanisms of simulated decadal variability, 2) which parts of simulated variability are more faithful representations of reality, and 3) the implications for climate predictions. Here, we investigate the decadal variability in the NA SPG in the state-of-the-art, high resolution (0.25◦ ocean resolution), climate model ‘HadGEM3’. We find a decadal mode with a period of 17 years that explains 30% of the annual variance in related indices. The mode arises due to the advection of heat content anomalies, and shows asymmetries in the timescale of phase reversal between positive and negative phases. A negative feedback from temperature-driven density anomalies in the Labrador Sea (LS) allows for the phase reversal. The North Atlantic Oscillation (NAO), which exhibits the same periodicity, amplifies the mode. The atmosphere-ocean coupling is stronger during positive rather than negative NAO states, explaining the asymmetry. Within the NA SPG, there is potential predictability arising partly from this mode for up to 5 years. There are important similarities between observed and simulated variability, such as the apparent role for the propagation of heat content anomalies. However, observations suggest interannual LS density anomalies are salinity-driven. Salinity control of density would change the temperature feedback to the south, possibly limiting real-world predictive skill in the southern NA SPG with this model. Finally, to understand the diversity of behaviours, we analyse 42 present-generation climate models. Temperature and salinity biases are found to systematically influence the driver of density variability in the LS. Resolution is a good predictor of the biases. The dependence of variability on the background state has important implications for decadal predictions.

Decadal predictions with the HiGEM high resolution global coupled climate model: description and basic evaluation

This paper describes the development and basic evaluation of decadal predictions produced using the HiGEM coupled climate model. HiGEM is a higher resolution version of the HadGEM1 Met Office Unified Model. The horizontal resolution in HiGEM has been increased to 1.25◦ × 0.83◦ in longitude and latitude for the atmosphere, and 1/3◦ × 1/3◦ globally for the ocean. The HiGEM decadal predictions are initialised using an anomaly assimilation scheme that relaxes anomalies of ocean temperature and salinity to observed anomalies. 10 year hindcasts are produced for 10 start dates (1960, 1965,..., 2000, 2005). To determine the relative contributions to prediction skill from initial conditions and external forcing, the HiGEM decadal predictions are compared to uninitialised HiGEM transient experiments. The HiGEM decadal predictions have substantial skill for predictions of annual mean surface air temperature and 100 m upper ocean temperature. For lead times up to 10 years, anomaly correlations (ACC) over large areas of the North Atlantic Ocean, the Western Pacific Ocean and the Indian Ocean exceed values of 0.6. Initialisation of the HiGEM decadal predictions significantly increases skill over regions of the Atlantic Ocean,the Maritime Continent and regions of the subtropical North and South Pacific Ocean. In particular, HiGEM produces skillful predictions of the North Atlantic subpolar gyre for up to 4 years lead time (with ACC > 0.7), which are significantly larger than the uninitialised HiGEM transient experiments.