Global Hybrid Coordinate Ocean Model Research Articles

Mesoscale eddies in the Gulf of Mexico: A three-dimensional characterization based on global HYCOM

The Gulf of Mexico (GoM) is characterized by strong mesoscale eddy activities that have been studied extensively, yet the comprehensive three-dimensional (3-D) kinematic properties of GoM eddies are still not well documented. In this study, the 3-D mesoscale eddy activities in the upper layer (0–800 m) of the GoM are characterized using 14-year (1997–2010) global Hybrid Coordinate Ocean Model (HYCOM) outputs. Most eddies in the upper layer (both cyclonic and anticyclonic) have radii of ∼30–60 km and lifespans shorter than 30 days. The spatial distributions of GoM eddies do not vary much with depth, while their intensity decreases with depth. The size of cyclonic eddies does not vary much with depth, while the size of anticyclonic eddies decreases slightly with depth. Cyclonic eddies are often found to be generated in the eastern GoM (especially in the Loop Current region), the Bay of Campeche, and on the continental slope of the Campeche Bank, while anticyclonic eddies are often generated on the northeastern and northwestern GoM continental slopes, and in the central GoM (near 24°N) and the Bay of Campeche (92–94°W). In addition, long-lived GoM eddies (e.g., lifespan >150 days) tend to have intermediate eddy intensity (e.g., 0.13–0.32 for cyclonic eddies at the 10 m level). Both cyclonic and anticyclonic eddies are found to play an important role in the horizontal and vertical transport of heat and salt, and eddy-induced anomalies of water temperature and salinity at both surface and subsurface are generally more pronounced in the eastern GoM than in the western GoM.

Nonlinear Internal Tides in a Realistically Forced Global Ocean Simulation

AbstractThe decay of the low‐mode internal tide due to the superharmonic energy cascade is investigated in a realistically forced global Hybrid Coordinate Ocean Model simulation with 1/25° (4 km) horizontal grid spacing. Time‐mean and depth‐integrated supertidal kinetic energy is found to be largest near low‐latitude internal tide generation sites, such as the Bay of Bengal, Amazon Shelf, and Mascarene Ridge. The supertidal kinetic energy can make up to 50% of the total internal tide kinetic energy several hundred kilometers from the generation sites. As opposed to the tidal flux divergence, the supertidal flux divergence does not correlate with the barotropic to baroclinic energy conversion. Instead, the time‐mean and depth‐integrated supertidal flux divergence correlates with the nonlinear kinetic energy transfers from (sub)tidal to supertidal frequency bands as estimated with a novel coarse‐graining approach. The regular spaced banding patterns of the surface‐intensified nonlinear energy transfers are attributed to semidiurnal mode 1 and mode 2 internal waves that interfere constructively at the surface. This causes patches where both surface tidal kinetic energy and nonlinear energy transfers are elevated. The simulated internal tide off the Amazon Shelf steepens significantly near these patches, generating solitary‐like waves in good agreement with Synthetic Aperture Radar imagery. Globally, we find that regions of high supertidal energy flux also show a high correlation with observed instances of internal solitary waves.

Physical Characteristics and Evolution of a Long-Lasting Mesoscale Cyclonic Eddy in the Straits of Florida

Ocean eddies along the Loop Current (LC)/Florida Current (FC) front have been studied for decades, yet studies of the entire evolution of individual eddies are rare. Here, satellite altimetry and ocean color observations, Argo profiling float records and shipborne acoustic Doppler current profiler (ADCP) measurements, together with high-resolution simulations from the global Hybrid Coordinate Ocean Model (HYCOM) are used to investigate the physical and biochemical properties, 3-dimensional (3-D) structure, and evolution of a long-lasting cyclonic eddy (CE) in the Straits of Florida (SoF) along the LC/FC front during April–August 2017. An Angular Momentum Eddy Detection Algorithm (AMEDA) is used to detect and track the CE during its evolution process. The long-lasting CE is found to form along the eastern edge of the LC on April 9th, and remained quasi-stationary for about 3 months (April 23 to July 15) off the Dry Tortugas (DT) until becoming much smaller due to its interaction with the FC and topography. This frontal eddy is named a Tortugas Eddy (TE) and is characterized with higher Chlorophyll (Chl) and lower temperature than surrounding waters, with a mean diameter of ∼100 km and a penetrating depth of ∼800 m. The mechanisms that contributed to the growth and evolution of this long-lasting TE are also explored, which reveal the significant role of oceanic internal instability.

