Cyclogenetic Perturbations and Analysis Errors Decomposed into Singular Vectors

2005 ◽  
Vol 62 (7) ◽  
pp. 2234-2247 ◽  
Author(s):  
Chris Snyder ◽  
Gregory J. Hakim

Abstract Singular vectors (SVs) have been applied to cyclogenesis, to initializing ensemble forecasts, and in predictability studies. Ideally, the calculation of the SVs would employ the analysis error covariance norm at the initial time or, in the case of cyclogenesis, a norm based on the statistics of initial perturbations, but the energy norm is often used as a more practical substitute. To illustrate the roles of the choice of norm and the vertical structure of initial perturbations, an upper-level wave with no potential vorticity perturbation in the troposphere is considered as a typical cyclogenetic perturbation or analysis error, and this perturbation is then decomposed by its projection onto each energy SV. All calculations are made, for simplicity, in the context of the quasigeostrophic Eady model (i.e., for a background flow with constant vertical shear and horizontal temperature gradient). Viewed in terms of the energy SVs, the smooth vertical structure of the typical perturbation, as well as its evolution, results from strong cancellation between the growing and decaying SVs, most of which are highly structured and tilted in the vertical. A simpler picture, involving less cancellation, follows from decomposition of the typical perturbation into SVs using an alternative initial norm, which is based on the relation between initial norms and the statistics of initial perturbations together with the empirical assumption that the initial perturbations are not dominated by interior potential vorticity. Differences between the energy SVs and those based on the alternative initial norm can be understood by noting that the energy norm implicitly assumes initial perturbations with second-order statistics given by the covariance matrix whose inverse defines the energy norm. Unlike the “typical” perturbation, perturbations with those statistics have large variance of potential vorticity in the troposphere and fine vertical structure. Finally, a brief assessment is presented of the extent to which the upper wave, and more generally the alternative initial norm, is representative of cyclogenetic perturbations and analysis errors. There is substantial evidence supporting deep perturbations with little vertical structure as frequent precursors to cyclogenesis, but surrogates for analysis errors are less conclusive: operational midlatitude analysis differences have vertical structure similar to that of the perturbations implied by the energy norm, while short-range forecast errors and analysis errors from assimilation experiments with simulated observations are more consistent with the alternative norm.

2003 ◽  
Vol 131 (8) ◽  
pp. 1741-1758 ◽  
Author(s):  
Thomas M. Hamill ◽  
Chris Snyder ◽  
Jeffrey S. Whitaker

Abstract Approximations to flow-dependent analysis-error covariance singular vectors (AEC SVs) were calculated in a dry, T31 L15 primitive-equation global model. Sets of 400-member ensembles of analyses were generated by an ensemble-based data assimilation system. A sparse network of simulated rawinsonde observations were assimilated, and a perfect model was assumed. Ensembles of 48-h forecasts were also generated from these analyses. The structure of evolved singular vectors was determined by finding the linear combination of the forecast ensemble members that resulted in the largest forecast-error variance, here measured in a total-energy norm north of 20°N latitude. The same linear combination of analyses specifies the initial-time structure that should evolve to the forecast singular vector under assumptions of linearity of error growth. The structures of these AEC SVs are important because they represent the analysis-error structures associated with the largest forecast errors. If singular vectors using other initial norms have very different structures, this indicates that these structures may be statistically unlikely to occur. The European Centre for Medium-Range Weather Forecasts currently uses singular vectors using an initial total-energy norm [“total-energy singular vectors” or (TE SVs)] to generate perturbations to initialize their ensemble forecasts. Approximate TE SVs were also calculated by drawing an initial random ensemble with perturbations that were white in total energy and applying the same approach as for AEC SVs. Comparing AEC SVs and approximate TE SVs, the AEC SVs had maximum amplitude in midlatitudes near the tropopause, both at the initial and evolved times. The AEC SVs were synoptic in scale, deep, and did not appear to be geographically localized nor tilted dramatically upshear. This contrasts with TE SVs, which started off relatively smaller in scale, were tilted upshear, and had amplitudes typically largest in the lower to midtroposphere. The difference between AEC SVs and TE SVs suggests that operational ensemble forecasts based on TE SVs could be improved by changing the type of singular vector used to generate initial perturbations. This is particularly true for short-range ensemble forecasts, where the structure of the forecast ensemble is more closely tied to the analysis ensemble.


