Skip to main content
Have a personal or library account? Click to login
A simpler formulation of forecast sensitivity to observations: application to ensemble Kalman filters Cover

A simpler formulation of forecast sensitivity to observations: application to ensemble Kalman filters

Open Access
|Dec 2012

Figures & Tables

Fig. 1. 

Schematic of the perceived forecast error verified against the analysis at the verification time t from two forecasts started from the analysis at t=0 h, and from the analysis at t=−6 h . Since the forecast started at t=−6 h serves as first guess for the analysis at t=0 h, the only difference between the two forecasts is the assimilation of the observations y 0 at t=0 h. Adapted from Langland and Baker (2004).

Table 1. Forecast error reduction of each formulation averaged over time and over 10 experiments. Truth and analysis show the average forecast error reduction verified against the truth, and against the analysis, respectively. Values in the parenthesis show the standard deviation of the 10 experiments. All values are multiplied by 100. No localization was used in the LETKF

6 hr 12 hr 1 d 2 d 3 d 5 d 7 d Truth −24 (0) −33 (0) −57 (0) −122 (3) −218 (8) −606 (29) −1171 (52) Analysis −27 (0) −32 (0) −53 (0) −116 (3) −211 (8) −593 (29) −1152 (51) (4-ADJ) −23 (0) −28 (0) −46 (0) −101 (2) −181 (6) −467 (22) −840 (82) (5-OLD) −27 (0) −33 (0) −54 (0) −121 (2) −218 (7) −421 (21) −260 (27) (6-NEW) −31 (0) −37 (0) −60 (1) −128 (3) −220 (9) −564 (41) −1167 (86)
Fig. 2. 

Skill score of the estimation of the total forecast error reduction verified against the true forecast error reduction. Black, red and blue lines show the results of equations (4-ADJ), (5-OLD) and (6-NEW). The vertical error bars show the standard deviation of the 10 experiments. No localization was used in the LETKF.

Fig. 3. 

Same as Fig. 2 but with only 10 members and using localization in the analysis. No localization was applied for the ensemble-based observation impact computations.

Fig. 4. 

As per Fig. 3 but now with localization used in the observation impact computations. Red and blue lines show the results of equations (5-OLD) and (6-NEW) with a fixed localization (Fxloc) function as used in the LETKF analysis. Dotted green and light blue lines show the results of equation (6-NEW) with localization functions moving with non-linear evolution (NL-loc) and constant speed equal to the group velocity (CL-loc).

Table 2. Same as Table 1 but with 10 ensemble members and localization when indicated

6 h 12 h 1 d 2 d 3 d 5 d 7 d Truth −23 (0) −32 (0) −56 (0) −121 (2) −213 (6) −582 (27) −1144 (47) Analysis −25 (0) −31 (0) −51 (0) −115 (2) −206 (6) −570 (27) −1126 (46) (4-ADJ) −21 (0) −25 (0) −39 (0) −77 (1) −132 (4) −344 (17) −627 (52) (5-OLD) Noloc −12 (0) −14 (0) −24 (1) −50 (3) −69 (4) 136 (22) 1163 (65) (6-NEW) Noloc −32 (0) −39 (0) −62 (1) −133 (4) −227 (6) −574 (41) −1220 (80) (5-OLD) Fxloc −25 (0) −30 (0) −47 (0) −89 (2) −117 (4) −9 (9) 486 (29) (6-NEW) Fxloc −28 (0) −33 (0) −50 (0) −86 (2) −106 (4) −125 (11) −162 (29) (6-NEW) NL-loc −24 (0) −24 (0) −47 (0) −107 (2) −178 (4) −422 (27) −788 (57) (6-NEW) CL-loc −28 (0) −33 (0) −54 (0) −99 (2) −151 (5) −284 (18) −462 (30)
Fig. 5. 

Example of the evolution of the localization function with two methods. (a) Localization function with the non-linear incremental evolution (NL-loc); and (b) localization function moving with constant group velocity (CL-loc).

Fig. 6. 

Time series of the total forecast error reduction at 5-d forecast from cycle 1000 to 1100. Pink lines show the true forecast error reduction, the black line is the adjoint estimation, the blue line is the ensemble estimation with fixed localization and the dotted green line is the ensemble estimation with non-linear evolution of the localization. A total of 40 observations are assimilated on each analysis time.

Fig. 7. 

Fraction of the observations that are estimated to decrease the forecast error in each method. Purple solid line shows the result of data denial experiments, denoted as an Observing System Experiment (OSE).

Fig. 8. 

Time average of observation impact estimates at each point on 5-d forecast with (a) 1000 and (b) 100 samples, using the adjoint approach eq. (4), the LLK10 formulation eq. (5) and the new ensemble formulation eq. (6), both with fixed localization. The observation at grid point 11 has an error variance of 0.8 rather than the assumed value of 0.2.

Fig. A1. 

Similarity indices between the non-linear and tangent linear evolution of the perturbations. Standard deviation of the perturbations are 0.2 (red), 0.5 (yellow), 1.0 (green) and 2.0 (blue). Values are averaged over 1000 realizations.

Table A1. Similarity indices between non-linear and tangent linear evolution of the perturbations (1 day evolution) with various perturbation amplitudes

Amplitude SI 10−5 0.9999999999953 10−4 0.9999999994920 10−3 0.9999999423196 10−2 0.9999937996438 10−1 0.9994241704626
Language: English
Page range: 18462 - 18462
Submitted on: Apr 1, 2012
Accepted on: Aug 25, 2012
Published on: Dec 1, 2012
Published by: Stockholm University Press
In partnership with: Paradigm Publishing Services

© 2012 Eugenia Kalnay, Yoichiro Ota, Takemasa Miyoshi, Junjie Liu, published by Stockholm University Press
This work is licensed under the Creative Commons Attribution 4.0 License.