The use of trajectory cluster analysis to interpret trace gas measurements at Mace Head, Ireland

A first step in interpreting the wide variation in trace gas concentrations measured over time at a given site is to classify the data according to the prevailing weather conditions. In order to classify measurements made during two intensive field campaigns at Mace Head, on the west coast of Ireland, an objective method of assigning data to different weather types has been developed. Air-mass back trajectories calculated using winds from ECMWF analyses, arriving at the site in 1995–1997, were allocated to clusters based on a statistical analysis of the latitude, longitude and pressure of the trajectory at 12 h intervals over 5 days. The robustness of the analysis was assessed by using an ensemble of back trajectories calculated for four points around Mace Head. Separate analyses were made for each of the 3 years, and for four 3-month periods. The use of these clusters in classifying ground-based ozone measurements at Mace Head is described, including the need to exclude data which have been influenced by local perturbations to the regional flow pattern, for example, by sea breezes. Even with a limited data set, based on 2 months of intensive field measurements in 1996 and 1997, there are statistically significant differences in ozone concentrations in air from the different clusters. The limitations of this type of analysis for classification and interpretation of ground-based chemistry measurements are discussed.

Sting jets in simulations of a real cyclone by two mesoscale models

The existence of sting jets as a potential source of damaging surface winds during the passage of extratropical cyclones has recently been recognized However, there are still very few published studies on the subject Furthermore, although ills known that other models are capable of reproducing sting jets, in the published literature only one numerical model [the Met Office Unified Model (MetUM)] has been used to numerically analyze these phenomena This article alms to improve our understanding of the processes that contribute to the development of sting jets and show that model differences affect the evolution of modeled sting jets A sting jet event during the passage of a cyclone over the United Kingdom on 26 February 2002 has been simulated using two mesoscale models namely the MetUM and the Consortium for Small Scale Modeling (COSMO) model to compare their performance Given the known critical importance of vertical resolution in the simulation of sting jets the vertical resolution of both models has been enhanced with respect to their operational versions Both simulations have been verified against surface measurements of maximum gusts, satellite imagery and Met Office operational synoptic analyses, as well as operational analyses from the ECMWF It is shown that both models are capable of reproducing sting jets with similar, though not identical. features Through the comparison of the results from these two models, the relevance of physical mechanisms, such as evaporative cooling and the release of conditional symmetric instability, in the generation and evolution of sting jets is also discussed

Equatorial waves in opposite QBO phases

A methodology for identifying equatorial waves is used to analyze the multilevel 40-yr ECMWF Re-Analysis (ERA-40) data for two different years (1992 and 1993) to investigate the behavior of the equatorial waves under opposite phases of the quasi-biennial oscillation (QBO). A comprehensive view of 3D structures and of zonal and vertical propagation of equatorial Kelvin, westward-moving mixed Rossby–gravity (WMRG), and n = 1 Rossby (R1) waves in different QBO phases is presented. Consistent with expectation based on theory, upward-propagating Kelvin waves occur more frequently during the easterly QBO phase than during the westerly QBO phase. However, the westward-moving WMRG and R1 waves show the opposite behavior. The presence of vertically propagating equatorial waves in the stratosphere also depends on the upper tropospheric winds and tropospheric forcing. Typical propagation parameters such as the zonal wavenumber, zonal phase speed, period, vertical wavelength, and vertical group velocity are found. In general, waves in the lower stratosphere have a smaller zonal wavenumber, shorter period, faster phase speed, and shorter vertical wavelength than those in the upper troposphere. All of the waves in the lower stratosphere show an upward group velocity and downward phase speed. When the phase of the QBO is not favorable for waves to propagate, their phase speed in the lower stratosphere is larger and their period is shorter than in the favorable phase, suggesting Doppler shifting by the ambient flow and a filtering of the slow waves. Tropospheric WMRG and R1 waves in the Western Hemisphere also show upward phase speed and downward group velocity, with an indication of their forcing from middle latitudes. Although the waves observed in the lower stratosphere are dominated by “free” waves, there is evidence of some connection with previous tropical convection in the favorable year for the Kelvin waves in the warm water hemisphere and WMRG and R1 waves in the Western Hemisphere, which is suggestive of the importance of convective forcing for the existence of propagating coupled Kelvin waves and midlatitude forcing for the existence of coupled WMRG and R1 waves.

