The Black Sea Upwelling System: Analysis on the Western Shallow Waters

: Upwelling is due to the combined effect of the coastal divergence process and Ekman pumping. The author aims to investigate two new upwelling indices for the Black Sea, derived from climate reanalysis models and the following in situ data: (a) the Coastal Upwelling Transport Index (CUTI) that estimates the rate of vertical volume transport and (b) the Biologically Effective Upwelling Transport Index (BEUTI) that estimates the nitrate flux into the surface mixed layer. Average monthly wind by the European Centre for Medium-Range Weather Forecasts (ECMWF) and Copernicus Marine Services for the Black Sea basin is used to calculate the CUTI and BEUTI Indexes for over 26 years (1993–2019) to analyse the sites along the North-Western Black Sea where changes in divergence phenomena occur. From 2000 to 2018, 31 divergence processes were observed based on daily in situ data from the coastal monitoring stations, with significant predominance in late spring and early summer. Nitrate supply by coastal upwelling has been estimated by combining sea surface temperature and salinity for the in situ data for the North-Western Black Sea shallow waters, and BEUTI indices were determined. Comparing 18 years of data results, the calculated indices and the observed upwelling events showed significant correlations


Introduction
Scaled by the baroclinic Rossby radius of deformation, coastal upwelling is essential in the world's ocean and semi-enclosed seas like the Black Sea as a driver of ecosystem productivity.The horizontal scale of geostrophic relaxation is the "natural" scale of baroclinic boundary currents, eddies and fronts.In addition, it sets the scale of the waves that grow most rapidly due to baroclinic instability and long Rossby wave speed [1][2][3].Along the western coast of the world's continents, the most productive regions are known as the eastern boundary currents of the Eastern Boundary Upwelling System (EBUS), representing only about 1% of the global ocean [4].The dynamics of wind-driven coastal upwelling were first described mathematically by V. Wilfrid Ekman [5].For over a century, the Ekman theory stood as foundational dynamical oceanography and the basis for understanding wind-driven ocean circulation [6].Ekman determined that the net transfer of water is at a 90-degree angle to the wind direction, and the total water transport in the Ekman layer displaced to the right of the prevailing winds is called Ekman transport.The Ekman transport results from a balance between the Coriolis effect and turbulent drag forces within surface waters and geographic features (the shoreline and seabed).
To describe the upwelling mechanisms and indicators, refs.[7,8] developed the Bakun Index as a simple time series to represent variations in upwelling along the coast, calculated based on Ekman's theory.Ref. [9] estimates the offshore Ekman transport driven by geostrophic wind stress (Environmental Research Division ERD or upwelling index) through calculations derived from mean surface atmospheric pressure fields provided by the US Navy Fleet Numerical Meteorological and Oceanographic Center (FNMOC).
The upwelled deep water on the sea's surface is generally cold and rich in nutrients like nitrogen and phosphorus.At the surface, nutrient-rich water is supplied to the sunlit surface layer by wind-driven upwelling, which stimulates phytoplankton growth and ultimately fuels diverse and productive marine ecosystems.As a result, plankton populations often multiply in the upwelling layer, which is a phenomenon called phytoplanktonic bloom.In addition to this bottom-up approach, forcing through the base of the food web, upwelling can influence higher trophic levels directly through exposure to the physical and chemical signatures of the deep ocean (e.g., lower temperature, oxygen concentration, and pH).Nitrate storage by phytoplankton cells during the early stages of upwelling development was studied in [10].In non-and upwelling events, they found that intracellular nitrate reserves account only for less than 2% of daily primary production during upwelling events.Impacts of upwelling variability can be felt on timescales ranging from single events (days) to decades and over a more extended time period.In addition, they propagate commercial and recreational activities that derive considerable socio-economic benefits [11].
Upwelling tends to be highly seasonal in temperate regions, where it peaks in the spring-summer but tends toward year-round continuity in most tropical areas [12].Located in a temperate region, the Black Sea, the world's most isolated sea, is connected to the Mediterranean Sea through the Marmara Sea and the Aegean Sea by the narrow and shallow Bosporus and Dardanelles Straits, respectively.As well as in other parts of the ocean, the Ekman pumping causes intense upwelling on the western coast [13], the western cyclonic gyre [14], and the southern part of the Black Sea basin.The upwelling season occurs on the North-Western shelf from April to September, strengthening the across-shelf exchange of rich-nutrient deep water within the surface layer.
Studies of upwelling events for the Black Sea have been based mainly on in situ and satellite data.Based on satellite data evaluation, refs.[15,16] summarise that the South Crimea coast, the Western part of the Turkish coast, and the Bulgarian and Romanian beaches have frequent upwelling events.For near-coastal North-Western Black Sea waters, refs.[13,14] concluded that strong southern wind or moderate-to-strong westerly wind generates offshore-oriented currents that favour surface water replenishments by deep waters (characterised by low temperature and high salinity).
This paper aims to explore the divergent processes and increase understanding of the physical-chemical-biological interaction dynamics in Black Sea basin upwelling regions, emphasising the North-Western Shelf.The analysis will consist of integrated data, including in situ and satellite data, and numerical models to understand the upwelling dynamics.

