Search

Search Results

• 1. [Article] On Intermediate Models for Barotropic Continental Shelf and Slope Flow Fields. Part II: Comparison of Numerical Model Solutions in Doubly Periodic Domains

  As part of a program to improve understanding of the dynamics of the complicated, vigorous eddy and jet flow fields recently observed over the continental shelf and slope, we investigate the potential ...

  Citation × Citation Citation Abstract Copy Title: On Intermediate Models for Barotropic Continental Shelf and Slope Flow Fields. Part II: Comparison of Numerical Model Solutions in Doubly Periodic Domains Author: Newberger, P. A., Allen, J. S., Barth, J. A. As part of a program to improve understanding of the dynamics of the complicated, vigorous eddy and jet flow fields recently observed over the continental shelf and slope, we investigate the potential of intermediate models for use in both process and data assimilation studies of these flows. Intermediate models incorporate physics simpler than that contained in the full primitive equations yet more complete than in the quasi-geostrophic equations, and are capable of representing subinertial flows over O(1) bottom topographic variations and/or with O(1) isopycnal slopes. In addition, intermediate models dynamically filter out high-frequency gravity- inertial motions leading, potentially, to higher computational efficiency and well-posed limited area forecast/hindcast models. Initial studies focus on single layer flows on an f-plane with a free surface, governed by the shallow-water equations. In Part I, various intermediate models are formulated and their accuracy assessed by comparing some exact nonlinear analytical solutions that exist for the shallow-water equations with corresponding analytical solutions of the intermediate models. Here in Part II, an extensive set of numerical finite-difference solutions to initial-value problems in doubly periodic domains (to isolate model differences from the influence of boundary condition implementation on solid walls) is used to determine the accuracy of various intermediate models by comparing their predictions with those of a shallow-water equation model that uses a potential enstrophy and energy conserving numerical scheme (SWE). Intermediate model results are also contrasted with those from a quasi-geostrophic (QG) model. The intermediate models considered are based on the geostrophic momentum (GM) approximation, the derivation of Salmon utilizing Hamilton's principle (HP), a geostrophic vorticity (GV) approximation, a combination of the quasi-geostrophic momentum and full continuity equations (IM), the linear balance equations (LBE), the balance equations (BE), the related balance-type (HBE, BEM, NBE) and modified linear balance equations (LQBE), and on Lynch's slow equations in their original form (SE) and in a modified form (MSE). In addition, a semi-implicit version (SEMI) of the shallow-water equations, which numerically filters high-frequency motions, is included in the study. The basic initial-value problem used to test the various intermediate models involves sinusoidal flow over a symmetric Gaussian-shaped bottom topographic feature. Comparisons are made for a range of the relevant dimensionless model parameters including the strength of the flow (as measured by the Rossby number), the square of the ratio of a characteristic horizontal length scale to the Rossby radius of deformation, and the height of the topographic feature. A second initial-value problem involves the evolution of a rotating elliptical vortex over both flat and variable bottom topography. Results show that in cases with low local Rossby number flow, but with large topographic height, most of the intermediate models are substantially more accurate than QG. Even for flows with O(1) local Rossby numbers, some of the intermediate models continue to give excellent results. Specifically, BE, BEM and SEMI consistently do the best in the comparisons. Although the relative ordering of the remaining models is somewhat parameter-dependent, in general the next most accurate models are MSE and LQBE followed by HP and NBE. These are followed in quality by GV and GM while the remaining models, HBE, IM and LBE, perform least accurately for the range of parameters studied. Generally, intermediate models that have an integral invariant corresponding to conservation of potential enstrophy do better than those without, with the best results coming from models which have this property and a nonlinear balance equation. Title: On Intermediate Models for Barotropic Continental Shelf and Slope Flow Fields. Part II: Comparison of Numerical Model Solutions in Doubly Periodic Domains Author: Newberger, P. A., Allen, J. S., Barth, J. A. As part of a program to improve understanding of the dynamics of the complicated, vigorous eddy and jet flow fields recently observed over the continental shelf and slope, we investigate the potential of intermediate models for use in both process and data assimilation studies of these flows. Intermediate models incorporate physics simpler than that contained in the full primitive equations yet more complete than in the quasi-geostrophic equations, and are capable of representing subinertial flows over O(1) bottom topographic variations and/or with O(1) isopycnal slopes. In addition, intermediate models dynamically filter out high-frequency gravity- inertial motions leading, potentially, to higher computational efficiency and well-posed limited area forecast/hindcast models. Initial studies focus on single layer flows on an f-plane with a free surface, governed by the shallow-water equations. In Part I, various intermediate models are formulated and their accuracy assessed by comparing some exact nonlinear analytical solutions that exist for the shallow-water equations with corresponding analytical solutions of the intermediate models. Here in Part II, an extensive set of numerical finite-difference solutions to initial-value problems in doubly periodic domains (to isolate model differences from the influence of boundary condition implementation on solid walls) is used to determine the accuracy of various intermediate models by comparing their predictions with those of a shallow-water equation model that uses a potential enstrophy and energy conserving numerical scheme (SWE). Intermediate model results are also contrasted with those from a quasi-geostrophic (QG) model. The intermediate models considered are based on the geostrophic momentum (GM) approximation, the derivation of Salmon utilizing Hamilton's principle (HP), a geostrophic vorticity (GV) approximation, a combination of the quasi-geostrophic momentum and full continuity equations (IM), the linear balance equations (LBE), the balance equations (BE), the related balance-type (HBE, BEM, NBE) and modified linear balance equations (LQBE), and on Lynch's slow equations in their original form (SE) and in a modified form (MSE). In addition, a semi-implicit version (SEMI) of the shallow-water equations, which numerically filters high-frequency motions, is included in the study. The basic initial-value problem used to test the various intermediate models involves sinusoidal flow over a symmetric Gaussian-shaped bottom topographic feature. Comparisons are made for a range of the relevant dimensionless model parameters including the strength of the flow (as measured by the Rossby number), the square of the ratio of a characteristic horizontal length scale to the Rossby radius of deformation, and the height of the topographic feature. A second initial-value problem involves the evolution of a rotating elliptical vortex over both flat and variable bottom topography. Results show that in cases with low local Rossby number flow, but with large topographic height, most of the intermediate models are substantially more accurate than QG. Even for flows with O(1) local Rossby numbers, some of the intermediate models continue to give excellent results. Specifically, BE, BEM and SEMI consistently do the best in the comparisons. Although the relative ordering of the remaining models is somewhat parameter-dependent, in general the next most accurate models are MSE and LQBE followed by HP and NBE. These are followed in quality by GV and GM while the remaining models, HBE, IM and LBE, perform least accurately for the range of parameters studied. Generally, intermediate models that have an integral invariant corresponding to conservation of potential enstrophy do better than those without, with the best results coming from models which have this property and a nonlinear balance equation. Close

