High-resolution Numerical Model Research Articles

Model-Based Study of Near-Surface Transport in and around Cape Cod Bay, Its Seasonal Variability, and Response to Wind

Abstract Output from a high-resolution numerical model is used to study near-surface transport in and around Cape Cod Bay using a Lagrangian approach. Key questions include the following: What are the dominant transport pathways? How do they vary in time on seasonal-to-interannual scales? What is the role of wind in driving this variability? Application to a possible release of wastewater into Cape Cod Bay from the recently closed Pilgrim Nuclear Power Station is discussed. Analysis reveals a seasonality in Cape Cod Bay transport patterns, with shorter residence times throughout the bay and an increased probability of outflow waters exiting the bay during spring and summer. Wind-induced Ekman currents are identified as a dominant driver of this variability. Significance Statement This study is motivated by a possible release of radioisotope-contaminated wastewater into Cape Cod Bay, a region important to fishing, aquaculture, and tourist industries. The specific aim is to better understand near-surface transport patterns and mechanisms in Cape Cod Bay both in general and within the context of a wastewater release from Pilgrim Nuclear Power Station.

Numerical investigation on the evolution process of hydraulic transport and overcurrent blockage of silted urban rainwater pipelines: An approach based on computational fluid dynamics and discrete element method

During urban flood events, the effect of urban rainwater pipeline siltation on overflow and stagflation intensifies the severity of flood disaster. However, the dynamic coupling mechanism of pipeline sedimentation and water flow is still unclear. To investigate the influence of two-phase flow on the hydraulic transport of siltation particles in rainwater pipelines, the numerical simulation model based on computational fluid dynamics (CFD) and discrete element method (DEM) is constructed. Then, the transient continuity governing equation and conservation equation of momentum are formulated to provide dynamic guidance and boundary constraint for CFD-DEM simulation. On this basis, the optimal drag force model and measurement method of equivalent siltation degree of pipeline are proposed and nested with CFD-DEM, and then, a high resolution numerical simulation model of pipeline sedimentation is formulated. The results show that the siltation degree affects the efficiency of drainage pipeline to a degree of 47%, which is much greater than the degree of influence of 33% for siltation length and 18% for slope. When the siltation degree is 0.1, the thickness of the silted bed surface under the influence of water flow scour is reduced by 33%. It revealed that the influence degree of siltation degree and flow rate was 168% and 20%, respectively, which was much larger than that of siltation length and slope. This study can provide technical support for subsequent pipeline cleaning and maintenance as well as flood prevention and mitigation.

Observation and Numerical Simulation of a Windshear Case at an Airport in the Qinghai-Tibet Plateau

This paper documents a windshear case for an airport in the Qinghai-Tibet Plateau and explores, for the first time, the capability for high-resolution numerical weather simulation of the wind shear features. The windshear appears to be associated with pulses of the wind speed in a low-level easterly jet. The features are basically reproduced quite well with the high-resolution numerical model, though some discrepancies are identified, such as the maximum wind speed of the easterly jet and the magnitude of the eddy dissipation rate as compared with the actual Doppler LIDAR observations. Statistical analysis has been performed between the observation and the simulation results. The sensitivity of the modeling result to the choice of turbulence parameterization scheme has also been studied. The study result shows that it is possible to forecast the windshear feature using a high-resolution numerical weather prediction model for an airport in the complex terrain of the Plateau.

Impact of altimeter-buoy data-pairing methods on the validation of Sentinel-3A coastal significant wave heights

Sea state information is critical for a broad range of human activities (e.g. shipping, marine energy, marine engineering) most of them being concentrated along the coastal zone. Satellite altimeter records of significant wave heights (SWH) represent the largest source of sea state observations available to date. However, the quality of altimeter observations is reduced in the coastal zone due to surface heterogeneity within the radar signal footprint. Major difficulties to assess the performance of coastal altimetry in the coastal zone are the reduced number of valid altimeter records and the increased sea state variability, which have recently fostered the development of new methods to pair and compare nearby altimeter and buoy data. In this study, we use a high-resolution numerical wave model implemented over the European coastal waters in order to characterize the spatial variability of sea states in the proximity of coastal in situ buoys, we explore different model-based data-pairing methods to account for coastal sea state variability and we assess their impact on the validation of Sentinel-3A 20Hz SWH measurements. Three Sentinel-3A processing modes are considered: the pseudo low rate mode processing, the SAR processing and the Low Resolution with Range Migration Correction (LR-RMC) processing. Our results indicate major impacts of data-pairing methods on the Sentinel-3A coastal validation and reveals the contribution of more frequent low SWH conditions, poorly resolved by radar altimeters, in the coastal zone as an additional source of errors in coastal altimetry.

