Application of an unstructured mesh model to storm surge propagation in the Mersey estuary region of the Irish Sea

The storm surge period of 13–16 November 1977 when there was a major positive surge followed by a negative surge in the Irish Sea is investigated using a two-dimensional unstructured mesh model of the west coast of Britain. The model accounts for tidal and external surge forcing across its open boundaries which are situated in the Celtic Sea and off the west coast of Scotland. Although this period has been examined previously using a uniform finite-difference model, and a finite element model, neither of these could resolve the Mersey estuary which is the focus of the present study. By using a finite element model with very high mesh resolution within the Mersey, the spatial variability of surge elevations and currents within the Mersey to rapidly changing surge dynamics can be examined. The mesh in the model varies from about 7 km in deep water, to the order of 100 m in the Mersey, with the largest mesh length reaching 17 km in deep offshore regions, and smallest of order 26 m occurring in shallow coastal regions of the Mersey estuary. The model accounts for wetting/drying which occurs in shallow water coastal areas. Calculations showed that during the positive surge period, the amplitude and speed of propagation of the surge was largest in the deep water channels. This gave rise to significant spatial variability of surge elevations and currents within the estuary. As wind stresses decreased over the Irish Sea, a negative surge occurred over Liverpool Bay and at the entrance to the Mersey. However, within the Mersey there was a local positive surge which continued to propagate down the estuary. This clearly showed that although the large scale response of the Irish Sea to changing wind fields occurred rapidly, the response in the Mersey was much slower. These calculations with a west coast variable mesh model that included a high-resolution representation of the Mersey revealed for the first time how elevations and currents within the Mersey responded to Irish Sea surges that rapidly changed from positive to negative.     

Modelling tide–surge interaction effects using finite volume and finite element models of the Irish Sea

An unstructured mesh model of the west coast of Britain, covering the same domain and using topography and open boundary forcing that are identical to a previous validated uniform grid finite difference model of the region, is used to compare the performance of a finite volume (FV) and a finite element (FE) model of the area in determining tide–surge interaction in the region. Initial calculations show that although qualitatively both models give comparable tidal solutions in the region, comparison with observations shows that the FV model tends to under-estimate tidal amplitudes and hence background tidal friction in the eastern Irish Sea. Storm surge elevations in the eastern Irish Sea due to westerly, northerly and southerly uniform wind stresses computed with the FV model tend to be slightly higher than those computed with the FE model, due to differences in background tidal friction. However, both models showed comparable non-linear tide–surge interaction effects for all wind directions, suggesting that they can reproduce the extensive tide–surge interaction processes that occur in the eastern Irish Sea. Following on from this model comparison study, the physical processes contributing to surge generation and tide–surge interaction in the region are examined. Calculations are performed with uniform wind stresses from a range of directions, and the balance of various terms in the hydrodynamic equations is examined. A detailed comparison of the spatial variability of time series of non-linear bottom friction and non-linear momentum advection terms at six adjacent nodes at two locations in water depths of 20 and 6 m showed some spatial variability from one node to another. This suggests that even in the near coastal region, where water depths are of the order of 6 m and the mesh is fine (of order 0.5 km), there is significant spatial variability in the non-linear terms. In addition, distributions of maximum bed stress due to tides and wind forcing in nearshore regions show appreciable spatial variability. This suggests that intensive measurement campaigns and very high-resolution mesh models are required to validate and reproduce the non-linear processes that occur in these regions and to predict extreme bed stresses that can give rise to sediment movement. High-resolution meshes will also be required in pollution transport problems.     

Surge modelling in the eastern Irish Sea: present and future storm impact

It is believed that, in the future, the intensity and frequency of extreme coastal flooding events may increase as a result of climate change. The Natural Environment Research Council (NERC) Flood Risk from Extreme Events (FREE) project, Coastal Flooding by Extreme Events and EU FP7 Morphological Impacts and Coastal Risks Induced by Extreme Storm Events project are investigating the flood risks in the eastern Irish Sea, an area that includes most of England’s coastal types. Using a previously modelled and validated historical extreme surge event, in November 1977, we now investigate the changes in peak surge as a result of possible future climate conditions. In order to simulate the surge, we have set up a one-way nested approach, using the Proudman Oceanographic Laboratory Coastal Ocean Modelling System 3D baroclinic model, from a domain covering the whole NW European continental shelf, through to a 1.85 km Irish Sea model; both areas are forced by tides, atmospheric pressure and winds. We use this modelling system to investigate the impact of enhanced wind velocities and increased sea levels on the peak surge elevation and residual current pattern. The results show that sea level rise has greater potential to increase surge levels than increased wind speeds.     

