Short-range prediction experiments of the Kuroshio path variabilities south of Japan

Assimilation and prediction experiments of the Kuroshio path variability south of Japan were conducted to investigate the predictability of the Kuroshio path. The assimilation and prediction system is composed of an eddy-resolving model and a three-dimensional variational analysis scheme with vertical coupled temperature–salinity empirical orthogonal function modes. The sea surface height (SSH) variability and the variations of the Kuroshio path of the assimilation fields are in good agreement with those observed. The results of the assimilation are then used as the initial conditions for 138 cases of 90-day prediction experiments conducted from 1993 to 2004. The predictive limit of our system is assessed by the SSH anomaly in the assimilation field and is found to be around 40–60 days, which is much longer than that of the persistence. The prediction results show good performance in the transition stage from a straight to a meandering path. For example, a large meandering event that occurred in August 2004 is successfully predicted in a 2-month forecast. Two types of failure cases are investigated. One is a case where the eastward propagation speed of the meander is faster than a real state. The dynamical response of the model to the assimilation revealed that an initial shock, caused by the dynamically unbalanced initial condition, induces the fast eastward propagation of the meander. The other case exhibits an unrealistic meander. In this case, a cold anomaly at an intermediate layer in the initial condition grows rapidly and results in the unrealistic meander. This implies that the Kuroshio path south of Japan has a chaotic nature. These facts revealed by the failure cases give us some insight for improving the predictive skill of the Kuroshio path variability.     

Open and coastal seas interactions south of Japan represented by an ensemble Kalman filter

We investigated the feasibility of the ensemble Kalman filter (EnKF) to reproduce oceanic conditions south of Japan. We have adopted the local ensemble transformation Kalman filter algorithm based on 20 members’ ensemble simulations of the parallelized Princeton Ocean Model (the Stony Brook Parallel Ocean Model) with horizontal resolution of 1/36°. By assimilating satellite sea surface height anomaly, satellite sea surface temperature, and in situ temperature and salinity profiles, we reproduced the Kuroshio variation south of Japan for the period from 8 to 28 February 2010. EnKF successfully reproduced the Kuroshio path positions and the water mass property of the Kuroshio waters as observed. It also detected the variation of the steep thermohaline front in the Kii Channel due to the intrusion of the Kuroshio water based on the observation, suggesting efficiency of EnKF for detection of open and coastal seas interactions with highly complicated spatiotemporal variability.     

Modeling fast ice formation and destruction in the Sea of Japan

The problems of parametric representation of the initial formation and the subsequent evolution of fast ice in freezing seas are discussed within the framework of the model of marine ice cover evolution. The mathematical form of this representation takes into account the processes of transformation of sea ice formations in open water areas in coastal regions into fast ice during its formation and inverse processes at the stage of its destruction. The parametric identification of the model was based on samples of long-term distributions of ice cover characteristics in the Sea of Japan, as well as the distributions of air temperature and wind speeds over sea water area. The model can be used to predict the state of fast ice in the ice cover of the Sea of Japan.     

Dynamical Control of Pacific Oceanic Low-frequency Variability on Western Boundary and Marginal Seas

The adjustment of sea surface height (SSH) around the coasts of the Japan/East Sea (JES) and the South China Sea (SCS) basins subjected to extratropical Pacific Oceanic low frequency variability is studied using a Kelvin-planetary wave model and a high resolution numerical model. It is found that the modulation of SSH around the coast of Japan is mainly determined by slow adjustment of planetary waves, which radiate from the west coast of Honshu and Hokkaido due to the coastal Kelvin wave. In contrast, the SSH modulation around the cost of the South China Sea basin is mainly determined by the coastal Kelvin wave, which transfers the anomalous SSH into the SCS via the Luzon Strait and out via the Mindoro Strait. The planetary waves radiating from the west coast of Palawan establish a nearly uniform SSH anomaly in the southern part of the SCS, bounded by an eastward jet at the latitude of the Mindoro Strait. Along the western boundary, SSH anomaly decreases almost linearly toward the south, in accordance with the changing local deformation radius. In these two marginal seas, the mean subtropical Pacific gyre circulation enhances SSH modulation induced by extratropical Pacific low frequency variability. Overall, the SSH adjustment in the JES and the SCS predicted by the analytical model agrees well with the numerical model simulation. Application of this model to interaction between these marginal seas and the open ocean is discussed.     

