Atmospheric blocking in a high resolution climate model: influences of mean state, orography and eddy forcing

An underestimate of atmospheric blocking occurrence is a well-known limitation of many climate models. This article presents an analysis of Northern Hemisphere winter blocking in an atmospheric model with increased horizontal resolution. European blocking frequency increases with model resolution, and this results from an improvement in the atmospheric patterns of variability as well as a simple improvement in the mean state. There is some evidence that the transient eddy momentum forcing of European blocks is increased at high resolution, which could account for this. However, it is also shown that the increase in resolution of the orography is needed to realise the improvement in blocking, consistent with the increase in height of the Rocky Mountains acting to increase the tilt of the Atlantic jet stream and giving higher mean geopotential heights over northern Europe. Blocking frequencies in the Pacific sector are also increased with atmospheric resolution, but in this case the improvement in orography actually leads to a decrease in blocking

Simulation of late-twenty-first-century changes in wintertime atmospheric circulation over Europe due to anthropogenic causes

Evidence is presented, based on an ensemble of climate change scenarios performed with a global general circulation model of the atmosphere with high horizontal resolution over Europe, to suggest that the end-of-century anthropogenic climate change over the North Atlantic--European region strongly projects onto the positive phase of the North Atlantic Oscillation during wintertime. It is reflected in a doubling of the residence frequency of the climate system in the associated circulation regime, in agreement with the nonlinear climate perspective. The strong increase in the amplitude of the response, compared to coarse-resolution coupled model studies, suggests that improved model representation of regional climate is needed to achieve more reliable projections of anthropogenic climate change on European climate.

An oceanic origin for the increase of atmospheric radiocarbon during the Younger Dryas

Variations in carbon-14 to carbon-12 ratio in the atmosphere (Δ14Catm) provide a powerful diagnostic for elucidating the timing and nature of geophysical and anthropological change. The (Atlantic) marine archive suggests a rapid Δ14Catm increase of 50‰ at the onset of the Younger Dryas (YD) cold reversal (12.9–11.7 kyr BP), which has not yet been satisfactorily explained in terms of magnitude or causal mechanism, as either a change in ocean ventilation or production rate. Using Earth-system model simulations and comparison of marine-based radiocarbon records from different ocean basins, we demonstrate that the YD Δ14Catm increase is smaller than suggested by the marine archive. This is due to changes in reservoir age, predominantly caused by reduced ocean ventilation.

The atmospheric impact of uncertainties in recent Arctic Sea ice reconstructions

There are significant discrepancies between observational datasets of Arctic sea ice concentrations covering the last three decades, which result in differences of over 20% in Arctic summer sea ice extent/area and 5%–10% in winter. Previous modeling studies have shown that idealized sea ice anomalies have the potential for making a substantial impact on climate. In this paper, this theory is further developed by performing a set of simulations using the third Hadley Centre Coupled Atmospheric Model (HadAM3). The model was driven with monthly climatologies of sea ice fractions derived from three of these records to investigate potential implications of sea ice inaccuracies for climate simulations. The standard sea ice climatology from the Met Office provided a control. This study focuses on the effects of actual inaccuracies of concentration retrievals, which vary spatially and are larger in summer than winter. The smaller sea ice discrepancies in winter have a much larger influence on climate than the much greater summer sea ice differences. High sensitivity to sea ice prescription was observed, even though no SST feedbacks were included. Significant effects on surface fields were observed in the Arctic, North Atlantic, and North Pacific. Arctic average surface air temperature anomalies in winter vary by 2.5°C, and locally exceed 12°C. Arctic mean sea level pressure varies by up to 5 mb locally. Anomalies extend to 45°N over North America and Eurasia but not to lower latitudes, and with limited changes in circulation above the boundary layer. No statistically significant impact on climate variability was simulated, in terms of the North Atlantic Oscillation. Results suggest that the uncertainty in summer sea ice prescription is not critical but that winter values require greater accuracy, with the caveats that the influences of ocean–sea ice feedbacks were not included in this study.

Impact of atmospheric forcing on Antarctic continental shelf water masses

