A numerical world ocean general circulation model Part I. Basic design and barotropic experiment

A new six-layer world ocean general circulation model based on the primitive system of equations is described in detail and its performance in the case of a homogeneous ocean is described. These test integrations show that the model is capable of reproducing the observed mean barotropic or vertically-integrated transport, as well as the seasonal variability of the major ocean gyres. The surface currents, however, are dominated by the Ekman transport, and such non-linear features as the western boundary currents and the equatorial countercurrents are poorly represented. The abyssal boundary countercurrents are also absent due to the lack of thermohaline forcing. The most conspicuous effect of the bottom topography on a homogeneous ocean is seen in the Southern ocean where the calculated Antarctic circumpolar transport through the Drake passage ( ≈ 10 Sv, with bathymetry included) greatly underestimates the observed transport (≈ 100 Sv).     

Modeled Sverdrup flow in the North Atlantic from 11 different wind stress climatologies

In studies of large-scale ocean dynamics, often quoted values of Sverdrup transport are computed using the Hellerman–Rosenstein wind stress climatology. The Sverdrup solution varies, however, depending on the wind set used. We examine the differences in the large-scale upper ocean response to different surface momentum forcing fields for the North Atlantic Ocean by comparing the different Sverdrup interior/Munk western boundary layer solutions produced by a 1/16° linear numerical ocean model forced by 11 different wind stress climatologies. Significant differences in the results underscore the importance of careful selection of a wind set for Sverdrup transport calculation and for driving nonlinear models. This high-resolution modeling approach to solving the linear wind-driven ocean circulation problem is a convenient way to discern details of the Sverdrup flow and Munk western boundary layers in areas of complicated geometry such as the Caribbean and Bahamas. In addition, the linear solutions from a large number of wind sets provide a well-understood baseline oceanic response to wind stress forcing and thus, (1) insight into the dynamics of observed circulation features, by themselves and in conjunction with nonlinear models, and (2) insight into nonlinear model sensitivity to the choice of wind-forcing product.The wind stress products are evaluated and insight into the linear dynamics of specific ocean features is obtained by examining wind stress curl patterns in relation to the corresponding high-resolution linear solutions in conjunction with observational knowledge of the ocean circulation. In the Sverdrup/Munk solutions, the Gulf Stream pathway consists of two branches. One separates from the coast at the observed separation point, but penetrates due east in an unrealistic manner. The other, which overshoots the separation point at Cape Hatteras and continues to flow northward along the continental boundary, is required to balance the Sverdrup interior transport. A similar depiction of the Gulf Stream is commonly seen in the mean flow of nonlinear, eddy-resolving basin-scale models of the North Atlantic Ocean. An O(1) change from linear dynamics is required for realistic simulation of the Gulf Stream pathway. Nine of the eleven Sverdrup solutions have a C-shaped subtropical gyre, similar to what is seen in dynamic height contours derived from observations. Three mechanisms are identified that can contribute to this pattern in the Sverdrup transport contours. Along 27°N, several wind sets drive realistic total western boundary current transport (within 10% of observed) when a 14 Sv global thermohaline contribution is added (COADS, ECMWF 10 m re-analysis and operational, Hellerman–Rosenstein and National Centers for Environmental Prediction (NCEP) surface stress re-analysis), a few drive transport that is substantially too high (ECMWF 1000 mb re-analysis and operational and Isemer–Hasse) and Fleet Numerical Meteorology and Oceanography Center (FNMOC) surface stresses give linear transport that is slightly weaker than observed. However, higher order dynamics are required to explain the partitioning of this transport between the Florida Straits and just east of the Bahamas (minimal in the linear solutions vs. 5 Sv observed east of the Bahamas). Part of the Azores Current transport is explained by Sverdrup dynamics. So are the basic path of the North Atlantic Current (NAC) and the circulation features within the Intra-Americas Sea (IAS), when a linear rendition of the northward upper ocean return flow of the global thermohaline circulation is added in the form of a Munk western boundary layer.     

