Seasonal and Interannual Variations of Heat Fluxes in the Barents Sea Region

Seasonal and interannual variations in adjective heat fluxes in the ocean (dQ_oc) and the convergence of advective heat fluxes in the atmosphere (dQ_atm) in the Barents Sea region have been investigated over the period of 1993–2012 using the results of the MIT regional eddy-permitting model and ERA-Interim atmospheric reanalysis. Wavelet analysis and singular spectrum analysis are used to reveal concealed periodicities. Seasonal 2- to 4- and 5- to 8-year cycles are revealed in the dQ_oc and dQ_atm data. It is also found that seasonal variations in dQ_oc are primarily determined by the integrated volume fluxes through the western boundary of the Barents Sea, whereas the 20-year trend is determined by the temperature variation of the transported water. A cross-wavelet analysis of dQ_oc and dQ_atm in the Barents Sea region shows that the seasonal variations in dQ_oc and dQ_atm are nearly in-phase, while their interannual variations are out-of-phase. It is concluded that the basin of the Barents Sea plays an important role in maintaining the feedback mechanism (the Bjerknes compensation) of the ocean–atmosphere system in the Arctic region.     

Vertical turbulent exchange in the Black Sea pycnocline and its relation to water dynamics

Estimates of vertical turbulent diffusion coefficient (K _t ) in the Black Sea pycnohalocline have been obtained from data of simultaneous observations of seawater temperature, salinity, density, and horizontal current velocity, obtained in the northeastern part of the Black Sea during 2013–2014 with a moored Aqualog profiler. A Munk and Andersson (1948) type parameterization, adapted for the Black Sea environment, is proposed for calculating K _t . Strong short-period (several days) variability of turbulent exchange is revealed, induced by vertical shear variations of the current velocity.     

Changes of cooling near mesopause under global warming from observations and model simulations

The results of joint analysis of temperature variations near mesopause from long-term measurements at the Zvenigorod Scientific Station of the Obukhov Institute of Atmospheric Physics RAS in 1960–2015 and variations of surface air temperature characterizing global climate change. Together with variations of temperature at the mesopause T _ms from measurements of the hydroxyl emissions we analyzed the temperature variations near mesopause T _m reduced to the same level of solar activity. The observed strong decrease in temperature near mesopause during last decades, particularly in winter, with its tendency to slow down since the 1980’s is was detected against the background of general increase in the surface air temperature of the Northern Hemisphere T _NHs and the Earth as a whole. It was revealed a sharp drop in winter temperature near mesopause in 1970s. and its synchronicity with the shift in climatic features at the surface associated with changes in formation of El Nino events and their impact on the global climate. The general significant negative correlation of temperature variations near mesopause and T _NHs detected from 56-year observational data was not accompanied by any significant coherence between the most long-period temperature variations from the cross-wavelet analysis. To assess the possible manifestation of this coherence the results of numerical simulations with a global climate model were used. According to model simulations for the 20–21 centuries taking into account anthropogenic forcings for significant coherence between long-term variations T _m and T _NHs the prolonged observations are required for temperature near mesopause–about a century or more.     

On the spatial–temporal variations in the chlorophyll-<Emphasis Type="Italic">a</Emphasis> concentration on the Peter the Great Bay shelf during the winter–spring phytoplankton bloom according to satellite and subsatellite data

Chlorophyll-a concentration (C _chl) variations in the cross section within and outside the Peter the Great Bay shelf during different stages of the winter–spring phytoplankton bloom in 2003–2005 has been considered based on a ship (obtained during the R/V Akademik M.A. Lavrent’ev voyage of February 26 to March 9, 2003) and MODIS-Aqua spectroradiometer and the SeaWiFS color-scanner satellite data. A comparison of the C _chl variability obtained from the ship and satellite data indicates that these data are inconsistent. According to satellite data obtained at the MUMM atmospheric correction, the C _chl variability is distorted less than the NIR-correction data. Studying the variations in the coefficients of light absorption by the detritus and yellow substance (a _dg) and light backscattering by suspended particles (b _bp), C _chl, chlorophyll-a fluorescence (F _chl) according to the satellite data allow us to state that the variations in the discrepancy between the satellite and ship C _chl values are mainly caused by the variations in the content of the detritus and yellow substance in water. Based on the satellite data, it has been revealed that the a _dg values increase with increasing wind mixing after the phytoplankton bloom (about 2–5 km areas where the a _dg, C _chl, F _chl, and bbp values abruptly increased in 2005, apparently due to eddy formation). It has been indicated that the F _chl characteristic, which is close to C _chl, increases when the favorable conditions for the phytoplankton bloom deteriorate. Therefore, this characteristic cannot be used to identify C _chl under the indicated conditions.     