Effects of wave-induced forcing on a circulation model of the North Sea

The effect of wind waves on water level and currents during two storms in the North Sea is investigated using a high-resolution Nucleus for European Modelling of the Ocean (NEMO) model forced with fluxes and fields from a high-resolution wave model. The additional terms accounting for wave-current interaction that are considered in this study are the Stokes-Coriolis force, the sea-state-dependent energy and momentum fluxes. The individual and collective role of these processes is quantified and the results are compared with a control run without wave effects as well as against current and water-level measurements from coastal stations. We find a better agreement with observations when the circulation model is forced by sea-state-dependent fluxes, especially in extreme events. The two extreme events, the storm Christian (25–27 October 2013), and about a month later, the storm Xaver (5–7 December 2013), induce different wave and surge conditions over the North Sea. Including the wave effects in the circulation model for the storm Xaver raises the modelled surge by more than 40 cm compared with the control run in the German Bight area. For the storm Christian, a difference of 20–30 cm in the surge level between the wave-forced and the stand-alone ocean model is found over the whole southern part of the North Sea. Moreover, the modelled vertical velocity profile fits the observations very well when the wave forcing is accounted for. The contribution of wave-induced forcing has been quantified indicating that this represents an important mechanism for improving water-level and current predictions.     

Storm surge computations in estuarine and near-coastal regions: the Mersey estuary and Irish Sea area

An unstructured grid storm surge model of the west coast of Britain, incorporating a high-resolution representation of the Mersey estuary is used to examine storm surge dynamics in the region. The focus of the study is the major surge that occurred during the period 11–14 November 1977, which has been investigated previously using uniform grid finite difference models and a finite element model of the west coast of Britain. However, none of these models included the Mersey estuary. Comparison of solutions in the eastern Irish Sea with those computed with these earlier models showed that, away from the Liverpool Bay region, the inclusion of the Mersey estuary had little effect. However, at the entrance to the Mersey, its inclusion did influence the solution. By including a detailed representation of the Mersey estuary within the model, it was possible to conduct a detailed study of storm surge propagation in the Mersey, which had never previously been performed. This detailed study showed for the first time that the surge’s temporal variability within the estuary is influenced by surge elevation at its entrance. This varies with time as a function of spatial and temporal variations of wind stress over the west coast of Britain. Within the Mersey, calculations show that the spatial variability is mainly determined by changes in bottom topography, which had not been included in earlier finite difference models. However, since water depth is influenced by variations in tidal elevation, this, together with tide surge interaction through bottom friction and momentum advection, influences the surge. The ability of the finite element model to vary the mesh in near-shore regions to such an extent that it can resolve the Mersey and hence the impact of the Mersey estuary upon the Liverpool Bay circulation shows that it has distinct advantages over earlier finite difference models. In the absence of detailed measurements within the Mersey at the time of the surge, it was not possible to validate predicted surge elevations within the Mersey. However, significant insight into physical processes influencing the surge propagation down the estuary, its reflection and spatial/temporal variability could be gained.     

Assessing changes in extreme sea levels: Application to the English Channel, 1900–2006

A recently extended and spatially rich English Channel sea level dataset has been used to evaluate changes in extreme still water levels throughout the 20th century. Sea level records from 18 tide gauges have been rigorously checked for errors and split into mean sea level, tidal and non-tidal components. These components and the interaction between surge and tide have been analysed separately for significant trends before determining changes in extreme sea level. Mean sea level is rising at 0.8–2.3 mm/year, depending on location. There is a small increase (0.1–0.3 mm/year) in the annual mean high water of astronomical tidal origin, relative to mean sea level, and an increase (0.2–0.6 mm/year) in annual mean tidal range. There is considerable intra- and inter-decadal variability in surge intensity with the strongest intensity in the late 1950s. Storm surges show a statistically significant weak negative correlation to the winter North Atlantic Oscillation index throughout the Channel and a stronger significant positive correlation at the boundary with the southern North Sea. Tide–surge interactions increase eastward along the English Channel, but no significant long-term changes in the distribution of tide–surge interaction are evident. In conclusion, extreme sea levels increased at all of the 18 sites, but at rates not statistically different from that observed in mean sea level.     

