Estimating the characteristics of the vertical turbulent viscosity in the upper 200-m layer of the Black Sea

The influence of high vertical velocity gradients in the Black Sea Rim Current on the intensity of the vertical turbulent exchange is demonstrated on the basis of numerical modeling based on CTD data. The vertical turbulent exchange is confirmed by the anomalous distribution of the hydrochemical parameters in the redox layer. A system of equations for the kinetic energy of the turbulence and dissipation rate (k-? model) is used for the calculation of the coefficient of the vertical turbulent viscosity (diffusivity).     

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

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.     

Evaluation of spatial distribution of turbulent mixing in the central Pacific

A long-term mean turbulent mixing in the depth range of 200–1000 m produced by breaking of internal waves across the middle and low latitudes (40°S–40°N) of the Pacific between 160°W and 140°W is examined by applying fine-scale parameterization depending on strain variance to 8-year (2005–2012) Argo float data. Results show that elevated turbulent dissipation rate (ε) is related to significant topographic regions, along the equator, and on the northern side of 20°N spanning to 24°N throughout the depth range. Two patterns of latitudinal variations of ε and the corresponding diffusivity (K_ρ) for different depth ranges are confirmed: One is for 200–450 m with significant larger ε and K_ρ, and the maximum values are obtained between 4°N and 6°N, where eddy kinetic energy also reaches its maximum; The other is for 350–1000 m with smaller ε and K_ρ, and the maximum values are obtained near the equator, and between 18°S and 12°S in the southern hemisphere, 20°N and 22°N in the northern hemisphere. Most elevated turbulent dissipation in the depth range of 350–1000 m relates to rough bottom roughness (correlation coefficient?=?0.63), excluding the equatorial area. In the temporal mean field, energy flux from surface wind stress to inertial motions is not significant enough to account for the relatively intensified turbulent mixing in the upper layer.     

Effect of the nonlinear interaction of tidal harmonics on their spatial structure as applied to the system of the white and Barents Seas

Using the three-dimensional nonlinear finite-element thermohydrodynamic model QUODDY-4, we obtain two solutions giving an idea of the role of the nonlinear interaction of tidal harmonics in the formation of their spatial structure. The first solution is induced by the total tide (M ₂ + S ₂ + K ₁ + O ₁) at the open boundary and by the total static tide inside the area under consideration; then, the solution obtained is subject to harmonic analysis. The second solution is obtained by specifying such tidal sea level elevations at the open boundary that meet individual tidal harmonics. These two solutions are compared. It is shown that the differences between the solutions for the S ₂, K ₁, and O ₁ waves can be significant, especially near the open boundary between the White and Barents seas. This conclusion remains valid also for the maximum velocities (major semiaxes of the ellipses) of the barotropic (depth-average) tidal current as well as for the average (over the tidal cycle) densities of the total tidal energy and components of the tidal-energy budget. The emergence of this feature indicates that there are resonance modes with frequencies that differ from those of the S ₂, K ₁, and O ₁ harmonics to a lesser extent than the M ₂ harmonic frequency. The same conclusion can be made by comparing the values of the amplification factor, which is defined as the ratio between the actual and static tides, for the system of the Barents and White seas as a whole and for Mezen Bay in the White Sea and Czech Bay in the Barents Sea, taken separately.     

Variability of Particulate Organic Phosphorus in the Northwestern Part of the Black Sea

Based on long-term (1985–1995) monitoring data, the paper considers the peculiarities of seasonal variability in the spatial and vertical distribution of particulate organic phosphorus (Р_POM) in the surface layer and in the photosynthetic zone in the northwestern Black Sea. Regression equations, experimental data, and satellite observations for the chlorophyll a concentration allowed us to evaluate the seasonal longterm (1979–1995) variability in Р_POM in the surface layer and photosynthesis zone. The ratios of the concentrations of particulate organic carbon, nitrogen, phosphorus, and chlorophyll a are calculated and statistical estimates of seasonal changes in the Р_POM in the areas with different degrees of influence of river runoff and water of open seas are obtained. The consistency of intra-annual changes in the concentrations of Р_POM, chlorophyll a, and phytoplankton biomass is shown, which indicates the role of phytoplankton in the formation of Р_POM and in its intra- and interannual variability in the northwestern part of the sea. It is shown that long-term seasonal variations in Р_POM and related changes in the concentration of chlorophyll a depend on the variability of bulk river runoff, the extent of its abundance in the northwestern shelf, and regional hydrometeorological conditions.     

