Large-scale context for the UK floods in summer 2007

This paper shows how the rainfall distribution over the UK, in the three major events on 13-15 June, 25 June and 20 July 2007, was related to troughs in the upper-level flow, and investigates the relationship of these features to a persistent large-scale flow pattern which extended around the northern hemisphere and its possible origins. Remote influences can be mediated by the propagation of large-scale atmospheric waves across the northern hemisphere and also by the origins of the air-masses that are wrapped into the developing weather systems delivering the rain to the UK. These dynamical influences are examined using analyses and forecasts produced by a range of atmospheric models.

Application of the Lighthill-Ford theory of spontaneous imbalance to clear-air turbulence forecasting

A new method of clear-air turbulence (CAT) forecasting based on the Lighthill–Ford theory of spontaneous imbalance and emission of inertia–gravity waves has been derived and applied on episodic and seasonal time scales. A scale analysis of this shallow-water theory for midlatitude synoptic-scale flows identifies advection of relative vorticity as the leading-order source term. Examination of leading- and second-order terms elucidates previous, more empirically inspired CAT forecast diagnostics. Application of the Lighthill–Ford theory to the Upper Mississippi and Ohio Valleys CAT outbreak of 9 March 2006 results in good agreement with pilot reports of turbulence. Application of Lighthill–Ford theory to CAT forecasting for the 3 November 2005–26 March 2006 period using 1-h forecasts of the Rapid Update Cycle (RUC) 2 1500 UTC model run leads to superior forecasts compared to the current operational version of the Graphical Turbulence Guidance (GTG1) algorithm, the most skillful operational CAT forecasting method in existence. The results suggest that major improvements in CAT forecasting could result if the methods presented herein become operational.

Decadal predictability of the Atlantic Ocean in a coupled GCM: forecast skill and optimal perturbations using Linear Inverse Modelling

The decadal predictability of three-dimensional Atlantic Ocean anomalies is examined in a coupled global climate model (HadCM3) using a Linear Inverse Modelling (LIM) approach. It is found that the evolution of temperature and salinity in the Atlantic, and the strength of the meridional overturning circulation (MOC), can be effectively described by a linear dynamical system forced by white noise. The forecasts produced using this linear model are more skillful than other reference forecasts for several decades. Furthermore, significant non-normal amplification is found under several different norms. The regions from which this growth occurs are found to be fairly shallow and located in the far North Atlantic. Initially, anomalies in the Nordic Seas impact the MOC, and the anomalies then grow to fill the entire Atlantic basin, especially at depth, over one to three decades. It is found that the structure of the optimal initial condition for amplification is sensitive to the norm employed, but the initial growth seems to be dominated by MOC-related basin scale changes, irrespective of the choice of norm. The consistent identification of the far North Atlantic as the most sensitive region for small perturbations suggests that additional observations in this region would be optimal for constraining decadal climate predictions.

The impact of the stratosphere on the troposphere during the southern hemisphere stratospheric sudden warming, September 2002

Recent research has established that a small but statistically significant link exists between the stratosphere and the troposphere in the northern hemisphere extratropics. In this paper it is shown that a similar link exists between the stratosphere and troposphere during the unprecedented September 2002 sudden warming in the southern hemisphere. Two ensemble forecasts of the stratospheric sudden warming are run which have different stratospheric initial conditions and identical tropospheric initial conditions. Stratospheric initial conditions have an impact on the tropospheric flow at the peak of the major warming (5 days into the run) and on longer time-scales (18 days into the run). The character of this influence is a localized, equatorward shift of the tropospheric storm track. The averaged impact of the change in the position of the storm-track maps strongly onto the Southern Annular Mode structure, but does not have an annular character.

Stratospheric circulation in seasonal forecasting models: implications for seasonal prediction

Accurate seasonal forecasts rely on the presence of low frequency, predictable signals in the climate system which have a sufficiently well understood and significant impact on the atmospheric circulation. In the Northern European region, signals associated with seasonal scale variability such as ENSO, North Atlantic SST anomalies and the North Atlantic Oscillation have not yet proven sufficient to enable satisfactorily skilful dynamical seasonal forecasts. The winter-time circulations of the stratosphere and troposphere are highly coupled. It is therefore possible that additional seasonal forecasting skill may be gained by including a realistic stratosphere in models. In this study we assess the ability of five seasonal forecasting models to simulate the Northern Hemisphere extra-tropical winter-time stratospheric circulation. Our results show that all of the models have a polar night jet which is too weak and displaced southward compared to re-analysis data. It is shown that the models underestimate the number, magnitude and duration of periods of anomalous stratospheric circulation. Despite the poor representation of the general circulation of the stratosphere, the results indicate that there may be a detectable tropospheric response following anomalous circulation events in the stratosphere. However, the models fail to exhibit any predictability in their forecasts. These results highlight some of the deficiencies of current seasonal forecasting models with a poorly resolved stratosphere. The combination of these results with other recent studies which show a tropospheric response to stratospheric variability, demonstrates a real prospect for improving the skill of seasonal forecasts.

