The study of ocean circulation dynamics with numerical models

Summarized is the author's study of the ocean circulation dynamics with numerical models, for which he was honored by the Okada Prize (1979). Cited topics are formation of the western boundary current and its variation associated with imposed wind stress, some effects of a marine ridge on the boundary current, coastal upwelling circulation and coastal thermohaline front formation. Recent modelling efforts in Japan, specifically on numerical study of ocean circulation dynamics are also reviewed.     

Trends of the Galician upwelling in the context of climate change

Coastal upwelling is a phenomenon of great importance both for the study of ocean dynamics and for the development of fish production in some coastal regions. Our study region, the Galician coast, lies at the northern end of the Canary–Iberian Peninsula upwelling system. Knowing the changes provoked by climate change on this upwelling system is particularly relevant for the future of this area taking into account the social and economic importance of fishing activities in this region. In this paper we study the trends in the intensity and frequency of upwelling in the Galician coast and the expected changes in this phenomenon for the next decades using three regional models implemented within the European project ENSEMBLES. As a main result, we observe that the models show a positive trend in both the intensity and frequency of upwelling phenomenon for the future, particularly significant in spring and summer which are the seasons favorable for upwelling. In autumn and winter there are no significant changes.     

地幔柱与洋脊分段成因探讨

地幔柱产生岩浆的核幔交界部位相当于"岩浆洋",下地幔与上地幔、上地幔与软流圈之间的岩浆流通道相当于"岩浆河流",岩浆由软流圈经岩石圈喷出地表的通道相当于"岩浆小溪",在它们之间都会有许多的"岩浆支流",其中有的"支流"将这些"河流"与"河流"或"河流"与"小溪"联系起来.洋脊分段现象是由于地幔柱中存在横向的"岩浆支流"...     

Observational and numerical modeling methods for quantifying coastal ocean turbulence and mixing

In this review paper, state-of-the-art observational and numerical modeling methods for small scale turbulence and mixing with applications to coastal oceans are presented in one context. Unresolved dynamics and remaining problems of field observations and numerical simulations are reviewed on the basis of the approach that modern process-oriented studies should be based on both observations and models. First of all, the basic dynamics of surface and bottom boundary layers as well as intermediate stratified regimes including the interaction of turbulence and internal waves are briefly discussed. Then, an overview is given on just established or recently emerging mechanical, acoustic and optical observational techniques. Microstructure shear probes although developed already in the 1970s have only recently become reliable commercial products. Specifically under surface waves turbulence measurements are difficult due to the necessary decomposition of waves and turbulence. The methods to apply Acoustic Doppler Current Profilers (ADCPs) for estimations of Reynolds stresses, turbulence kinetic energy and dissipation rates are under further development. Finally, applications of well-established turbulence resolving particle image velocimetry (PIV) to the dynamics of the bottom boundary layer are presented. As counterpart to the field methods the state-of-the-art in numerical modeling in coastal seas is presented. This includes the application of the Large Eddy Simulation (LES) method to shallow water Langmuir Circulation (LC) and to stratified flow over a topographic obstacle. Furthermore, statistical turbulence closure methods as well as empirical turbulence parameterizations and their applicability to coastal ocean turbulence and mixing are discussed. Specific problems related to the combined wave-current bottom boundary layer are discussed. Finally, two coastal modeling sensitivity studies are presented as applications, a two-dimensional study of upwelling and downwelling and a three-dimensional study for a marginal sea scenario (Baltic Sea). It is concluded that the discussed methods need further refinements specifically to account for the complex dynamics associated with the presence of surface and internal waves.     

The Long-term Ecosystem Observatory: an integrated coastalobservatory

An integrated ocean observatory has been developed and operated in the coastal waters off the central coast of New Jersey, USA. One major goal for the Long-term Ecosystem Observatory (LEO) is to develop a real-time capability for rapid environmental assessment and physical/biological forecasting in coastal waters. To this end, observational data are collected from satellites, aircrafts, ships, fixed/relocatable moorings and autonomous underwater vehicles. The majority of the data are available in real-time allowing for adaptive sampling of episodic events and are assimilated into ocean forecast models. In this observationally rich environment, model forecast errors are dominated by uncertainties in the model physics or future boundary conditions rather than initial conditions. Therefore, ensemble forecasts with differing model parameterizations provide a unique opportunity for model refinement and validation. The system has been operated during three annual coastal predictive skill experiments from 1998 through 2000. To illustrate the capabilities of the system, case studies on coastal upwelling and small-scale biological slicks are discussed. This observatory is one part of the expanding network of ocean observatories that will form the basis of a national observation network     