Area of Interest
The Black Sea (BS) is an inland sea in southeastern Europe, bordered by six countries: Bulgaria, Romania, Ukraine, Russia, Georgia, and Turkey.Unlike many coastal regions worldwide that experience significant upwelling, the Black Sea does not have strong, continuous upwelling events due to its unique physical and hydrographic characteristics.Instead, the Black Sea is considered a strongly stratified sea, meaning it has distinct layers of water that do not mix as thoroughly as in a typical ocean.Rivers play a substantial role in pollution and eutrophication processes, leading to profound alterations in the Black Sea ecosystem.
Four major rivers, namely the Danube, Dnieper, Dniester, and Southern Bug, which collectively possess a water basin size 23 times larger than that of the Northwest Shelf (NWS) area [17], serve as conduits for industrial, agricultural, and domestic wastewater from several European countries into the Black Sea.However, these waterborne effluents more frequently propagate southward along the western coast (Figure 1).In either scenario, the coastal water region, extending from the Danube mouth to the Dnieper-Bug estuary, bears the brunt of the highest anthropogenic load.In some instances, this area even receives the combined flow of all four rivers.This is particularly true when the Danube's flow diverts southward, resulting in an influx of eutrophicated waters from the north. north.
These waterways are major contributors to the influx of pollutants and exces nutrients, which trigger significant transformations in the delicate balance of life within the Black Sea.Between the 1970s and 1980s decades, the Black Sea, especially its northwes shelf, underwent a dramatic shift in ecological balance.This change led to persistent and extensive hypoxia, causing significant disruption to the pelagic (water column) and benthic (sediment) ecosystems, leaving enduring impacts in its wake.As a result, th ecosystem of the Black Sea has undergone dramatic shifts with far-reaching consequences Nonetheless, the most environmentally precarious conditions in the Northwest Shel (NWS) manifest during the warmer months.At the onset of the spring-summer period typically from April to May [13,14], rising water temperatures induce a stratification phenomenon within the NWS, encompassing temperature and salinity variations.Thi stratification facilitates the dispersion of altered river runoff across the surface, resulting in a notable elevation of nutrient concentrations within the warm, brackish surface waters Consequently, this triggers a robust proliferation of phytoplankton, which gradually consume the abundant nutrients, simultaneously enhancing the dissolved oxygen level through photosynthesis in the surface waters.
Regarding its physical characteristics, the NWS Black Sea stretches from the shorelin to the shelf's edge, reaching depths of 130 to 200 m, with its widest point spanning 22 kilometres.This expansive region encompasses the majority of the western Black Sea defined in this context as the area west of the narrowest basin width, approximately 25 kilometres wide, between the Crimea peninsula and the Turkish coast.The seabed terrain exhibits a smooth, uniform profile shaped by the continuous sediment deposition carried by the major rivers that discharge into this sector of the Black Sea.However, the shelf' overall topography gains complexity due to specific geographical features, including th presence of the Danube Delta, the coastline's northward curvature beyond the delta irregularities along the coastline caused by prominent capes, and a reduction in the shelf' width immediately south of the Danube River.These waterways are major contributors to the influx of pollutants and excess nutrients, which trigger significant transformations in the delicate balance of life within the Black Sea.Between the 1970s and 1980s decades, the Black Sea, especially its northwest shelf, underwent a dramatic shift in ecological balance.This change led to persistent and extensive hypoxia, causing significant disruption to the pelagic (water column) and benthic (sediment) ecosystems, leaving enduring impacts in its wake.As a result, the ecosystem of the Black Sea has undergone dramatic shifts with far-reaching consequences.
Nonetheless, the most environmentally precarious conditions in the Northwest Shelf (NWS) manifest during the warmer months.At the onset of the spring-summer period, typically from April to May [13,14], rising water temperatures induce a stratification phenomenon within the NWS, encompassing temperature and salinity variations.This stratification facilitates the dispersion of altered river runoff across the surface, resulting in a notable elevation of nutrient concentrations within the warm, brackish surface waters.Consequently, this triggers a robust proliferation of phytoplankton, which gradually consume the abundant nutrients, simultaneously enhancing the dissolved oxygen levels through photosynthesis in the surface waters.
Regarding its physical characteristics, the NWS Black Sea stretches from the shoreline to the shelf's edge, reaching depths of 130 to 200 m, with its widest point spanning 220 km.This expansive region encompasses the majority of the western Black Sea, defined in this context as the area west of the narrowest basin width, approximately 250 km wide, between the Crimea peninsula and the Turkish coast.The seabed terrain exhibits a smooth, uniform profile shaped by the continuous sediment deposition carried by the major rivers that discharge into this sector of the Black Sea.However, the shelf's overall topography gains complexity due to specific geographical features, including the presence of the Danube Delta, the coastline's northward curvature beyond the delta, irregularities along the coastline caused by prominent capes, and a reduction in the shelf's width immediately south of the Danube River.