Variability in clear-sky longwave radiative cooling of the atmosphere

The longwave radiative cooling of the clear-sky atmosphere (Q(LWc)) is a crucial component of the global hydrological cycle and is composed of the clear-sky outgoing longwave radiation to space (OLRc) and the net downward minus upward clear-sky longwave radiation to the surface (SNLc). Estimates of QLWc from reanalyses and observations are presented for the period 1979-2004. Compared to other reanalyses data sets, the European Centre for Medium-range Weather Forecasts 40-year reanalysis (ERA40) produces the largest Q(LWc) over the tropical oceans (217 W m(-2)), explained by the least negative SNLc. On the basis of comparisons with data derived from satellite measurements, ERA40 provides the most realistic QLWc climatology over the tropical oceans but exhibits a spurious interannual variability for column integrated water vapor (CWV) and SNLc. Interannual monthly anomalies of QLWc are broadly consistent between data sets with large increases during the warm El Nino events. Since relative humidity ( RH) errors applying throughout the troposphere result in compensating effects on the cooling to space and to the surface, they exert only a marginal effect on QLWc. An observed increase in CWV with surface temperature of 3 kg m(-2) K-1 over the tropical oceans is important in explaining a positive relationship between QLWc and surface temperature, in particular over ascending regimes; over tropical ocean descending regions this relationship ranges from 3.6 to 4.6 +/- 0.4 W m(-2) K-1 for the data sets considered, consistent with idealized sensitivity tests in which tropospheric warming is applied and RH is held constant and implying an increase in precipitation with warming.

Regional Differences in the Prediction of Extratropical Cyclones by the ECMWF Ensemble Prediction System

A regional study of the prediction of extratropical cyclones by the European Centre for Medium-Range Weather Forecasts (ECMWF) Ensemble Prediction System (EPS) has been performed. An objective feature-tracking method has been used to identify and track the cyclones along the forecast trajectories. Forecast error statistics have then been produced for the position, intensity, and propagation speed of the storms. In previous work, data limitations meant it was only possible to present the diagnostics for the entire Northern Hemisphere (NH) or Southern Hemisphere. A larger data sample has allowed the diagnostics to be computed separately for smaller regions around the globe and has made it possible to explore the regional differences in the prediction of storms by the EPS. Results show that in the NH there is a larger ensemble mean error in the position of storms over the Atlantic Ocean. Further analysis revealed that this is mainly due to errors in the prediction of storm propagation speed rather than in direction. Forecast storms propagate too slowly in all regions, but the bias is about 2 times as large in the NH Atlantic region. The results show that storm intensity is generally overpredicted over the ocean and underpredicted over the land and that the absolute error in intensity is larger over the ocean than over the land. In the NH, large errors occur in the prediction of the intensity of storms that originate as tropical cyclones but then move into the extratropics. The ensemble is underdispersive for the intensity of cyclones (i.e., the spread is smaller than the mean error) in all regions. The spatial patterns of the ensemble mean error and ensemble spread are very different for the intensity of cyclones. Spatial distributions of the ensemble mean error suggest that large errors occur during the growth phase of storm development, but this is not indicated by the spatial distributions of the ensemble spread. In the NH there are further differences. First, the large errors in the prediction of the intensity of cyclones that originate in the tropics are not indicated by the spread. Second, the ensemble mean error is larger over the Pacific Ocean than over the Atlantic, whereas the opposite is true for the spread. The use of a storm-tracking approach, to both weather forecasters and developers of forecast systems, is also discussed.

Variational data assimilation for parameter estimation: application to a simple morphodynamic model

Data assimilation is a sophisticated mathematical technique for combining observational data with model predictions to produce state and parameter estimates that most accurately approximate the current and future states of the true system. The technique is commonly used in atmospheric and oceanic modelling, combining empirical observations with model predictions to produce more accurate and well-calibrated forecasts. Here, we consider a novel application within a coastal environment and describe how the method can also be used to deliver improved estimates of uncertain morphodynamic model parameters. This is achieved using a technique known as state augmentation. Earlier applications of state augmentation have typically employed the 4D-Var, Kalman filter or ensemble Kalman filter assimilation schemes. Our new method is based on a computationally inexpensive 3D-Var scheme, where the specification of the error covariance matrices is crucial for success. A simple 1D model of bed-form propagation is used to demonstrate the method. The scheme is capable of recovering near-perfect parameter values and, therefore, improves the capability of our model to predict future bathymetry. Such positive results suggest the potential for application to more complex morphodynamic models.

