We document the configuration and emergent simulation features from the Geophysical Fluid Dynamics Laboratory (GFDL) OM4.0 ocean/sea ice model. OM4 serves as the ocean/sea ice component for the GFDL climate and Earth system models. It is also used for climate science research and is contributing to the Coupled Model Intercomparison Project version 6 Ocean Model Intercomparison Project. The ocean component of OM4 uses version 6 of the Modular Ocean Model and the sea ice component uses version 2 of the Sea Ice Simulator, which have identical horizontal grid layouts (Arakawa C‐grid). We follow the Coordinated Ocean‐sea ice Reference Experiments protocol to assess simulation quality across a broad suite of climate‐relevant features. We present results from two versions differing by horizontal grid spacing and physical parameterizations: OM4p5 has nominal 0.5° spacing and includes mesoscale eddy parameterizations and OM4p25 has nominal 0.25° spacing with no mesoscale eddy parameterization. Modular Ocean Model version 6 makes use of a vertical Lagrangian‐remap algorithm that enables general vertical coordinates. We show that use of a hybrid depth‐isopycnal coordinate reduces the middepth ocean warming drift commonly found in pure z* vertical coordinate ocean models. To test the need for the mesoscale eddy parameterization used in OM4p5, we examine the results from a simulation that removes the eddy parameterization. The water mass structure and model drift are physically degraded relative to OM4p5, thus supporting the key role for a mesoscale closure at this resolution.

Internal gravity waves, the subsurface analogue of the familiar surface gravity waves that break on beaches, are ubiquitous in the ocean. Because of their strong Internal gravity waves are propagating disturbances of the ocean's density stratification. Their physics resembles that of surface gravity waves but with buoyancy rather than gravity providing their restoring force -making them much larger (10's to 100's of meters instead of 1 to 10 meters) and slower (hours instead of seconds). Generated primarily by tidal flow past seafloor topography and winds blowing on the sea surface, and typically having multi-kilometer-scale horizontal wavelengths, their estimated 1 TW of deep-sea dissipation is understood to play a crucial role in the ocean's global redistribution of heat and momentum 12 . A major challenge is to improve understanding of internal wave generation, propagation, steepening and dissipation, so that the role of internal waves can be more accurately incorporated in climate models.The internal waves that originate from the Luzon Strait on the eastern margin of the South China Sea (SCS) are the largest documented in the global oceans ( Figure 1).As the waves propagate west from the Luzon Strait they steepen dramatically ( Figure 1a), producing distinctive solitary wave fronts evident in sun glint and synthetic aperture radar (SAR) images from satellites ( Figure 1b). When they shoal onto the continental slope to the west, the downward displacement of the ocean's layers associated with these solitary waves can exceed 250 m in 5 minutes 8 . On such a scale, these waves pose hazards for underwater navigation and offshore drilling 4 , and supply nutrients from the deep ocean that nourish coral reefs 1 and pilot whale populations that forage in their wakes 13 .Over the past decade a number of field studies have been conducted in the region; this work has been comprehensively reviewed 10,11 . All of these studies, however, focused on the propagation of the internal waves across the SCS and their interactions with the continental shelf of China. Until the present study there had been no substantial in situ data gathered at the generation site of the Luzon Strait, in large part because of the extremely challenging operating conditions. A consequence has been persistent 5 confusion regarding the nature of the generation mechanism 11 ; an underlying cause being the sensitivity of the models employed to the system parameters, such as the chosen transect for a two-dimensional model, the linear internal wave speed or the assumed location of the waves' origin within the Luzon Strait. Furthermore, the lack of in situ data from the Luzon Strait has meant an inability to test numerical predictions of energy budgets 9 and no knowledge of the impact of the Kuroshio on the emergence of internal solitary waves 11 .The goal of IWISE is to obtain the first comprehensive in situ data set from the Luzon Strait, which in combination with high-resolution three-dimensional numerical modeling supports a cradle-to-grave picture ...

Turbulent Boussinesq convection under the influence of rapid rotation (i.e. with comparable characteristic rotation and convection timescales) is studied. The transition to turbulence proceeds through a relatively simple bifurcation sequence, starting with unstable convection rolls at moderate Rayleigh (Ra) and Taylor numbers (Ta) and culminating in a state dominated by coherent plume structures at high Ra and Ta. Like non-rotating turbulent convection, the rapidly rotating state exhibits a simple power-law dependence on Ra for all statistical properties of the flow. When the fluid layer is bounded by no-slip surfaces, the convective heat transport (Nu − 1, where Nu is the Nusselt number) exhibits scaling with Ra2/7 similar to non-rotating laboratory experiments. When the boundaries are stress free, the heat transport obeys ‘classical’ scaling (Ra1/3) for a limited range in Ra, then appears to undergo a transition to a different law at Ra ≈ 4 × 107. Important dynamical differences between rotating and non-rotating convection are observed: aside from the (expected) differences in the boundary layers due to Ekman pumping effects, angular momentum conservation forces all plume structures created at flow-convergent sites of the heated and cooled boundaries to spin-up cyclonically; the resulting plume/cyclones undergo strong vortex-vortex interactions which dramatically alter the mean state of the flow and result in a finite background temperature gradient as Ra → ∞, holding Ra/Ta fixed.