Inter-Annual Variability of the Seawater Light Absorption in Surface Layer of the Northeastern Black Sea in Connection with Hydrometeorological Factors

Data on the light absorption by seawater and its components are needed in many theoretical and practical aspects of marine science and engineering. However, up to now, there is a lack of such data for the northeastern part of the Black Sea. This article presents the data on light absorption measured by a portable integrated cavity absorption meter (ICAM) spectrophotometer in the Gelendzhik region of the Black Sea during field studies in June 2017 and 2018, together with other bio-optical and oceanographic data from in situ measurements and satellite observations. In 2018, the elevated values of the colored dissolved organic matter (CDOM) absorption in the surface layer were observed concurrently with high values of salinity, contradicting the idea of river runoff being the main CDOM source. The vertical profiles of salinity differed in 2017 and 2018, especially in shallow waters; in the upper layer, the salinity increased from 17.1 psu in 2017 to 17.8 psu in 2018, while the values of CDOM absorption increased from 0.10 to 0.16 m−1. The analysis of available hydrometeorological data pointed to intensive vertical mixing due to the strong wind forcing as a main factor in increasing values of both salinity and the CDOM absorption in the surface layer in 2018.


Introduction
This work is aimed at studying the inter-annual variability of the light absorption properties of seawater in the surface layer of the northeastern part of the Black Sea. Absorption, along with scattering, is a basic process governing the light radiation propagation in water from both natural and artificial sources; the absorption parameters are included (directly or indirectly) in the radiation transfer equation and most approximations. Knowledge of the absorption coefficients is needed to assess both spectral transmission and reflectance of a given water layer or the whole water column.
The seawater absorption properties depend on the content of optically active components of seawater, such as the colored dissolved organic matter (CDOM), phytoplankton pigments, detritus and other particles, and the seawater itself. The content of these components varies depending on water types, as well as seasonally and inter-annually. The spectral values of seawater absorption coefficient are very sensitive to variability of these components and reflect their changes. Specialists in different fields need such data in their studies: marine biologists for assessment of the photosynthetically available radiation (PAR) at different depths and assimilation of solar radiation for primary production; physical oceanographers for estimating ocean albedo and volume absorption in water column; marine ecologists for tracking river runoff and propagation of various impurities.
Satellite observations focused on the eastern open part and the eastern and southern shelves, which correspond to sub-regions #7 and #8 (see Figure 1) according to the demarcation below [19,20].  [19,20]; (b) location of the drift stations in 2017 (blue) and 2018 (red); the blue line is the border between sub-regions #7 and #8.
The analysis of the obtained data confirmed the assumption that it is possible to use ag for estimating the contribution of particles brought by river runoff to the particle backscattering coefficient bbp and for estimating the terrigenous component TSM_trg of the suspended matter. According to data from 2005-2007, the coefficient of determination between TSM_trg and ag was 0.86 [14]. In the coastal zone, the coefficient of determination between bbp and ag was 0.82, whereas, in the open part, it was close to 0. The average value of ag in the open part was 0.098 m −1 as compared to 0.13 m −1 in the coastal zone. As for the background values, the lowest monthly mean of ag, derived from satellite data over 2003-2010, was equal to 0.047 m −1 with the standard deviation of 0.013 m −1 (range: 0.021-0.060 m −1 ) [14].
Direct measurements of the spectral absorption coefficients ag(λ) were conducted since 2017 with a portable spectrophotometer in the integrated cavity absorption meter (ICAM) configuration, developed at the Department of Biophysics, Biological Faculty, Moscow State University [21]. Its description is given in the next section. Measurements were carried out on samples taken from different horizons at drift stations. The stations were located at depths of 25-100 m in the shallow-water part and 500-1500 m in the deep-water part ( Figure 1b); for more details, see Reference [22].
In this paper, we focus on the results of 2017-2018; the analysis is based on both satellite and in situ data (see the next section). These years are quite different in their oceanographic and hydrometeorological conditions; in particular, in 2017, the most intensive CB within the whole period of studies was observed, whereas there was no pronounced CB in 2018. The research highlight of our work is an explanation of the unusual situation concerning the CDOM absorption in 2018, namely, high CDOM absorption under elevated values of salinity. This means that the river runoff is not a main source of CDOM as mentioned above, and an alternative source should be found. To be sure of the appropriate quality of the measurement data, we compared the values of ag derived from different measurements. To obtain a more complete picture of the influence of various factors, the available data from previous years were also analyzed.  [19,20]; (b) location of the drift stations in 2017 (blue) and 2018 (red); the blue line is the border between sub-regions #7 and #8.

