Jeon, C. - H., Buijsman, M. C., Wallcraft, A. J., Shriver, J. F., Arbic, B. K., Richman, J. G., et al. (2019). Improving surface tidal accuracy through two-way nesting in a global ocean model. Ocean Modelling , 137 , 98–113.
Abstract: In global ocean simulations, forward (non-data-assimilative) tide models generally feature large sea-surface-height errors near Hudson Strait in the North Atlantic Ocean with respect to altimetry-constrained tidal solutions. These errors may be associated with tidal resonances that are not well resolved by the complex coastal-shelf bathymetry in low-resolution simulations. An online two-way nesting framework has been implemented to improve global surface tides in the HYbrid Coordinate Ocean Model (HYCOM). In this framework, a high-resolution child domain, covering Hudson Strait, is coupled with a relatively low-resolution parent domain for computational efficiency. Data such as barotropic pressure and velocity are exchanged between the child and parent domains with the external coupler OASIS3-MCT. The developed nesting framework is validated with semi-idealized basin-scale model simulations. The M2 sea-surface heights show very good accuracy in the one-way and two-way nesting simulations in Hudson Strait, where large tidal elevations are observed. In addition, the mass and tidal energy flux are not adversely impacted at the nesting boundaries in the semi-idealized simulations. In a next step, the nesting framework is applied to a realistic global tide simulation. In this simulation, the resolution of the child domain (1/75°) is three times as fine as that of the parent domain (1/25°). The M2 sea-surface-height root-mean-square errors with tide gauge data and the altimetry-constrained global FES2014 and TPXO9-atlas tidal solutions are evaluated for the nesting and no-nesting solutions. The better resolved coastal bathymetry and the finer grid in the child domain improve the local tides in Hudson Strait and Bay, and the back-effect of the coastal tides induces an improvement of the barotropic tides in the open ocean of the Atlantic.
Nyadjro, E. S., Rydbeck, A. V., Jensen, T. G., Richman, J. G., & Shriver, J. F. (2020). On the Generation and Salinity Impacts of Intraseasonal Westward Jets in the Equatorial Indian Ocean. J. Geophys. Res. Oceans , 125 (6), e2020JC016066.
Abstract: While westerly winds dominate the equatorial Indian Ocean and generate the well‐known eastward flowing Wyrtki Jets during boreal spring and fall, there is evidence of a strong westward surface jet during winter that is swifter than eastward currents during that season. A weaker westward jet is found in summer. In this study, we report the occurrence, characteristics, and intraseasonal variability of this westward jet and its impact on mixed layer salinity in the equatorial Indian Ocean using the HYbrid Coordinate Ocean Model (HYCOM) reanalysis with the Navy Coupled Ocean Data Assimilation (NCODA). The westward jet typically occurs in the upper 50 m, above an eastward flowing equatorial undercurrent, with peak westward volume transport of approximately −8 Sv. The westward jet builds up gradually, decays rapidly, and is primarily forced by local intraseasonal wind stress anomalies generated by atmospheric intraseasonal convection. Westward acceleration of the jet occurs when the dominant intraseasonal westward wind anomaly is not balanced by the zonal pressure gradient (ZPG) force. The intraseasonal westward jet generates strong horizontal advection and is the leading cause of mixed layer freshening in the western equatorial Indian Ocean. Without it, a saltier mixed layer would persist and weaken any barrier layers. Existing barrier layers are strengthened following the passage of freshwater‐laden westward jets. Deceleration of the westward jet occurs when the eastward ZPG becomes increasingly important and the westward intraseasonal wind anomalies weaken. A rapid reversal of atmospheric intraseasonal convection‐driven surface winds eventually terminates the westward jet.
Nelson, A. D., Arbic, B. K., Zaron, E. D., Savage, A. C., Richman, J. G., Buijsman, M. C., et al. (2019). Toward Realistic Nonstationarity of Semidiurnal Baroclinic Tides in a Hydrodynamic Model. J. Geophys. Res. Oceans , 124 (9), 6632–6641.
Abstract: Semidiurnal baroclinic tide sea surface height (SSH) variance and semidiurnal nonstationary variance fraction (SNVF) are compared between a hydrodynamic model and altimetry for the low- to middle-latitude global ocean. Tidal frequencies are aliased by similar to 10-day altimeter sampling, which makes it impossible to unambiguously identify nonstationary tidal signals from the observations. In order to better understand altimeter sampling artifacts, the model was analyzed using its native hourly outputs and by subsampling it in the same manner as altimeters. Different estimates of the semidiurnal nonstationary and total SSH variance are obtained with the model depending on whether they are identified in the frequency domain or wave number domain and depending on the temporal sampling of the model output. Five sources of ambiguity in the interpretation of the altimetry are identified and briefly discussed. When the model and altimetry are analyzed in the same manner, they display qualitatively similar spatial patterns of semidiurnal baroclinic tides. The SNVF typically correlates above 80% at all latitudes between the different analysis methods and above 60% between the model and altimetry. The choice of analysis methodology was found to have a profound effect on estimates of the semidiurnal baroclinic SSH variance with the wave number domain methodology underestimating the semidiurnal nonstationary and total SSH variances by 68% and 66%, respectively. These results produce a SNVF estimate from altimetry that is biased low by a factor of 0.92. This bias is primarily a consequence of the ambiguity in the separation of tidal and mesoscale signals in the wave number domain.
Timko, P. G., Arbic, B. K., Hyder, P., Richman, J. G., Zamudio, L., O'Dea, E., et al. (2019). Assessment of shelf sea tides and tidal mixing fronts in a global ocean model. Ocean Modelling , 136 , 66–84.
Abstract: Tidal mixing fronts, which represent boundaries between stratified and tidally mixed waters, are locations of enhanced biological activity. They occur in summer shelf seas when, in the presence of strong tidal currents, mixing due to bottom friction balances buoyancy production due to seasonal heat flux. In this paper we examine the occurrence and fidelity of tidal mixing fronts in shelf seas generated within a global 3-dimensional simulation of the HYbrid Coordinate Ocean Model (HYCOM) that is simultaneously forced by atmospheric fields and the astronomical tidal potential. We perform a first order assessment of shelf sea tides in global HYCOM through comparison of sea surface temperature, sea surface tidal elevations, and tidal currents with observations. HYCOM was tuned to minimize errors in M2 sea surface heights in deep water. Over the global coastal and shelf seas (depths <200 m) the area-weighted root mean square error of the M2 sea surface amplitude in HYCOM represents 35% of the 50 cm root mean squared M2 sea surface amplitude when compared to satellite constrained models TPXO8 and FES2014. HYCOM and the altimeter constrained tidal models TPXO8 and FES2014 exhibit similar skill in reproducing barotropic tidal currents estimated from in-situ current meter observations. Through comparison of a global HYCOM simulation with tidal forcing to a global HYCOM simulation with no tides, and also to previous regional studies of tidal mixing fronts in shelf seas, we demonstrate that HYCOM with embedded tides exhibits quite high skill in reproducing known tidal mixing fronts in shelf seas. Our results indicate that the amount of variability in the location of the tidal mixing fronts in HYCOM, estimated using the Simpson-Hunter parameter, is consistent with previous studies when the differences in the net downward heat flux, on a global scale, are taken into account. We also provide evidence of tidal mixing fronts on the North West Australian Shelf for which we have been unable to find references in the existing scientific literature.
Jeon, C. - H., Buijsman, M. C., Wallcraft, A. J., Shriver, J. F., Arbic, B. K., Richman, J. G., et al. (2019). Improving surface tidal accuracy through two-way nesting in a global ocean model. Ocean Modelling , 137 , 98–113.
Abstract: In global ocean simulations, forward (non-data-assimilative) tide models generally feature large sea-surface-height errors near Hudson Strait in the North Atlantic Ocean with respect to altimetry-constrained tidal solutions. These errors may be associated with tidal resonances that are not well resolved by the complex coastal-shelf bathymetry in low-resolution simulations. An online two-way nesting framework has been implemented to improve global surface tides in the HYbrid Coordinate Ocean Model (HYCOM). In this framework, a high-resolution child domain, covering Hudson Strait, is coupled with a relatively low-resolution parent domain for computational efficiency. Data such as barotropic pressure and velocity are exchanged between the child and parent domains with the external coupler OASIS3-MCT. The developed nesting framework is validated with semi-idealized basin-scale model simulations. The M2 sea-surface heights show very good accuracy in the one-way and two-way nesting simulations in Hudson Strait, where large tidal elevations are observed. In addition, the mass and tidal energy flux are not adversely impacted at the nesting boundaries in the semi-idealized simulations. In a next step, the nesting framework is applied to a realistic global tide simulation. In this simulation, the resolution of the child domain (1/75°) is three times as fine as that of the parent domain (1/25°). The M2 sea-surface-height root-mean-square errors with tide gauge data and the altimetry-constrained global FES2014 and TPXO9-atlas tidal solutions are evaluated for the nesting and no-nesting solutions. The better resolved coastal bathymetry and the finer grid in the child domain improve the local tides in Hudson Strait and Bay, and the back-effect of the coastal tides induces an improvement of the barotropic tides in the open ocean of the Atlantic.
Jeon, C. - H., Buijsman, M. C., Wallcraft, A. J., Shriver, J. F., Arbic, B. K., Richman, J. G., et al. Improving surface tidal accuracy through two-way nesting in a global ocean model. Ocean Modelling , .
Smedstad, O. M., Hurlburt, H. E., Metzger, E. J., Rhodes, R. C., Shriver, J. F., Wallcraft, A. J., et al. (2003). An operational Eddy resolving 1/16° global ocean nowcast/forecast system. Journal of Marine Systems , 40-41 , 341–361.
Tilburg, C. E., Hurlburt, H. E., O'Brien, J. J., & Shriver, J. F. (2002). Remote Topographic Forcing of a Baroclinic Western Boundary Current: An Explanation for the Southland Current and the Pathway of the Subtropical Front East of New Zealand*. J. Phys. Oceanogr. , 32 (11), 3216–3232.
Tilburg, C. E., Hurlburt, H. E., O'Brien, J. J., & Shriver, J. F. (2001). The Dynamics of the East Australian Current System: The Tasman Front, the East Auckland Current, and the East Cape Current. J. Phys. Oceanogr. , 31 (10), 2917–2943.
Ansong, J. K., Arbic, B. K., Simmons, H. L., Alford, M. H., Buijsman, M. C., Timko, P. G., et al. (2018). Geographical Distribution of Diurnal and Semidiurnal Parametric Subharmonic Instability in a Global Ocean Circulation Model. J. Phys. Oceanogr. , 48 (6), 1409–1431.
Abstract: The evidence for, baroclinic energetics of, and geographic distribution of parametric subharmonic instability (PSI) arising from both diurnal and semidiurnal tides in a global ocean general circulation model is investigated using 1/12.5° and 1/25° simulations that are forced by both atmospheric analysis fields and the astronomical tidal potential. The paper examines whether PSI occurs in the model, and whether it accounts for a significant fraction of the tidal baroclinic energy loss. Using energy transfer calculations and bispectral analyses, evidence is found for PSI around the critical latitudes of the tides. The intensity of both diurnal and semidiurnal PSI in the simulations is greatest in the upper ocean, consistent with previous results from idealized simulations, and quickly drops off about 5° from the critical latitudes. The sign of energy transfer depends on location; the transfer is positive (from the tides to subharmonic waves) in some locations and negative in others. The net globally integrated energy transfer is positive in all simulations and is 0.5%�10% of the amount of energy required to close the baroclinic energy budget in the model. The net amount of energy transfer is about an order of magnitude larger in the 1/25° semidiurnal simulation than the 1/12.5° one, implying the dependence of the rate of energy transfer on model resolution.