2012年夏季海南岛东岸上升流区的混合观测

The turbulent mixing in the upwelling region east of Hainan Island in the South China Sea is analyzed based on in situ microstructure observations made in July 2012. During the observation, strong upwelling appears in the coastal waters, which are 3℃ cooler than the offshore waters and have a salinity 1.0 greater than that of the offshore waters. The magnitude of the dissipation rate of turbulent kinetic energy ε in the upwelling region is O(10–9 W/kg), which is comparable to the general oceanic dissipation. The inferred eddy diffusivity K_ρ is O(10–6 m~2/s), which is one order of magnitude lower than that in the open ocean. The values are elevated to K_ρ≈O(10–4 m~2/s) near the boundaries. Weak mixing in the upwelling region is consistent with weak instability as a result of moderate shears versus strong stratifications by the joint influence of surface heating and upwelling of cold water.The validity of two fine-scale structure mixing parameterization models are tested by comparison with the observed dissipation rates. The results indicate that the model developed by Mac Kinnon and Gregg in 2003 provides relatively better estimates with magnitudes close to the observations. Mixing parameterization models need to be further improved in the coastal upwelling region.     

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.     

Vertical Structure of Currents in Northern Black Sea,October 2016

The paper presents the results of an analysis of hydrological data obtained on cruise 89 of the R/V Professor Vodyanitskii from September 30 through October 19, 2016, in the central area of the northern Black Sea. Currents were measured using a lowered acoustic Doppler current profiler (LADCP) to a depth of 350 m. Examples of the current velocity profiles and conditional density at individual stations are presented. The averaged characteristics of the vertical current structure in the seasonal and permanent pycnoclines are discussed. The isopycnally averaged kinetic energy (KE) demonstrates an almost linear dependence on the conditional density (σ_θ). In the seasonal pycnocline, EK = 0.2–0.013σ_θ; in the permanent pycnocline, EK = 0.12–0.007σ_θ. An interpretation is given for the variability of the kinetic energy with depth. It is shown that the maxima of the averaged vertical shear profile correspond to the occurrence depth of the seasonal and permanent pycnoclines. The mean Richardson number has a value of ~10 in the upper sea layer and decreases with depth, assuming a value of ~3 at the 300 m horizon. The estimate for the vertical turbulent diffusion coefficient from the simplified model shows its increase with depth.     

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.     

Features of the spatial distribution of phytoplankton in Nhatrang Bay of the South China Sea during the rainy season

The species composition, phytoplankton abundance, and relative yield of the variable fluorescence (F _v /F _m ) were determined in the mesotrophic Nhatrang Bay in October–November of 2004. The species diversity (250 taxonomic units) and heterogeneity of the phytoplankton structure were high. With respect to the number of species and their abundance, diatoms prevailed. In selected parts of the bay, dinoflagellates dominated. The mean biomass in the water column under 1 m² (B _t ) varied from 2.3 to 64.4 mg C/m³ being 31.0 mg C/m³ on average. The values of B _t were the lowest at the stations nearest to the river mouth. Seaward, B _t increased. The values of B _t increased with depth at some stations and decreased at others. In the surface sea layers, the biomass was lower than that in the underlying waters. The values of F _v /F _m ranged from 0.10 to 0.64 (at a mean value of 0.49). The lowest values of F _v /F _m were observed in the area close to the seaport. Over the greater part of the bay, the values of F _v /F _m were higher than 0.47. Such values are indicative of the relatively high potential photosynthetic activity of the phytoplankton. The abundance and species diversity were higher than those in the dry season (March–April).     