2007 ◽  
Vol 135 (7) ◽  
pp. 2754-2777 ◽  
Author(s):  
Jean-François Caron ◽  
M. K. Yau ◽  
Stéphane Laroche

Abstract This paper presents a diagnostic study of the evolution of initial corrections obtained from the key analysis error algorithm that minimizes the short-range (24 h) forecast errors for four specific events poorly forecasted over the eastern part of North America. A potential vorticity (PV) perspective is employed. It is shown that the modification to the low-level structure at the initial time is mainly attributed to the modification of the low-level PV distribution, while changes in the upper-level structure are attributed to the modification of the upper-level PV distribution. The low-level corrections grow mainly through background surface potential temperature advection by the wind corrections attributable to the interior PV corrections. Changes in the diabatic processes and the vertical alignment of low-level PV corrections by differential PV advection also increase the magnitude of the low-level corrections with time. The upper-level corrections grow by advection of background PV from wind corrections. However, the cause of these latter wind corrections responsible for upper-level background PV advection varies from case to case. An investigation of the relative importance of the low-level and of the upper-level initial corrections to produce the final-time corrections also reveals strong variability between cases. Finally, comparison of two cases in which the key analysis errors propagate vertically with two others without significant vertical propagation shows how the relative position of the key analysis errors with respect to the structure of the background flow can influence the evolution of the initial corrections.


2015 ◽  
Vol 143 (10) ◽  
pp. 4064-4087 ◽  
Author(s):  
Ryan D. Torn ◽  
Glen S. Romine

Abstract The role of upstream subsynoptic forecast errors on forecasts of two different central Oklahoma severe convection events (19 and 31 May 2013) characterized by strong synoptic forcing during the Mesoscale Predictability Experiment (MPEX) are evaluated by applying the ensemble-based sensitivity technique to WRF ensemble forecasts with explicit convection. During both cases, the forecast of the timing and intensity of convection over central Oklahoma is modulated by the southward extent of upstream midtropospheric potential vorticity anomalies that are moving through the base of a larger-scale upstream trough but pass by central Oklahoma prior to convective initiation. In addition, the convection forecasts are also sensitive to the position of lower-tropospheric boundaries, such that moving the boundaries in a manner that would lead to increased equivalent potential temperature over central Oklahoma prior to convective initiation leads to more precipitation. Statistical PV inversion and correlation calculations suggest that the midtropospheric PV and near-surface boundary sensitivities are not independent; the winds associated with the PV error can modulate the position of the lower-tropospheric boundary through advection in a manner consistent with the implied sensitivity. As a consequence, it appears that reducing the uncertainty in specific upstream subsynoptic features prior to convective initiation could improve subsequent forecasts of severe convection.


2008 ◽  
Vol 136 (5) ◽  
pp. 1565-1581 ◽  
Author(s):  
David A. Gold ◽  
John W. Nielsen-Gammon

Abstract Observational and modeling studies have shown that shear and instability are powerful predictors of the likelihood of severe weather and tornadoes. To the extent that upper-tropospheric forecast errors can be described as potential vorticity (PV) anomalies on the forecasted PV field, knowing (and being able to quantify) the effects of such errors on shear and instability would allow forecasters to anticipate the effects of those errors on the likely mode of severe weather. To test the sensitivity of the severe convective environment to PV fluctuations, a PV inversion framework is adopted that utilizes nonlinear balance. The observed PV field is modified in a way that mimics realistic perturbations of trough intensity, location, or shape. Soundings, including moisture profiles, are reconstructed from the balanced geopotential height field assuming that air parcels conserve mixing ratio while their isentropic surfaces are displaced upward or downward by the addition of anomalous PV. Unperturbed balanced soundings agree reasonably well with full, unbalanced soundings, and differences are attributable to departures from nonlinear balance in areas of strong vorticity or acceleration. Balanced vertical wind profiles do not include the effects of friction, so the vertical shear of the balanced wind departs unacceptably from total shear within the lowest 1 km of the troposphere. The balanced wind perturbations are added to the total analyzed shear profile to estimate the effect of PV perturbations on shear and storm-relative helicity. By this process, the importance of typical or hypothesized upper-tropospheric forecast errors may be addressed in an idealized, case-study, or operational context.