A zonally averaged,three-basin ocean circulation model for climate studies

A two-dimensional, three-basin ocean model suitable for long-term climate studies is developed. The model is based on the zonally averaged form of the primitive equations written in spherical coordinates. The east-west density difference which arises upon averaging the momentum equations is taken to be proportional to the meridional density gradient. Lateral exchanges of heat and salt between the basins are explicitly resolved. Moreover, the model includes bottom topography and has representations of the Arctic Ocean and of the Weddell and Ross seas. Under realistic restoring boundary conditions, the model reproduces the global conveyor belt: deep water is formed in the Atlantic between 60 and 70°N at a rate of about 17 Sv (1 Sv=10⁶ m3 s^–1) and in the vicinity of the Antarctic continent, while the Indian and Pacific basins show broad upwelling. Superimposed on this thermohaline circulation are vigorous wind-driven cells in the upper thermocline. The simulated temperature and salinity fields and the computed meridional heat transport compare reasonably well with the observational estimates. When mixed boundary conditions (i.e., a restoring condition on sea-surface temperature and flux condition on sea-surface salinity) are applied, the model exhibits an irregular behavior before reaching a steady state characterized by self-sustained oscillations of 8.5-y period. The conveyor-belt circulation always results at this stage. A series of perturbation experiments illustrates the ability of the model to reproduce different steady-state circulations under mixed boundary conditions. Finally, the model sensitivity to various factors is examined. This sensitivity study reveals that the bottom topography and the presence of a submarine meridional ridge in the zone of the Drake Passage play a crucial role in determining the properties of the model bottom-water masses. The importance of the seasonality of the surface forcing is also stressed.     

Impacts of the South China Sea Throughflow on seasonal and interannual variations of the Indonesian Throughflow

Impacts of the South China Sea Throughflow (SCST) on seasonal and interannual variations of the Indonesian Throughflow are studied by comparing outputs from ocean general circulation model (OGCM) experiments with and without the SCST. The observed subsurface maximum in the southward flow through the Makassar Strait is simulated only when the SCST, which is driven by the large-scale wind, is allowed in the model. The mean volume and heat transport by the Makassar Strait Throughflow are reduced by 1.7 Sv and 0.19 PW, respectively, by the existence of the SCST in the model. The difference is particularly remarkable during boreal winter when the SCST reaches its seasonal maximum. Furthermore, the SCST is strengthened during El Niño, leading to the weakening in the southward volume and heat transport through the Makassar Strait by 0.37 Sv and 0.05 PW, respectively. These findings from the OGCM experiments suggest that the SCST may play an important role in climate variability of the Indo-Pacific Ocean.     

Interocean circulation and heat and freshwater budgets of the South China Sea based on a numerical model

The South China Sea (SCS) interocean circulation and its associated heat and freshwater budgets are examined using the results of a variable-grid global ocean model. The ocean model has a 1/6° resolution in the SCS and its adjacent oceans. The model results from 1982 to 2003 show that the western Pacific waters enter the SCS through the Luzon Strait with an annual mean volume transport of 4.80 Sv, of which 1.71 Sv returns to the western Pacific through the Taiwan Strait and East China Sea and 3.09 Sv flows toward the Indian Ocean. The heat in the western Pacific is transported to the SCS with a rate of 0.373 PW (relative to a reference temperature 3.72 °C), while the total heat transport through the outflow straits is 0.432 PW. The net heat transport out of the SCS is thus 0.059 PW, which is balanced by a mean net downward heat flux of 17 W/m² across the SCS air–sea interface. Therefore, the interocean circulation acts as an “air conditioner”, cooling the SCS and its overlaying atmosphere. The SCS contributes a heat transport of 0.279 PW to the Indian Ocean, of which 0.240 PW is from the Pacific Ocean through the Luzon Strait and 0.039 PW is from the SCS interior gained from the air–sea exchange. The Luzon Strait salt transport is greater than the total salt transport leaving the SCS by 3.97 Gg/s, implying a mean freshwater flux of 0.112 Sv (or 3.54 × 10¹² m³/year) from the land discharge and P − E (precipitation minus evaporation). The total annual land discharge to the SCS is estimated to be 1.60 × 10¹² m³/year, the total annual P − E over the SCS is thus 1.94 × 10¹² m³/year, equivalent to a mean P − E of 0.55 m/year. The SCS freshwater contribution to the Indian Ocean is 0.096 Sv. The pattern of the SCS interocean circulation in winter differs greatly from that in summer. The SCS branch of the Pacific-to-Indian Ocean throughflow exists in winter, but not in summer. In winter this branching flow starts at the Luzon Strait and extends to the Karimata Strait. In summer the interocean circulation is featured by a north-northeastward current starting at the Karimata Strait and extending to the Taiwan and Luzon Straits, and a subsurface inflow from the Luzon Strait that upwells into the surface layer in the SCS interior to supply the outward transports.     

