Search

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. Plain Language Summary: Oceanic mixing impacts circulation, stratification (layering by density), and the uptake and transport of heat and nutrients. Over most of the ocean, mixing is caused by the breaking (turnover) of internal waves lying on the interfaces of density layers. Most ocean models do not contain a resolved internal wavefield, and therefore must parameterize internal wave (IW) mixing based upon external information. Recently developed high‐resolution ocean models with credible representations of internal waves may make it possible to map and understand global IW mixing without use of external information. Here we compare vertical profiles (profiles in depth) of IW dissipation in a regional model, which can be used to understand sensitivities to numerical schemes and grid spacings. With grid spacings that are attainable in global models, modeled dissipation profiles lie somewhat close to observed profiles, as long as certain choices are made within the numerical schemes. One numerical dissipation scheme is designed to realistically remove energy from eddy fields, which are the non‐wavelike motions in the ocean, and we have adapted this scheme to act less strongly on internal waves. Using this modified scheme, we find that high‐resolution global models may already be able to map IW dissipation. Key Points: Vertical profiles of internal wave (IW) dissipation in high‐resolution regional ocean simulations are compared with observed profilesProfiles in runs with a restricted form of horizontal viscosity and resolutions attainable in global models are close to observationsResults suggest high‐resolution global models can be used to map IW dissipation after numerical sensitivities are tested in regional models [ABSTRACT FROM AUTHOR]

Through interactions with the continental margins, incident low‐mode internal tides (ITs) can be reflected, scattered to high modes, transmitted onto the shelf and dissipated. We investigate the fate of remotely generated mode‐1 ITs in the U.S. West Coast (USWC) continental margin using two 4‐km horizontal resolution regional simulations. These 1‐year long simulations have realistic stratification, and atmospheric, tidal, and sub‐tidal forcings. In addition, one of these simulations has remote internal wave (IW) forcing at the open boundaries while the other does not. To compute the IT reflectivity of the USWC margin, we separate the IT energy fluxes into onshore and offshore propagating components using a Discrete Fourier Transform in space and time. Overall, ∼20% of the remote mode‐1 semidiurnal IT energy fluxes reflect off the USWC margin, 40% is scattered to modes 2–5, and 7% is transmitted onto the shelf while the remaining is dissipated on the continental slope. Furthermore, our results reveal that differences in stratification, slope criticality, topographic roughness and angle of incidence cause these fractions to vary spatially and temporally along the USWC margin. However, there is no clear seasonal variability in these estimates. Remote IWs enhance the advection and diffusion of heat in the continental margin, resulting in cooling at the surface and warming at depth, and a reduction in the thermocline stratification. These results suggest that low‐mode ITs can cause water mass transformation in continental margins that are far away from their generation sites. Plain Language Summary: Ocean dynamics that exist over a wide range of spatial scales, such as the internal tides (ITs), are best resolved in numerical simulations with fine grid spacings. To accomplish this, model solutions are downscaled by nesting domains with increasingly higher resolution. However, the long propagation distances of low‐mode ITs away from their generation sites necessitate their inclusion at the open boundaries of regional simulations. Here, we show that these remotely generated ITs, for example, from Hawaii, modulate turbulence in the U.S. West Coast (USWC) continental margin. Approximately 80% of the remote mode‐1 semidiurnal IT energy is available for mixing in the USWC margin. In the absence of an accurate IT forcing at the open boundaries of coastal regional models, IT mixing may be underestimated. Key Points: The bulk of the internal tide (IT) energy approaching the U.S. West Coast comes from far‐field sources, for example, HawaiiOn average, 20% of the remote mode‐1 IT energy is reflected while the remaining energy is scattered to higher modes and dissipatedThe remote ITs advect and diffuse heat in the continental margin, causing surface cooling and subsurface warming [ABSTRACT FROM AUTHOR]

It is generally understood that the origin of ocean diapycnal diffusivity is primarily associated with the stratified turbulence produced by breaking internal (gravity) waves (IW). However, it requires significant effort to verify diffusivity values in ocean general circulation models in any particular geographical region of the ocean due to the scarcity of microstructure measurements. Recent analyses of downscaled IW fields from an internal‐wave‐admitting global ocean simulation into higher‐resolution regional configurations northwest of Hawaii have demonstrated a much‐improved fit of the simulated IW spectra to the in‐situ profiler measurements such as the Garrett‐Munk (GM) spectrum. Here, we employ this dynamically downscaled ocean simulation to directly analyze the nature of the IW‐breaking and the wave‐turbulence cascade in this region. We employ a modified version of the Kappa Profile Parameterization (KPP) to infer what the horizontally averaged vertical profile of diapycnal diffusivity should be, and compare this to the background profile that would be employed in the ocean component of a low‐resolution coupled climate model such as the Community Earth System Model (CESM) of the US National Center for Atmospheric Research (NCAR). In pursuing this goal, we also demonstrate that the wavefield in the high‐resolution regional domain is dominated by a well‐resolved spectrum of low‐mode IWs that are predictable by solving an appropriate eigenvalue problem for stratified flow. We finally suggest a new tentative approach to improve the KPP parameterization. Plain Language Summary: A much‐improved spectrum of the simulated internal wave (IW) field has recently been obtained by downscaling a global ocean model into a higher‐resolution regional configuration. The global simulation is based on the Massachusetts Institute of Technology general circulation model (MITgcm) forced by both astronomical tidal potential and surface atmospheric processes. By employing a mathematical framework to predict the structure of IWs, we first demonstrate that the interior wavefield of the high‐resolution regional domain is well dominated by a series of low‐order IW modes. Then, we address the issue as to whether the component of the K‐Profile Parameterization (KPP) associated with IW shear might be able to explain the physical origins of the background depth dependence of diapycnal diffusivity that would normally be employed in the ocean component of a modern coupled climate model. Finally, we suggest a tentative approach to further improve KPP. Key Points: A representation of internal wave modes that includes dissipation is derived that explains spectra in a tidally forced numerical simulationEliminating the background component of KPP leads to a representation of the physics of internal wave breakingRealistic diapycnal diffusivity profiles can be obtained by minor adjustments to the shear component of KPP [ABSTRACT FROM AUTHOR]