An Intercomparison of Global Reanalysis Products for Southern Africa’s Major Oceanographic Features

A regional assessment of three global ocean reanalysis products is presented for southern Africa’s major oceanographic features. The reanalyses include Mercator Ocean’s Global Reanalysis (GLORYS), the Commonwealth Scientific and Industrial Research Organisation’s (CSIRO) Bluelink Reanalysis (BRAN) and the Fleet Numerical Meteorology and Oceanography Center’s (FNMOC) global Hybrid Coordinate Ocean Model (HYCOM) reanalysis. The aim is to provide modelers with sufficient information for selecting the appropriate product for use as boundary conditions to force their regional ocean models, as well as to provide marine industries, relevant government agencies and academics with insight into the optimal reanalysis product for their purposes. The reanalyses are compared to both assimilated and independent observational datasets spanning various regions within the southern African marine environment. While all reanalysis products reproduce the eastern and western boundary current systems surrounding southern Africa, limitations exist. BRAN outperforms the other reanalyses in its representation of the Mixed Layer Depth, contributing to its good representation of coastal SSTs in the Benguela upwelling system, whereas GLORYS and HYCOM’s misrepresented MLD result in significant warm biases in this region. The Angola-Benguela Frontal Zone and it’s variability is best reproduced by BRAN and HYCOM. The Agulhas Current system’s major components are well reproduced by both GLORYS and BRAN. HYCOM, however, simulates considerably more early retroflections than are observed which have resulted in its mean eastward location. While all the reanalyses overestimate the occurrence of Agulhas meanders, GLORYS and BRAN resolve the associated variability best. Agulhas Current transport is best resolved by GLORYS, unlike BRAN and HYCOM which largely overestimate the magnitude of its south-westward flow, linked to their misrepresentation of the Current’s vertical structure. The bay-scale and nearshore evaluations highlighted issues pertaining to the resolution of the reanalyses and their use at such a small scale. The reanalyses are limited by their resolution, as well as by their misrepresentation of submesoscale processes or lack thereof, prompting the need for the development of regional downscaled models in and around the southern African oceans based on the global ocean reanalysis products.

Intraseasonal variability of the Indonesian throughflow associated with the Madden-Julian Oscillation

Intraseasonal variability (20–90 days) of the Indonesian Throughflow (ITF), which is primarily forced by the Madden-Julian Oscillation (MJO), is investigated using a high-resolution global ocean reanalysis and satellite altimeter data. Previous studies show that during the MJO active phase, downwelling Kelvin waves generated in the central equatorial Indian Ocean propagate along the coast of Sumatra and Java islands, affecting the ITF transport at exit passages in the Indonesian Seas and Makassar Strait. However, the intraseasonal variation of the ITF transport through these straits over the MJO life cycle, especially during the suppressed phase, has not been quantified. To quantify the ITF transport associated with the MJO, composites of ITF transport through major straits in the Indonesian Seas are constructed using a 0.08° global HYbrid Coordinate Ocean Model (HYCOM) reanalysis. A prominent reduction of ITF transport through major straits is found during the MJO active phase, and a transport enhancement comparable to the reduction is evident during the suppressed phase. As a result, the net effect of the MJO on the mean ITF transport is very small due to the cancellation of the enhancement and reduction. The magnitude of the MJO-associated ITF transport variation through all ITF exit passages is about 6 Sv, which is about 50% of the total transport. While the propagation of coastal Kelvin waves during the MJO active phase and their impact on ITF transport is clearly evident in the composite, such Kelvin wave influence on ITF transport is not clearly detected during the suppressed phase. This suggests that anomalous winds over the Maritime Continent (MC) area are mostly responsible for the ITF variation during the suppressed phase in many MJO events. Yet, during some MJO events, remotely-forced upwelling Kelvin waves and their significant impact on the ITF transport are evident during the MJO suppressed phase.