Tide–surge interactions and their effects on total sea levels in Irish coastal waters

The research presented in this paper involves the application of the joint probability method to the estimation of extreme water levels resulting from astronomical tides and surge residuals and the investigation of the effects of tide–surge interactions on extreme water levels. The distribution of tide peaks was analysed from field records (<20 years) and a 46-year dataset of monthly maximum tidal amplitudes. Large surges were extracted from both field records and a numerical model hindcast covering the 48 largest storm events in the Irish Sea over the period 1959–2005. Extreme storm surges and tides were independently modelled using the generalised extreme value statistical model, and derived probability distributions were used to compute extreme water levels. An important, and novel, aspect of this research is an analysis of tide–surge interactions and their effects on total water level; where interactions exist, they lead to lower total water levels than in the case of independency. The degree of decrease varies with interaction strength, magnitude of surge peak at a particular phase of tide and the distribution of peaks over a tidal cycle. Therefore, including interactions in the computation of extreme levels may provide very useful information at the design stage of coastal protection systems.     

Post-glacial sediment dynamics in the Irish Sea and sediment wave morphology: Data–model comparisons

The irregular seafloor of the narrow Irish Sea on the NW European Shelf has been documented over several decades. From recently collected swath bathymetry data, very large trochoidal, nearly symmetrical sediment waves are observed in many parts of the Irish Sea and appear similar to those described from other continental shelf seas in North America that were covered by glacigenic sediments during the Last Glacial Maximum. Swath multibeam and single beam bathymetry data, backscatter intensity, shallow seismic imagery, video footage and sediment cores from the Irish Sea high sediment waves have been integrated to identify their genesis with reference to present and past hydrodynamic variability. From cross-sectional profiles over asymmetrical sediment waves in the Irish Sea the direction of asymmetry is used to map residual bed stress directions and associated bedload transport paths. Irish Sea peak bed stress vectors were generated using a two-dimensional palaeo-tidal model for the NW European shelf seas and compare well with the observations. Tidally induced bed stresses are modelled to have increased between 7–10 ka BP, to be nearly symmetrical in magnitude and to have reversed in dominant direction on a millennial scale. These environmental conditions during the post-glacial marine transgression are suggested here to help comprehend the construction of the very large sediment waves, with local variations due to differences in sediment grain size, sediment supply, water depth and intensified currents due to seafloor slopes. Model parameterisation using an open ocean boundary with time-dependent tidal changes and the implementation of high-resolution bathymetric information will improve future models of small-scale bed shear stress patterns and improve the predictive value of such modelling efforts.     

On the sensitivity of computed surges to open-boundary formulation

A coarse-grid (resolution of order 7 km) model of the west coast of Britain is used to examine the sensitivity of computed storm-surge elevations and currents to a range of open-boundary conditions. The storm-surge period 1 to 26 March 1994 is used for this comparison, as it is a time of significant wind activity. Also current measurements in the North Channel of the Irish Sea together with coastal elevation measurements are available for model validation. Elevations and currents previously computed with a coarse-grid shelf-wide model can also be incorporated into the open-boundary condition to examine the influence of far-field effects. Initial model calculations with no far-field input show the importance of including shelf-wide effects from either the external shelf model, or by using observations from coastal gauges interpolated along the open boundary of the west-coast model. Provided the west-coast models open boundary is taken sufficiently far away from the region of interest, in this case the Irish Sea, then either a radiation condition or an elevation-specified condition is appropriate provided far-field effects are taken into account. If these are not included, then neither boundary condition is successful. For the radiation condition it is necessary to include both elevations and currents from a far-field model in order to reproduce the surge. In the case of an elevation-specified boundary condition far-field effects can be incorporated in hindcast calculations by including observed sea-level changes. In a storm-surge prediction calculation the radiation condition with a far-field model is required. Calculations show that computed elevations are spatially more coherent than currents, with flows through the western Irish Sea showing the greatest sensitivity to open-boundary formulation during storm events.Responsible Editor: Phil Dyke     

A two-way nested-grid ocean-circulation model for the Meso-American Barrier Reef System

