Are there inescapable issues prohibiting the use of terrain-following coordinates in climate models?
| Autor: | Alexander F. Shchepetkin, James C. McWilliams, M. Jeroen Molemaker, Florian Lemarié, François Colas, Jaison Kurian |
|---|---|
| Přispěvatelé: | Department of Atmospheric Sciences [Los Angeles], University of California [Los Angeles] (UCLA), University of California (UC)-University of California (UC), University of California-University of California |
| Rok vydání: | 2012 |
| Předmět: |
Atmospheric Science
Geopotential 010504 meteorology & atmospheric sciences Meteorology Scale (ratio) Ocean modeling Terrain Oceanography Sigma coordinates K-Profile Parameterization 01 natural sciences Diapycnal mixing [SPI.MECA.MEFL]Engineering Sciences [physics]/Mechanics [physics.med-ph]/Fluids mechanics [physics.class-ph] [MATH.MATH-MP]Mathematics [math]/Mathematical Physics [math-ph] Computer Science (miscellaneous) Applied mathematics [PHYS.MECA.MEFL]Physics [physics]/Mechanics [physics]/Fluid mechanics [physics.class-ph] 14. Life underwater Physics::Atmospheric and Oceanic Physics 0105 earth and related environmental sciences Pacific Ocean 010505 oceanography Advection Geotechnical Engineering and Engineering Geology 13. Climate action Biharmonic equation Climate model Temporal discretization [MATH.MATH-NA]Mathematics [math]/Numerical Analysis [math.NA] Geology Smoothing |
| Zdroj: | Ocean Modelling Ocean Modelling, 2011, 42, pp.57-79. ⟨10.1016/j.ocemod.2011.11.007⟩ Ocean Modelling, Elsevier, 2011, 42, pp.57-79. ⟨10.1016/j.ocemod.2011.11.007⟩ |
| ISSN: | 1463-5003 |
| DOI: | 10.1016/j.ocemod.2011.11.007 |
| Popis: | International audience; The use of $\sigma$-coordinate ocean models has historically been considered a disadvantage for large-scale climate studies. The main reason resides in the non-alignment of the vertical coordinate isosurfaces with either geopotential surfaces or isopycnals making it harder to accurately compute the horizontal pressure gradient, advection, and isoneutral tracer diffusion. Moreover, this class of model requires a vertical mixing parameterization robust to large changes in the vertical resolution between shallow and deep areas. In this paper we show that, with some adjustments of the tracer advection, the surface boundary layer parameterization and the vertical grid, a $\sigma$-coordinate model can achieve an accurate representation of the oceanic interior and mixed-layer dynamics. To do so, a new way of handling the temporal discretization of the rotated biharmonic operator is used to achieve tracer variance dissipation in an adiabatic and computationally efficient way. Furthermore, a redesign of the K-Profile surface layer Parameterization (KPP) to improve the regularity of the solution and the overall numerical efficiency of the scheme is introduced. To validate the new algorithmic developments, we perform a set of coarse-resolution realistic basin-scale Pacific simulations. Besides improving the conservation of water mass properties, the use of an isoneutral tracer hyperdiffusion is shown to have a negative feedback on the circulation error growth rate, thus significantly reducing the sensitivity of the model solution to the degree of topographic smoothing. The overall validation of our simulations, focusing on the key characteristics of the circulation in the Pacific Ocean, provides some evidence of the efficacy of a terrain-following coordinate for large scale applications. |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|
Full Text from ScienceDirect