Measurements of the Spectral Absorption Parameters
Direct measurements of the spectral absorption coefficients a g (λ) were conducted since 2017 with a portable spectrophotometer in the integrated cavity absorption meter (ICAM) configuration, developed at the Department of Biophysics, Biological Faculty, Moscow State University [21]. Its description is given in the next section. Measurements were carried out on samples taken from different horizons at drift stations. The stations were located at depths of 25-100 m in the shallow-water part and 500-1500 m in the deep-water part ( Figure 1b); for more details, see Reference [22].
In this paper, we focus on the results of 2017-2018; the analysis is based on both satellite and in situ data (see the next section). These years are quite different in their oceanographic and hydrometeorological conditions; in particular, in 2017, the most intensive CB within the whole period of studies was observed, whereas there was no pronounced CB in 2018. The research highlight of our work is an explanation of the unusual situation concerning the CDOM absorption in 2018, namely, high CDOM absorption under elevated values of salinity. This means that the river runoff is not a main source of CDOM as mentioned above, and an alternative source should be found. To be sure of the appropriate quality of the measurement data, we compared the values of a g derived from different measurements. To obtain a more complete picture of the influence of various factors, the available data from previous years were also analyzed.

Measurements of the Spectral Absorption Parameters
Spectral seawater absorption coefficients were obtained using a portable spectrophotometer in the ICAM configuration, developed at the Department of Biophysics, Biological Faculty, Moscow State University, Moscow, Russia [21].
The so-called ICAM (integrated cavity absorption meter) technologies, named so because the studied water is placed inside the integrating sphere, were actively developed in the past two decades. The use of such an approach allows avoiding problems associated with light scattering (see Section 1) and increasing sensitivity due to the multiple reflections of light inside a sphere. However, for the determination of the absolute values of the absorption coefficient, it is necessary to know the effective path length of photons, taking into account their multiple reflections. This problem was studied, and a special algorithm with a processing program was developed. The Monte Carlo calculations performed showed that the results of the determination of the absorption coefficient are independent of the seawater scattering properties when the scattering coefficient b varies from 0 to 5 m −1 [23].
Three measurements of the absorption spectrum were taken for each water sample: (1) seawater; (2) after filtration through the nuclear filter of 0.4 µm (absorption spectra of the filtrate); (3) the empty sphere. The spectral absorption coefficient of seawater (or the filtrate) was calculated through the ratio of the measured spectral dependencies for the sphere filled with the sample to the empty sphere [23]. Distilled water spectra were recorded in the control tests. The volume of the water sample was 330 mL; replicate samples were collected and measured for each series of measurements. The absorption spectra of the particles were not measured but determined as a difference between the spectral values of the seawater absorption and of the filtrate absorption. A detailed description of the measurement procedure and calculation is given in Reference [23].

Measurements of the Vertical Structure
A conductivity, temperature, and depth (CTD) profiler with the "SBE 19+, Seacat" sampling system was used to measure the vertical distributions of the following hydrological parameters: salinity, temperature, density, and chlorophyll fluorescence (F_chl). To measure the vertical distribution of the light attenuation coefficient, a submersible transparency meter [24] was used.
In accordance with the features of the vertical structure, water samples at each station were taken from 4-5 selected depths for study in an onshore laboratory. The laboratory measurements included the absorption spectra by the ICAM instrument, determination of concentrations of chlorophyll, suspended matter, phytoplankton content, and species composition. These studies were performed according to standard procedures [14].
Direct measurements of the depth dependence of spectral underwater irradiance E d (z, λ) created by the downwelling flux of solar radiation were performed using a RAMSES hyperspectral radiometer (model ACC-VIS), manufactured by TriOS Optical Sensors, Germany. The radiometer is designed to instantaneously measure spectral irradiance at wavelengths from 320 to 950 nm with a 3-nm resolution down to a depth of 100 m. The dynamic range of the irradiance channel is 3-0.003 W·m −2 ·nm −1 . A special frame was constructed for properly submersing the instrument; the data were transmitted via cable to a laboratory unit. A notebook computer was used to record and process the data. The radiometer was calibrated by the TriOS company in compliance with National Institute of Standards and Technology (NIST) specifications. The measured E d (z, λ) values allowed us to compute the diffuse attenuation coefficients for underwater irradiance, K d (λ) under the assumption that irradiance attenuation with depth obeys the exponential law E d (z) = E d (0 − ) exp[−K d z], where wavelength λ is omitted for simplicity.

