The steady‐state structure of the natural stratosphere and ozone distribution in a 2‐D model incorporating radiation and O‐H‐N photochemistry and the effects of stratospheric pollutants

A mean meridional circulation model of the stratosphere, incorporating radiative heating and photochemistry of the oxygen‐hydrogen‐nitrogen atmosphere, is used to simulate the meridional distributions of O₃, HO_X, N₂O,NO_X, temperature and the three components of mean motion for the summer and winter seasons under steady‐state conditions. The results are generally in good agreement with the available observations in the normal stratosphere. The model has been applied to assess the effects of water vapour and nitrogen oxide perturbations resulting from aircraft emissions in the stratosphere. It is found that a fleet of 500 Boeing‐type sst's, flying at 20 km and 45°N in the summer hemisphere and inserting NO_x at a rate of 1.8 megatons per year, has the effect of reducing the global total ozone by 14.7%. Similar calculations for 342 Concorde/TU‐114's, cruising at 17 km and injecting NO_x at a rate of 0.35 megatons per year, show a global‐average total‐ozone reduction of 1.85%. Although water vapour is considered important, because of its ability to convert NO₂ into HNO₃, the direct effect on global‐average total‐ozone reduction resulting from the 100% increase in the stratospheric water content is less than 1%. The changes in the chemical structure (HO^NO^), temperature, and mean motions associated with the ozone reduction are also investigated in the case of the 1.8‐megaton‐per‐year NO_X perturbation. It is shown that the reduced meridional temperature gradient in the middle and upper stratosphere resulting from the NO_x perturbation leads to the weakening of the tropical easterly jet in the summer hemisphere and mid‐latitude westerlies in the winter season.

The sensitivity of the model solutions to an alternate choice of input parameters (diffusion coefficients and solar photodissociation data) is tested and the main deficiency of the model is pointed out.     

On the efficiency of horizontal diffusion and numerical filtering in an Arakawa‐type model

Abstract

In a series of 5‐day forecasts with a 3‐layer/2.8° hemispheric model, horizontal diffusion schemes of the second and fourth degree are compared with a numericalfilter technique. The results, which are discussed mainly in terms of spectral energetics in zonal wavenumber space, indicate that fourth‐degree diffusion is more scale selective than second‐degree and equivalent to filtering. The seventh‐order filter applied only intermittently to the prognostic variables is superior to fourth‐degree diffusion from the viewpoint of computational economy. Excessive dissipation of the long waves may inhibit the production of eddy kinetic energy from eddy available potential energy.     

SPOM: A regional model of the sub‐polar north Atlantic

Abstract

A regional model of the sub‐polar North Atlantic has been developed for use in process and variability studies of this important high‐latitude area. Open boundary conditions handle connections with the rest of the Atlantic Ocean at 38°N, while buffer zones are used in the northern boundary regions. Extensive testing and experimentation has led to a model which can reproduce major elements of the hydrography and circulation in the region, although limitations exist. A key model feature is the inclusion of a finite volume partial cell topographic representation that significantly improves the structure of the underlying bottom topography. Improvements include a tighter and sharper gyre structure, increased transports, sub‐polar mode water formation sites linked to the topographic slope along the outside of the gyre and a more reasonable representation of Labrador Sea water properties and dispersal pathways. The choice of inflow conditions for the open southern boundary affects the deep western boundary current, as well as the representation of Mediterranean Water, which has a significant effect on Labrador Sea water in the eastern basin.     

Sea‐ice interaction and the stability of the thermohaline circulation

Abstract

The role of sea‐ice in affecting the stability and long‐term variability of the oceanic thermohaline circulation (THC) is studied in this paper. The emphasis is placed on studying how sea‐ice might affect the stability and the long‐term variability of the THC through modulations of the surface heat and freshwater fluxes. A simple box model is analyzed to elucidate qualitatively the distinct physical meanings of these two processes. The analytical solution of this simple model indicates that, for the long timescales considered here, the thermal insulation stabilizes the THC while the freshwater feedback increases the effective inertia of the coupled ice‐ocean system. Sea‐ice insulation lessens the negative feedback between heat flux and the SST, and therefore, allows the SST to play a greater role in counteracting changes of the THC and high latitude salinity field. The freshwater feedback effectively links the surface heat flux to a freshwater reservoir, and thus, increases the effective inertia of the coupled ocean‐ice system. A two‐dimensional ocean model coupled with a thermodynamic sea‐ice model is used to estimate quantitatively the magnitudes of these two feedbacks. The numerical experiments involve the model's responses both to initial anomalies and to changes of forcing fields. For the free response cases (model responses to initial anomalies without changing the forcing fields), the model shows that the decay rate of an initial anomaly is greater when sea‐ice is included. For small perturbations the thermal insulation effect dominates over the freshwater feedback. The latter becomes increasingly more important for larger perturbations. In response to a change of external forcing, the presence of sea‐ice reduces the magnitude and the pace of the model's response. The numerical results are qualitatively consistent with the analytical solution of the box model.     