珠江口溶解有机物的遥感动态监测

河口有色溶解有机物（colored dissolved organic matter,CDOM）的分布是各种物理-生物地球化学过程共同作用的结果。为实现河口高动态变化CDOM的监测,遥感是一种重要的手段。由珠江口四个不同季节的航次获得的实测数据,本文构建了一个遥感算法以反演CDOM在400 nm的吸收系数（a_CDOM （400））。该算法使用以波段反射率比值R_rs （667）/R_rs （443）和R_rs （748）/R_rs （412）为自变量。将构建的算法应用于2002-2014年的MODIS/Aqua数据,本文计算了珠江口不同季节的a_CDOM （400）气候态分布。CDOM的分布主要受珠江径流量和区域水下地形特征的影响。沿着垂直于水深梯度的断面,气候态a_CDOM （400）呈指数减少（y=ae^bx,b<0）,但不同季节差异很大。珠江口CDOM主要是河流淡水输运而来。其中,富里酸比例随盐度的增加而降低。基于构建的算法、CDOM保守混合方程和径流量,本文由MODIS/Aqua数据进一步估算了2002-2014年夏季和冬季珠江DOC的有效入海浓度和有效入海通量。珠江的有效入海浓度和有效入海通量都与流量存在正相关关系,且在夏季的相关性更明显,R2分别为0.698和0.9657。     

Features of the atmospheric boundary layer block for three versions of the WAM wind wave model

The estimated characteristics of the atmospheric boundary layer, obtained by the simulation of wind wave fields using three versions of the WAM numerical model are compared with the well-known empirical dependences of drag coefficient C _d on wind speed U ₁₀ and wave age A, as well as with the dependence of dimensionless roughness height z _n on inverse wave age u*/с р. Calculations carried out for several years in the areas of the Pacific and Indian oceans, based on the ERA-interim and CFSR wind reanalyses have shown good agreement between the model and empirical dependences C _d (U ₁₀) and C _d (A). The range of estimated variability for z _n (u*/с _р ) has been found to be significantly less than empirical. It has been also found that estimated values of wind speed U _10W (t) are overestimated from 5 to 10% in all versions of WAM models compared with the input wind reanalysis U _10R (t) at the moments of appearance maximum values of wind U _10R (t). The reasons for the established features of the WAM model and their dependence on the model version are discussed.     

Long-period variations in the optical and microphysical parameters of surface aerosol in a Moscow suburb

The mass concentration M of surface submicron aerosol and its condensation activity χ were measured at the Zvenigorod Scientific Station of the A.M. Obukhov Institute of Atmospheric Physics, Russian Academy of Sciences, in 2001–2006. Variations in these quantities on scales of several days to several months are analyzed. It is shown that the occurrence of dense and superdense hazes over central Russia became more frequent during these years. This resulted in a change in the character of the seasonal dependence of the mass aerosol concentration M: ordered oscillations whose amplitude increases from year to year appeared in the seasonal cycle of M. Due to more frequent cases of dense and superdense hazes in 2002–2006, a significant monotonic increase is observed in the annual means of M. The annual mean of M increased more than twofold in 2005 (56 μm/m³) compared to 2001 (23 μm/m³). The time regimes of M and χ are characterized by variations with periods in the vicinity of 10, 20, and 40–50 days, as well as variations with longer periods within 70–100 days. In addition, the variability of mass concentration is of a strongly intermittent character. A cross-spectral analysis has shown that, on the whole, variations in M are coherent with those in χ. If the parameters M and χ vary on the same time scales, then, as a rule, variations in M lag behind variations in χ by several days. To understand the relationship between the parameters M and χ, special investigations, including studies of aerosol composition, are needed.     

Cenoses of phototrophic algae of ultrasaline lakes in the Kulunda steppe (Altai krai,Russian Federation)