The impact of internal tides and gravity waves on key acoustic parameters and acoustic transmission loss near the Amazon shelf

Upper ocean variability depends on tides and wind forcing. The barotropic tide creates baroclinic internal tides and internal gravity waves (IGWs) that can propagate long distances and dissipate, impacting mixed layer and thermocline depths, vertical thermal gradients, and mixing. In the last decade, there has been an effort to include tidal motions in global HYCOM (HYbrid Coordinate Ocean Model). Here, we look at the impacts of tidally driven motions on acoustic propagation along the Amazon shelf by comparing acoustic transmission loss and key acoustic parameters in 1/25 deg Global HYCOM model runs with and without tides. Acoustic runs were performed at mid-frequencies using Bellhop3D at depths above and below the sonic-layer depth (SLD). Key acoustic parameters assessed include the SLD and vertical sound-speed gradients. Although the impacts of the internal waves and tides on acoustics are intermingled with the effects of other stochastic ocean processes (e.g., mesoscale eddies, fronts, and atmospheric forcing), timeseries analysis indicates that there is acoustic sensitivity to tidally forced variability. These ocean-acoustic results can be used for future development of machine learning algorithms to model coupled oceanographic and acoustical processes.

Simulated zonal current characteristics in the southeastern tropical Indian Ocean (SETIO)

Abstract. Detailed ocean currents in the southeastern tropical Indian Ocean adjacent to southern Sumatran and Javan coasts have not been fully explained because of limited observations. In this study, zonal current characteristics in the region have been studied using simulation results of a 1/8∘ global hybrid coordinate ocean model from 1950 to 2013. The simulated zonal currents across three meridional sections were then investigated using an empirical orthogonal function (EOF), where the first three modes account for 75 %–98 % of the total variance. The first temporal mode of EOF is then investigated using ensemble empirical mode decomposition (EEMD) to distinguish the signals. This study has revealed distinctive features of currents in the South Java Current (SJC) region, the Indonesian Throughflow (ITF)–South Equatorial Current (SEC) region, and the transition zone between these regions. The vertical structures of zonal currents in southern Java and offshore Sumatra are characterized by a one-layer flow. Conversely, a two-layer flow is observed in the nearshore and transition regions of Sumatra. Current variation in the SJC region has peak energies that are sequentially dominated by semiannual, intraseasonal, and annual timescales. Meanwhile, the transition zone is characterized by semiannual and intraseasonal periods with pronounced interannual variations. In contrast, interannual variability associated with El Niño–Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD) modulates the prominent intraseasonal variability of current in the ITF–SEC region. ENSO has the strongest influence at the outflow ITF, while the IOD's strongest influence is in southwestern Sumatra, with the ENSO (IOD) leading the current by 4 months (1 month). Moreover, the contributions (largest to smallest) of each EEMD mode at the nearshore of Java and offshore Sumatra are intraseasonal, semiannual, annual, interannual, and long-term fluctuations. The contribution of long-term variation (19.2 %) in the far offshore eastern Indian Ocean is larger than the interannual (16.3 %) and annual (14.7 %) variations. Future studies should be conducted to investigate this long-term variation.

The features of the coastal fronts in the Eastern Guangdong coastal waters during the downwelling-favorable wind period

