A comparison of GOCE and drifter-based estimates of the North Atlantic steady-state surface circulation

Over the last decade, due to the Gravity Recovery And Climate Experiment (GRACE) mission and, more recently, the Gravity and steady state Ocean Circulation Explorer (GOCE) mission, our ability to measure the ocean’s mean dynamic topography (MDT) from space has improved dramatically. Here we use GOCE to measure surface current speeds in the North Atlantic and compare our results with a range of independent estimates that use drifter data to improve small scales. We find that, with filtering, GOCE can recover 70% of the Gulf Steam strength relative to the best drifter-based estimates. In the subpolar gyre the boundary currents obtained from GOCE are close to the drifter-based estimates. Crucial to this result is careful filtering which is required to remove small-scale errors, or noise, in the computed surface. We show that our heuristic noise metric, used to determine the degree of filtering, compares well with the quadratic sum of mean sea surface and formal geoid errors obtained from the error variance–covariance matrix associated with the GOCE gravity model. At a resolution of 100 km the North Atlantic mean GOCE MDT error before filtering is 5 cm with almost all of this coming from the GOCE gravity model.

A Lagrangian analysis of ice-supersaturated air over the North Atlantic

Understanding the nature of air parcels that exhibit ice-supersaturation is important because they are the regions of potential formation of both cirrus and aircraft contrails, which affect the radiation balance. Ice-supersaturated air parcels in the upper troposphere and lower stratosphere over the North Atlantic are investigated using Lagrangian trajectories. The trajectory calculations use ERA-Interim data for three winter and three summer seasons, resulting in approximately 200,000 trajectories with ice-supersaturation for each season. For both summer and winter, the median duration of ice-supersaturation along a trajectory is less than 6 hours. 5% of air which becomes ice-supersaturated in the troposphere, and 23% of air which becomes ice-supersaturated in the stratosphere will remain ice-supersaturated for at least 24 hours. Weighting the ice-supersaturation duration with the observed frequency indicates the likely overall importance of the longer duration ice-supersaturated trajectories. Ice-supersaturated air parcels typically experience a decrease in moisture content while ice-supersaturated, suggesting that cirrus clouds eventually form in the majority of such air. A comparison is made between short-lived (less than 24 h) and long-lived (greater than 24 h) ice-supersaturated air flows. For both air flows, ice-supersaturation occurs around the northernmost part of the trajectory. Short-lived ice-supersaturated air flows show no significant differences in speed or direction of movement to subsaturated air parcels. However, long-lived ice-supersaturated air occurs in slower moving air flows, which implies that they are not associated with the fastest moving air through a jet stream.

Predictability of the western North Pacific summer climate associated with different ENSO phases by ENSEMBLES multi-model seasonal forecasts

Predictability of the western North Pacific (WNP) summer climate associated with different El Niño–Southern Oscillation (ENSO) phases is investigated in this study based on the 1-month lead retrospective forecasts of five state-of-the-art coupled models from ENSEMBLES. During the period from 1960 to 2005, the models well capture the WNP summer climate anomalies during most of years in different ENSO phases except the La Niña decaying summers. In the El Niño developing, El Niño decaying and La Niña developing summers, the prediction skills are high for the WNP summer monsoon index (WNPMI), with the prediction correlation larger than 0.7. The high prediction skills of the lower-tropospheric circulation during these phases are found mainly over the tropical western Pacific Ocean, South China Sea and subtropical WNP. These good predictions correspond well to their close teleconnection with ENSO and the high prediction skills of tropical SSTs. By contrast, for the La Niña decaying summers, the prediction skills are considerably low with the prediction correlation for the WNPMI near to zero and low prediction skills around the Philippines and subtropical WNP. These poor predictions relate to the weak summer anomalies of the WNPMI during the La Niña decaying years and no significant connections between the WNP lower-tropospheric circulation anomalies and the SSTs over the tropical central and eastern Pacific Ocean in observations. However, the models tend to predict an apparent anomalous cyclone over the WNP during the La Niña decaying years, indicating a linearity of the circulation response over WNP in the models prediction in comparison with that during the El Niño decaying years which differs from observations. In addition, the models show considerable capability in describing the WNP summer anomalies during the ENSO neutral summers. These anomalies are related to the positive feedback between the WNP lower-tropospheric circulation and the local SSTs. The models can capture this positive feedback but with some uncertainties from different ensemble members during the ENSO neutral summers.

The drivers of projected North Atlantic sea level change

