• Bennett, A. F., 1985: Array design by inverse methods. Progress in Oceanography, Vol. 15, Pergamon, 129–156.

    • Crossref
    • Export Citation
  • ——, 1990: Inverse methods for assessing ship-of-opportunity networks and estimating circulation and winds from tropical expendable bathythermograph data. J. Geophys. Res.,95 (C9), 16 111–16 148.

  • ——, 1992: Inverse Methods in Physical Oceanography. Cambridge University Press, 346 pp.

  • ——, and P. C. McIntosh, 1982: Open ocean modeling as an inverse problem: Tidal theory. J. Phys. Oceanogr.,12, 1004–1018.

    • Crossref
    • Export Citation
  • Csanady, G. T., 1978: The arrested topographic wave. J. Phys. Oceanogr.,8, 47–62.

    • Crossref
    • Export Citation
  • Dowd, M. G., and K. R. Thompson, 1996: Extraction of tidal currents from a shipborne ADCP. J. Geophys. Res.,101, 8943–8956.

    • Crossref
    • Export Citation
  • Draper, N. R. 1981: Applied Regression Analysis. John Wiley and Sons, Inc., 736 pp.

  • Foreman, M. G. G., M. L. Delves, I. Barrodale, and R. F. Henry, 1982: On the use of the Proudman–Heaps tidal theorem. Geophys. J. Roy Astron. Soc.,63, 467–478.

    • Crossref
    • Export Citation
  • Gekeler, J, 1995: Assimilating data into a tidal model of the Irish and Celtic Seas. Cont. Shelf Res.,15, 1381–1408.

    • Crossref
    • Export Citation
  • Gilbert, J. R., C. Moler, and R. Schreiber, 1992: Sparse matrices in Matlab: Design and implementation. SIAM J. Matrix Anal.,13, 333–356.

    • Crossref
    • Export Citation
  • Greenberg, D., F. Werner, and D. R. Lynch, 1998: A diagnostic finite element ocean circulation model in spherical-polar coordinates. J. Atmos. Oceanic Technol.,15, 942–958.

    • Crossref
    • Export Citation
  • Han, G., C. G. Hannah, J. W. Loder, and P. C. Smith, 1997: Seasonal variation of the three-dimensional mean circulation over the Scotian Shelf. J. Geophys. Res.,102 (C1), 1011–1025.

    • Crossref
    • Export Citation
  • Hannah, C., and D. G. Wright, 1995: Depth dependent analytical and numerical solutions for wind-driven flow in the coastal ocean. Quantitative Skill Assessment for Coastal Ocean Models, D. R. Lynch and D. R. Davies, Eds., Amer. Geophys. Union, 125–152.

    • Crossref
    • Export Citation
  • ——, J. W. Loder, and D. G. Wright, 1996: Seasonal variation in the baroclinic circulation in the Scotian marine region. Buoyancy Effects on Coastal Dynamics, D. Aubrey, Ed., Amer. Geophys. Union.

  • Lardner, R. W., 1992: Optimal control of open boundary conditions for a numerical tidal model. Comp. Methods Appl. Mech. Eng.,102, 367–387.

    • Crossref
    • Export Citation
  • Le Dimet, F. X., and O. Talagrand, 1986: Variational algorithms for analysis and assimilation of meteorological observations. Tellus,38, 97–110.

    • Crossref
    • Export Citation
  • Lewis, J. M., and J. C. Derber, 1985: The use of adjoint equations to solve a variational adjustment problem with advective constraints. Tellus,37, 309–322.

    • Crossref
    • Export Citation
  • Lough, R. G., and Coauthors, 1994: Influence of wind driven advection on the interannual variability in cod egg and larval distribution on Georges Bank: 1982 vs 1985. ICES Mar. Sci. Symp.,198, 356–378.

  • Lynch, D. R., and F. E. Werner, 1987: Three-dimensional hydrodynamics on finite elements, Part I: Linearized harmonic model. Int. J. Numer. Methods Fluids,7, 871–909.

    • Crossref
    • Export Citation
  • ——, and C. E. Naimie, 1993: The m2 tide and its residual on the outer banks of the Gulf of Maine. J. Phys. Oceanogr.,23, 2222–2253.

    • Crossref
    • Export Citation
  • ——, F. E. Werner, D. Greenberg, and J. W. Loder, 1992: Diagnostic model for baroclinic, wind-driven and tidal circulation in shallow seas. Cont. Shelf Res.,12, 37–64.

    • Crossref
    • Export Citation
  • ——, J. Ip, C. Naimie, and F. Werner, 1996: Comprehensive coastal circulation model with application to the Gulf of Maine. Cont. Shelf Res.,16, 875–906.

    • Crossref
    • Export Citation
  • Mathworks, 1992: Matlab Reference Guide 548 pp.

  • McIntosh, P. C., and A. F. Bennett, 1984: Open ocean modeling as an inverse problem: M2 tides in bass strait. J. Phys. Oceanogr.,14, 601–614.

    • Crossref
    • Export Citation
  • Naimie, C. E., and D. R. Lynch, 1993: FUNDY5 User’s Manual. Thayer School of Engineering, 40 pp.

  • Pedlosky, J., 1979: Geophysical Fluid Dynamics. Springer-Verlag, 624 pp.

    • Crossref
    • Export Citation
  • Sasaki, Y., 1955: A variational study of the numerical prediction based on the variational principle. J. Meteor. Soc. Japan,33, 262–275.

    • Crossref
    • Export Citation
  • ——, 1970: Some basic formalisms in numerical variational analysis. Mon. Wea. Rev.,98, 875–883.

    • Crossref
    • Export Citation
  • Seber, G. A. F., 1977: Linear Regression Analysis. John Wiley and Sons, 496 pp.

  • Seiler, U., 1993: Estimation of open boundary conditions with the adjoint method. J. Geophys. Res.,98 (C12), 22 855–22 870.

    • Crossref
    • Export Citation
  • Thacker, W. C., 1988: Three lectures of fitting numerical models to observations. GKSS-Forschungszentrum Geesthacht GmbH Geethacht, GKSS 87/E/65, 64 pp.

  • ——, and R. B. Long, 1988: Fitting dynamics to data. J. Geophys. Res.,93, 1227–1240.

    • Crossref
    • Export Citation
  • Tremblay, M. J., J. W. Loder, F. E. Werner, C. E. Naimie, F. H. Page, and M. M. Sinclair, 1994: Drift of sea scallop larval placopecten magellanicus on Georges Bank: A model study of the roles of mean advection, larval origin. Deep-Sea Res.,41, 7–49.

    • Crossref
    • Export Citation
  • Werner, F. E., F. H. Page, D. R. Lynch, J. W. Loder, R. G. Lough, R. I. Perry, and D. A. Greenberg, 1993: Influence of mean advection and simple behavior on the distribution of cod and haddock early life stages on Georges Bank. Fish. Oceanogr.,2, 184–193.

    • Crossref
    • Export Citation
  • Wright, D. G., and J. W. Loder, 1988: On the influences on nonlinear bottom friction on the topographic rectification of tidal currents. Geophys. Astrophys. Fluid Dyn.,42, 227–245.

    • Crossref
    • Export Citation
  • Wunsch, C., 1996: The Ocean Circulation Inverse Problem. Cambridge University Press, 442 pp.

    • Crossref
    • Export Citation
  • Xu, Z., 1995: A transport approach to the convolution method for numerical modeling of linearized 3-D circulation. Int. J. Numer. Methods Fluids,20, 363–391.

    • Crossref
    • Export Citation
  • Zahel, W., 1991: Modeling ocean tides with and without assimilating data. J. Geophys. Res.,96, (20), 379–391.

    • Crossref
    • Export Citation
  • ——, 1995: Assimilating ocean tide determined data into global tidal models. J. Mar. Syst.,6, 3–13.

    • Crossref
    • Export Citation
  • Zhu, K., M. Navon, and X. Zou, 1994: Variational data assimilation with a variable resolution finite-element shallow-water equations model. Mon. Wea. Rev.,122, 946–965.

    • Crossref
    • Export Citation
  • Zou, J., W. W. Hsieh, and I. M. Navon, 1995: Sequential open-boundary control by data assimilation in a limited-area model. Mon. Wea. Rev.,123, 2899–2909.

    • Crossref
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 138 138 3
PDF Downloads 10 10 0

