1 / 24

Advantages of data assimilation in coastal ocean circulation models: Oregon perspective

Advantages of data assimilation in coastal ocean circulation models: Oregon perspective Alexander L. Kurapov, J. S. Allen, G. D. Egbert, R. N. Miller COAS/Oregon State University In cooperation with P. M. Kosro, M. D. Levine, T. Boyd, J. A. Barth, J. N. Moum, P. T. Strub, S. Erofeeva.

avari
Download Presentation

Advantages of data assimilation in coastal ocean circulation models: Oregon perspective

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Advantages of data assimilation in coastal ocean circulation models: Oregon perspective Alexander L. Kurapov, J. S. Allen, G. D. Egbert, R. N. Miller COAS/Oregon State University In cooperation with P. M. Kosro, M. D. Levine, T. Boyd, J. A. Barth, J. N. Moum, P. T. Strub, S. Erofeeva 29 January 2004, AGU/Ocean Sciences http://www.coas.oregonstate.edu/po/research/kurapov/main.html

  2. HF radars Moorings (ADP, T, S) currents: 3D+time density:3D+time Summer circulation on the Oregon shelf: • wind stress (upwelling favorable) is dominant forcing • strong effects of flow-topography interactions • energetic internal tide Summer 2001: DA system is implemented with data from COAST observational program • Data assimilation: • improves prediction of the ocean state, • provides solution error estimates, • is used as a tool for data synthesis, • helps to design an observational system (e.g., suggests optimal observational locations)

  3. Dual approach: • Objectives: • to develop practical, but still nearly optimal methods for the assimilation of data into coastal circulation models • to apply these methods to measurements from the Oregon shelf • to utilize DA to increase scientific understanding of shelf circulation

  4. Model of of M2 internal tide [Kurapov et al., JPO33, 2003] - linearized, primitive eqns, 3D, periodic in time [~exp(iwt)] - terrain following coordinates e.g., momentum equations: • Model domain: 40  60 km, Dx=1 km, 21 s-layers • Zone of coverage of 2 HF radars (May-July 1998) • Efficient model solver (direct factorization of the model operator) • Address open boundary issues HF (P. M. Kosro) ADP HF Most internal tide comes from outside the computational domain DA: corrects open boundary baroclinic flux

  5. Generalized Inverse Method (GIM): • Solution minimizes a cost function: • Cost Function = || Model error ||2 + || BCond error ||2 + || Obs error ||2 min • Explicit statistical assumptions about errors in the inputs • Statistics in the output (prior model and inverse solutions) are computed [Bennett, 1992, 2002] • State vector: v = {velocity, sea surface elevation, density} • Model+BCond: S v = f + em • Data: L v = d + ed errors in model forcing and data specified prior to assimilation

  6. Use of Representers: Model+BCond: S v = f + em Data: L v = d + ed vo Strongly constrained dynamics: Fwd solver Adjoint solver HF radars: K=900 locations where radial velocity components are available • Reduce burden of representer computation with: • reduced basis representer approach • indirect representer approach [Egbert et al., JGR, 1994] Standard feature in Inverse Ocean Modeling system [IOM, Chua and Bennett, Ocean Modeling, 2001]

  7. Solution sensitivity to the choice of model error covariance COB (in an experiment with synthetic data) • -”true” solution: forced at open boundary (OB) with a significantly baroclinic flux • synthetic data (velocity harmonic constants) are sampled from true solution • prior model: forced with depth-averaged OB current • DA: corrects OB baroclinic fluxes Depth-ave RMS error with respect to true solution DA, COB(Type I) DA, COB(Type II) Prior these two solutions allow for OB b/clinic correction of the same magnitude (but different correlation structure)

  8. DA COB (Type I) is obtained by nesting approach: representer  column of prior solution error covariance matrix • In a large domain, compute representers for small domain boundary data • then sample these representers along the OB of small domain • COB (covariance for the errors on the OB of the small domain, witha dynamically consistent spatial structure) • COBcontrols radiation at an open boundary COB (Type II): our best guess w/out nesting

  9. A series of M2 tidal solutions, May-July 1998 Internal tide intermittence: analysis in 2-week overlapping time windows DA: in each time window depth-ave (rotating CCW) Validation ADP Tidal ellipses of horizontal currents at ADP location, vs depth: (a) observed, (b) prior model, (c) DA. deviations from depth-ave. (CW) No DA ADP DA solution Assimilation of HF surface currents improves prediction at depth

  10. M2 tidal ellipses on the surface: internal tide velocities can be twice as large as barotropic tidal velocities Deviations from depth-ave (time window centered on day 139) Depth-ave CCW rotation CW rotation

  11. day, 1998 Energy balance is closed : Data assimilation corrects only boundary inputs Some persistent features are found: e.g., baroclinic phase and energy propagation is from NW. Most baroclinic signal comes into the computational domain from outside 40 W m-1 Terms in the baroclinic energy equation (time and space averaged) Baroclinic energy flux (depth-integrated and time-ave.)

  12. Baroclinic KE averaged over a series of days 139-167: a) surface, b) bottom, c) cross-section north of Stonewall Bank, d) cross-section through Stonewall Bank. Zones of higher KE variability are aligned along the coast, consistent with energetic of a internal Poincare wave interaction with bathymetry beams over Stonewall B Dominance of 1st baroclinic mode A series of tidal solutions (constrained by HF radar data) provides a uniquely detailed description of spatial and temporal variability of M2 internal tide

  13. AVHRR SST, oC [courtesy P.T. Strub] Model of wind-driven circulation: • Princeton Ocean Model: 220350 km, • periodic OB conditions (south-north), • Dx~2 km, 31 s-layers • Forcing: alongshore wind stress, heat flux • Data assimilation: Optimal Interpolation • Initial implementation (summer 1998): assimilation ofHF radar data improves modeled circulation at depth [Oke et al., JGR-Oceans, 2002] • Data from COAST program (summer 2001): assimilate moored ADP currents

  14. Optimal Interpolation (3DVAR): matrix matching observations to state vector model w/out DA ||Error|| forecast DA Incremental approach: correction is applied gradually over the analysis time window (1/4 of inertial period) analysis Time Forecast error covariance (stationary in OI):Pf= Pm  F(laggedPm, Cd) where Pm is the covariance of errors in the model solution not constrained by the data (in contrast, Pfis conditioned upon previously assimilated data) [Kurapov et al., Mon. Wea Rev., 2002] Pfhas a shorter horizontal scale in the alongshore direction thanPm (effect of propagation) Pm:could be obtained as representer calculation, if an adjoint model were available Presently, Pm is computed from an ensemble of model solutions

  15. Spatial structure of Pf: NMS, 12m SSB, 16m [cm2 s-2]

  16. Time- and depth-ave terms in the momentum eqn. (along-jet direction) no DA DA (ADPs in south) Smooth, large scale correction (in this case, DA tends to reduce upwelling intensity) Dominant dynamical balance is preserved

  17. Assimilation of moored ADP velocities (May-Aug 2001): • Study is focused on: • Distant effect of data assimilation • Multivariate capabilities (effect on SSH, isopycnals, temperature, salinity transport, turbulent dissipation rate) 90 km Central part of model domain with mooring locations, Bathymetry each 100 m (black) and 10 m (half-tone, from 0 to 200 m) Moorings: Lines N and S – COAST (Kosro, Levine, Boyd), NH10 – GLOBEC (Kosro)

  18. Case 1: assimilate currents at Northern Line  improve currents at NH10, SSB Correction can be advected by a predominantly southward current Alongshore depth-ave current: obs, no DA, DA rmse: 7.85.8 cm s1, corr:0.18  0.71 NH10 90 km rmse: 9.67.1 cm s1, corr:0.36  0.70 SSB ADP sites, May-Aug 2001 Assimilated ADP sites Sites where DA is better than model only solution (smaller model-data rms error, larger correlation)

  19. ADP sites, May-Aug 2001 Assimilated ADP sites Sites where DA is better than model only solution (smaller model-data rms error, larger correlation) Case 2: assimilate ADP currents at Southern Line  improve currents up North Correction can be propagated northward with coastal trapped waves Alongshore depth-ave current:obs, no DA, DA NMS rmse: 11.37.9 cm s1, corr:0.46  0.79 NH10 rmse: 7.86.9 cm s1, corr:0.18  0.63

  20. Assimilated site DA is better than model only solution DA is worse than model only solution Posterior error statistics analysis E.g., compare expected and actual analysis rms error as a consistency test for Pf Expected performance diag (Pm) and (Pa) are compared, where Pa = Pf – G H Pf is the analysis error covariance Actual performance Discrepancy between expected and actual outcome when assimilating inner-shelf data : artificially large decorrelation scale in Pf  inclusion of a more realistic spatially varying wind stress is a necessity

  21. Multivariate capabilities SSH: obs, model only, DA (Lines N+S) e.g., effect on SSH (validation - tide gauge data): Model-data Corr.: 0.51  0.78, rmse: 5.4  3.8 cm effect on isopycnal slope: SeaSoar measurements (Barth et al.) no DA DA (South) (white contours are measuredsq= 24, 25, and 26 kg m-3) + improvement in temperature correlations, surface salinity transport

  22. Turbulent Dissipation rate (e): Microstructure data [J.Moum, A. Perlin] No DA DA (North) 12 transects on Line N Time series ofeaveraged near bottom (in box area) yearday, 2001 DA correction in near-bottom velocity field yields improvement in e Analysis of BBL dynamics is extended for the whole study period – presentation OS52I-08

  23. SUMMARY: • Progress has been made on both aspects of the dual approach to coastal ocean DA • Linearized dynamics, variational DA (internal tides) • has provided unique information on spatial and temporal variability of internal tide from HF radar measurements of surface currents • has given us experience in open boundary DA • Nonlinear dynamics, sequential OI DA (wind-driven circulation) • - has shown the value of assimilation of currents from HF radar and from moored ADPs (distant effect, multivariate capabilities, BBL analysis) • has provided information on optimal ADP mooring locations and on effective alongshore scales of ADP current measurements • In both cases, formulation of error hypotheses isthe science and art of DA • DA is utilized to increase scientific understanding of shelf circulation http://www.coas.oregonstate.edu/po/research/kurapov/main.html

  24. PLANNED RESEARCH: • Merger of approaches: use tangent linear and adjoint codes for a fully non-linear ocean circulation model (ROMS) • Use data assimilation to help provide open boundary conditions for high-resolution limited-area coastal models • Tidal research: study effect of wind-forced subinertial flows on internal tide propagation • Study of wind-forced upwelling circulation: analyze cross-shelf transport, bottom boundary layer processes, dynamical balances

More Related