A coastal front was detected in the eastern Guangdong (EGD) coastal waters during a downwelling-favorable wind period by using the diffuse attenuation coefficient at 490 nm (Kd(490)). Long-term satellite data, meteorological data and hydrographic data collected from 2003 to 2017 were jointly utilized to analyze the environmental factors affecting coastal fronts. The intensities of the coastal fronts were found to be associated with the downwelling intensity. The monthly mean Kd(490) anomalies in shallow coastal waters less than 25 m deep along the EGD coast and the monthly mean Ekman pumping velocities retrieved by the ERA5 dataset were negatively correlated, with a Pearson correlation of − 0.71. The fronts started in October, became weaker and gradually disappeared after January, extending southwestward from the southeastern coast of Guangdong Province to the Wanshan Archipelago in the South China Sea (SCS). The cross-frontal differences in the mean Kd(490) values could reach 3.7 m−1. Noticeable peaks were found in the meridional distribution of the mean Kd(490) values at 22.5°N and 22.2°N and in the zonal distribution of the mean Kd(490) values at 114.7°E and 114.4°E. The peaks tended to narrow as the latitude increased. The average coastal surface currents obtained from the global Hybrid Coordinate Ocean Model (HYCOM) showed that waters with high nutrient and sediment contents in the Fujian and Zhejiang coastal areas in the southern part of the East China Sea could flow into the SCS. The directions and lengths of the fronts were found to be associated with the flow advection.

Impact of mesoscale eddies on the source funnel of sediment trap measurements in the South China Sea

The mooring tethered time-series sediment traps (TS-traps) collect sinking particles in the ocean, enabling the estimation of biological pump (BP) efficiency based on the assumption of vertical settling. However, the advection of sea water can disperse particles over long distance during the sinking process, introducing uncertainties into the estimates of BP efficiency. In the South China Sea, mesoscale eddies generate significant lateral transport above the mesopelagic zone, which has a large impact on the collecting area and source funnel of TS-traps. We studied the role of eddies in the lateral transport of sinking particles using a Lagrangian Particle Tracking Model (LPTM) which was forced by three-dimensional currents of eddy composites derived from the global Hybrid Coordinate Ocean Model (HYCOM). The eddy reanalysis using the HYCOM outputs was validated using occurrence, genesis and dissipation statistics of eddies. Backward particle tracking was performed at three stations located in the northern, central and southwestern SCS, respectively. The results indicate that the vertical sinking assumption is invalid when the particle sinking velocity is relatively low. Furthermore, the fraction of surface sources from the edge of the eddies depends on both the sinking velocity and seasonal variability of the background current. The eddy current controls the movement of particles within the eddy, and many of the particles with a sinking velocity <80 m d−1 can be transported from the edge to the eddy center. The source funnel within the mesoscale eddy is tilted vertically, leading to an asymmetrical distribution of the particle source around the trap location. Finally, for slow-sinking particles, trap-collected particles could originate from the subsurface located tens to hundreds of kilometers away from the trap location. Overall, our results highlight the role of mesoscale eddies in interpreting trap measurements.

A modelling synthesis of the volume flux through Cook Strait

ABSTRACT We present hindcasts of currents in Cook Strait from a baroclinic, tide-resolving ocean model. Three hindcasts are run, forced respectively at the lateral boundaries with data from the Bluelink Reanalysis (BRAN), the global Hybrid Coordinate Ocean Model (HYCOM) and a New Zealand Shelf Seas baroclinic hindcast (NZROMS). The hindcast forced by BRAN is evaluated against measured currents at locations in the Cook Strait Narrows and South Taranaki Bight. The model reproduces the observed tidal, subtidal and mean currents well. The three hindcasts produce a 3-year-mean volume flux southward through the Narrows of 0.360, 0.734 and 0.230 Sv, respectively. When corrected for the observed discrepancies between model and measurements in the Narrows, the estimates are 0.430, 0.498 and 0.337 Sv. We conclude that true 3-year mean volume flux through Cook Strait is 0.42 ± 0.08 Sv. The tidal volume flux amplitude is 4.68 Sv (M2). The subtidal fluctuations in volume flux have a standard deviation of 0.62 Sv. The effect of wind is investigated with a wind-forced barotropic model, which produces volume flux fluctuations that are highly correlated (r = 0.93) with the subtidal fluctuations from the baroclinic model, implying that the latter are largely wind-generated.

On the interplay between horizontal resolution and wave drag and their effect on tidal baroclinic mode waves in realistic global ocean simulations