Some properties and comparative performance of the semi‐Lagrangian method of Robert in the solution of the advection‐diffusion equation

Abstract

A numerical method for solving the advection‐diffusion equation, based on the semi‐Lagrangian algorithm of Robert (1981, 1982) is described, analysed and evaluated in comparison with other methods through a series of test problems. It is found that this method is generally better than other semi‐Lagrangian methods, and is a viable alternative to existing methods for LRTAP and other meteorological modelling because of its flexibility in application, its computational stability and its accuracy.     

A decision‐analytic study of the joint value of seasonal precipitation and temperature forecasts in a choice‐of‐crop problem

Abstract

A static decision‐analytic method is used to investigate the economic value of bivariate ‐ precipitation and temperature ‐ seasonal forecasts of the form currently issued by the U.S. National Weather Service. This method is applied to a corn versus spring wheat choice‐of‐crop decision‐making problem by considering a transect of four counties across the northwestern margin of the North American corn belt. Numerical results indicate that seasonal forecasts of current quality can be of appreciable value (≥$1/ha) for some locations when the optimal action chosen on the basis of climatological information is only marginally preferred to another action. Increases in forecast value follow from hypothetical increases in the quality of both the precipitation and temperature components of the forecasts in the spring wheat region, whereas forecast value increases primarily as a function of the quality of the precipitation forecasts alone in the corn belt region. The results are very sensitive to absolute and relative crop prices.     

Integration of a semi‐implicit model with time‐dependent boundary conditions

A semi‐implicit barotropic primitive equations model is integrated over a limited area with time dependent boundary conditions using the standard mesh and a finer mesh. Following a theorem by Charney, a minimum number of variables are specified as boundary conditions for the limited area integrations in order to avoid mathematical over‐specification. The comparison of coarse mesh limited area forecasts with the corresponding forecasts made over a much larger domain demonstrates that the essential features, namely the Rossby type perturbations, are handled almost perfectly. The fine mesh forecasts over the same limited area are also very good. Finally, the effect of specifying inaccurate boundary conditions, in the form of twelve‐hour forecasts, is briefly illustrated.     

Reply to comments on and addenda to “some properties and comparative performance of the semi‐lagrangian method of Robert in the solution of the advection‐diffusion equation”

Abstract

Droughts are major natural disasters for many parts of the world. Dry areas, where the precipitation pattern is markedly seasonal, or is otherwise highly variable, are the most susceptible. The Canadian Prairies, together with the U.S. Great Plains, are one such area. While immediate loss of life is seldom a feature of most droughts, malnutrition and even starvation do follow severe droughts in some parts of the world. In Canada, economic losses, particularly in the agricultural sector, may reach several hundred millions of dollars in a drought year, with major socio‐economic repercussions affecting the entire region. Environmental damages include soil degradation and erosion, vegetation damage, slough and lake deterioration and wildlife loss. Unlike most other natural disasters, drought onset is difficult to identify. Droughts develop slowly, and until human activity begins to be affected by an on‐going reduction of precipitation, their existence is unrecognized. Development and application of specific soil moisture and drought indices based on cumulative precipitation deficits have enhanced drought monitoring programs. These in turn provide guidance on the need for mitigative measures that can be initiated early in the course of a drought. Any improvement in the timely application of these measures is, however, strongly dependent on being able to determine the drought's course, extent and likely severity at the earliest stage possible. The identification of precursor conditions for past drought has raised the possibility that the likelihood of a drought occurring in a particular year or growing season might be predictable. Teleconnections between North American precipitation patterns and ENSO events and other surface boundary conditions in the North Pacific have been detected. Forecasting seasonal temperature and perhaps precipitation anomalies appears to be potentially feasible using a suitable merging of precursor parameters and modelling methodologies. Clearly, future research must focus both on those precursors that have been identified and on a search for possible new ones. Development of better forecasting methodology is also essential. Research activity to identify and evaluate new mitigative measures should also be increased to keep pace with the prospects of drought predictability.     