Improved Orographic Gravity Wave Drag Parameterization Accounting for the Nonhydrostatic Effect in the Weather Research and Forecasting Model: Tests for Short-Range Forecast of Northeast China Cold Vortices

Abstract The parameterization of orographic gravity wave drag (OGWD) is essential for accurate numerical weather prediction in regions of complex terrain. Current OGWD schemes assume hydrostatic OGWs but the parameterized OGWs in fine-resolution models with narrow subgrid-scale orography can be significantly affected by the nonhydrostatic effects (NHE). In our recent work, the OGWD scheme in the Model for Prediction Across Scales (MPAS) was revised by accounting for the NHE on the surface wave momentum flux of upward-propagating OGWs. Herein, the revised OGWD scheme is implemented in the Weather Research and Forecasting (WRF) model to evaluate its performance in short-range weather forecast. Two sets of 36-hr WRF simulations are conducted for nine Northeast China cold vortices (NECVs) that occurred in the warm season of 2011 using the original and revised OGWD schemes. Results show that the WRF model tends to underestimate the intensity of the NECVs, producing too high geopotential height. When accounting for the NHE in the OGWD scheme, the NECV intensity biases are significantly reduced. Analyses reveal that the NHE act to weaken the lower-tropospheric OGWD by decreasing the surface wave momentum flux, which strengthens the NECV in the lower troposphere. Consequently, the strengthened low-level cyclonic circulation increases the post-trough cold advection to the southwest of the NECV which in turn enhances the NECV in the mid-upper troposphere with reduced geopotential height. The NHE are found to increase as the model horizontal resolution increases, suggesting greater importance of the NHE in the OGWD parameterization of high-resolution numerical models.

High-resolution numerical modelling reveals tsunami risk hotspots in Xiamen City, China

High-resolution numerical modelling reveals tsunami risk hotspots in Xiamen City, China

Trench-parallel mid-ocean ridge subduction driven by along-strike transmission of slab pull

Abstract Trench-parallel mid-ocean ridge (MOR) subduction is observed and/or predicted on the present Earth and during geological history. Slab break-off normally occurs after such a MOR subduction, leading to an absence of local slab pull. The driving force of such MOR subduction is a puzzling issue. We realize that the MOR is generally dislocated by transform faults, which means that the MOR does not enter the trench simultaneously along strike. In this case, the slab pull does not vanish simultaneously along the entire subduction zone. Consequently, the trench-parallel MOR subduction may be driven by along-strike transmission of neighboring slab pull. We tested this idea using a series of 3-D, high-resolution numerical models. The results indicate that the transform fault (TF) and fracture zone (FZ) should not be too weak in order for the lateral transmission of slab pull. Such rheological strength of a TF/FZ is consistent with observation-based inferences and rheological analyses. In addition, the thermal structure and strength of oceanic plates neighboring the MOR also affect the MOR subduction: faster spreading MOR in a younger plate leads to easier subduction. Based on the model results and geological constraints, we propose a self-driven MOR subduction model, which highlights the role of along-strike transmission of slab pull during diachronic entry of MOR into the trench.

Mesoscale ocean eddies determine dispersal and connectivity of corals at the RMS Titanic wreck site

The sinking of the RMS Titanic on 15 April 1912 remains one of most iconic maritime disasters in history. Today, the wreck site lies in waters 3800 m deep approximately 690 km south southeast of Newfoundland, Atlantic Canada. The wreck and debris field have been colonized by many marine organisms including the octocoral Chrysogorgia agassizii. Because of the rapid deterioration of the Titanic and the vulnerability of natural deep-sea coral populations to environmental changes, it is vital to understand the role the Titanic as well as other such structures could play in connecting ecosystems along the North American slope. Based on Lagrangian experiments with more than one million virtual particles and different scenarios for larval behavior, given the uncertainties around the biology of chrysogorgiids, the dispersal of larvae spawned at the Titanic wreck is studied in a high-resolution numerical ocean model. While the large-scale bathymetry shields the Titanic from a strong mean flow, mesoscale ocean eddies can considerably affect the deep circulation and cause a significant speed up, or also a reversal, of the circulation. As a consequence, the position of upper and mid-ocean eddies in the model largely controls the direction and distance of larval dispersal, with the impact of eddies outweighing the importance of active larval swimming in our experiments. Although dependent on larval buoyancy and longevity, we find that the Titanic could be reached by larvae spawned on the upper slope east of the Grand Banks. Therefore, the Titanic could act as a stepping stone connecting the upper to the deep continental slope off Newfoundland. From the Titanic, larvae then spread into deep Canadian waters and areas beyond national jurisdiction.