Calculation of Hydrometeorological Parameters
Calculations of precipitation and wind characteristics were carried out on the basis of data from the ERA-Interim (European Interim Reanalysis) reanalysis archive using the ECMWF model (European Center for Medium-Range Weather Forecasts) with a 12-hour time step and a spatial resolution of 0.125 × 0.125 degrees, available at the website [25]. These models were used to calculate daily spatial distributions from 1 April to 11 June (the expedition end date) To assess the river discharge in this part of the Black Sea, data from an automated flood situation monitoring system in the Krasnodar region, available at the website [26], were used. Unfortunately, most of the posts provided only data on the water level in the rivers in the Baltic height system. We used data for the Ashamba River (Yashamba), the closest to the area of the field studies, flowing into Rybachya (Blue) Bay, and Mzymta, a deep river in the region under consideration, fasting station AGK-0160 (43.517 • N, 39.996 • E). The measurement data were obtained with an interval of 10 min; in constructing the graphs, the running average method with an averaging interval of one day was used.

Satellite Observation Data
On the National Aeronautics and Space Administration (NASA) website [27], the parameters of our interest are freely available, but most of them are calculated using empirical (regression) or semi-analytical algorithms [19,28]. Most of these algorithms were developed by using data from field measurements mainly for ocean waters. In waters that are strongly influenced by river runoff, such as the waters of the Black Sea, the standard algorithms can produce large errors [19]. Thus, in such seas, it is necessary to use regional algorithms to obtain correct data on bio-optical parameters. These algorithms were developed on the basis of data from the field measurements in the studied regions, taking into account their features. Validation studies performed in SIO RAS expeditions showed that the use of regional algorithms can significantly reduce errors of the satellite estimation, compared to standard algorithms.
The algorithms used in this work are described in detail in Reference [19] and on the website [20]. They allow to calculate the chlorophyll concentration Chl, the particle backscattering coefficient b bp , and the absorption coefficient of yellow substance a g ; the special regional algorithms were developed to derive concentrations of the chlorophyll, suspended matter, and coccolithophore cells during the mass blooming period which usually occurs in June [14,19,20].
We used the Level 2 data of satellite spectroradiometers Moderate-Resolution Imaging Spectroradiometer (MODIS)-Aqua and MODIS-Terra, available through the NASA website [27]. These are the values of geophysical parameters calculated after atmospheric correction, corresponding to the initial pixels, i.e., with the same coordinate and time values as the original values. The five-day averages were calculated for the spatial distribution of the values of a g , Chl, and TSM in May and June 2017-2018. The monthly average values of these characteristics and sea surface temperature were estimated for the period of 2012-2018. The calculation was performed with a spatial resolution of 2 × 2 km for five-day averages and 3 × 3 km for monthly averages by using software [29] developed in the Laboratory of Ocean Optics SIO RAS, Moscow, Russia. Data marked as land, clouds, and stray light flags were excluded.
The above-mentioned package includes programs for maintaining the database, which contains files from satellite and ship measurements, programs for batch processing of large amounts of satellite data, and a host program for visualizing satellite and ship data, as well as for launching other components of the program complex. A set of programs for batch processing of satellite files contain modules for calculating new products and recording results in HDF (Hierarchical Data Format), and NetCDF (Network Common Data Form) formats, averaging data, constructing maps in a given projection, creating time series, etc. Features of the programs of this system are the dynamic loading of libraries that can be developed by users and implemented without recompilation of programs, and the presence of a convenient and fairly universal language for task management [19]. Figure 2 shows the changes in monthly mean values of the yellow substance absorption coefficients a g and the particle backscattering coefficients b bp in the #7 and #8 sub-regions (see Figure 1) derived from MODIS-Aqua data of 2012-2018 by using the SIO RAS regional algorithms. The black lines show the seasonal variability of the "climatic" values which are averaged values over the period of 1998-2011. One can see the pronounced inter-annual variability of the bio-optical parameters compared with the seasonal variability of the climatic values. Our purpose is to identify the causes of the observed variability, primarily for the yellow substance absorption coefficients a g . The solution to this problem is possible only on the basis of in situ measured data and an analysis of factors that determine the seasonal and inter-annual variability of bio-optical and hydrometeorological characteristics.