Cloud-resolving simulations of deep convection over a heated mountain

Cloud-resolving numerical simulations of airflow over a diurnally heated mountain ridge are conducted to explore the mechanisms and sensitivities of convective initiation under high pressure conditions. The simulations are based on a well-observed convection event from the Convective and Orographically Induced Precipitation Study (COPS) during summer 2007, where an isolated afternoon thunderstorm developed over the Black Forest mountains of central Europe, but they are idealized to facilitate understanding and reduce computational expense. In the conditionally unstable but strongly inhibited flow under consideration, sharp horizontal convergence over the mountain acts to locally weaken the inhibition and moisten the dry midtroposphere through shallow cumulus detrainment. The onset of deep convection occurs not through the deep ascent of a single updraft but rather through a rapid succession of thermals that are vented through the mountain convergence zone into the deepening cloud mass. Emerging thermals rise through the saturated wakes of their predecessors, which diminishes the suppressive effects of entrainment and allows for rapid glaciation above the freezing level as supercooled cloud drops rime onto preexisting ice particles. These effects strongly enhance the midlevel cloud buoyancy and enable rapid ascent to the tropopause. The existence and vigor of the convection is highly sensitive to small changes in background wind speed U0, which controls the strength of the mountain convergence and the ability of midlevel moisture to accumulate above the mountain. Whereas vigorous deep convection develops for U0 = 0 m s−1, deep convection is completely eliminated for U0 = 3 m s−1. Although deep convection is able to develop under intermediate winds (U0 = 1.5 m s−1), its formation is highly sensitive to small-amplitude perturbations in the initial flow.

Ensemble prediction for nowcasting with a convection-permitting model—I: description of the system and the impact of radar-derived surface precipitation rates

A key strategy to improve the skill of quantitative predictions of precipitation, as well as hazardous weather such as severe thunderstorms and flash floods is to exploit the use of observations of convective activity (e.g. from radar). In this paper, a convection-permitting ensemble prediction system (EPS) aimed at addressing the problems of forecasting localized weather events with relatively short predictability time scale and based on a 1.5 km grid-length version of the Met Office Unified Model is presented. Particular attention is given to the impact of using predicted observations of radar-derived precipitation intensity in the ensemble transform Kalman filter (ETKF) used within the EPS. Our initial results based on the use of a 24-member ensemble of forecasts for two summer case studies show that the convective-scale EPS produces fairly reliable forecasts of temperature, horizontal winds and relative humidity at 1 h lead time, as evident from the inspection of rank histograms. On the other hand, the rank histograms seem also to show that the EPS generates too much spread for forecasts of (i) surface pressure and (ii) surface precipitation intensity. These may indicate that for (i) the value of surface pressure observation error standard deviation used to generate surface pressure rank histograms is too large and for (ii) may be the result of non-Gaussian precipitation observation errors. However, further investigations are needed to better understand these findings. Finally, the inclusion of predicted observations of precipitation from radar in the 24-member EPS considered in this paper does not seem to improve the 1-h lead time forecast skill.

Exploring the land–ocean contrast in convective vigor using islands