Urbisphere-Berlin Campaign: Investigating Multiscale Urban Impacts on the Atmospheric Boundary Layer

Abstract For next-generation weather and climate numerical models to resolve cities, both higher spatial resolution and subgrid parameterizations of urban canopy–atmosphere processes are required. The key is to better understand intraurban variability and urban–rural differences in atmospheric boundary layer (ABL) dynamics. This includes upwind–downwind effects due to cities’ influences on the atmosphere beyond their boundaries. To address these aspects, a network of >25 ground-based remote sensing sites was designed for the Berlin region (Germany), considering city form, function, and typical weather conditions. This allows investigation of how different urban densities and human activities impact ABL dynamics. As part of the interdisciplinary European Research Council Grant urbisphere, the network was operated from autumn 2021 to autumn 2022. Here, we provide an overview of the scientific aims, campaign setup, and results from 2 days, highlighting multiscale urban impacts on the atmosphere in combination with high-resolution numerical modeling at 100-m grid spacing. During a spring day, the analyses show systematic upwind-city-downwind effects in ABL heights, largely driven by urban–rural differences in surface heat fluxes. During a heatwave day, ABL height is remarkably deep, yet spatial differences in ABL heights are less pronounced due to regionally dry soil conditions, resulting in similar observed surface heat fluxes. Our modeling results provide further insights into ABL characteristics not resolved by the observation network, highlighting synergies between both approaches. Our data and findings will support modeling to help deliver services to a wider community from citizens to those managing health, energy, transport, land use, and other city infrastructure and operations.

Modulation of Internal Solitary Waves by One Mesoscale Eddy Pair West of the Luzon Strait

Abstract The characteristics of modulated internal solitary waves (ISWs) under the influence of one mesoscale eddy pair in the Luzon Strait, involving one anticyclonic eddy (AE) and one cyclonic eddy (CE) induced by the Kuroshio intrusion, were investigated using a nested high-resolution numerical model in the northeastern South China Sea (SCS). The presence of mesoscale eddies greatly impacts the nonlinear evolution of type-a and type-b ISWs. The eddy pair contributes to distinct wave properties and energy evolutions. Compared to type-b waves, type-a waves display more pronounced modulatory characteristics with a larger spatial scale. CE currents and horizontal inhomogeneous stratification are crucial in modulating the wave behaviors, which induce extremely large-amplitude depression ISWs. The AE thereafter yields retardation effects on the wave energy evolution. The average depth-integrated available potential and kinetic energy showed relative growth rates of −66.12% and −46.07%, respectively, for type-a waves, and −24.26% and −20.15%, respectively, for type-b waves along the propagation path up to the AE core. The deformed and distorted ISW crest lines propagating further northward exhibit a more dramatic shoaling evolution. The maximum total energies of type-a and type-b waves at the north station are approximately 13.5 and 3.5 times, respectively, greater than those at the south station on the continental shelf of the Dongsha Atoll. This work provides essential insights into modulated ISW dynamics under the mesoscale eddy pair within the northeastern SCS deep basin.

High-resolution numerical modelling of seasonal volume, freshwater, and heat transport along the Indian coast