Sea level change predicted by the CMIP5 atmosphere–ocean general circulation models (AOGCMs) is not spatially homogeneous. In particular, the sea level change in the North Atlantic is usually characterised by a meridional dipole pattern with higher sea level rise north of 40°N and lower to the south. The spread among models is also high in that region. Here we evaluate the role of surface buoyancy fluxes by carrying out simulations with the FAMOUS low-resolution AOGCM forced by surface freshwater and heat flux changes from CO2-forced climate change experiments with CMIP5 AOGCMs, and by a standard idealised surface freshwater flux applied in the North Atlantic. Both kinds of buoyancy flux change lead to the formation of the sea level dipole pattern, although the effect of the heat flux has a greater magnitude, and is the main cause of the spread of results among the CMIP5 models. By using passive tracers in FAMOUS to distinguish between additional and redistributed buoyancy, we show that the enhanced sea level rise north of 40°N is mainly due to the direct steric effect (the reduction of sea water density) caused by adding heat or freshwater locally. The surface buoyancy forcing also causes a weakening of the Atlantic meridional overturning circulation, and the consequent reduction of the northward ocean heat transport imposes a negative tendency on sea level rise, producing the reduced rise south of 40°N. However, unlike previous authors, we find that this indirect effect of buoyancy forcing is generally less important than the direct one, except in a narrow band along the east coast of the US, where it plays a major role and leads to sea level rise, as found by previous authors.

An anatomy of the cooling of the North Atlantic Ocean in the 1960s and 1970s

In the 1960s and early 1970s sea surface temperatures in the North Atlantic Ocean cooled rapidly. There is still considerable uncertainty about the causes of this event, although various mechanisms have been proposed. In this observational study it is demonstrated that the cooling proceeded in several distinct stages. Cool anomalies initially appeared in the mid-1960s in the Nordic Seas and Gulf Stream Extension, before spreading to cover most of the Subpolar Gyre. Subsequently, cool anomalies spread into the tropical North Atlantic before retreating, in the late 1970s, back to the Subpolar Gyre. There is strong evidence that changes in atmospheric circulation, linked to a southward shift of the Atlantic ITCZ, played an important role in the event, particularly in the period 1972-76. Theories for the cooling event must account for its distinctive space-time evolution. Our analysis suggests that the most likely drivers were: 1) The “Great Salinity Anomaly” of the late 1960s; 2) An earlier warming of the subpolar North Atlantic, which may have led to a slow-down in the Atlantic Meridional Overturning Circulation; 3) An increase in anthropogenic sulphur dioxide emissions. Determining the relative importance of these factors is a key area for future work.

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.

Cluster analysis of downscaled and explicitly simulated North Atlantic tropical cyclone tracks

A realistic representation of the North Atlantic tropical cyclone tracks is crucial as it allows, for example, explaining potential changes in US landfalling systems. Here we present a tentative study, which examines the ability of recent climate models to represent North Atlantic tropical cyclone tracks. Tracks from two types of climate models are evaluated: explicit tracks are obtained from tropical cyclones simulated in regional or global climate models with moderate to high horizontal resolution (1° to 0.25°), and downscaled tracks are obtained using a downscaling technique with large-scale environmental fields from a subset of these models. For both configurations, tracks are objectively separated into four groups using a cluster technique, leading to a zonal and a meridional separation of the tracks. The meridional separation largely captures the separation between deep tropical and sub-tropical, hybrid or baroclinic cyclones, while the zonal separation segregates Gulf of Mexico and Cape Verde storms. The properties of the tracks’ seasonality, intensity and power dissipation index in each cluster are documented for both configurations. Our results show that except for the seasonality, the downscaled tracks better capture the observed characteristics of the clusters. We also use three different idealized scenarios to examine the possible future changes of tropical cyclone tracks under 1) warming sea surface temperature, 2) increasing carbon dioxide, and 3) a combination of the two. The response to each scenario is highly variable depending on the simulation considered. Finally, we examine the role of each cluster in these future changes and find no preponderant contribution of any single cluster over the others.

Sources, sinks, and transatlantic transport of North African dust aerosol: a multimodel analysis and comparison with remote sensing data