Phase‐locked and asymmetric correlations of the wintertime atmospheric patterns with the ENSO

Abstract

Teleconnections between sea surface temperature (SST) anomalies over the Pacific and the dominant patterns of wintertime Northern Hemisphere 500‐hPa height are examined by applying statistical techniques such as rotated principal component analysis and composite analysis. It is shown that the Pacific/North American (PNA) patterns in December through March are correlated most significantly with the ENSO‐related SST anomalies in the previous October, while the western Pacific (WP) patterns in December through February are most closely linked to the ENSO‐related SST anomalies in the same season. In addition, the PNA response to the ENSO signal during La Niña events is more significant than that during El Niño events, while the WP response is stronger during El Niño events than during La Niña events. A composite analysis shows that in the El Niño winters the North Pacific centre of the PNA pattern is located about 10 degrees east of its normal position, leading to a less significant correlation between the ENSO signal and the PNA pattern in these winters.

The ENSO‐related SST anomalies include a large centre of action over the tropical Pacific and an oppositely signed anomaly centre over the North Pacific. The North Pacific centre appears to the west of the dateline in September and October. This ENSO‐related seed of SST anomalies slowly moves eastward in the following months, gradually cutting off its connection with SST anomalies over the tropical Pacific and being coupled with the PNA pattern. It is pointed out that, although the wintertime SST anomaly over the North Pacific may appear as a mode linearly independent of the ENSO signal in the same season, it is partially related to the ENSO signal in the preceding autumn.

Possible dynamical explanations of the above results are discussed. It is suggested that the WP pattern can be linked to the tropical Pacific heat source via advection of vorticity by the upper‐tropospheric divergent/convergent flow, and the intensification of vorticity gradients associated with a stronger east Asian jet is likely to be responsible for a more significant WP pattern response to the ENSO signal in the El Niño winters. On the other hand, the ENSO‐related PNA pattern could be considered a manifestation of the eastward extension (El Niño) or westward withdrawal (La Niña) of the east Asian jet stream due to the local Hadley cell over the Pacific. In addition, the ENSO‐related seed of extratropical SST anomaly over the western Pacific in autumn may also play an important role in the development of the PNA pattern in the following winter.     

2‐D model simulations of anomalous spring variations of ozone,temperature and other minor trace constituents in the southern polar region

Abstract

A 2‐D time‐dependent Eulerian model that takes into account the major interaction between chemistry and dynamical transports in the stratosphere is used to simulate, in a self‐consistent manner, the anomalous spring variations of ozone, temperature and other relevant minor constituents in the southern polar region. The normal transport coefficients and the standard chemistry have been altered to represent the anomalous dynamical and chemical conditions in the extremely cold lower stratosphere of Southern Hemisphere (south of 65° S) spring. The results show that it is necessary to invoke both dynamical and chemical mechanisms in order to explain the observed rapid spring decline of total ozone in the southern polar region.     

Non‐linear coupling between modes in a low‐dimensional model of ENSO

Abstract

An intermediate coupled model of the tropical Pacific ocean‐atmosphere system was reduced by projecting the non‐linear model onto a truncated basis set of its own empirical orthogonal functions (EOFs). For moderate coupling strengths, the simulated El Niño/Southern Oscillation (ENSO) variability consists of a dominant quasi‐quadrennial mode with a period of approximately four years and a smaller quasi‐biennial mode at a period of approximately two years. In the absence of a seasonal cycle, the leading two EOFs capture the dynamics of the leading interannual mode, with a further two EOFs being required to capture the secondary oscillation. The presence of seasonal forcing increases the EOF requirement by two, the leading pair of EOFs being dominated by the annual cycle. Normal mode analysis of the reduced models indicates that the quasi‐biennial mode manifests itself, even though it is linearly stable, by non‐linear coupling to the quasi‐quadrennial mode. The nonlinearity does not produce the quasi‐biennial signal unless the spatial degrees of freedom associated with the linear quasi‐biennial mode are present. Other linearly stable modes also couple non‐linearly to the leading interannual mode and to the seasonal cycle, but the quasi‐biennial mode is favoured over other, less‐damped linear modes because of its proximity to a multiple of the quasi‐quadrennial frequency.     