2007 ◽  
Vol 135 (2) ◽  
pp. 249-266 ◽  
Author(s):  
Jean-François Caron ◽  
M. K. Yau ◽  
Stéphane Laroche ◽  
Peter Zwack

Abstract The characteristics of the initial corrections obtained from the Canadian Meteorological Centre (CMC) energy-norm-based key analysis error algorithm that minimizes short-range (24 h) forecast errors were investigated for four specific CMC operational analyses. The results show that both the rotational and the divergent components of the initial corrections are strongly out of balance. Some dispersive modes are also present in the mass component of the initial corrections. The results from one experiment where the initial state errors were known suggest that the current algorithm always selects a set of unbalanced initial corrections with more mass correction than wind correction, regardless of the characteristics of the real initial condition errors. Comparison with observational data showed that the corrected analysis is systematically farther away from the observations than the control analysis even in large forecast error events where most of the forecast errors are believed to have originated from errors in the initial state.


2007 ◽  
Vol 135 (2) ◽  
pp. 267-280 ◽  
Author(s):  
Jean-François Caron ◽  
M. K. Yau ◽  
Stéphane Laroche ◽  
Peter Zwack

Abstract This study examines a few approaches to isolate the balanced component of the initial corrections from the Canadian Meteorological Centre energy-norm-based key analysis error algorithm, in an attempt to capture the part of the key analysis errors responsible for short-range forecast errors. The best results were obtained with the nonlinear balance potential vorticity (PV) inversion technique. It was shown that the PV component of the initial corrections contains the essential information for reducing short-range forecast errors. The remaining imbalance part of the initial corrections does not grow in time and does not contribute to the improvement of the forecast. The removal of the imbalance part of the initial corrections makes the corrected analysis slightly closer to the observations, but remains systematically farther away as compared with the original analysis. Thus the balanced part of the key analysis errors cannot justifiably be associated to analysis errors. A methodology to balance the divergent part of the initial corrections, which reduces significantly the spinup in the vertical motion corrections, is also presented. Finally, in light of the results presented in this paper, some recommendations to improve the key analysis error algorithm are proposed.


2005 ◽  
Vol 133 (3) ◽  
pp. 567-578 ◽  
Author(s):  
Gregory J. Hakim

Abstract The dominant vertical structures for analysis and forecast errors are estimated in midlatitudes using a small ensemble of operational analyses. Errors for fixed locations in the central North Pacific and eastern North America are selected for comparing errors in regions with relatively low and high observation density, respectively. Results for these fixed locations are compared with results for zonal wavenumber 9, which provides a representative sample of baroclinic waves. This study focuses on deviations from the ensemble mean for meridional wind and temperature at 40°N; these quantities are chosen for simplicity and because they capture dynamical and thermodynamical aspects of midlatitude baroclinic waves. Results for the meridional wind show that analysis and forecast errors share the same dominant vertical structure as the analyses. This structure peaks near the tropopause and decays smoothly toward small values in the middle and lower troposphere. The dominant vertical structure for analysis errors exhibits upshear tilt and peaks just below the tropopause, suggesting an asymmetry in errors of the tropopause location, with a bias toward greater errors for downward tropopause displacements. The dominant vertical structure for temperature analysis errors is distinctly different from temperature analyses. Analysis errors have a sharp peak in the lower troposphere, with a secondary structure near the tropopause, whereas forecast errors and analyses show a dipole straddling the tropopause and smooth vertical structure, consistent with potential vorticity anomalies due to variance in tropopause position. Linear regression of forecast errors onto analysis errors for the western North Pacific is used to assess the nonseparable zonal-height structure of errors and their propagation. Analysis errors near the tropopause rapidly develop into a spreading wave packet, with a group speed that matches the mean zonal wind speed of 31 m s−1. A complementary calculation for the regression of 24-h forecast errors onto analysis errors shows that forecast errors originate from analysis errors in the middle and upper troposphere. These errors rapidly expand in the vertical to span the troposphere, with a peak at the tropopause.


2018 ◽  
Vol 35 (9) ◽  
pp. 1737-1752 ◽  
Author(s):  
Dae-Hui Kim ◽  
Hyun Mee Kim