This paper examines spectra of horizontal kinetic energy (HKE) in the surface and sub‐surface ocean, with an emphasis on internal gravity wave (IGW) motions, in global high‐resolution ocean simulations. Horizontal wavenumber‐frequency spectra of surface HKE are computed over seven oceanic regions from two global simulations of the HYbrid Coordinate Ocean Model (HYCOM) and three global simulations of the Massachusetts Institute of Technology general circulation model (MITgcm). In regions with high IGW activity, high surface HKE variance in the horizontal wavenumber‐frequency spectra is aligned along IGW linear dispersion curves. For both HYCOM and MITgcm, and in almost all regions, finer horizontal resolution yields more energetic supertidal IGW continuum spectra. The ratio of high‐horizontal‐wavenumber variance in semi‐diurnal and supertidal motions relative to lower‐frequency motions, a quantity of great interest for swath altimetry, depends on the model employed and the horizontal resolution within the model, implying that quantitative predictions of the partition between low‐ and high‐frequency motions taken from particular simulations should be treated with care. The frequency‐vertical wavenumber spectra, frequency spectra, and vertical wavenumber spectra from the models are compared to spectra computed from McLane profilers at nine locations. In general, MITgcm spectra match the McLane profiler spectra more closely at high frequencies (|ω| > 4.5 cpd). In both models, vertical wavenumber spectra roll off more steeply than observations at high vertical wavenumbers (m > 10−2 cpm). The vertical wavenumber spectra in such models is an important target for improvement, due to turbulence production and dissipation that takes place at high vertical wavenumbers. Plain Language Summary: Recently, a small but growing number of global ocean models have begun to employ simultaneous tidal and atmospheric forcing. At the same time, increasing supercomputer power has allowed for simulations of oceanic motions with increasing accuracy, increasing feature (spatial) resolution, and more frequent time slices. Global ocean models with fine grid spacing, and simultaneous tidal and atmospheric forcing, host a vigorous spectrum of high‐frequency waves that control mixing over most of the ocean water column, and are important for many operational oceanography challenges. As an example of the latter, high‐resolution global internal wave models have been used to study the relative partition of high‐frequency versus low‐frequency motions at the small horizontal scales that will be measured by the new Surface Water Ocean Topography mission. The partition described above depends on the model employed and the grid spacing employed within that model, meaning that conclusions about the partition are dependent on the model used to estimate it. Comparisons between the models and vertically profiling instruments indicate that resolving fine scale motions in the vertical direction, where ocean mixing takes place, is not yet handled well by the models. Modeling of fine‐vertical scale motions is therefore an important future research direction. Key Points: Vertical wavenumber spectra of internal gravity wave kinetic energy in two high‐resolution global models are compared to observed spectraModels under‐estimate motions at high vertical wavenumbers (small vertical scales), flagging this as a target for model improvementThe ratio of high‐ versus low‐frequency surface kinetic energy at small horizontal scales is dependent on the model and grid spacings employed [ABSTRACT FROM AUTHOR]