Satellite Data
period of 1998-2011. One can see the pronounced inter-annual variability of the bio-optical parameters compared with the seasonal variability of the climatic values. Our purpose is to identify the causes of the observed variability, primarily for the yellow substance absorption coefficients ag. The solution to this problem is possible only on the basis of in situ measured data and an analysis of factors that determine the seasonal and inter-annual variability of bio-optical and hydrometeorological characteristics. We are mainly interested in data for June when our field studies were carried out, as well as, of course, the above area itself. Thus, we limited the representation of satellite data mainly to this area (see Figure 1) and considered the dynamics of the changes during the period when these studies were conducted. Figure 2 clearly demonstrates that the seasonal and inter-annual changes in the monthly means of ag and bbp are quite different. In particular, the June monthly means of ag were larger than the climatic values in 2015, 2016, and 2018, whereas the bbp values in these years were lower. This certainly indicates that the variability of ag and bbp is controlled by different factors. Figure 3 shows the spatial distributions of the yellow substance absorption coefficient ag and the particle backscattering coefficient bbp from MODIS data averaged over five days from 1 June to 15 June 2017 and 2018. One can see that the years 2017 and 2018 sharply differed from each other. In 2017, high values of bbp but low values of the yellow substance absorption coefficient ag were observed. In 2018, on the contrary, markedly higher values of the yellow substance absorption coefficients were observed, especially near the coast, along with lower bbp values.
The reason for elevated values of the particle backscattering coefficients bbp is known; CB of very high intensity was observed in June 2017. Figure 4 shows the spatial distributions of the coccolithophore concentration, calculated from MODIS-Aqua data on 5 June and 8 June 2017, using the SIO RAS algorithm [14]. As seen, in both days, the concentration of cells in the coastal zone in the area of our work exceeded 8 × 10 6 cells/L. We are mainly interested in data for June when our field studies were carried out, as well as, of course, the above area itself. Thus, we limited the representation of satellite data mainly to this area (see Figure 1) and considered the dynamics of the changes during the period when these studies were conducted. Figure 2 clearly demonstrates that the seasonal and inter-annual changes in the monthly means of a g and b bp are quite different. In particular, the June monthly means of a g were larger than the climatic values in 2015, 2016, and 2018, whereas the b bp values in these years were lower. This certainly indicates that the variability of a g and b bp is controlled by different factors. Figure 3 shows the spatial distributions of the yellow substance absorption coefficient a g and the particle backscattering coefficient b bp from MODIS data averaged over five days from 1 June to 15 June 2017 and 2018. One can see that the years 2017 and 2018 sharply differed from each other. In 2017, high values of b bp but low values of the yellow substance absorption coefficient a g were observed.
In 2018, on the contrary, markedly higher values of the yellow substance absorption coefficients were observed, especially near the coast, along with lower b bp values.
The reason for elevated values of the particle backscattering coefficients b bp is known; CB of very high intensity was observed in June 2017. Figure 4 shows the spatial distributions of the coccolithophore concentration, calculated from MODIS-Aqua data on 5 June and 8 June 2017, using the SIO RAS algorithm [14]. As seen, in both days, the concentration of cells in the coastal zone in the area of our work exceeded 8 × 10 6 cells/L.
The difference in the values of the particle backscattering coefficients b bp in 2017 and 2018 was strongly pronounced in the satellite observation under clear sky conditions. It is clearly seen from the comparison between the values of the spectral remote-sensing reflectance R rs (λ) measured by the MODIS-Aqua sensor in 2017 and 2018 ( Figure 5). It is also seen that the difference between shallow-water and deep-water stations was rather small in 2018 and higher in 2017.
On 8 June, the full-scale in situ measurements were carried out and the samples were taken to determine the quantitative and species composition of phytoplankton. The coccolithophore cell concentration in the coastal zone turned out to be approximately the same as that from the satellite data, i.e., about 8 × 10 6 cells/L.   The difference in the values of the particle backscattering coefficients bbp in 2017 and 2018 was strongly pronounced in the satellite observation under clear sky conditions. It is clearly seen from the comparison between the values of the spectral remote-sensing reflectance Rrs (λ) measured by the MODIS-Aqua sensor in 2017 and 2018 ( Figure 5). It is also seen that the difference between shallow-water and deep-water stations was rather small in 2018 and higher in 2017. On 8 June, the full-scale in situ measurements were carried out and the samples were taken to determine the quantitative and species composition of phytoplankton. The coccolithophore cell concentration in the coastal zone turned out to be approximately the same as that from the satellite data, i.e., about 8 × 10 6 cells/L.