Development of a 2-km resolution ocean model covering the coastal seas around Japan for operational application

Ocean Dynamics - In order to expand the coastal ocean monitoring and forecasting system of the Japan Meteorological Agency from the Seto Inland Sea to the entire coastal seas of Japan, a 2-km...     

The predictability of near-coastal currents using a baroclinic unstructured grid model

A limited domain, coastal ocean forecast system consisting of an unstructured grid model, a meteorological model, a regional ocean model, and a global tidal database is designed to be globally relocatable. For such a system to be viable, the predictability of coastal currents must be well understood with error sources clearly identified. To this end, the coastal forecast system is applied at the mouth of Chesapeake Bay in response to a Navy exercise. Two-day forecasts are produced for a 10-day period from 4 to 14 June 2010 and compared to real-time observations. Interplay between the temporal frequency of the regional model boundary forcing and the application of external tides to the coastal model impacts the tidal characteristics of the coastal current, even contributing a small phase error. Frequencies of at least 3 h are needed to resolve the tidal signal within the regional model; otherwise, externally applied tides from a database are needed to capture the tidal variability. Spatial resolution of the regional model (3 vs 1 km) does not impact skill of the current prediction. Tidal response of the system indicates excellent representation of the dominant M ₂ tide for water level and currents. Diurnal tides, especially K ₁, are amplified unrealistically with the application of coarse 27-km winds. Higher-resolution winds reduce current forecast error with the exception of wind originating from the SSW, SSE, and E. These winds run shore parallel and are subject to strong interaction with the shoreline that is poorly represented even by the 3-km wind fields. The vertical distribution of currents is also well predicted by the coastal model. Spatial and temporal resolution of the wind forcing including areas close to the shoreline is the most critical component for accurate current forecasts. Additionally, it is demonstrated that wind resolution plays a large role in establishing realistic thermal and density structures in upwelling prone regions.     

雷达卫星SAR与防卫气象卫星SSM/I对渤海海冰的观测研究

用雷达卫星(RADARSAT)合成孔径雷达(SAR)和防卫气象卫星(DMSP)特别微波成像辐射计(SSM/I)在1999年1月23日同一天对中国渤海区域海冰的观测数据，进行了主动SAR与被动SSM/I的组合研究.用一层海冰电磁散射辐射的建模与数值模拟，分析了中纬度渤海沿岸海冰的特征性变化，提出用SSM/I辐射亮度温度的散射指数、极化指数、极化比率来识别渤海海冰；用RADARSAT高分辨率SAR数据(水平极化后向散射系数)对渤海海冰物理特征进行识别与分类.这些特征指数在时间与空间尺度上的相关变化可有效地应用于渤海海冰的监察.     

Predictability of storm wave heights in the ice-free Beaufort Sea

A predictability study on wave forecast of the Arctic Ocean is necessary to help identify hazardous areas and ensure sustainable shipping along the trans-Arctic routes. To assist with validation of the Arctic Ocean wave model, two drifting wave buoys were deployed off Point Barrow, Alaska for two months in September 2016. Both buoys measured significant wave heights exceeding 4 m during two different storm events on 19 September and 22 October. The NOAA-WAVEWATCH III^? model with 16-km resolution was forced using wind and sea ice reanalysis data and obtained general agreement with the observation. The September storm was reproduced well; however, model accuracy deteriorated in October with a negative wave height bias of around 1 m during the October storm. Utilising reanalysis data, including the most up-to-date ERA5, this study investigated the cause: grid resolution, wind and ice forcing, and in situ sea level pressure observations assimilated for reanalysis. The analysis has found that there is a 20% reduction of in situ SLP observations in the area of interest, presumably due to fewer ships and deployment options during the sea ice advance period. The 63-member atmospheric ensemble reanalysis, ALERA2, has shown that this led to a larger ensemble spread in the October monthly mean wind field compared to September. Since atmospheric physics is complex during sea ice advance, it is speculated that the elevated uncertainty of synoptic-scale wind caused the negative wave model bias. This has implications for wave hindcasts and forecasts in the Arctic Ocean.     