In 2012, expeditions of the Institute of Microbiology, Russian Academy of Sciences, delivered samples of algo-bacterial mats from Kulunda steppe alkaline lakes (Petukhovskoe alkaline lake, Tanatar VI, and Gorchina III). The filamentous alga Ctenocladus circinnatus (Chlorophyta) acted as an edificator of the mats. The composition of cenoses algocomponents also included chlorophytes Dunaliella viridis and Picocystis salinarum as well as diatoms Anomeoneis sphaerophora, Brachysira brebissonii, B. zellensis, Mastogloia pusilla var. subcapitata, Nitzschia amphibia, N. cf. communis, and Nitzschia sp. 1. The composition and structure of phototrophic algae cenoses (including diatom taxocenes) were described for the investigated lakes for the first time. For the period from 2011 to 2012, the total mineralization significantly increased in lakes. This involved sensible alterations of cenoses. B. zellensis was the most permanent component of diatom taxocenes in both seasons. In the summer of 2011, it was often accompanied by A. sphaerophora and B. brebissonii. In the summer of 2012, A. sphaerophora was found only singularly in Lake Gorchina III, and some biotopes of Lake Tanatar VI were massively inhabited by N. cf. communis, including colonies that had not been previously described for the species. The genetic analysis of three diatoms, which are markedly different from each other in their appearance and were sampled from different lakes but were all determined as Nitzschia cf. communis, showed their complete similarity to each other with the 18S rRNA gene fragment and the highest similarity of all the three diatoms with the species Nitzschia communis.     

Feeding,egg production,and respiration rate of pteropods <Emphasis Type="Italic">Limacina</Emphasis> in Arctic seas

The feeding, egg production, and respiration rate of the dominant pteropod Limacina helicina have been studied in Russia’s Arctic seas. The sinking rates of fecal pellets and dead individuals have been measured to estimate their role in vertical carbon flux. As has been shown, the rate of ecophysiological processes taking place in the pteropods is higher than that of copepods, the main consumers of phytoplankton. The gut pigment content in Limacina (3084 ng ind^–1 as a maximum) was two orders of magnitude higher than in copepods. The egg production rate in Limacina even without feeding reached 4000 eggs ind^–1 versus 350–450 egg ind^–1 typical of the dominant copepods even with excess food. A close correlation between the pteropod feeding rate and individual body weight was observed for Limacina rather than a correlation with food concentration. The experimentally estimated sinking rate of Limacina fecal pellets was 270 m day^–1, higher than for most copepods. The sinking rate of dead pteropods reaches 2000 m day^–1. According to the literature, discarded mucous feeding nets sink at a rate of 80 to 1080 m day^–1. Evidently, pteropods play a significant role in biogeochemical cycles by accelerating sedimentation. High rates of all studied processes suggest that Limacina are an important component of plankton communities and play the most important role in trophodynamics at sites of their accumulation.     

Using the acoustic-pulse conservation law in estimating the energy of surface acoustic sources by remote sounding

The atmospheric effect on the characteristics of infrasonic signals from explosions has been studied. New methods have been proposed to remotely estimate the energy of explosions using the data of infrasonic wave registration. One method is based on the law of conservation of acoustic pulse I, which is equal to the product of the wave profile area S/2 of the studied infrasonic signal and the distance to the source E_I [kt] = 1.38 × 10^–10 (I [kg/s])^1.482. The second method is based on the relationship between the explosion energy and the dominant period T of the recorded signal, EТ [kt] =1.02 × (Т [s]²/σ)^3/2, where σ is a dimensionless distance used for determining the degree of manifestation of nonlinear effects in the propagation of sound along ray trajectories. When compared to the conventional E_W (Whitaker’s) relation, the advantage of the EI relation is that it can be used for pulsed sources located at an arbitrary height over the land surface and having an arbitrary form of the initial-pulse profile and for any type of infrasonic arrivals. A distinctive feature of the expression for E_Т is that the atmospheric effect on the characteristics of recorded infrasonic signals is explicitly taken into account. These methods have been tested using infrasonic data recorded at a distance of 322 km from the sources (30 explosions caused by a fire that occurred at the Pugachevo armory in Udmurtia on June 2, 2011). For the same explosion, empirical relations have been found between energy values obtained by different methods: E_I = 1.107 × E _W , E _Т = 2.201 × E _I .     

自然死亡率变化对黄海蓝点马鲛渔业管理的影响