The accuracy of Doppler radar wind retrievals using insects as targets

Insect returns from the UK's Doppler weather radars were collected in the summers of 2007 and 2008, to ascertain their usefulness in providing information about boundary layer winds. Such observations could be assimilated into numerical weather prediction models to improve forecasts of convective showers before precipitation begins. Significant numbers of insect returns were observed during daylight hours on a number of days through this period, when they were detected at up to 30 km range from the radars, and up to 2 km above sea level. The range of detectable insect returns was found to vary with time of year and temperature. There was also a very weak correlation with wind speed and direction. Use of a dual-polarized radar revealed that the insects did not orient themselves at random, but showed distinct evidence of common orientation on several days, sometimes at an angle to their direction of travel. Observation minus model background residuals of wind profiles showed greater bias and standard deviation than that of other wind measurement types, which may be due to the insects' headings/airspeeds and to imperfect data extraction. The method used here, similar to the Met Office's procedure for extracting precipitation returns, requires further development as clutter contamination remained one of the largest error contributors. Wind observations derived from the insect returns would then be useful for data assimilation applications.

The ERA-40 re-analysis

ERA-40 is a re-analysis of meteorological observations from September 1957 to August 2002 produced by the European Centre for Medium-Range Weather Forecasts (ECMWF) in collaboration with many institutions. The observing system changed considerably over this re-analysis period, with assimilable data provided by a succession of satellite-borne instruments from the 1970s onwards, supplemented by increasing numbers of observations from aircraft, ocean-buoys and other surface platforms, but with a declining number of radiosonde ascents since the late 1980s. The observations used in ERA-40 were accumulated from many sources. The first part of this paper describes the data acquisition and the principal changes in data type and coverage over the period. It also describes the data assimilation system used for ERA-40. This benefited from many of the changes introduced into operational forecasting since the mid-1990s, when the systems used for the 15-year ECMWF re-analysis (ERA-15) and the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) re-analysis were implemented. Several of the improvements are discussed. General aspects of the production of the analyses are also summarized. A number of results indicative of the overall performance of the data assimilation system, and implicitly of the observing system, are presented and discussed. The comparison of background (short-range) forecasts and analyses with observations, the consistency of the global mass budget, the magnitude of differences between analysis and background fields and the accuracy of medium-range forecasts run from the ERA-40 analyses are illustrated. Several results demonstrate the marked improvement that was made to the observing system for the southern hemisphere in the 1970s, particularly towards the end of the decade. In contrast, the synoptic quality of the analysis for the northern hemisphere is sufficient to provide forecasts that remain skilful well into the medium range for all years. Two particular problems are also examined: excessive precipitation over tropical oceans and a too strong Brewer-Dobson circulation, both of which are pronounced in later years. Several other aspects of the quality of the re-analyses revealed by monitoring and validation studies are summarized. Expectations that the second-generation ERA-40 re-analysis would provide products that are better than those from the firstgeneration ERA-15 and NCEP/NCAR re-analyses are found to have been met in most cases. © Royal Meteorological Society, 2005. The contributions of N. A. Rayner and R. W. Saunders are Crown copyright.

Water vapor variability in the tropics and its links to dynamics and precipitation

The distribution and variability of water vapor and its links with radiative cooling and latent heating via precipitation are crucial to understanding feedbacks and processes operating within the climate system. Column-integrated water vapor (CWV) and additional variables from the European Centre for Medium-Range Weather Forecasts (ECMWF) 40-year reanalysis (ERA40) are utilized to quantify the spatial and temporal variability in tropical water vapor over the period 1979–2001. The moisture variability is partitioned between dynamical and thermodynamic influences and compared with variations in precipitation provided by the Climate Prediction Center Merged Analysis of Precipitation (CMAP) and the Global Precipitation Climatology Project (GPCP). The spatial distribution of CWV is strongly determined by thermodynamic constraints. Spatial variability in CWV is dominated by changes in the large-scale dynamics, in particular associated with the El Niño–Southern Oscillation (ENSO). Trends in CWV are also dominated by dynamics rather than thermodynamics over the period considered. However, increases in CWV associated with changes in temperature are significant over the equatorial east Pacific when analyzing interannual variability and over the north and northwest Pacific when analyzing trends. Significant positive trends in CWV tend to predominate over the oceans while negative trends in CWV are found over equatorial Africa and Brazil. Links between changes in CWV and vertical motion fields are identified over these regions and also the equatorial Atlantic. However, trends in precipitation are generally incoherent and show little association with the CWV trends. This may in part reflect the inadequacies of the precipitation data sets and reanalysis products when analyzing decadal variability. Though the dynamic component of CWV is a major factor in determining precipitation variability in the tropics, in some regions/seasons the thermodynamic component cancels its effect on precipitation variability.