A Direct Inverse Method for Inferring Open Boundary Conditions of a Finite-Element Linear Harmonic Ocean Circulation Model

View More View Less
  • 1 Ocean Sciences Division, Fisheries and Oceans Canada, Bedford Institute of Oceanography, Dartmouth, Nova Scotia, Canada
© Get Permissions Rent on DeepDyve
Restricted access

Abstract

A direct inverse method is presented for inferring numerical model open boundary conditions from interior observational data. The dynamical context of the method is the frequency-domain 3D linear shallow water equations. A set of weight matrices is derived via finite-element discretization of the dynamical equations. The weight matrices explicitly express any interior solution of elevations or velocities as a weighted sum of boundary elevations. The interior data assimilation is then cast as a regression problem.

The weight matrix may be singular, which implies there may be an infinite set of boundary conditions that fit the data equally well. With the singular value decomposition technique, a general solution is provided for this infinite set of minimum-squared-misfit boundary conditions. Among them, a particular boundary condition, which minimizes potential energy on the boundary (hence the whole domain), is studied in detail: its confidence interval is defined and a way to smooth it is discussed.

Green’s function maps for the weight matrix provide insights into the dynamics inherent to the model domain. Such maps should be useful for many purposes. One of their usages demonstrated in this paper is to provide a physical explanation for the singularity of the system. Also discussed are how to assess the compatibility between data and the model and how to design a null-space smoothing device for smoothing the potential energy minimum boundary condition. While maintaining the goodness of fit between data and model, smoothing of the boundary condition may improve the interior solutions at the locations where the data have poor control.