Moist convection is well known to be generally more intense over continental than maritime regions, with larger updraft velocities, graupel, and lightning production. This study explores the transition from maritime to continental convection by comparing the trends in Tropical Rainfall Measuring Mission (TRMM) radar and microwave (37 and 85 GHz) observations over islands of increasing size to those simulated by a cloud-resolving model. The observed storms were essentially maritime over islands of <100 km2 and continental over islands >10 000 km2, with a gradual transition in between. Equivalent radar and microwave quantities were simulated from cloud-resolving runs of the Weather Research and Forecasting model via offline radiation codes. The model configuration was idealized, with islands represented by regions of uniform surface heat flux without orography, using a range of initial sounding conditions without strong horizontal winds or aerosols. Simulated storm strength varied with initial sounding, as expected, but also increased sharply with island size in a manner similar to observations. Stronger simulated storms were associated with higher concentrations of large hydrometeors. Although biases varied with different ice microphysical schemes, the trend was similar for all three schemes tested and was also seen in 2D and 3D model configurations. The successful reproduction of the trend with such idealized forcing supports previous suggestions that mesoscale variation in surface heating—rather than any difference in humidity, aerosol, or other aspects of the atmospheric state—is the main reason that convection is more intense over continents and large islands than over oceans. Some dynamical storm aspects, notably the peak rainfall and minimum surface pressure low, were more sensitive to surface forcing than to the atmospheric sounding or ice scheme. Large hydrometeor concentrations and simulated microwave and radar signatures, however, were at least as sensitive to initial humidity levels as to surface forcing and were more sensitive to the ice scheme. Issues with running the TRMM simulator on 2D simulations are discussed, but they appear to be less serious than sensitivities to model microphysics, which were similar in 2D and 3D. This supports the further use of 2D simulations to economically explore modeling uncertainties.

On the response of the Antarctic Circumpolar Current transport to climate change in coupled climate models

The CMIP3 (IPCC AR4) models show a consistent intensification and poleward shift of the westerly winds over the Southern Ocean during the 21st century. However, the responses of the Antarctic Circumpolar Currents (ACC) show great diversity in these models, with many even showing reductions in transport. To obtain some understanding of diverse responses in the ACC transport, we investigate both external atmospheric and internal oceanic processes that control the ACC transport responses in these models. While the strengthened westerlies act to increase the tilt of isopycnal surfaces and hence the ACC transport through Ekman pumping effects, the associated changes in buoyancy forcing generally tend to reduce the surface meridional density gradient. The steepening of isopycnal surfaces induced by increased wind forcing leads to enhanced (parameterized) eddy-induced transports that act to reduce the isopycnal slopes. There is also considerable narrowing of the ACC that tends to reduce the ACC transport, caused mainly by the poleward shifts of the subtropical gyres and to a lesser extent by the equatorward expansions of the subpolar gyres in some models. If the combined effect of these retarding processes is larger than that of enhanced Ekman pumping, the ACC transport will be reduced. In addition, the effect of Ekman pumping on the ACC is reduced in weakly stratified models. These findings give insight into the reliability of IPCC-class model predictions of the Southern Ocean circulation, and into the observed decadal-scale steady ACC transport.

Boundary layer dynamics over London, UK, as observed using Doppler lidar during REPARTEE-II

Urban boundary layers (UBLs) can be highly complex due to the heterogeneous roughness and heating of the surface, particularly at night. Due to a general lack of observations, it is not clear whether canonical models of boundary layer mixing are appropriate in modelling air quality in urban areas. This paper reports Doppler lidar observations of turbulence profiles in the centre of London, UK, as part of the second REPARTEE campaign in autumn 2007. Lidar-measured standard deviation of vertical velocity averaged over 30 min intervals generally compared well with in situ sonic anemometer measurements at 190 m on the BT telecommunications Tower. During calm, nocturnal periods, the lidar underestimated turbulent mixing due mainly to limited sampling rate. Mixing height derived from the turbulence, and aerosol layer height from the backscatter profiles, showed similar diurnal cycles ranging from c. 300 to 800 m, increasing to c. 200 to 850 m under clear skies. The aerosol layer height was sometimes significantly different to the mixing height, particularly at night under clear skies. For convective and neutral cases, the scaled turbulence profiles resembled canonical results; this was less clear for the stable case. Lidar observations clearly showed enhanced mixing beneath stratocumulus clouds reaching down on occasion to approximately half daytime boundary layer depth. On one occasion the nocturnal turbulent structure was consistent with a nocturnal jet, suggesting a stable layer. Given the general agreement between observations and canonical turbulence profiles, mixing timescales were calculated for passive scalars released at street level to reach the BT Tower using existing models of turbulent mixing. It was estimated to take c. 10 min to diffuse up to 190 m, rising to between 20 and 50 min at night, depending on stability. Determination of mixing timescales is important when comparing to physico-chemical processes acting on pollutant species measured simultaneously at both the ground and at the BT Tower during the campaign. From the 3 week autumnal data-set there is evidence for occasional stable layers in central London, effectively decoupling surface emissions from air aloft.