On tidal resonance in the global ocean and the back‐effect of coastal tides upon open‐ocean tides

Abstract

The resonance of semi‐diurnal tidal elevations is investigated with a forward numerical forced damped global tide model and an analytical model of forced‐damped tides in a deep ocean basin coupled to a shelf. The analytical model contains the classical half‐wavelength and quarter‐wavelength resonances in the deep ocean and shelf, respectively, as well as a forcing‐scale dependence which depends on the ratio of the phase speed of open‐ocean gravity waves to that of the astronomical forcing. In the analytical model, when the deep ocean and shelf resonate separately at the same frequency, the resonance in the coupled system shifts to frequencies slightly higher and lower than the original frequency, such that a ‘double bump’ is seen in plots of elevation amplitude versus frequency. The addition of a shelf to a resonant open ocean tends to reduce open‐ocean tides, especially when the shelf is also near resonance. The magnitude of this ‘back‐effect’ is controlled by shelf friction. A weakly damped resonant shelf has a larger back‐effect on the open‐ocean tide than does a strongly damped shelf. Numerical simulations largely bear out the analytical model predictions, at least qualitatively. Idealized simulations show that continents enhance tides by enabling the half‐wavelength resonance. Simulations with realistic geometry and topography but varying longitudinal structure in the astronomical forcing display an influence of the forcing scale on tidal amplitudes somewhat similar to that seen in the analytical model. A frequency sweep in the semi‐diurnal band in experiments with realistic geometry and topography reveals weakly resonant peaks in the amplitudes of several shelf regions and in the globally averaged open‐ocean amplitudes. Finally, the back‐effect of the shelf upon the open ocean is seen in simulations in which locations of resonant coastal tides are blocked out and open‐ocean tidal elevations are significantly altered (increased, generally) as a result.     

Comparison between measured and modelled cloudless‐sky solar irradiances at the surface and at the top of the atmosphere

Abstract

Cloudless‐sky solar fluxes calculated by the radiative transfer algorithm used in the Canadian Climate Centre's general circulation climate model are compared with measurements of upwelling radiation at the top of the atmosphere (TOA) and downwelling radiation at the surface. The 12‐layer model partitions the solar spectrum into two broad wavebands (0.25–0.68 and 0.68–4.00 μm). The comparison utilized TOA fluxes estimated from Nimbus‐ 7 measurements and measured downwelling fluxes at the surface for Kalgoorlie, West Australia, and downwelling fluxes at the surface for Woodbridge, Ontario. Model estimates and measurements agreed to within experimental error for most solar zenith angles. Estimates improved, especially at Woodbridge, when aerosol effects were included. The mean bias error was less than 4% for surface irradiance and less than 6% for upwelling TOA irradiance, which produces a TOA albedo error of about 0.01.     

Decadal‐to‐interdecadal fluctuations of arctic sea‐ice cover and the atmospheric circulation during 1954–1994

Abstract

The relationship between the Arctic and subarctic sea‐ice concentration (SIC) anomalies, particularly those associated with the decadal‐scale Greenland and Labrador Seas “Ice and Salinity Anomalies (ISAs) “, and the overlying atmospheric circulation fluctuations is investigated using the singular value decomposition (SVD) and composite map analysis methods. The data analyzed are monthly SIC and sea level pressure (SLP) anomalies, which cover the northern hemisphere poleward of 45°N and extend over the 41‐year period 1954–1994.

The SVD₁ (first) mode of the coupled variability, which accounts for 57% of the square covariance, is for the most part an atmosphere‐to‐ice forcing mode characterized by the decadal timescale. The aforementioned ISA anomalies are clearly captured by this mode whose SIC anomalies are dominated by a strong dipole across Greenland. However, as part of the same mode, there is also a weaker SIC dipole in the northern North Pacific which has opposite‐signed anomalies in the Sea of Okhotsk and the Bering Sea. It is also shown that there exists a significant negative correlation between the decadal SIC variability in the Greenland‐Barents Seas region associated with this mode and the North Atlantic Oscillation, whose spectrum also exhibits a quasi‐decadal signal.

The SVD₂ mode accounts for 12% of the square covariance and shows no evidence of a dominant forcing field of either SIC or SLP. This SVD mode exhibits very low frequency (interdecadal) variability, and its co‐variability is mainly concentrated in the northern North Pacific. It appears to be a high‐latitude extension of the recently investigated interdecadal North Pacific Oscillation. The spatial structure of the second mode complements the case of the first SVD mode whose co‐variability mainly occurs in the northern North Atlantic.     

