Changes between Version 5 and Version 6 of EnsembleGeneration


Ignore:
Timestamp:
Jan 29, 2021, 10:00:55 AM (4 months ago)
Author:
lnerger
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • EnsembleGeneration

    v5 v6  
    77To start ensemble-based data assimilation one has to generate an ensemble of model states that represents both the state estimate (provided by the ensemble mean state) and the uncertainty of the state estimate, given by the ensemble spread. In case of ensemble Kalman filters, the uncertainty is estimate by the sample covariance matrix of the ensemble, i.e. the covariance matrix of the deviations of the ensemble states from the ensemble mean state.
    88
    9 There are different possibilities to generate an ensemble. With PDAF we currently provide support for one sampling method, the co-called second-order exact sampling from a model trajectory. This method was introduced by Pham et al. (see e.g. Pham, Monthly Weather Review 129 (2001) 1194-1207). The method is as follows:
     9There are different possibilities to generate an ensemble. With PDAF we currently provide support for one sampling method, the so-called second-order exact sampling from a model trajectory. This method was introduced by Pham et al. (see e.g. Pham, Monthly Weather Review 129 (2001) 1194-1207). The method is as follows:
    10101. Run the model with which the data assimilation will be performed and write a trajectory of model snapshots into files.
    11111. In a separate program read the model snapshots. Compute the temporal mean state and subtract it from all snapshots. Now compute the singular value decomposition of the resulting matrix. This is known as the EOF-decomposition. It provides the singular vectors (EOFS, empirical orthogonal functions) and singular values (EOF coefficients). The number of EOFs is usually '''r''' if the trajectory contains '''r+1''' model snapshots. Store the singular vectors and singular values in a file.