The Kuroshio near the Luzon Strait and circulation in the northern South China Sea during August and September 1994

Wind data from NCEP and hydrographic data obtained from August 28 to September 10, 1994 have been used to compute circulation in the northern South China Sea and near Luzon Strait using three-dimensional diagnostic models with a modified inverse method. The numerical results are as follows: the main Kuroshio is located above 400 m levels near Taiwan’s eastern coast and above 800 m levels away from it. Near Luzon Strait above 400 m levels a branch of the Kuroshio joins with a part of the northward current, which comes from an area west of Luzon’s western coast and intrudes northwestward, then it branchs into western and eastern parts near 20°30′ N. The eastern part flows northward into an area east of Taiwan, while its western part continues to intrude northwestward, flowing through an area southwest of Taiwan. Net westward intruded volume transport through longitude Section AB at 121°00′ E from 19°00′ N to 21° 43′ N is about 3.5 × 10⁶ m³s⁻¹ in a layer above 400 m levels. The anticyclonic eddies W1 and W3 exist above 700 m levels east of Dongsha Islands and below 200 m levels in the eastern part of the region, respectively. The circulation in the middle region is dominated mainly by a basin-scale cyclonic gyre, and consists of three cyclonic eddies. Strong upwelling occurs in the middle region. The joint effect of baroclinity and relief and interaction between wind stress and relief both are important for real forcing of flow across contours of fH ⁻¹ in effecting the circulation pattern.     

Modeling the global circulation response and the regional response of the Arctic Ocean to the external forcing anomalies

The problem of numerical modeling and analysis of the large-scale World Ocean circulation variability under variations of the external forcing is considered. A numerical model was developed in the INM RAS and is based on the primitive equations of the ocean circulation written in a spherical generalized σ-coordinate system. The model’s equations are approximated on a grid with resolution of 2.5° × 2° × 33, and the North Pole is displaced to the continental point (60°E, 60.5°N). There are two stages for the numerical experiments. The quasi-equilibrium circulation of the World Ocean under the climatological atmospheric forcing is simulated at the first stage. The run is carried out over a period of 3000 years during which a quasi-equilibrium model regime is formed. At the second stage, the sensitivity of the model ocean circulation to the atmospheric forcing perturbations in the Southern Hemisphere is studied. According to the results, the strongest regional changes in the hydrography take place in the Arctic Ocean. Substantial changes of sea’s surface height and local anomalies of the temperature and salinity are formed there.     

Technique for Simulation of Black Sea Circulation with Increased Resolution in the Area of the IO RAS Polygon

A numerical technique is presented for simulating the hydrophysical fields of the Black Sea on a variable-step grid with refinement in the area of IO RAS polygon. Model primitive equations are written in spherical coordinates with an arbitrary arrangement of poles. In order to increase the horizontal resolution of the coastal zone in the area of the IO RAS polygon in the northeastern part of the sea near Gelendzhik, one of the poles is placed at a land point (38.35° E, 44.75° N). The model horizontal resolution varies from 150 m in the area of the IO RAS polygon to 4.6 km in the southwestern part of the Black Sea. The numerical technique makes it possible to simulate a large-scale structure of Black Sea circulation as well as the meso- and submesoscale dynamics of the coastal zone. In order to compute the atmospheric forcing, the results of the regional climate model WRF with a resolution of about 10 km in space and 1 h in time are used. In order to demonstrate the technique, Black Sea hydrophysical fields for 2011–2012 and a passive tracer transport representing self-cleaning of Gelendzhik Bay in July 2012 are simulated.     

Investigation on the Cyclonic Seastate along Southeast Coast of India