Sea-ice anomalies in the Sea of Okhotsk and the relationship with storm tracks in the Northern Hemisphere during winter

Previous studies have shown that sea-ice in the Sea of Okhotsk can be affected by local storms; in turn, the resultant sea-ice changes can affect the downstream development of storm tracks in the Pacific and possibly dampen a pre-existing North Atlantic Oscillation (NAO) signal in late winter. In this paper, a storm tracking algorithm was applied to the six hourly horizontal winds from the National Centers for Environmental Prediction (NCEP) reanalysis data from 1978(9) to 2007 and output from the atmospheric general circulation model (AGCM) ECHAM5 forced by sea-ice anomalies in the Sea of Okhotsk. The life cycle response of storms to sea-ice anomalies is investigated using various aspects of storm activity—cyclone genesis, lysis, intensity and track density. Results show that, for enhanced positive sea-ice concentrations in the Sea of Okhotsk, there is a decrease in secondary cyclogenesis, a westward shift in cyclolysis and changes in the subtropical jet are seen in the North Pacific. In the Atlantic, a pattern resembling the negative phase of the NAO is observed. This pattern is confirmed by the AGCM ECHAM5 experiments driven with above normal sea-ice anomalies in the Sea of Okhotsk

A comparison of extra-tropical cyclones in recent re-analyses; ERA-Interim, NASA MERRA, NCEP CFSR and JRA-25

Extra-tropical cyclones are identified and compared using data from four recent re-analyses for the winter periods in both hemispheres. Results show the largest differences occur between the older lower resolution JRA25 re-analysis when compared with the newer high resolution re-analyses, in particular in the Southern Hemisphere (SH). Spatial differences between the newest re-analyses are small in both hemispheres and generally not significant except some common regions associated with cyclogenesis close to orography. Intensities are generally related to spatial resolution except NASA-MERRA which has larger intensities for several different measures. Matching storms between re-analyses shows the number matched between ERA-Interim and the other re-analyses are similar in the Northern Hemisphere (NH). In the SH the number matched between JRA25 and ERA-Interim is lower than in the NH, but for NASA-MERRA and NCEP-CFSR the number matched is similar to the NH. The mean separation of the identically same cyclones is typically less than 20 geodesic in both hemispheres for the latest re-analyses, whereas JRA25 compared with the other re-analyses has a broader distribution in the SH indicating greater uncertainty. The instantaneous intensity differences for matched storms shows narrow distributions for pressure while for winds and vorticity the distributions are much broader indicating larger uncertainty typical of smaller scale fields. Composite cyclone diagnostics show that cyclones are very similar between the re-analyses, with differences being related to the intensities, consistent with the intensity results. Overall, results show NH cyclones correspond well between re-analyses, with a significant improvement in the SH for the latest re-analyses, indicating a convergence between re-analyses for cyclone properties.

Northern Hemisphere extratropical cyclones in a warming climate in the HiGEM high-resolution climate model