The effects of horizontal resolution and wave drag damping on the semidiurnal M2 tidal energetics are studied for two realistically-forced global HYbrid Coordinate Ocean Model (HYCOM) simulations with 41 layers and horizontal resolutions of 8 km (1∕12.5∘; H12) and 4 km (1∕25∘; H25). In both simulations, the surface tidal error is minimized by tuning the strength of the linear wave drag, which is a parameterization of the surface-tide energy conversion to the unresolved baroclinic wave modes. In both simulations the M2 surface tide error with TPXO8-atlas, an altimetry constrained model, is 2.6 cm. Compared to H12, the surface tide energy conversion to the resolved vertical modes is increased by 50% in H25. This coincides with an equivalent reduction in the tuned loss of energy from the surface tide to the wave drag. For the configurations studied here, the horizontal and not the vertical resolution is the factor limiting the number of vertical modes that are resolved in most of the global ocean: modes 1–2 in H12 and modes 1–5 in H25. The wave drag also dampens the resolved internal tides. The 40% reduction in wave-drag strength does not result in a proportional increase in the mode-1 energy density in H25. In the higher-resolution simulations, topographic mode-scattering and wave–wave interactions are better resolved. This allows for an energy flux out of mode 1 to the higher modes, mitigating the need for an internal tide damping term. The HYCOM simulations are validated with analytical conversion models and altimetry-inferred sea-surface height, fluxes, and surface tide dissipation. H25 agrees best with these data sets to within ∼10%. To facilitate the comparison of stationary tide signals extracted from time series with different durations, we successfully apply a spatially-varying correction factor.

Origin and formation of the Ryukyu Current revealed by HYCOM reanalysis

The origin of the Ryukyu Current (RC) and the formation of its subsurface velocity core were investigated using a 23-year (1993–2015) global Hybrid Coordinate Ocean Model (HYCOM) dataset. The volume transport of the RC comes from the Kuroshio eastward branch (KEB) east of Taiwan and part of the North Pacific Subtropical Gyre (p-NPSG). From the surface to 2 000 m depth, the KEB (p-NPSG) transport contributes 41.5% (58.5%) to the mean total RC transport. The KEB originally forms the subsurface velocity core of the RC east of Taiwan due to blockage of the subsurface Kuroshio by the Ilan Ridge (sill depth: 700 m). Above 700 m, the Kuroshio can enter the East China Sea (ECS) over the Ilan Ridge, meanwhile, the blocked Kuroshio below 700 m turns to the right and flows along the Ryukyu Islands. With the RC flowing northeastward, the p-NPSG contribution strengthens the subsurface maximum structure of the RC owing to the blockage of the Ryukyu Ridge. In the surface layer, the p-NPSG cannot form a stable northeastward current due to frequent disturbance by mesoscale eddies and water exchange through the gaps (with net volume transport into ECS) between the Ryukyu Islands.

Long-Term Trend of Green and Golden Tides in the Eastern Yellow Sea

Kim, K.; Shin, J.; Kim, K.Y., and Ryu, J.H., 2019. Long-term trend of green and golden tide in the eastern Yellow Sea. In: Jung, H.-S.; Lee, S.; Ryu, J.-H., and Cui, T. (eds.), Advances in Remote Sensing and Geoscience Information Systems of Coastal Environments. Journal of Coastal Research, Special Issue No. 90, pp. 317-323. Coconut Creek (Florida), ISSN 0749-0208.Since 2008, floating green tides (Ulva sp.) have been occurring continuously in the Yellow Sea (YS), and after 2013 floating golden tides (Sargassum sp.) have also occurred. The distribution, areal coverage, and migration of floating green tides have been actively studied, but most research has focused only on the western YS. Little is known about the floating golden tides in the eastern YS. The purpose of this study was to determine the long-term distribution of floating green and golden tides in the eastern YS using Geostationary Ocean Color Imager (GOCI), Moderate Resolution Imaging Spectroradiometer (MODIS), and Landsat satellite images from 2008 to 2017. In addition, the migration of floating macroalgae with Global Hybrid Coordinate Ocean Model (HYCOM) surface current data were compared. Green tides were observed in 2008, 2009, 2011, 2015, and 2016 in the eastern YS. From a satellite image backtracking analysis, it was confirmed that the green tides observed in the eastern YS were supplied from the western YS. When the maximum areal coverage of green tide was compared between the eastern and western YS, the coverage in the eastern YS was found to be about 4 % of that of the western YS. However, in 2011, the largest amount of floating macroalgae was found in the eastern YS and it accounted for about 45 % of the amount in the western part of the YS. In the eastern YS, floating golden tides were found in 2013, 2015, 2016, and 2017, with the largest amount of floating macroalgae occurring in 2017. Although there were no long-term golden tide data for the western YS, such that it was not compared to the areal coverage of the eastern YS, it was confirmed that the amount of golden tide supplied to the eastern YS gradually increased. A comparison between the migration of floating macroalgae and HYCOM surface current data suggested that the migration and flow directions were not identical, and were considerably affected by surface ocean currents during their passage into the eastern YS. From this study, the long-term distribution and changes in areal coverage of green and golden tides in the eastern YS were obtained for the first time. This information will be useful for understanding the long-term patterns of green and golden tides, and provides basic data for predicting the occurrence and migration of floating macroalgae.

