Abstract
The sigma (SIG) coordinate system in ocean circulation simulation models results inevitably in horizontal pressure gradient error. This problem also emerges in models of deep lakes or reservoirs with the same characteristics of underwater terrain mutation. SIG coordinates reflect vertical relative stratification but cannot be used to calculate horizontal pressure gradient force in places with drastic topographic changes; this results in vertical water temperature and circulation errors. In deep lakes or reservoirs, differences in water density caused by the temperature difference between upper and lower water bodies is the primary cause of thermal stratification phenomena. Lake Mead was used as a case study on steep topography based on Environmental Fluid Dynamics Code (EFDC) model in this study. SIG coordinates result in close agreement between the calibrated temperature time series at the top and middle water layers, but disparity in the bottom water layer. The error emerges in the horizontal pressure gradient error due to the SIG coordinate transformation. Neither increasing the vertical resolution nor adjusting the horizontal viscosity coefficient resolve this error. We test the sigma-zed (SGZ) coordinate which combines Z coordinate and SIG coordinate as a replacement for the SIG coordinate to find that they effectively reduce the model’s runtime and simulation efficiency. The vertical temperature distribution in SGZ coordinate mode is more accurate than the distribution in SIG coordinate mode. The Navier-Stokes horizontal gradient and advection diffusion equation results under SIG coordinates are very sensitive to the pressure gradient. The replacement also enhances resolution near the thermocline, facilitates reclosing of the water bottom and the equal sigma surface, lends significant advantages in terms of vertical temperature in the simulation for local deep water with steep terrain, and shortens runtime for 0.14 h. SGZ mixed coordinates are recommended in the simulation of deep lakes or reservoirs wherein the underwater topography is large (with abundant continuous deep trenches or reefs).
Similar content being viewed by others
References
Arifin RR, James SC, Pitts DADA, Hamlet AF, Sharma A, Fernando HJS (2016) Simulating the thermal behavior in Lake Ontario using EFDC. J Great Lakes Res 42(3):511–523
Barron CN, Kara AB, Martin PJ, Rhodes RC, Smedstad LF (2006) Formulation, implementation and examination of vertical coordinate choices in the Global Navy Coastal Ocean Model (NCOM). Ocean Model 11(3):347–375
Berntsen J, Thiem Ø, Avlesen H (2015) Internal pressure gradient errors in σ -coordinate ocean models in high resolution fjord studies. Ocean Model 92:42–55
Eckermann S (2008) Hybrid σ-p coordinate choices for a global model. Mon Weather Rev 137(1)
Ezer T (2005) Entrainment, diapycnal mixing and transport in three-dimensional bottom gravity current simulations using the Mellor–Yamada turbulence scheme. Ocean Model 9(2):151–168
Halliwell GR (2004) Evaluation of vertical coordinate and vertical mixing algorithms in the HYbrid-Coordinate Ocean Model (HYCOM). Ocean Model 7(3–4):285–322
Hamrick JM(1992) A three-dimensional environmental fluid dynamics computer code: theoretical and computational aspects. The College ofWilliam and Mary, Virginia Institute of Marine Science, Special Report 317, 63
Haney RL (1991) On the pressure gradient force over steep topography in sigma coordinate ocean models. J Phys Oceanogr 21(4):610–618
Hayat T, Shehzad SA, Al-Sulami HH, Asghar S (2013) Influence of thermal stratification on the radiative flow of Maxwell fluid. J Braz Soc Mech Sci 35(4):381–389
Hodges BR, Laval B, Wadzuk BM (2006) Numerical error assessment and a temporal horizon for internal waves in a hydrostatic model. Ocean Model 13(1):44–64
Huu Chung D, Craig PM (2014) An algorithm on vertical layering in EFDC with sigma-zed.Proceedings of the 11th International Conference on Hydrodynamics (ICHD 2014) October 19-24, 2014 Singapore
Laval B, Hodges BR, Imberger J (2003) Reducing numerical diffusion effects with pycnocline filter. J Hydraul Eng 129(3):215–224
Li YP, Acharya K, Chen D, Stone M, Tietjen TE, Holdren GC (2010) Modeling water ages and thermal structure of lake mead under changing water levels. Lake Reserv Manage 26(4):258–272
Li YP, Acharya K, Yu Z (2011) Modeling impacts of Yangtze River water transfer on water ages in Lake Taihu, China. Ecol Eng 37(2):325–334
Li Y, Wang D, Wang B (2012) A new approach to implement sigma coordinate in a numerical model. Commun Comput Phys 12(4):1033–1050
Li Y, Huang TL, Zhou ZZ, Long SH, Zhang HH (2018) Effects of reservoir operation and climate change on thermal stratification of a canyon-shaped reservoir, in Northwest China. Water Sci Tech Water Supply 18(2):ws2017068
Liu X, Huang W (2008) An effective algorithm to reduce horizontal pressure gradient errors in σ-coordinate in EFDC hydrodynamic model. J Coast Res 10052:193–204
Marti C L, Imberger J, Garibaldi L, Leoni B (2016) Using time scales to characterize phytoplankton assemblages in a deep subalpine lake during the thermal stratification period: Lake Iseo, Italy Water Resour Res 52 (3):n/a-n/a
Mellor GL, Ezer T, Oey LY (1994) The pressure gradient conundrum of sigma coordinate ocean models. J Atmos Ocean Technol 11(4):1126–1134
Phillips NA (1957) A coordinate system having some special advantages for numerical forecasting. J Atmos Sci 14(2):184–185
Preusse M, Freistuehler H, Peeters F (2012) Seasonal variation of solitary wave properties in Lake Constance. J Geophys Res Oceans 117(C04026)
Sahoo GB, Forrest AL, Schladow SG, Reuter JE, Coats R, Dettinger M (2016) Climate change impacts on lake thermal dynamics and ecosystem vulnerabilities. Limnol Oceanogr 61(2):496–507
Smagorinsky J (1963) General circulation experiments with the primitive equations. Mon Weather Rev 91(3):99–164
Twichell DC, Cross VAA, Belew SD (2003) Mapping the floor of Lake Mead (Nevada and Arizona): preliminary discussion and GIS data release
Vilhena LC, Marti CL, Imberger J (2013) The importance of nonlinear internal waves in a deep subalpine lake: Lake Iseo, Italy. Limnol Oceanogr 58(5):1871–1891
Winder M, Schindler DE (2004) Climate change uncouples trophic interactions in an aquatic ecosystem. Ecology 85(8):2100–2106
Acknowledgments
We would also like to thank Ronald Veley with the US Geological Survey (USGS) for providing monitoring data, and Paul Craig for valuable assistance in developing the Lake Mead model.
Funding
This research was supported by the National Key Research and Development Program of China (2017YFC0405203), the Fundamental Research Funds for the Central Universities, and the World Class Universities (Disciplines), the Characteristic Development Guidance Funds for the Central Universities, the Fundamental Research Funds for the Central Universities (2018B48214, 2017B20514) and PAPD, and the Chinese National Science Foundation (51779072, 51579071, 41323001, 51539003), National Science Funds for Creative Research Groups of China (51421006).
Author information
Authors and Affiliations
Corresponding author
Additional information
Responsible editor: Marcus Schulz
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Shi, Y., Li, Y., Cheng, Y. et al. A novel sigma coordinate system to simulate abrupt changes of underwater terrain in a hydrodynamic model: application to Lake Mead, USA. Environ Sci Pollut Res 27, 9261–9273 (2020). https://doi.org/10.1007/s11356-019-06994-y
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11356-019-06994-y