A model study of factors influencing projected changes in regional sea level over the 21st century

In addition to projected increases in global mean sea level over the 21st century, model simulations suggest there will also be changes in the regional distribution of sea level relative to the global mean. There is a considerable spread in the projected patterns of these changes by current models, as shown by the recent Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment (AR4). This spread has not reduced from that given by the Third Assessment models. Comparison with projections by ensembles of models based on a single structure supports an earlier suggestion that models of similar formulation give more similar patterns of sea level change. Analysing an AR4 ensemble of model projections under a business-as-usual scenario shows that steric changes (associated with subsurface ocean density changes) largely dominate the sea level pattern changes. The relative importance of subsurface temperature or salinity changes in contributing to this differs from region to region and, to an extent, from model-to-model. In general, thermosteric changes give the spatial variations in the Southern Ocean, halosteric changes dominate in the Arctic and strong compensation between thermosteric and halosteric changes characterises the Atlantic. The magnitude of sea level and component changes in the Atlantic appear to be linked to the amount of Atlantic meridional overturning circulation (MOC) weakening. When the MOC weakening is substantial, the Atlantic thermosteric patterns of change arise from a dominant role of ocean advective heat flux changes.

Processes governing the predictability of the Atlantic meridional overturning circulation in a coupled GCM

The processes that govern the predictability of decadal variations in the North Atlantic meridional overturning circulation (MOC) are investigated in a long control simulation of the ECHO-G coupled atmosphere–ocean model. We elucidate the roles of local stochastic forcing by the atmosphere, and other potential ocean processes, and use our results to build a predictive regression model. The primary influence on MOC variability is found to come from air–sea heat fluxes over the Eastern Labrador Sea. The maximum correlation between such anomalies and the variations in the MOC occurs at a lead time of 2 years, but we demonstrate that the MOC integrates the heat flux variations over a period of 10 years. The corresponding univariate regression model accounts for 74.5% of the interannual variability in the MOC (after the Ekman component has been removed). Dense anomalies to the south of the Greenland-Scotland ridge are also shown to precede the overturning variations by 4–6 years, and provide a second predictor. With the inclusion of this second predictor the resulting regression model explains 82.8% of the total variance of the MOC. This final bivariate model is also tested during large rapid decadal overturning events. The sign of the rapid change is always well represented by the bivariate model, but the magnitude is usually underestimated, suggesting that other processes are also important for these large rapid decadal changes in the MOC.

Mountain ranges favour vigorous Atlantic Meridional Overturning

We use a global Ocean-Atmosphere General Circulation Model (OAGCM) to show that the major mountain ranges of the world have a significant role in maintenance of the Atlantic Meridional Overturning Circulation (AMOC). A simulation with mountains has a maximum AMOC of 18 Sv (1 Sv=106 m3 s-1) compared with ~0 Sv for a simulation without mountains. Atlantic heat transport at 25N is 1.1 PW with mountains compared to 0.2 PW without. The difference in AMOC is due to major changes in surface heat and freshwater (FW) fluxes over the Atlantic. In the Pacific changed surface fluxes lead to a meridional overturning circulation of 10 Sv. Our results suggest that the effects of mountains on the large-scale atmospheric circulation is to force the ocean towards a state with a vigorous AMOC and with no overturning in the Pacific.

Wind-driven mixing below the oceanic mixed layer

This study describes the turbulent processes in the upper ocean boundary layer forced by a constant surface stress in the absence of the Coriolis force using large-eddy simulation. The boundary layer that develops has a two-layer structure, a well-mixed layer above a stratified shear layer. The depth of the mixed layer is approximately constant, whereas the depth of the shear layer increases with time. The turbulent momentum flux varies approximately linearly from the surface to the base of the shear layer. There is a maximum in the production of turbulence through shear at the base of the mixed layer. The magnitude of the shear production increases with time. The increase is mainly a result of the increase in the turbulent momentum flux at the base of the mixed layer due to the increase in the depth of the boundary layer. The length scale for the shear turbulence is the boundary layer depth. A simple scaling is proposed for the magnitude of the shear production that depends on the surface forcing and the average mixed layer current. The scaling can be interpreted in terms of the divergence of a mean kinetic energy flux. A simple bulk model of the boundary layer is developed to obtain equations describing the variation of the mixed layer and boundary layer depths with time. The model shows that the rate at which the boundary layer deepens does not depend on the stratification of the thermocline. The bulk model shows that the variation in the mixed layer depth is small as long as the surface buoyancy flux is small.

Precise calculations of the existence of multiple AMOC equilibria in coupled climate models Part I: equilibrium states