Changes to the Northern Hemisphere winter (December, January and February) extratropical storm tracks and cyclones in a warming climate are investigated. Two idealised climate change experiments with HiGEM1.1, a doubled CO2 and a quadrupled CO2 experiment, are compared against a present day control run. An objective feature tracking method is used and a focus given to regional changes. The climatology of extratropical storm tracks from the control run is shown to be in good agreement with ERA-40, while the frequency distribution of cyclone intensity also compares well. In both simulations the mean climate changes are generally consistent with the simulations of the IPCC AR4 models, with a strongly enhanced surface warming at the winter pole and the reduced lower tropospheric warming over the North Atlantic Ocean associated with the slowdown of the Meridional Overturning Circulation. The circulation changes in the North Atlantic are different between the two idealised simulations with different CO2 forcings. In the North Atlantic the storm tracks are influenced by the slowdown of the MOC, the enhanced surface polar warming, and the enhanced upper tropical troposphere warming, giving a north eastward shift of the storm tracks in the 2XCO2 experiment, but no shift in the 4XCO2 experiment. Over the Pacific, in the 2XCO2 experiment, changes in the mean climate are associated with local temperature changes, while in the 4XCO2 experiment the changes in the Pacific are impacted by the weakened tropical circulation. The storm track changes are consistent with the shifts in the zonal wind. Total cyclone numbers are found to decrease over the Northern Hemisphere with increasing CO2 forcing. Changes in cyclone intensity are found using 850hPa vorticity, mean sea level pressure, and 850hPa winds. The intensity of the Northern Hemisphere cyclones is found to decrease relative to the control.

Impact of increasing resolution and a warmer climate on extreme weather from NH extra-tropical cyclones

The effect of a warmer climate on the properties of extra-tropical cyclones is investigated using simulations of the ECHAM5 global climate model at resolutions of T213 (60 km) and T319 (40 km). Two periods representative of the end of the 20th and 21st centuries are investigated using the IPCC A1B scenario. The focus of the paper is on precipitation for the NH summer and winter seasons, however results from vorticity and winds are also presented. Similar number of events are identified at both resolutions. There are, however, a greater number of extreme precipitation events in the higher reso- lution run. The difference between maximum intensity distributions are shown to be statistically significant using a Kolmogorov-Smirnov test. A Generalised Pareto Distribution is used to analyse changes in extreme precipitation and wind events. In both resolutions, there is an increase in the number of ex- treme precipitation events in a warmer climate for all seasons, together with a reduction in return period. This is not associated with any increased verti- cal velocity, or with any increase in wind intensity in the winter and spring. However, there is an increase in wind extremes in the summer and autumn associated with tropical cyclones migrating into the extra-tropics.

Ultra-distal tephra deposits from super-eruptions: examples from Toba, Indonesia and Taupo Volcanic Zone, New Zealand

Voluminous rhyolitic eruptions from Toba, Indonesia, and Taupo Volcanic Zone (TVZ), New Zealand, have dispersed volcanic ash over vast areas in the late Quaternary. The ~74 ka Youngest Toba Tuff (YTT) eruption deposited ash over the Bay of Bengal and the Indian subcontinent to the west. The ~340 ka Whakamaru eruption (TVZ) deposited the widespread Rangitawa Tephra, dominantly to the southeast (in addition to occurrences northwest of vent), extending across the landmass of New Zealand, and the South Pacific Ocean and Tasman Sea, with distal terrestrial exposures on the Chatham Islands. These super-eruptions involved ~2500 km^3 and ~1500 km3 of magma (dense-rock equivalent; DRE), respectively. Ultra-distal terrestrial exposures of YTT at two localities in India, Middle Son Valley, Madhya Pradesh, and Jurreru River Valley, Andhra Pradesh, at distances of >2000 km from the source caldera, show a basal ‘primary’ ashfall unit ~4 cm thick, although deposits containing reworked ash are up to ~3 m in total thickness. Exposures of Rangitawa Tephra on the Chatham Islands, >900 km from the source caldera, are ~15-30 cm thick. At more proximal localities (~200 km from source), Rangitawa Tephra is ~55-70 cm thick and characterized by a crystal-rich basal layer and normal grading. Both distal tephra deposits are characterized by very-fine ash (with high PM10 fractions) and are crystal-poor. Glass chemistry, stratigraphy and grain-size data for these distal tephra deposits are presented with comparisons of their correlation, dispersal and preservation. Using field observations, ash transport and deposition were modeled for both eruptions using a semi-analytical model (HAZMAP), with assumptions concerning average wind direction and strength during eruption, column shape and vent size. Model outputs provide new insights into eruption dynamics and better estimates of eruption volumes associ- ated with tephra fallout. Modeling based on observed YTT distal tephra thicknesses indicate a relatively low (<40 km high), very turbulent eruption column, consistent with deposition from a co-ignimbrite cloud extending over a broad region. Similarly, the Whakamaru eruption was modeled as producing a predominantly Plinian column (~45 km high), with dispersal to the southeast by strong prevailing winds. Significant ash fallout of the main dispersal direction, to the northwest of source, cannot be replicated in this modeling. The widespread dispersal of large volumes of fine ash from both eruptions may have had global environmental consequences, acutely affecting areas up to thousands of kilometers from vent.