The Antarctic continental shelf seas feature a bimodal distribution of water mass temperature, with the Amundsen and Bellingshausen Seas flooded by Circumpolar Deep Water that is several degrees Celsius warmer than the cold shelf waters prevalent in the Weddell and Ross Seas. This bimodal distribution could be caused by differences in atmospheric forcing, ocean dynamics, ocean and ice feedbacks, or some combination of these factors. In this study, a highly simplified coupled sea ice–mixed layer model is developed to investigate the physical processes controlling this situation. Under regional atmospheric forcings and parameter choices the 10-yr simulations demonstrate a complete destratification of the Weddell Sea water column in winter, forming cold, relatively saline shelf waters, while the Amundsen Sea winter mixed layer remains shallower, allowing a layer of deep warm water to persist. Applying the Weddell atmospheric forcing to the Amundsen Sea model destratifies the water column after two years, and applying the Amundsen forcing to the Weddell Sea model results in a shallower steady-state winter mixed layer that no longer destratifies the water column. This suggests that the regional difference in atmospheric forcings alone is sufficient to account for the bimodal distribution in Antarctic shelf-sea temperatures. The model prediction of mixed layer depth is most sensitive to the air temperature forcing, but a switch in all forcings is required to prevent destratification of the Weddell Sea water column.

The integrated WRF/urban modelling system: development, evaluation, and applications to urban environmental problems

To bridge the gaps between traditional mesoscale modelling and microscale modelling, the National Center for Atmospheric Research, in collaboration with other agencies and research groups, has developed an integrated urban modelling system coupled to the weather research and forecasting (WRF) model as a community tool to address urban environmental issues. The core of this WRF/urban modelling system consists of the following: (1) three methods with different degrees of freedom to parameterize urban surface processes, ranging from a simple bulk parameterization to a sophisticated multi-layer urban canopy model with an indoor–outdoor exchange sub-model that directly interacts with the atmospheric boundary layer, (2) coupling to fine-scale computational fluid dynamic Reynolds-averaged Navier–Stokes and Large-Eddy simulation models for transport and dispersion (T&D) applications, (3) procedures to incorporate high-resolution urban land use, building morphology, and anthropogenic heating data using the National Urban Database and Access Portal Tool (NUDAPT), and (4) an urbanized high-resolution land data assimilation system. This paper provides an overview of this modelling system; addresses the daunting challenges of initializing the coupled WRF/urban model and of specifying the potentially vast number of parameters required to execute the WRF/urban model; explores the model sensitivity to these urban parameters; and evaluates the ability of WRF/urban to capture urban heat islands, complex boundary-layer structures aloft, and urban plume T&D for several major metropolitan regions. Recent applications of this modelling system illustrate its promising utility, as a regional climate-modelling tool, to investigate impacts of future urbanization on regional meteorological conditions and on air quality under future climate change scenarios. Copyright © 2010 Royal Meteorological Society

Ice Shelf Water plume flow beneath Filchner-Ronne Ice Shelf, Antarctica

[1] A two-dimensional plume model is used to study the interaction between Filchner-Ronne Ice Shelf, Antarctica and its underlying ocean cavity. Ice Shelf Water (ISW) plumes are initiated by the freshwater released from a melting ice shelf and, if they rise, may become supercooled and deposit marine ice due to the pressure increase in the in situ freezing temperature. The aim of this modeling study is to determine the origin of the thick accretions of marine ice at the base of Filchner-Ronne Ice Shelf and thus improve our understanding of ISW flow paths. The model domain is defined from measurements of ice shelf draft, and from this ISW the model is able to predict plumes that exit the cavity in the correct locations. The modeled plumes also produce basal freezing rates that account for measured marine ice thicknesses in the western part of Ronne Ice Shelf. We find that the freezing rate and plume properties are significantly influenced by the confluence of plumes from different meltwater sources. We are less successful in matching observations of marine ice under the rest of Filchner-Ronne Ice Shelf, which we attribute primarily to this model’s neglect of circulations in the ocean outside the plume.

Northern Hemisphere atmospheric blocking as simulated by 15 atmospheric general circulation models in the period 1979-1988