Response of the overturning circulation to high-latitude fresh-water perturbations in the North Atlantic

Studies have suggested that sea-ice cover east and west of Greenland fluctuates out-of phase as a part of the Atlantic decadal climate variability, and greater changes are possible under global warming conditions. In this study, the response of the Atlantic meridional overturning circulation (MOC) to the distribution of surface fresh-water flux is explored using a global isopycnal ocean model. An Arctic ice related fresh-water flux of 0.1 Sv entering the Nordic Seas is shown to reduce the maximum overturning by 1 to 2 Sv (10⁶ m³ s^–1). A further decrease of 3 to 5 Sv in the MOC is observed when the fresh-water flux is shifted from the Fram Strait to the southern Baffin Bay area. Surprisingly, the salinity in much of the upper Nordic Seas actually increases when the Arctic fresh-water source is the strongest there, as a result of enhanced global overturning. It reflects the great influence of Labrador Sea convection on this models MOC. By applying a weaker surface fresh-water transport perturbation (0.02 Sv) on the Baffin Bay area and therefore perturbing the Labrador Sea Water (LSW) formation, we have also investigated the interaction between the overflows across the Greenland–Scotland Ridge and the LSW and find that, with the same surface forcing conditions in the Nordic Seas, volume transport of the overflows weakens when the LSW formation intensifies.     

涡分辨率全球海洋环流模式LICOM模拟的吕宋海峡流场的季节变化

分析了一个1/10°的涡分辨率全球环流模式LICOM(LASG/IAP Climate system Ocean Model)对吕宋海峡附近海洋环流的模拟能力。结果表明,模拟的吕宋海峡附近上层环流及输运具有明显的季节变化特征,除6月是东向净流出外,其余月份均为西向流入,冬季流量最大。年平均流量在-3.76 Sv(1 Sv=106 m3/s),其中上层(600 m以上)流量起主要贡献,为-3.60 Sv,与目前已有的研究结果基本一致。南海通过6个海峡完成与外界的水交换,其中吕宋海峡和巴拉巴克海峡是大洋水进入南海的主要通道,其余海峡均以流出为主,流出量最大的是台湾海峡(1.99 Sv),其次是卡里玛塔海峡(1.03 Sv)。进一步分析表明,由季风引起的埃克曼输送量约占吕宋海峡流量的11%,而由季风引起的吕宋海峡压力梯度形成的西向的地转流对吕宋海峡的输运起支配作用。作为黑潮源头的太平洋北赤道流流量对吕宋海峡输运的季节变化也有一定影响。     

The influence of the Bering Strait on the circulation in a coarse resolution global ocean model

