350 likes | 503 Views
Observations and Ocean State Estimation: Impact, Sensitivity and Predictability. Andy Moore University of California Santa Cruz Hernan Arango Rutgers University. Outline. State estimation Observation impact Information content Observation sensitivity
E N D
Observations and Ocean State Estimation: Impact, Sensitivity and Predictability Andy Moore University of California Santa Cruz Hernan Arango Rutgers University
Outline • State estimation • Observation impact • Information content • Observation sensitivity • Forecast error and predictability • Examples from the California Current
may include ocean dynamics as in 4D-Var Ocean State Estimation and Data Assimilation Ocean state vector: From Bayes’ Theorem: Posterior Prior Gain Innovation Obs Obs operator Kalman gain: Obs error cov Prior error cov TL obs operator
identified iteratively using 4D-Var Practicalities of 4D-Var Typically done using conjugate gradients (CG) Practical Gain Matrix Alternatively, view CG and 4D-Var as a function: 4D-Var
Two Views of 4D-Var OR Posterior Prior Innovation
The Role of Observations Q: What is the influence of the observations on the analysis? Observation impact: Observation sensitivity: Predictability:
Observation Impact Q: How does each obs contribute to the analysis? Consider a function of the circulation:
Observation Sensitivity Q: How does the analysis change if the observations change? Consider again the function Let
Observations and Predictability Q: How does each obs contribute to the forecast predictability? Consider now an ensemble of forecast function values for obtained by perturbing priors and obs. Expected forecast error variance: Forecast error covariance
Observations and Predictability Forecast error covariance Tangent linear model Control priors where: Posterior error covariance Tangent Linear 4D-Var Adjoint Linear 4D-Var
37N 500m Mesoscale eddies The California Current
fb(t), Bf bb(t), Bb xb(0), Bx Previous assimilation cycle The California Current The Regional Ocean Modeling System (ROMS) COAMPS forcing ECCO open boundary conditions 30km, 10 km & 3 km grids, 30- 42 levels Veneziani et al (2009) Broquet et al (2009) Moore et al (2010)
Observations (y) CalCOFI & GLOBEC ~90% SST & SSH EN3 Ingleby and Huddleston (2007) TOPP Elephant Seals (APB) ~10% ARGO Data from Dan Costa
Sequential Data Assimilation: July 2002-Dec2004 Observations Observations Observations prior prior prior Data Assimilation Data Assimilation Data Assimilation Posterior Posterior Posterior 7 days Forecast Forecast Forecast
XBT CTD All in situ data: July 2002 – Dec 2004 ARGO TOPP
Observation Impact Q: How does each obs contribute to the analysis?
Observation Impacts on Transport 7day average transport 10km ROMS Transport increment = (Posterior-Prior)
Example: 37N Transport CUC CC No assim JAS time mean alongshore Flow (10km, 42 lev) CC = California Current CUC = California Under Current CC CUC With Data Assimilation
Poleward Equatorward Offshore Onshore Prior alongshore transport (CC+CUC+CJ) Prior cross-shore transport
rms Analysis Cycle – Observation Impacts 10km ROMS Poleward Alongshore transport Equatorward Cross-shore transport Offshore Onshore (Moore et al, 2011c)
IGW Adjoint CTW (GT) IGW CTW (G) IGW Gyre Circulation Alongshore Transport Impacts Sv (10-5) SSH
Degrees of Freedom (dof) y y Observation vector Observation vector x x “Perfect World” dof ~ Nobs Redundancy dof < Nobs
Information Content of Obs Degrees of freedom of obs (dof): Cardinali et al, 2004; Desroziers et al., 2009 Bennett & McIntosh, 1982 Tr{KH}/Nobs vs assimilation cycle upper & lower bounds 30km Only ~10% of obs contain independent info (Moore et al, 2011b)
Observation Sensitivity Q: How does the analysis change if the observations change? Q: How does the analysis change if the observation array changes?
Impact vs Sensitivity Single 4D-Var cycle obs Impact on 37N transport Sensitivity of 37N transport to removing observations
Observation Sensitivity: An OSSE 4D-Var change in 37N transport when all SSH removed Change in 37N transport predicted by obs sensitivity Change in 37N transport predicted by obs impact
Sequential Data Assimilation: July 2002-Dec2004 Observations Observations Observations prior prior prior Data Assimilation Data Assimilation Data Assimilation Posterior Posterior Posterior 7 days Forecast Forecast Forecast
Observations and Predictability Q: How does each obs contribute to the forecast predictability?
Forecast Ensembles Expected uncertainty of analysis Small spread Predictable Large spread Unpredictable
Predictability Forecast ensemble t0 t0-7d t0+14d t0+7d Analysis cycle ending t0 Forecast cycle ending t0+14 Analysis cycle ending t0+7d Forecast cycle ending t0+14
Predictability positive impact of obs on predictability Alongshore transport Cross-shore transport
Summary • In situ observations have a large impact on • circulation estimates, despite small number. • Adjoint operators provide considerable utility • for quantifying the impact and value of ocean • observations. • Routine monitoring of adjoint-based • diagnostics → real-time monitoring of • observing array. • Quantification of the true value of observations.