What controls the extreme flow through the Kerama Gap: a global HYbrid Coordinate Ocean Model reanalysis point of view

The temporal variability of volume transport from the North Pacific Ocean to the East China Sea (ECS) through the Kerama Gap (between Okinawa Island and Miyakojima Island—a part of Ryukyu Island Arc) is investigated using a 20-year global HYbrid Coordinate Ocean Model reanalysis with the Navy Coupled Ocean Data Assimilation from 1993 to 2012. We study the causes of extreme flow events defined as the volume transport through the Kerama Gap exceeding one standard deviation (4.50 Sv) above (inflow, from the North Pacific into the ECS) or below (outflow, from the ECS to the North Pacific) the mean volume transport (1.95 Sv). In our results, the impinging mesoscale eddies from the Pacific Ocean are the most important factor in controlling the extreme flow occurrence, explaining 70% (68%) of the mean extreme inflow (outflow)-mean volume transport anomaly. Extreme inflow (outflow) through the Kerama Gap is associated with the impinging anticyclonic (cyclonic) eddies along the eastern side of Okinawa and/or cyclonic (anticyclonic) eddies to the southwest of the gap. The Kuroshio frontal meanders cause the mean Kuroshio axis to shift away from (toward) the Kerama Gap during the extreme inflow (outflow)-mean events and play a secondary role in determining the extreme flow occurrence, representing 15% (14%) of the mean extreme inflow (outflow)-mean volume transport anomaly. The interannual variability of the flow through the gap also plays a subsidiary role, explaining another 13% (14%) of the mean extreme inflow (outflow)-mean volume transport anomaly.

Assessment of shelf sea tides and tidal mixing fronts in a global ocean model

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 <200 m) the area-weighted root mean square error of the M2 sea surface amplitude in HYCOM represents 35% of the 50 cm root mean squared M2 sea surface amplitude when compared to satellite constrained models TPXO8 and FES2014. HYCOM and the altimeter constrained tidal models TPXO8 and FES2014 exhibit similar skill in reproducing barotropic tidal currents estimated from in-situ current meter observations. Through comparison of a global HYCOM simulation with tidal forcing to a global HYCOM simulation with no tides, and also to previous regional studies of tidal mixing fronts in shelf seas, we demonstrate that HYCOM with embedded tides exhibits quite high skill in reproducing known tidal mixing fronts in shelf seas. Our results indicate that the amount of variability in the location of the tidal mixing fronts in HYCOM, estimated using the Simpson-Hunter parameter, is consistent with previous studies when the differences in the net downward heat flux, on a global scale, are taken into account. We also provide evidence of tidal mixing fronts on the North West Australian Shelf for which we have been unable to find references in the existing scientific literature.

Assimilation of Sea Surface Temperature in a Global Hybrid Coordinate Ocean Model