This study evaluates model-simulated dust aerosols over North Africa and the North Atlantic from five global models that participated in the Aerosol Comparison between Observations and Models phase II model experiments. The model results are compared with satellite aerosol optical depth (AOD) data from Moderate Resolution Imaging Spectroradiometer (MODIS), Multiangle Imaging Spectroradiometer (MISR), and Sea-viewing Wide Field-of-view Sensor, dust optical depth (DOD) derived from MODIS and MISR, AOD and coarse-mode AOD (as a proxy of DOD) from ground-based Aerosol Robotic Network Sun photometer measurements, and dust vertical distributions/centroid height from Cloud Aerosol Lidar with Orthogonal Polarization and Atmospheric Infrared Sounder satellite AOD retrievals. We examine the following quantities of AOD and DOD: (1) the magnitudes over land and over ocean in our study domain, (2) the longitudinal gradient from the dust source region over North Africa to the western North Atlantic, (3) seasonal variations at different locations, and (4) the dust vertical profile shape and the AOD centroid height (altitude above or below which half of the AOD is located). The different satellite data show consistent features in most of these aspects; however, the models display large diversity in all of them, with significant differences among the models and between models and observations. By examining dust emission, removal, and mass extinction efficiency in the five models, we also find remarkable differences among the models that all contribute to the discrepancies of model-simulated dust amount and distribution. This study highlights the challenges in simulating the dust physical and optical processes, even in the best known dust environment, and stresses the need for observable quantities to constrain the model processes.

Improving climate change detection through optimal seasonal averaging: the case of the North Atlantic jet and European precipitation

The detection of anthropogenic climate change can be improved by recognising the seasonality in the climate change response. This is demonstrated for the North Atlantic jet (zonal wind at 850 hPa, U850) and European precipitation responses projected by the CMIP5 climate models. The U850 future response is characterised by a marked seasonality: an eastward extension of the North Atlantic jet into Europe in November-April, and a poleward shift in May-October. Under the RCP8.5 scenario, the multi-model mean response in U850 in these two extended seasonal means emerges by 2035-2040 for the lower--latitude features and by 2050-2070 for the higher--latitude features, relative to the 1960-1990 climate. This is 5-15 years earlier than when evaluated in the traditional meteorological seasons (December--February, June--August), and it results from an increase in the signal to noise ratio associated with the spatial coherence of the response within the extended seasons. The annual mean response lacks important information on the seasonality of the response without improving the signal to noise ratio. The same two extended seasons are demonstrated to capture the seasonality of the European precipitation response to climate change and to anticipate its emergence by 10-20 years. Furthermore, some of the regional responses, such as the Mediterranean precipitation decline and the U850 response in North Africa in the extended winter, are projected to emerge by 2020-2025, according to the models with a strong response. Therefore, observations might soon be useful to test aspects of the atmospheric circulation response predicted by some of the CMIP5 models.

Exploring the impact of CMIP5 model biases on the simulation of North Atlantic decadal variability

Instrumental observations, palaeo-proxies, and climate models suggest significant decadal variability within the North Atlantic subpolar gyre (NASPG). However, a poorly sampled observational record and a diversity of model behaviours mean that the precise nature and mechanisms of this variability are unclear. Here, we analyse an exceptionally large multi-model ensemble of 42 present-generation climate models to test whether NASPG mean state biases systematically affect the representation of decadal variability. Temperature and salinity biases in the Labrador Sea co-vary and influence whether density variability is controlled by temperature or salinity variations. Ocean horizontal resolution is a good predictor of the biases and the location of the dominant dynamical feedbacks within the NASPG. However, we find no link to the spectral characteristics of the variability. Our results suggest that the mean state and mechanisms of variability within the NASPG are not independent. This represents an important caveat for decadal predictions using anomaly-assimilation methods.

Future changes in the Western North Pacific tropical cyclone activity projected by a multidecadal simulation with a 16-km global atmospheric GCM

How tropical cyclone (TC) activity in the northwestern Pacific might change in a future climate is assessed using multidecadal Atmospheric Model Intercomparison Project (AMIP)-style and time-slice simulations with the ECMWF Integrated Forecast System (IFS) at 16-km and 125-km global resolution. Both models reproduce many aspects of the present-day TC climatology and variability well, although the 16-km IFS is far more skillful in simulating the full intensity distribution and genesis locations, including their changes in response to El Niño–Southern Oscillation. Both IFS models project a small change in TC frequency at the end of the twenty-first century related to distinct shifts in genesis locations. In the 16-km IFS, this shift is southward and is likely driven by the southeastward penetration of the monsoon trough/subtropical high circulation system and the southward shift in activity of the synoptic-scale tropical disturbances in response to the strengthening of deep convective activity over the central equatorial Pacific in a future climate. The 16-km IFS also projects about a 50% increase in the power dissipation index, mainly due to significant increases in the frequency of the more intense storms, which is comparable to the natural variability in the model. Based on composite analysis of large samples of supertyphoons, both the development rate and the peak intensities of these storms increase in a future climate, which is consistent with their tendency to develop more to the south, within an environment that is thermodynamically more favorable for faster development and higher intensities. Coherent changes in the vertical structure of supertyphoon composites show system-scale amplification of the primary and secondary circulations with signs of contraction, a deeper warm core, and an upward shift in the outflow layer and the frequency of the most intense updrafts. Considering the large differences in the projections of TC intensity change between the 16-km and 125-km IFS, this study further emphasizes the need for high-resolution modeling in assessing potential changes in TC activity.