Shelf-basin exchange in the Laptev Sea in the warming climate: a model study

GCM-based forecast simulations predict continuously increasing seasonality of the sea ice cover and an almost ice-free, summer-time, Arctic Ocean within several decades from the present. In this study we use a primitive equation ocean model: NEMO, coupled with the sea ice model LIM2, to test the hypothesis that under such an increased range in seasonal ice cover the intensity of shelf-basin water exchange will significantly increase. We use the simulated results for the Laptev Sea from a global model run 1958–2007 and compare results for two years with anomalously high and low summer sea ice extents: 1986–1987 and 2006–2007. The shelf–basin fluxes of volume, heat and salt during specific seasons are evaluated and attributed to plausible driving processes, with particular attention to dense water cascading. Analyses of the model temperature distribution at the depth of the intermediate maximum, associated with Atlantic Water, have shown a marked increase of the amount of the local origin cold water in late winter 2007 in the region, where dense water typically appears as a result of its formation on the shelf and subsequent downslope leakage. Calculation of the shelf-basin exchange during March-May in both years confirmed a substantial increase (a factor of two) of fluxes in “ice-free” 2007 compared to the “icy” 1987. According to several past model studies, dense water production on Arctic shelves in winter driven by ice freezing and brine rejection is not likely to cease in a warmer climate, but rather to increase. There is also observational evidence that cascading in the seasonally ice covered seas (e.g. the Barents Sea) is much more efficient than it is in the permanently ice covered Arctic Ocean, which supports these model results.     

秋季北极海冰与欧亚冬季气温在年代际和年际尺度上的不同联系

北极海冰的急剧消融在近年来欧亚大陆频发的低温事件中扮演着关键角色.秋季北极海冰的偏少对应着冬季欧亚大陆的低温天气,然而二者的联系在年代际和年际两种时间尺度上存在显著区别.本文运用1979—2012年哈德莱中心第一套海冰覆盖率(HadISST1)、欧洲中心(ERA_Interim)的2m温度、风场、海平面气压场、高度场等资料,分别研究了年代际和年际时间尺度上前期秋季北极海冰与欧亚冬季气温的联系.结果表明,欧亚和北极地区(0°—160°E,15°N—90°N)的冬季气温具有显著的年代际和年际变化.在年代际尺度上,温度异常分布在21世纪初由北极冷-大陆暖转为北极暖-大陆冷.这一年代际转折与前期秋季整个北极地区的海冰年代际减少联系密切.秋季北极全区海冰年代际偏少对应冬季欧亚大陆中高纬地区的高压异常,有利于北大西洋的暖湿气流北上和北极的冷空气南侵,造成北极暖-大陆冷的温度分布;在年际时间尺度上,温度异常分布主要由第一模态的年际变化部分和第二模态组成,且第一模态包含的年际变率信号也存在显著的年代际变化.年际尺度上全区北极海冰对欧亚冬季气温的影响远不及位于北冰洋西南边缘的巴伦支海、喀拉海和拉普捷夫海西部(30°E—120°E,75°N—85°N)的关键区海冰影响显著.关键区内海冰的偏少会引发冬季的北大西洋涛动负位相,导致北大西洋吹往欧亚大陆的暖湿气流减弱和欧亚大陆中高纬地区的气温偏低.     

Dissolved barium as a tracer of Kuroshio incursion in the Kuroshio region east of Taiwan Island and the adjacent East China Sea