A two-way nested-grid ocean-circulation model is developed for the Meso-American Barrier Reef System (MBRS), using a newly developed two-way interactive nesting technique. The unique feature of this new nesting technique is its use of the semi-prognostic method (Sheng et al. 2001) to exchange information between different grids through the model momentum equations. The nested-grid model for the MBRS has a fine-resolution inner model embedded in a coarse-resolution outer model. The outer model is the western Caribbean Sea model developed by Sheng and Tang (2003), with a horizontal resolution of roughly 19 km. The inner model domain covers the northwest Caribbean Sea (NWCS) between 79°W and 89°W and between 15.5°N and 22°N, with a horizontal resolution of roughly 6 km. The nested-grid ocean model is initialized with the January mean temperature and salinity and forced by the monthly mean COADS (comprehensive ocean-atmosphere data set) wind stress and surface heat flux. The model sea-surface salinity is restored to the monthly mean climatology. The nested-grid model is integrated for 2 years and the second-year model results are presented in this paper. The model-calculated annual-mean near-surface currents over the NWCS agree reasonably well with the time-mean near-surface currents inferred by Fratantoni (2001) from trajectories of the satellite-tracked 15-m drogued drifters in the 1990s. The two-way nested model is also used to quantify the role of local wind stress, local density gradients and boundary forcings of the outer model in driving the annual-mean circulation in the region.Responsible Editor: Phil Dyke     

Influence of sea surface wind wave turbulence upon wind-induced circulation, tide–surge interaction and bed stress

A three-dimensional finite volume unstructured mesh model of the west coast of Britain, with high resolution in the coastal regions, is used to investigate the role of wind wave turbulence and wind and tide forced currents in producing maximum bed stress in the eastern Irish Sea. The spatial distribution of the maximum bed stress, which is important in sediment transport problems, is determined, together with how it is modified by the direction of wind forced currents, tide–surge interaction and a surface source of wind wave turbulence associated with wave breaking. Initial calculations show that to first order the distribution of maximum bed stress is determined by the tide. However, since maximum sediment transport occurs at times of episodic events, such as storm surges, their effects upon maximum bed stresses are examined for the case of strong northerly, southerly and westerly wind forcing. Calculations show that due to tide–surge interaction both the tidal distribution and the surge are modified by non-linear effects. Consequently, the magnitude and spatial distribution of maximum bed stress during major wind events depends upon wind direction. In addition calculations show that a surface source of turbulence due to wind wave breaking in shallow water can influence the maximum bed stress. In turn, this influences the wind forced flow and hence the movement of suspended sediment. Calculations of the spatial variability of maximum bed stress indicate the level of measurements required for model validation.     

Modelling travel and residence times in the eastern Irish Sea

The Irish Sea, which lies between 51 degrees N-56 degrees N and 2 degrees 50'W-7 degrees W, provides a sheltered environment to exploit valuable fisheries resource. Anthropogenic activity is a real threat to its water quality. The majority of freshwater input down rivers flows into the eastern Irish Sea. The structure of the water circulation was not well understood during the planning of Sellafield nuclear plant outfall site in the eastern Irish Sea. A three-dimensional primitive equation numerical model was applied to the Irish Sea to simulate both barotropic and baroclinic circulation within the region. High accuracy was achieved with regard to the prediction of both tidal circulation and surface and nearbed water temperatures across the region. The model properly represented the Western Irish Sea Gyre, induced by thermal stratification and not known during planning Sellafield. Passive tracer simulations based on the developed hydrodynamic model were used to deliver residence times of the eastern Irish Sea region for various times of the year as well as travel times from the Sellafield outfall site to various locations within the Irish Sea. The results indicate a strong seasonal variability of travel times from Sellafield to the examined locations. Travel time to the Clyde Sea is the shortest for the autumnal tracer release (90 days); it takes almost a year for the tracer to arrive at the same location if it is released in January. Travel times from Sellafield to Dublin Bay fall within the range of 180-360 days. The average residence time of the entire eastern Irish Sea is around 7 months. The areas surrounding the Isle of Man are initially flushed due to a predominant northward flow; a backwater is formed in Liverpool Bay. Thus, elevated tracer concentrations are predicted in Liverpool Bay in the case of accidental spills at the Sellafield outfall site.     

Modelling of waves in the Irish Sea: effects of oceanic wave and wind forcing

