Intense Cyclones in the Black Sea Region: Change, Variability, Predictability and Manifestations in the Storm Activity

: Cyclonic activity in the midlatitudes is a form of general atmospheric circulation, and the most intense cyclones are the cause of hydrometeorological anomalies that lead to economic damage, casualties and human losses. This paper examines the features of variability of intense cyclonic activity in the Black Sea region and the examples of their regional manifestations in the storm types. Based on 6-hourly NCEP / NCAR reanalysis data on 1000 hPa geopotential height ﬁelds with 2 ◦ × 2 ◦ spatial resolution and using the methodology by M.Yu. Bardin, objective data were obtained for the identiﬁcation and estimation of the frequency of deep cyclones (reaching 0.75 and 0.95 quantiles by intensity and depth—intense and extreme cyclones, respectively) for the Black Sea region during the period 1951–2017. Additionally, a speciﬁc methodology of more precise cyclone identiﬁcation based on spherical spline interpolation was successfully applied, and then the two methodologies were compared. The key point of the study is the following: In the background of negative signiﬁcant linear trends and interdecadal variability (period of about 35 years), typical scales of their interannual variability on the periods of about 2.5–3.5 and 6–8 years were identiﬁed. These periods coincide with the time scales of the North Atlantic Oscillation and El Nino–Southern Oscillation, providing an outlook for further study of the patterns of their connection. Besides, seasonal forecasts of frequency of intense cyclones in the Black Sea region were successfully modeled using an artiﬁcial neural network technique. Finally, the case studies of regional manifestations of deep cyclones in the types of storms in the northern Black Sea coast revealed substantial di ﬀ erences in the location of deep centers of cyclones and storm tracks associated with the large-scale pressure ﬁelds.


Introduction
Knowledge of the periods and patterns of regional climate variability is important for the development of long-term forecasts of weather and climate anomalies and for solving applied problems of long-term planning of regional sustainable development under climate change [1][2][3]. Extratropical cyclones in the middle latitudes are responsible for the transfer of heat and moisture, and especially deep ones can have great destructive force due to strong winds [4,5].
According to reanalysis, an increase of global surface temperature is generally accompanied by a decrease in the total number of extratropical cyclones and their packing density in extratropical latitudes of the northern hemisphere, for example in the North Atlantic and European region [6].

Data and Methods
One of the complexities in the field of cyclonic activity research is a wide variety of applied identification and tracking methods (minimum sea level pressure/geopotential height or maximum vorticity) [43] and data sets, mainly reanalysis (NCEP/NCAR, ERA, JRA, MERRA, 20CR, see acronyms in Table A1) or measurements of atmospheric pressure at stations forming triangles and analyzed in terms of percentiles [21,22,44,45].
Comparison between different cyclone identification and tracking methods and input data was made for example in [32,43,46], in particular for the Mediterranean region [47]. Different methods show different cyclone numbers using the same input data, but there is a good agreement in location of cyclogenesis regions, annual cycles and trends of cyclones. If a standardized value, for example anomalies, is used, statistical regime characteristics are lost, but instead a consensus between different methods on variability, sign and significance of trends is reached. Accordingly, different methods do not constrain the analysis of variability of anomalies of cyclones. As for different input data used, the differences are not only in absolute values of cyclonic activity due to the resolution of the model, but even in the trend estimates both between the measurement and reanalysis data, and between different reanalyses, mainly because of small cyclones [32].
A cyclone was identified as a minimum in the 1000 hPa geopotential height field surrounded by closed isohypses with a step of 10 m. The cyclone area was determined by the figure bounded by the Sustainability 2020, 12, 4468 4 of 24 last closed isohypse. Cyclone depth was determined by the difference in geopotential heights between the minimum found at the grid node and the value at the external bounding isohypse. To account for the contribution of the cyclone area to its "power" characteristics, the cyclone intensity (or average depth) was obtained, which is the ratio of the integral cyclone depth to its area.

Deep, Intense and Extreme Cyclones
Deep cyclones were identified by exceeding the depth or intensity thresholds of 0.75 (upper quartile) and 0.95 quantiles by analogy with the method used for calculating extreme precipitation [50] and storminess by pressure data [21,44]. Correspondingly, the series of deep cyclones were obtained: Rd75 (≥40 gpm) and Rd95 (≥90 gpm) by the depth, and Ri75 (≥16.8 gpm) and Ri95 (≥33.7 gpm) by the intensity of cyclones, which were called intense (Rd75 and Ri75) and extreme (Rd95 and Ri95) cyclones. It is necessary to note that the variability of intense cyclones defined as the 75th percentile includes the variability of extreme ones defined as the 95th percentile.
The frequency (dimensionless quantity) of intense and extreme cyclones was determined as the number of centers of cyclonic systems over the Black Sea region exceeding a given threshold (0.75 and 0.95 quantiles) by depth or intensity for a certain period to the total number of times analyzed (equal to the number of days in a period multiplied by 4 times because 6-hourly daily data were used).
Variation of deep cyclones was analyzed in the following terms: -Comparison of seasonal and annual mean values; -Annual cycle of long-term monthly mean and dispersion of the frequency of deep cyclones; -Long-term variability of seasonal and annual frequency of deep cyclones and their anomalies normalized on standard deviation σ; -Linear trends and their p-value with contribution to dispersion (determination coefficient R 2 ); -Variability of low-frequency (≥14 years) and high-frequency (<14 years) components of seasonal and annual series with the assessments of contribution to dispersion; -Spectral Fourier analysis of high-frequency component (<14 years) of seasonal and annual frequency of deep cyclones.