The Hybrid Coordinate Ocean Model (HYCOM) uses different vertical coordinate choices in different regions. In HYCOM, the prognostic variables include not only the seawater temperature, salinity and current fields, but also the layer thickness. All prognostic variables are usually adjusted in the assimilation when multivariate data assimilation methods are used to assimilate sea surface temperature (SST). This paper investigates the effects of SST assimilation in a global HYCOM model using the Ensemble Optimal Interpolation multivariate assimilation method. Three assimilation experiments are conducted from 2006-08. In the first experiment, all model variables are adjusted during the assimilation process. In the other two experiments, the temperature alone is adjusted in the entire water column and in the mixed layer. For comparison, a control experiment without assimilation is also conducted. The three assimilation experiments yield notable SST improvements over the results of the control experiment. Additionally, the experiments in which all variables are adjusted and the temperature alone in all model layers is adjusted, produce significant negative effects on the subsurface temperature. Also, they yield negative effects on the subsurface salinity because it is associated with temperature and layer thickness. The experiment adjusting the temperature alone in the mixed layer yields positive effects and outperforms the other experiments. The heat content in the upper 300 m and 300-700 m layers further suggests that it yields the best performance among the experiments.

The wind work input into the global ocean revealed by a 17-year global HYbrid coordinate ocean model reanalysis

A 17-year (1999–2015) HYbrid Coordinate Ocean Model (HYCOM) reanalysis is used to calculate the wind work to the global oceanic circulation with a focus on the global integral of the wind work on the ageostrophic currents. While the wind work on the geostrophic currents estimate of the study is in good agreement with previous research, our estimation of the wind work on the ageostrophic currents (5.44 TW) is significantly larger than previous estimates. The two main reasons for this difference are: (1) the thickness of HYCOM's first vertical layer is 1 m and thus the surface ageostrophic currents are better represented than those in numerical models that use a thicker top layer; and (2) the reanalysis is forced with hourly surface stresses and thus contains wind work contributions from high frequency wind stress variations compared against previous estimates using daily wind stress.

Observation and prediction of flotsam trajectories in the California Current System based on surface drift of RAFOS floats

Surface drift of eighty-nine undrogued RAFOS floats in the California Current System (CCS) has been studied to describe patterns of flotsam drift, its seasonal variability and predictability limitations. The floats were launched in the California Undercurrent during 1992–2010 and were tracked by the ARGOS system when they surfaced at the end of their subsurface missions. The duration of surface trajectories varied from as short as 11 to as long as 280 days. The surface drift of these floats was typically equatorward in the California Current. However, some floats moved poleward into the Subpolar Gyre, and others drifted westward into the North Equatorial Current.Usually, observations of surface currents use drifters which are coupled to the surface layer by drogues located at 15 m depth. While drogued observations are useful for studies of circulation of the upper layer of the ocean, a more typical operational problem involves trying to find flotsam that has fallen off the deck of a ship or to predict the path of an abandoned vessel. To better understand the behavior of these surface drifting objects, observations of the surface drifts of RAFOS floats in the California Current system were used to compare the floats’ motions to wind-induced drift and to evaluate the drift prediction by three ocean models: Ocean Surface Current Simulator (OSCURS), Global Navy Coastal Model (gNCOM) and Hybrid Coordinate Ocean Model (HYCOM).Cross-Calibrated Multi-Platform Ocean Surface Wind Vector L3.0 First-Look Analysis wind data were used to analyze the surface drift of the RAFOS floats. The best correlation between float drift speed and wind speed was observed during summer and fall months, when the wind regime over the Northeast Pacific is dominated by the North Pacific High. The predominant winds associated with this regime are steady and contribute to the surface drift in the same direction as the large-scale geostrophic component of the California Current.Evaluation of the drift prediction by three ocean models was conducted by comparing observed drifter trajectories with model-simulated trajectories at 7-day time scales. The model-simulated trajectories were initially collocated with RAFOS positions and restarted every 15 days. The OSCURS model was the only one of three analyzed models which allowed for tuning based on the surface RAFOS trajectories. Two main parameters, Wind Angle Deflection (WAD) and Geostrophic Current Factor (GCF), were adjusted to minimize the separation between observed and modeled trajectories. At shorter time scales, from weeks to several months, the tuned OSCURS model outperformed gNCOM and HYCOM models in reproducing the observed surface drift. For seasonal and longer time scales, both gNCOM and HYCOM can reproduce the seasonal variability, flow direction, and eddy-like structures in the California Current system, which may be useful for long-term forecasts of the Northeast Pacific circulation.