传统的渔业资源评估方法对数据有一定的要求,而现有的许多渔业缺乏相关数据,造成了渔业评估和管理上的障碍。特别是自然死亡率,由于观测和计算的困难,研究中一般通过经验公式等方法进行假话,而不同来源的数据和估算方法的结果存在较大差异,可能导致在渔业资源评估和管理中很大的不确定性。本研究以蓝点马鲛（Scomberomorus niphonius）为例,利用MSE比较了22种基于有限数据的渔业管理规程（MP）,筛选出DCAC、SPMSY、curE75、minlenLopt1 4种对于M不确定性较为稳健的的管理方法,并针对不同年代的自然死亡率以及通过不同计算方法得出的自然死亡率对这4种MP进行了管理策略评价。结果表明：M自1980s至今的不同年代间有显著降低趋势,该变化对MP管理效果产生了较强影响。随着M降低,种群B>0.5BMSY的概率显著上升,而过度捕捞可能性在不同MP间有方向性的差异。M变化对对curE75、minlenLopt1的长期产量的影响较小,而对DCAC、SPMSY的AAVY<15%的可能性变化较大,即降低产量的波动。不同M估算方法的估计值相差不大,对MP管理影响不明显。综合结果来看,Cur75受自然死亡率不确定性影响较小,对于渔业资源的开发和保护两个方面有着最好的权衡,最适用于我国蓝点马鲛渔业的管理。     

温度、盐度和光照对浒苔生长和光合作用的影响

Intensive Pyropia aquaculture in the coast of southwestern Yellow Sea and its subsequent waste, including disposed Ulva prolifera, was speculated to be one of the major sources for the large-scale green tide proceeding in the Yellow Sea since 2007. It was, however, unclear how the detached U. prolifera responded and resumed growing after they detached from its original habitat. In this study, we investigated the growth and photosynthetic response of the detached U. prolifera to various temperature, salinity and irradiance in the laboratory. The photosynthetic rate of the detached U. prolifera was significantly higher at moderate temperature levels(14–27°C)and high salinity(26–32), with optimum at 23°C and 32. Both low(14°C) and highest temperature(40°C), as well as low salinity(8) had adverse effects on the photosynthesis. Compared with the other Ulva species, U. prolifera showed higher saturated irradiance and no significant photoinhibition at high irradiance, indicating the great tolerance of U. prolifera to the high irradiance. The dense branch and complex structure of floating mats could help protect the thalli and reduce photoinhibition in field. Furthermore, temperature exerted a stronger influence on the growth rate of the detached U. prolifera compared to salinity. Overall, the high growth rate of this detached U. prolifera(10.6%–16.7% d~(–1)) at a wide range of temperature(5–32°C) and salinity(14–32) implied its blooming tendency with fluctuated salinity and temperature during floating. The environmental parameters in the southwestern Yellow Sea at the beginning of green tide were coincident with the optimal conditions for the detached U. prolifera.     

南极西北威德尔海上层海洋湍流混合研究

Turbulent mixing in the upper ocean(30-200 m) of the northwestern Weddell Sea is investigated based on profiles of temperature,salinity and microstructure data obtained during February 2014.Vertical thermohaline structures are distinct due to geographic features and sea ice distribution,resulting in that turbulent dissipation rates(ε) and turbulent diffusivity(K) are vertically and spatially non-uniform.On the shelf north of Antarctic Peninsula and Philip Ridge,with a relatively homogeneous vertical structure of temperature and salinity through the entire water column in the upper 200 m,both ε and K show significantly enhanced values in the order of O(10~(-7))-O(10~(-6)) W/kg and O(10~(-3))-O(10~(-2)) m~2/s respectively,about two or three orders of magnitude higher than those in the open ocean.Mixing intensities tend to be mild due to strong stratification in the Powell Basin and South Orkney Plateau,where s decreases with depth from O(10~(-8)) to O(10~(-9)) W/kg,while K changes vertically in an inverse direction relative to s from O(10~(-6)) to O(10~(-5)) m~2/s.In the marginal ice zone,K is vertically stable with the order of10~(-4) m~2/s although both intense dissipation and strong stratification occur at depth of 50-100 m below a cold freshened mixed layer.Though previous studies indentify wind work and tides as the primary energy sources for turbulent mixing in coastal regions,our results indicate weak relationship between K and wind stress or tidal kinetic energy.Instead,intensified mixing occurs with large bottom roughness,demonstrating that only when internal waves generated by wind and tide impinge on steep topography can the energy dissipate to support mixing.In addition,geostrophic current flowing out of the Weddell Sea through the gap west of Philip Passage is another energy source contributing to the local intense mixing.     

The spatial and seasonal variation of trace metals in coastal seawater and soft tissue of <Emphasis Type="Italic">Chthamalus montagui</Emphasis> around the northwest coast of Portugal