Dynamic of ENSO towards upwelling and thermal front zone in the east coast of Peninsular Malaysia

The El Ni?o Southern Oscillation(ENSO) is a natural phenomenon that relates to the fluctuation of temperatures over the Pacific Ocean. The ENSO significantly affects the ocean dynamics including upwelling event and coastal front. A recent study discovered the seasonal upwelling in the east coast of Peninsular Malaysia(ECPM), which is significant to the fishery industry in this region. Thus, it is vital to have a better understanding of the influence of ENSO towards the coastal upwelling and thermal front in the ECPM. The sea surface temperature(SST) data achieved from moderate resolution imaging spectroradiometer(MODIS) aboard Aqua satellite are used in this study to observe the SST changes from 2005 to 2015. However, due to cloud cover issue, a reconstruction of data set is applied to MODIS data using the data interpolating empirical orthogonal function(DINEOF) to fill in the missing gap in the dataset based on spatial and temporal available data. Besides, a wavelet transformation analysis is done to determine the temperature fluctuation throughout the time series. The DINEOF results show the coastal upwelling in the ECPM develops in July and reaches its peak in August with a clear cold water patch off the coast. There is also a significant change of SST distribution during the El Ni?o years which weaken the coastal upwelling event along the ECPM. The wavelet transformation analysis shows the highest temperature fluctuation is in 2009–2010 which indicates the strongest El Ni?o throughout the time period. It is suggested that the El Ni?o is favourable for the stratification in water column thus it is weakening the upwelling and thermal frontal zone formation in ECPM waters.     

Modeling the seasonal variability of marine ecosystem in the region of the Central-Eastern Atlantic

A 3D eco-hydrodynamical model of high resolution (0.25° × 0.25°, 27 σ-levels) is used to simulate the seasonal variability of the ocean circulation and marine ecosystem in the Central-Eastern Basin of the North Atlantic including the Canary upwelling system. According to the model results, in the winter period, the “patches” of maximal phytoplankton and zooplankton biomass are often located in upwelling zones in the open ocean on the periphery of cyclonic eddies rather than in the coastal upwelling zones. In the summer period, when the phytoplankton biomass reaches maximal (in the annual cycle) values, the maxima of the phytoplankton are located in the coastal upwelling zones. As shown, there is no simple relationship between the nitrate distributions, on the one hand, and the phytoplankton and zooplankton ones, on the other hand.     

Patterns and processes in the California Current System

The California Current System (CCS) is forced by the distribution of atmospheric pressure and associated winds in relation to the west coast of North America. In this paper, we begin with a simplified case of winds and a linear coast, then consider variability characteristic of the CCS, and conclude by considering future change. The CCS extends from the North Pacific Current (～50°N) to off Baja California, Mexico (～15–25°N) with a major discontinuity at Point Conception (34.5°N). Variation in atmospheric pressure affects winds and thus upwelling. Coastal, wind-driven upwelling results in nutrification and biological production and a southward coastal jet. Offshore, curl-driven upwelling results in a spatially large, productive habitat. The California Current flows equatorward and derives from the North Pacific Current and the coastal jet. Dominant modes of spatial and temporal variability in physical processes and biological responses are discussed. High surface production results in deep and bottom waters depleted in oxygen and enriched in carbon dioxide. Fishing has depleted demersal stocks more than pelagic stocks, and marine mammals, including whales, are recovering. Krill, squid, and micronekton are poorly known and merit study. Future climate change will differ from past change and thus prediction of the CCS requires an understanding of its dynamics. Of particular concern are changes in winds, stratification, and ocean chemistry.     

2006年夏季琼东、粤西沿岸上升流研究