Abstract. Seasonal reversal of winds and equatorial remote forcing influences the circulation of the Arabian Sea (AS) and Bay of Bengal (BoB) basins in the northern Indian Ocean. In this study, we numerically modelled the physical characteristics of the AS and BoB, using the Massachusetts Institute of Technology general circulation model (MITgcm) at a high spatial resolution of 1/20° forced with climatological initial and boundary conditions. The simulated temperature, salinity, and flow fields were validated with satellite and in situ datasets. We then studied the exchange of coastal waters by evaluating transports computed from the model simulations. The alongshore volume transport on the eastern coast is stronger with high seasonal variability due to the poleward-flowing western boundary current and equatorward-flowing East Indian Coastal Current. West coast transport is influenced by large intraseasonal oscillations. The alongshore freshwater transport is an order less than the alongshore volume transport. Out of the net volume transport, freshwater accounts for a maximum of 6.03 % during the southwestern monsoon season, followed by 4.85 % in the post-monsoon season. We observe an inverse relationship between alongshore freshwater and volume transport on the western coast and a direct relationship on the eastern coast. The contribution of eddy-induced heat and freshwater transport was also examined. The relation between net heat transport and net heat flux illustrates the role of coastal currents and equatorial forcing in dissipating heat within the coastal waters. We observed that meridional heat transport over the AS is stronger than over the BoB. Both basins act as a heat source during the summer monsoon and heat sink during the winter. This high-resolution model setup simulates all the important physical climatological patterns, leading to a better understanding of the state of the northern Indian Ocean.

OceanNet: a principled neural operator-based digital twin for regional oceans

While data-driven approaches demonstrate great potential in atmospheric modeling and weather forecasting, ocean modeling poses distinct challenges due to complex bathymetry, land, vertical structure, and flow non-linearity. This study introduces OceanNet, a principled neural operator-based digital twin for regional sea-suface height emulation. OceanNet uses a Fourier neural operator and predictor-evaluate-corrector integration scheme to mitigate autoregressive error growth and enhance stability over extended time scales. A spectral regularizer counteracts spectral bias at smaller scales. OceanNet is applied to the northwest Atlantic Ocean western boundary current (the Gulf Stream), focusing on the task of seasonal prediction for Loop Current eddies and the Gulf Stream meander. Trained using historical sea surface height (SSH) data, OceanNet demonstrates competitive forecast skill compared to a state-of-the-art dynamical ocean model forecast, reducing computation by 500,000 times. These accomplishments demonstrate initial steps for physics-inspired deep neural operators as cost-effective alternatives to high-resolution numerical ocean models.

A Postlaunch Update on the Effects of Instrumental Measurement Errors on SWOT Estimates of Sea Surface Height, Velocity, and Vorticity

Abstract The Ka-band radar interferometer (KaRIn) on the Surface Water and Ocean Topography (SWOT) satellite that was launched in December 2022 is providing the first two-dimensional altimetric views of sea surface height (SSH). Measurements are made across two parallel swaths of 50-km width separated by a 20-km gap. In the data product that will be used for most oceanographic applications, SSH estimates with a footprint diameter of about 3 km are provided on a 2 km × 2 km grid. Early analyses of in-flight KaRIn data conclude that the instrumental noise for this footprint diameter has a standard deviation less than σ3km = 0.40 cm for conditions of 2-m significant wave height. This is a factor of 2.3 better than the prelaunch expectation based on the science requirement specification. The SSH fields measured by KaRIn allow the first satellite estimates of essentially instantaneous surface current velocity and vorticity computed geostrophically from SSH. The effects of instrumental noise on smoothed estimates of velocity and vorticity based on early postlaunch assessments are quantified here as functions of the half-power filter cutoff wavelength of the smoothing. Signal-to-noise ratios for smoothed estimates of velocity and vorticity are determined from simulated noisy KaRIn data derived from a high-resolution numerical model of the California Current System. The wavelength resolution capabilities for σ3km = 0.40 cm are found to be about 17 and 35 km for velocity and vorticity, respectively, which correspond to feature diameters of about 8.5 and 17.5 km, and are better than the prelaunch expectations by about 45% and 35%.

Numerical simulations of the interactions between counterflowing density currents

The collision and interactions between counterflowing density currents frequently occur in nature and engineering practices along the coastal regions. However, the mechanisms and implications of these interactions remain an active research area. In this study, we developed a high-resolution numerical model to simulate the hydrodynamic characteristics of symmetric and asymmetric colliding saline density currents. The aim was to understand the interaction mechanisms, quantify collision parameters, and provide baseline physical insights into the implications of such interactions. The direct numerical simulations reveal that when equal height counterflowing saline density currents collide, the dynamic characteristics of the resulting flow depend on the density difference between the colliding currents, collaborating experimental results. The interaction process initiated mixing and significant vertical motion in the resulting mixed fluid, which attained a maximum rise height of approximately twice that of the colliding currents, thereby enhancing the vertical motion of hypersaline fluid in the ambient column. While the extent of mixing at the collision boundary increases with increasing asymmetry between the colliding currents, a relationship is observed between the buoyancy ratio of the colliding currents and the maximum rise height of the vertically displaced fluid. Given that desalination operations require effective dispersion of brine into the marine environment to minimize local impacts, understanding the collision dynamics of counterflowing saline currents would facilitate improved engineering design and implementation of environmentally compliant coastal management strategies.