Metal concentrations (Cd, Cr, Cu, Fe, Mn and Zn) were determined in seawaters and soft tissues of Chthamalus montagui from the northwest coast of Portugal during the four seasons of 2011. The main objectives of this work were to assess seasonal and spatial variations of metals in order to detect hot spots of contamination, to establish correlations between metals in coastal seawaters and C. montagui and to calculate metal bioaccumulation factors (BAFs) in each season. Metal concentrations in coastal seawaters ranged within Cd: 1.2–35 ng L^?1; Cr: 15–87 ng L^?1; Mn: 77–1763 ng L^?1; Cu: 126–1819 ng L^?1; Fe: 430–4048 ng L^?1 and Zn: 2889–16867 ng L^?1 and in C. montagui ranged for Cd: 0.39–1.98 mg kg^?1; Cr: 0.45–3.13 mg kg^?1; Cu: 0.93–5.70 mg kg^?1; Mn: 2.2–20.4 mg kg^?1; Fe: 135–707 mg kg^?1 and Zn: 119–782 mg kg^?1. Significant spatial and seasonal variations were found between: (i) metal concentrations in seawaters and C. montagui tissues; (ii) the distribution of metal concentrations in C. montagui tissues were Fe > Zn > Mn > Cu > Cr > Cd and (iii) C. montagui showed higher bioaccumulation factors for Fe and Cd than for Cu, Mn and Zn in all seasons. Regarding the metal concentrations accumulated in C. montagui tissues during each season of 2011, the ecological quality classifications of the NW coast of Portugal varied from “Class I–Unpolluted” to “Class III–Remarkably Polluted”.     

温度和铜离子浓度对近江牡蛎超氧化物岐化酶活力的联合效应

Superoxide dismutase(SOD) is a crucial antioxidant enzyme playing the first defense line in antioxidant pathways against reactive oxygen species in various organisms including marine invertebrates. There exist mainly two specific forms, Cu/Zn-SOD(SOD1) and Mn-SOD(SOD2), in eukaryotes. SODs are known to be concurrently modulated by a variety of environmental stressors. By using central composite experimental design and response surface method, the joint effects of water temperature(18–34°C) and copper ion concentration(0.1–1.5 mg/L) on the total SOD activity in the digestive gland of Crassostrea ariakensis were studied. The results showed that the linear effect of temperature was highly significant(P0.01), the quadratic effect of temperature was significant(P0.05); the linear effect of copper ion concentration was not significant(P0.05), while the quadratic effect of copper ion concentration was highly significant(P0.01); the interactive effect of temperature and copper ion concentration was not significant(P0.05); the effect of temperature was greater than that of copper ion concentration. The model equation of digestive gland SOD enzyme activity towards the two factors of interest was established, with R2 and predictive R2 as high as 0.961 6 and 0.820 7, respectively, suggesting that the goodness-offit to experimental data be very satisfactory, and could be applied to prediction of digestive gland SOD activity in C. ariakensis under the conditions of the experiment. Our results would be conducive to addressing the health of aquatic animals and/or to detecting environmental problems by taking SOD as a potential bioindicator.     

Upper ocean heat content and atmospheric anomaly fields in the off-equatorial North Pacific related to ENSO

We have investigated interannual-scale variations of oceanic and atmospheric anomaly fields, such as upper ocean heat content (OHC), sea surface temperature (SST), latent heat flux (LHF) through the sea surface, sea level pressure (SLP) and wind stress curl (WSC) in the tropical Pacific and their relationships to El Niño/Southern Oscillation (ENSO) events. The results reported here show that the OHC and SST anomalies are almost in phase and lead LHF anomalies in the western tropical Pacific (WTP) region, which are preferable to the generation of subsequent atmospheric anomalies in the WTP. We also describe linear relationships between the amplitudes of these variables in the WTP. In addition, the results show that the both WSC and LHF anomalies are in phase with the temporal trend of OHC anomalies in the WTP, and suggest a combined effect of the local WSC and LHF anomaly in the WTP and ENSO-related, off-equatorial, westward propagating OHC anomaly to generate a large OHC anomaly in the WTP. In contrast to the WTP, OHC and SST anomalies are not in phase to the east of the WTP. The results also indicate that OHC anomalies in the WTP have a potential effect on the generation of an equatorial OHC anomaly via both a reflection of waves at the western boundary and atmospheric variations, which force the enhancement of western equatorial OHC anomaly. Therefore, the WTP is a key region where ENSO events are significantly modulated, and OHC anomalies in the WTP play an important role in the subsequent ENSO event.     

