A singular vector perspective of 4D-Var: Filtering and interpolation

Four-dimensional variational data assimilation (4D-Var) combines the information from a time sequence of observations with the model dynamics and a background state to produce an analysis. In this paper, a new mathematical insight into the behaviour of 4D-Var is gained from an extension of concepts that are used to assess the qualitative information content of observations in satellite retrievals. It is shown that the 4D-Var analysis increments can be written as a linear combination of the singular vectors of a matrix which is a function of both the observational and the forecast model systems. This formulation is used to consider the filtering and interpolating aspects of 4D-Var using idealized case-studies based on a simple model of baroclinic instability. The results of the 4D-Var case-studies exhibit the reconstruction of the state in unobserved regions as a consequence of the interpolation of observations through time. The results also exhibit the filtering of components with small spatial scales that correspond to noise, and the filtering of structures in unobserved regions. The singular vector perspective gives a very clear view of this filtering and interpolating by the 4D-Var algorithm and shows that the appropriate specification of the a priori statistics is vital to extract the largest possible amount of useful information from the observations. Copyright © 2005 Royal Meteorological Society

Moist singular vectors and the predictability of some high impact European cyclones

The ECMWF full-physics and dry singular vector (SV) packages, using a dry energy norm and a 1-day optimization time, are applied to four high impact European cyclones of recent years that were almost universally badly forecast in the short range. It is shown that these full-physics SVs are much more relevant to severe cyclonic development than those based on dry dynamics plus boundary layer alone. The crucial extra ingredient is the representation of large-scale latent heat release. The severe winter storms all have a long, nearly straight region of high baroclinicity stretching across the Atlantic towards Europe, with a tongue of very high moisture content on its equatorward flank. In each case some of the final-time top SV structures pick out the region of the actual storm. The initial structures were generally located in the mid- to low troposphere. Forecasts based on initial conditions perturbed by moist SVs with opposite signs and various amplitudes show the range of possible 1-day outcomes for reasonable magnitudes of forecast error. In each case one of the perturbation structures gave a forecast very much closer to the actual storm than the control forecast. Deductions are made about the predictability of high-impact extratropical cyclone events. Implications are drawn for the short-range forecast problem and suggestions made for one practicable way to approach short-range ensemble forecasting. Copyright © 2005 Royal Meteorological Society.

The influence of physical processes on extratropical singular vectors

An investigation is made of the impact of a full linearized physical (moist) parameterization package on extratropical singular vectors (SVs) using the ECMWF integrated forecasting system (IFS). Comparison is made for one particular period with a dry physical package including only vertical diffusion and surface drag. The crucial extra ingredient in the full package is found to be the large-scale latent heat release. Consistent with basic theory, its inclusion results in a shift to smaller horizontal scales and enhanced growth for the SVs. Whereas, for the dry SVs, T42 resolution is sufficient, the moist SVs require T63 to resolve their structure and growth. A 24-h optimization time appears to be appropriate for the moist SVs because of the larger growth of moist SVs compared with dry SVs. Like dry SVs, moist SVs tend to occur in regions of high baroclinicity, but their location is also influenced by the availability of moisture. The most rapidly growing SVs appear to enhance or reduce large-scale rain in regions ahead of major cold fronts. The enhancement occurs in and ahead of a cyclonic perturbation and the reduction in and ahead of an anticyclonic perturbation. Most of the moist SVs for this situation are slightly modified versions of the dry SVs. However, some occur in new locations and have particularly confined structures. The most rapidly growing SV is shown to exhibit quite linear behavior in the nonlinear model as it grows from 0.5 to 12 hPa in 1 day. For 5 times this amplitude the structure is similar but the growth is about half as the perturbation damps a potential vorticity (PV) trough or produces a cutoff, depending on its sign.

Changes in the global sulfate burden due to perturbations in global CO2 concentrations