This study examines criteria for the existence of two stable states of the Atlantic Meridional Overturning Circulation (AMOC) using a combination of theory and simulations from a numerical coupled atmosphere–ocean climate model. By formulating a simple collection of state parameters and their relationships, the authors reconstruct the North Atlantic Deep Water (NADW) OFF state behavior under a varying external salt-flux forcing. This part (Part I) of the paper examines the steady-state solution, which gives insight into the mechanisms that sustain the NADW OFF state in this coupled model; Part II deals with the transient behavior predicted by the evolution equation. The nonlinear behavior of the Antarctic Intermediate Water (AAIW) reverse cell is critical to the OFF state. Higher Atlantic salinity leads both to a reduced AAIW reverse cell and to a greater vertical salinity gradient in the South Atlantic. The former tends to reduce Atlantic salt export to the Southern Ocean, while the latter tends to increases it. These competing effects produce a nonlinear response of Atlantic salinity and salt export to salt forcing, and the existence of maxima in these quantities. Thus the authors obtain a natural and accurate analytical saddle-node condition for the maximal surface salt flux for which a NADW OFF state exists. By contrast, the bistability indicator proposed by De Vries and Weber does not generally work in this model. It is applicable only when the effect of the AAIW reverse cell on the Atlantic salt budget is weak.

The link between the Barents Sea and ENSO events simulated by NEMO model

An analysis of observational data in the Barents Sea along a meridian at 33°30' E between 70°30' and 72°30' N has reported a negative correlation between El Niño/La Niña Southern Oscillation (ENSO) events and water temperature in the top 200 m: the temperature drops about 0.5 °C during warm ENSO events while during cold ENSO events the top 200 m layer of the Barents Sea is warmer. Results from 1 and 1/4-degree global NEMO models show a similar response for the whole Barents Sea. During the strong warm ENSO event in 1997–1998 an anomalous anticyclonic atmospheric circulation over the Barents Sea enhances heat loses, as well as substantially influencing the Barents Sea inflow from the North Atlantic, via changes in ocean currents. Under normal conditions along the Scandinavian peninsula there is a warm current entering the Barents Sea from the North Atlantic, however after the 1997–1998 event this current is weakened. During 1997–1998 the model annual mean temperature in the Barents Sea is decreased by about 0.8 °C, also resulting in a higher sea ice volume. In contrast during the cold ENSO events in 1999–2000 and 2007–2008, the model shows a lower sea ice volume, and higher annual mean temperatures in the upper layer of the Barents Sea of about 0.7 °C. An analysis of model data shows that the strength of the Atlantic inflow in the Barents Sea is the main cause of heat content variability, and is forced by changing pressure and winds in the North Atlantic. However, surface heat-exchange with the atmosphere provides the means by which the Barents sea heat budget relaxes to normal in the subsequent year after the ENSO events.

The changing energy balance of the polar regions in a warmer climate

Energy fluxes for polar regions are examined for two 30-year periods, representing the end of the 20th and 21st centuries, using data from high resolution simulations with the ECHAM5 climate model. The net radiation to space for the present climate agrees well with data from the Clouds and the Earth’s Radiant Energy System (CERES) over the northern polar region but shows an underestimation in planetary albedo for the southern polar region. This suggests there are systematic errors in the atmospheric circulation or in the net surface energy fluxes in the southern polar region. The simulation of the future climate is based on the Intergovernmental Panel on Climate Change (IPCC) A1B scenario. The total energy transport is broadly the same for the two 30 year periods, but there is an increase in the moist energy transport of the order of 6 W m−2 and a corresponding reduction in the dry static energy. For the southern polar region the proportion of moist energy transport is larger and the dry static energy correspondingly smaller for both periods. The results suggest a possible mechanism for the warming of the Arctic that is discussed. Changes between the 20th and 21st centuries in the northern polar region show the net ocean surface radiation flux in summer increases ~18W m−2 (24%). For the southern polar region the response is different as there is a decrease in surface solar radiation. We suggest that this is caused by changes in cloudiness associated with the poleward migration of the storm tracks.

Tropical Atlantic variability modes (1979–2002). Part I: Time-evolving SST modes related to West African rainfall

