Skip to main content
Have a personal or library account? Click to login
Impacts on sea ice analyses from the assumption of uncorrelated ice thickness observation errors: Experiments using a 1D toy model Cover

Impacts on sea ice analyses from the assumption of uncorrelated ice thickness observation errors: Experiments using a 1D toy model

Open Access
|Jan 2018

Figures & Tables

Figure 1.

Spectral variances of the Kalman gain matrix for a fully observed system. Line with markers is the case where the length scale of observation error correlations is greater than that of the background error correlations; solid line is for the case where the length scale of the background error correlations is greater than that of the observation error correlations.

Figure 2.

Temporal evolution of spread in the five-hundred-member ensemble during the thirty-day freeze-up period. The thickness spread (a) increased gradually to a value of approximately 0.23 m. The concentration (b) and velocity (c) spreads increased abruptly at the beginning of the simulation, but gradually decreased as the mean concentration approached 100% and the ice strengthened.

Figure 3.

Ensemble of sea ice model states with the true state following the thirty-day freeze-up period. The mean and the 10th/90th percentiles of the ensemble are plotted to illustrate the spread. Note that where the true ice concentration was near 100%, the spread of the ice concentration and the velocity errors were the least but there was little correlation between the state and the thickness spread.

Figure 4.

Prescribed true observation error standard deviation as a function of ice thickness. As the true ice thickness increases beyond 0.3m, the observation error standard deviation increases quadratically, meaning that there is low observation accuracy for regions of thick i.e.

Figure 5.

Evaluation of the ice thickness observation error standard deviation function at the true state mapped to the observation locations. The background error standard deviations (solid black) can be compared to the observation error standard deviations (grey), revealing that the background error standard deviations are relatively constant across the domain while the observation error standard deviations are more variable. We can expect greater impact from the assimilation in two isolated regions where the observation error standard deviations are less than the background error standard deviations.

Table 1.

Correlation structure of each observation error covariance matrix.

MatrixCorrelation StructureDecorrelation Length (km)R150Exponential150R50Exponential50RdiagDiagonal0
Table 2.

Four variants of the analysis error covariance matrix that were computed for each twin experiment.

ACalculationAidealAideal=(I-KH)B(I-KH)T+KRtrueKT      K=BHT(HBHT+Rtrue)-1AestAest=(I-KH)B(I-KH)T+KRestKT     K=BHT(HBHT+Rest)-1AactualAactual=(I-KH)B(I-KH)T+KRtrueKT         K=BHT(HBHT+Rest)-1BNo assimilation
Figure 6.

Correlation structure of ice thickness errors extracted from the background error covariance matrix compared to the prescribed error correlation structures for the observation error covariance matrices.

Table 3.

Analysis error standard deviations for ice thickness (h), ice concentration (a) and ice velocity (u) for the case where Rtrue=R150 (left column), Rtrue=R50 (centre column), and Rtrue=Rdiag (right column).

Rtrue=R150Rtrue=R50Rtrue=RdiagResthauhauhau(10-1m)a (10-2)(10-2m/s)(10-1m)a (10-2)(10-2m/s)(10-1m)a (10-2)(10-2m/s)σb2.237.321.382.237.321.382.237.321.38Rtrueσesta1.996.831.272.026.951.291.976.941.28σactuala1.996.831.272.026.951.291.976.941.28Rdiagσesta1.976.941.281.976.941.281.976.941.28σactuala2.236.981.322.076.991.301.976.941.28Rinflσesta2.097.111.322.097.111.322.097.111.32σactuala2.097.051.312.057.041.312.037.021.30
Figure 7.

Spectral densities of thickness analysis error variances for the experiment where the true observation error covariance matrix had an exponential correlation structure with decorrelation length of 150 km. Panel (a) Rest=R150, (b) Rest=R50, (c) Rest=Rinfl, (d) Rest=Rdiag.

Figure 8.

Spectral densities of (a) ice concentration and (b) ice velocity, analysis error variances for the experiment where Rtrue=R150 and Rest=Rdiag.

Figure 9.

Estimated and actual thickness analysis error standard deviation as a function of the inflation factor for (a) Rtrue=Rdiag(b) Rtrue=R50 and (c) Rtrue=R150. Each experiment had an optimal inflation factor between 1 and 4. Underestimation of the inflation factor resulted in a divergence between the estimated and actual analysis error standard deviation. Overestimation of the inflation factor had a lesser effect.

Table A1.

Default model parameter values.

ParameterUnitsDefault valueCa–0.0015Cw–0.0015C–20e–2h0m0.05PNm-15000ρikgm-3910ρwkgm-31035ρakgm-31.3Δxm1000Δth0.5
Language: English
Page range: 1445379 - 1445379
Submitted on: Sep 25, 2017
Accepted on: Feb 20, 2018
Published on: Jan 1, 2018
Published by: Stockholm University Press
In partnership with: Paradigm Publishing Services

© 2018 Graham Stonebridge, K. Andrea Scott, Mark Buehner, published by Stockholm University Press
This work is licensed under the Creative Commons Attribution 4.0 License.