Methods to Calculate the Upwelling Indexes
Here, we briefly review the upwelling indices and their key differences.Details of Coastal Upwelling Transport Index (CUTI) and Biologically Effective Upwelling Transport Index (BEUTI) calculations and an exhaustive comparison with the Bakun Index are provided in [6].
CUTI is derived from the surface wind stress, sea surface height, and mixed layer depth.CUTI is an estimation of vertical transport into (or out of) the surface mixed layer (SML) as the sum of two components: vertical transport in near-surface Ekman transport (U Ek ) and vertical transport associated with an alongshore pressure gradient (U Geo ) CUTI = U Ek + U Geo . ( The Bakun Index can be obtained using the sea level pressure fields from an atmospheric reanalysis to derive estimated near-surface winds.Details of the peculiarity of the Bakun Index calculation are available in technical reports in [7][8][9]. For the Biologically Effective Upwelling Transport Index (BEUTI), the volume transport of upwelling/downwelling and the amount of nitrate upwelled/downwelled is derived by multiplying the vertical transport at the mixed layer base: In contrast to the Bakun Index, CUTI and BEUTI use winds directly from atmospheric reanalyses that assimilate satellite and in situ wind measurements.As a result, they account for cross-shore geostrophic transport, which can suppress or enhance upwelling relative to the wind-driven Ekman transport and fully consider the upwelling driven by wind stress curls.In contrast, the Bakun Index omits upwelling associated with alongshore changes in the wind.
Another method proposed in the specialised literature is represented by the simple thermal difference (∆T) between the coastal and the offshore SST acquired at the same latitude to estimate the latitudinal and temporal variability of the intensity of the upwelling, using the following formulation: The maximum temperature, SST max , is chosen as the offshore temperature, which used to be the coastal SST, where the influence of the upwelling is predictably insignificant.The SST max determines the SST-based indices from various locations at arbitrary distances from the coast.

Data and Computation
In the present study, the Coastal Upwelling Transport Index (CUTI) of the National Oceanic and Administrative Administration's Environmental Research Division (NOAA-ERD) was used to derive the time series of the coastal upwelling index in four locations on the north-western Black Sea coast.The coastal upwelling index time series was calculated using monthly average wind fields from the European Centre for Medium-Range Weather Forecasts (ECMWFs) reanalysis [19] and MATLAB to compute the CUTI (Figure 2).
values of UI as a response to upwelling or downwelling favourable winds.
To compute the BEUTI, Copernicus Marine Service [20] data are used for dedicated locations.For in situ records of the sea surface temperature (SST), sea surface salinity (SSS), and nitrates can be used to extract the upwelling events, and we also used longterm daily data from the monitoring station of the NIMRD at the Mamaia Casino site (44°14′ N and 28°38′ E) [13,14].

Results
In this paper, the sea surface temperature and salinity data are used as indicators of the upwelling phenomena.Data were collected with YSI CastAway CTD at Mamaia Casino (44°14′ N and 28°38′ E).Referring to the occurrence statistics, from 18 years of analysis (2000 to 2018), 31 divergence processes were observed based on daily in situ data from the coastal monitoring station, with significant predominance in late spring and early summer (Figures 3 and 4).
Figure 3 illustrates the occurrences of recorded upwelling events in the nearshore waters of the Western Black Sea, specifically at Constanta, during the period spanning from May to September from 2000 to 2019.Notably, divergence processes occur in May and June nearly every year, with one event per period being the norm and sometimes extending into July.Ekman transport calculations involve the determination of zonal and meridional Ekman transport components by computing the wind stress values τ x and τ y (in N m −2 ).These wind stress values are derived from the monthly mean values of the zonal (u) and meridional (v) components of "pseudo wind stress" (m 2 s −2 ): where ρ a = 1.2 kg/m 3 (air density), C d = 1.4 × 10 3 (drag coefficient), x corresponds to the zonal component, and y corresponds to the meridional wind stress component.Finally, the upwelling index (UI) is computed using CUTI Formula (1), defined as CUTI (m 3 •s −1 •100 m −1 ), representing the volume transport per distance unit of an alongshore section.The sign of Ekman transport is changed to define positive or negative values of UI as a response to upwelling or downwelling favourable winds.
To compute the BEUTI, Copernicus Marine Service [20] data are used for dedicated locations.For in situ records of the sea surface temperature (SST), sea surface salinity (SSS), and nitrates can be used to extract the upwelling events, and we also used long-term daily data from the monitoring station of the NIMRD at the Mamaia Casino site (44 • 14 ′ N and 28 • 38 ′ E) [13,14].

Results
In this paper, the sea surface temperature and salinity data are used as indicators of the upwelling phenomena.Data were collected with YSI CastAway CTD at Mamaia Casino (44 • 14 ′ N and 28 • 38 ′ E).Referring to the occurrence statistics, from 18 years of analysis (2000 to 2018), 31 divergence processes were observed based on daily in situ data from the coastal monitoring station, with significant predominance in late spring and early summer (Figures 3 and 4).
Figure 3 illustrates the occurrences of recorded upwelling events in the nearshore waters of the Western Black Sea, specifically at Constanta, during the period spanning from May to September from 2000 to 2019.Notably, divergence processes occur in May and June nearly every year, with one event per period being the norm and sometimes extending into July.The intensity of upwelling and downwelling varies significantly among the stations.Some stations experience stronger upwelling or downwelling than others.For example, the ECTA5 station tends to have stronger upwelling than other stations, while the DANUBE station experiences stronger downwelling than the rest.In addition, some stations, such as ECTA1 and ECTA3, exhibit considerable interannual variability in upwelling intensity.There seems to be a general trend of decreasing upwelling intensity.The intensity of upwelling and downwelling varies significantly among the stations.Some stations experience stronger upwelling or downwelling than others.For example, the ECTA5 station tends to have stronger upwelling than other stations, while the DANUBE station experiences stronger downwelling than the rest.In addition, some stations, such as ECTA1 and ECTA3, exhibit considerable interannual variability in upwelling intensity.There seems to be a general trend of decreasing upwelling intensity.Some stations experience stronger upwelling or downwelling than others.For example, the ECTA5 station tends to have stronger upwelling than other stations, while the DANUBE station experiences stronger downwelling than the rest.In addition, some stations, such as ECTA1 and ECTA3, exhibit considerable interannual variability in upwelling intensity.There seems to be a general trend of decreasing upwelling intensity.The intensity of upwelling and downwelling varies significantly among the stations.Some stations experience stronger upwelling or downwelling than others.For example, the ECTA5 station tends to have stronger upwelling than other stations, while the DANUBE station experiences stronger downwelling than the rest.In addition, some stations, such as ECTA1 and ECTA3, exhibit considerable interannual variability in upwelling intensity.There seems to be a general trend of decreasing upwelling intensity.For ECTA1, nearshore water, Figure 7 shows a clear seasonal pattern in the CUTI, with positive values occurring predominantly in the spring and summer and negative values occurring more commonly in autumn and winter.In contrast, Figure 8 reveals how, for the offshore waters, the CUTI values mainly fluctuate around zero, indicating a relatively balanced state between upwelling and downwelling events throughout the observed period.Moreover, considerable variability in the intensity of upwelling can be observed yearly, with some years experiencing much stronger upwelling events than others (Figure 7).However, Figure 8 appears to be a weak seasonal pattern, with slightly more positive CUTI values during the spring and summer than autumn and winter.While there are fluctuations in the CUTI over the 26 years, there does not appear to be a consistent trend towards increasing or decreasing upwelling intensity (Figure 7) or no exceptionally strong upwelling or downwelling events observed in the time series (Figure 8).Year-to-year fluctuations in upwelling intensity and frequency, as revealed by the CUTI, underscore the dynamic nature of this process in the nearshore zone.These fluctuations, coupled with the pronounced seasonal variations in upwelling, highlight the significant influence of seasonal factors on driving upwelling events in this region.The offshore area of Constanta experiences a relatively moderate and balanced upwelling regime, with weak seasonal variations and notable interannual fluctuations.This information is valuable for understanding the oceanographic processes and potential impacts on marine ecosystems in this region.
Figure 9 shows that the SST exhibits a clear seasonal pattern, with higher temperatures in the summer months and lower temperatures in the winter months.This pattern is consistent for both nearshore and offshore locations.The nearshore SST tends to be slightly higher than the offshore SST, particularly during the warmer months.This could be due to the influence of shallower waters and land-based heat sources.The CUTI shows negative values during periods of upwelling.These events appear to occur more frequently during spring and summer, corresponding with periods of higher SST.There does not appear to be a significant long-term trend in either SST or the frequency and intensity of upwelling events over the time period covered by the graph.Figure 9 illustrates the complex interplay between SST, upwelling, and seasonality in the study For ECTA1, nearshore water, Figure 7 shows a clear seasonal pattern in the CUTI, with positive values occurring predominantly in the spring and summer and negative values occurring more commonly in autumn and winter.In contrast, Figure 8 reveals how, for the offshore waters, the CUTI values mainly fluctuate around zero, indicating a relatively balanced state between upwelling and downwelling events throughout the observed period.Moreover, considerable variability in the intensity of upwelling can be observed yearly, with some years experiencing much stronger upwelling events than others (Figure 7).However, Figure 8 appears to be a weak seasonal pattern, with slightly more positive CUTI values during the spring and summer than autumn and winter.While there are fluctuations in the CUTI over the 26 years, there does not appear to be a consistent trend towards increasing or decreasing upwelling intensity (Figure 7) or no exceptionally strong upwelling or downwelling events observed in the time series (Figure 8).Year-to-year fluctuations in upwelling intensity and frequency, as revealed by the CUTI, underscore the dynamic nature of this process in the nearshore zone.These fluctuations, coupled with the pronounced seasonal variations in upwelling, highlight the significant influence of seasonal factors on driving upwelling events in this region.The offshore area of Constanta experiences a relatively moderate and balanced upwelling regime, with weak seasonal variations and notable interannual fluctuations.This information is valuable for understanding the oceanographic processes and potential impacts on marine ecosystems in this region.
Figure 9 shows that the SST exhibits a clear seasonal pattern, with higher temperatures in the summer months and lower temperatures in the winter months.This pattern is consistent for both nearshore and offshore locations.The nearshore SST tends to be slightly higher than the offshore SST, particularly during the warmer months.This could be due to the influence of shallower waters and land-based heat sources.The CUTI shows negative values during periods of upwelling.These events appear to occur more frequently during spring and summer, corresponding with periods of higher SST.There does not appear to be a significant long-term trend in either SST or the frequency and intensity of upwelling events over the time period covered by the graph.Figure 9 illustrates the complex interplay between SST, upwelling, and seasonality in the study area.This information is crucial for understanding the region's marine ecosystem dynamics, as upwelling plays a critical role in nutrient supply and primary productivity.
area.This information is crucial for understanding the region's marine ecosystem dynamics, as upwelling plays a critical role in nutrient supply and primary productivity.Figure 10 illustrates the time series of three variables: CUTI ECTA for nearshore and offshore waters and the calculated CUTI using sea surface temperature (SST) data.Figure 10 spans from January 1993 to September 2019 and shows the monthly variability of these indices.The nearshore values generally exhibit a larger magnitude than the offshore values, suggesting that upwelling is more intense in the nearshore region.Moreover, there is a clear seasonal pattern in all three variables, with the strongest upwelling events occurring during the spring and summer months (April to September).These three variables appear to be highly correlated, indicating that they all capture the same underlying phenomenon of upwelling.Furthermore, there is considerable interannual area.This information is crucial for understanding the region's marine ecosystem dynamics, as upwelling plays a critical role in nutrient supply and primary productivity.Figure 10 illustrates the time series of three variables: CUTI ECTA for nearshore and offshore waters and the calculated CUTI using sea surface temperature (SST) data.Figure 10 spans from January 1993 to September 2019 and shows the monthly variability of these indices.The nearshore values generally exhibit a larger magnitude than the offshore values, suggesting that upwelling is more intense in the nearshore region.Moreover, there is a clear seasonal pattern in all three variables, with the strongest upwelling events occurring during the spring and summer months (April to September).These three variables appear to be highly correlated, indicating that they all capture the same underlying phenomenon of upwelling.Furthermore, there is considerable interannual area.This information is crucial for understanding the region's marine ecosystem dynamics, as upwelling plays a critical role in nutrient supply and primary productivity.Figure 10 illustrates the time series of three variables: CUTI ECTA for nearshore and offshore waters and the calculated CUTI using sea surface temperature (SST) data.Figure 10 spans from January 1993 to September 2019 and shows the monthly variability of these indices.The nearshore values generally exhibit a larger magnitude than the offshore values, suggesting that upwelling is more intense in the nearshore region.Moreover, there is a clear seasonal pattern in all three variables, with the strongest upwelling events occurring during the spring and summer months (April to September).These three variables appear to be highly correlated, indicating that they all capture the same underlying phenomenon of upwelling.Furthermore, there is considerable interannual Figure 10 illustrates the time series of three variables: CUTI ECTA for nearshore and offshore waters and the calculated CUTI using sea surface temperature (SST) data.Figure 10 spans from January 1993 to September 2019 and shows the monthly variability of these indices.The nearshore values generally exhibit a larger magnitude than the offshore values, suggesting that upwelling is more intense in the nearshore region.Moreover, there is a clear seasonal pattern in all three variables, with the strongest upwelling events occurring during the spring and summer months (April to September).These three variables appear to be highly correlated, indicating that they all capture the same underlying phenomenon of upwelling.Furthermore, there is considerable interannual variability in the intensity of upwelling, with some years experiencing much stronger events than others.Overall, Figure 10 provides valuable insights into the spatial and temporal dynamics of upwelling in the study area, highlighting the importance of both nearshore and offshore processes and the significant role of seasonality in driving upwelling events.variability in the intensity of upwelling, with some years experiencing much stronger events than others.Overall, Figure 10 provides valuable insights into the spatial and temporal dynamics of upwelling in the study area, highlighting the importance of both nearshore and offshore processes and the significant role of seasonality in driving upwelling events.For the evaluation, homogeneity analysis (the Pettitt test, Table 1) was applied to the sea surface temperature (SST), CUTI (m 3 /s/100 m coastline) and CUTI SST (°C).Widely used in hydrological studies, the Pettitt test detects the significant changes in the mean of a time series, with performance analysis using the Monte Carlo simulation.The change point was identified as being in June 2010 and May 2014 (as shown in Figure 10) with a fixed significance level of α = 0.05.The seasonal variability of SST varies strongly by longitude depending on the upwelling phenomena intensity and the energy exchange between the sea and the atmosphere.The climatological SST displays a seasonality that varies sinusoidally for the Constanta transect (Figure 10).The CUTI SST index describes the seasonal variability of the upwelling intensity for the north-western Black Sea system from 28.6 to 32.85 East transect.For the evaluation, homogeneity analysis (the Pettitt test, Table 1) was applied to the sea surface temperature (SST), CUTI (m 3 /s/100 m coastline) and CUTI SST ( • C).Widely used in hydrological studies, the Pettitt test detects the significant changes in the mean of a time series, with performance analysis using the Monte Carlo simulation.The change point was identified as being in June 2010 and May 2014 (as shown in Figure 10) with a fixed significance level of α = 0.05.The seasonal variability of SST varies strongly by longitude depending on the upwelling phenomena intensity and the energy exchange between the sea and the atmosphere.The climatological SST displays a seasonality that varies sinusoidally for the Constanta transect (Figure 10).The CUTI SST index describes the seasonal variability of the upwelling intensity for the north-western Black Sea system from 28.6 to 32.85 East transect.
Figure 11 shows considerable interannual variability in the CUTI values, indicating that the intensity and frequency of upwelling events varied from year to year at the offshore Mangalia (MG6) station from January 1993 to March 2020.The CUTI values mostly fluctuate around zero, indicating a relatively balanced state between upwelling and downwelling events throughout the observed period.Moreover, there seems to be a slight tendency for upwelling to be more pronounced during the spring and summer months, as indicated by more frequent positive CUTI values during those periods, but no exceptionally strong upwelling or downwelling events can be observed in the time series.Figure 11 suggests that the MG6 station experiences a relatively moderate and balanced upwelling regime, with slight seasonal variations and notable interannual fluctuations.Figure 11 shows considerable interannual variability in the CUTI values, indicating that the intensity and frequency of upwelling events varied from year to year at the offshore Mangalia (MG6) station from January 1993 to March 2020.The CUTI values mostly fluctuate around zero, indicating a relatively balanced state between upwelling and downwelling events throughout the observed period.Moreover, there seems to be a slight tendency for upwelling to be more pronounced during the spring and summer months, as indicated by more frequent positive CUTI values during those periods, but no exceptionally strong upwelling or downwelling events can be observed in the time series.Figure 11 suggests that the MG6 station experiences a relatively moderate and balanced upwelling regime, with slight seasonal variations and notable interannual fluctuations.Figure 12 illustrates the comparison between the number of observed in situ upwelling events and the number of upwelling events predicted by the Coastal Upwelling Transport Index (CUTI) model over the period 2000-2019.The x-axis represents the months from May to September for each year, while the y-axis indicates the number of events, ranging from 0 to 3.5.Both in situ and CUTI-predicted upwelling events show a strong seasonal pattern, with the majority of events occurring during the summer months (June-August).This is consistent with the typical upwelling season in the Black Sea.In most cases, the number of in situ and CUTI-predicted events are in close agreement, Figure 12 illustrates the comparison between the number of observed in situ upwelling events and the number of upwelling events predicted by the Coastal Upwelling Transport Index (CUTI) model over the period 2000-2019.The x-axis represents the months from May to September for each year, while the y-axis indicates the number of events, ranging from 0 to 3.5.Both in situ and CUTI-predicted upwelling events show a strong seasonal pattern, with the majority of events occurring during the summer months (June-August).This is consistent with the typical upwelling season in the Black Sea.In most cases, the number of in situ and CUTI-predicted events are in close agreement, suggesting that the CUTI model is a reliable tool for predicting upwelling occurrences in this region.However, there are a few instances where the two datasets diverge, particularly in 2009 and 2018, where the CUTI model seems to overpredict the number of upwelling events compared to the in situ observations.suggesting that the CUTI model is a reliable tool for predicting upwelling occurrences in this region.However, there are a few instances where the two datasets diverge, particularly in 2009 and 2018, where the CUTI model seems to overpredict the number of upwelling events compared to the in situ observations.Overall, Figure 12 highlights the utility of the CUTI model in predicting upwelling events in the Black Sea while also emphasising the need for continuous validation with in situ data to improve the model's accuracy.Additionally, the graph reveals the distinct seasonal pattern of upwelling in this region, providing valuable information for understanding the region's marine ecosystem dynamics.Overall, Figure 12 highlights the utility of the CUTI model in predicting upwelling events in the Black Sea while also emphasising the need for continuous validation with in situ data to improve the model's accuracy.Additionally, the graph reveals the distinct seasonal pattern of upwelling in this region, providing valuable information for understanding the region's marine ecosystem dynamics.
Figure 13 displays the time series of two variables, CUTI ECTA1 (Coastal Upwelling Transport Index for the nearest coastal station) and BEUTI (Biologically Effective Upwelling Transport Index), from January 1993 to December 2019.The x-axis represents time (in months and years), and the y-axis shows the magnitude of both indices.Using Formula (2), we can calculate the BEUTI index that quantifies the amount of nitrate brought to the surface by upwelling.

Discussion
The Black Sea is an inland sea located in southeastern Europe, bordered by six countries: Bulgaria, Romania, Ukraine, Russia, Georgia, and Turkey.Unlike many coastal regions around the world that experience significant upwelling, the Black Sea does not have strong, continuous upwelling events due to its unique physical and hydrographic characteristics.Instead, the Black Sea is considered a "strongly stratified" sea, which means that it has distinct layers of water that do not mix as thoroughly as in a typical ocean.
Several notable features of upwelling and water circulation in the Black Sea basin can be briefly described:

−
Stratified Nature: The Black Sea is stratified into layers of water with different properties.The surface layer, known as the "epilimnion," is relatively warm and well-oxygenated.Below this layer is the "metalimnion," which acts as a transition zone.Beneath the metalimnion is the "hypolimnion," which is cold, stagnant, and lacks oxygen.The presence of the hypolimnion inhibits the significant vertical mixing of waters, which is characteristic of upwelling in other seas.-Limited Wind-Driven Upwelling: While the Black Sea does experience some winddriven vertical mixing, it is generally not strong enough to bring up significant nutrient-rich water from the depths to the surface regularly.The prevailing winds in the Black Sea are typically light, and the depth of the Black Sea is relatively shallow compared to large ocean basins.-Nutrient Distribution: The nutrients in the Black Sea are primarily supplied by river discharge, particularly from the Danube River, which carries nutrients from agricultural runoff.These nutrients are distributed throughout the water column but tend to accumulate in the anoxic (oxygen-depleted) hypolimnion.As a result, the surface waters are generally nutrient-poor compared to regions with frequent upwelling.-Anoxic Conditions: The Black Sea's unique stratification and limited vertical mixing can lead to the development of anoxic conditions in the deep hypolimnion.These anoxic conditions result in the accumulation of hydrogen sulphide, which can be toxic to marine life and have significant ecological implications.
In the realm of marine environmental conservation, comprehending upwelling is Both CUTI ECTA1 and BEUTI exhibit a distinct seasonal pattern, with positive values (upwelling) occurring predominantly in the spring and summer months, with negative values (downwelling) more common in autumn and winter.These two indices show a strong positive correlation, indicating that higher rates of upwelling are associated with increased nitrate transport to the surface.Also, there is considerable interannual variability in the magnitude of both indices, suggesting that the intensity of upwelling and associated nitrate transport vary from year to year.While there are fluctuations in both indices over time, there does not appear to be a clear long-term trend towards increasing or decreasing upwelling intensity in this location.
Overall, Figure 13 explains the dynamic nature of upwelling and its associated biological productivity in the Black Sea, emphasising the importance of considering both physical and biogeochemical processes when studying this phenomenon.

Discussion
The Black Sea is an inland sea located in southeastern Europe, bordered by six countries: Bulgaria, Romania, Ukraine, Russia, Georgia, and Turkey.Unlike many coastal regions around the world that experience significant upwelling, the Black Sea does not have strong, continuous upwelling events due to its unique physical and hydrographic characteristics.Instead, the Black Sea is considered a "strongly stratified" sea, which means that it has distinct layers of water that do not mix as thoroughly as in a typical ocean.
Several notable features of upwelling and water circulation in the Black Sea basin can be briefly described: -Stratified Nature: The Black Sea is stratified into layers of water with different properties.The surface layer, known as the "epilimnion," is relatively warm and welloxygenated.Below this layer is the "metalimnion," which acts as a transition zone.Beneath the metalimnion is the "hypolimnion," which is cold, stagnant, and lacks oxygen.The presence of the hypolimnion inhibits the significant vertical mixing of waters, which is characteristic of upwelling in other seas.
-Limited Wind-Driven Upwelling: the Black Sea does experience some winddriven vertical mixing, it is generally not strong enough to bring up significant nutrientrich water from the depths to the surface regularly.The prevailing winds in the Black Sea are typically light, and the depth of the Black Sea is relatively shallow compared to large ocean basins.-Nutrient Distribution: The nutrients in the Black Sea are primarily supplied by river discharge, particularly from the Danube River, which carries nutrients from agricultural runoff.These nutrients are distributed throughout the water column but tend to accumulate in the anoxic (oxygen-depleted) hypolimnion.As a result, the surface waters are generally nutrient-poor compared to regions with frequent upwelling.-Anoxic Conditions: The Black Sea's unique stratification and limited vertical mixing can lead to the development of anoxic conditions in the deep hypolimnion.These anoxic conditions result in the accumulation of hydrogen sulphide, which can be toxic to marine life and have significant ecological implications.
In the realm of marine environmental conservation, comprehending upwelling is essential for safeguarding this region against potential sources of pollution, whether they originate from industrial activities or ship waste [21].The occurrence, timing, and intensity of upwelling events are pivotal considerations in this context.Assessing the presence or absence of upwelling can be achieved through vertical temperature measurements of seawater using tools such as moorings or Argo floats [22].However, it is important to note that these methods are associated with high costs, localised data collection, and the inability to access real-time or near-real-time data due to the need for extensive post-processing.
Conversely, a more expedient and cost-effective approach widely adopted in the study of upwelling involves the examination of upwelling indicators [23].These indicators offer a swift and economical means of ascertaining the presence or absence of upwelling events in a given marine environment.
While the Coastal Upwelling Transport Index (CUTI) and Biologically Effective Upwelling Transport Index (BEUTI) have proven valuable in open ocean settings, their application to the semi-enclosed and highly stratified Black Sea necessitates a nuanced interpretation.Directly applying these indices, originally designed for large-scale, wind-driven upwelling systems, may not fully capture the complex upwelling dynamics influenced by factors such as riverine input and unique bathymetry.Therefore, validation with in situ observations of upwelling events, as presented in this study, is crucial for assessing the accuracy and applicability of CUTI and BEUTI within the Black Sea context.Despite these limitations, these indices remain valuable tools for understanding the relative variability and seasonal patterns of upwelling along the Romanian Black Sea coast.By focusing on comparative analyses and trend identification rather than absolute quantification, we can utilise CUTI and BEUTI to gain insights into the drivers and potential ecological impacts of upwelling in this unique marine environment.
The study of upwelling is, therefore, important in the NW Black Sea for advances in the knowledge of primary productivity processes.The present paper extends the statistical investigation of the Black Seas' upwelling events.It introduces the CUTI index covering the North-Western Black Sea shelf, based on free atmospheric and oceanographic data portals, such as Copernicus Marine Service data and ERDDAP (the Environmental Research Division's Data Access Program), validated with observations of in situ data.The results confirm the previously reported seasonality of upwelling in this region, with events primarily occurring between May and September ( Figures 3,4 and 12).The spatial variability of upwelling intensity is evident across different stations, with ECTA5 experiencing the most pronounced events (Figure 4a).This study's findings align with previous research highlighting the Black Sea's unique characteristics that influence upwelling patterns.The strong stratified nature of the sea, limited wind-driven upwelling due to weak winds and shallow depths, and the significant influence of riverine input, particularly from the Danube, all contribute to the observed dynamics.
The strong correlation between calculated CUTI and observed in situ upwelling events (Figure 12) validates the effectiveness of the CUTI as a cost-effective tool for assessing upwelling in the Black Sea.This is particularly important for a region where traditional methods like moorings or Argo floats are not feasible due to financial constraints.The successful application of CUTI in this study opens up avenues for its wider use in other regions with similar limitations.However, the discrepancies observed between in situ and CUTI-predicted events in some years (Figure 12) underscore the need for further refinement of the model.Incorporating additional parameters, such as local bathymetry and coastline configuration, can enhance the model's accuracy and predictive capabilities.
The implications of these findings extend beyond the Black Sea.This study contributes to a broader understanding of upwelling dynamics in semi-enclosed seas, which are often characterised by complex interactions between physical, chemical, and biological processes.Moreover, the demonstrated effectiveness of the CUTI model has potential applications for other regions facing similar challenges in monitoring and managing upwelling events.
To further enhance our understanding of upwelling processes and their role in shaping marine ecosystems, future research could delve into high-resolution modelling to investigate the fine-scale spatial and temporal dynamics of upwelling, potentially revealing underlying mechanisms and their ecological consequences.Additionally, coupling the CUTI model with biogeochemical models could help assess the impact of upwelling on nutrient transport, phytoplankton blooms, and overall ecosystem productivity in the Black Sea.Examining the potential effects of climate change on upwelling patterns in the Black Sea is crucial for predicting future changes in this region's marine ecosystem and developing appropriate adaptation strategies.Furthermore, comparative studies with other semi-enclosed seas could reveal commonalities and differences in upwelling mechanisms and their ecological implications, ultimately contributing to more effective conservation and management strategies in the face of ongoing environmental challenges.