Spectral Absorption Measurements on the ICAM
In Section 1, we mentioned an unusual situation concerning the CDOM absorption in 2018, namely, high CDOM absorption under elevated values of salinity. Here, we try to show this situation in detail and understand which additional information is required to explain the results obtained.   On 8 June, the full-scale in situ measurements were carried out and the samples were taken to determine the quantitative and species composition of phytoplankton. The coccolithophore cell concentration in the coastal zone turned out to be approximately the same as that from the satellite data, i.e., about 8 × 10 6 cells/L.

Spectral Absorption Measurements on the ICAM
In Section 1, we mentioned an unusual situation concerning the CDOM absorption in 2018, namely, high CDOM absorption under elevated values of salinity. Here, we try to show this situation in detail and understand which additional information is required to explain the results obtained.

Spectral Absorption Measurements on the ICAM
In Section 1, we mentioned an unusual situation concerning the CDOM absorption in 2018, namely, high CDOM absorption under elevated values of salinity. Here, we try to show this situation in detail and understand which additional information is required to explain the results obtained. Figure 6a shows the absorption spectra of the filtrate (a f ), and Figure 6b shows the suspended particles (a p ), measured in the surface layer at three shallow-water stations in 2017, 2018, and 2019. As seen, the highest absorption by the filtrate was observed in 2018, and by suspended particles in 2017, where the latter is obviously associated with intensive CB.
Recall that the absorption spectra of particles were obtained as a difference between the spectral values of the absorption coefficient of seawater before and after filtration through a filter with a pore size of 0.4 µm. As seen from Figure 6b, the absorption spectra of particles in all years expose the maxima caused by phytoplankton pigments: the blue maximum at 430-440 nm, the shoulder near 480 nm due to absorption by carotenoids, and the red maximum at 675-680 nm caused by chlorophyll absorption. This suggests that most suspended particles, the sizes of which exceed 0.4 µm, are of biological origin. Figure 6a shows the absorption spectra of the filtrate (af), and Figure 6b shows the suspended particles (ap), measured in the surface layer at three shallow-water stations in 2017, 2018, and 2019. As seen, the highest absorption by the filtrate was observed in 2018, and by suspended particles in 2017, where the latter is obviously associated with intensive CB. Recall that the absorption spectra of particles were obtained as a difference between the spectral values of the absorption coefficient of seawater before and after filtration through a filter with a pore size of 0.4 μm. As seen from Figure 6b, the absorption spectra of particles in all years expose the maxima caused by phytoplankton pigments: the blue maximum at 430-440 nm, the shoulder near 480 nm due to absorption by carotenoids, and the red maximum at 675-680 nm caused by chlorophyll absorption. This suggests that most suspended particles, the sizes of which exceed 0.4 μm, are of biological origin.
The average values of contributions to the seawater light absorption from the CDOM and suspended particles at the wavelengths corresponding to the spectral channels of MODIS in the layers at 0-5 m and 6-35 m are shown in Figure 7. It is seen that the highest values of the seawater absorption coefficients were recorded in 2018, due to the increased CDOM absorption. The higher (about twofold) values of the particle absorption coefficient ap were recorded in 2017 under conditions of strong CB. In 2019, absorption values were close to those observed in 2017, but with less particulate absorption than in 2017 and less CDOM absorption than in 2018. The inter-annual differences are also observed in the underlying layers, even more pronounced than in the surface layer.
One should pay attention to absorption in the underlying layer of 6-35 m; in all years, it was noticeably higher than in the subsurface layer of 0-5 m. Next, we turn to this observation.