From May to June 2014, the geochemical characteristics of dissolved barium(Ba) in sea water and its influx from the Kuroshio into the East China Sea(ECS) were studied by investigation of the Kuroshio mainstream east of Taiwan Island and the adjacent ECS. This allowed for the scope and extent of the Kuroshio incursion to be quantitatively described for the first time by using Ba as a tracer. The concentration of Ba in the Kuroshio mainstream increased gradually downward from the surface in the range 4.91–19.2 μg L.1. In the surface layer of the ECS, the Ba concentration was highest in coastal water and gradually decreased seaward, while it was higher in coastal and offshore water but lowest in middle shelf for bottom layer. The influx of Ba from Kuroshio into the ECS during May to October was calculated to be 2.19×108 kg by a water exchange model, in which the subsurface layer had the largest portion. The distribution of Ba indicated that Kuroshio upwelled in the sea area northeast of Taiwan Island. The north-flowing water in the Taiwan Strait restrained the incursion of Kuroshio surface water onto the ECS shelf, while Kuroshio subsurface water gradually affected the bottom of the ECS from outside. The results of end member calculation, using Ba as a parameter, showed that the Kuroshio surface water had little impact on the ECS, while the Kuroshio subsurface water formed an intrusion current by climbing northwest along the bottom of the middle shelf from the sea area northeast of Taiwan Island into the Qiantang Estuary, of which the volume of Kuroshio water was nearly 65%. Kuroshio water was the predominant part of the water on the outer shelf bottom and its proportion in areas deeper than the 100 m isobath could reach more than 95%. In the DH9 section(north of Taiwan Island), Kuroshio subsurface water intruded westward along the bottom from the shelf edge and then rose upward(in lower proportion). Kuroshio water accounted for 95% of the ocean volume could reach as far as 122°E. Ba was able to provide detailed tracing of the Kuroshio incursion into the ECS owing to its geochemical characteristics and became an effective tracer for revealing quantitative interactions between the Kuroshio and the ECS.     

Numerical simulation of ice-ocean variability in the Barents Sea region

A dynamic–thermodynamic sea ice model has been coupled to a three-dimensional ocean general circulation model for the purpose of conducting ocean climate dynamical downscaling experiments for the Barents Sea region. To assess model performance and suitability for such an application, the coupled model has been used to conduct a hindcast for the period 1990–2002. A comparison with available observations shows that the model successfully tracks seasonal and inter-annual variability in the ocean temperature field and that the simulated horizontal and vertical distribution of temperature are in good agreement with observations. The model results follow the seasonal and inter-annual variability in sea ice cover in the region, with the exception that the model results show too much ice melting in the northern Barents Sea during summer. The spatial distribution of the winter simulated sea ice cover is in close agreement with observations. Modelled temperatures and ice concentrations in the central Barents Sea are biased too high and too low, respectively. The probable cause is too high inflow of Atlantic Water into the Barents. The seasonal and inter-annual fluctuations in temperature and sea ice cover in the central Barents are, however, in excellent agreement with observations. Salt release during the freezing process in the numerical simulation exhibits considerable inter-annual variability and tends to vary in an opposite manner to the net inflow volume flux at the western entrance of the Barents Sea. Overall, the model produces realistic ice-ocean seasonal and inter-annual variability and should prove to be a useful tool for dynamical downscaling applications.     

Ensemble numerical forecasts of the sporadic Kuroshio water intrusion (kyucho) into shelf and coastal waters

The finite volume coastal ocean model downscaling ocean reanalysis and forecast data provided by the Japan Coastal Ocean Predictability Experiment (JCOPE2) are used to forecast sudden Kuroshio water intrusion events (kyucho) induced by frontal waves amplified south of the Bungo Channel in 2010. Two-month hindcast computations give initial conditions of the following 3-month forecasts computations which consist of ten ensemble members. The temperature time series computed by these ten members are averaged to compare with that actually observed in the Bungo Channel, where sudden temperature rises related to kyucho events are remarkable in February, August, and September. Overall, the intense kyucho events actually observed in these months are predicted successfully. However, intense kyucho events are forecasted frequently during the period of May through June even though intense kyucho events are absent during this period in the actual ocean. It is suggested that the present downscaling forecast model requires reliable lateral boundary conditions provided by JCOPE2 data to which numerous Argo data are assimilated to enhance the accuracy. In addition, it seems likely that the model accuracy is reduced by small eddies moving along the shelf break.     

Investigation of three-dimensional circulation and hydrography over the Pearl River Estuary of China using a nested-grid coastal circulation model

The Pearl River Estuary (PRE) in South China's Guangdong Province is a subtropical estuary with highly irregular topography and dynamically complicated circulations. A nested-grid coastal circulation modelling system is used in this study to examine dynamic responses of the PRE to tides, meteorological forcing and buoyancy forcing. The nested-grid modelling system is based on the Princeton Ocean Model and consists of three downscaling subcomponents: including an outer-most model with a coarse horizontal resolution of ~7 km for simulating tidally forced and wind-driven surface elevations and depth-mean currents over the China Seas from Bohai Sea to the northern South China Sea and an innermost model with a fine resolution of ~1.2 km for simulating the 3D coastal circulation and hydrography over the PRE and adjacent coastal waters. Model results during the winter northeast monsoon surge in January and super typhoon Koryn in June of 1993 are used to demonstrate that the 3D coastal circulation and hydrographic distributions in the PRE are affected by tides, winds and buoyancy forcing associated with river discharge from the Pearl River with significant seasonal and synoptic variabilities.     