Figure 1 .
Figure 1.The Danube's propagation altered summer freshwater flow patterns, as evidenced by th distribution of chlorophyll in surface waters.Reproduced with permission from [Kubryakov, A.A published by J. Mar.Syst.2018 [18].

Figure 1 .
Figure 1.The Danube's propagation altered summer freshwater flow patterns, as evidenced by the distribution of chlorophyll in surface waters.Reproduced with permission from [Kubryakov, A.A] published by J. Mar.Syst.2018 [18].

Figure 2 .
Figure 2. European Center for Medium-Range Weather Forecasts (ECMWFs) grid for the North-Western Black Sea shelf (Sfantu Gheorghe-abbreviation in the text; Station: Danube; Portitaabbreviation in the text; Station: PO; Constanta-abbreviation in the text; Station: ECTA; Mangaliaabbreviation in the text; Station: MG.Stations are numbered from the nearshore/near the coast from 1 to offshore station 6 or 7).

Figure 2 .
Figure 2. European Center for Medium-Range Weather Forecasts (ECMWFs) grid for the North-Western Black Sea shelf (Sfantu Gheorghe-abbreviation in the text; Station: Danube; Portita-abbreviation in the text; Station: PO; Constanta-abbreviation in the text; Station: ECTA; Mangalia-abbreviation in the text; Station: MG.Stations are numbered from the nearshore/near the coast from 1 to offshore station 6 or 7).