AbstractIn this study, the effect of assimilating Himawari-8 (HIMA-8) atmospheric motion vectors (AMVs) on forecast errors in East Asia is evaluated using observation system experiments based on the Weather Research and Forecasting Model and three-dimensional variational data assimilation system. The experimental period is from 1 August to 30 September 2015, during which both HIMA-8 and Multifunctional Transport Satellite-2 (MTSAT-2) AMVs exist. The energy-norm forecast error based on the analysis of each experiment as reference was reduced more by replacing MTSAT-2 AMVs with HIMA-8 AMVs than by adding HIMA-8 AMVs to the MTSAT-2 AMVs. When the HIMA-8 AMVs replaced or were added to MTSAT-2 AMVs, the observation impact was reduced, which implies the analysis–forecast system was improved by assimilating HIMA-8 AMVs. The root-mean-square error (RMSE) of the 500-hPa geopotential height forecasts based on the analysis of each experiment decreases more effectively when the region lacking in upper-air wind observations is reduced by assimilating both MTSAT-2 and HIMA-8 AMVs. When the upper-air radiosonde (SOUND) observations are used as reference, assimilating more HIMA-8 AMVs decreases the forecast error. Based on various measures, the assimilation of HIMA-8 AMVs has a positive effect on the reduction of forecast errors. The effects on the energy-norm forecast error and the RMSE based on SOUND observations are greater when HIMA-8 AMVs replaced MTSAT-2 AMVs. However, the effects on the RMSE of the 500-hPa geopotential height forecasts are greater when both HIMA-8 and MTSAT-2 AMVs were assimilated, which implies potential benefits of assimilating AMVs from several satellites for forecasts over East Asia depending on the choice of measurement.


2018 ◽  
Vol 57 (4) ◽  
pp. 1011-1019 ◽  
Author(s):  
H. F. Dacre ◽  
N. J. Harvey

ABSTRACTVolcanic ash poses an ongoing risk to safety in the airspace worldwide. The accuracy with which volcanic ash dispersion can be forecast depends on the conditions of the atmosphere into which it is emitted. In this study, meteorological ensemble forecasts are used to drive a volcanic ash transport and dispersion model for the 2010 Eyjafjallajökull eruption in Iceland. From analysis of these simulations, the authors determine why the skill of deterministic-meteorological forecasts decreases with increasing ash residence time and identify the atmospheric conditions in which this drop in skill occurs most rapidly. Large forecast errors are more likely when ash particles encounter regions of large horizontal flow separation in the atmosphere. Nearby ash particle trajectories can rapidly diverge, leading to a reduction in the forecast accuracy of deterministic forecasts that do not represent variability in wind fields at the synoptic scale. The flow‐separation diagnostic identifies where and why large ensemble spread may occur. This diagnostic can be used to alert forecasters to situations in which the ensemble mean is not representative of the individual ensemble‐member volcanic ash distributions. Knowledge of potential ensemble outliers can be used to assess confidence in the forecast and to avoid potentially dangerous situations in which forecasts fail to predict harmful levels of volcanic ash.


2013 ◽  
Vol 7 (1) ◽  
pp. 37-50
Author(s):  
Masanori Yamasaki

This paper describes results from numerical experiments which have been made toward a better understanding of tropical cyclone formation. This study uses a nonhydrostatic version of the author’s mesoscale-convection-resolving model that was developed in the 1980s to improve paramerization schemes of moist convection. In this study the horizontal grid size is taken to be 20 km in an area of 6,000 km x 3,000 km, and a non-uniform coarse grid is used in two areas to its north and south. Results from two numerical experiments are presented; one (case 1) without any environmental flow, and the other (case 2) with an easterly flow without low-level vertical shear. Three circular buoyancy perturbations are placed in the west-east direction at the initial time. Convection is initiated in the imposed latently unstable (positive CAPE) area. In both cases, a vortex with a pressure low is formed, and two band-shaped convective systems are formed to the north and the south of the vortex center. The vortex and two convective systems are oriented in the westsouthwest – eastnortheast direction, and their horizontal scales are nearly 2,000 km. In case 1, the band-shaped convective system on the southern side is stronger, and winds are stronger just to its south. In contrast, in case 2, the northern convective system is stronger, and winds are stronger just to its north. Therefore, the distributions of the equivalent potential temperature in the boundary layer and latent instability (positive buoyancy of the rising air) are also quite different between cases 1 and 2. The TC formation processes in these different cases are discussed, with an emphasis on the importance of examining the time change of latent instability field.


Sign in / Sign up

Export Citation Format

Share Document