Interdecadal changes on the seasonal prediction of the western North Pacific summer climate around the late 1970s and early 1990s

Identifying predictability and the corresponding sources for the western North Pacific (WNP) summer climate in the case of non-stationary teleconnections during recent decades benefits for further improvements of long-range prediction on the WNP and East Asian summers. In the past few decades, pronounced increases on the summer sea surface temperature (SST) and associated interannual variability are observed over the tropical Indian Ocean and eastern Pacific around the late 1970s and over the Maritime Continent and western–central Pacific around the early 1990s. These increases are associated with significant enhancements of the interannual variability for the lower-tropospheric wind over the WNP. In this study, we further assess interdecadal changes on the seasonal prediction of the WNP summer anomalies, using May-start retrospective forecasts from the ENSEMBLES multi-model project in the period 1960–2005. It is found that prediction of the WNP summer anomalies exhibits an interdecadal shift with higher prediction skills since the late 1970s, particularly after the early 1990s. Improvements of the prediction skills for SSTs after the late 1970s are mainly found around tropical Indian Ocean and the WNP. The better prediction of the WNP after the late 1970s may arise mainly from the improvement of the SST prediction around the tropical eastern Indian Ocean. The close teleconnections between the tropical eastern Indian Ocean and WNP summer variability work both in the model predictions and observations. After the early 1990s, on the other hand, the improvements are detected mainly around the South China Sea and Philippines for the lower-tropospheric zonal wind and precipitation anomalies, associating with a better description of the SST anomalies around the Maritime Continent. A dipole SST pattern over the Maritime Continent and the central equatorial Pacific Ocean is closely related to the WNP summer anomalies after the early 1990s. This teleconnection mode is quite predictable, which is realistically reproduced by the models, presenting more predictable signals to the WNP summer climate after the early 1990s.

Impact of the North Atlantic Oscillation on transatlantic flight routes and clear-air turbulence

The variation of wind-optimal transatlantic flight routes and their turbulence potential is investigated to understand how upper-level winds and large-scale flow patterns can affect the efficiency and safety of long-haul flights. In this study, the wind-optimal routes (WORs) that minimize the total flight time by considering wind variations are modeled for flights between John F. Kennedy International Airport (JFK) in New York, New York, and Heathrow Airport (LHR) in London, United Kingdom, during two distinct winter periods of abnormally high and low phases of North Atlantic Oscillation (NAO) teleconnection patterns. Eastbound WORs approximate the JFK–LHR great circle (GC) route following northerly shifted jets in the +NAO period. Those WORs deviate southward following southerly shifted jets during the −NAO period, because eastbound WORs fly closely to the prevailing westerly jets to maximize tailwinds. Westbound WORs, however, spread meridionally to avoid the jets near the GC in the +NAO period to minimize headwinds. In the −NAO period, westbound WORs are north of the GC because of the southerly shifted jets. Consequently, eastbound WORs are faster but have higher probabilities of encountering clear-air turbulence than westbound ones, because eastbound WORs are close to the jet streams, especially near the cyclonic shear side of the jets in the northern (southern) part of the GC in the +NAO (−NAO) period. This study suggests how predicted teleconnection weather patterns can be used for long-haul strategic flight planning, ultimately contributing to minimizing aviation’s impact on the environment

On the ability of statistical wind-wave models to capture the variability and long-term trends of the North Atlantic winter wave climate

A dynamical wind-wave climate simulation covering the North Atlantic Ocean and spanning the whole 21st century under the A1B scenario has been compared with a set of statistical projections using atmospheric variables or large scale climate indices as predictors. As a first step, the performance of all statistical models has been evaluated for the present-day climate; namely they have been compared with a dynamical wind-wave hindcast in terms of winter Significant Wave Height (SWH) trends and variance as well as with altimetry data. For the projections, it has been found that statistical models that use wind speed as independent variable predictor are able to capture a larger fraction of the winter SWH inter-annual variability (68% on average) and of the long term changes projected by the dynamical simulation. Conversely, regression models using climate indices, sea level pressure and/or pressure gradient as predictors, account for a smaller SWH variance (from 2.8% to 33%) and do not reproduce the dynamically projected long term trends over the North Atlantic. Investigating the wind-sea and swell components separately, we have found that the combination of two regression models, one for wind-sea waves and another one for the swell component, can improve significantly the wave field projections obtained from single regression models over the North Atlantic.