利用2006年夏季广东、海南、广西近海的海洋水文调查资料和卫星遥感QuikSCAT风场资料分析琼东、粤西沿岸上升流的空间结构特征, 探讨风场、风应力旋度对上升流的影响以及上升流区水温、海流、海平面对上升流的响应。结果表明:琼东、粤西沿岸上升流区并非相互独立, 从10 m层以下已经连成一片。琼东沿岸上升流主要由夏季西南季风驱动而产生, 风应力旋度也有一定贡献。琼东沿岸上升流的强度比粤西强。琼东沿岸海域的上层海水(18 m以浅)以离岸运动为主, 中下层海水以向岸运动为主。上层的离岸流速大于中下层的向岸流速。琼东沿岸的上升流现象是间歇性的, 与沿岸风速强弱有关。琼东沿岸海域海平面的升降与上升流的强弱有良好的关系, 上升流的强弱滞后于海平面的升降约1～2 d。     

The physical structure of a cold filament in a Chilean upwelling zone (Península de Mejillones, Chile, 23°S)

Cold filaments associated with Eastern Boundary Currents are typically narrower than 100 km but can be several hundred kilometers long, extending from the coast to the open ocean in upwelling areas. One such structure, observed off Península de Mejillones (23°S, Chile), was studied with both satellite images and two 5-days hydrographic cruises carried out during January 1997. The study used a coastal grid of 31 stations in an area of 165 ×155 km², approximately. The spatial distribution of the filament and its change between cruises are described from the horizontal distributions of dynamic height, temperature, salinity and dissolved oxygen. The filament was a shallow feature (thickness <100 m) and extended at least 165 km toward the open ocean. A meandering northward current flowed at the borders of the filament, separating oceanic and coastal waters of different physical properties. Comparisons of cross sections of the filament near the coast and in the oceanic zone show the ascent of the shallow salinity minimum (SSM), and its extension toward the ocean, bound to the filament. It is concluded that Subantarctic Water ((SAAW) distinguish by low salinity, high dissolved oxygen) and Equatorial Subsurface Water ((ESSW) high salinity, low dissolved oxygen, high nutrient content) form this filament, and that their relative proportions depend on the strength of the coastal upwelling. Thus, the knowledge of the dynamics of these structures is fundamental to better understanding of the spatial distribution of important biological variables, such as nutrients and chlorophyll, in the coastal ecosystem.     

Coastal upwelling in the California current system

Coastal upwelling in the California Current system has been the subject of large scale studies off California and Baja California, and of small scale studies off Oregon. Recent studies of the winds along the entire coast from 25°N to 50°N indicate that there are significant along-shore variations in the strength of coastal upwelling, which are reflected in the observed temperature distribution. Active upwelling appears to be restricted to a narrow coastal band (about 10–25 km wide) along the entire coast, but the region influenced by coastal upwelling may be much wider. Intensive observations of the upwelling zone during summer off Oregon show the presence of a southward coastal jet at the surface, a mean vertical shear, a poleward undercurrent along the bottom, and persistently sloping isopycnals over the continental shelf; most of the upwelling there occurs during relatively short periods (several days long) of upwelling-favorable winds. During the upwelling season off Oregon, the offshore Ekman transport is carried by the surface Ekman layer, and the onshore return flow occurs through a quasi-geostrophic interior. It is not known whether the structure and dynamics observed off Oregon are typical of the upwelling zone along the entire coast, though some of the same features have been observed off Baja California. Current and future research will eventually show whether the Oregon results are also applicable in the region of persistently strong upwelling-favorable winds off northern California, and in the region of complex bathymetry off central and southern California.     

Sea Level Variations Associated with Upwelling/ Downwelling Along the West Coast of India

A three-dimensional numerical model is developed and used to study the coastal upwelling processes and corresponding seasonal changes in the sea level along the west coast of India. The upwelling and associated sea level variations are seen as a response of coastal ocean to pure wind stress forcing. The model is designed to represent coastal ocean physics by resolving surface and bottom Ekman layers as realistically as possible. The prognostic variables are the three components of the velocity field, temperature, salinity and turbulent energy. The governing equations together with their boundary conditions are solved by finite-difference techniques. Experiments are performed to investigate sea level fluctuations associated with the thermal response and alongshore currents of the coastal waters. The model is forced with mean monthly wind stress forcing of January, May, July and September representing northeast monsoon and different phases of the southwest monsoon. It is known from the observational study that the upwelling process reaches to the surface waters by May along the coastal waters of the extreme southwest peninsular region. The process is more intense in July compared to May and September and its strength decreases from south to north. However, during the northeast monsoon season, which is represented by January wind stress forcing in the model, downwelling is simulated along the coast. The model simulations of the coastal response are compared with the observations and are found to be in good agreement. The maximum computed vertical velocity of about 2.0 ×10 -3 cm s -1 is predicted in July in the southern region off the coast.     