Cross‐shore variations of near‐surface wind velocity and atmospheric turbulence at the land‐sea boundary during CASP

Abstract

Airborne measurements of mean wind velocity and turbulence in the atmospheric boundary layer under wintertime conditions of cold offshore advection suggest that at a height of 50 m the mean wind speed increases with offshore distance by roughly 20% over a horizontal scale of order 10 km. Similarly, the vertical gust velocity and turbulent kinetic energy decay on scales of order 3.5 km by factors of 1.5 and 3.2, respectively. The scale of cross‐shore variations in the vertical fluxes of heat and downwind momentum is also 10 km, and the momentum flux is found to be roughly constant to 300 m, whereas the heat flux decreases with height. The stability parameter, z/L (where z = 50 m and L is the local Monin‐Obukhov length), is generally small over land but may reach order one over the warm ocean. The magnitude and horizontal length scales associated with the offshore variations in wind speed and turbulence are reasonably consistent with model results for a simple roughness change, but a more sophisticated model is required to interpret the combined effects of surface roughness and heat flux contrasts between land and sea.

Comparisons between aircraft and profile‐adjusted surface measurements of wind speed indicate that Doppler biases of 1–2 m s^?1 in the aircraft data caused by surface motions must be accounted for. In addition, the wind direction measurements of the Minimet anemometer buoy deployed in CASP are found to be in error by 25 ± 5°, possibly due to a misalignment of the anemometer vane. The vertical fluxes of heat and momentum show reasonably good agreement with surface estimates based on the Minimet data.     

The influence of Hudson Bay runoff and ice‐melt on the salinity of the inner Newfoundland Shelf

Abstract

The present study examines sources of the interannual variability in salinity on the Newfoundland continental shelf observed in a 40‐year time series from an oceanographic station known as Station 27. Specifically, we investigate, through lag‐correlation analysis, the a priori hypotheses that the salinity anomalies at Station 27 are determined by freshwater runoff anomalies from Hudson and Ungava bays and by ice‐melt anomalies in Hudson Bay and on the Labrador Shelf. Interannual variations of summer runoff into Hudson Bay were significantly negatively correlated with salinity anomalies on the Newfoundland Shelf with a lag (9 months) that is consistent with expected travel times based on known current velocities in Hudson Bay and along the Labrador Shelf. Sea‐ice extent over the Labrador and northern Newfoundland shelves was significantly negatively correlated with salinity at a lag of 3 to 4 months, corresponding to the time of minimum salinity at Station 27. It appears that ice‐melt over the Labrador‐northern Newfoundland Shelf is primarily responsible for the seasonal salinity minimum over the Newfoundland Shelf. Interannual variability in runoff into Ungava Bay and ice‐melt in Hudson Bay were not correlated with interannual salinity variations on the Newfoundland Shelf.     

The dynamics of moist frontogenesis in a semi‐geostrophic model

Abstract

The effects of condensational heating on the semi‐geostrophic dynamics of frontogenesis are studied using a two‐dimensional deformation model. The model includes water vapour and allows the formation of stratiform clouds. Analysis and numerical results show that heating due to stratiform clouds has the effect of reducing stability to slantwise convection, as found in previous studies (Thorpe and Emanuel, 1985). In addition, heating‐induced potential vorticity and temperature anomalies play a very important role in the frontal circulation. The ageostrophic flow induced by these anomalies tends to reinforce the effect of heating and increases the strength of frontal cloud. The model is also able to produce the low‐level jet maximum ahead of a cold front at an elevated level, in agreement with observations, owing to the explicit condensation scheme used in the model.     

Initialisation and predictability of the AMOC over the last 50 years in a climate model