Figure 4 .
Figure 4. (a) Monthly and (b) yearly upwelling events from September for the 2000-2019 period at Constanta, Western Black Sea nearshore waters.

Figure 5
Figure 5 presents the monthly variability of the Ekman Transport Index (which is proportional to the wind stress) at different stations in the Black Sea throughout the year.The stations are represented by different coloured lines (ECTA and PO), and the x-axis shows the months.Most stations have clear seasonal patterns in upwelling and downwelling trends.Upwelling (positive values) is more common during spring and summer, while downwelling (negative values) is more prevalent during autumn and winter.The intensity of upwelling and downwelling varies significantly among the stations.Some stations experience stronger upwelling or downwelling than others.For example, the ECTA5 station tends to have stronger upwelling than other stations, while the DANUBE station experiences stronger downwelling than the rest.In addition, some stations, such as ECTA1 and ECTA3, exhibit considerable interannual variability in upwelling intensity.There seems to be a general trend of decreasing upwelling intensity.

Figure 5 .
Figure 5. Monthly variability in the Ekman Transport Index for the Western Black Sea.

Figure 6
Figure 6 illustrates the monthly variability of the Ekman Coastal Transport Index for East Constanta (ECTA) at eight different locations from the coast to deep waters (ECTA1

Figure 4 .
Figure 4. (a) Monthly and (b) yearly upwelling events from September for the 2000-2019 period at Constanta, Western Black Sea nearshore waters.

Figure 5
Figure 5 presents the monthly variability of the Ekman Transport Index (which is proportional to the wind stress) at different stations in the Black Sea throughout the year.The stations are represented by different coloured lines (ECTA and PO), and the x-axis shows the months.Most stations have clear seasonal patterns in upwelling and downwelling trends.Upwelling (positive values) is more common during spring and summer, while downwelling (negative values) is more prevalent during autumn and winter.The intensity of upwelling and downwelling varies significantly among the stations.Some stations experience stronger upwelling or downwelling than others.For example, the ECTA5 station tends to have stronger upwelling than other stations, while the DANUBE station experiences stronger downwelling than the rest.In addition, some stations, such as ECTA1 and ECTA3, exhibit considerable interannual variability in upwelling intensity.There seems to be a general trend of decreasing upwelling intensity.

Figure 5 .
Figure 5. Monthly variability in the Ekman Transport Index for the Western Black Sea.

Figure 6
Figure 6 illustrates the monthly variability of the Ekman Coastal Transport Index for East Constanta (ECTA) at eight different locations from the coast to deep waters (ECTA1

Figure 4 .
Figure 4. (a) Monthly and (b) yearly upwelling events from September for the 2000-2019 period at Constanta, Western Black Sea nearshore waters.

Figure 5
Figure 5 presents the monthly variability of the Ekman Transport Index (which is proportional to the wind stress) at different stations in the Black Sea throughout the year.The stations are represented by different coloured lines (ECTA and PO), and the x-axis shows the months.Most stations have clear seasonal patterns in upwelling and downwelling trends.Upwelling (positive values) is more common during spring and summer, while downwelling (negative values) is more prevalent during autumn and winter.The intensity of upwelling and downwelling varies significantly among the stations.Some stations experience stronger upwelling or downwelling than others.For example, the ECTA5 station tends to have stronger upwelling than other stations, while the DANUBE station experiences stronger downwelling than the rest.In addition, some stations, such as ECTA1 and ECTA3, exhibit considerable interannual variability in upwelling intensity.There seems to be a general trend of decreasing upwelling intensity.

15 Figure 3 .
Figure 3. Recorded upwelling [events no.] phenomena for May-September during the 2000-2019 period at Constanta, Western Black Sea nearshore waters.

Figure 4 .
Figure 4. (a) Monthly and (b) yearly upwelling events from September for the 2000-2019 period at Constanta, Western Black Sea nearshore waters.

Figure 5
Figure 5 presents the monthly variability of the Ekman Transport Index (which is proportional to the wind stress) at different stations in the Black Sea throughout the year.The stations are represented by different coloured lines (ECTA and PO), and the x-axis shows the months.Most stations have clear seasonal patterns in upwelling and downwelling trends.Upwelling (positive values) is more common during spring and summer, while downwelling (negative values) is more prevalent during autumn and winter.The intensity of upwelling and downwelling varies significantly among the stations.Some stations experience stronger upwelling or downwelling than others.For example, the ECTA5 station tends to have stronger upwelling than other stations, while the DANUBE station experiences stronger downwelling than the rest.In addition, some stations, such as ECTA1 and ECTA3, exhibit considerable interannual variability in upwelling intensity.There seems to be a general trend of decreasing upwelling intensity.

Figure 5 .
Figure 5. Monthly variability in the Ekman Transport Index for the Western Black Sea.

Figure 6 Figure 5 .
Figure 6 illustrates the monthly variability of the Ekman Coastal Transport Index for East Constanta (ECTA) at eight different locations from the coast to deep waters (ECTA1

Figure 6
Figure 6 illustrates the monthly variability of the Ekman Coastal Transport Index for East Constanta (ECTA) at eight different locations from the coast to deep waters (ECTA1 through ECTA8) throughout the year.Positive values indicate upwelling, while negative values indicate downwelling.Upwelling is generally more pronounced in the spring and summer (April to September), particularly at locations ECTA4, ECTA5, ECTA7, and ECTA8.Upwelling intensity varies significantly between locations, with ECTA5 and ECTA8 showing the strongest upwelling events.Some locations, such as ECTA1 and ECTA3, exhibit considerable interannual variability in upwelling intensity.Atmosphere 2024, 15, x FOR PEER REVIEW 7 of 15

Figure 6 .
Figure 6.Spatial and temporal dynamics of upwelling along the coast (Constanta), highlighting the complex interactions between wind patterns and coastal topography.

Figure 6 .
Figure 6.Spatial and temporal dynamics of upwelling along the coast (Constanta), highlighting the complex interactions between wind patterns and coastal topography.

Figure 7 .
Figure 7. Monthly variability in the Coastal Upwelling Transport Index (CUTI) near the shore of Constanta, Romania, from January 1993 to April 2019.

Figure 8 .
Figure 8. Coastal Upwelling Transport Index (CUTI) for the offshore region of Constanta, Romania, from January 1993 to December 2019.

Figure 9 .
Figure 9. Evolution of the sea surface temperature (SST) in both nearshore and offshore areas, along with the Coastal Upwelling Transport Index (CUTI) over time, likely in the Black Sea.

Figure 7 .
Figure 7. Monthly variability in the Coastal Upwelling Transport Index (CUTI) near the shore of Constanta, Romania, from January 1993 to April 2019.

Figure 7 .
Figure 7. Monthly variability in the Coastal Upwelling Transport Index (CUTI) near the shore of Constanta, Romania, from January 1993 to April 2019.

Figure 8 .
Figure 8. Coastal Upwelling Transport Index (CUTI) for the offshore region of Constanta, Romania, from January 1993 to December 2019.

Figure 9 .
Figure 9. Evolution of the sea surface temperature (SST) in both nearshore and offshore areas, along with the Coastal Upwelling Transport Index (CUTI) over time, likely in the Black Sea.

Figure 8 .
Figure 8. Coastal Upwelling Transport Index (CUTI) for the offshore region of Constanta, Romania, from January 1993 to December 2019.

Figure 7 .
Figure 7. Monthly variability in the Coastal Upwelling Transport Index (CUTI) near the shore of Constanta, Romania, from January 1993 to April 2019.

Figure 8 .
Figure 8. Coastal Upwelling Transport Index (CUTI) for the offshore region of Constanta, Romania, from January 1993 to December 2019.

Figure 9 .
Figure 9. Evolution of the sea surface temperature (SST) in both nearshore and offshore areas, along with the Coastal Upwelling Transport Index (CUTI) over time, likely in the Black Sea.

Figure 9 .
Figure 9. Evolution of the sea surface temperature (SST) in both nearshore and offshore areas, along with the Coastal Upwelling Transport Index (CUTI) over time, likely in the Black Sea.

Figure 10 .
Figure 10.CUTI for nearshore and offshore waters the for East Constanta transect (CUTI ECTA).

Figure 10 .
Figure 10.CUTI for nearshore and offshore waters the for East Constanta transect (CUTI ECTA).

Figure 12 .
Figure 12.Comparison between observed upwelling events (in situ) and those calculated CUTI for the 1993-2019 period at Constanta for the favourable upwelling month.

Figure 13
Figure 13 displays the time series of two variables, CUTI ECTA1 (Coastal Upwelling Transport Index for the nearest coastal station) and BEUTI (Biologically Effective Upwelling Transport Index), from January 1993 to December 2019.The x-axis represents time (in months and years), and the y-axis shows the magnitude of both indices.Using Formula (2), we can calculate the BEUTI index that quantifies the amount of nitrate

Figure 12 .
Figure 12.Comparison between observed upwelling events (in situ) and those calculated CUTI for the 1993-2019 period at Constanta for the favourable upwelling month.

Figure 13 .
Figure 13.CUTI ECTA1 (Coastal Upwelling Transport Index for the nearest coastal station) and BEUTI (Biologically Effective Upwelling Transport Index), from January 1993 to December 2019.