Feeding of the Dominant Herbivorous Plankton Species in the Black Sea and Their Role in Coccolithophorid Consumption

The feeding of abundant herbivore plankton species has been investigated during different stages of coccolithophorid development in the northeastern Black Sea, namely, at the initial stage of development (March 2009) and at the massive bloom stage (June 2011 and May 2013). The role of coccolithophorids as a food source for Black Sea copepods Calanus euxinus, Acartia clausi, and Pseudocalanus elongatus; the heterotrophic dinoflagellate Noctiluca scintillans; and the larvacean Oikopleura dioica has been characterized for the first time. The contribution of coccolithophorids to copepod ration was low (0.2 to 13% of total carbon intake) during all study periods. The content of these organisms in the diet of N. scintillans ranged from 17 to 100%. The coccolithophorid bloom apparently had a stronger positive effect on the nutrition of fine filter feeders, such as N. scintillans and O. dioica, than on copepod nutrition. Daily coccolithophorid consumption by zooplankton ranged from 0.7 to 39.4% of the biomass in different study periods. The grazing of N. scintillans and O. dioica populations made the greatest contribution to coccolithophorid consumption (up to 26 and 23% of the coccolithophorid biomass, respectively).     

Age determination and growth rate of <Emphasis Type="Italic">Mactra chinensis</Emphasis> (Bivalvia: Mactridae) by external rings and chondrophore growth bands

Age, growth and mortality of Mactra chinensis were investigated during the period from October 2012 to September 2013 in Busan, South Korea. The monthly variation of the marginal index (MI) of the shell and chondrophore showed that the ring of this species was formed once a year during July. We estimated the age of M. chinensis by reading the external rings on the shell and the growth bands of the chondrophore to compare growth parameters between the two growth characters. The age of this species ranged from 0 to 8 years (shell-based age reading) and from 0 to 10 years (chondrophore-based age reading). Based on external rings and growth bands of chondrophore for the same period, the von Bertalanffy growth functions were expressed by the equation, L _t = 101.53[1-exp {-0.15(t + 0.75)}] and L _t = 90.03[1-exp {-0.20(t + 0.50)}], respectively. The likelihood test showed that there was a significant difference in L _∞ (P < 0.001), K (P < 0.001), to (P < 0.001) estimated from non-linear regression between the two growth characters.     

A new formula for saturated water steam pressure within the temperature range −25 to 220°C

Instead of approximation formula ln(E(t)/E(0)) = [(a ? bt)t/(c + T)] commonly used at present for representing dependence of pressure of saturated streams of liquid water E upon temperature we suggested new approximation formula of greater accuracy in the form ln(E(t)/E(0)) = [(A ? Bt + Ct ²)t/T], where t and T are temperature in °C and K respectively. For this formula with parameters A = 19.846, B = 8.97 × 10^?3, C = 1.248 × 10^?5 and E(0) = 6.1121 GPa with ITS-90 temperature scale and for temperature range from 0°C to 110°C relative difference of approximation applying six parameter formula by W. Wagner and A. Pruß 2002, developed for positive temperatures, is less than 0.005%, that is approximately 15 times less than accuracy obtained with the firs formula. Increase of temperature range results in relative difference increasing, but for even temperature range from 0°C to 220°C it does not higher than 0.1%. For negative temperatures relative difference between our formula and a formula of D. M. Murphy and T. Koop, 2005, is less than 0.1% for temperatures higher than ?25°C. This paper also presents values of coefficients for approximation of Goff and Grach formula recommended by IMO. The procedure of finding dew point T _d for known water steam pressure e _n based on our formula adds up to solving an algebraic equation of a third degree, which coefficients are presented in this paper. For simplifying this procedure this paper also includes approximation ratio applying a coefficient A noted above, in the form T _d (e _n ) = \(\frac{{AT_0 }}{{A - \varepsilon }}\) + 0.0866?² + 0.0116?^10/3, where ? = ln(e _n /E(T ₀)). Error of dew point recovery in this ratio is less than 0.005 K within the range from 0 to 50°C.     

Waves and diffusion on the sea surface

