Search

Oceanic tides are seldom represented in Earth system models (ESMs) owing to the need for high horizontal resolution to accurately represent the associated barotropic waves close to coasts. This paper presents results of tides implemented in the Model for Prediction Across Scales–Ocean or MPAS-Ocean, which is the ocean component within the U.S. Department of Energy developed Energy Exascale Earth System Model (E3SM). MPAS-Ocean circumvents the limitation of low resolution using unstructured global meshing. We are at this stage simulating the largest semidiurnal (M2, S2, N2) and diurnal (K1, O1) tidal constituents in a single-layer version of MPAS-O. First, we show that the tidal constituents calculated using MPAS-Ocean closely agree with the results of the global tidal prediction model TPXO8 when suitably tuned topographic wave drag and bottom drag coefficients are employed. Thereafter, we present the sensitivity of global tidal evolution due to the presence of Antarctic ice shelf cavities. The effect of ice shelves on the amplitude and phase of tidal constituents are presented. Lower values of complex errors (with respect to TPXO8 results) for the M2 tidal constituents are observed when the ice shelf is added in the simulations, with particularly strong improvement in the Southern Ocean. Our work points towards future research with varying Antarctic ice shelf geometries and sea ice coupling that might lead to better comparison and prediction of tides and thus better prediction of sea-level rise and also the future climate variability.

The autocovariance of the semidiurnal internal tide (IT) is examined in a 32 d segment of a global run of the HYbrid Coordinate Ocean Model (HYCOM). This numerical simulation, with 41 vertical layers and 1/25∘ horizontal resolution, includes tidal and atmospheric forcing, allowing for the generation and propagation of ITs to take place within a realistic eddying general circulation. The HYCOM data are in turn compared with global observations of the IT around 1000 dbar, from Argo float park-phase data and mooring records. HYCOM is found to be globally biased low in terms of the IT variance and decay of the IT autocovariance over timescales shorter than 32 d. Except in the Southern Ocean, where limitations in the model cause the discrepancy with in situ measurements to grow poleward, the spatial correlation between the Argo and HYCOM tidal variance suggests that the generation of low-mode semidiurnal ITs is globally well captured by the model.

