An official website of the United States government
Abstract The mechanisms and geographic distribution of global tidal dissipation in barotropic tidal models are examined using a high resolution unstructured mesh finite element model. Mesh resolution varies between 2 and 25 km and is especially focused on inner shelves and steep bathymetric gradients. Tidal response sensitivities to bathymetric changes are examined to put into context response sensitivities to frictional processes. We confirm that the Ronne Ice Shelf dramatically affects Atlantic tides but also find that bathymetry in the Hudson Bay system is a critical control. We follow a sequential frictional parameter optimization process and use TPXO9 data‐assimilated tidal elevations as a reference solution. From simulated velocities and depths, dissipation within the global model is estimated and allows us to pinpoint dissipation at high resolution. Boundary layer dissipation is extremely focused with 1.4% of the ocean accounting for 90% of the total. Internal tide friction is much more distributed with 16.7% of the ocean accounting for 90% of the total. Often highly regional dissipation can impact basin‐scale and even ocean wide tides. Optimized boundary layer friction parameters correlate very well with the physical characteristics of the locality with high friction factors associated with energetic tidal regions, deep ocean island chains, and ice covered areas. Global complex M2tide errors are 1.94 cm in deep waters. Total global boundary layer and internal tide dissipation are estimated, respectively, at 1.83 and 1.49 TW. This continues the trend in the literature toward attributing more dissipation to internal tides.
more »
« less
This content will become publicly available on December 1, 2025
Dissipation Scaled Internal Wave Drag in a Global Heterogeneously Coupled Internal/External Mode Total Water Level Model
Abstract This study showcases a global, heterogeneously coupled total water level system wherein salinity and temperature outputs from a coarser‐resolution (12 km) ocean general circulation model are used to calculate density‐driven terms within a global, higher‐resolution (2.5 km) depth‐averaged total water level model. We demonstrate that the inclusion of baroclinic forcing in the barotropic model requires modification of the internal wave drag term to prevent excess degradation of tidal results compared to the barotropic model. By scaling the internal tide dissipation by an easy to calculate dissipation ratio, the resulting heterogeneously coupled model has complex root mean square errors (RMSE) of 2.27 cm in the deep ocean and 12.16 cm in shallow waters for the tidal constituent. While this represents a 10%–20% deterioration as compared to the barotropic model, the improvements in total water level prediction more than offset this degradation. Global median RMSE compared to observations of total water levels, 30‐day sea levels, and non‐tidal residuals improve by 1.86 (18.5%), 2.55 (42.5%), and 0.36 (5.3%) cm respectively. The drastic improvement in model performance highlights the importance of including density‐driven effects within global hydrodynamic models and will help to improve the results of both hindcasts and forecasts in modeling extreme and nuisance flooding. With only an 11% increase in model run time compared to the fully barotropic total water level model, this approach paves the way for high resolution coastal water level and flood models to be used alongside climate models, improving operational forecasting of total water levels.
more »
« less
- Award ID(s):
- 1855047
- PAR ID:
- 10565931
- Publisher / Repository:
- AGU
- Date Published:
- Journal Name:
- Journal of Advances in Modeling Earth Systems
- Volume:
- 16
- Issue:
- 12
- ISSN:
- 1942-2466
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Regional ocean general circulation models are generally forced at the boundaries by mesoscale ocean dynamics and barotropic tides. In this work we provide evidence that remotely forced internal waves can be a significant source of energy for the dynamics. We compare global and regional model solutions within the California Current System. Both models have similar inputs, forcings, and identical grids and numerics. The global model has a steric height power spectrum consistent with mooring observations at superinertial frequencies, while the regional model spectrum is weaker. The regional model also has less sea surface height variance at high wavenumber than the global model. The vertical velocity variance is significantly larger in the global model, except in the sheltered Southern California Bight. While the regional model has roughly equal high‐pass baroclinic and barotropic kinetic energy levels, the global model high‐pass baroclinic kinetic energy is 28% (0.39 PJ) greater than the barotropic energy. An internal wave energy flux analysis reveals that the regional model domain boundaries act as a sink of 183 MW, while in the global model the analysis domain boundaries act as a source of 539 MW. This 722 MW difference can account for the relative increase of 0.39 PJ high‐pass baroclinic energy in the global model, assuming a baroclinic kinetic energy dissipation time in the domain of approximately 6.3 days. The results here imply that most regional ocean models will need to account for internal wave boundary fluxes in order to reproduce the observed internal wave continuum spectrum.more » « less
-
Abstract Motivated by the importance of mixing arising from dissipating internal waves (IWs), vertical profiles of internal‐wave dissipation from a high‐resolution regional ocean model are compared with finestructure estimates made from observations. A horizontal viscosity scheme restricted to only act on horizontally rotational modes (such as eddies) is introduced and tested. At lower resolutions with horizontal grid spacings of 2 km, the modeled IW dissipation from numerical model agrees reasonably well with observations in some cases when the restricted form of horizontal viscosity is used. This suggests the possibility that if restricted forms of horizontal viscosity are adopted by global models with similar resolutions, they could be used to diagnose and map IW dissipation distributions. At higher resolutions with horizontal grid spacings of ∼250 m, the dissipation from vertical shear and horizontal viscosity act much more strongly resulting in dissipation overestimates; however, the vertical‐shear dissipation itself is found to agree well with observations.more » « less
-
The effects of horizontal resolution and wave drag damping on the semidiurnal M2 tidal energetics are studied for two realistically-forced global HYbrid Coordinate Ocean Model (HYCOM) simulations with 41 layers and horizontal resolutions of 8 km (H12) and 4 km (H25). In both simulations, the surface tidal error is minimized by tuning the strength of the linear wave drag, which is a parameterization of the surface-tide energy conversion to the unresolved baroclinic wave modes. In both simulations the M2 surface tide error with TPXO8-atlas, an altimetry constrained model, is 2.6 cm. Compared to H12, the surface tide energy conversion to the resolved vertical modes is increased by 50% in H25. This coincides with an equivalent reduction in the tuned loss of energy from the surface tide to the wave drag. For the configurations studied here, the horizontal and not the vertical resolution is the factor limiting the number of vertical modes that are resolved in most of the global ocean: modes 1–2 in H12 and modes 1–5 in H25. The wave drag also dampens the resolved internal tides. The 40% reduction in wave-drag strength does not result in a proportional increase in the mode-1 energy density in H25. In the higher-resolution simulations, topographic mode-scattering and wave–wave interactions are better resolved. This allows for an energy flux out of mode 1 to the higher modes, mitigating the need for an internal tide damping term. The HYCOM simulations are validated with analytical conversion models and altimetry-inferred sea-surface height, fluxes, and surface tide dissipation. H25 agrees best with these data sets to within 10%. To facilitate the comparison of stationary tide signals extracted from time series with different durations, we successfully apply a spatially-varying correction factor.more » « less
-
In this study, we investigate the transition of semidiurnal Kelvin waves into Hybrid Kelvin-Edge (HKE) waves and associated generation of internal tides at widening shelves using theory, a realistic global baroclinic ocean model simulation, and quasi-realistic regional barotropic model simulations. Using the global model simulation, we identify several areas where a tidal HKE wave transition co-exists with internal wave generation. Of all areas considered, the Celtic Sea/Bay of Biscay shelf has the widest shelf and the strongest internal tide generation. We find that the global simulation agrees better with the theoretical Kelvin modes on the narrow than with the hybrid edge modes on the wide shelves. To help us understand the effect of complex, realistic bathymetry on the HKE wave transition, we perform quasi-realistic 1/25◦ barotropic simulations of the Celtic Sea/Bay of Biscay shelf areas. In these simulations, we gradually change the realistic bathymetry to a more idealized bathymetry. The idealized simulations show that the complex bathymetry steers the barotropic energy flux and causes standing wave patterns, which mask the HKE wave transition. Based on this analysis, we conclude that the HKE wave transition in the Celtic Sea/Bay of Biscay and other shelf areas in the global ocean is most likely masked by the effects of complex bathymetry and that offshelf baroclinic fluxes cannot be exclusively attributed to the HKE wave transition.more » « less