The frequency spectrum of surface elevations in the presence of wind waves is well known. On this basis, one can estimate the frequency spectrum of vertical velocities in sea-surface waves. Owing to liquid incompressibility, the spectrum of horizontal velocities should have the same frequency dependence. The use of the dispersion equation for waves on the surface of a heavy liquid allows one to obtain to the spatial spectrum of velocities. Therefore, one can estimate the spatial structure function of the velocity field. For short waves and large depths, the structure function increases as r ^1/2, where r is the distance between the points of observations. For long waves and shallow depths h, this increase is proportional to r. The coefficient of turbulent mixing K(r) of pollution spots of size r on the sea surface is now estimated as the product of the spot size and the rms difference of velocities. As a result, depending on r and h, the exponent in the r ⁿ dependence of K(r) may vary between 1.25 and 1.5. This outcome provides an explanation for a scatter in the values of the exponent n, a phenomenon that has been observed by many experimentalists.     

Variation of the photosynthetic electron transfer rate and electron requirement for daily net carbon fixation in Ariake Bay,Japan

Fast repetition rate fluorometry (FRRf) provides a potential means to examine marine primary productivity; however, FRRf-based productivity estimations require knowledge of the electron requirement (K) for carbon (C) uptake (K _C) to scale an electron transfer rate (ETR) to the CO₂ uptake rate. Most previous studies have derived K _C from parallel measurements of ETR and CO₂ uptake over relatively short incubations, with few from longer-term daily-integrated periods. Here we determined K _C by comparing depth-specific, daily ETRs and CO₂-uptake rates obtained from 24-h on-deck incubation experiments undertaken on seven cruises in Ariake Bay, Japan, from 2008 to 2010. The purpose of this study was to determine the extent of variability of K _C and to what extent this variability could be reconciled with the prevailing environmental conditions and ultimately to develop a method for determining net primary productivity (NPP) based on FRRf measurements. Both daily ETR and K _C of the upper layer varied considerably, from 0.5 to 115.7 mmol e^? mg Chl-a ^?1 day^?1 and 4.1–26.6 mol e^? (mol C)^?1, respectively, throughout the entire data set. Multivariate analysis revealed a strong correlation between daily photosynthetically active radiation (PAR) and K _C (r ² = 0.94). A simple PAR-dependent relationship derived from the data set was used for generating K _C, and this relationship was validated by comparing the FRRf-predicted NPP with the ¹³C uptake measured in 2007. These new observations demonstrate the potential application of FRRf for estimating regional NPP from ETR.     

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.     

High mortality and poor growth of green mussels, <Emphasis Type="Italic">Perna viridis</Emphasis>, in high chlorophyll-<Emphasis Type="Italic">a</Emphasis> environment

The current study was carried out from May 2014 to April 2015 to estimate the stock status of P. viridis in Marudu Bay. The gonad development was monitored by histological examination, while the population parameters including asymptotic length (L_∞), growth coefficient (K), mortality rate (Z, F and M), exploitation level (E) and recruitment of P. viridis were estimated using the lengthfrequency data. Results of the current study demonstrated that P. viridis in Marudu Bay spawned throughout the year with two major peaks, one in April to May and another one in October to December. The recruitment pattern was continuous with the peak in May to June 2014, which corresponded to the first spawning peak in April. However, no significant recruitment was observed from the second spawning peak due to the difference in spawning timing between male and female populations. The estimated asymptotic length (L_∞), growth coefficient (K), total mortality (Z), natural mortality (M), fishing mortality (F) and growth performance (ф) of P. viridis in Marudu Bay were estimate to be 117 mm, 0.97 yr^-1, 4.39 yr^-1, 1.23 yr^-1, 3.16 yr^-1 and 4.123, respectively. The exponent b of the lengthweight relationship was 2.4 and exploitation level (E) was 0.72. The high mortality, low condition indices and negative allometric of P. viridis in Marudu Bay is caused by a lack of suitable food in the surrounding water.     

Estimating the parameters of pulsed sources from data on acoustic waves recorded in the atmosphere

