Skip to main content
Have a personal or library account? Click to login
Data compression in the presence of observational error correlations Cover

Data compression in the presence of observational error correlations

By:   
Open Access
|Jan 2019

Figures & Tables

Fig. 1.

Top: Eigenvalues of MMT (20) (yellow line), R (blue triangles) and B (red circles) for the case of SOAR correlation functions with LB=5 and LR=0.1. Also plotted are the eigenvalues of Pa (purple crosses) when compressed observations retaining 75% of the total ER are assimilated (13 compressed observations in this case). Bottom: ERc (left), DFSc (middle) and trace(BPa)c (right) as a function of the number of compressed observations ordered according to the eigenvalues of MMT. The dashed lines indicate the number of compressed observations needed to achieve 75% of the total value. These numbers are also given in Table 1.

Fig. 2.

As in Fig. 1 but for LR=5. The black line of the last panel shows tr(BPa)c as a function of the number of compressed observations when the assimilation is performed using large-scale observations first, and the red line when the assimilation is performed using small-scale observations first.

Fig. 3.

As in Fig. 1 but for LR=10.

Table 1.

75% of the value of ER, DFS, and tr(BP) when all observations are assimilated for different values of LR. In brackets are the number of compressed observations need to achieve these values.

Criterion0.75 × ERall0.75 × DFSall0.75 × tr(BPa)allLR=0.17.08 (13 obs)9.48 (15 obs)14.5 (8 obs)LR=58.32 (23 obs)12.0 (24 obs)12.0 (10 large-scale obs, 29 small-scale obs)LR=1016.9 (18 obs)17.2 (20 obs)13.4 (27 obs)
Fig. 4.

Analysis error correlation structure for the case illustrated in Fig. 2. The number of compressed observations assimilated are chosen to conserve 75% of ER. Blue: compressed observations favour large scales and red: compressed observations favour small scales.

Fig. 5.

Rows of the observation operator matrix for the five strategies for reducing the observation data detailed in Section 4. The optimal strategies are illustrated for the first observation time.

Table 2.

The wave-numbers selected by the optimal Fourier data compression method at each observation time for the two different types of observations. The first observation time corresponds to the observation operator plotted in Fig. 5.

Ob timeLR=0.1LR = 2110, 11, 9, 12, 812, 13, 11, 10, 15210, 11, 12, 13, 914, 12, 13, 20, 15312, 9, 11, 10, 1315, 14, 17, 20, 18412, 10, 11, 9, 817, 19, 18, 16, 20513, 14, 11, 12, 1019, 17, 15, 18, 16
Fig. 6.

Top left: pseudo inverse of Xf(Xf)T illustrated for the first observation time. Other panels: HcTRc1Hc (the accuracy of the compressed observation in state space) for the 5 different thinning strategies when observation errors are uncorrelated (middle panels) and correlated (right panels) illustrated for the first observation time.

Fig. 7.

(a) Ensemble spread, (b) entropy computed using the ensemble estimate of (11), (c) log10ER, and (d) log10 of the condition number of the analysis error covariance matrix for the five different thinning strategies detailed in Section 3.1. Solid lines represent results when the observation errors are uncorrelated, dashed lines represent results when the observation errors are correlated. Results are averaged over 200 experiments with different realisations of the observation and model error.

Language: English
Page range: 1634937 - 1634937
Published on: Jan 1, 2019
Published by: Stockholm University Press
In partnership with: Paradigm Publishing Services

© 2019 A.M. Fowler, published by Stockholm University Press
This work is licensed under the Creative Commons Attribution 4.0 License.