This work presents a description of the 1979–2002 tropical Atlantic (TA) SST variability modes coupled to the anomalous West African (WA) rainfall during the monsoon season. The time-evolving SST patterns, with an impact on WA rainfall variability, are analyzed using a new methodology based on maximum covariance analysis. The enhanced Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP) dataset, which includes measures over the ocean, gives a complete picture of the interannual WA rainfall patterns for the Sahel dry period. The leading TA SST pattern, related to the Atlantic El Niño, is coupled to anomalous precipitation over the coast of the Gulf of Guinea, which corresponds to the second WA rainfall principal component. The thermodynamics and dynamics involved in the generation, development, and damping of this mode are studied and compared with previous works. The SST mode starts at the Angola/Benguela region and is caused by alongshore wind anomalies. It then propagates westward via Rossby waves and damps because of latent heat flux anomalies and Kelvin wave eastward propagation from an off-equatorial forcing. The second SST mode includes the Mediterranean and the Atlantic Ocean, showing how the Mediterranean SST anomalies are those that are directly associated with the Sahelian rainfall. The global signature of the TA SST patterns is analyzed, adding new insights about the Pacific– Atlantic link in relation to WA rainfall during this period. Also, this global picture suggests that the Mediterranean SST anomalies are a fingerprint of large-scale forcing. This work updates the results given by other authors, whose studies are based on different datasets dating back to the 1950s, including both the wet and the dry Sahel periods.

Changes in tropical Atlantic interannual variability from a substantial weakening of the meridional overturning circulation

In response to a substantial weakening of the Atlantic Meridional Overturning Circulation (AMOC)— from a coupled ocean–atmosphere general circulation model experiment—significant changes in the interannual variability are found over the tropical Atlantic, characterized by an increase of variance (by ~150 %) in boreal late spring-early summer and a decrease of variance (by ~60 %) in boreal autumn. This study focuses on understanding physical mechanisms responsible for these changes in interannual variability in the tropical Atlantic. It demonstrates that the increase of variability in spring is a consequence of an increase in the variance of the El Niño-Southern Oscillation, which has a large impact on the tropical Atlantic via anomalous surface heat fluxes. Winter El Niño (La Niña) affects the eastern equatorial Atlantic by decreasing (increasing) cloud cover and surface wind speed which is associated with anomalous downward (upward) short wave radiation and reduced (enhanced) upward latent heat fluxes, creating anomalous positive (negative) sea surface temperature (SST) anomalies over the region from winter to spring. On the other hand, the decrease of SST variance in autumn is due to a deeper mean thermocline which weakens the impact of the thermocline movement on SST variation. The comparison between the model results and observations is not straightforward owing to the influence of model biases and the lack of a major MOC weakening event in the instrumental record. However, it is argued that the basic physical mechanisms found in the model simulations are likely to be robust and therefore have relevance to understanding tropical Atlantic variability in the real world, perhaps with modified seasonality.

The transient atmospheric response to midlatitude SST anomalies

To study the transient atmospheric response to midlatitude SST anomalies, a three-layer quasigeostrophic (QG) model coupled to a slab oceanic mixed layer in the North Atlantic is used. As diagnosed from a coupled run in perpetual winter conditions, the first two modes of SST variability are linked to the model North Atlantic Oscillation (NAO) and eastern Atlantic pattern (EAP), respectively, the dominant atmospheric modes in the Atlantic sector. The two SST anomaly patterns are then prescribed as fixed anomalous boundary conditions for the model atmosphere, and its transient responses are established from a large ensemble of simulations. In both cases, the tendency of the air–sea heat fluxes to damp the SST anomalies results in an anomalous diabatic heating of the atmosphere that, in turn, forces a baroclinic response, as predicted by linear theory. This initial response rapidly modifies the transient eddy activity and thus the convergence of eddy momentum and heat fluxes. The latter transforms the baroclinic response into a growing barotropic one that resembles the atmospheric mode that had created the SST anomaly in the coupled run and is thus associated with a positive feedback. The total adjustment time is as long as 3–4 months for the NAO-like response and 1–2 months for the EAP-like one. The positive feedback, in both cases, is dependent on the polarity of the SST anomaly, but is stronger in the NAO case, thereby contributing to its predominance at low frequency in the coupled system. However, the feedback is too weak to lead to an instability of the atmospheric modes and primarily results in an increase of their amplitude and persistence and a weakening of the heat flux damping of the SST anomaly.

Eddy-driven exchange between the open ocean and a sub-ice shelf cavity