The possibility of estimating the parameters of surface pulsed sources from data on acoustic waves recorded in the atmosphere is studied. Experimental values are given for peak pressure P ₊ of recorded acoustic signals, wave-profile area S ₊ in their positive phase, and length t ₊ of this phase, and the approximations of these parameters are obtained within wide ranges of source energy 10^–3 < E < 10¹⁰ kg TNT and scaled distances 1 < R/E ^1/3 < 4 × 10⁴ m/kg^1/3. Conventional methods of estimating the acoustic energy E according to data obtained from acoustic measurements in the atmosphere are analyzed, and ways to improve their accuracy are proposed. The influence of the type of explosions on the parameters P ₊, S ₊, and t ₊ of acoustic signals at long distances R/E ^1/3 > 500 m/kg^1/3 from explosions is shown.     

Oceanic processes of upper ocean heat content associated with two types of ENSO

The spatiotemporal variability of equatorial Pacific upper ocean heat content (HC) and subsurface heat during two types of El Niño-Southern Oscillation (ENSO), namely eastern and central Pacific (EP and CP) types, is investigated using subsurface ocean heat budget analysis. Results show that HC tendencies during both types of ENSO are mainly controlled by oceanic heat advection beneath the mixed layer to the thermocline, and the role of net surface heat flux can be neglected. The most important three terms are the zonal and vertical advections of anomalous heat by climatological currents (QU ₀ T′, QW ₀ T′) and zonal advection of climatological heat by anomalous current (QU′T ₀). The large contribution of QU ₀ T′ extends from west to east along the equatorial Pacific. The considerable contribution of QU′T ₀ is confined to the east of 160°W, and that of the QW ₀ T′ is observed in the central Pacific between 180°E and 120°W. In particular, a major contribution of QW ₀ T′ is also observed in the far eastern Pacific east of 100°W during EP ENSO. There is also a small contribution from meridional advection of climatological heat by anomalous current (QV′T ₀). In contrast, the meridional advection of anomalous heat by climatological currents (QV ₀ T′) and vertical advection of climatological heat by anomalous current (QW′T ₀) are two damping factors in the HC tendency, with the former dominating. Differences in spatial distribution of the heat advection associated with the two types of ENSO are also presented. We define a warm water heat index (WWH) as integrated heat content above 26 kg m^?3 potential density (26σ _? ) isopycnal depth within 130°E–80°W and 5°S–5°N. Further examination suggests that the recharge–discharge of WWH is involved in both types of El Niño, though with some differences. First, it takes about 42 (55) months for the evolution of a recharge–discharge cycle during an EP (CP) ENSO. Second, the EP El Niño event peaks during the discharge phase, 7–8 months after the recharge time. The CP El Niño peaks during the recharge phase, 4–5 months before the recharge time. The locations of HC anomalies in the El Niño mature phase relative to those at recharged time explain why the EP and CP El Niño peak in different stages of the recharge–discharge process.     

黄海小黄鱼生物学特征对多重压力的响应

Temporal changes in biological characteristics of small yellow croaker Larimichthys polyactis in the Yellow Sea were examined for the period of 1960–2008. The body size and age of small yellow croaker decreased substantially, in particular, average length of fish in 2008 was reduced by ~85% than those occurring in 1985, and at that time ~93% of the total catch was dominated by one-year-old individuals. Correspondingly, growth parameters also varied significantly over the years, i.e., k(growth coefficient) and t_0(zero-length age) gradually increased from 0.26 and –0.58 year in 1960 to 0.56 and –0.25 year in 2008, respectively. Although, L∞(body length)sharply decreased from 34.21 cm in 1960 to 24.06 cm in 2008, and t_r(inflexion age) decreased from 3.78 year in1960 to 1.61 year in 2008. There was a great increase both in natural mortality coefficient and fishing mortality coefficient. However, according to the gray correlation analysis, changes in the biological characteristics of small yellow croaker were induced by different stressors ranked as: fishing vessel powerfeeding gradesea surface temperature. This study suggests that the active fishery management measures for biological characters of fish populations should be considered.