Utilizing high-resolution 3D Voronoi meshing to analyze field data from the Brine Availability Test in Salt (BATS)

Salt is an attractive disposal medium for radioactive waste because intact salt is essentially impermeable and non-porous. However, upon drift or borehole excavation a damaged region develops surrounding the excavation which causes increased permeability and porosity creating potential flow paths for brine. Brine leads to corrosion of waste forms and waste packages and is a possible transport vector for radionuclides, so it is important to better understand the early-time behavior and evolution of brine flow in a salt. As a result, this study is part of Task E of DECOVALEX-2023 which focuses on understanding the evolution of thermal, two-phase hydrological, and mechanical processes in the excavation damaged zone in salt. Field measurements from The Brine Availability Test in Salt (BATS) 1a heater experiment are analyzed by implementing a high-resolution three-dimensional numerical model. This salt heater experiment consists of 28 days of heating and 13 days of cooling in a central borehole within bedded salt at the Waste Isolation Pilot Plant (WIPP). Here, the flow simulator PFLOTRAN is utilized; simulations are run on a Voronoi mesh, with temperature-dependent thermal conductivity, permeability and porosity decay away from excavations. The temperature-dependency of permeability is done to match field measurements. Results from the simulation match temperature measured in the field within + /- 0.1 °C and the total brine inflow over the 41-day experiment. This study illustrates that the accuracy of the temperature evolution within salt is critically important when analyzing and modeling experimental data by simulating three heating scenarios of the BATS 1a experiment showing that temperature has a direct effect on total brine inflow.

Active and passive salt diapirs: a numerical study

Summary Salt diapirs dominate the structure in many sedimentary basins and control the preservation and migration of hydrocarbon. The formation of salt diapirs generally falls into two endmember models: active (up-building) and passive (down-building) diapirism. In the active model, salt diapirs rise from salt buoyancy to pierce through the sedimentary overburden, whereas in the passive model, salt diapirs result from differential loading of sediments during deposition. These endmember models are mostly conceptual or kinematic, the mechanics of active and passive diapirism, and their relative roles and interactions in the formation of salt diapirs, remain uncertain. Here, we use two-dimensional high-resolution numerical models to investigate the primary factors and critical conditions for active and passive diapirism. Our results indicate that it is improper to use driving mechanisms to classify salt diapirs, because the buoyancy-driven active salt diapirism involves differential loading, while the passive diapirism requires salt buoyancy. The rise of salt diapirs is more sensitive to the effective viscosity of the overburden than to the salt viscosity. Stiff overburdens could prevent the rise of salt diapirs, but they could be pierced by salt diapirs if plastic yield of the overburden is allowed. During deposition, the coupled salt-sediment deformation, driven by both salt buoyancy and differential loading of sediments, can lead to various diapiric salt structures and minibasins. Regional tectonic stress generally promotes salt diapirism by enhancing strain weakening of salts and overburdens. We suggest that the classification of active and passive salt diapirism is an oversimplification in most cases. We propose a general model of the formation of salt diapirs that usually begins with dome initiation driven by salt buoyancy, followed by syndepositional down-building controlled by sedimentation and differential loading, and ends with canopy formation when sedimentation stops.

Modelling water temperature dynamics for eelgrass (Zostera marina) areas in the nearshore Scotian Shelf