Comparison between the ag values derived from different measurements
Comparison of the values of the yellow substance absorption coefficient measured by the ICAM and calculated from MODIS-Aqua and MODIS-Terra satellite data using the SIO RAS algorithm showed good agreement between them (Figure 8): coefficient of determination r 2 = 0.63 (correlation The average values of contributions to the seawater light absorption from the CDOM and suspended particles at the wavelengths corresponding to the spectral channels of MODIS in the layers at 0-5 m and 6-35 m are shown in Figure 7. It is seen that the highest values of the seawater absorption coefficients were recorded in 2018, due to the increased CDOM absorption. The higher (about twofold) values of the particle absorption coefficient a p were recorded in 2017 under conditions of strong CB. In 2019, absorption values were close to those observed in 2017, but with less particulate absorption than in 2017 and less CDOM absorption than in 2018. The inter-annual differences are also observed in the underlying layers, even more pronounced than in the surface layer.
One should pay attention to absorption in the underlying layer of 6-35 m; in all years, it was noticeably higher than in the subsurface layer of 0-5 m. Next, we turn to this observation.  Table 1 presents results of the comparison between the values of ag(443) measured by ICAM and those calculated using MODIS data on the remote-sensing reflectance Rrs and from data of Kd measurements in 2018.

Comparison between the a g values derived from different measurements
The spectral diffuse attenuation coefficient Kd in the near-surface layer (0-10 m) was calculated from the measured values of the downwelling spectral irradiance. Furthermore, the sum of the absorption and backscattering coefficients was computed using Gordon's formula [30]. The particle backscattering coefficient bbp(555) was determined from satellite data using the algorithm in Reference [31]; then, the spectral particle backscattering coefficient was calculated using the following formula: bbp(λ) = bbp(555) (555/λ) 0.5 . The data in Reference [32] were used for backscattering by pure seawater. As a result, the estimates for the spectral coefficient of seawater absorption were obtained, which are in good agreement with the measurement data. For 2017, no calculations were carried out, since Gordon's formula gives big errors for conditions of strong CB (this formula was derived with quasi-single scattering approximation [30]).  Table 1 presents results of the comparison between the values of a g (443) measured by ICAM and those calculated using MODIS data on the remote-sensing reflectance R rs and from data of K d measurements in 2018. The spectral diffuse attenuation coefficient K d in the near-surface layer (0-10 m) was calculated from the measured values of the downwelling spectral irradiance. Furthermore, the sum of the absorption and backscattering coefficients was computed using Gordon's formula [30]. The particle backscattering coefficient b bp (555) was determined from satellite data using the algorithm in Reference [31]; then, the spectral particle backscattering coefficient was calculated using the following formula: b bp (λ) = b bp (555) (555/λ) 0.5 . The data in Reference [32] were used for backscattering by pure seawater. As a result, the estimates for the spectral coefficient of seawater absorption were obtained, which are in good agreement with the measurement data. For 2017, no calculations were carried out, since Gordon's formula gives big errors for conditions of strong CB (this formula was derived with quasi-single scattering approximation [30]).
From Table 1, it is seen that the K d measurements give smaller errors than those from MODIS data. In the first case, the difference is within 9-22%; in the second, the error varies from −6% to +30%. Table 2 presents the average values of the main parameters measured in situ in 2017 and 2018 and averaged over three layers of different depths; these are the yellow substance absorption coefficient, salinity, temperature, chlorophyll concentration, and the beam attenuation coefficient c for shallow-water (station depth 25-80 m) and deep-water stations (more than 300 m), averaged over three layers. The values of the yellow substance absorption coefficient a g and the salinity S in the subsurface layer are highlighted in red. It is seen that, in 2018, the selected values of a g were larger than in 2017 both in the shelf zone and in the open sea. It is noteworthy that the values of salinity were also greater. This is an unusual situation because it is generally accepted that the CDOM concentration is associated with river runoff, and its increase must be accompanied by a decrease in salinity. Figure 9 shows the daily averaged profiles of salinity, temperature, and chlorophyll from the CTD fluorescence sensor, and the beam attenuation coefficient from a transparency meter, separately for shallow-water and deep-water stations in 2017 and 2018. thermocline, like the maximum fluorescence, was deeper in 2018 than in 2017. There is a sharp distinction observed between the vertical salinity profiles in 2017 and 2018 both in shallow-water and deep-water stations. If, in 2017, the differences in salinity on the surface and below 15 m were about 1 psu, then, in 2018, these differences were only about 0.4 psu. It can be assumed that the observed differences are associated with different hydrometeorological conditions in 2017 and 2018, and, to test this assumption, we examined precipitation, wind, and river levels. , and, to test this assumption, we examined precipitation, wind, and river levels.

Discussion
The results of our field studies in 2017 and 2018 with the involvement of satellite data revealed an unusual situation in 2018 as compared with 2017, namely, the elevated values of the CDOM absorption coefficient a g concurrently with higher values of salinity in the surface layer. This is in contradiction with the accepted idea that the main source of CDOM is the river runoff and that, therefore, the increased values of the CDOM absorption should be accompanied by lower values of salinity. In this section, the various factors are discussed to understand which of them can potentially affect the situation.

Interannual Salinity Measurements
Variability of salinity in the upper layer in Gelendzhik region of the Black Sea during the warm season in 2010-2013 was studied, together with analysis of coastal precipitation and wind forcing, by Podymov and Zatsepin [22].

Discussion
The results of our field studies in 2017 and 2018 with the involvement of satellite data revealed an unusual situation in 2018 as compared with 2017, namely, the elevated values of the CDOM absorption coefficient ag concurrently with higher values of salinity in the surface layer. This is in contradiction with the accepted idea that the main source of CDOM is the river runoff and that, therefore, the increased values of the CDOM absorption should be accompanied by lower values of salinity. In this section, the various factors are discussed to understand which of them can potentially affect the situation.

Interannual Salinity Measurements
Variability of salinity in the upper layer in Gelendzhik region of the Black Sea during the warm season in 2010-2013 was studied, together with analysis of coastal precipitation and wind forcing, by Podymov and Zatsepin [22].  It is seen from Table 2 and Figure 9 in the previous section that, in 2018, the salinity in the upper 5-m layer in the shelf zone was even higher than in 2013-about 17.8 psu (in 2017, it was slightly higher than in 2011). Looking ahead, we note that precipitation in 2018 was lower than in 2013 (from 1 April to 11 June, 0.68 mm and 0.83 mm, respectively); in 2017, it was also lower than in 2011-1.89 and 2.00 mm, respectively (Table 3). It is seen from Table 2 and Figure 9 in the previous section that, in 2018, the salinity in the upper 5-m layer in the shelf zone was even higher than in 2013-about 17.8 psu (in 2017, it was slightly higher than in 2011). Looking ahead, we note that precipitation in 2018 was lower than in 2013 (from 1 April to 11 June, 0.68 mm and 0.83 mm, respectively); in 2017, it was also lower than in 2011-1.89 and 2.00 mm, respectively (Table 3).
According to the results obtained in [22], high salinity values in June 2013 are primarily associated with precipitation or, rather, their absence; in 2013, the total precipitation at the beginning of May was 1.5 times less than in 2010-2011 (214 mm), and this trend continued in June. According to Reference [22], the effect of wind variability on surface salinity was not noticeable; however, generally speaking, this effect can be significant, since its impact is proportional to the third power of wind speed. At a depth of 30 m, salinity was approximately 18.2 psu in both 2017 and 2018 (see Figure 9)  Figure 11 shows the daily precipitation from 1 April to 11 June (the end of the expedition) 2010-2018 and the monthly mean precipitation, summarized over the catchment area (see Figure 12); numerical values are given in Table 3. These values were obtained from data of the ERA-Interim reanalysis archive (European Interim Reanalysis), using the ECMWF model (European Center for Medium-Range Weather Forecasts).  According to the results obtained in [22], high salinity values in June 2013 are primarily associated with precipitation or, rather, their absence; in 2013, the total precipitation at the beginning of May was 1.5 times less than in 2010-2011 (214 mm), and this trend continued in June. According to Reference [22], the effect of wind variability on surface salinity was not noticeable; however, generally speaking, this effect can be significant, since its impact is proportional to the third power of wind speed.

Precipitation and Wind
At a depth of 30 m, salinity was approximately 18.2 psu in both 2017 and 2018 (see Figure 9)  Figure 11 shows the daily precipitation from 1 April to 11 June (the end of the expedition) 2010-2018 and the monthly mean precipitation, summarized over the catchment area (see Figure 12); numerical values are given in Table 3. These values were obtained from data of the ERA-Interim reanalysis archive (European Interim Reanalysis), using the ECMWF model (European Center for Medium-Range Weather Forecasts).    According to the results obtained in [22], high salinity values in June 2013 are primarily associated with precipitation or, rather, their absence; in 2013, the total precipitation at the beginning of May was 1.5 times less than in 2010-2011 (214 mm), and this trend continued in June. According to Reference [22], the effect of wind variability on surface salinity was not noticeable; however, generally speaking, this effect can be significant, since its impact is proportional to the third power of wind speed.

Precipitation and Wind
At a depth of 30 m, salinity was approximately 18.2 psu in both 2017 and 2018 (see Figure 9), but there are noticeable differences in the upper 15-m layer between 2017 and 2018. Pay attention to the change in the lower depth of the halocline; in 2017, it was about 13 m; in 2018, it increased to 15 m. Figure 11 shows the daily precipitation from 1 April to 11 June (the end of the expedition) 2010-2018 and the monthly mean precipitation, summarized over the catchment area (see Figure 12); numerical values are given in Table 3. These values were obtained from data of the ERA-Interim reanalysis archive (European Interim Reanalysis), using the ECMWF model (European Center for Medium-Range Weather Forecasts).     Comparing Figures 11 and 13, one can see that the maxima of the water level in the rivers often correspond to the maximum precipitation. For example, in 2017 three well-seen maxima of precipitation (9,16,and 23 May) are also displayed on the graph of changes in water level in the Mzymta River. For the Ashamba River during this period, an increase in water level is seen only on 16-17 May, which is consistent with the distribution of precipitation in the region.

Main Results and Conclusion
In 2017 and 2018, comprehensive bio-optical studies were carried out in the northeastern part of the Black Sea (near Gelendzhik), including direct measurements of the seawater spectral absorption coefficients by means of a portable spectrophotometer. An unusual situation was observed in 2018, when elevated values of the CDOM absorption in the surface layer were accompanied by a high Comparing Figures 11 and 13, one can see that the maxima of the water level in the rivers often correspond to the maximum precipitation. For example, in 2017 three well-seen maxima of precipitation (9, 16, and 23 May) are also displayed on the graph of changes in water level in the Mzymta River. For the Ashamba River during this period, an increase in water level is seen only on 16-17 May, which is consistent with the distribution of precipitation in the region. Figure 14 presents the values of the third degree of wind speed according to ECMWF ERA-Interim reanalysis, averaged over the area at a height of 10 m, which is proportional to the magnitude of wind exposure. The average wind speed from 1 May to 11 June was 3.5 m/s in both 2017 and 2018, but the highest value (10.7 m/s) was observed on 29 May 2018 ( Figure 14) and the total wind exposure in 2018 was significantly higher than in 2017-more than 4000 m 3 /s 3 and less than 3000 m 3 /s 3 , respectively. The weakest wind effect was observed in 2011-about 1650 m 3 /s 3 . by an average of 20-30 cm, which corresponds to the estimates of the total precipitation for the area under consideration, which was more in 2017. Comparing Figures 11 and 13, one can see that the maxima of the water level in the rivers often correspond to the maximum precipitation. For example, in 2017 three well-seen maxima of precipitation (9,16,and 23 May) are also displayed on the graph of changes in water level in the Mzymta River. For the Ashamba River during this period, an increase in water level is seen only on 16-17 May, which is consistent with the distribution of precipitation in the region.

Main Results and Conclusion
In 2017 and 2018, comprehensive bio-optical studies were carried out in the northeastern part of the Black Sea (near Gelendzhik), including direct measurements of the seawater spectral absorption coefficients by means of a portable spectrophotometer. An unusual situation was observed in 2018, when elevated values of the CDOM absorption in the surface layer were accompanied by a high

Main Results and Conclusion
In 2017 and 2018, comprehensive bio-optical studies were carried out in the northeastern part of the Black Sea (near Gelendzhik), including direct measurements of the seawater spectral absorption coefficients by means of a portable spectrophotometer. An unusual situation was observed in 2018, when elevated values of the CDOM absorption in the surface layer were accompanied by a high value of salinity. This means that the river runoff is not the only source of yellow substance in the surface layer, as it is considered at present.
Comparison between vertical profiles of salinity in 2017 and 2018 found not only an increase in salinity in the upper 5-m layer but also a change in the shape of a vertical profile-a violation of the pronounced vertical stratification, which was in 2017, and erosion of the halocline, lowering its lower boundary. If the increase in salinity in the upper 5-m layer can be explained by a decrease in precipitation in 2018 (the smallest amount from 1 April to 11 June in the period 2010-2018), the change in underlying layer suggests that it was due to vertical mixing.
A thorough analysis of changes in salinity in a given area under the influence of various factors was presented in Reference [22]. One of the factors could be wind forcing. However, in 2010-2013 (the period considered in Reference [22]), the wind was weak and the wind forcing was hardly noticeable. In 2018, the wind forcing increased dramatically-about 1.5 times compared to 2013 (Table 4), and this suggests that, in 2018, a strong wind forcing was the main factor behind the increasing values of salinity and a g in the surface layer. We can assume that such a situation was observed not only in 2018, but also in 2015. In the latter, the wind impact was strongest for the entire period 2010-2018 (about 1.7 times stronger than in 2018-see Table 4), and the salinity in the 0-5-m layer was slightly higher than in 2018 ( Figure 10). This was not even prevented by quite strong precipitation (total for the period about twice as much as in 2018-see Table 3). Unfortunately, in 2015, there were no direct measurements of the CDOM absorption, but the estimate from satellite data shows that, in June, the values of a g were not lower than in 2018-see Figure 2).