Sea Level Variations Associated with Upwelling/ Downwelling Along the West Coast of India

A three-dimensional numerical model is developed and used to study the coastal upwelling processes and corresponding seasonal changes in the sea level along the west coast of India. The upwelling and associated sea level variations are seen as a response of coastal ocean to pure wind stress forcing. The model is designed to represent coastal ocean physics by resolving surface and bottom Ekman layers as realistically as possible. The prognostic variables are the three components of the velocity field, temperature, salinity and turbulent energy. The governing equations together with their boundary conditions are solved by finite-difference techniques. Experiments are performed to investigate sea level fluctuations associated with the thermal response and alongshore currents of the coastal waters. The model is forced with mean monthly wind stress forcing of January, May, July and September representing northeast monsoon and different phases of the southwest monsoon. It is known from the observational study that the upwelling process reaches to the surface waters by May along the coastal waters of the extreme southwest peninsular region. The process is more intense in July compared to May and September and its strength decreases from south to north. However, during the northeast monsoon season, which is represented by January wind stress forcing in the model, downwelling is simulated along the coast. The model simulations of the coastal response are compared with the observations and are found to be in good agreement. The maximum computed vertical velocity of about 2.0 2 10 -3 cm s -1 is predicted in July in the southern region off the coast.     

Impact of Estuarine Outflow on Coastal Upwelling

A two-dimensional numerical model is developed facilitating the locationwise study of coastal upwelling. The coastal rigid boundary in the model is replaced with an open boundary to understand the dynamical response of the coastal ocean in the presence of an estuary. The model is applied to the east coast of India in a plane perpendicular to the coast of Kakinada where the Godavari river joins the Bay of Bengal. The model is driven, starting from a state of rest, by the combined effect of the wind stress forcing and the freshwater discharge from the estuary. Two numerical experiments were conducted to study the effect of the variation in the freshwater discharge on upwelling. It is found that the freshwater discharge from the Godavari estuary suppresses the upwelling off Kakinada.     

Impact of Estuarine Outflow on Coastal Upwelling

A two-dimensional numerical model is developed facilitating the locationwise study of coastal upwelling. The coastal rigid boundary in the model is replaced with an open boundary to understand the dynamical response of the coastal ocean in the presence of an estuary. The model is applied to the east coast of India in a plane perpendicular to the coast of Kakinada where the Godavari river joins the Bay of Bengal. The model is driven, starting from a state of rest, by the combined effect of the wind stress forcing and the freshwater discharge from the estuary. Two numerical experiments were conducted to study the effect of the variation in the freshwater discharge on upwelling. It is found that the freshwater discharge from the Godavari estuary suppresses the upwelling off Kakinada.     

Responses of nutrient biogeochemistry and nitrogen cycle to seasonal upwelling in coastal waters of the eastern Hainan Island

The coastal upwelling has profound influence on the surrounding ecosystem by supplying the nutrient-replete water to the euphotic zone. Nutrient biogeochemistry was investigated in coastal waters of the eastern Hainan Island in summer 2015 and autumn 2016. From perspectives of nutrient dynamics and physical transport, the nutrient fluxes entered the upper 50 m water depth（between the mixed layer and the euphotic zone） arisen from the upwelling were estimated to be 2.5-5.4 mmol/（m²·d）,...     