Seasonal variation of upper layer circulation in the northern part of the East/Japan Sea

Seasonal variation of upper layer circulation in the northern part of the East/Japan Sea and its mechanism were investigated using empirical orthogonal function (EOF) analysis with satellite sea surface heights over the northern East/Japan Sea and a three-dimensional circulation model. The spatial structure and temporal variation of first EOF mode, which explains about 64% of the total variance, indicate that a large cyclonic circulation in the northern East/Japan Sea shows a semi-annual variation with maximum strength in summer and winter. According to numerical model result, the Liman Cold Current, accepted as a major current in the northern East/Japan Sea, is well mixed vertically by the winter monsoon and the current in the upper layer has a relatively deep structure, with a maximum westward speed of about 20 cm/s in winter. On the other hand, in summer the current has a stronger baroclinic structure of velocity than in winter. Numerical experiments showed that in summer the temporal variation of upper layer circulation is controlled by thermal forcing, such as sea surface heat flux and inflow of heat transport into the East/Japan Sea through the Korea/Tsushima Strait. Moreover, the cyclonic circulation in the upper layer of the northern East/Japan Sea is also generated and strengthened by the positive wind stress curl occupying most of the East/Japan Sea during the winter. The seasonal variation of each forcing that drives the circulation is responsible for the strength or weakness of the upper layer circulation in the northern East/Japan Sea. The contribution of each forcing to the seasonal variation of the upper layer circulation is examined through sensitivity experiments. According to these numerical experiments, the upper layer circulation in the northern East/Japan Sea is strengthened twice a year, in winter and summer, and this semi-annual variation is determined by a combination of wind (winter) and thermal (summer) forcing.     

Korean East Coast wave predictions by means of ensemble Kalman filter data assimilation

To respond to the need for preventing offshore and coastal accidents, damage and flooding, a state-of-the-art coastal wave forecast system for the East Coast of Korea waters is being developed. Given that the quality of the input wind has been identified as the main factor influencing the quality of the wave results, the effectiveness of adjusting the wind fields by means of data assimilation using the ensemble Kalman filter technique has been explored. In this article the model setup, the data assimilation parameters and the validation of the predictions during stormy periods is described. The validation shows that the model is able to provide predictions of coastal waves fulfilling available benchmarks; especially, the data assimilation analysis and forecast predictions are judged to be of high quality.     

An altimetric transport index for Kuroshio inflow northeast of Taiwan Island

The Kuroshio inflow northeast of Taiwan Island plays an important role in the heat and nutrient balances over the East China Sea(ECS). Based on merged satellite altimeter data and the PCM-1 mooring observation at the East Taiwan Channel(ETC), the study employs a correlation iteration scheme to find the optimal transport index for the Kuroshio inflow. The sea level difference with the highest correlation to the ETC transport is across the ECS shelf break rather than along the PCM-1 line. The counter-intuitive result is caused by large signal noise and poor track coverage of altimeters near the Taiwan coast. The optimal altimetric index is highly correlated with the two-year in-situ measurements as well as the ten-year output of the global assimilation model. It serves as a better estimator of Kuroshio inflow than those using tidal gauge data, and helps pinpoint a 5 cm mismatch of mean sea level in the Keelung tidal record. The mean transport of Kuroshio inflow based on the twenty-year altimetric index is 20.55 Sv with a standard deviation of 3.05 Sv. Wavelet spectrum of the index reveals that semi-annual period dominates the Kuroshio variation northeast of Taiwan Island.     

Argo data assimilation in ocean general circulation model of Northwest Pacific Ocean