This study uses a series of scenarios of wave (boundary) and wind (local) forcing to examine the sensitivity and to quantify the effects associated with nesting ProWAM and POLCOMS models for downscaling predictions of waves in the Irish Sea. The model results show that the response of the modelling system to the wave and wind forcing during the downscaling varies widely depending on wind conditions. Generally, the wave forcing has a greater effect on overall wave prediction in most of the Irish Sea, except for the eastern Irish Sea/Liverpool Bay. The study also suggests detailed look-up tables at specific locations to quantify the impacts of the different forcing scenarios over the Irish Sea, which can be readily extended to the location on any other sites.     

Storm surge computations for the west coast of Britain using a finite element model (TELEMAC)

An unstructured mesh finite element model of the sea region off the west coast of Britain is used to examine the storm surge event of November 1977. This period is chosen because accurate meteorological data to drive the model and coastal observations for validation purposes are available. In addition, previous published results from a coarse-grid (resolution 7 km) finite difference model of the region and high-resolution (1 km) limited area (namely eastern Irish Sea) model are available for comparison purposes. To enable a “like with like” comparison to be made, the finite element model covers the same domain and has the same meteorological forcing as these earlier finite difference models. In addition, the mesh is based on an identical set of water depths. Calculations show that the finite element model can reproduce both the “external” and “internal” components of the surge in the region. This shows that the “far field” (external) component of the surge can accurately propagate through the irregular mesh, and the model responds accurately, without over- or under-damping, to local wind forcing. Calculations show significant temporal and spatial variability in the surge in close agreement with that found in earlier finite difference calculations. In addition, root mean square errors between computed and observed surge are comparable to those found in previous finite different calculations. The ability to vary the mesh in nearshore regions reveals appreciable small-scale variability that was not found in the previous finite difference solutions. However, the requirement to perform a “like with like” comparison using the same water depths means that the full potential of the unstructured grid model to improve resolution in the nearshore region is inhibited. This is clearly evident in the Mersey estuary region where a higher resolution unstructured mesh model, forced with uniform winds, had shown high topographic variability due to small-scale variations in topography that are not resolved here. Despite the lack of high resolution in the nearshore region, the model showed results that were consistent with the previous storm surge models of the region. Calculations suggest that to improve on these earlier results, a finer nearshore mesh is required based upon accurate nearshore topography.     

Climate change and North Sea storm surge extremes: an ensemble study of storm surge extremes expected in a changed climate projected by four different regional climate models

The coastal zones are facing the prospect of changing storm surge statistics due to anthropogenic climate change. In the present study, we examine these prospects for the North Sea based on numerical modelling. The main tool is the barotropic tide-surge model TRIMGEO (Tidal Residual and Intertidal Mudflat Model) to derive storm surge climate and extremes from atmospheric conditions. The analysis is carried out by using an ensemble of four 30-year atmospheric regional simulations under present-day and possible future-enhanced greenhouse gas conditions. The atmospheric regional simulations were prepared within the EU project PRUDENCE (Prediction of Regional scenarios and Uncertainties for Defining EuropeaN Climate change risks and Effects). The research strategy of PRUDENCE is to compare simulations of different regional models driven by the same global control and climate change simulations. These global conditions, representative for 1961–1990 and 2071–2100 were prepared by the Hadley Center based on the IPCC A2 SRES scenario. The results suggest that under future climatic conditions, storm surge extremes may increase along the North Sea coast towards the end of this century. Based on a comparison between the results of the different ensemble members as well as on the variability estimated from a high-resolution storm surge reconstruction of the recent decades it is found that this increase is significantly different from zero at the 95% confidence level for most of the North Sea coast. An exception represents the East coast of the UK which is not affected by this increase of storm surge extremes.     

Influence of seasonal circulation on flushing of the Irish Sea

We applied a three-dimensional general ocean and coastal circulation model to the Irish Sea in order to determine water renewal time scales in the region. The model was forced with meteorological data for 1995, a year with relatively warm summer and when extensive hydrographic surveys were conducted in the Irish Sea. We investigated intra-annual variability in the rates of net flow through the Irish Sea and carried out several flushing simulations based on conservative tracer transport. The results indicate that the net northward flow of 2.50 km³/d is seasonally highly variable and under certain conditions is reversed to southward. The variability in obtained residence times is high; baroclinic effects are significant. Obtained results point at the importance of spatial and temporal consideration for transport of pollutants in the shelf seas. Implications for management are numerous and involve activities such as transport, fishing, use of resources, nature conservation, monitoring, tourism and recreation.     