Additional Methodology and Intercomparison of Data Sets
At the same time, a specific methodology [51] of cyclone identification on the basis of spherical spline interpolation was involved in the research for the case studies in Section 3.5. It was done after we had obtained the following result: about 30% of extreme storms in the region were caused by shallow local cyclones occurring over the Black Sea, some of which could not be identified using the method [49] because of a coarser isohypse step. The cyclone tracking algorithm consisted of 3 stages: finding a minima for baric formations in a coarse grid; transferring the pole of the spherical coordinate system into the low pressure center and spherical spline interpolation [52] in order to obtain geometry characteristics of cyclone; combining the low-pressure centers into a single cyclone event at a given distance. To find centers of low pressure, the geopotential data at levels of 1000 hPa (Z 1000 ) and 500 hPa (Z 500 ) were recalculated into sea level pressure (SLP) according to Equation (1): The procedure of low-pressure center search in a coarse grid was as follows: if during sequential considerations of each grid node SLP increased in 4 directions, the starting point was identified as a minimum of SLP field. Next, the pole of the spherical coordinate system was transferred into the center of low pressure, after that subsequent interpolation was performed using spherical splines. Then a closed isobar was found sequentially in 36 directions from the center of a cyclone. Thus, the average radius, pressure gradient and the cyclone area were calculated. The depth of the cyclone was estimated as the SLP difference between its center and the last closed isobar. Then, the identified low-pressure Sustainability 2020, 12, 4468 5 of 24 centers were combined into a single cyclone if the distance between their centers did not exceed 500 km at the same time (every 6 h).
Good quality of both methods used was determined by their comparison in Figure 1. The curve of density probability of deviations between baric systems in the data sets is shown in Figure 1a for thousands of kilometers. The probability density tends to zero for deviations of more than 500 km. Figure 1b shows the number of deviations between the two data sets up to 200 km in more detail. According to it, the maximum repeatability was observed in the deviation range of 20-25 km.
Sustainability 2020, 12, x FOR PEER REVIEW 5 of 25 (a) (b) Figure 1. Results of statistical comparison of baric formations using data set #1 after [49] and the specific data set #2 after [51]: (a) the curve of density probability of deviations between baric systems in compared data sets; (b) histogram of repeatability of deviations between baric formations in compared data sets on 0-200 km spatial scale.

Spectral Analysis
To determine the periods of variability of deep cyclones in the Black Sea region, a spectral Fourier analysis of seasonal and annual values of the frequency of cyclones was carried out [53]. To obtain the spectra estimates, the Parzen window was used for averaging periodograms [54]. The window width was taken equal to 20 with the number of frequency points being 100. A linear trend for the period 1951-2017 was previously removed from the series of frequencies of deep cyclones. Then, a high-frequency filtering with periods of up to 14 years was applied to the residuals of the detrending procedure because it was not possible to determine reliable spectral peaks during longer periods for a 67-year time interval.

Neural Network Model
An approach of modeling and forecasting time variability of intense cyclones in the Black Sea region was tested using the method of artificial neural network. A set of monthly average climate indices was used as input parameters in the neural network model consisting of the indices of largescale signals of the ocean-atmosphere system in the northern hemisphere [18,55], southern hemisphere [56] and equatorial zone (Southern Oscillation index) [57]. These indices were calculated as a difference in pressure or geopotential height between the centers of action of these signals using the monthly average NCEP/NCAR reanalysis fields [48] on the 500 hPa geopotential height and SLP for the period 1948-2019 with a spatial resolution of 2.5° × 2.5°. The indices were calculated and calculations were carried out for each month separately, thereby eliminating seasonal variability. Linear trends of the indices were excluded for the whole period.
To simulate the frequency of intense cyclones in the Black Sea region, a unidirectional heteroassociative single-layer (one hidden layer) neural network with a teacher was used [58]. The neurons of the input, hidden and output layer were represented by a sigmoidal bipolar function tanh(βx). The model assimilated data until 2006 within the 1951-1975 training period and 1976-2006 testing period. After that, there was a period of modeling/predicting cyclones and the validation was made for the 2007-2017 control period.
The learning was based on the backpropagation algorithm according to [59]. The aim of the Results of statistical comparison of baric formations using data set #1 after [49] and the specific data set #2 after [51]: (a) the curve of density probability of deviations between baric systems in compared data sets; (b) histogram of repeatability of deviations between baric formations in compared data sets on 0-200 km spatial scale.

Spectral Analysis
To determine the periods of variability of deep cyclones in the Black Sea region, a spectral Fourier analysis of seasonal and annual values of the frequency of cyclones was carried out [53]. To obtain the spectra estimates, the Parzen window was used for averaging periodograms [54]. The window width was taken equal to 20 with the number of frequency points being 100. A linear trend for the period 1951-2017 was previously removed from the series of frequencies of deep cyclones. Then, a high-frequency filtering with periods of up to 14 years was applied to the residuals of the detrending procedure because it was not possible to determine reliable spectral peaks during longer periods for a 67-year time interval.

Neural Network Model
An approach of modeling and forecasting time variability of intense cyclones in the Black Sea region was tested using the method of artificial neural network. A set of monthly average climate indices was used as input parameters in the neural network model consisting of the indices of large-scale signals of the ocean-atmosphere system in the northern hemisphere [18,55], southern hemisphere [56] and equatorial zone (Southern Oscillation index) [57]. These indices were calculated as a difference in pressure or geopotential height between the centers of action of these signals using the monthly average NCEP/NCAR reanalysis fields [48] on the 500 hPa geopotential height and SLP for the period 1948-2019 with a spatial resolution of 2.5 • × 2.5 • . The indices were calculated and calculations were carried out for each month separately, thereby eliminating seasonal variability. Linear trends of the indices were excluded for the whole period.
To simulate the frequency of intense cyclones in the Black Sea region, a unidirectional heteroassociative single-layer (one hidden layer) neural network with a teacher was used [58]. The neurons of the input, hidden and output layer were represented by a sigmoidal bipolar function tanh(βx). The model assimilated data until 2006 within the 1951-1975 training period and 1976-2006 testing period. After that, there was a period of modeling/predicting cyclones and the validation was made for the 2007-2017 control period.
The learning was based on the backpropagation algorithm according to [59]. The aim of the learning was to determine the neuron weights of each layer in such a way that for a given input vector it could be possible to obtain output values coinciding with the expected values of the training sample with the required accuracy [58]. In the classical statement of the task, the maximum learning is determined by the objective function calculated at each step of the training sample and is defined usually as a quadratic sum of the differences between the actual and expected values of the output signal [58,59]. When it is impossible to achieve coincidence between model calculations and known values with the required accuracy, the maximum learning is defined as the minimum of the objective function [58]. In this study, the training sample consisted of only 30 values, which complicated the determination of the maximum learning by classical method. In this regard, to determine the maximum learning, the method from the paper [60] was used, which was as follows: Correction of the weights of neurons was carried out using the classical objective function, but the training sample was repeated many times, while the values of the weights were saved. At each iteration of the training, the correlation coefficients with the training (r study ) and test (r test ) samples were calculated. The learning maximum had to satisfy 2 conditions: determining the iteration of the learning with the maximum r study and the inequality |r study − r test | < 0.15 had to be satisfied. If the second condition was not fulfilled, the closest iteration corresponding to both conditions was taken as the maximum learning. Since a test sample was used in the search for learning maximum, it was necessary to use an independent control sample for an objective assessment of the model capabilities.
The adaptation of the model was carried out by analogy with the paper [61] and included the stages of preprocessing, modeling and postprocessing. At the stage of preliminary data processing, a search was made for the correlation relationships between the frequency of intense cyclones and the values of the indices of large-scale signals in the previous 2, 4 and 6 months, after that the input parameters having the greatest relationship with the prognostic parameter were selected for the model operational work. At the modeling stage, the model was launched using, at the input, all possible combinations of indices selected at the preliminary stage with different numbers of hidden layer neurons. After training for each such neural network, correlations of model calculations with the results of the test and training samples were remembered. At the final stage, the best 20 neural network combinations were selected based on a ranking of correlation coefficients at a test sample. Then, their ensemble was calculated, the quality of which was estimated in comparison with the control sample.