As a part of the Atmospheric Model Intercomparison Project (AMIP), the behaviour of 15 general circulation models has been analysed in order to diagnose and compare the ability of the different models in simulating Northern Hemisphere midlatitude atmospheric blocking. In accordance with the established AMIP procedure, the 10-year model integrations were performed using prescribed, time-evolving monthly mean observed SSTs spanning the period January 1979–December 1988. Atmospheric observational data (ECMWF analyses) over the same period have been also used to verify the models results. The models involved in this comparison represent a wide spectrum of model complexity, with different horizontal and vertical resolution, numerical techniques and physical parametrizations, and exhibit large differences in blocking behaviour. Nevertheless, a few common features can be found, such as the general tendency to underestimate both blocking frequency and the average duration of blocks. The problem of the possible relationship between model blocking and model systematic errors has also been assessed, although without resorting to ad-hoc numerical experimentation it is impossible to relate with certainty particular model deficiencies in representing blocking to precise parts of the model formulation.

Some results from an intercomparison of the climates simulated by 14 atmospheric general circulation models

Some climatological information from 14 atmospheric general circulation models is presented and compared in order to assess the ability of a broad group of models to simulate current climate. The quantities considered are cross sections of temperature, zonal wind, and meridional stream function together with latitudinal distributions of mean sea level pressure and precipitation rate. The nature of the deficiencies in the simulated climates that are common to all models and those which differ among models is investigated; the general improvement in the ability of models to simulate certain aspects of the climate is shown; consideration is given to the effect of increasing resolution on simulated climate; and approaches to understanding and reducing model deficiencies are discussed. The information presented here is a subset of a more voluminous compilation which is available in report form (Boer et al., 1991). This report contains essentially the same text, but results from all 14 models are presented together with additional results in the form of geographical distributions of surface variables and certain difference statistics.

An intercomparison of the climates simulated by 14 atmospheric general circulation models

Climatological information from fourteen atmospheric general circulation models is presented and compared in order to assess the ability of a broad group of models to simulate current climate. The quantities considered are cross sections of temperature, zonal wind and meridional stream function together with latitudinal distributions of mean sea-level pressure and precipitation rate. The nature of the deficiencies in the simulated climates that are common to all models and those which differ among models is investigated, general improvement in the ability of models to simulate certain aspects of the climate is shown, consideration is given to the effect of increasing resolution on simulated climate and approaches to the understanding and reduction of model deficiencies are discussed.

Northern Hemisphere atmospheric blocking as simulated by 15 atmospheric general circulation models in the period 1979-1988 (Results from an AMIP diagnostic subproject)

As a part of the Atmospheric Model Intercomparison Project (AMIP), the behaviour of 15 general circulation models has been analysed in order to diagnose and compare the ability of the different models in simulating Northern Hemisphere midlatitude atmospheric blocking. In accordance with the established AMIP procedure, the 10-year model integrations were performed using prescribed, time-evolving monthly mean observed SSTs spanning the period January 1979–December 1988. Atmospheric observational data (ECMWF analyses) over the same period have been also used to verify the models results. The models involved in this comparison represent a wide spectrum of model complexity, with different horizontal and vertical resolution, numerical techniques and physical parametrizations, and exhibit large differences in blocking behaviour. Nevertheless, a few common features can be found, such as the general tendency to underestimate both blocking frequency and the average duration of blocks. The problem of the possible relationship between model blocking and model systematic errors has also been assessed, although without resorting to ad-hoc numerical experimentation it is impossible to relate with certainty particular model deficiencies in representing blocking to precise parts of the model formulation.

Horizontal and vertical structure of the Eyjafjallajökull ash cloud over the UK: a comparison of airborne lidar observations and simulations