A large ensemble of general circulation model (GCM) integrations coupled to a fully interactive sulfur cycle scheme were run on the climateprediction.net platform to investigate the uncertainty in the climate response to sulfate aerosol and carbon dioxide (CO2) forcing. The sulfate burden within the model (and the atmosphere) depends on the balance between formation processes and deposition (wet and dry). The wet removal processes for sulfate aerosol are much faster than dry removal and so any changes in atmospheric circulation, cloud cover, and precipitation will feed back on the sulfate burden. When CO2 is doubled in the Hadley Centre Slab Ocean Model (HadSM3), global mean precipitation increased by 5%; however, the global mean sulfate burden increased by 10%. Despite the global mean increase in precipitation, there were large areas of the model showing decreases in precipitation (and cloud cover) in the Northern Hemisphere during June–August, which reduced wet deposition and allowed the sulfate burden to increase. Further experiments were also undertaken with and without doubling CO2 while including a future anthropogenic sulfur emissions scenario. Doubling CO2 further enhanced the increases in sulfate burden associated with increased anthropogenic sulfur emissions as observed in the doubled CO2-only experiment. The implications are that the climate response to doubling CO2 can influence the amount of sulfate within the atmosphere and, despite increases in global mean precipitation, may act to increase it.

Coriolis perturbations in the infrared spectrum of ethylene

Rotational structure has been resolved and analyzed in two of the infrared‐active perpendicular bands of C2H4 vapor: the Type b fundamental band, ν10, at 826 cm—1, and the Type c fundamental band, ν7, at 949 cm—1. Many of the individual PP and RR branch lines have been observed. The analysis has been confined to values of the quantum number K≥3, for which energy levels ethylene shows no detectable deviations from a symmetric‐top rotational structure. The analysis reveals a Coriolis interaction between ν7 and ν10, and between ν4 and ν10, and values of the Coriolis constants ζ7,10z and ζ4,10y are obtained; these are related to normal coordinate calculations for the appropriate symmetry species, and force constants are derived to fit the observed zeta constants. The band center of ν10 has been revised from the original figure of 810 cm—1 to the new value, 826 cm—1, and the inactive frequency ν4 is estimated to lie at 1023±3 cm—1, in good agreement with the previous estimate of 1027 cm—1. The change in the value of ν10 leads to a suggested change in the value of the Raman‐active fundamental ν6 from 1236 to 1222 cm—1. New combination bands have been observed at 2174 cm—1, assigned as ν3+ν10; and at 2252 cm—1, assigned as ν4+ν6; also rotational structure has been resolved and analyzed in the ν6+ν10 band at 2048 cm—1. The new data obtained for the C2H4 molecule are summarized in Table XII, with all of the other data presently available on the vibrational and rotational constants.

Coriolis perturbations in the infrared spectrum of allene

High resolution infrared spectra of the ν9 and ν10 perpendicular fundamentals of the allene molecule are reported, in which the J structure in the sub-bands has been partially resolved. Analysis of the latter shows that the vibrational origin ν9 = 999 cm−1, some 35 cm−1 below previous assignments. The pronounced asymmetry in the intensity distribution of the rotational structure which this assignment implies is shown to be expected theoretically, due to the Coriolis perturbations involved, and it is interpreted in terms of the sign and magnitude of the ratio of the dipole moment derivatives in the two fundamentals. The results of this analysis are shown to be in good agreement with observations on allene-1.1-d2, where similar intensity perturbations are observed, and with an independent analysis of the ν8 band of allene-h4. The A rotational constant of allene-h4 is found to have the value 4.82 ± 0.01 cm−1, and for the molecular geometry we obtain r(CH) = 1.084 A, r(CC) = 1.308 A, and HCH = 118.4°. A partial analysis of the rotational structure of the hot bands (ν9 + ν11 − ν11) and (ν10 + ν11 − ν11) is presented; these provide an example of a strong Coriolis interaction between nearly degenerate A1A2 and B1B2 pairs of vibrational levels. Some localized rotational perturbations in the ν9 and ν10 fundamentals are also noted, and their possible interpretations are discussed.

Coriolis perturbations in the infrared spectrum of diborane