An ocean general circulation model of global domain, full continental geometry and bottom topography, is used to study the influence of the Bering Strait on the general circulation by comparing equilibrium solutions obtained with and without a land-bridge between Siberia and Alaska. The model is integrated with restoring boundary conditions (BC) on temperature and salinity, and later, with mixed BC in which a restoring BC on temperature is maintained but a specified flux condition on salinity is imposed. In both cases, the effect of the Bering Strait is to allow a flow of about 1.25–1.5 Sv from the North Pacific to the Arctic Ocean and, ultimately, back to the North Pacific along the western boundary current regions of the Atlantic and Indian Oceans. When a restoring BC on salinity is used, the overturning associated with North Atlantic Deep Water and Antarctic Intermediate Water formation are increased if the Bering Strait is present in the model geometry. The result of switching to a specified flux BC on salinity is to cause a transition in the THC in which the overturning associated with North Atlantic Deep Water formation increases from about 12 Sv to about 22 Sv. This transition occurs in an essentially smooth fashion with no significant variability and is about 12% smaller in magnitude if the Bering Strait is present in the model geometry. Because the Bering Strait appears to exert some influence on the general circulation and the formation of deep water masses, it is recommended that this Strait be included in the geometry of similar resolution models designed to study the deep ocean and potential changes in climate. Correspondence to: CJC Reason     

The effect of Antarctic Circumpolar Current transport on the frontal variability in Drake Passage

The Antarctic Circumpolar Current (ACC) is composed of three major fronts: the Sub-Antarctic Front (SAF), the Polar Front (PF), the Southern ACC Front (SACCF). The locations of these fronts are variable. The PF can shift away from its historical (mean) location by as much as 100 km. The transport of the ACC in Drake Passage varies from its mean (134 Sv) by as much as 60 Sv. A regional numerical circulation model is used to study frontal variability in Drake Passage as affected by a range of volume transports (from 95 Sv to 155 Sv with an interval of 10 Sv). Large transport shifts the fronts northward while the smaller transport causes a southward shift. The mean shifting distance of the PF from the historical mean location is minimum with 135 Sv transport. The SAF and the SACCF are confined by northern and southern walls, respectively, while the PF is loosely controlled by the topography. Due to impact of the eddies and meanders on the PF at several regions in Drake Passage, the PF may move northward to join the SAF or move southward to combine with the SACCF, especially in central Scotia Sea. The SAF and PF are more stable with higher transport. The SAF behaves as a narrow, strong frontal jet with large transport while displaying meanders with smaller transport. In the model simulations, the Ertel Potential Vorticity (EPV) is strongly correlated with the volume transport stream function. EPV at depths between 1000 and 2500 m is correlated with the transport stream function with a coefficient above 0.9. Near the bottom, the correlation is about 0.6 due to the disruptive influence of bottom topography. Within 750 m of the surface, the correlation is much reduced due to the effect of K-Profile Parameterization (KPP) mixing and eddy mixing.     

Decadal oscillations in a simplified coupled model due to unstable interactions between zonal winds and ocean gyres

A simplified coupled ocean–atmosphere model, consisting of a one-layer bidimensional ocean model and a one-layer unidimensional energy balance atmospheric model [J. Clim. 13 (2000) 232] is used to study the unstable interactions between zonal winds and ocean gyres. In a specific range of parameters, decadal variability is found. Anomalies, quite homogeneous zonally, show small-scale wavelength in latitude: perturbations emerge and grow at the southern limb of the intergyre boundary and propagate southward before decaying. The wind stress anomalies are proportional to the meridional gradient of the atmospheric temperature anomalies: this ratio acts as a positive amplification factor, as confirmed by a parameter sensitivity analysis. Assuming zonally-averaged anomalies harmonic in the meridional direction, a very simple analytical model for the perturbations is derived, based on forced Rossby wave adjustment of the western boundary current and its associated anomalous heat transport: it accounts for the scale selection, the growth and the southward propagation of sea surface temperature anomalies in the subtropical gyre. The latter is not only due to the slow advection by the mean current, but to a prevailing mechanism of self-advecting coupled oceanic and atmospheric waves, out of phase in latitude. Relevance to the observational record is discussed.     