The method is tested in a realistic domain but with synthetic data. The test yields very satisfactory results. Application of the method to Lardner’s open bay tidal problem demonstrates its advantages in computational accuracy and inexpensiveness.

Corresponding author address: Dr. Zhigang Xu, Ocean Sciences Division, Fisheries and Ocean Canada, Bedford Institute of Oceanography, P.O. Box 1006, Dartmouth, NS B2Y 4A2 Canada.

Abstract

A direct inverse method is presented for inferring numerical model open boundary conditions from interior observational data. The dynamical context of the method is the frequency-domain 3D linear shallow water equations. A set of weight matrices is derived via finite-element discretization of the dynamical equations. The weight matrices explicitly express any interior solution of elevations or velocities as a weighted sum of boundary elevations. The interior data assimilation is then cast as a regression problem.

The weight matrix may be singular, which implies there may be an infinite set of boundary conditions that fit the data equally well. With the singular value decomposition technique, a general solution is provided for this infinite set of minimum-squared-misfit boundary conditions. Among them, a particular boundary condition, which minimizes potential energy on the boundary (hence the whole domain), is studied in detail: its confidence interval is defined and a way to smooth it is discussed.

Green’s function maps for the weight matrix provide insights into the dynamics inherent to the model domain. Such maps should be useful for many purposes. One of their usages demonstrated in this paper is to provide a physical explanation for the singularity of the system. Also discussed are how to assess the compatibility between data and the model and how to design a null-space smoothing device for smoothing the potential energy minimum boundary condition. While maintaining the goodness of fit between data and model, smoothing of the boundary condition may improve the interior solutions at the locations where the data have poor control.

The method is tested in a realistic domain but with synthetic data. The test yields very satisfactory results. Application of the method to Lardner’s open bay tidal problem demonstrates its advantages in computational accuracy and inexpensiveness.

Corresponding author address: Dr. Zhigang Xu, Ocean Sciences Division, Fisheries and Ocean Canada, Bedford Institute of Oceanography, P.O. Box 1006, Dartmouth, NS B2Y 4A2 Canada.

Save