The impact of targeted observations made during the Greenland Flow Distortion Experiment

The impact of targeted sonde observations on the 1-3 day forecasts for northern Europe is evaluated using the Met Office four-dimensional variational data assimilation scheme and a 24 km gridlength limited-area version of the Unified Model (MetUM). The targeted observations were carried out during February and March 2007 as part of the Greenland Flow Distortion Experiment, using a research aircraft based in Iceland. Sensitive area predictions using either total energy singular vectors or an ensemble transform Kalman filter were used to predict where additional observations should be made to reduce errors in the initial conditions of forecasts for northern Europe. Targeted sonde data was assimilated operationally into the MetUM. Hindcasts show that the impact of the sondes was mixed. Only two out of the five cases showed clear forecast improvement; the maximum forecast improvement seen over the verifying region was approximately 5% of the forecast error 24 hours into the forecast. These two cases are presented in more detail: in the first the improvement propagates into the verification region with a developing polar low; and in the second the improvement is associated with an upper-level trough. The impact of cycling targeted data in the background of the forecast (including the memory of previous targeted observations) is investigated. This is shown to cause a greater forecast impact, but does not necessarily lead to a greater forecast improvement. Finally, the robustness of the results is assessed using a small ensemble of forecasts.

Voronoi, Delaunay, and Block-Structured Mesh Refinement for Solution of the Shallow-Water Equations on the Sphere

Alternative meshes of the sphere and adaptive mesh refinement could be immensely beneficial for weather and climate forecasts, but it is not clear how mesh refinement should be achieved. A finite-volume model that solves the shallow-water equations on any mesh of the surface of the sphere is presented. The accuracy and cost effectiveness of four quasi-uniform meshes of the sphere are compared: a cubed sphere, reduced latitude–longitude, hexagonal–icosahedral, and triangular–icosahedral. On some standard shallow-water tests, the hexagonal–icosahedral mesh performs best and the reduced latitude–longitude mesh performs well only when the flow is aligned with the mesh. The inclusion of a refined mesh over a disc-shaped region is achieved using either gradual Delaunay, gradual Voronoi, or abrupt 2:1 block-structured refinement. These refined regions can actually degrade global accuracy, presumably because of changes in wave dispersion where the mesh is highly nonuniform. However, using gradual refinement to resolve a mountain in an otherwise coarse mesh can improve accuracy for the same cost. The model prognostic variables are height and momentum collocated at cell centers, and (to remove grid-scale oscillations of the A grid) the mass flux between cells is advanced from the old momentum using the momentum equation. Quadratic and upwind biased cubic differencing methods are used as explicit corrections to a fast implicit solution that uses linear differencing.

Surface wave processes in air-sea interaction

We review briefly recent progress on understanding the role of surface waves on the marine atmospheric boundary layer and the ocean mixed layer and give a global perspective on these processes by analysing ERA-40 data. Ocean surface waves interact with the marine atmospheric boundary layer in two broad regimes: (i) the conventional wind-driven wave regime, when fast winds blow over slower moving waves, and (ii) a wave-driven wind regime when long wavelength swell propagates under low winds, and generates a wave-driven jet in the lower part of the marine boundary layer. Analysis of ERA-40 data indicates that the wave-driven wind regime is as prevalent as the conventional wind-driven regime. Ocean surface waves also change profoundly mixing in the ocean mixed layer through generation of Langmuir circulation. Results from large-eddy simulation are used here to develop a scaling for the resulting Langmuir turbulence, which is a necessary step in developing a parametrization of the process. ERA-40 data is then used to show that the Langmuir regime is the predominant regime over much of the global ocean, providing a compelling motivation for parameterising this process in ocean general circulation models.

Variability in the summer season hydrological cycle over the Atlantic-Europe region 1979-2007

Variability in aspects of the hydrological cycle over the Europe-Atlantic region during the summer season is analysed for the period 1979-2007, using observational estimates, reanalyses and climate model simulations. Warming and moistening trends are evident in observations and models although decadal changes in water vapour are not well represented by reanalyses, including the new European Centre for Medium Range Weather Forecasts (ECMWF) Interim reanalysis. Over the north Atlantic and northern Europe, observed water vapour trends are close to that expected from the temperature trends and Clausius-Clapeyron equation (7% K-1), larger than the model simulations. Precipitation over Europe is dominated by large-scale dynamics with positive phases of the North Atlantic Oscillation coinciding with drier conditions over north Europe and wetter conditions over the Mediterranean region. Evaporation trends over Europe are positive in reanalyses and models, especially for the Mediterranean region (1-3% per decade in reanalyses and climate models). Over the north Atlantic, declining precipitation combined with increased moisture contributed to an apparent rise in water vapour residence time. Maximum precipitation minus evaporation over the north Atlantic occurred during summer 1991, declining thereafter.