Three-dimensional perspectives of the Florida current: transport, potential vorticity, and related dynamical properties

The Florida Current flows through the Straits of Florida, which starts as a zonal channel and turns to become a meridional channel. The spatial structure of the Florida Current and its transport, potential vorticity, and related dynamical properties are investigated using a three-dimensional, baroclinic, primitive equation model with a mesoscale-admitting (5.6 km) horizontal resolution and 25 vertical (sigma: terrain-following) levels. At 83°W, the Florida Current fills only a portion of the channel; however, due to the interaction with the shoaling bottom topography (from a maximum depth of over 2000 m at 83°W to less than 800 m at 27°N) and the narrowing Straits of Florida (from a maximum width of about 170 km at 83°W to about 110 km at 27°N), the Florida Current fills the entire channel at 27°N, and the potential vorticity distribution is altered. The specified transport of 28.6 Sverdrup (1 Sv = 10⁶ m³ s⁻¹) from the Loop Current at the western boundary and the inflow from the Old Bahama Channel of 1.9 Sv converge into the meridional channel. With an additional inflow of 1.2 Sv from the Northwest Providence Channel, the simulated total transport of 31.8 Sv at 27°N is comparable to the STACS (Subtropical Atlantic Climate Studies) mean transport of 31.7 Sv. Both vertically and laterally integrated subsectional transports are examined at transects 83°W, 82°W, 81°W, 25°N, 26°N, and 27°N. The potential vorticity increases (decreases) on the cyclonic (anticyclonic) side of the Florida Current at 27°N compared to 83°W. The downstream variation of static stability, relative vorticity, and Froude number is also examined. While the vertical shear is strong only on the northern side at 83°W it is comparable on the both western and eastern sides downstream at 27°N, reaching to the bottom of the meridional channel. Large values of the Froude number exist only in the upper 300 m of the zonal channel, but they reach to the bottom of the meridional channel.     

The effect of clouds on aerosol growth in the rural atmosphere

Measurements of accumulation mode aerosol in the atmospheric boundary layer under cloudy and cloud-free conditions, and in the lower free troposphere under cloud-free conditions, were conducted over the rural northwest of England. Normalised size distributions in the cloud-free boundary layer (CFBL) and the cloud-free free troposphere (CFFT) exhibited almost identical spectral similarities with both size distributions possessing a concentration peak mode-radius of ≈0.05 μm or less. By comparison, aerosol distributions observed in cloudy air exhibited a distinctive log-normal distribution with mode-radii occurring at ≈0.1 μm concomitant with a local minimum at ≈0.05 μm. The consistent and noticeable difference in spectral features observed between cloudy and cloud-free conditions suggest that a greater amount of gas-to-particle conversion occurs on cloudy days, presumably through in-cloud aqueous phase oxidation processes, leading to larger sized accumulation mode particles. Apart from the distinct difference between cloudy and cloud-free aerosol spectra on cloudy days, aerosol concentration and mass were observed to be significantly enhanced above that of the ambient background in the vicinity of clouds. Volatility analysis during one case of cloud processing indicated an increase in the relative contribution of aerosol mass volatile at temperatures characteristic of sulphuric acid, along with a smaller fraction of more volatile material (possibly nitric acid and/or organic aerosol). Growth-law analysis of possible growth mechanisms point to aqueous phase oxidation of aerosol precursors in cloud droplets as being the only feasible mechanism capable of producing the observed growth. The effect of cloud processing is to alter the cloud condensation nuclei (CCN) supersaturation spectrum in a manner which increases the availability of CCN at lower cloud supersaturations.     

A diagnostic calculation of the circulation in the upper and middle layers of the Luzon Strait and the northern South China Sea during March 1992