Storm surge simulation along the U.S. East and Gulf Coasts using a multi-scale numerical model approach

The effectiveness of simulating surge inundation using the Eulerian–Lagrangian circulation (ELCIRC) model over multi-scale unstructured grids was examined in this study. The large domain model grid encompasses the western North Atlantic Ocean, the Gulf of Mexico, and the Caribbean Sea to appropriately account for remote and resonance effects during hurricane events and simplify the specification of the open boundary condition. The U.S. East and Gulf Coasts were divided into 12 overlapping basins with fine-resolution (up to 30 × 30 m) grids to model overland surge flooding. These overlapping basins have different fine-resolution grids near the coastal region, but have an identical coarse-resolution grid in the offshore region within the large model domain. Thus, the storm surge prediction can be conducted without reducing computation efficiency by executing multiple model runs with local fine-resolution grids where potential hurricane landfalls may occur. The capability of the multi-scale approach was examined by simulating storm surge caused by Hurricanes Andrew (1992) and Isabel (2003) along the South Florida coast and in the Chesapeake Bay. Comparisons between simulated and observed results suggest that multi-scale models proficiently simulated storm surges in the Biscayne Bay and the Chesapeake Bay during two hurricanes. A series of sensitivity tests demonstrated that the simulation of surge flooding was improved when LiDAR topographic data and special bottom drag coefficient values for mangrove forests were employed. The tests also showed that appropriate representation of linear hydrologic features is important for computing surge inundation in an urban area.     

Spatial and temporal variability in nutrient concentrations in Liverpool Bay,a temperate latitude region of freshwater influence

This paper presents data for the temporal and spatial distribution of nutrients in Liverpool Bay between 2003 and 2009 and an analysis of inputs of nutrients from the major rivers. The spatial distribution of winter nutrient concentrations are controlled by the region of freshwater influence (ROFI) in Liverpool Bay through the mixing of riverine freshwater and Irish Sea water, with strong linear relationships between nutrient concentration and salinity between December and February. The location of highest spring and summer phytoplankton biomass reflects the nutrient distributions as controlled by the ROFI. Analysis of 7 years of data showed that the seasonal cycle of winter maximum nutrient concentrations in February and drawdown in April/May is a recurrent feature of this location, with the timing of the drawdown varying by several weeks between years. A comparison of observed nutrient concentrations in Liverpool Bay with those predicted from inputs from rivers has been presented. Nutrient concentrations in the rivers flowing into Liverpool Bay were highly variable and there was reasonable agreement between predicted freshwater nutrient concentrations using data from this study and riverine nutrient concentrations weighted on the basis of river flow, although the exact nature of mixing between the rivers could not be determined. Predicted Irish Sea nutrient concentrations in the winter were lower than those reported for the input waters of the North Atlantic, supporting findings from previous work that nitrogen is lost through denitrification in the Irish Sea.     

Observation of the seasonal evolution of the Yellow Sea Cold Water Mass in 1996–1998

We use the hydrographic data obtained during the joint survey of the Yellow Sea by the First Institute of Oceanography, China and the Korea Ocean Research and Development Institute, Korea, to quantify the spatial structures and temporal evolution of the southern Yellow Sea Cold Water Mass (YSCWM). It is indicated that the southern YSCWM is a water mass that develops in summer and decays in fall. In winter, due to the intrusion of the Yellow Sea Warm Current (YSWC), the central area (approximately between 34°N and 35°N, 122°E and 124°E) of the Yellow Sea is mainly occupied by relatively high temperature water (T>10 °C). By contrast, from early summer to fall, under the seasonal thermocline, the central area of Yellow Sea is occupied by cold water (T<10 °C). In summer, the southern YSCWM has two cold cores. One is formed locally southeast of Shandong Peninsula, and the other one has a tongue-like feature occupying the area approximately between 34°N and 37°N, 123°E and 126°E. The bottom layer temperature anomalies from February to July in the cold tongue region, along with the trajectories of the bottom floaters, suggest that the cold water mass in the northeast region has a displacement from the north to the central area of the Yellow Sea during the summer.