The mixing of ocean waters on continental shelves, which is mainly driven by waves, tides, and currents, plays a key role in the physics, biogeochemistry, and ecology of coastal regions. This study focuses on four months of continuous data recorded along a telecommunication cable offshore Oregon, USA, with Distributed Acoustic Sensing (DAS). We apply a cross-correlation approach to the continuous DAS data to infer the propagation of ocean surface gravity waves in the 3 to 100 s period range and estimate near-surface ocean flows. We observe strong spatio-temporal variations of ocean flows along the cable over four months, with strong impacts from a series of storms in late October 2021. We find that our measurements capture oceanic surface motions as those measured by nearby traditional oceanographic instruments. This study demonstrates that ocean-bottom DAS can be used to infer the dynamic properties of near-shore oceans with an unprecedented spatio-temporal resolution.

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 (

SUMMARY Observations of the ψ1 earth tide yield valuable insights into the earth’s free core nutation, especially if the effects of the ψ1 ocean tide can be removed. The ocean tide is extremely small, with amplitudes rarely more than a few millimetres, and developing an accurate model is challenging. Direct observations are inadequate to support a global model. The alternative—numerical simulation—must account for a multitude of possible effects. The ocean tide is forced by the gravitational tidal potential, by pressure loading from atmospheric tides, by seasonal modulation of the nearby K1 constituent, and possibly by non-linear interactions among several other constituents. Here we construct a model of the ψ1 ocean tide which accounts for (or attempts to bound) each of these effects. The radiational component (from atmospheric pressure loading), although relatively small, is complicated by the presence of non-tidal atmospheric variability in the diurnal band. The ocean’s response is dynamic, but there is also high-wavenumber pressure forcing with a near-isostatic response. A general circulation model, forced by both winds and the tidal potential, suggests that annual variability in K1 leads to pronounced ψ1 amplitudes in some marginal seas, especially in the western Pacific.

Marine business and resources play a major role in the economics and way of life in coastal West African countries. Such countries see great profitability from their marine resources while also facing challenges that come with a bordering sea. Despite this fact, there has been limited research into the optimal way for West African Coastal States to coexist with, and sustainably use their marine resources, a research deficit that is mainly due to a lack of infrastructure for in-situ work, lack of capacity development, and comprehensive datasets to undertake oceanographic research. The Coastal Ocean Environment Summer School in Ghana (COESSING; www.coessing.org ) was developed to help meet some of these challenges. Each summer since 2015, ocean scientists (e.g., biologists, chemists, physicists, hydrologists) from the USA and Europe have collaborated with West African colleagues to lead a week-long intensive summer school in Accra, Ghana, alternating in location between the Regional Maritime University and the University of Ghana. The school receives in excess of 100 participants drawn from universities, government agencies, and the private sector organizations, mainly from Ghana and neighboring Liberia, Nigeria, Togo, and Benin, among others. The format of the school includes morning lectures, afternoon field trips, and hands-on laboratory exercises and one-on-one coaching of students. Important to the COESSING program is the satellite oceanography component which introduces participants to the extensive and often free, remotely sensed oceanographic datasets. Participants develop skills that allow them to access, process, and analyze these datasets in order to better understand regional oceanographic phenomena, such as upwelling, pollution, habitat characterization, sea level rise, and coastal erosion. Following the school, facilitators keep in touch with program participants, helping them acquire and analyze data for their studies, dissertations, and often graduate school applications, etc. In summary, schools such as COESSING are critical not only for science in the region but for the global ocean community as such training develops eager, bright minds while leading to improved regional observing and modeling strategies in severely under-sampled seas. Here, we describe a unique case in which satellite oceanography has led to such outcomes for countries bordering the Gulf of Guinea, West Africa.

Ocean–atmosphere coupling modifies the variability of Earth’s climate over a wide range of time scales. However, attribution of the processes that generate this variability remains an outstanding problem. In this article, air–sea coupling is investigated in an eddy-resolving, medium-complexity, idealized ocean–atmosphere model. The model is run in three configurations: fully coupled, partially coupled (where the effect of the ocean geostrophic velocity on the sea surface temperature field is minimal), and atmosphere-only. A surface boundary layer temperature variance budget analysis computed in the frequency domain is shown to be a powerful tool for studying air–sea interactions, as it differentiates the relative contributions to the variability in the temperature field from each process across a range of time scales (from daily to multidecadal). This method compares terms in the ocean and atmosphere across the different model configurations to infer the underlying mechanisms driving temperature variability. Horizontal advection plays a dominant role in driving temperature variance in both the ocean and the atmosphere, particularly at time scales shorter than annual. At longer time scales, the temperature variance is dominated by strong coupling between atmosphere and ocean. Furthermore, the Ekman transport contribution to the ocean’s horizontal advection is found to underlie the low-frequency behavior in the atmosphere. The ocean geostrophic eddy field is an important driver of ocean variability across all frequencies and is reflected in the atmospheric variability in the western boundary current separation region at longer time scales.

Oceanic tides are seldom represented in Earth System Models (ESMs) owing to the need for high horizontal resolution to accurately represent the associated barotropic waves close to coasts. This paper presents results of tides implemented in the Model for Prediction Across Scales–Ocean or MPAS-Ocean, which is the ocean component within the US Department of Energy developed Energy Exascale Earth System Model (E3SM). MPAS-Ocean circumvents the limitation of low resolution using unstructured global meshing. We are at this stage simulating the largest semidiurnal (M2, S2, N2) and diurnal (K1, O1) tidal constituents in a single layer version of MPAS-O. First, we show that the tidal constituents calculated using MPAS-Ocean closely agree with TPXO8 results when suitably tuned topographic wave drag and bottom drag coefficients are employed. Thereafter, we present the sensitivity of global tidal evolution due to the presence of Antarctic ice shelf cavities. The effect of ice shelves on the amplitude and phase of tidal constituents are presented. Lower values of complex errors (with respect to TPX08 results) for the M2 tidal constituents are observed when ice shelf is added in the simulations, with particularly strong improvement in the Southern Ocean. Our work points towards future research with varying Antarctic ice shelf geometries and sea ice coupling that might lead to better comparison and prediction of tides, and thus better prediction of sea-level rise and also the future climate variability.

The interaction of a barotropic flow with topography generates baroclinic motion that exerts a stress on the barotropic flow. Here, explicit solutions are calculated for the spatial-mean flow (i.e., the barotropic tide) resulting from a spatially uniform but time-varying body force (i.e., astronomical forcing) acting over rough topography. This approach of prescribing the force contrasts with that of previous authors who have prescribed the barotropic flow. It is found that the topographic stress, and thus the impact on the spatial-mean flow, depend on the nature of the baroclinic motion that is generated. Two types of stress are identified: (i) a “wave drag” force associated with propagating wave motion, which extracts energy from the spatial-mean flow, and (ii) a topographic “spring” force associated with standing motion at the seafloor, including bottom-trapped internal tides and propagating low-mode internal tides, which significantly damps the time-mean kinetic energy of the spatial-mean flow but extracts no energy in the time-mean. The topographic spring force is shown to be analogous to the force exerted by a mechanical spring in a forced-dissipative harmonic oscillator. Expressions for the topographic stresses appropriate for implementation as baroclinic drag parameterizations in global models are presented.

We consider the power-law spectra of internal gravity waves in a rotating and stratified ocean. Field measurements have shown considerable variability of spectral slopes compared to the high-wavenumber, high-frequency portion of the Garrett–Munk (GM) spectrum. Theoretical explanations have been developed through wave turbulence theory (WTT), where different power-law solutions of the kinetic equation can be found depending on the mechanisms underlying the nonlinear interactions. Mathematically, these are reflected by the convergence properties of the so-called collision integral (CL) at low- and high-frequency limits. In this work, we study the mechanisms in the formation of the power-law spectra of internal gravity waves, utilizing numerical data from the high-resolution modeling of internal waves (HRMIW) in a region northwest of Hawaii. The model captures the power-law spectra in broad ranges of space and time scales, with scalings ω−2.05±0.2 in frequency and m−2.58±0.4 in vertical wavenumber. The latter clearly deviates from the GM76 spectrum but is closer to a family of induced-diffusion-dominated solutions predicted by WTT. Our analysis of nonlinear interactions is performed directly on these model outputs, which is fundamentally different from previous work assuming a GM76 spectrum. By applying a bicoherence analysis and evaluations of modal energy transfer, we show that the CL is dominated by nonlocal interactions between modes in the power-law range and low-frequency inertial motions. We further identify induced diffusion and the near-resonances at its spectral vicinity as dominating the formation of power-law spectrum.

This data repository includes the model and altimetry-based sea surface height wavenumber spectra that arepresented in the manuscript submitted to JGR-Ocean, titled"On the Spatial variability of the sea surface heightWavenumber Spectra in the Atlantic Ocean"&nbsp

With the increase in computational power, ocean models with kilometer-scale resolution have emerged over the last decade. These models have been used for quantifying the energetic exchanges between spatial scales, informing the design of eddy parametrizations, and preparing observing networks. The increase in resolution, however, has drastically increased the size of model outputs, making it difficult to transfer and analyze the data. It remains, nonetheless, of primary importance to assess more systematically the realism of these models. Here, we showcase a cloud-based analysis framework proposed by the Pangeo project that aims to tackle such distribution and analysis challenges. We analyze the output of eight submesoscale-permitting simulations, all on the cloud, for a crossover region of the upcoming Surface Water and Ocean Topography (SWOT) altimeter mission near the Gulf Stream separation. The cloud-based analysis framework (i) minimizes the cost of duplicating and storing ghost copies of data and (ii) allows for seamless sharing of analysis results amongst collaborators. We describe the framework and provide example analyses (e.g., sea-surface height variability, submesoscale vertical buoyancy fluxes, and comparison to predictions from the mixed-layer instability parametrization). Basin- to global-scale, submesoscale-permitting models are still at their early stage of development; their cost and carbon footprints are also rather large. It would, therefore, benefit the community to document the different model configurations for future best practices. We also argue that an emphasis on data analysis strategies would be crucial for improving the models themselves.

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.5Ga 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 252Ma 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 ∼3Ga , 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.5Ga, as indicated by most lunar formation models. Additional processes, missing from our current efforts, are discussed as topics for future investigation.

Oceanography is by nature a global science, and thus requires a global trained workforce. Yet in many coastal nations, the number of trained professionals working in ocean science fields is lacking. Global Ocean Corps and Conveyor (GOCC), an endorsed capacity development programme of the UN Decade of Ocean Science for Sustainable Development, aims to increase the geographical and cultural diversity of the ocean science workforce through facilitating and building sustained long-term education and research collaborations between scientists around the globe. Based upon our collective experience with schools and workshops held in Ghana, Malaysia, University of Rhode Island Coastal Resources Center, and elsewhere, we are confident that a well-funded Ocean Corps would inspire large numbers of scientists, especially early-career scientists, into its ranks, thus molding many of them into champions for international capacity development for the remainder of their careers, and fostering truly global ocean science collaborations worldwide.

Climate variability is investigated by identifying the energy sources and sinks in an idealized, coupled, ocean–atmosphere model, tuned to mimic the North Atlantic region. The spectral energy budget is calculated in the frequency domain to determine the processes that either deposit energy into or extract energy from each fluid, over time scales from one day up to 100 years. Nonlinear advection of kinetic energy is found to be the dominant source of low-frequency variability in both the ocean and the atmosphere, albeit in differing layers in each fluid. To understand the spatial patterns of the spectral energy budget, spatial maps of certain terms in the spectral energy budget are plotted, averaged over various frequency bands. These maps reveal three dynamically distinct regions: along the western boundary, the western boundary current separation, and the remainder of the domain. The western boundary current separation is found to be a preferred region to energize oceanic variability across a broad range of time scales (from monthly to decadal), while the western boundary itself acts as the dominant sink of energy in the domain at time scales longer than 50 days. This study paves the way for future work, using the same spectral methods, to address the question of forced versus intrinsic variability in a coupled climate system.

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.

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

The biotic and abiotic controls on major shifts in atmospheric oxygen and the persistence of low-oxygen periods over a majority of Earth’s history remain under debate. Explanations of Earth’s stepwise pattern of oxygenation have mostly neglected the effect of changing diel illumination dynamics linked to daylength, which has increased through geological time due to Earth’s rotational deceleration caused by tidal friction. Here we used microsensor measurements and dynamic modelling of interfacial solute fluxes in cyanobacterial mats to investigate the effect of changing daylength on Precambrian benthic ecosystems. Simulated increases in daylength across Earth’s historical range boosted the diel benthic oxygen export, even when the gross photosynthetic production remained constant. This fundamental relationship between net productivity and daylength emerges from the interaction of diffusive mass transfer and diel illumination dynamics, and is amplified by metabolic regulation and microbial behaviour. We found that the resultant daylength-driven surplus organic carbon burial could have shaped the increase in atmospheric oxygen that occurred during the Great and Neoproterozoic Oxidation Events. Our suggested mechanism, which links the coinciding increases in daylength and atmospheric oxygen via enhanced net productivity, reveals a possible contribution of planetary mechanics to the evolution of Earth’s biology and geochemistry. Rotational deceleration has increased daylength on Earth, potentially linking the increased burial of organic carbon by cyanobacterial mats and planetary oxygenation, according to experiments and modelling of Precambrian benthic ecosystems.

Internal waves generated at the seafloor propagate through the interior of the ocean, driving mixing where they break and dissipate. However, existing theories only describe these waves in two limiting cases. In one limit, the presence of an upper boundary permits bottom-generated waves to reflect from the ocean surface back to the seafloor, and all the energy flux is at discrete wavenumbers corresponding to resonant modes. In the other limit, waves are strongly dissipated such that they do not interact with the upper boundary and the energy flux is continuous over wavenumber. Here, a novel linear theory is developed for internal tides and lee waves that spans the parameter space in between these two limits. The linear theory is compared with a set of numerical simulations of internal tide and lee wave generation at realistic abyssal hill topography. The linear theory is able to replicate the spatially-averaged kinetic energy and dissipation of even highly non-linear wave fields in the numerical simulations via an appropriate choice of the linear dissipation operator, which represents turbulent wave breaking processes.

Quite a handful of past studies have reported lack of energy near the tidal bands in high-resolution, regional model simulations’ frequency-wave number spectra when compared to observations. A plausible reason for this discrepancy could be the lack of remotely generated internal tides in regional simulations. In this study, we consider the impact of remote internal tides on the energetics in regional simulations of the California Current System (CCS). The CCS is an eddy-rich mid-latitude region, where energetic NIWs and internal tidal waves coexist. We run high-resolution realistic regional simulations using the Regional Ocean Modelling System (ROMS). The ROMS simulations are boundary-forced with high-frequency offline data from the Hybrid Coordinate Ocean Model (HYCOM). We consider a year-long HYCOM expt_06.1 simulation with 8-km horizontal grid resolution and 41 depth layers. The HYCOM simulation is realistically forced with tides and atmospheric forcing.Time-mean and depth-integrated internal tidal flux computed for the parent HYCOM domain shows radiation of remotely generated internal tide beams into the ROMS domain. These beams comprise mainly of modes 1 & 2. To ensure that we provide satisfactory open boundary conditions (OBCs) for our regional simulation, we conduct sensitivity runs using two main types of OBCs (clamped and adaptive OBCs). For the runs with clamped OBCs, we varied the sponge layer viscosities at the open boundaries from 100 to 800 m2/s. Both nudging parameters and sponge layer viscosities are varied for simulations with the adaptive OBCs. Although, we observe remotely generated internal tides in all our simulations, we find that the amount of internal tidal energy that is transmitted through the open boundaries and the internal tidal energetics in the interior of the domain depend on the nudging time scales, sponge layer width and/or viscosity values.In the future, we plan to nest down to increasing high-resolution horizontal and vertical grids and perform simulations with different boundary forcings e.g. with total internal tides, stationary internal tides, and no internal tides. We will also force the ROMS model with unidirectional internal tides. We will evaluate the impacts of these scenarios on the internal tide energetics in the ROMS domain.

Wind generated near-inertial waves (NIWs) are a major source of energy for deep-ocean mixing by transmitting wind energy from the ocean surface into the interior. Recently, it has been established that the NIW energy transmission to ocean depths is significantly modulated by background mesoscale vorticity. Thus, understanding NIW energetics in the presence of mesoscale eddies on a global scale is crucial.We study the generation, propagation and dissipation of NIWs in global 1/25o Hybrid Coordinate Ocean Model (HYCOM) simulations with realistic tidal forcing. The model has 41 layers with uniform vertical coordinates in the mixed layer and isopycnal coordinates in the ocean interior. The model is forced by 1/3hr wind from the NAVGEM atmospheric model. We analyze one month of model data for May-June 2019. The 3D HYCOM fields are projected on vertical normal modes to compute the wind input, wave kinetic energy (KE), flux divergence and dissipation per mode.We find that the globally integrated wind input in surface near-inertial motions is 0.21 TW for the 30-day period and is consistent with previous studies. The sum of the wind input to the first 5 modes accounts to only 31% of the total wind input while the sum of the NIW kinetic energy in the first 5 modes adds up to 60% of the total NIW KE. The difference in the fraction of the total between the wind input and NIW KE (31% and 60%) suggests that a significant portion of wind-induced near-inertial motions is dissipated close to the surface without being projected onto modes. We also find that NIW horizontal fluxes diverge from areas with cyclonic vorticity and converge in areas with anticyclonic vorticity, i.e., anticyclonic eddies are a sink for NIW energy in the global ocean.The residual NIW KE that does not project onto modes is found to be largely trapped in anticyclonic eddies. In a next step, we will study the fate of this energy, which most likely propagate downward as beam-like features with large wave numbers. We will compute the near-inertial wave energy balance for fixed subsurface layers and consider the energy exchange between these layers to understand the vertical structure of NIW energy dissipation. We find that the downward NIW radiation to the ocean interior at 500 m depth is 19% of the surface near-inertial wind input for the 30-day period.