The exchange between the open ocean and sub-ice shelf cavities is important to both water mass transformations and ice shelf melting. Here we use a high-resolution (500 m) numerical model to investigate to which degree eddies produced by frontal instability at the edge of a polynya are capable of transporting dense High Salinity Shelf Water (HSSW) underneath an ice shelf. The applied surface buoyancy flux and ice shelf geometry is based on Ronne Ice Shelf in the southern Weddell Sea, an area of intense wintertime sea ice production where a flow of HSSW into the cavity has been observed. Results show that eddies are able to enter the cavity at the southwestern corner of the polynya where an anticyclonic rim current intersects the ice shelf front. The size and time scale of simulated eddies are in agreement with observations close to the Ronne Ice Front. The properties and strength of the inflow are sensitive to the prescribed total ice production, flushing the ice shelf cavity at a rate of 0.2–0.4 × 106 m3 s−1 depending on polynya size and magnitude of surface buoyancy flux. Eddy-driven HSSW transport into the cavity is reduced by about 50% if the model grid resolution is decreased to 2-5 km and eddies are not properly resolved.

A nonlinear oscillator describing storm track variability

We construct a two-variable model which describes the interaction between local baroclinicity and eddy heat flux in order to understand aspects of the variance in storm tracks. It is a heuristic model for diabatically forced baroclinic instability close to baroclinic neutrality. The two-variable model has the structure of a nonlinear oscillator. It exhibits some realistic properties of observed storm track variability, most notably the intermittent nature of eddy activity. This suggests that apparent threshold behaviour can be more accurately and succinctly described by a simple nonlinearity. An analogy is drawn with triggering of convective events.

Transient atmospheric response to interactive SST anomalies

The transient atmospheric response to interactive SST anomalies in the midlatitudes is investigated using a three-layer QG model coupled in perpetual winter conditions to a slab oceanic mixed layer in the North Atlantic. The SST anomalies are diagnosed from a coupled run and prescribed as initial conditions, but are free to evolve. The initial evolution of the atmospheric response is similar to that obtained with a prescribed SST anomaly, starting as a quasi-linear baroclinic and then quickly evolving into a growing equivalent barotropic one. Because of the heat flux damping, the SST anomaly amplitude slowly decreases, albeit with little change in pattern. Correspondingly, the atmospheric response only increases until it reaches a maximum amplitude after about 1–3.5 months, depending on the SST anomaly considered. The response is similar to that at equilibrium in the fixed SST case, but it is 1.5–2 times smaller, and then slowly decays away.

The lifecycle of the North Atlantic storm track

The North Atlantic eddy-driven jet exhibits latitudinal variability, with evidence of three preferred latitudinal locations: south, middle and north. Here we examine the drivers of this variability and the variability of the associated storm track. We investigate the changes in the storm track characteristics for the three jet locations, and propose a mechanism by which enhanced storm track activity, as measured by upstream heat flux, is responsible for cyclical downstream latitudinal shifts in the jet. This mechanism is based on a nonlinear oscillator relationship between the enhanced meridional temperature gradient (and thus baroclinicity) and the meridional high-frequency (periods of shorter than 10 days) eddy heat flux. Such oscillations in baroclinicity and heat flux induce variability in eddy anisotropy which is associated with the changes in the dominant type of wave breaking and a different latitudinal deflection of the jet. Our results suggest that high heat flux is conducive to a northward deflection of the jet, whereas low heat flux is conducive to a more zonal jet. This jet deflecting effect was found to operate most prominently downstream of the storm track maximum, while the storm track and the jet remain anchored at a fixed latitudinal location at the beginning of the storm track. These cyclical changes in storm track characteristics can be viewed as different stages of the storm track’s spatio-temporal lifecycle.

Structure of the planetary boundary layer over southeast England: modeling and measurements

The Weather Research and Forecasting model was applied to analyze variations in the planetary boundary layer (PBL) structure over Southeast England including central and suburban London. The parameterizations and predictive skills of two nonlocal mixing PBL schemes, YSU and ACM2, and two local mixing PBL schemes, MYJ and MYNN2, were evaluated over a variety of stability conditions, with model predictions at a 3 km grid spacing. The PBL height predictions, which are critical for scaling turbulence and diffusion in meteorological and air quality models, show significant intra-scheme variance (> 20%), and the reasons are presented. ACM2 diagnoses the PBL height thermodynamically using the bulk Richardson number method, which leads to a good agreement with the lidar data for both unstable and stable conditions. The modeled vertical profiles in the PBL, such as wind speed, turbulent kinetic energy (TKE), and heat flux, exhibit large spreads across the PBL schemes. The TKE predicted by MYJ were found to be too small and show much less diurnal variation as compared with observations over London. MYNN2 produces better TKE predictions at low levels than MYJ, but its turbulent length scale increases with height in the upper part of the strongly convective PBL, where it should decrease. The local PBL schemes considerably underestimate the entrainment heat fluxes for convective cases. The nonlocal PBL schemes exhibit stronger mixing in the mean wind fields under convective conditions than the local PBL schemes and agree better with large-eddy simulation (LES) studies.