Implementation of a Second-Order TVD Transport Algorithm in the General Ocean Model (GOM)

This study presents the implementation of a scalar transport algorithm in the recently developed General Ocean Model (GOM), a three-dimensional, unstructured grid, finite volume/finite difference model. Solving the advection–diffusion transport equation is an essential part of any ocean circulation model since the baroclinic density gradient distinguishes saline water from freshwater. To achieve both high accuracy and computational efficiency, we adopted a second-order semi-implicit Total Variation Diminishing (TVD) scheme. The TVD approach, known for its ability to suppress non-physical oscillations near steep gradients, provides a higher-fidelity representation of salinity fronts without introducing significant numerical artifacts. The TVD algorithm is constructed with the first-order Upwind scheme, which is known for suffering from excessive numerical diffusion, and the higher-order anti-diffusive flux term. The implemented transport algorithm is evaluated using two standard test cases, an ideal lock exchange problem and a U-shaped channel problem, and it is further applied to simulate salinity dynamics in Mobile Bay, Alabama. The model results from both the first-order Upwind and second-order TVD schemes are compared. The results indicate that the TVD scheme marginally improves the resolution of salinity fronts while maintaining computational stability and efficiency. The implementation enables a flexible and straightforward transition between the first-order scheme, which is faster than the second-order scheme, and the second-order scheme, which is less diffusive than the first-order scheme, enhancing the GOM’s capability for realistic and efficient salinity simulations in a tidally driven estuarine system.