In recent years, efforts have been made to include tides in both operational ocean models as well as climate and earth system models. The accuracy of the barotropic tides is often limited by the model topography, which is in turn limited by model horizontal resolution. In this work, we explore the reduction of barotropic tidal errors in an ocean general circulation model (Modular Ocean Model version 6; MOM6) using sub‐grid scale topography representation. We follow the methodology from Adcroft (2013, https://doi.org/10.1016/j.ocemod.2013.03.002), which utilizes statistics from finer resolution topographic data sets to represent sub‐grid scale features with a light computational cost in a structured finite volume formulation. The geometric effect from sub‐grid scale topography can be introduced to the model with only a few parameters at each grid cell. The porous barriers, which are implemented at the walls of the grid cells, are used to modify transport between grid cells. Our results show that the globally averaged tidal error in lower‐resolution simulations is significantly reduced with the use of porous barriers. We argue this method is a potentially useful tool to improve simulations of tides (and other flows) in low‐resolution simulations. Plain Language Summary: In recent years, significant attention has been given to simulating tides in diverse ocean and climate models, because of the roles of tides in impacting sea levels. In order to better predict future sea level changes, it is therefore essential to improve the accuracy of surface tides simulated in numerical ocean models. One known source of errors that could impact the accuracy of modeled tides stems from the accuracy of ocean topography in the models, which relies on how many grid boxes (i.e., model resolution) are used to represent the complicated seafloor elevations. The model's resolution, however, is often constrained by the efficiency of computations. Our objective in this study is to enhance the portrayal of ocean topography at a given resolution to minimize tidal errors. To achieve this, we employ a technique known as "porous barriers," which uses approximations to mimic the effects of ocean topography that would be omitted in the models, without imposing significant computational burdens. Our findings demonstrate that this technique of porous barriers can substantially reduce tidal errors in our ocean model. We also expect a wider application of this technique to numerous other physical processes within the ocean. Key Points: The accuracy of tides in ocean models is affected by the representation of topography, which is constrained by model horizontal resolutionThe geometric effect of the sub‐grid scale topographic features can be approximated in the model with the porous barriersThe porous barriers can significantly reduce globally averaged errors for barotropic tides [ABSTRACT FROM AUTHOR]

Internal tides (ITs) play a critical role in ocean mixing, and have strong signatures in ocean observations. Here, global IT sea surface height (SSH) in nadir altimetry is compared with an ocean forecast model that assimilates de‐tided SSH from nadir altimetry. The forecast model removes IT SSH variance from nadir altimetry at skill levels comparable to those achieved with empirical analysis of nadir altimetry. Accurate removal of IT SSH is needed to fully reveal lower‐frequency mesoscale eddies and currents in altimeter data. Analysis windows of order 30–120 days, made possible by the frequent (hourly) outputs of the forecast model, remove more IT SSH variance than longer windows. Forecast models offer a promising new approach for global internal tide mapping and altimetry correction. Because they provide information on the full water column, forecast models can also help to improve understanding of the underlying dynamics of ITs. Plain Language Summary: Tidal flow over topographic features on the seafloor generates vertical displacements along the interfaces of ocean layers that have different densities. These vertical displacements at tidal frequencies are known as internal tides. Internal tide displacements are largest well below the sea surface, but also display a sea surface height (SSH) signature that is large enough to be measured by satellite altimeters. Removing internal tide signals from satellite altimeter SSH allows for a more accurate accounting of non‐tidal features, including slowly evolving ocean currents and eddies, that are also measured by altimeters. Here, we show that supercomputer ocean forecast simulations of the global internal tide field are able to remove internal tide SSH from satellite altimeter measurements with a skill level that is comparable to the skill of internal tide SSH removal based upon analysis of the satellite altimeter data itself. Thus, forecast models offer a complementary method for this important task. In addition, forecast models provide information on the entire ocean water column, not just the sea surface. Finally, the hourly outputs of forecast models allow for a greater variety of tidal analysis record lengths than can be achieved with altimeter outputs, which report sea surface height fields much less frequently. Key Points: Global ocean forecast models can accurately simulate both long‐term (phase‐locked) internal tides and their short‐term modulationsOcean forecast models offer a useful complement to empirical models for mapping internal tides and correcting altimetry for internal tidesIn regions of strong internal tides, optimal variance reduction in nadir altimetry is attained through short‐term tidal analyses (∼60 days) [ABSTRACT FROM AUTHOR]

The decay of the low‐mode internal tide due to the superharmonic energy cascade is investigated in a realistically forced global Hybrid Coordinate Ocean Model simulation with 1/25° (4 km) horizontal grid spacing. Time‐mean and depth‐integrated supertidal kinetic energy is found to be largest near low‐latitude internal tide generation sites, such as the Bay of Bengal, Amazon Shelf, and Mascarene Ridge. The supertidal kinetic energy can make up to 50% of the total internal tide kinetic energy several hundred kilometers from the generation sites. As opposed to the tidal flux divergence, the supertidal flux divergence does not correlate with the barotropic to baroclinic energy conversion. Instead, the time‐mean and depth‐integrated supertidal flux divergence correlates with the nonlinear kinetic energy transfers from (sub)tidal to supertidal frequency bands as estimated with a novel coarse‐graining approach. The regular spaced banding patterns of the surface‐intensified nonlinear energy transfers are attributed to semidiurnal mode 1 and mode 2 internal waves that interfere constructively at the surface. This causes patches where both surface tidal kinetic energy and nonlinear energy transfers are elevated. The simulated internal tide off the Amazon Shelf steepens significantly near these patches, generating solitary‐like waves in good agreement with Synthetic Aperture Radar imagery. Globally, we find that regions of high supertidal energy flux also show a high correlation with observed instances of internal solitary waves. Plain Language Summary: Ocean tides generate internal waves at steep ridges and continental shelves in a stratified ocean. In some regions with strong tidal forcing, surface intensified stratification, and weak background rotation, these waves are observed to steepen into solitary‐like waves. In this process, energy is transferred from the tidal band to higher (supertidal) frequencies and wavenumbers. This study investigates "internal tide steepening" in a realistically forced global Hybrid Coordinate Ocean Model (HYCOM) simulation with a 4‐km horizontal grid spacing. In our simulation, supertidal internal wave energy accounts for about 5% of the total energy available to internal tides equatorward of ±25°. However in some regions, such as the Amazon Shelf and Bay of Bengal, this ratio can be as much as 50%. Along the internal tide beams, the supertidal energy transfers are enhanced in regular spaced "patches" where longer and faster mode 1 semidiurnal internal tides overtake the shorter and slower mode 2 internal tides. In these patches the surface velocities of these waves constructively superpose, enhancing the surface kinetic energy and nonlinear energy transfers. At the Amazon Shelf, sharp‐crested solitary‐like waves that emerge in these patches in the HYCOM simulation are in good agreement with satellite observations. Although nonlinear internal tides are simulated in our global ocean model simulation, higher resolutions are needed to better resolve the energy cascade to smaller scales. Key Points: Supertidal energy (>2.67 cycles per day [cpd]) is elevated at low latitudes, making up to 50% of total tidal energy (>0.9 cpd) in some areasSupertidal flux divergence and surface tidal energy reveal banding patterns due to the interaction between mode 1 and 2 internal tidesSupertidal flux divergence is due to energy transfers from internal tides to higher‐harmonic frequencies as computed with coarse graining [ABSTRACT FROM AUTHOR]

We examine ocean tides in the barotropic version of the Model for Prediction Across Scales (MPAS‐Ocean), the ocean component of the Department of Energy Earth system model. We focus on four factors that affect tidal accuracy: self‐attraction and loading (SAL), model resolution, details of the underlying bathymetry, and parameterized topographic wave drag. The SAL term accounts for the tidal loading of Earth's crust and the self‐gravitation of the ocean and the load‐deformed Earth. A common method for calculating SAL is to decompose mass anomalies into their spherical harmonic constituents. Here, we compare a scalar SAL approximation versus an inline SAL using a fast spherical harmonic transform package. Wave drag accounts for energy lost by breaking internal tides that are produced by barotropic tidal flow over topographic features. We compare a series of successively finer quasi‐uniform resolution meshes (62.9, 31.5, 15.7, and 7.87 km) to a variable resolution (45 to 5 km) configuration. We ran MPAS‐Ocean in a single‐layer barotropic mode forced by five tidal constituents. The 45 to 5 km variable resolution mesh obtained the best total root‐mean‐square error (5.4 cm) for the deep ocean (> $ > $1,000 m) M2 ${\mathrm{M}}_{2}$ tide compared to TPXO8 and ran twice as fast as the quasi‐uniform 8 km mesh, which had an error of 5.8 cm. This error is comparable to those found in other forward (non‐assimilative) ocean tide models. In future work, we plan to use MPAS‐Ocean to study tidal interactions with other Earth system components, and the tidal response to climate change. Plain Language Summary: Over the next century, climate change impacts on coastal regions will include floods, droughts, erosion, and severe weather events. The Department of Energy (DoE) is funding the Integrated Coastal Modeling Project to understand these potential risks better. In this paper, we implement tides in the DoE ocean model. Tides themselves respond to climate change, altering coastal flooding risk assessments. We explore the sensitivity of tides to model resolution (the spacing of model gridpoints), ocean‐floor topography, and the so‐called "self‐attraction and loading" (SAL) effect. Self‐attraction and loading occurs as the mass of water in a location fluctuates, causing a deformation of the Earth's crust and changes in the gravitational potential, which must be accounted for when modeling tides. We present a computationally efficient method of calculating the SAL effects and show that it is more accurate than other commonly used approximations. In future work we will examine interactions of tides with other components of the climate system, including sea ice, floating ice shelves, rivers, and current systems. Key Points: We calculate the full self‐attraction and loading (SAL) term inline, in a barotropic configuration of Model for Prediction Across Scales (MPAS‐Ocean)Inclusion of the inline SAL and higher resolution meshes yield improved tidal accuracy in stand‐alone barotropic MPAS‐Ocean configurationsA 45 to 5‐km variable‐resolution mesh provides reduced tidal errors with better computational performance than a quasi‐uniform 8‐km mesh [ABSTRACT FROM AUTHOR]

The wavenumber spectral slope of sea surface height (SSH) computed within the mesoscale range (70–250 km) from satellite altimetry exhibits a large spatial variability which, until now, has not been reproduced in numerical ocean models. This study documents the impacts of including internal tides, high‐resolution bathymetry, and high‐frequency atmospheric variability on the SSH wavenumber spectra in the Atlantic Ocean, using a series of 1/50° North and equatorial Atlantic simulations with a realistic representation of barotropic/baroclinic tides and mesoscale‐to‐submesoscale variability. The results show that the inclusion of internal tides does increase high frequency SSH variability (with clear peaks near 120 and 70 km) and flattens the spectra slope in the mesoscale range in a good agreement with observations. The surface signature of internal tides, mostly in the equatorial Atlantic but also in subtropical regions in the eastern North Atlantic, is the primary reason behind the observed large spatial variability of the spectral slope in the Atlantic. Internal tides are stronger in the tropical regions when compared to higher latitudes because of the stronger barotropic tides and stronger stratification in the upper layer of the water column. High‐resolution bathymetry does play an important role in the internal tide generation on a local scale, but its impact on the modeled large‐scale SSH variability and SSH wavenumber spectra is quite small. High‐frequency wind variability plays only a minor role on the generation of the simulated high‐frequency SSH variability. Impact of altimetry sampling on the computation of the SSH wavenumber spectra is also discussed. Plain Language Summary: The sea surface height (SSH) wavenumber spectrum is a quantification of the magnitude of SSH variability on different spatial scales. The slope of SSH wavenumber spectra is essential to understanding how kinetic energy is transferred across different scales in the turbulent ocean. The most striking feature of the SSH spectral slope as derived from satellite observations is that it has a wide range of values over different parts of the Atlantic Ocean: from an order of −5 in the western boundary current, to an order of −4 to −3 in the interior of subtropical North Atlantic, and to an order of −1 to −0.5 in the tropical Atlantic. In this study, we reproduced the observed spatial variability of the SSH spectral slope in numerical ocean models, and further documented that internal tide is the main driver behind the scenes. The surface manifestation of internal tides flattens the spectral slope of the SSH. This flattening is most significant in the tropical Atlantic, but also occurs in the subtropical eastern North Atlantic. The impacts from other factors such as high‐resolution bathymetry and high‐frequency atmospheric forcing are relatively small. Key Points: The slope of the sea surface height (SSH) power spectra as derived from satellite observations exhibits a large spatial variability, which is now reproduced in our numerical simulationsThe surface manifestation of internal tides with energy peaks near 120 and 70 km is the main driver for spatial variability in the slope of the SSH power spectraHigh‐resolution bathymetry plays an important role in the internal tide generation locally but has a relatively small impact on the SSH power spectra on a broader scale [ABSTRACT FROM AUTHOR]

The geographical variability, frequency content, and vertical structure of near‐surface oceanic kinetic energy (KE) are important for air‐sea interaction, marine ecosystems, operational oceanography, pollutant tracking, and interpreting remotely sensed velocity measurements. Here, KE in high‐resolution global simulations (HYbrid Coordinate Ocean Model; HYCOM, and Massachusetts Institute of Technology general circulation model; MITgcm), at the sea surface (0 m) and at 15 m, are compared with KE from undrogued and drogued surface drifters, respectively. Global maps and zonal averages are computed for low‐frequency (<0.5 cpd), near‐inertial, diurnal, and semidiurnal bands. Both models exhibit low‐frequency equatorial KE that is low relative to drifter values. HYCOM near‐inertial KE is higher than in MITgcm, and closer to drifter values, probably due to more frequently updated atmospheric forcing. HYCOM semidiurnal KE is lower than in MITgcm, and closer to drifter values, likely due to inclusion of a parameterized topographic internal wave drag. A concurrent tidal harmonic analysis in the diurnal band demonstrates that much of the diurnal flow is nontidal. We compute simple proxies of near‐surface vertical structure—the ratio 0 m KE/(0 m KE + 15 m KE) in model outputs, and the ratio undrogued KE/(undrogued KE + drogued KE) in drifter observations. Over most latitudes and frequency bands, model ratios track the drifter ratios to within error bars. Values of this ratio demonstrate significant vertical structure in all frequency bands except the semidiurnal band. Latitudinal dependence in the ratio is greatest in diurnal and low‐frequency bands. Plain Language Summary: It is important to map and understand ocean surface currents because they affect climate and marine ecosystems. Recent advances in global ocean models include the addition of astronomical tidal forcing alongside atmospheric forcing and the usage of more powerful computers that can resolve finer features. Here, we evaluate ocean surface currents in high‐resolution simulations of two different ocean models through comparison with observations from surface drifting buoys. We examine near‐inertial motions, forced by fast‐changing winds; semidiurnal tides, forced by the astronomical tidal potential; diurnal motions, arising from tidal and other sources; and low‐frequency currents and eddies, forced by atmospheric fields. Global patterns in the models and drifters are broadly consistent. The two models differ in their degree of proximity to drifter measurements in the near‐inertial band, most likely due to different update intervals for atmospheric forcing and in the semidiurnal band, most likely due to different damping schemes. A simple proxy for vertical structure of the currents, measured by differences in drifter flows at the surface versus 15 m depth, is tracked reasonably well by the models. Discrepancies between models and observations motivate future improvements in the models. Key Points: We examine frequency content of ocean kinetic energy (KE) at the sea surface (0 m) and 15 m depth with global drifter data and two modelsNear‐surface near‐inertial and tidal KE in numerical models are sensitive to atmospheric forcing frequency and dampingThe ratio 0 m KE/(0 m KE + 15 m KE) in models lies within error bars of the observational ratios over some latitudes and frequency bands [ABSTRACT FROM AUTHOR]

We present improvements in the modeling of the vertical wavenumber spectrum of the internal gravity wave continuum in high‐resolution regional ocean simulations. We focus on model sensitivities to mixing parameters and comparisons to McLane moored profiler observations in a Pacific region near the Hawaiian Ridge, which features strong semidiurnal tidal beams. In these simulations, the modeled continuum exhibits high sensitivity to the background mixing components of the K‐Profile Parameterization (KPP) vertical mixing scheme. Without the KPP background mixing, stronger vertical gradients in velocity are sustained in the simulations and the modeled kinetic energy and shear spectral slopes are significantly closer to the observations. The improved representation of internal wave dynamics in these simulations makes them suitable for improving ocean mixing estimates and for the interpretation of satellite missions such as the Surface Water and Ocean Topography mission. Plain Language Summary: Internal waves (IWs) exist in the ocean interior due to differences in fluid densities. Breaking IWs cause mixing, which has important effects on ocean temperatures and nutrients. Interactions between internal tides generated by tidal flow over bathymetric features and near‐inertial waves generated by wind yield a spectrum of IWs at many frequencies. Here, we compare the IW spectrum in high‐resolution numerical simulations of a region in the North Pacific with observations from moored instruments. We study the effects of the "background" mixing components of the widely used K‐Profile Parameterization (KPP) vertical mixing scheme on the vertical structure of the IW field. The KPP background parameterizes the mixing action of IWs, which is not resolved in coarser‐resolution global ocean models. In our high‐resolution simulations, the IW field is highly active, and the KPP background components turn out to be mostly redundant in this setting. The modeled IW field lies closer to observations when we turn off the KPP background. Improved IW representation in ocean models can play an important role in the accurate representation of IW‐driven mixing in ocean simulations and interpretation of IW signatures from the upcoming Surface Water and Ocean Topography mission. Key Points: Regional ocean simulations are ideal for examining sensitivity of internal gravity wave (IGW) spectra to model mixing parametersTurning off the background components of K‐Profile Parameterization yields more realistic IGW vertical structure in high‐resolution regional modelsIGW spectra are most correctly estimated in models away from tidal generation sites and lateral boundaries [ABSTRACT FROM AUTHOR]

Tides and Earth‐Moon system evolution are coupled over geological time. Tidal energy dissipation on Earth slows Earth′s rotation rate, increases obliquity, lunar orbit semi‐major axis and eccentricity, and decreases lunar inclination. Tidal and core‐mantle boundary dissipation within the Moon decrease inclination, eccentricity and semi‐major axis. Here we integrate the Earth‐Moon system backwards for 4.5 Ga with orbital dynamics and explicit ocean tide models that are "high‐level" (i.e., not idealized). To account for uncertain plate tectonic histories, we employ Monte Carlo simulations, with tidal energy dissipation rates (normalized relative to astronomical forcing parameters) randomly selected from ocean tide simulations with modern ocean basin geometry and with 55, 116, and 252 Ma reconstructed basin paleogeometries. The normalized dissipation rates depend upon basin geometry and Earth′s rotation rate. Faster Earth rotation generally yields lower normalized dissipation rates. The Monte Carlo results provide a spread of possible early values for the Earth‐Moon system parameters. Of consequence for ocean circulation and climate, absolute (un‐normalized) ocean tidal energy dissipation rates on the early Earth may have exceeded today′s rate due to a closer Moon. Prior to ∼3 Ga, evolution of inclination and eccentricity is dominated by tidal and core‐mantle boundary dissipation within the Moon, which yield high lunar orbit inclinations in the early Earth‐Moon system. A drawback for our results is that the semi‐major axis does not collapse to near‐zero values at 4.5 Ga, as indicated by most lunar formation models. Additional processes, missing from our current efforts, are discussed as topics for future investigation. Plain Language Summary: Tidal dissipation in Earth′s oceans and solid body cause the distance to the Moon and the length of day to increase over time. Tides also change the eccentricity and tilt of the lunar orbit, and Earth′s obliquity (the tilt between the equator plane and the ecliptic plane of our orbit around the Sun). This paper attempts to calculate the evolution of the Earth‐Moon system over the whole of Earth′s history using sophisticated ocean tide and orbit models. Over long time scales, the rate at which tidal energy is being dissipated is affected by the geometrical configuration of the continents, the length of day, and mean sea level, which is affected by plate tectonic forces and the presence or absence of large ice caps. The faster rotating Earth of the past was less efficient at dissipating energy and the present placement of the continents enhances some tides due to resonances. In addition, tidal dissipation in the Moon slows the orbit evolution by absorbing energy from the orbit and there was a time in the distant past when the Moon′s tidal dissipation was large. The evolution of the Earth‐Moon system is complex and uncertain, but it can be addressed with advanced models. Key Points: Long‐term Earth‐Moon system evolution is estimated with backwards‐in‐time integrations using high‐level orbit and ocean tide modelsRapid Earth rotation reduces paleotidal energy dissipation rate relative to paleotidal forcing. Ocean basin geometry is another key factorTidal and core/mantle boundary dissipation within the Moon significantly impact the orbital evolution from about 3–4.5 Ga in the past [ABSTRACT FROM AUTHOR]

The location of mixing due to internal tides is important for both the ocean circulation as well as local biogeochemical processes. Numerous observations and modeling studies have shown that submarine canyons may be regions of enhanced internal tide‐driven mixing, but there has not yet been a systematic study of all submarine canyons resolved in bathymetric datasets. Here, we parameterize the internal tide‐driven dissipation from a suite of simulations and pair this with a global high‐resolution, internal tide‐resolving model and bathymetric dataset to estimate the internal‐tide‐driven dissipation that occurs in all documented submarine canyons. We find that submarine canyons dissipate a significant fraction of the incoming internal tide's energy, which is consistent with observations. When globally integrated, submarine canyons are responsible for dissipating 30.8–75.3 GW, or 3.2%–7.8% of the energy input into the M2‐frequency internal tides. This percentage of the internal tide energy that is dissipated in submarine canyons is comparable to or larger than previous calculations using extrapolations from observations of single canyons. Plain Language Summary: Internal waves, or waves that propagate in density layers below the surface of the ocean, are responsible for transporting a significant amount of energy throughout the ocean. When these waves break, they deposit their energy in local mixing events. Submarine canyons have been identified as a type of topography that leads to significant internal wave‐driven mixing. In this study, we use a global map of canyons, together with a high‐resolution ocean model, to calculate the percentage of internal wave energy that is lost to mixing in submarine canyons. We find that a significant fraction of the incident internal wave energy is lost within each submarine canyon. We then sum the energy loss for all submarine canyons to calculate the amount of internal wave‐driven mixing that occurs within canyons over the global ocean. We find that approximately 5%, a non‐negligible amount, of the energy in the global M2‐frequency internal wavefield is lost in canyons. This percentage of the internal wave energy that is lost in submarine canyons is comparable to or larger than previous calculations using extrapolations from observations of single canyons. Key Points: We present a framework for estimating the range of global internal tide‐driven dissipation in submarine canyonsSubmarine canyons may dissipate approximately 5% of the energy input into the global internal tidesPrior estimates of canyon‐driven dissipation may be underestimates, as regions of max canyon‐induced mixing have been un‐ or under‐observed [ABSTRACT FROM AUTHOR]

The flow of tides over rough bathymetry in the deep ocean generates baroclinic motion including internal waves and bottom‐trapped tides. The stresses generated by this motion feedback on the amplitude and phase of the large‐scale tide. Quantifying the stresses associated with tidal flow over abyssal hills is especially important, as this scale of bathymetry is often unresolved in global baroclinic tide models, and the stresses must therefore be parameterized. Here, we extend the previous theoretical work of the authors to determine the amplitude, phasing, and vertical location of the stresses exerted on a flow driven by a time‐periodic body force when it encounters rough bathymetry. The theory compares favorably with a suite of fully nonlinear numerical simulations. It is shown that all topographic stresses are applied directly above the bathymetry, leading to a two‐layer baroclinic flow, with the near‐bottom spatial‐mean flow (the benthic tide) strongly modified by topographic stresses, and the flow at height unperturbed by the presence of topography. Our results provide a framework to improve baroclinic tide models by (i) providing a simple parameterization for the subinertial stress which is currently not included in any models, (ii) establishing that parameterized stresses should be applied in the diffusive boundary layer directly above the topography, independent of where internal tides may dissipate, and (iii) identifying a minimum resolution of ∼10 km for baroclinic tidal models to adequately capture wave resonance effects that can significantly impact the magnitude of the benthic tide. Plain Language Summary: The ocean tide is a periodic motion of the water column at the scale of ocean basins. The amplitude of the ocean tide is determined by a balance between the gravitational attraction of the moon and the forces exerted when the tide flows over the rough seafloor. Much of the roughness in the deep ocean is due to mountains on the seafloor, 10–50 km in width, known as "abyssal hills." In this study, mathematical descriptions are developed for the strength and distribution of the forces exerted on the ocean tide due to its interaction with these hills. These forces always act below the crest of the hills, leading to a sub‐crest or "benthic tide" that differs from the tide above the crest of the hills. Our results provide a framework for improving the representation of tides in numerical ocean models. Key Points: The interaction of the spatial‐mean tide with rough bathymetry generates baroclinic motion including internal tidesThese motions exert stresses on the spatial‐mean tide near the seafloor which modifies the benthic tidal amplitudeThe benthic tidal amplitude may be decreased or increased depending on the nature of the interaction [ABSTRACT FROM AUTHOR]

Temperature variance and kinetic energy (KE) from two global simulations of the HYbrid Coordinate Ocean Model (HYCOM; 1/12° and 1/25°) and three global simulations of the Massachusetts Institute of Technology general circulation model (MITgcm; 1/12°, 1/24°, and 1/48°), all of which are forced by atmospheric fields and the astronomical tidal potential, are compared with temperature variance and KE from a database of about 2,000 moored historical observations (MHOs). The variances are computed across frequencies ranging from supertidal, dominated by the internal gravity wave continuum, to subtidal, dominated by currents and mesoscale eddies. The most important qualitative difference between HYCOM and MITgcm, and between simulations of different resolutions, is in the supertidal band, where the 1/48° MITgcm lies closest to observations. Across all frequency bands examined, the HYCOM simulations display higher spatial correlation with the MHO than do the MITgcm simulations. The supertidal, semidiurnal, and diurnal velocities in the HYCOM simulations also compare more closely with observations than do the MITgcm simulations in a number of specific continental margin/marginal sea regions. To complement the model‐MHO comparisons, this paper also compares the surface ocean geostrophic eddy KE in HYCOM, MITgcm, and a gridded satellite altimeter product. Consistent with the model‐MHO comparisons, the HYCOM simulations have a higher spatial correlation with the altimeter product than the MITgcm simulations do. On the other hand, the surface ocean geostrophic eddy KE is too large, relative to the altimeter product, in the HYCOM simulations. [ABSTRACT FROM AUTHOR]

Scientists and engineers have observed for some time that tidal amplitudes at many locations are shifting considerably due to nonastronomical factors. Here we review comprehensively these important changes in tidal properties, many of which remain poorly understood. Over long geological time scales, tectonic processes drive variations in basin size, depth, and shape and hence the resonant properties of ocean basins. On shorter geological time scales, changes in oceanic tidal properties are dominated by variations in water depth. A growing number of studies have identified widespread, sometimes regionally coherent, positive, and negative trends in tidal constituents and levels during the 19th, 20th, and early 21st centuries. Determining the causes is challenging because a tide measured at a coastal gauge integrates the effects of local, regional, and oceanic changes. Here, we highlight six main factors that can cause changes in measured tidal statistics on local scales and a further eight possible regional/global driving mechanisms. Since only a few studies have combined observations and models, or modeled at a temporal/spatial resolution capable of resolving both ultralocal and large‐scale global changes, the individual contributions from local and regional mechanisms remain uncertain. Nonetheless, modeling studies project that sea level rise and climate change will continue to alter tides over the next several centuries, with regionally coherent modes of change caused by alterations to coastal morphology and ice sheet extent. Hence, a better understanding of the causes and consequences of tidal variations is needed to help assess the implications for coastal defense, risk assessment, and ecological change.Plain Language Summary: Tides are one of the most persistent and dominant forces that shape our planet. The regular and predictable daily, fortnightly, monthly, annual, interannual, and longer‐term changes in tides, driven by astronomical forces, are well understood. However, scientists and engineers have observed for some time that tides at many locations are shifting considerably due to nonastronomical factors. Here, we carry out a review of these important changes in tides, many of which remain poorly understood. We highlight that over long geological time scales, changes in tides are driven by tectonic processes, which alter the size, depth, and shape of the ocean. Over shorter geological time scales, changes in tides are mainly driven by changes in water depth. In recent decades, a growing number of studies have identified widespread, and sometimes regionally coherent, changes in tides during the last 150 years. However, determining exactly what has caused these more recent changes in tides has proven difficult. We discuss the local and regional/global mechanisms that might be responsible for the observed changes. Modeling studies predict that sea level rise and climate change will continue to alter tides over the next several centuries. Therefore, a better understanding of what is causing changes in tides is needed.Key Points: Tidal properties have changed, and continue to evolve, due to nonastronomical factorsAttributing causation remains challengingRegionally coherent increases/decreases in tidal properties are likely to occur over the next centuries [ABSTRACT FROM AUTHOR]

Results from a range of Earth System and climate models of various resolution run under high‐CO2 emission scenarios challenge the paradigm that seafloor CaCO3 dissolution will grow in extent and intensify as ocean acidification develops over the next century. Under the "business as usual," RCP8.5 scenario, CaCO3 dissolution increases in some areas of the deep ocean, such as the eastern central Pacific Ocean, but is projected to decrease in the Northern Pacific and abyssal Atlantic Ocean by the year 2100. The flux of CaCO3 to the seafloor and bottom‐current speeds, both of which are expected to decrease globally through the 21st century, govern changes in benthic CaCO3 dissolution rates over 53% and 31% of the dissolving seafloor, respectively. Below the calcite compensation depth, a reduced CaCO3 flux to the CaCO3‐free seabed modulates the amount of CaCO3 material dissolved at the sediment‐water interface. Slower bottom‐water circulation leads to thicker diffusive boundary layers above the sediment bed and a consequent stronger transport barrier to CaCO3 dissolution. While all investigated models predict a weakening of bottom current speeds over most of the seafloor by the end of the 21st century, strong discrepancies exist in the magnitude of the predicted speeds. Overall, the poor performance of most models in reproducing modern bottom‐water velocities and CaCO3 rain rates coupled with the existence of large disparities in predicted bottom‐water chemistry across models hampers our ability to robustly estimate the magnitude and temporal evolution of anthropogenic CaCO3 dissolution rates and the associated anthropogenic CO2 neutralization. Plain language summary: Carbon dioxide (CO2), produced and released to the atmosphere by human activities, has been accumulating in the oceans for two centuries and will continue to do so well beyond the end of this century if emissions are not curbed. One direct consequence of CO2 buildup in the ocean is the acidification of seawater. Calcite, a mineral secreted by many organisms living in the surface ocean to produce their shells and skeletons, covers a large part of the seafloor and acts as a natural antacid, neutralizing this excess CO2. Model projections for the 21st century, under a "business as usual" scenario, reveal that seawater will become more corrosive to this mineral, but calcite dissolution at the seafloor will only increase slightly due to reductions in bottom‐current speeds and in the amount of calcite particles delivered to the seafloor over that period. These results indicate that the neutralization of human‐made CO2 by calcite dissolution at the seafloor may take longer than previously anticipated. Key Points: Reduced CaCO3 flux to the seafloor and weaker bottom‐current speeds curtail benthic CaCO3 dissolution over the 21st centuryModeled bottom currents underestimate current meter observations by up to 90%Under RCP8.5, the mean calcite compensation depth may rise by ~800 m by the end of this century [ABSTRACT FROM AUTHOR]

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 ∼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. Key Points: Hydrodynamic models incorporating mesoscale dynamics and tides are beginning to resolve stationary and nonstationary baroclinic tidesThe ratio of nonstationary to total semidiurnal variance computed from altimetry and HyCOM simulations agrees at low and middle latitudesComparisons of analysis methodologies show that total and nonstationary semidiurnal variances are underestimated in altimetry on average [ABSTRACT FROM AUTHOR]

Two global ocean models ranging in horizontal resolution from 1/12° to 1/48° are used to study the space and time scales of sea surface height (SSH) signals associated with internal gravity waves (IGWs). Frequency-horizontal wavenumber SSH spectral densities are computed over seven regions of the world ocean from two simulations of the HYbrid Coordinate Ocean Model (HYCOM) and three simulations of the Massachusetts Institute of Technology general circulation model (MITgcm). High wavenumber, high-frequency SSH variance follows the predicted IGW linear dispersion curves. The realism of high-frequency motions (> [ABSTRACT FROM AUTHOR]

The jets in the equatorial Pacific Ocean of a realistically forced global circulation model with a horizontal resolution of [ABSTRACT FROM AUTHOR]