New records of three <Emphasis Type="Italic">Paraprionospio</Emphasis> species (Polychaeta: Spionidae) from Korean waters

The spionid polychaete Paraprionospio pinnata (Ehlers 1901) has been widely reported from Korean waters. We examined some specimens belonging to the genus Paraprionospio that had been collected from Korean waters, and identified them as Paraprionospio coora Wilson, 1990, Paraprionospio cordifolia Yokoyama, 2007 and Paraprionospio patiens Yokoyama, 2007, which are new to the fauna of Korea. The present study suggests misidentification of the specimens that have been previously reported as P. pinnata from Korean waters.     

Copepods <Emphasis Type="Italic">Oithona similis</Emphasis> and <Emphasis Type="Italic">Oithona davisae</Emphasis>: Two strategies of ecological–physiological adaptation in the Black Sea

Salinity tolerance, energy metabolism, buoyancy, and passive sinking and swimming speeds have been studied for comparative assessment of the adaptive potential of two species of cyclopoid copepods in the Black Sea, the native Oithona similis and new invader Oithona davisae. Both species were considered marine euryhaline copepods, but the range of salinity tolerance of O. davisae was much broader (5–55‰). The energy metabolism, locomotor activity, mean body mass density, and speed of passive sinking at the same temperature were significantly higher in O. davisae than in O. similis. The relationship between the physiological and behavioral parameters and ecological characteristics of the species are discussed.     

Seasonal dynamics and estimation of the annual primary production of phytoplankton in the Black Sea

Seasonal variations in the surface chlorophyll a concentrations (Chl _s) and the integrated primary production (PP _inf) were investigated for ten regions of the Black Sea based on long term observations (1973–1997). Two or three maximums of both Chl _s and PP _inf were registered in most of the shelf regions (SR, <200 m), the continental slope (CS, 200–1500 m), and the deep regions (DSR >1500 m) in February–March, June–August, and October–November. Such a pattern suggests that the seasonal dynamics of PP _inf strongly depend on the Chl _s variability. The mean annual values of the PP _inf comprised 130–420, 130–150, and 140–150 g C m^?2 in the SR, CS, and DSR, respectively. These values are mainly typical of the eutrophic layer and the transition between the eutrophic and mesotrophic waters (SR) or for the upper boundary of the mesotrophic waters (CS and DSR). The maximal contribution of the wintertime (December–March) to the total PP _inf values (40–42%) was observed in the DSR. In the SR and the adjacent eastern CS areas, the proportion of the PP _inf summertime production (June–September) reaches 40–60% and is higher than the wintertime production. The lowest values of PP _inf (9–17%) were produced in the spring and autumn periods. The total annual values of PP _inf in the Black Sea are close to 50–70 Mt C.     

Seasonal and interannual variation in mesozooplankton community structure off Tongyeong,southeastern coast of Korea,from 2011 to 2014

Mesozooplankton community structure and environmental factors were monitored monthly at a fixed station off Tongyeong, southeastern coast of Korea, from 2011 to 2014 to better understand the variability of the mesozooplankton community in relation to changes in the marine environment. Total mesozooplankton density varied from 747 to 8,945 inds. m^-3 with peaks in summer. The surface water temperature (r = 0.338, p < 0.05) and chlorophyll-a (Chl-a) concentration (r = 0.505, p < 0.001) were parts of the factors that may have induced the mesozooplankton peaks in summer. Copepods accounted for 71% of total mesozooplankton. Total copepod density, particularly cyclopoid copepods, increased during the study period. Cumulative sum plots and anomalies of the cyclopoid copepod density revealed a change of the cyclopoid density from negative to positive in June 2013. A positive relationship between cyclopoid copepods and the Chl-a concentration (r = 0.327, p < 0.05) appeared to be one of the reasons for the increase in cyclopoids. Dominant mesozooplankton species such as Paracalanus parvus s.l., Oikopleura spp., Evadne tergestina, Cirripedia larvae, Corycaeus affinis, Calanus sinicus, and Oithona similis accounted for 60% of total mesozooplankton density. Based on cluster analysis of the mesozooplankton community by year, the seasonal distinction among groups was different in 2014 compared to other years. P. parvus s.l. and its copepodites contributed most in all groups in all four years. Our results suggest that the high Chl-a concentration since 2013 may have caused the changes in mesozooplankton community structure in the study area.