During April and May 2010 the ash cloud from the eruption of the Icelandic volcano Eyjafjallajökull caused widespread disruption to aviation over northern Europe. The location and impact of the eruption led to a wealth of observations of the ash cloud were being obtained which can be used to assess modelling of the long range transport of ash in the troposphere. The UK FAAM (Facility for Airborne Atmospheric Measurements) BAe-146-301 research aircraft overflew the ash cloud on a number of days during May. The aircraft carries a downward looking lidar which detected the ash layer through the backscatter of the laser light. In this study ash concentrations derived from the lidar are compared with simulations of the ash cloud made with NAME (Numerical Atmospheric-dispersion Modelling Environment), a general purpose atmospheric transport and dispersion model. The simulated ash clouds are compared to the lidar data to determine how well NAME simulates the horizontal and vertical structure of the ash clouds. Comparison between the ash concentrations derived from the lidar and those from NAME is used to define the fraction of ash emitted in the eruption that is transported over long distances compared to the total emission of tephra. In making these comparisons possible position errors in the simulated ash clouds are identified and accounted for. The ash layers seen by the lidar considered in this study were thin, with typical depths of 550–750 m. The vertical structure of the ash cloud simulated by NAME was generally consistent with the observed ash layers, although the layers in the simulated ash clouds that are identified with observed ash layers are about twice the depth of the observed layers. The structure of the simulated ash clouds were sensitive to the profile of ash emissions that was assumed. In terms of horizontal and vertical structure the best results were obtained by assuming that the emission occurred at the top of the eruption plume, consistent with the observed structure of eruption plumes. However, early in the period when the intensity of the eruption was low, assuming that the emission of ash was uniform with height gives better guidance on the horizontal and vertical structure of the ash cloud. Comparison of the lidar concentrations with those from NAME show that 2–5% of the total mass erupted by the volcano remained in the ash cloud over the United Kingdom.

The effect of a new drag-law parameterization on ice shelf water plume dynamics

A drag law accounting for Ekman rotation adjacent to a flat, horizontal bou ndary is proposed for use in a plume model that is written in terms of the depth-mean velocity. The drag l aw contains a variable turning angle between the mean velocity and the drag imposed by the turbulent bound ary layer. The effect of the variable turning angle in the drag law is studied for a plume of ice shelf wat er (ISW) ascending and turning beneath an Antarctic ice shelf with draft decreasing away from the groundi ng line. As the ISW plume ascends the sloping ice shelf–ocean boundary, it can melt the ice shelf, wh ich alters the buoyancy forcing driving the plume motion. Under these conditions, the typical turning ang le is of order 10° over most of the plume area for a range of drag coefficients (the minus sign arises for th e Southern Hemisphere). The rotation of the drag with respect to the mean velocity is found to be signifi cant if the drag coefficient exceeds 0.003; in this case the plume body propagates farther along and across the b ase of the ice shelf than a plume with the standard quadratic drag law with no turning angle.

Validity of the ice shelf water plume concept beneath Filchner-Ronne Ice Shelf

In winter, brine rejection from sea ice formation and export in the Weddell Sea, offshore of Filchner-Ronne Ice Shelf (FRIS), leads to the formation of High Salinity Shelf Water (HSSW). This dense water mass enters the cavity beneath FRIS by sinking southward down the sloping continental shelf towards the grounding line. Melting occurs when the HSSW encounters the ice shelf, and the meltwater released cools and freshens the HSSW to form a water mass known as Ice Shelf Water (ISW). If this ISW rises, the ‘ice pump’ is initiated (Lewis and Perkin, 1986), whereby the ascending ISW becomes supercooled and deposits marine ice at shallower locations due to the pressure increase in the in-situ freezing temperature. Sandh¨ager et al. (2004) were able to infer the thickness patterns of marine ice deposits at the base of FRIS (figure 1), so the primary aim of this work is to try to understand the ocean flows that determine these patterns. The plume model we use to investigate ISW flow is described fully by Holland and Feltham (accepted) so only a relatively brief outline is presented here. The plume is simulated by combining a parameterisation of ice shelf basal interaction and a multiplesize- class frazil dynamics model with an unsteady, depth-averaged reduced-gravity plume model. In the model an active region of ISW evolves above and within an expanse of stagnant ambient fluid, which is considered to be ice-free and has fixed profiles of temperature and salinity. The two main assumptions of the model are that there is a well-mixed layer underneath the ice shelf and that the ambient fluid outside the plume is stagnant with fixed properties. The topography of the ice shelf that the plume flows beneath is set to the FRIS ice shelf draft calculated by Sandh¨ager et al. (2004) masked with the grounding line from the Antarctic Digital Database (ADD Consortium, 2002). To initiate the plumes, we assume that the intrusion of dense HSSW initially causes melting at the points on the grounding line where the glaciological tributaries feeding FRIS go afloat.