The physical structure of a cold filament in a Chilean upwelling zone (Penı&#x0301;nsula de Mejillones,Chile, 23°S)

Cold filaments associated with Eastern Boundary Currents are typically narrower than 100 km but can be several hundred kilometers long, extending from the coast to the open ocean in upwelling areas. One such structure, observed off Penı&#x0301;nsula de Mejillones (23°S, Chile), was studied with both satellite images and two 5-days hydrographic cruises carried out during January 1997. The study used a coastal grid of 31 stations in an area of 165 ×155 km², approximately. The spatial distribution of the filament and its change between cruises are described from the horizontal distributions of dynamic height, temperature, salinity and dissolved oxygen. The filament was a shallow feature (thickness <100 m) and extended at least 165 km toward the open ocean. A meandering northward current flowed at the borders of the filament, separating oceanic and coastal waters of different physical properties. Comparisons of cross sections of the filament near the coast and in the oceanic zone show the ascent of the shallow salinity minimum (SSM), and its extension toward the ocean, bound to the filament. It is concluded that Subantarctic Water ((SAAW) distinguish by low salinity, high dissolved oxygen) and Equatorial Subsurface Water ((ESSW) high salinity, low dissolved oxygen, high nutrient content) form this filament, and that their relative proportions depend on the strength of the coastal upwelling. Thus, the knowledge of the dynamics of these structures is fundamental to better understanding of the spatial distribution of important biological variables, such as nutrients and chlorophyll, in the coastal ecosystem.     

Influence of bottom topography on an upwelling current: Generation of long trapped filaments

We investigate the influence of bottom topography on the formation and trapping of long upwelling filaments using a 2-layer shallow water model on the f-plane. A wind forced along-shore current, associated with coastal upwelling along a vertical wall, encounters a promontory of finite width and length, perpendicular to the coast.In the lower layer, topographic eddies form, which are shown to drive the formation of a filament on the front. Indeed, as the upwelling current and front develop along the coast, the along shore flow crosses the promontory, re-arranging the potential vorticity structure and generating intense vortical structures: water columns with high potential vorticity initially localized upon the promontory are advected into the deep ocean, forming cyclonic eddies, while water columns from the deep ocean with low potential vorticity climb on the topography forming a trapped anticyclonic circulation. These topographic eddies interact with the upper layer upwelling front and form an elongated, trapped and narrow filament.Sensitivity tests are then carried out and it is shown that:

• •baroclinic instability of the front does not play a major role on the formation of long trapped filaments;
• •increasing the duration of the wind forcing increases the upwelling current and limits the offshore growth of the filament;
• •modifying the promontory characteristics (width, length, height and slopes) has strong impact on the filament evolution, sometimes leading to a multipolarisation of the potential vorticity anomaly structure which results in much more complicated patterns in the upper layer (numerous shorter and less coherent filaments). This shows that only specific promontory shapes can lead to the formation of well defined filaments;
• •adding bottom friction introduces a slight generation of potential vorticity in the bottom layer over the promontory, but does not significantly alter significantly the formation of the filament along the outcropped front in the present configuration;
• •modifying the stratification characteristics, in particular the density jump between the layers, has only a weak influence on the dynamics of topographic eddies and on filament formation;
• •the influence of capes is also modest in our simulations, showing that topography plays the major role in the formation of long and trapped upwelling filaments.

     

Observations of a cycle of intense coastal upwelling and downwelling at the research site of the Shirshov Institute of Oceanology in the Black Sea

The paper presents the results of joint analysis of the response of vertical temperature and current velocity profile distributions in the coastal zone of the Gelendzhik region of the Black Sea to strong wind forcing in the third ten-day period of September 2013. This forcing was caused by the propagation of an atmospheric cyclone, which first initiated coastal upwelling that was later replaced by downwelling. We formulate a criterion for the development of full coastal upwelling and demonstrate its efficiency. We assume that frequent events of incomplete coastal upwelling and downwelling are associated with changes in the water dynamics (variations in the intensity and direction of the alongshore current) generally not related to local wind forcing.     

Modelling of upwelling and downwelling in the ocean

Numerical experiments on the reconstruction of upwelling and downwelling at the eastern boundary of the ocean were carried out in the framework of a multilayer model of the ocean involving the upper mixed layer (UML). The peculiarities of these phenomena when they are formed and attenuated owing to the strong intensification and abatement of the longshore wind have been studied. It is shown that cold waters are always involved from the thermocline to the UML during upwelling. In downwelling, this occurs as a rule. However, during upwelling the abatement of the wind may result in subduction —the inflow of warm waters from the UML to the thermocline.Translated by Mikhail M. Trufanov.