Wind data from NCEP and hydrographic data obtained during 8–27 March 1992 have been used to compute circulation in the Luzon Strait and the northern South China Sea using three-dimensional diagnostic models with a modified inverse method. Numerical results are as follows: the main Kuroshio is located above 800 m levels. It has two intrusive branches of the Kuroshio in the areas above 400 m. One part intrudes anti-cyclonically northwestward, then flows through the area above 200 m southwest of Taiwan and into the Taiwan Strait. The other part intrudes westward and flows cyclonically in the areas north of the cyclonic eddies, then flows southward through the southern boundary of the region. The net westward volume transport (VT) through Section at 120°15′E between Luzon Island and Taiwan Island is about 3.0 Sv, net northward VT through northern boundaries into the Taiwan Strait is about 1.4 Sv and net southward VT through southern boundaries is about 1.6 Sv, which finally flows into the Karimata and Mindoro Straits. In the areas above 400 m east of 117°15′E, the circulation is mainly dominated by the basin-scale cyclonic gyre, which consists of two cyclonic eddies. However, in the areas below 400 m east of 119°00′E, the circulation is mainly dominated by basin-scale anti-cyclonic gyre. The joint effect of baroclinity and relief and interaction between wind stress and relief are important in different area respectively for the pattern of the depth-averaged flow across contours of fH⁻¹.     

A Note on the South China Sea Shallow Interocean Circulation

1. IntroductionThe South China Sea (SCS) has many channelsconnecting with the outer oceans/seas (Fig. 1). Thewidest and deepest channel is the Luzón Strait, whichis the main entrance to the SCS from the WesternPacific Ocean, having a sill depth of about 2500 m.On the north, the Taiwan Strait connects with theEast China Sea, with a sill depth of about 70 m. Inthe vicinity of Mindoro Island, there are a numberof channels connecting the SCS with the Sulu Sea.The main channel is the M…     

The seasonal variation of the upper layers of the South China Sea (SCS) circulation and the Indonesian through flow (ITF): An ocean model study

The upper layer, wind-driven circulation of the South China Sea (SCS), its through-flow (SCSTF) and the Indonesian through flow (ITF) are simulated using a high resolution model, FVCOM (finite volume coastal ocean model) in a regional domain comprising the Maritime Continent. The regional model is embedded in the MIT global ocean general circulation model (ogcm) which provides surface forcing and boundary conditions of all the oceanographic variables at the lateral open boundaries in the Pacific and Indian oceans. A five decade long simulation is available from the MITgcm and we choose to investigate and compare the climatologies of two decades, 1960–1969 and 1990–1999.The seasonal variability of the wind-driven circulation produced by the monsoon system is realistically simulated. In the SCS the dominant driving force is the monsoon wind and the surface circulation reverses accordingly, with a net cyclonic tendency in winter and anticyclonic in summer. The SCS circulation in the 90s is weaker than in the 60s because of the weaker monsoon system in the 90s. In the upper 50 m the interaction between the SCSTF and ITF is very important. The southward ITF can be blocked by the SCSTF at the Makassar Strait during winter. In summer, part of the ITF feeds the SCSTF flowing into the SCS through the Karimata Strait. Differently from the SCS, the ITF is primarily controlled by the sea level difference between the western Pacific and eastern Indian Ocean. The ITF flow, consistently southwestward below the surface layer, is stronger in the 90s.The volume transports for winter, summer and yearly are estimated from the simulation through all the interocean straits. On the annual average, there is a ∼5.6 Sv of western Pacific water entering the SCS through the Luzon Strait and ∼1.4 Sv exiting through the Karimata Strait into the Java Sea. Also, ∼2 Sv of SCS water enters the Sulu Sea through the Mindoro Strait, while ∼2.9 Sv flow southwards through the Sibutu Strait merging into the ITF. The ITF inflow occurs through the Makassar Strait (up to ∼62%) and the Lifamatola Strait (∼38%). The annual average volume transport of the ITF inflow from the simulation is ∼15 Sv in the 60s and ∼16.6 Sv in the 90s, very close to the long term observations. The ITF outflow through the Lombok, Ombai and Timor straits is ∼16.8 Sv in the 60s and 18.9 Sv in the 90s, with the outflow greater by 1.7 Sv and 2.3 Sv respectively. The transport estimates of the simulation at all the straits are in rather good agreement with the observational estimates.We analyze the thermal structure of the domain in the 60s and 90s and assess the simulated temperature patterns against the SODA reanalysis product, with special focus on the shallow region of the SCS. The SODA dataset clearly shows that the yearly averaged temperatures of the 90s are overall warmer than those of the 60s in the surface, intermediate and some of the deep layers and the decadal differences (90s − 60s) indicate that the overall warming of the SCS interior is a local effect. In the simulation the warm trend from the 60s to the 90s in well reproduced in the surface layer. In particular, the simulated temperature profiles at two shallow sites at midway in the SCSTF agree rather well with the SODA profiles. However, the warming trend in the intermediate (deep) layers is not reproduced in the simulation. We find that this deficiency is mostly due to a deficiency in the initial temperature fields provide by the MITgcm.     