The equatorial response to higher latitude forcing

The classic view, following Charney and Webster and Holton, is that significant midlatitude forcing of the Tropics can be expected only in regions with westerly winds in the upper troposphere because it is only in these regions that stationary Rossby waves will be able to propagate toward the equator. Here it is shown that higherlatitude forcing can project directly onto equatorial waves and give a significant tropical response in both easterly and westerly tropical flow. The equatorial response to higher-latitude forcing is considered in the context of a dry atmosphere and a localized higher-latitude forcing with eastward or westward phase speed. Previous ideas of the Doppler shifting of equatorial waves by zonal flows are extended to include consideration of a forcing involving a range of zonal wavenumbers. A Gill-type model suggests that there can be significant forcing of equatorial waves by either vorticity forcing or heating in higher latitudes. In agreement with the theory, the Kelvin wave response to eastward forcing is peaked at high frequencies/short periods but reduces only slowly with decreasing frequency. Primitive-equation experiments confirm the strong equatorial response associated with a deep Kelvin wave for forcing in midlatitudes. The response is strongest in the Eastern Hemisphere with its equatorial, upper-tropospheric easterlies. The possible importance of this equatorial response in the organization of large-scale, deep tropical convection and the initiation of the Madden–Julian oscillation is discussed. The ability of westward forcing in higher latitudes to trigger Rossby–gravity and Rossby waves is found in the primitive-equation model to be significant but rather less robust. These wave signatures are clearest in the lower troposphere. For shorter periods the Rossby–gravity wave dominates, and for upper-tropospheric forcing, downward and eastward wave activity propagation is seen. Upper-tropospheric westerlies are found to enhance the response.

Atmospheric conservation properties in ERA-Interim

We study the global atmospheric budgets of mass, moisture, energy and angular momentum in the latest reanalysis from the European Centre for Medium-Range Weather Forecasts (ECMWF), ERA-Interim, for the period 1989–2008 and compare with ERA-40. Most of the measures we use indicate that the ERA-Interim reanalysis is superior in quality to ERA-40. In ERA-Interim the standard deviation of the monthly mean global dry mass of 0.7 kg m−2 (0.007%) is slightly worse than in ERA-40, and long time-scale variations in dry mass originate predominately in the surface pressure field. The divergent winds are improved in ERA-Interim: the global standard deviation of the time-averaged dry mass budget residual is 10 kg m−2 day−1 and the quality of the cross-equatorial mass fluxes is improved. The temporal variations in the global evaporation minus precipitation (E − P) are too large but the global moisture budget residual is 0.003 kg m−2 day−1 with a spatial standard deviation of 0.3 kg m−2 day−1. Both the E − P over ocean and P − E over land are about 15% larger than the 1.1 Tg s−1 transport of water from ocean to land. The top of atmosphere (TOA) net energy losses are improved, with a value of 1 W m−2, but the meridional gradient of the TOA net energy flux is smaller than that from the Clouds and the Earth's Radiant Energy System (CERES) data. At the surface the global energy losses are worse, with a value of 7 W m−2. Over land however, the energy loss is only 0.5 W m−2. The downwelling thermal radiation at the surface in ERA-Interim of 341 W m−2 is towards the higher end of previous estimates. The global mass-adjusted energy budget residual is 8 W m−2 with a spatial standard deviation of 11 W m−2, and the mass-adjusted atmospheric energy transport from low to high latitudes (the sum for the two hemispheres) is 9.5 PW