The east coast of the Indian Peninsula experiences the effects of a devastating cyclone at least annually. The Thane cyclone of 29–30 December 2011 has been once such event that resulted in significant damages along the coastline of Tamil Nadu on the southeast coast of India (13° 9′ 10′ N and 80° 21′E). Waves as high as 8–12 m in a water depth of 20 m have been measured. Such huge waves, combined with a storm surge of 0.5 m, lead to severe damages to coastal structures during the passage of the cyclone. As a part of an exercise in assessing the sediment transport rates through measurements of the hydrodynamic driving parameters along the coast of major port of Chennai instruments were deployed for the measurement of waves and flow field. The measurement campaign was carried out at a location of about 120 km north of the cyclone made landfall. The ENCEP wind data formed the input for executing the WAM model for the simulation of wave characteristics, which are compared with the measured wave data. The agreement between them is found to be good. The details of the analysis of the results are presented and discussed in this paper.     

Upper Ocean Four-Dimensional Variational Data Assimilation in the Arabian Sea and Bay of Bengal

A regional ocean circulation model with four-dimensional variational data assimilation scheme is configured to study the ocean state of the Indian Ocean region (65°E–95°E; 5°N–20°N) covering the Arabian Sea (AS) and Bay of Bengal (BoB). The state estimation setup uses 10 km horizontal resolution and 5 m vertical resolution in the upper ocean. The in-situ temperature and salinity, satellite-derived observations of sea surface height, and blended (in-situ and satellite-derived) observations of sea surface temperature alongwith their associated uncertainties are used for data assimilation with the regionally configured ocean model. The ocean state estimation is carried out for 61 days (1 June to 31 July 2013). The assimilated fields are closer to observations compared to other global state estimates. The mixed layer depth (MLD) of the region shows deepening during the period of assimilation with AS showing higher MLD compared to the BoB. An empirical forecast equation is derived for the prediction of MLD using the air–sea forcing variables as predictors. The surface and sub-surface (50 m) heat and salt budget tendencies of the region are also investigated. It is found that at the sub-surface, only the advection and diffusion temperature and salt tendencies are important.     

Hydrography and circulation of the bay of Bengal during withdrawal phase of the southwest Monsoon

Hydrographic data were collected from 3 to 10 September 1996 along two transects; one at 18° N and the other at 90° E. The data were used to examine the thermohaline, circulation and chemical properties of the Bay of Bengal during the withdrawal phase of the southwest monsoon. The surface salinity exhibited wide spatial variability with values as low as 25.78 at 18° N / 87° E and as high as 34.79 at 8° N / 90° E. Two high salinity cells (S > 35.2) were noticed around 100 m depth along the 90° E transect. The wide scatter in T-S values between 100 and 200 m depth was attributed to the presence of the Arabian Sea High Salinity (ASHS) water mass. Though the warm and low salinity conditions at the sea surface were conducive to a rise in the sea surface topography at 18° N / 87° E, the dynamic height showed a reduction of 0.2 dyn.m. This fall was attributed to thermocline upwelling at this location. The geostrophic currents showed alternating flows across both the transects. Relatively stronger and mutually opposite currents were noticed around 25 m depth across the 18° N transect with velocity slightly in excess of 30 cm s⁻¹. Similar high velocity (> 40 cm s⁻¹) pockets were also noticed to extend up to 30 m depths in the southern region of the 90° E transect. However, the currents below 250 m were weak and in general < 5 cm s⁻¹. The net geostrophic volume transports were found to be of the order of 1.5 × 10⁶ m³ s⁻¹ towards the north and of 6 × 10⁶ m³ s⁻¹ towards west across the 18° N and 90° E transects respectively. The surface circulation patterns were also investigated using the trajectories of drifting buoys deployed in the eastern Indian Ocean around the same observation period. Poleward movement of the drifting buoy with the arrival of the Indian Monsoon Current (IMC) at about 12° N along the eastern rim of the Bay of Bengal has been noticed to occur around the beginning of October. The presence of an eddy off the southeast coast of India and the IMC along the southern periphery of the Bay of Bengal were also evident in the drifting buoy data.     

The Kuroshio Extension Front from satellite sea surface temperature measurements