High‐resolution infrared spectra of B2H6 vapor are reported. The sample was prepared from the naturally occurring 11B☒10B isotopic mixture. The rotational structure of the infrared bands has been analysed for Coriolis perturbations due to rotation about the axis of least moment of inertia (the B⋅⋅⋅B axis). The following results have been obtained: (a) interaction between the Type A fundamental ν18 and the inactive fundamental ν5 has been observed, thus confirming the assignment of ν5 at 833 cm—1, giving ∣ ζ5,18Z ∣=0.55±0.05; (b) interaction observed between the Type A combination band (ν10+ν12) at 1283 cm—1 and the inactive combination (ν10+ν7) gives an estimate of the unobserved fundamental ν7 as 850±30 cm—1, and an estimate of ∣ ζ7,12Z ∣=0.6±0.1; (c) the absence of any observed perturbation of the Type C fundamental ν14 at 973 cm—1, suggests, by negative arguments, that either the unobserved fundamental ν9 does not lie in the frequency range 900 to 1100 cm—1, or ∣ ζ9,14Z ∣<0.2. The assignment of the unobserved fundamental vibrations of diborane is discussed in the light of this evidence.

Perturbations in the infrared spectrum of monofluoroacetylene, and their relationship to intramolecular vibrational energy redistribution

Several high-order vibration-rotation perturbations in the high-resolution infrared spectrum of monofluoroacetylene, HCCF, are assigned and analyzed in detail. They result in avoided crossings in the rotational structure of several bands, and precise values for the effective high-order terms in the Hamiltonian have been determined. The significance of these results for intramolecular vibrational redistribution is discussed.

Can 4D-Var use dynamical information from targeted observations of a baroclinic structure?

Targeted observations are generally taken in regions of high baroclinicity, but often show little impact. One plausible explanation is that important dynamical information, such as upshear tilt, is not extracted from the targeted observations by the data assimilation scheme and used to correct initial condition error. This is investigated by generating pseudo targeted observations which contain a singular vector (SV) structure that is not present in the background field or routine observations, i.e. assuming that the background has an initial condition error with tilted growing structure. Experiments were performed for a single case-study with varying numbers of pseudo targeted observations. These were assimilated by the Met Office four-dimensional variational (4D-Var) data assimilation scheme, which uses a 6 h window for observations and background-error covariances calculated using the National Meteorological Centre (NMC) method. The forecasts were run using the operational Met Office Unified Model on a 24 km grid. The results presented clearly demonstrate that a 6 h window 4D-Var system is capable of extracting baroclinic information from a limited set of observations and using it to correct initial condition error. To capture the SV structure well (projection of 0.72 in total energy), 50 sondes over an area of 1×106 km2 were required. When the SV was represented by only eight sondes along an example targeting flight track covering a smaller area, the projection onto the SV structure was lower; the resulting forecast perturbations showed an SV structure with increased tilt and reduced initial energy. The total energy contained in the perturbations decreased as the SV structure was less well described by the set of observations (i.e. as fewer pseudo observations were assimilated). The assimilated perturbation had lower energy than the SV unless the pseudo observations were assimilated with the dropsonde observation errors halved from operational values. Copyright © 2010 Royal Meteorological Society

l-resonance perturbations in infrared perpendicular bands

The effects of ℓ-type resonance on rovibrational bands in infrared spectra are reviewed. Observed spectra are compared with computer-simulated spectra obtained by solving the Hamiltonian matrix numerically and calculating the true (perturbed) wavenumber and intensity of each line in the band. The most obvious effects in the spectra are shown to result from intensity perturbations rather than line-shifts; in oblate symmetric tops the Q branch structure near the band center may show anomalies due to ℓ-resonance even at quite low resolution and even when the accidental resonance is not very exact. Numerical values of ℓ-doubling constants are obtained for several cyclopropane bands by comparing observed band contours at about 0.2-cm−1 resolution with computed contours. Although the constants are not determined with great precision, the sign of the ℓ-doubling constants is determined unambiguously.

Intensity perturbations due to the v3, v4 Coriolis interaction in methane

Formulas are obtained for the intensity asymmetry (Herman-Wallis) factors in the ν3 and ν4 fundamentals of methane due to the ζ34 Coriolis interaction. The results are also applicable to the ν3 and ν4 bands of SF6.

Singular vector ensemble forecasting systems and the prediction of flow dependent uncertainty

The ECMWF ensemble weather forecasts are generated by perturbing the initial conditions of the forecast using a subset of the singular vectors of the linearised propagator. Previous results show that when creating probabilistic forecasts from this ensemble better forecasts are obtained if the mean of the spread and the variability of the spread are calibrated separately. We show results from a simple linear model that suggest that this may be a generic property for all singular vector based ensemble forecasting systems based on only a subset of the full set of singular vectors.