A homogeneous model of the Deacon cell in a circumpolar ocean

Dynamics of the three-dimensional structure of the wind-driven Deacon cell in a β-plane channel are discussed in a homogeneous model in the presence of a sufficiently high ridge. The emphasis is on the water mass balance: how the northward surface Ekman drift is returned. It is demonstrated that a sufficiently high ridge can break up the geostrophic constraint and a net geostrophic volume flux across the open latitude band is allowed. It is found that: (1) the Deacon cell is a fundamentally three-dimensional structure, (2) wind forcing can drive an inter-basin water mass exchange in the Southern Ocean, and (3) zonal through-channel transport in the circumpolar ocean varies at different longitudes.     

Steering of upper ocean currents and fronts by the topographically constrained abyssal circulation

A two-layer theory is used to investigate (1) the steering of upper ocean current pathways by topographically constrained abyssal currents that do not impinge on the bottom topography and (2) its application to upper ocean – topographic coupling via flow instabilities where topographically constrained eddy-driven deep mean flows in turn steer the mean pathways of upper ocean currents and associated fronts. In earlier studies the two-layer theory was applied to ocean models with low vertical resolution (2–6 layers). Here we investigate its relevance to complex ocean general circulation models (OGCMs) with high vertical resolution that are designed to simulate a wide range of ocean processes. The theory can be easily applied to models ranging from idealized to complex OGCMs, provided it is valid for the application. It can also be used in understanding some persistent features seen in observed ocean frontal pathways (over deep water) derived from satellite imagery and other data. To facilitate its application, a more thorough explanation of the theory is presented that emphasizes its range of validity. Three regions of the world ocean are used to investigate its application to eddy-resolving ocean models with high vertical resolution, including one where an assumption of the two-layer theory is violated. Results from the OGCMs with high vertical resolution are compared to those from models with low vertical resolution and to observations. In the Kuroshio region upper ocean – topographic coupling via flow instabilities and a modest seamount complex are used to explain the observed northward mean meander east of Japan where the Kuroshio separates from the coast. The Japan/East Sea (JES) is used to demonstrate the impact of upper ocean – topographic coupling in a relatively weak flow regime. East of South Island, New Zealand, the Southland Current is an observed western boundary current that flows in a direction counter to the demands of Sverdrup flow and counter to the direction simulated in nonlinear global flat bottom and reduced gravity models. A model with high vertical resolution (and topography extending through any number of layers) and a model with low vertical resolution (and vertically compressed but otherwise realistic topography confined to the lowest layer) both simulate a Southland Current in the observed direction with dynamics depending on the configuration of the regional seafloor. However, the dynamics of these simulations are very different because the Campbell Plateau and Chatham Rise east and southeast of New Zealand are rare features of the world ocean where the topography intrudes into the stratified water column over a relatively broad area but lies deeper than the nominal 200 m depth of the continental shelf break, violating a limitation of the two-layer theory. Observations confirm the results from the high vertical resolution model. Overall, the model simulations show increasingly widespread upper ocean – topographic coupling via flow instabilities as the horizontal resolution of the ocean models is increased, but fine resolution of mesoscale variability and the associated flow instabilities are required to obtain sufficient coupling. As a result, this type of coupling is critical in distinguishing between eddy-resolving and eddy-permitting ocean models in regions where it occurs.     