The position and strength of the surface Kuroshio Extension Front (KEF), defined as the sea surface temperature (SST) gradient maximum adjacent to the Kuroshio Extension (KE) axis (approximated by a specific SSH contour consistently located at, or near, the maximum of the SSH gradient magnitude), have been studied using weekly, microwave SST measurements from the later 1997 to early 2008. The mean KEF meanders twice around ∼36°N between the east coast of Japan and 153°E. It then migrates southeast to ∼34°N, just before reaching the Shatsky Rise (∼158°E), then progresses mostly eastward. Spatially, the KEF is strongest near the Japan coast, while it is seasonally strongest in winter and weakest in summer. Low-frequency variations of its strength, most notably in its upstream region, can be related to the known bimodal states of the KE. During 2003–2005, when the KE was in its stable state, the winter KEF SST gradient exceeded 10°C/100 km.     

Bathymetry Inversion with the Gravity-Geologic Method: A Study of Long-Wavelength Gravity Modeling Based on Adaptive Mesh

Modeling of long-wavelength gravity anomaly is crucial for bathymetry inversion with a gravity-geologic method. We propose a new method, named as iGGM, to approximate the long-wavelength gravity anomalies by using a finite element method based on an adaptive triangular mesh which is constructed by uneven control points. The mesh size is suitably controlled to ensure that there are several control points in each grid. By using iGGM, the bathymetry in the South China Sea (Test Area #1: 112°E–119°E, 12°N–20°N) and East China Sea (Test Area #2: 125°E–130°E, 25°N–30°N) is estimated. The performance of the method was evaluated by comparing the predictions with Earth topographical database 1 (ETOPO1) model and shipborne depths in the test points. Results show that the depths derived by iGGM have a strong correlation with the shipborne depths. In the test points, the mean values of their differences are smaller than 10 m. The standard deviations of their differences are smaller than 120 m and their correlation is stronger than 0.98. Meanwhile, the results provided by the iGGM model are comparable with that obtained by the ETOPO1 model, e.g., the differences between iGGM and ETOPO1 models in test points for Test Areas 1 and 2 are 116 and 70 m in standard deviation, respectively.     

Numerical study of sea waves created by tropical cyclone Jelawat

A numerical study of sea waves generated by tropical cyclone Jelawat is carried out using the cycle 4 version of the WAM.The model domain currently covers the latitudes 20-45 N and longitudes 115-135 E,and the model spatial resolution reaches 0.25 ×0.25.Comparison of the model results with buoy observations reveals that the model can fairly reproduce the temporal variation of observed waves.Two-dimensional comparison is also made against the satellite altimeter significant wave heights derived from TOPEX/PO...     

Major variability modes and wind patterns over the Sea of Japan and adjacent areas

The expansion of wind fields observed at fixed times (four times daily) in complex empirical orthogonal functions is performed for the Japan Sea area (34°–53° N, 127°–143° E). The wind fields are taken from the 1998–2004 NCEP/NCAR Reanalysis data with better spatial resolution (1° × 1°) than the standard product, which are publicly available on the Internet. Major modes of wind variability in the Japan Sea area are identified. The modes determine a general direction of air-mass transport throughout a year, zonal and meridional modulation, and a cyclonic and an anticyclonic eddy component. Objective classification of wind fields with respect to the prevailing flow direction is performed, and wind stress and wind-curl patterns are obtained for major events in the cold and warm periods of the year. The pattern obtained can be used in hydrodynamic numerical models of the general circulation of the Japan Sea.     

Numerical Simulation of North Indian Ocean State Prior to the Onset of SW Monsoon Using SSM/I Winds

The Sea Surface Temperatures (SST) and currents are simulated over the north Indian Ocean, during the onset phase of southwest monsoon for the three years 1994, 1995, and 1996, using daily Special Sensor Microwave/Imager (SSM/I) winds and National Center for Environmental Prediction (NCEP) heat fluxes as forcings in the 2½ layer thermodynamic numerical ocean model. The results are discussed for the 30-day period from 16 May to 13 June for all the three years, to determine the ocean state during the onset phase of SW monsoon. The maximum variability in the simulated SST is found along the Somali coast, Indian coasts, and equatorial regions. The maximum SST in the North Arabian Sea is found to be greater than 30°C and minimum SST in the west equatorial region is 25°C during the onset phase of all three years. Model SSTs are in agreement with Reynolds SST. SST gradients in the north-south as well as in the east-west directions, west of 80°E are found to change significantly prior to the onset. It can be inferred from the study that the SST gradient of 2.5°C/2000 km is seen due north and due west of the region 2° - 7°S, 60° - 65°E, about 8 to 10 days prior to the arrival of SW monsoon near Kerala coast. Upper and lower layer circulation fields do not show prominent interannual variability.     

Tide and surge dynamics along the Iberian Atlantic coast

The eastern Atlantic barotropic dynamics (in a region spanning from 20° N to 48° N and 34° W to 0°) are studied through numerical modelling and in situ measurements. The main source of data is the tidal gauge network REDMAR, managed by Clima Marítimo (Puertos del Estado). The numerical model employed is the HAMSOM, developed both by the Institut für Meereskunde (Hamburg University) and Clima Marítimo.In this paper, tidal and storm surge dynamics are studied for the region, focusing particularly on the nonlinear transfer of energy between the different forcings.The results of tidal simulations show good agreement between semidiurnal harmonic components and the values obtained from the tidal gauges (both coastal and pelagic) and current metres. The nonlinear transfers of energy from semidiurnal to higher order harmonics, such as M4 and M6, were mapped. Those transfers were found to be important in only two areas: The French continental shelf in the Bay of Biscay and the widest part of the African shelf, south of Cabo Bojador. The modelled diurnal constituents show larger relative differences with measurements than semidiurnal harmonics, especially in data concerning the phase.A method to isolate the nonlinear transfers of energy between tidal and atmospheric forcing during a storm surge was developed. These transfers were found to be significant in the same areas where tidal nonlinear activity was present. The effect of short period wind generated waves on sea surface elevation was also investigated. The magnitude of the spatial derivatives of radiation stress was compared with wind-induced stress. As a result of this comparison, we found the inclusion of a forcing term that depends on radiation stress in ocean model simulations at this scale and resolution to be not essential. The effect of computing wind-induced stresses, with a formulation that explicitly depends on sea state, was also explored by means of a coupled run of the HAMSOM and the spectral wave model WAM for a storm surge event in the Spanish coast. This formulation was not found to be an improvement over a classical parameterisation which only depends on wind fields.     

南海中部上层海洋湍流混合的空间分布特征及参数化模型

通过对2010年5月南海16°N和14.5°N断面的湍流微结构剖面观测资料分析,给出了南海海盆上层湍流混合空间分布特征:在16°N断面上,上层10~400m垂向平均湍动能耗散率ερ在东侧略大于西侧;相反,在14.5°N断面上,西侧ερ均值约是东侧ερ的4倍,其中,西侧110.5°~111°E的ερ的平均值为2.6×10-6 W/m3,东侧118.5°E的ερ仅为5.89×10-7 W/m3。通过分析细结构剪切和湍流混合的相关性,发现剪切是南海中部上层强湍流混合的主要驱动力,揭示了高模态内波破碎可能是湍流混合的主要机制。另外,研究了大洋中的3种参数化模型,发现适用于大洋近海的参数化MacKinnon-Gregg(MG)模型能较好地用浮频和剪切估算南海中部深海区上层湍流耗散率。     

A description of short period temperature fluctuations in the upper ocean

Dynamical properties of short-period temperature fluctuations are studied. Water temperature was measured continuously at several depths at the following stations: at 38°29.5′N, 141°35.8′E (100 m depth) on the continental shelf off Miyagi Prefecture in the summer of 1967, at 35°01.8′N, 139°0.8.5′E (100 m depth) in Sagami Bay in the summer of 1968, and at 32°32.2′N, 129°53.7′E (74 m depth) in Tachibana Bay in the summer of 1970. These measurements were made with a thermistor array laid down from the R. V.Tanseimaru (Ocean Research Institute, University of Tokyo) which was fixed with bow and stern anchors. Significant temperature fluctuations found at the first and the third stations are thought to be due to first mode internal waves having amplitude 3 to 5 m and period 5 to 20 minutes. The wave length of the waves is estimated to be 25 m to 400 m from the observed density structure. At the second station, we found second-mode internal waves. The period, amplitude and wave length of the waves are about 30 minutes, 1.3 m and 600 m, respectively. In all cases, the spectral density of the temperature fluctuations decreases with increase in frequency. However, the decrease obey neither the ?3 power law nor the ?5/3 power law. Coherences in the temperature fluctuations between two depths of measurement in the seasonal thermocline are significantly high in the range of frequencies lower than the local Brunt-Väisälä frequency, but are low in the higher frequency range. At the first and the third stations, the difference in the level of coherences between the lower frequencies and the higher frequencies are large. Phase differences between two depths in the thermocline are small in the lower frequency range. This suggests that the first-mode internal waves are predominant over higher-mode internal waves and over other disturbances.     

Intercomparison of High-Resolution Bay of Bengal Circulation Models Forced with Different Winds

The high-resolution Bay of Bengal circulation modeling in the region [80E–95E; 5N–22N] is performed with a horizontal resolution of 10 km and the highest vertical resolution of 5 m near the surface. The intercomparison experiments, with ocean model forced with the near-surface (1) National Centers for Environmental Prediction (NCEP) reanalysis winds and (2) blended seawinds data (a combination of remotely sensed scatterometer and in situ observations) are carried out for a period of 17 years during 1998–2014. The seasonal variability of the realistically simulated surface hydrographic (temperature and salinity) and circulation (currents) variables from both the experiments is compared and contrasted with the observational data. The mixed layer depth seasonal variability of the region is also studied. The mesoscale features of currents at 50 and 100 m are also studied. The volume transport across different sections in the Bay of Bengal is computed and its relation with summer monsoon rainfall is investigated. The results suggest that there is no real advantage of using high-resolution blended seawinds over the much coarser NCEP reanalysis winds.     

Discretization of Complex Geometric Domain Through Computer-Based Stair-Step Representation Method for Estimating Water Levels Associated with A Storm

In this paper, an improvement has been made to the approximation technique of a complex domain through the stairstep approach to have a considerable accuracy, minimize computational cost, and avoid the hardship of manual work. A novel stair-step representation algorithm is used in this regard, where the entire procedure is carried out through our developed MATLAB routine. Arakawa C-grid is used in our approximation with (1/120)° grid resolution. As a test case, the method is applied to approximate the domain covering the area between 15°–23°N latitudes and 85°–95° E longitudes in the Bay of Bengal. Along with the approximation of the land-sea interface, coastal stations are also identified. Approximated land-sea interfaces and coastal stations are found to be in good agreement with the actual ones based on the similarity index, overlap fraction, and extra fraction criteria. The method can be used for approximating an irregular geometric domain to employ the finite difference method in solving problems related to long waves. As a test case, shallow water equations in Cartesian coordinates are solved on the domain of interest for simulating water levels due to the nonlinear tide-surge interaction associated with the storms April 1991 and AILA, 2009 along the coast of Bangladesh. The same input except for the discretized domain and bathymetry as that of Paul et al. (2016) is used in our simulation. The results are found to be in reasonable agreement with the observed data procured from Bangladesh Inland Water Transport Authority.

     

Seasonal Variability of Thermohaline Front in the Central South China Sea

An upper layer thermohaline front across the South China Sea (SCS) basin from the South Vietnamese coast (around 15°N) to Luzon Island (around 19°N) has been identified using the Navy's open domain Generalized Digital Environmental Model (GDEM) monthly mean temperature and salinity data on a 0.5° × 0.5° grid. This front does not occur at the surface in summer. The strength of this front is around 1°C/100 km at the surface and 1.4°C/100 km at the subsurface (50 m deep). A cross-basin current, inverted using the P-vector method, is associating with the front. Meandering and eddies have been identified along this current. Seasonal and vertical variabilities of the thermohaline structure across this front are reported in this paper. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the seasonal variation of transport through the Tokara Strait

Results are described from a limited area barotropic model of the North Pacific with 1/3°×1/3° resolution and bounded by latitudes 10°N and 50°N and longitudes 120°E and 160°E. The model employs dynamics linearised about a state of rest and incorporates realistic bottom topography. It is driven by the seasonally varying part of the Hellerman and Rosenstein wind stress field and by inflow along the northern and eastern boundaries specified from a 1°×1° version of the model applied to the whole North Pacific. The model-calculated transport variations through the Tokara Strait are similar to those of the observed seasonal sea level differences across the Strait, although the model appears to underestimate the amplitude of the signal by a factor of at least 2. The inclusion of realistic bottom topography is shown to be crucial in determining the model response.     

Observed Variability of Thermohaline Fields,Currents and Eddies in the Western Bay of Bengal During BOBMEX-99

Field measurements during the Bay of Bengal Monsoon Experiment (BOBMEX-99), data from a deep sea moored buoy, and satellite altimeter were used to describe variability in the hydrographic and meso-scale features in the Bay of Bengal (BoB) during the summer monsoon of 1999. The thermohaline fields showed two regions of upsloping of isopleths centered at 82°E and 84.75°E, ～110 km and 450 km away from the coast, respectively, followed by downsloping. The upsloping/downsloping of isopleths and the alternating currents was part of cyclonic and anti-cyclonic circulation patterns in the western BoB. In this region, both wind and current were important in the dynamics of coastal upwelling. The observations showed a relationship between the propagating waves and eddy on variability of thermohaline fields. On an annual cycle, four Kelvin waves were observed in the BoB, but only the downwelling Kelvin wave formed during October entered the Arabian Sea. During the monsoon season, four eddies were formed in the western BoB, of which the anticyclonic eddy centered at 15°N, 84°E and the cyclonic eddy centered at 17.5°N, 84.5°E were prominent. The baroclinic instability caused by the opposing currents along the east coast and the wind stress curl favored the formation of eddies. Okhubo-Weiss and Isern-Fontanet parameter confirmed the presence of eddies in the BoB.     

Ocean temperature climate off North‐East New Zealand

The ocean temperature field off the north‐east coast of New Zealand is studied to quantify the annual cycle and reveal the intra‐ and inter‐annual variability. The data used are repeat expendable bathythermograph (XBT) sections between Auckland and either Suva or Honolulu which have been collected quarterly since 1986. These sections give temperature measurements between the surface and 800 m and Auckland and 30°S from 1986 to August 1999. The mean and annual cycle are compared with those from the NOAA World Ocean Atlas (WOA98). The results are similar; however WOA98 lacks the horizontal resolution to fully discern the East Auckland Current and North Cape Eddy, while the XBT analysis lacks the temporal resolution to discern higher frequency intra‐annual signals. The temperature variability in the mixed layer is dominated by the annual cycle, which accounts for 80–90% of the variance. The amplitude of the annual cycle diminishes rapidly with depth, from 2.8°C at the surface, to c. 0.1°C at 180 m. The phase of the annual cycle is retarded with depth, with peak temperatures occurring in February at the surface and in June/July at 180 m. Removing the annual cycle from the time series reveals the more subtle inter‐ and intra‐annual variability. This variability is of the order of 1°C in the upper 50 m, decreasing to 0.3°C at 400–500 m. The surface layer was cold between 1991 and 1994 (c. 0.7°C cooler than average), and 0.7°C warmer than average in 1999. The deeper ocean shows a different signal, being up to 0.3°C cooler in 1990–92, 0.3°C warmer in 1998, and c. 0.2°C warmer than average in 1999. The inter‐annual mixed layer variability is highly correlated with the Southern Oscillation Index and also with inter‐annual terrestrial air temperature and wind measurements from northern New Zealand. In contrast, at higher intra‐annual frequencies, the mixed layer variability is not correlated with air and wind measurements. At these higher frequencies, the air temperature is better correlated with the sea surface temperature (SST) than with the bulk mixed layer temperature.