• 2. [Article] On Intermediate Models for Barotropic Continental Shelf and Slope Flow Fields. Part III: Comparison of Numerical Model Solutions in Periodic Channels

  The study of intermediate models for barotropic continental shelf and slope flow fields initiated in Parts I and II is continued. The objective is to investigate the possible use of intermediate models ...

  Citation × Citation Citation Abstract Copy Title: On Intermediate Models for Barotropic Continental Shelf and Slope Flow Fields. Part III: Comparison of Numerical Model Solutions in Periodic Channels Author: Newberger, P. A., Barth, J. A., Allen, J. S. The study of intermediate models for barotropic continental shelf and slope flow fields initiated in Parts I and II is continued. The objective is to investigate the possible use of intermediate models for process and data assimilation studies of nonlinear mesoscale eddy and jet current fields over the continental shelf and slope. Intermediate models contain physics between that in the primitive equations and that in the quasi-geostrophic equations and are capable of representing subinertial frequency motion over the O(1) topographic variations typical of the continental slope while filtering out high-frequency gravity–inertial waves. We concentrate on the application of intermediate models to the f-plane shallow-water equations. The accuracy of several intermediate models is evaluated here by a comparison of numerical finite-difference solutions with those of the primitive shallow-water equations (SWE) and with those of the quasi-geostrophic equations (QG) for flow in a periodic channel. The intermediate models that we consider are based on the balance equations (BE), the balance equations derived from momentum equations (BEM), the potential vorticity conserving linear balance equations (LQBE), the hybrid balance equations (HBE), the near balance equation (NBE), a geostrophic vorticity (GV) approximation, the geostrophic momentum (GM) approximation, and the quasi-geostrophic momentum and full continuity equations (IM). The periodic channel provides a basic geometry for the study of physical flow processes over the continental shelf and slope. Wall boundary conditions are formulated for the intermediate models and implemented in the numerical finite-difference approximations. The ability of intermediate models to represent linear ageostrophic coastally trapped waves, i.e., Kelvin and continental shelf waves, is verified by numerical experiments. The results of numerical solution intercomparisons for initial-value problems involving O(1) topographic variations are as follows. For flow at small local Rossby number |ϵ[subscript]L| < 0.2, where ϵ[subscript]L is given by the magnitude of the vorticity divided by f, all of the intermediate models do well, while the QG model does poorly. For flows with larger values of |ϵ[subscript]L|, e.g., |ϵ[subscript]L| ≈ 0.5, the performance of the different intermediate models varies. BEM and BE consistently give extremely accurate solutions while the solutions from LQBE are almost as good. The other models are substantially less accurate with errors generally increasing in the order NBE, HBE, GV, GM, IM. The QG solution always has the largest errors. Consistent with the results from the studies in Part II in a doubly periodic domain, the balance equations BE and BEM, followed closely by LQBE, appear to be the most accurate intermediate model Title: On Intermediate Models for Barotropic Continental Shelf and Slope Flow Fields. Part III: Comparison of Numerical Model Solutions in Periodic Channels Author: Newberger, P. A., Barth, J. A., Allen, J. S. The study of intermediate models for barotropic continental shelf and slope flow fields initiated in Parts I and II is continued. The objective is to investigate the possible use of intermediate models for process and data assimilation studies of nonlinear mesoscale eddy and jet current fields over the continental shelf and slope. Intermediate models contain physics between that in the primitive equations and that in the quasi-geostrophic equations and are capable of representing subinertial frequency motion over the O(1) topographic variations typical of the continental slope while filtering out high-frequency gravity–inertial waves. We concentrate on the application of intermediate models to the f-plane shallow-water equations. The accuracy of several intermediate models is evaluated here by a comparison of numerical finite-difference solutions with those of the primitive shallow-water equations (SWE) and with those of the quasi-geostrophic equations (QG) for flow in a periodic channel. The intermediate models that we consider are based on the balance equations (BE), the balance equations derived from momentum equations (BEM), the potential vorticity conserving linear balance equations (LQBE), the hybrid balance equations (HBE), the near balance equation (NBE), a geostrophic vorticity (GV) approximation, the geostrophic momentum (GM) approximation, and the quasi-geostrophic momentum and full continuity equations (IM). The periodic channel provides a basic geometry for the study of physical flow processes over the continental shelf and slope. Wall boundary conditions are formulated for the intermediate models and implemented in the numerical finite-difference approximations. The ability of intermediate models to represent linear ageostrophic coastally trapped waves, i.e., Kelvin and continental shelf waves, is verified by numerical experiments. The results of numerical solution intercomparisons for initial-value problems involving O(1) topographic variations are as follows. For flow at small local Rossby number |ϵ[subscript]L| < 0.2, where ϵ[subscript]L is given by the magnitude of the vorticity divided by f, all of the intermediate models do well, while the QG model does poorly. For flows with larger values of |ϵ[subscript]L|, e.g., |ϵ[subscript]L| ≈ 0.5, the performance of the different intermediate models varies. BEM and BE consistently give extremely accurate solutions while the solutions from LQBE are almost as good. The other models are substantially less accurate with errors generally increasing in the order NBE, HBE, GV, GM, IM. The QG solution always has the largest errors. Consistent with the results from the studies in Part II in a doubly periodic domain, the balance equations BE and BEM, followed closely by LQBE, appear to be the most accurate intermediate model Close