Seasonal and intra-seasonal variations of suspended-sediment distribution in the Yellow Sea

Seasonal and intra-seasonal variability of hydrodynamics and suspended sediment concentration (SSC) in the Yellow Sea were investigated using combined ship measurements and physical fields extracted from the global Hybrid Coordinate Ocean Model (HyCOM). Model-data comparison indicates that the HyCOM solutions fairly agree with the observed temperature fields, both showing strong horizontal and vertical seasonality. High SSC was observed in waters near the Chengshantou Cape and above the old Yellow River subaqueous delta in November 2012, and was coincident with enhanced velocity shear. Thermal fronts were identified in waters less than 60m deep with their prominent appearances in November 2012, inhibiting the offshore sediment transport. The previously reported depo-center of distal Yellow River mud-patch off the Chengshantou Cape was confined between the positions of the strong velocity shear in a shoreward direction and the thermal fronts in an offshore direction. We concluded that the nearshore mixing front and the offshore thermal front play critical roles in both capture of the fine-grained suspended sediment and formation of the Omega-shaped mud patch.

Development, implementation, and validation of a California coastal ocean modeling, data assimilation, and forecasting system

A three-dimensional, near real-time data-assimilative modeling system for the California coastal ocean is presented. The system consists of a Regional Ocean Modeling System (ROMS) forced by the North American Mesoscale Forecast System (NAM). The ocean model has a horizontal resolution of approximately three kilometers and utilizes a multi-scale three-dimensional variational (3DVAR) data assimilation methodology. The system is run in near real-time to produce a nowcast every six hours and a 72-hour forecast every day. The performance of this nowcast system is presented using results from a six-year period of 2009–2015.The ROMS results are first compared with the assimilated data as a consistency check. RMS differences in observed satellite infrared sea surface temperatures (SST) and vertical profiles of temperature between observations and ROMS nowcasts were found to be mostly less than 0.5°C, while the RMS differences in vertical profiles of salinity between observations and ROMS nowcasts were found to be 0.09 or less. The RMS differences in SST show a distinct seasonal cycle that mirrors the number of observations available: the nowcast is less skillful with larger RMS differences during the summer months when there are less infrared SST observations due to the presence of low-level clouds. The larger differences during summer were found primarily along the northern and central coasts in upwelling regions where strong gradients exist between colder upwelled waters nearshore and warmer offshore waters. RMS differences between HF radar surface current observations and ROMS nowcasts were approximately 7–8cms−1, which is about 30% of the time mean current speeds in this region. The RMS differences in sea surface height (SSH) between the AVISO (Archiving, Validation and Interpretation of Satellite Oceanographic) altimetric satellite observations and ROMS nowcasts were about 2cm. In addition, the system realistically reproduces the interannual variability in temperatures at the M1 mooring (122.03°W, 36.75°N) in Monterey Bay, including the strong warming of the California coastal ocean during 2014.The ROMS nowcasts were then validated against independent observations. A comparison of the ROMS nowcast with independent profile observations of temperature and salinity shows RMS differences of 0.7 to 0.92°C and 0.13 to 0.17, which are larger (by up to a factor of 2) than the differences found in the comparisons with assimilated data. Validation of the depth-averaged currents derived from Spray gliders shows that the flow patterns associated with California Current and California Undercurrent/Davidson current systems and their seasonal variations are qualitatively reproduced by the ROMS modeling system.Lastly, the impact of two recent upgrades to the system is quantified. Switching the lateral boundary conditions from a U.S. west coast regional model to the global HYCOM (HYbrid Coordinate Ocean Model) model results in an improvement in the simulation of the seasonal and interannual variations in the SSH, especially south of Pt. Conception (120.47°W, 34.45°N). The assimilation of altimetric satellite SSH data also results in an improvement in the model surface currents when compared to independent surface drifter observations.