Material for Case Studies Involving the Types of Storms
Regional manifestations of deep cyclones were studied for the selected large-scale stormy atmospheric situations in the northern Black Sea coast identified earlier in [15]. On the basis of the objective method of cluster analysis applied to observational hydrometeorological data, four types of stormy atmospheric situations (including 2 subtypes) associated with different wave heights in the northern Black Sea coast were identified. They represented four types of meridional circulation: type 1 (western), type 2 (mixed), type 3 (central) and type 4 (eastern). In addition, two subtypes were identified for the types 1 and 2 (western and mixed), for which the macrosynoptic conditions were significantly different. The northwest transfer of air masses was carried out for the subtype 1a of the western type and the northeastern/eastern for the subtype 1b. Similarly, the mixed type 2, dependent on the location of the pressure ridges relative to the Azov-Black Sea basin, was represented by two subtypes: the transfer of air masses carried out from the south to the Black Sea for the subtype 2a and from the north/northwest for the subtype 2b.
Although these types of storminess are based on the features of regional climatic manifestations, they do not contradict the types of large scale atmospheric circulation patterns over midlatitudes of the Sustainability 2020, 12, 4468 7 of 24 northern hemisphere and Central Europe originating from Katz's classification [62], Wangenheim-Girs's classification [63,64] with modifications [65], Grossweterlagen's classification [66][67][68] or types of synoptic situations with deep cyclones in central Europe [69].
Only those types with wave height higher than or equal to 5 m were selected and used as a range of data for daily composite analysis, they were 1b with 13 events, 2a with 17 events, 2b with 10 events and type 3 with 5 events (Table 1). Table 1. Dates of storms (day.month.year) with a wave height higher than or equal to 5 m associated with one of the types of atmospheric processes according to [15] used for composite analysis. Due to the limited number of storm events for each type of large-scale circulation regime, a case study approach was used to analyze regional manifestations of deep cyclones. For the range of dates for each selected type of stormy atmospheric situations with a wave height higher than or equal to 5 m, composite maps of the following parameters were obtained:

1.
Centers of cyclones in the days of storms based on data set #1 after [49]; 2.
Full storm tracks of cyclones based on data set #2 after [51], the centers of which were over the Black Sea region on the composite dates; 3-6. Sea level pressure; 1000 hPa and 500 hPa geopotential height; 1000 hPa wind vector.
Two data sets based on the methodologies by M.Y. Bardin [49] (data set #1) and the current authors' methodology of objective tracking [51] (data set #2) were compared and used to construct composite maps of cyclones to ensure the reliability of results. Other composite maps were obtained using an interactive tool for plotting daily mean composites on the basis of NCEP/NCAR reanalysis data, available at https://www.esrl.noaa.gov/psd/data/composites/day/.

Annual Cycle of Long-Term Mean Characteristics of the Frequency of Deep Cyclones
First, let us consider and compare the mean values and the spread of values (standard deviation σ) of the frequency of the identified series of deep cyclones by seasons and months. As for the total number of cyclones in the region, an increase was expected in winter and a decrease in summer [42].
Comparison of the seasonal (conventional seasons) and annual frequency of the series of deep cyclones can be seen in Figure 2. The frequency of deep cyclones determined by the depth was higher annually and in all seasons than that identified by the intensity, because the depth of cyclones does not take into account their area (and consequently a gradient) in contrast to intensity. As expected, the highest values of the frequency of deep cyclones were in winter (with a maximum in February) while the lowest ones were in summer (with the minimum in August). Mean winter values of the frequency of intense cyclones were 6-8 times higher than summer values, 2.3-2.5 times higher than autumn values and 1.2-1.3 times higher than spring values. The frequency of extreme cyclones in winter was 2.4-2.8 times higher than that in spring and more than 3.5 times higher than that in autumn. Long-term mean of the summer frequency of extreme cyclones in the 1951-2017 period was close to zero (from May to August) because such cyclones are observed only in certain years.
feature of the annual cycle of extreme cyclones (Figure 3b,d) was a steeper descending slope from winter to spring months, resulting in the highest activity of these cyclones from November to March (while for intense cyclones this period was from November to May). As for intense cyclones ( Figure  3a,c), September violated the monotonic increase from summer to winter, characterized by a peak of long-term mean values and standard deviations. The long-term mean values of intense cyclones were higher than the standard deviation from November to April, which indicates that the values tended to be closer to the mean these months, and lower from May to October, when the values were spread out over a wider range.
A feature of intense cyclones was that the largest spread of the values (σ) relative to the mean was observed in January (Figure 3a,c), while the mean values (M) of the time series of both intense and extreme cyclones and standard deviation of extreme cyclones were in February ( Figure 3). That is, in January, the mean values can be less predictable for intense cyclones than, for example, in February, as well as in September for the same reasons described above. This feature is also noted in the Section 3.4. As for extreme cyclones (Figure 3b,d), the standard deviation was higher than the mean values during all months. In this regard and due to low values of frequencies, we did not try to simulate extreme cyclones using the neural network model. The maximum frequency of the deep cyclones of the Northern hemisphere was also shown in February in [8]. In addition, in the paper [8], the maximum frequency of deep cyclones also shifted to the spring months in recent decades, i.e., the deepest cyclones were observed in January and February before the most recent period of global temperature rise (after 1980), and after that the deepest cyclones were observed in February and March. In this research, conducted throughout the entire period 1951-2017, it is noted that mean long-term frequency values of intense cyclones (Ri75 and Rd75) were similar in January and March, with the dominance in March; while for the frequency of extreme cyclones (Ri95 and Rd95), January substantially prevailed over March.
Due to the strong correlation (0.9-0.99) between the frequency of deep cyclones identified by the depth and intensity (Rd75/95 and Ri75/95), and taking into account that intensity is a more complex parameter concerning both depth and area of cyclones, we further understand intense and extreme cyclones as those deep cyclones which were identified by the intensity (respectively, Ri75 and Ri95). The annual cycle on a monthly basis was similar to both intense (Ri75 and Rd75) and extreme (Ri95 and Rd95) cyclones by their long-term mean values and standard deviations ( Figure 3). The feature of the annual cycle of extreme cyclones (Figure 3b,d) was a steeper descending slope from winter to spring months, resulting in the highest activity of these cyclones from November to March (while for intense cyclones this period was from November to May). As for intense cyclones (Figure 3a,c), September violated the monotonic increase from summer to winter, characterized by a peak of long-term mean values and standard deviations. The long-term mean values of intense cyclones were higher than the standard deviation from November to April, which indicates that the values tended to be closer to the mean these months, and lower from May to October, when the values were spread out over a wider range.
A feature of intense cyclones was that the largest spread of the values (σ) relative to the mean was observed in January (Figure 3a,c), while the mean values (M) of the time series of both intense and extreme cyclones and standard deviation of extreme cyclones were in February ( Figure 3). That is, in January, the mean values can be less predictable for intense cyclones than, for example, in February, as well as in September for the same reasons described above. This feature is also noted in the Section 3.4. As for extreme cyclones (Figure 3b,d), the standard deviation was higher than the mean values during all months. In this regard and due to low values of frequencies, we did not try to simulate extreme cyclones using the neural network model.
The maximum frequency of the deep cyclones of the Northern hemisphere was also shown in February in [8]. In addition, in the paper [8], the maximum frequency of deep cyclones also shifted to the spring months in recent decades, i.e., the deepest cyclones were observed in January and February before the most recent period of global temperature rise (after 1980), and after that the deepest cyclones were observed in February and March. In this research, conducted throughout the entire period 1951-2017, it is noted that mean long-term frequency values of intense cyclones (Ri75 and Rd75) were similar in January and March, with the dominance in March; while for the frequency of extreme cyclones (Ri95 and Rd95), January substantially prevailed over March.
Due to the strong correlation (0.9-0.99) between the frequency of deep cyclones identified by the depth and intensity (Rd75/95 and Ri75/95), and taking into account that intensity is a more complex parameter concerning both depth and area of cyclones, we further understand intense and extreme cyclones as those deep cyclones which were identified by the intensity (respectively, Ri75 and Ri95).

Variability, Trends and Anomalies of Seasonal and Annual Frequency of Deep Cyclones
This subsection is devoted to long-term variability of the seasonal (Figure 4a Variability of annual extreme cyclones (Figure 4d) is characterized by a greater contribution of interannual amplitude than that of intense cyclones (Figure 4c). Figure 4c,d shows negative linear trends of the annual frequency of both intense and extreme cyclones, which are significant at the level of 97.5% for intense cyclones and of 92.4% for extreme cyclones (according to Table 2 described below). Variability of 14-year low-frequency filtration residues of annual values (Figure 4c,d) demonstrates an interdecadal period of about 35 years. The significance of this period is not possible to determine using the 67-year data set (1951-2017); however, it is clear that contribution of lowfrequency variability to the variance reached significant values, so it was evaluated further using the coefficient of determination R 2 (see Table 2).

Variability, Trends and Anomalies of Seasonal and Annual Frequency of Deep Cyclones
This subsection is devoted to long-term variability of the seasonal (Figure 4a The most prominent feature of variability of annual frequencies of intense and extreme cyclones (Figure 4c,d) was a peak in the mid-1960s, mainly due to winter and spring contributions (Figure 4a,b), followed by a descent until the 1990s and a period of increase after that up to the end of the time interval. A secondary increase in winter frequency of both intense and extreme cyclones (Figure 4a,b) was found at the beginning of the time interval and after 2000. Autumn frequencies of deep cyclones were low or negligible in the 1980s.
Variability of annual extreme cyclones (Figure 4d) is characterized by a greater contribution of interannual amplitude than that of intense cyclones (Figure 4c). Figure 4c,d shows negative linear trends of the annual frequency of both intense and extreme cyclones, which are significant at the level of 97.5% for intense cyclones and of 92.4% for extreme cyclones (according to Table 2 described below). Variability of 14-year low-frequency filtration residues of annual values (Figure 4c,d) demonstrates an interdecadal period of about 35 years. The significance of this period is not possible to determine using the 67-year data set (1951-2017); however, it is clear that contribution of low-frequency variability to the variance reached significant values, so it was evaluated further using the coefficient of determination R 2 (see Table 2).  This clearly anomalous period approximately from the mid-1960s to the mid-1970s, with a notable maximum in the number of cyclones, was found both in depth and intensity of the total number of cyclones over the Black Sea region on shorter time series in earlier papers [23,42,70], as well as in the frequency of wave heights [15] and water temperature (on the surface and 100 m horizon) [71]. This and other periods of intensification of deep cyclones are consistent with the periods of increasing and decreasing storm activity in the Black Sea region identified in the papers [15,[72][73][74][75][76]. The authors of the paper [72] found the period of positive standardized frequency anomalies of the northern hemisphere cyclones at low latitudes from the mid-1960s to the mid-1970s  Table 2. Coefficients of linear trends of seasonal and annual frequency of deep cyclones (k lin ) and the significance of error (p-value). Estimates (in %) of contribution (determination coefficients R 2 ) to the dispersion of seasonal and annual frequency of intense (Ri75) and extreme (Ri95) cyclones: contribution of the linear trend (R 2 lin ); contribution of the low-frequency component without the linear trend (R 2 ≥14 ); contribution of the low-frequency component taking into account the linear trend (R 2 low ); contribution of the high-frequency component (R 2 high ).   Now, let us specify the intervals and relative values of the periods of low-frequency variability identified in Figure 4 in terms of standardized anomalies (normalized by standard deviation σ) of intense and extreme cyclones referred to as Ria75 and Ria95 ( Figure 5). A period of positive anomalies of annual frequency of the intense cyclones (Figure 5a), with some of them exceeding 3σ, was observed in 1960-1975, while the rest of the time interval was characterized by negative anomalies, with the strongest ones more than 1σ. A slightly shifted period of positive anomalies of annual extreme cyclones (Figure 5b), with a maximum reaching 2.5σ, was observed in 1962-1977, and the rest was characterized by negative anomalies with the module up to about 1σ, except for a couple of years at the beginning of the time interval (1951) that had weak positive anomalies. The character of interdecadal variability of winter anomalies of extreme cyclones (Figure 5d) and spring anomalies of intense cyclones (Figure 5e) was in good agreement with the corresponding annual anomalies (Figure 5a,b). As for the winter anomalies of intense cyclones (Figure 5c This clearly anomalous period approximately from the mid-1960s to the mid-1970s, with a notable maximum in the number of cyclones, was found both in depth and intensity of the total number of cyclones over the Black Sea region on shorter time series in earlier papers [23,42,70], as well as in the frequency of wave heights [15] and water temperature (on the surface and 100 m horizon) [71]. This and other periods of intensification of deep cyclones are consistent with the periods of increasing and decreasing storm activity in the Black Sea region identified in the papers [15,[72][73][74][75][76]. The authors of the paper [72] found the period of positive standardized frequency anomalies of the northern hemisphere cyclones at low latitudes from the mid-1960s to the mid-1970s and the period of negative anomalies in the 1990s. The greatest number of storms caused by cyclones was observed in the west of the Black Sea region in the 1960s and after the year 2000 [74], and off the coast of the Danube Delta in the late 1960s, 1970s and 1990s [75]. Besides that, we found earlier that cyclonic activity (total number of cyclones) was characterized by the following patterns over the Black Sea region: an increase in the frequency of winter cyclones over the central and southern parts of the sea in the 1980s and 1990s, while in the 1960s and 2000s the maximum frequency of cyclones moved to the land to the northern part of the Black sea region, while an increase of frequency of winter anticyclones over the southwestern part of the sea was observed. The results obtained further in Sections 3.4 and 3.5 confirmed that these anomalous episodes shown in the frequency of intense cyclones can be explained by atmospheric dynamics on a regional and global scale.
It can be supposed that the end of the time interval (2017) was characterized by the transition from negative to positive annual anomalies of the frequency of deep cyclones (Figure 5a,b) on the background of the sustainable growth of the interdecadal variability of anomalies. Besides, if recent years have been characterized by an increase in the surface temperature of the Black Sea, then the paper [77] based on numerical experiments on the sensitivity of the regional atmospheric circulation to the sea surface temperature can be referred to. It showed that additional cyclonic vorticity originates in both winter and summer central months for the case of a warm sea [77]. The analyses of long-term variability of the Black Sea surface temperature and intermediate 100 m layer water temperature for the interval 1951-2017 [71] confirmed the presence of interdecadal temperature variations with a period of 40-50 years and an amplitude of up to 0.8 • C, while a common trend representing a sharp rise in temperature over the past 20 years was identified for both horizons. Table 2 shows the coefficients of linear trends of seasonal and annual frequency of deep cyclones (k lin ) with the significance of error (p-value), as well as the contribution of high-frequency (up to 14 years) and low-frequency (≥14 years) components of variability, taking into account and without taking into account the linear trend, to the dispersion of the frequency of deep cyclones. As mentioned above, significant linear trends of the annual frequency of deep cyclones have a negative sign. As for the seasons, significant negative coefficients were identified for the frequency of intense cyclones at the level of more than 98% in winter and spring, and for the frequency of extreme cyclones only in winter at the level of 94.6%. Estimates of the contribution of linear trends (R 2 lin ) to the dispersion of the annual frequencies of intense and extreme cyclones were 7.3% and 4.8%, in particular due to winter and spring: respectively, 9.6% and 9% for intense cyclones and 5.8% and 1% for extreme cyclones.
The trends found for deep cyclones support the tendencies identified in earlier works for the total number of cyclones in the Black Sea region. For example, the paper [23] showed that the frequency of cyclones decreased in all seasons (except summer) as a result of the intensification of the North Atlantic Oscillation in the 1960-1990s and the displacement of the predominant paths of synoptic disturbances to the north.
The contribution of the low-frequency component to the annual values, without taking into account the linear trend (R 2 ≥14 ), was more than 53% for the frequency of intense cyclones (mainly due to winter and spring) and more than 33% for the frequency of extreme cyclones (due to winter).
Taking into account the linear trend, the low-frequency variability (R 2 low = R 2 ≥14 + R 2 lin ) explained more than 60% and 38% of the dispersion of annual frequency of intense and extreme cyclones, respectively, due to the winter season, as well as spring in the case of intense cyclones. The revealed linear and quasiperiodic tendencies could be both an indicator of global warming and part of the lower-frequency variability in the ocean-atmosphere system determined by climate signals of the Atlantic and Pacific oceans [72,78].
As for the contribution of high-frequency variability (R 2 high = 100 − R 2 low ) to the standard deviation of the frequency of intense and extreme cyclones, a summer-autumn maximum of 81-96%, as well as a spring one for extreme cyclones, was noted, while the contribution of high-frequency variability to the annual values reached 39.4% and 61.5%, respectively. For annual values of extreme cyclones, spectral peaks stood out at 5.6, 3.4 and 2.4 years. Herewith, a descending spectrum was observed for intense cyclones, and an ascending one for extreme cyclones, which confirms the previously mentioned feature: for extreme cyclones, more energy is concentrated at higher frequencies (shorter periods), and for intense cyclones at lower frequencies (longer periods).

Spectral Analysis of the High-Frequency Component of Deep Cyclones
The main contribution to the energy of the annual spectra of the high-frequency component of intense cyclones was made by spring, as well as winter, peaks and to a lesser extent by autumn peaks (Figure 6a). As for extreme cyclones, winter spectra made the main contribution to the energy of annual spectra, and the energy of autumn and spring spectra was approximately equal. Close spectral peaks of the high-frequency component of deep cyclones by seasons were identified in the following range of years: • For the frequency of intense cyclones (Figure 6b), 7.4-8.3 in spring and autumn, 5.12-5.7 in winter and summer, from 3.4-3.5 in winter and spring to 3.8 in autumn, and also 2.7 in summer (notably, lower-frequency peaks in the winter-spring period and higher-frequency peaks in the summer-autumn period have higher energy).

•
For the frequency of extreme cyclones (Figure 6c), 5.6, 6.3 and 7.7 in winter, spring and autumn, respectively, and from 3.4 in winter and autumn to 2.5 in spring (summer spectral peaks are insignificant due to the small number of extreme cyclones in this season).
The periods of variability of the frequency of deep cyclones in the Black Sea region identified by spectral analysis were also typical for the global processes in the ocean-atmosphere system, such as the North Atlantic Oscillation-NAO (2.4; 4.8; 7.7 years for the winter index) [79] and El Nino-Southern Oscillation-ENSO (3-4 and 7-8 years) [78]. Thus, the results obtained confirmed that the indicated climatic processes can be responsible for the interannual variability of the frequency of deep cyclones in the Black Sea region. Earlier, we showed in [42] that correlation coefficients of the winter-spring parameters of the total number of cyclones in the Black Sea region exceeded 0.6 with NAO indices, while they reached 0.49 with Southern Oscillation indices at time lags of 4-6 months, and the joint influence of NAO and ENSO significantly exceeded 50% of the interannual variability of cyclones in the studied region. That is why a similar correlation is expected with deep cyclones identified by quantiles. It can be a purpose of our further study.

Approach to Forecasting Intense Cyclones in the Black Sea Region Using the Neural Network Model
Long-term forecast (with a lead time from a month to a year) is an unresolved recent problem. The requirements for the quality of such forecasts still cannot be achieved using modern global climate models. At the same time, the need for high-quality long-term forecasts not only persists, but also becomes more relevant.
A model based on the method of artificial neural network to forecast the frequency of intense cyclones was tested in this paper. Previously, such a modeling approach was successfully applied to predict anomalies in the Nino3.4 region, with a lead time from 1 to 9 months [61]. SLP fluctuations in the centers of action of the global atmosphere are used as predictors in the model.
The ability of the neural network model to project the frequency of intense cyclones for each month separately with different lead times of 2, 4, and 6 months was tested in this paper. For the control period (2007-2017), standard deviations, correlation coefficients and root mean square error (RMSE) were calculated, and Taylor diagrams normalized to standard deviations were obtained (Figure 7). The critical value of the statistically significant correlation for a series of 11 values (from 2007 to 2017) being 0.45 (at p = 0.01) was considered as a threshold for the reliable result-all

Approach to Forecasting Intense Cyclones in the Black Sea Region Using the Neural Network Model
Long-term forecast (with a lead time from a month to a year) is an unresolved recent problem. The requirements for the quality of such forecasts still cannot be achieved using modern global climate models. At the same time, the need for high-quality long-term forecasts not only persists, but also becomes more relevant.
A model based on the method of artificial neural network to forecast the frequency of intense cyclones was tested in this paper. Previously, such a modeling approach was successfully applied to predict anomalies in the Nino3.4 region, with a lead time from 1 to 9 months [61]. SLP fluctuations in the centers of action of the global atmosphere are used as predictors in the model.
The ability of the neural network model to project the frequency of intense cyclones for each month separately with different lead times of 2, 4, and 6 months was tested in this paper. For the control period (2007-2017), standard deviations, correlation coefficients and root mean square error (RMSE) were calculated, and Taylor diagrams normalized to standard deviations were obtained (Figure 7). The critical value of the statistically significant correlation for a series of 11 values (from 2007 to 2017) being 0.45 (at p = 0.01) was considered as a threshold for the reliable result-all calculations with a correlation below this threshold were considered unreliable.  (Figure 7c), for the above mentioned November, February and July, also had high simulation quality for June and September (r = 0.74-0.88, RMSE/σ = 0.64-0.7). A satisfactory quality of modeling was noted for March (r = 0.61), and unsatisfactory for January, April, May, August and December (r < 0.45, RMSE/σ = 0.64). Months with an unsatisfactory forecast can be explained by a different set of predictors (climatic signals) for different lead times, as well as by a large standard deviation of the values, as mentioned in Section 3.1.
A comparison of the quality of the modeling of the annual time series of the frequency of intense cyclones for different lead times was drawn from Figure 8. The correlations between calculated and simulated frequency of intense cyclones for lead times 2 (see Figure 8a), 4, (see Figure 8b) and 6 (see Figure 8c) months were 0.66, 0.66 and 0.69, respectively, and RMSEs were 0.039, 0.039, and 0.037 with an average deviation of 0.061.
According to Figures 7 and 8, the model with a greater forecast lead time (6 months) appears to have been slightly better than models with lower lead times (2 and 4 months). It requires further study because the found feature may be associated with a typical period of manifestation of the global climate signal El Nino-Southern Oscillation, which is about 6 months for the Atlantic-European region [30,80,81]. After this time, the mature phase of ENSO and the responding rearrangement of the global atmospheric circulation develop [82,83].  (Figure 7c), for the above mentioned November, February and July, also had high simulation quality for June and September (r = 0.74-0.88, RMSE/σ = 0.64-0.7). A satisfactory quality of modeling was noted for March (r = 0.61), and unsatisfactory for January, April, May, August and December (r < 0.45, RMSE/σ = 0.64). Months with an unsatisfactory forecast can be explained by a different set of predictors (climatic signals) for different lead times, as well as by a large standard deviation of the values, as mentioned in Section 3.1.
A comparison of the quality of the modeling of the annual time series of the frequency of intense cyclones for different lead times was drawn from Figure 8. The correlations between calculated and simulated frequency of intense cyclones for lead times 2 (see Figure 8a), 4, (see Figure 8b) and 6 (see Figure 8c) months were 0.66, 0.66 and 0.69, respectively, and RMSEs were 0.039, 0.039, and 0.037 with an average deviation of 0.061.
According to Figures 7 and 8, the model with a greater forecast lead time (6 months) appears to have been slightly better than models with lower lead times (2 and 4 months). It requires further study because the found feature may be associated with a typical period of manifestation of the global climate signal El Nino-Southern Oscillation, which is about 6 months for the Atlantic-European region [30,80,81]. After this time, the mature phase of ENSO and the responding rearrangement of the global atmospheric circulation develop [82,83].

Case Studies of Regional Manifestations of Deep Cyclones in the Types of Storms in the Northern Black Sea Coast
We conducted the analysis of composite maps of centers of cyclones in the days of storms from data set #1, full storm tracks from data set #2, sea level pressure, 1000 hPa and 500 hPa geopotential height and 1000 hPa wind vector for the selected four types (1b, 2a, 2b and 3) of storminess according to [15] associated with wave heights of no less than 5 m at the northern Black Sea coast.
The storm type 1b (13 events) was characterized mainly by storm tracks crossing the Black Sea from southwest (SW) to northeast (NE) with deep centers of cyclones on the dates of storms located along this line and focused in the northern sector of the region. Minimum SLP less than 1010 hPa was located to the north of Crimea, as was minimal 1000 hPa geopotential height lower than 85 hPa, albeit slightly shifted to northeast. The west of the northern Black Sea coast was occupied by a trough in the field of 500 hPa geopotential height and southeastward 1000 hPa wind vector.
For the storm type 2a (17 events), storm tracks were located along the same line from SW to NE, but shifted to the north/northwest of the northern Black Sea coast, which was especially noticeable for deep cyclonic centers on the dates of storms. Locations of the minimal SLP and 1000 hPa geopotential height were generally the same as for the previous type 1b, but with a substantially lower central value: ˂1004 hPa for SLP and ˂10 m for geopotential height. A southern edge of this center over the northern Black Sea coast was characterized by the NW-SE direction of a 1000 hPa wind vector and 500 hPa geopotential height gradient.
While full storm tracks for the type 2b (10 events) were identified to the north and east of the Black Sea, majority of them were located to the northeast of the northern Black Sea coast. Most deep centers of cyclones on the dates of storms were located in the northeastern sector, unlike the two previous types (1b and 2a). Areas of minimal values of SLP, 1000 hPa and 500 hPa geopotential height, were located to the NE of the northern Black Sea coast, falling below 1004 hPa for the SLP, 30 m for the 1000 hPa geopotential height and 5200 m for the 500 hPa geopotential height (the lowest values

Case Studies of Regional Manifestations of Deep Cyclones in the Types of Storms in the Northern Black Sea Coast
We conducted the analysis of composite maps of centers of cyclones in the days of storms from data set #1, full storm tracks from data set #2, sea level pressure, 1000 hPa and 500 hPa geopotential height and 1000 hPa wind vector for the selected four types (1b, 2a, 2b and 3) of storminess according to [15] associated with wave heights of no less than 5 m at the northern Black Sea coast.
The storm type 1b (13 events) was characterized mainly by storm tracks crossing the Black Sea from southwest (SW) to northeast (NE) with deep centers of cyclones on the dates of storms located along this line and focused in the northern sector of the region. Minimum SLP less than 1010 hPa was located to the north of Crimea, as was minimal 1000 hPa geopotential height lower than 85 hPa, albeit slightly shifted to northeast. The west of the northern Black Sea coast was occupied by a trough in the field of 500 hPa geopotential height and southeastward 1000 hPa wind vector.
For the storm type 2a (17 events), storm tracks were located along the same line from SW to NE, but shifted to the north/northwest of the northern Black Sea coast, which was especially noticeable for deep cyclonic centers on the dates of storms. Locations of the minimal SLP and 1000 hPa geopotential height were generally the same as for the previous type 1b, but with a substantially lower central value: <1004 hPa for SLP and <10 m for geopotential height. A southern edge of this center over the northern Black Sea coast was characterized by the NW-SE direction of a 1000 hPa wind vector and 500 hPa geopotential height gradient.
While full storm tracks for the type 2b (10 events) were identified to the north and east of the Black Sea, majority of them were located to the northeast of the northern Black Sea coast. Most deep centers of cyclones on the dates of storms were located in the northeastern sector, unlike the two previous types (1b and 2a). Areas of minimal values of SLP, 1000 hPa and 500 hPa geopotential height, were located to the NE of the northern Black Sea coast, falling below 1004 hPa for the SLP, 30 m for the 1000 hPa geopotential height and 5200 m for the 500 hPa geopotential height (the lowest values in comparison with the previous types 1b and 2a). The direction of 1000 hPa wind vector was southeastward along the entire northern Black Sea coast.
The last storm type described here, type 3, included the least number of cases (5 events) but was still characterized by relatively dense frequency of full storm tracks and deep centers of cyclones on the dates of storms concentrated to the southeast of the Black Sea, where the areas of minimum values of the SLP (lower 1011 hPa) and 1000 hPa geopotential height (lower 90 m) were located as well. The minimal 500 hPa geopotential height was located over the NW Black Sea shelf. The 1000 hPa wind vector was mostly westward over the northern Black Sea coast.
Taking into account the geometry of cyclones (their entire area) from the data set #2, the composite maps of the frequency of cyclones as closed systems were obtained (Figure 9). The maps, grouped for these four storm types, show the frequency of passage of the entire area of a cyclone through a specific space point. Maximum passage frequency of the total number of cyclones for type 1b was located over the Crimean peninsula (Figure 9a), for type 2a to the north-northwest of the Black Sea coast (Figure 9b), for type 2b to the northeast (Figure 9c) and for type 3 to the south-southeast (Figure 9d).
These results obtained for the groups of severest storms (with wave heights equal to or higher than 5 m) agree with the early results by Chernyakova et al. [84,85] for separate storms in 1965-69, where baric fields and position of synoptic systems were shown for different stormy wind directions over the Black Sea. Two types of storm fields of atmospheric pressure with the position of the cyclonic center to the north or south of the Black Sea were distinguished recently in paper [86]. Our case studies showed that types of storms with wave height higher than or equal to 5 m over the northern Black Sea coast were associated with groups of centers of deep cyclones located relative to this coast in the north (type 1b), northwest (type 2a), northeast (type 2b) and southeast (type 3).
Sustainability 2020, 12, x FOR PEER REVIEW 17 of 25 in comparison with the previous types 1b and 2a). The direction of 1000 hPa wind vector was southeastward along the entire northern Black Sea coast. The last storm type described here, type 3, included the least number of cases (5 events) but was still characterized by relatively dense frequency of full storm tracks and deep centers of cyclones on the dates of storms concentrated to the southeast of the Black Sea, where the areas of minimum values of the SLP (lower 1011 hPa) and 1000 hPa geopotential height (lower 90 m) were located as well. The minimal 500 hPa geopotential height was located over the NW Black Sea shelf. The 1000 hPa wind vector was mostly westward over the northern Black Sea coast.
Taking into account the geometry of cyclones (their entire area) from the data set #2, the composite maps of the frequency of cyclones as closed systems were obtained (Figure 9). The maps, grouped for these four storm types, show the frequency of passage of the entire area of a cyclone through a specific space point. Maximum passage frequency of the total number of cyclones for type 1b was located over the Crimean peninsula (Figure 9a), for type 2a to the north-northwest of the Black Sea coast (Figure 9b), for type 2b to the northeast (Figure 9c) and for type 3 to the southsoutheast (Figure 9d).
These results obtained for the groups of severest storms (with wave heights equal to or higher than 5 m) agree with the early results by Chernyakova et al. [84,85] for separate storms in 1965-69, where baric fields and position of synoptic systems were shown for different stormy wind directions over the Black Sea. Two types of storm fields of atmospheric pressure with the position of the cyclonic center to the north or south of the Black Sea were distinguished recently in paper [86]. Our case studies showed that types of storms with wave height higher than or equal to 5 m over the northern Black Sea coast were associated with groups of centers of deep cyclones located relative to this coast in the north (type 1b), northwest (type 2a), northeast (type 2b) and southeast (type 3). (a)

Conclusions
This study was carried out as a part of ongoing research of regional climatic-environmental anomalies in the Black Sea region. Earlier papers were devoted to the total number of cyclones in the region [23,42,70] or individual cases of intense cyclones [40], but not to their climatology. The focus of this study was on the long-term variability of the frequency of deep cyclones, which have an important socioeconomic impact. Intense cyclones were identified as the 0.75 quantile and extreme cyclones as the 0.95 quantile of the total number of cyclones (while the variability of intense cyclones includes the variability of extreme cyclones). We studied the characteristics and described the peculiarities of their seasonal and interannual-interdecadal variability in the terms of conventional statistical and spectral analysis, tested the possibility of their predictability using the developed model of neural networks and studied their relationship with the types of storminess on the northern Black Sea coast.
It was shown that the frequency of deep cyclones identified by the depth was higher in all seasons than that identified by intensity, because the depth of cyclones in our study did not take into account their area (and consequently a gradient) in contrast to intensity. The highest values of the frequency of deep cyclones were in winter, with a maximum in February (except for the maximum of standard deviation of the frequency of intense cyclones in January), and the lowest (or close to zero for extreme cyclones) were in summer, with a minimum in August, with the 6-8 times difference between the maximum and minimum for intense cyclones. A feature of the annual cycle of extreme cyclones is a steeper descending frequency from winter to spring months. Variability of extreme cyclones has a higher contribution of interannual amplitude to the total variance than that of intense cyclones.
The most prominent feature of variability of annual values of deep cyclones was a peak in the mid-1960s (up to 3σ standardized anomalies for intense cyclones and 2.5σ for extreme cyclones), mostly pronounced in winter and also in spring for the intense cyclones. It was followed by a decrease until the 1990s (lower 1σ), during which autumn frequencies of deep cyclones were low or negligible, and the period of increase after that up to the end of the time interval (mostly due to winter). Low frequency variability after 14 years filtering was characterized by the periodicity of about 35 years. It is possible that the end of the time interval (2017) was characterized by the transition from negative to positive annual anomalies of the frequency of deep cyclones on the background of sustainable growth of the interdecadal residuals (≥14 years) of variability of anomalies.
Negative linear trends for the annual frequency of deep cyclones were estimated as significant at the level of 97.5% for intense cyclones (coefficient of determination R 2 = 7.3%) due to winter and spring and at the level of 92.4% for extreme cyclones (R 2 = 4.8%) due to winter. The contribution of the low-frequency component (≥14 years) to variance of the annual values (R 2 ) reached more than 53% (60% with the contribution of linear trend) for the frequency of intense cyclones (mainly due to winter and spring) and more than 33% (38%, respectively) for the frequency of extreme cyclones (due to winter). As for the contribution of high-frequency variability, a summer-autumn maximum of 81-96% for intense cyclones, as well as a spring one for extreme cyclones, was observed, while the contribution of the high-frequency component to the annual variance reached 39.4% and 61.5%, respectively.
Spectral analysis of the high-frequency component (up to 14 years) of the seasonal and annual frequency of deep cyclones showed a descending spectrum for intense cyclones and an ascending one for extreme cyclones, which indicates that for extreme cyclones more energy was concentrated at higher frequencies, while for intense cyclones at lower frequencies. Typical periods of interannual variability were identified as 2.5/2.4 and 7.7/5.6, 3.4 years for intense cyclones (due to spring contribution)/extreme cyclones (due to winter contribution).
Modeling using neural networks allows successful forecasting of the frequency of intense cyclones in the Black Sea region, with a lead time of up to 6 months for all months except December-January and July-August (probably due to the different set of predictors and higher standard deviations as compared with the mean values these months). Consequently, there is a potential of the neural network modeling for forecasting cyclones in the Black Sea region. Case studies of the role of deep cyclones in the generation of storms with wave heights above 5 m in the northern Black Sea coast showed that the types of storms were associated with specific baric situations when deep cyclone centers were concentrated in the north (type 1b), northwest (type 2a), northeast (type 2b) and southeast (type 3) of the Black Sea coast.
The key findings, which are of the greatest practical interest for the sustainable development of the studied region, can be outlined as follows: - The highest frequency of intense cyclones is from November to May, while for extreme cyclones this period is from November to March with the maximum in February for the both; -When characterizing the regime and variability of deep cyclones in the region, it is necessary to take into account significant linear trends for some seasons (which are responsible for up to 10% contribution to dispersion in winter), as well as the low-frequency component of variability (which makes up to half of the contribution to dispersion for extreme cyclones in winter and more than half of the contribution to dispersion for intense cyclones in winter and spring); -Typical periods of interannual variability of winter-spring deep cyclones are about 2.5-3.5 and 6-8 years, the same as for the North Atlantic Oscillation and El Nino-Southern Oscillation; - The neural network modeling for forecasting cyclones in the region is a promising approach and requires a separate study; - The plausible pattern of the distribution of deep cyclone centers associated with the types of storms on the northern Black Sea coast is shown.
The outlook for further research of intense cyclones in the Black Sea region is important, not only in terms of developing the theory of long-term forecasts of weather and climate anomalies (based on the patterns of regional manifestations of large-scale processes in the ocean-atmosphere system), but also for solving the applied problems of long-term planning of regional sustainable development under climate change. Examples include a number of developed projects and activities for solving urgent regional problems, taking into account trends and natural patterns of change of dangerous hydrometeorological phenomena. They concern protecting and strengthening coasts, risk management of the coastal zone, prevention of emergencies of technogenic and environmental problems, long-term perspective planning of safe activities of the recreational, agricultural, mining sectors, the city's fresh water sector (taking into account the level of reservoir filling), seaports and engineering structures.