The mechanisms involved in Atlantic meridional overturning circulation (AMOC) decadal variability and predictability over the last 50 years are analysed in the IPSL–CM5A–LR model using historical and initialised simulations. The initialisation procedure only uses nudging towards sea surface temperature anomalies with a physically based restoring coefficient. When compared to two independent AMOC reconstructions, both the historical and nudged ensemble simulations exhibit skill at reproducing AMOC variations from 1977 onwards, and in particular two maxima occurring respectively around 1978 and 1997. We argue that one source of skill is related to the large Mount Agung volcanic eruption starting in 1963, which reset an internal 20-year variability cycle in the North Atlantic in the model. This cycle involves the East Greenland Current intensity, and advection of active tracers along the subpolar gyre, which leads to an AMOC maximum around 15 years after the Mount Agung eruption. The 1997 maximum occurs approximately 20 years after the former one. The nudged simulations better reproduce this second maximum than the historical simulations. This is due to the initialisation of a cooling of the convection sites in the 1980s under the effect of a persistent North Atlantic oscillation (NAO) positive phase, a feature not captured in the historical simulations. Hence we argue that the 20-year cycle excited by the 1963 Mount Agung eruption together with the NAO forcing both contributed to the 1990s AMOC maximum. These results support the existence of a 20-year cycle in the North Atlantic in the observations. Hindcasts following the CMIP5 protocol are launched from a nudged simulation every 5 years for the 1960–2005 period. They exhibit significant correlation skill score as compared to an independent reconstruction of the AMOC from 4-year lead-time average. This encouraging result is accompanied by increased correlation skills in reproducing the observed 2-m air temperature in the bordering regions of the North Atlantic as compared to non-initialized simulations. To a lesser extent, predicted precipitation tends to correlate with the nudged simulation in the tropical Atlantic. We argue that this skill is due to the initialisation and predictability of the AMOC in the present prediction system. The mechanisms evidenced here support the idea of volcanic eruptions as a pacemaker for internal variability of the AMOC. Together with the existence of a 20-year cycle in the North Atlantic they propose a novel and complementary explanation for the AMOC variations over the last 50 years.     

Sensitivity of the CRCM atmospheric and the Gulf of St. Lawrence Ocean‐Ice models to each other

Abstract

The sensitivity of the Canadian Regional Climate Model (CRCM), developed at the Université du Québec à Montréal, and the Gulf of St. Lawrence Ocean Model (GOM), developed at the Institut Maurice‐ Lamontagne, to each other is tested with an ensemble of simulations over eastern Canada from 1 November 1989 to 31 March 1990. The goal of this study is to investigate the interaction of the CRCM and GOM with respect to each other's forcing fields. In the first part of the experiment, a series of simulations were performed using an iterative strategy, where both models run separately and alternately, using variables from the other model to supply the needed forcing fields for the computation of surface fluxes. The runs are iterated several times over the same period from the output of the previous run to allow the atmosphere and the ocean to interact several times with each other and to study the evolution of the solutions from one iteration to the next. In the second part of the experiment, a two‐way coupled simulation is performed over the same period. The results indicate that on a monthly or longer timescale, the CRCM is not very sensitive to the details of the oceanic fields from GOM, except locally over the Gulf of St. Lawrence (GSL). However, GOM is quite sensitive to the differences in atmospheric fields from the CRCM. The results of several iterations converge to a unique solution, suggesting that the CRCM and GOM reach equilibrium with respect to each other's forcing fields. Furthermore, the results of the coupled run also converge to this same solution.     

The economics of decarbonizing the energy system—results and insights from the RECIPE model intercomparison

This paper synthesizes the results from the model intercomparison exercise among regionalized global energy-economy models conducted in the context of the RECIPE project. The economic adjustment effects of long-term climate policy are investigated based on the cross-comparison of the intertemporal optimization models ReMIND-R and WITCH as well as the recursive dynamic computable general equilibrium model IMACLIM-R. A number of robust findings emerge. If the international community takes immediate action to mitigate climate change, the costs of stabilizing atmospheric CO₂ concentrations at 450?ppm (roughly 530?C550?ppm-e) discounted at 3% are estimated to be 1.4% or lower of global consumption over the twenty-first century. Second best settings with either a delay in climate policy or restrictions to the deployment of low-carbon technologies can result in substantial increases of mitigation costs. A delay of global climate policy until 2030 would render the 450?ppm target unachievable. Renewables and CCS are found to be the most critical mitigation technologies, and all models project a rapid switch of investments away from freely emitting energy conversion technologies towards renewables, CCS and nuclear. Concerning end use sectors, the models consistently show an almost full scale decarbonization of the electricity sector by the middle of the twenty-first century, while the decarbonization of non-electric energy demand, in particular in the transport sector remains incomplete in all mitigation scenarios. The results suggest that assumptions about low-carbon alternatives for non-electric energy demand are of key importance for the costs and achievability of very low stabilization scenarios.