Water temperature is an important environmental factor for many ecological processes in coastal ecosystems. Here, we study water temperature dynamics at a set of study sites on the Atlantic coast of Nova Scotia where eelgrass beds are found. The central emphasis is to predict temperature on scales relevant to coastal ecosystem processes using a high-resolution nearshore oceanographic model based on the Finite Volume Community Ocean Model (FVCOM). The model predictions were evaluated against observed temperature time series at six sites for three years from 2017-2019; the evaluation indicates that the model was able to replicate the temperature variation on time scales from hours to seasonal. We also used various biologically tailored temperature metrics relevant to eelgrass condition, including mean seasonal values and variability, daily ranges, growing degree day (GDD), and warm events, to validate the model against time series observations to better understand the temperature regime at the study sites. Frequency resolved Willmott skill scores were >0.7, and the temperature metrics were well predicted with the exception of a bias in GDD at some of the shallow sites. The eelgrass sites have a wide range of temperature conditions. Mean water temperature in the summer differed by more than 7°C between the shallowest and the deepest sites, and the rate of heat accumulation was fastest at shallow sites which had ≥ 12 extreme warm events per year. While the amplitude of the temperature variations within the high frequency band (<48 hr) was greater in shallower sites, temperature changes on meteorological time scales (48 hr to 60 days) were coherent at all sites, suggesting the importance of coast-wide processes. The results of this study demonstrated that our high resolution numerical model captured biologically relevant temperature dynamics at different time scales and over a large spatial region, and yet still accurately predicted detailed temperature dynamics at specific nearshore sites. Thus, the model can provide important insights into coastal temperature dynamics that are potentially useful for conservation planning and understanding the implications of future change.

Construction of a Real-Time Forecast Model with Deep Learning Techniques for Coastal Engineering and Processes: Nested in a Basin Scale Suite of Models

This study aims to develop a robust and adaptable real-time forecasting system for coastal and estuarine regions, considering the challenges posed by the limited or unavailable forecast data. A real-time forecast system has been designed to handle three distinct scenarios: forecast data from global models are available, unavailable, or intermittently unavailable. To address the challenge of data unavailability, this research proposes the application of a deep learning model (DLM). This study involves the construction of a high-resolution numerical model nested within global models. The performance of the DLM is evaluated by comparing the simulation results of the nested model generated using DLM-predicted data against those obtained using data from global models. The results demonstrate a high correlation and around 90% accuracy for up to the initial 5 days of the forecast. Even in the absence of hindcast data up to 4 days prior to the forecast beginning time, the real-time forecast with DLM can achieve an accuracy exceeding 70% of the forecast data from global models. However, direct application of the developed DLM to an unknown domain results in significantly lower performance, highlighting the importance of retraining the DLM with local data. Overall, this study presents a comprehensive approach to constructing a real-time forecasting system for coastal and estuarine regions that adeptly handles various data availability scenarios through the integration of machine learning techniques.

Submesoscale Kinetic Energy Induced by Vertical Buoyancy Fluxes During the Tropical Cyclone Haitang

AbstractSubmesoscale process is an important part in the kinetic energy cascade from large‐scale circulation to turbulent dissipation, and a key component of the global heat budget. Its dynamic response to weather event is an important process in forecasting marine bio‐chemical matter transport. So how will submesoscale instabilities response to tropical cyclones (TCs) is worth studying. Based on underwater glider observations and 1‐km high resolution numerical modeling, we investigated two TCs (Roke and Haitang)‐induced submesoscale baroclinic instabilities and their dynamic mechanisms in the Northern South China Sea. The TC Haitang induced significant surface cooling, mixed layer deepening, front sharpening, and enhanced the mixed layer baroclinic and symmetric instabilities. The submesoscale kinetic energy also enhanced sharply after TC Haitang, which was higher correlated with increased mesoscale strain rates. The submesoscale energetics analysis revealed that the enhanced frontal submesoscale kinetic energy after TC Haitang was mainly from potential energy release via baroclinic energy conversion. Four groups of sensitivity numerical experiments revealed that the turbulent heat buoyancy flux and the Ekman buoyancy flux contributed equally to the positive baroclinic energy conversion during the TC Haitang. This study helps us to understand the multiscale oceanic energy transfers and submesoscale air‐sea interaction processes.

Role of tidal mixing on ocean exchange through the Strait of Hormuz

We investigate the influence of tides on the exchange of water between the Arabian Gulf and the Sea of Oman through the Strait of Hormuz using a high-resolution numerical model. Two numerical simulations are contrasted, one with and one without tidal forcing. We find that tides suppress exchange through the Strait, by ∼20% in the annual mean, being largest in the summer (∼30%) and diminishing in the winter (∼13%). Tides enhance the parameterised shear-driven vertical mixing inside the Gulf and Strait, mixing warm, relatively fresh surface waters downward thus reducing the density of bottom waters flowing outwards. This reduces the lateral difference of density between Gulf and Sea of Oman and hence the exchange through the Strait. Maximum reductions occur in summer when both the vertical stratification and mixing is the largest.