The Argo temperature and salinity profiles in 2005–2009 are assimilated into a coastal ocean general circulation model of the Northwest Pacific Ocean using the ensemble adjustment Kalman filter (EAKF). Three numerical tests, including the control run (CTL) (without data assimilation, which serves as the reference experiment), ensemble free run (EnFR) (without data assimilation), and EAKF experiment (with Argo data assimilation using EAKF), are carried out to examine the performance of this system. Using the restarts of different years as the initial conditions of the ensemble integrations, the ensemble spreads from EnFR and EAKF are all kept at a finite value after a sharp decreasing in the first few months because of the sensitive of the model to the initial conditions, and the reducing of the ensemble spread due to Argo data assimilation is not much. The ensemble samples obtained in this way can well represent the probabilities of the real ocean states, and no ensemble inflation is necessary for this EAKF experiment. Different experiment results are compared with satellite sea surface temperature (SST) data and the Global Temperature-Salinity Profile Program (GTSPP) data. The comparison of SST shows that modeled SST errors are reduced after data assimilation; the error reduction percentage after assimilating the Argo profiles is about 10?% on average. The comparison against the GTSPP profiles, which are independent of the Argo profiles, shows improvements in both temperature and salinity. The comparison results indicated a great error reduction in all vertical layers relative to CTL and the ensemble mean of EnFR; the maximum value for temperature and salinity reaches to 85?% and 80?%, respectively. The standard deviations of sea surface height are employed to examine the simulation ability, and it is shown that the mesoscale variability is improved after Argo data assimilation, especially in the Kuroshio extension area and along the section of 10°N. All these results suggest that this system is potentially useful for improving the simulation ability of oceanic numerical models.     

Structure and dynamics of the sudden acceleration of Kuroshio off Cape Shionomisaki

A sudden acceleration of the Kuroshio jet appears off Cape Shionomisaki in the high-resolution (horizontal resolution of 1/36°) JCOPE 2 ocean reanalysis data. Using this dataset, we investigated the structure of the Kuroshio acceleration. The increase in the velocity of the current is accompanied by a downstream flow separation from the coast and an outcrop of cold temperature inshore. The acceleration of Kuroshio appears when it takes a near-shore path. Cape Shionomisaki amplifies the responses to the Kuroshio flow by creating the zonal velocity acceleration toward the downstream region when the Kuroshio flows closer to the coast. The Kuroshio acceleration coincided with the topographic ridge on the continental shelf near Cape Shionomisaki. This relation suggests that the dynamics of the acceleration is linked to the topographic feature. We proposed an explanation of the Kuroshio acceleration using a hydraulic control theory. An analytical solution was applied to the coastal topography around the Kii Peninsula. The solution captured some aspects of the Kuroshio acceleration.     

Effects of deforestation and weather on diurnal frost heave processes on the steep mountain slopes in south central Japan

Freezing and thawing processes play an important role for the gravitational transport of surface materials on steep mountain slopes in Japan. The effects of deforestation on frost heave activity were observed through the 2012/2013 winter season in Ikawa University Forest, a southern mountainous area in central Japan (1180–1310 m above sea level). During periods without snow cover, needle ice development prevailed at a clear‐cut site, and the downslope sediment movement of upper soil was 10 to 15 cm through the winter season. At a non‐cut site, rise and fall in the ground surface level prevailed on a weekly scale, with no evident downslope movements at the surface; ice lens formation in the soil layer is assumed. Abrupt changes in the radiation budget, such as the strengthening of nighttime radiative cooling and increases in daytime direct insolation, induced frequent development/deformation of needle ice at the clear‐cut site. In snow‐free periods, the day‐to‐day variability in needle ice growth length and in nighttime averaged net radiation showed significant correlations; cloudy weather with warmer and moist air intrusion associated with synoptic disturbances prevented the occurrence of needle ice. Namely, day‐to‐day weather changes directly affected the mass movement of the upper soil after deforestation. Shallow snow cover occurred discontinuously through the winter and is likely an important factor in keeping the soil moisture sufficiently high in the upper soil layer for initiating needle ice during snow‐free periods. We also discuss contributions of coastal extratropical cyclone activities providing both snow cover and cloudy weather in the southern mountain areas of central Japan to the intra‐seasonal variability in frost heave and its indirect effect on soil creep and landslides on the deforested steep slopes. Copyright © 2015 John Wiley & Sons, Ltd.