Steady-state frictional geostrophic circulation in a one-layer ocean model with thermodynamics

The governing equations are developed for a steady-state frictional geostrophic inhomogeneous 1.5-layer ocean model, with horizontal velocity field that is linearly sheared in the vertical coordinate. We show that in the adiabatic, thermally non-diffusive limit there are an infinite number of solutions for the temperature and depth fields of the subtropical gyre even with the constraint of identical mass within each temperature range. In the non-adiabatic case, a unique subtropical gyre solution exists that can exhibit a temperature front, containing an unbounded meridional gradient, in the northwest corner of the solution domain. The role of mixing of enthalpy in the western boundary layer (WBL) region was investigated by comparing the two extreme cases of no mixing and complete mixing of enthalpy in this region. Also investigated was the dependence of the meridional heat transport on the air–sea heat exchange coefficient, κ. The temperature field was found to be strongly influenced by mixing. However, both qualitatively and quantitatively, the heat transport is similar in the model with and without mixing. The heat transport attains a single local maximum at κ=κ_c, that lies within values that are oceanographically relevant.     

Coupled oceanic and atmospheric mixed layer model

A coupled, one-dimensional atmospheric-oceanic boundary layer model based on a single station assessment has been formulated from independent oceanic and atmospheric bulk boundary layer models. Sensitivity analyses are conducted to determine major differences in the responses of the coupled model compared to those of the separate oceanic and atmospheric models. The general behavior of the coupled model atmosphere is not significantly different from that of the atmospheric model over short term simulations (12–24 h). However, large differences may occur under certain limited conditions when winds are light and the lifting condensation level is close to the height of the inversion. Major differences between the predicted evolution of the ocean boundary layer by the ocean model and coupled model are more common, and the short term predictive ability of the ocean model in coupled form is enhanced.     

Microhydrodynamics of flotation processes in the sea surface layer

The uppermost surface of the ocean forms a peculiarly important ecosystem, the sea surface microlayer (SML). Comprising the top 1–1000 μm of the ocean surface, the SML concentrates many chemical substances, particularly those that are surface active. Important economically as a nursery for fish eggs and larvae, the SML unfortunately is also especially vulnerable to pollution. Contaminants that settle out from the air, have low solubility, or attach to floatable matter tend to accumulate in the SML.Bubbles contribute prominently to the dynamics of air–sea exchanges, playing an important role in geochemical cycling of material in the upper ocean and SML. In addition to the movement of bubbles, the development of a bubble cloud interrelates with the single particle dynamics of all other bubbles and particles. In the early sixties, several in situ oceanographic techniques revealed an “unbelievably immense” number of coastal bubbles of radius 15–300 μm. The spatial and temporal variation of bubble numbers were studied; acoustical oceanographers now use bubbles as tracers to determine ocean processes near the ocean surface. Sea state and rain noises have both been definitively ascribed to the radiation from huge numbers of infant micro bubbles [The Acoustic Bubble. Academic Press, San Diego].Our research programme aims at constructing a hydrodynamic model for particle transport processes occurring at the microscale, in multi-phase flotation suspensions. Current research addresses bubble and floc microhydrodynamics as building blocks for a microscale transport model. This paper reviews sea surface transport processes in the microlayer and the lower atmosphere, and identifies those amenable to microhydrodynamic modelling and simulation. It presents preliminary simulation results including the multi-body hydrodynamic mobility functions for the modelling of “dynamic bubble filters” and floc suspensions. Hydrodynamic interactions versus spatial anisotropy and size of particle clouds are investigated.