Assessment of the Hydrochemical Characteristics of the Carbon Observational Site ‘Carbon-Sakhalin’ (Aniva Bay, Sea of Okhotsk)

By inergency On Apr 5, 2024

[ad_1]

1. Introduction

Aniva Bay, situated in the northern part of La Perouse Strait, washes the southern coast of Sakhalin Island and belongs to the Sea of Okhotsk water basin (Figure 1). Aniva Bay is important for transport and fishery activities. Large hydrotechnical structures are located on the Bay’s coastline: the sea port of Korsakov (annual cargo turnover of >1.5 mln tonnes) and the oil and liquefied natural gas terminal in the port of Prigorodnoye (cargo turnover of >1.5 mln tons) [1], which are potential sources of marine waters technogenic pollution. Moreover, Salmon Bay (Bukhta Lososey), in the southwestern part of Aniva Bay, is exposed to pollutants through river runoff from the coastal cities of Aniva and Korsakov and the densely populated Susunai Valley. Littoral eelgrass meadows (Zostera, Phyllospadix) and subtidal kelps of brown seaweeds (Saccharina japonica) create areas of intense primary productivity and large biomass of commercial species of fish and invertebrates. The total stock of the most harvested species of brown algae Saccharina japonica is estimated here at more than 40 thousand tons [2]. The water area suitable for extensive mariculture is 177,455 ha [3].

Aniva Bay is a part of the Sea of Okhotsk with the widest range of variability of natural conditions and difficult ice conditions. It makes it impossible to conduct seasonal field observations of greenhouse gas absorption and constant monitoring of the carbon balance of the water area. The biological pump is the main part of the ocean’s carbon pump [4,5,6]. Many environmental factors need the phytoplankton blooms as the main force of the biological carbon pump (sufficient illumination, nutrient content, acceptable temperature, a certain regime of salinity and hydrodynamics, and much more). The range of environmental conditions is specific for each climatic zone and it determines the value of the primary productivity of the water area and carbon flux. Thus, the hydrological and hydrochemical conditions of the water area determine the net absorption of greenhouse gases by the marine ecosystem. Therefore, it is important to have a full-scale basis for modeling the annual cycle of carbon change in a particular ecosystem.

Aniva Bay is the marine part of the carbon observational site ‘Carbon-Sakhalin’ which is organized on the base of Sakhalin State University. ‘Carbon-Sakhalin’ is a part of the carbon polygons project, which is supported by the Government of the Russian Federation. This project is part of the national plan for adaptation to climate change [7].

The Program of Carbon Polygon in the Sakhalin Region was approved on 23 September 2021 at the Expert Council on Scientific Support for the Development of Carbon Balance Control Technologies of the Ministry of Education and Science of the Russian Federation.

The observational site ‘Carbon-Sakhalin’ is one of the first marine polygons in Russia. In 2022, the first all-Russian experiment to limit greenhouse gas emissions was launched in accordance with Federal Law No. 34-FZ, dated 6 March 2022, “On conducting an experiment to limit greenhouse gas emissions in certain subjects of the Russian Federation”. To solve the problem of estimating and reducing CO₂ emissions, as well as possible ways of sequestering carbon, the most studied water area of the Sea of Okhotsk was selected. Aniva Bay is located close to the scientific centers and has developed infrastructure that made it possible to ensure a successful start of the project.

The general project concept includes several aims:

–: Creating a network of carbon observational sites, which include terrestrial and marine areas in different climatic zones of Russia;
–: Investigation of the physical and biochemical processes that define the carbon cycle;
–: Measuring of emission and absorption of greenhouse gases (GHGs) using the standard equipment and methods;
–: Assessments of the regional biogeochemical carbon cycle and human impacts;
–: Investigation of the potential of GHG sequestration by natural ecosystems;
–: Development and testing of carbon balance control and technologies of GHG sequestration, and anthropogenic increase in the sequestration capacity of natural ecosystems within the framework of the “green” and “blue” economy.

Russia is still at the very beginning of the way to assess carbon emissions in the coastal zone of our seas. Research at the carbon observational site will certainly contribute to the sustainable development of the Sakhalin region.

The network of carbon polygons in Russia is created on a competitive basis and is funded by the Ministry of Science and Higher Education of the Russian Federation and businesses.

Financing is possible only after an assessment of the regional legislation in limiting greenhouse gas emissions, the regional system of carbon regulation, an inventory of all types of emissions and integral estimates of greenhouse gas emissions and absorption, the capabilities of scientific organizations, the complexity of the approach, and the presence of a variety of non-urbanized ecosystems.

The total greenhouse gas emissions in the Sakhalin Region were 12.333 million tons e-CO₂; total net greenhouse gas absorption was 11.068 million tons e-CO₂ (only by forest lands) in the Sakhalin Region in 2019 [8]. The total net greenhouse gas emissions in the Sakhalin Region amounted to 1.265 million tons of e-CO₂ in 2019 [8].

Marine ecosystems and blue carbon were not included in the accounting. Therefore, the study of greenhouse gas fluxes in coastal marine ecosystems and the assessment of net absorption and emissions for subsequent inclusion in the national plan adaptation and mitigation of Russia is our strategic mission. The adaptation of the world’s Ocean Carbon Dioxide Removal (CDR) technologies through the development of mariculture farms will help the Sakhalin region quickly achieve carbon neutrality.

A marine carbon site is a sea area where we carried out a complex of measurements of meteorological, hydrological, hydrochemical (HCh) (including fluxes of carbon dioxide in the ocean-atmosphere system), hydrobiological parameters in both field and remote mode. Several important points about the development of the network of carbon polygons in Russia are presented in [9].

The absorption of carbon by marine ecosystems, based on both physical and biological processes, is most active in natural eelgrass meadows and kelps of seaweed (including maricultural), foreshores, and areas with high phytoplankton density (estuaries, lagoons, upwelling zones, areas of organic matter accumulation, etc.). However, simultaneously with carbon absorption is its emission, the so-called ‘breath of the sea’. This breath consists of CO₂ emission through breathing hydrobionts, decomposition of organic matter by microorganisms, CO₂ emission from seawater during storms, ice melting, etc. [10,11,12,13,14,15,16]. The processes described above make the baseline of absorption, which is unique for every body of water.

The combination of natural and anthropogenic conditions makes Aniva Bay essential for studying the biogeochemical carbon cycle in the marine ecosystem. The Bay is shallow (up to 100 m), convenient for field research by boats, its ice regime has been studied, and there are abundant natural macrophyte populations, large littoral tidal flats, and a flux of large rivers with wastewater. An interesting point that complicates the carbon balance of this reservoir is the natural and artificial salmon populations, which additionally bring in a large amount of organic matter. Therefore, Aniva Bay was chosen as the marine part of the carbon measurements test areas.

It is planned to install several stationary oceanographic platforms with a set of meteorological and hydro-hydrochemical sensors, conduct regular field sampling of water, ice, sediments, plankton, and benthos, including the winter ice period to calculate the baseline of CO₂ absorption at the carbon site in Aniva Bay.

Another task of the ‘Carbon-Sakhalin’ site is to search for technologies for additional extraction of GHGs by the marine ecosystem. One of the most natural-like technologies of GHG sequestration, prevention of acidification and eutrophication in marine environments is the creation of seaweed farms, mussel farms, and integrated multi-trophic aquaculture farms. Such farms, with organisms of different trophic levels (autotrophic and heterotrophic), can significantly increase carbon absorption and decrease the concentration of biogenic elements to prevent eutrophication [13,17,18,19,20,21,22,23]. There is active development of algorithms to assess the intake capacity of bodies of water for such farms, risks of eutrophication, and interaction with natural populations [24,25,26,27]. Large kelp farms may decrease phytoplankton productivity due to the decrease in available nitrogen and phosphorus forms, which impacts the overall level of biological productivity in the area [28,29].

So, studies of the hydrochemical (nutrient) regime in Aniva Bay will identify suitable areas for the deployment of farms with essential hydrological and hydrochemical conditions and will help us plan the optimal number of carbon mariculture seaweed–bivalve farms. Management strategies for carbon farm deployment are aimed at limiting the area for farms, recommendations for placement sites and species composition. It is very important to improve the health and productivity of the Aniva Bay ecosystem and prevent eutrophication and damage of local ecosystems and biodiversity.

In general, modern hydrochemical data are the basis for assessments of the absorption capacity of marine areas for GHGs and help the regional government to make science-based decisions regarding economic and mitigation strategies for climate change in the region.

Understanding the hydrophysical variability is of utter importance for the research of processes of organic matter transformation and GHG emissions [30,31]. Since biogeochemical models [32,33,34,35] are used to assess GHG emissions on site, the hydrochemical characteristics of water areas are set as initial and boundary conditions. Taking into account climate change, it is essential to understand the seasonal and interannual variability of hydrochemical parameters of bay water, the tendencies of its changes, and the evaluation of its current state. Monitoring of various parameters at carbon sites using autonomous devices is very prospective research [36].

Aniva Bay is part of the Sea of Okhotsk and, at the same time, is influenced by the Sea of Japan, both of them play significant climate effects on the northern Pacific region. The interaction of subarctic and subtropical water masses in Aniva Bay creates a unique pattern of greenhouse gas flux from the northern Pacific region and back.

Therefore, the analysis of hydrochemical parameters in Aniva Bay provides information about the carbon cycle and fluxes of greenhouse gas emissions in part of the northern Pacific region.

In 1975, the National System for Environmental Observations and Control (NSEOC) was initiated in the Russian Far East. Within this system, seawater quality was monitored in Far Eastern seas and Aniva Bay, and this monitoring has been carried out with various regularity up to today [37].

The first fundamental overview of studies on the Sea of Okhotsk hydrochemistry was presented in a monograph Hydrometeorology and Hydrochemistry of Seas in the “Seas of the USSR” series [38]. However, the coarse spatial scale of the research (1 × 1°) did not permit us to provide details of the hydrochemical regime of Aniva Bay. The next large-scale review of hydrochemical research, conducted within the “Ecology of Russian Seas” project [39], was devoted to the Sea of Okhotsk and the Pacific Ocean water territories around the Kuril Straits, so Aniva Bay was not considered separately.

In 1998, the first version of the Hydrological and Hydrochemical Atlas of the Sakhalin shelf was created (“ATLAS” software № 930078) [40]. These applications permitted to plot maps of the seasonal distribution of hydrological sections in a 100-mile zone adjacent to Sakhalin, as well as the vertical distribution of parameters for intermediate oceanological stations in eight selected regions [41]. The program and paper version of the “ATLAS” has been widely used in the last decades by research institutions in the Far East as a practical reference guide during oceanological, hydrobiological, and numerical studies for the assessment of regime parameters of the marine environment in the Sakhalin shelf zone.

The first overall seasonal spatial distribution of hydrological and hydrochemical characteristics in Aniva Bay, based on the data from the “ATLAS” program, was calculated in a study guide [42] aimed at undergraduate students, PhD students, and researchers of Sakhalin State University. A more detailed edition, supplemented with results of the marine environment state modeling, was published in 2020 [43].

An analysis of separate characteristics of the hydrochemical regime of Aniva Bay is presented in a series of publications after hydrobiological studies were conducted, typically on micropolygons, in order to solve practical tasks [44,45,46]. Various aspects of the hydrochemical regime in Aniva Bay were discussed after the expedition results in works by the Sakhalin Scientific Research Institute of the Fish Facilities and Oceanography (SakhNIRO) [47,48,49].

There are a few significant articles devoted to the hydrochemical regime of the Sea of Okhotsk as a GHG absorber and the only source of intermediate water circulation in the northern Pacific Ocean. The available data from research vessels and satellites were summarized and they showed decreased water salinity, decreased dissolved oxygen, and increased transparency from a long-term perspective [50,51,52,53].

Biohydrochemical processes in the Sea of Okhotsk are still understudied, especially during ice season. It is known that cold conditions create a community of low-thermophilic nitrogen-fixing and photosynthetic microorganisms [54]. Concentrations of chlorophyll-a in the Sea of Okhotsk ice may reach 3.5 mg/m², which indicates active vegetation of sea-ice microcoenosis [55] resulting in the formation of a specific oligotrophic under-ice layer of water [56]. Research demonstrated that during ice melt in the Sea of Okhotsk, sea-ice microcoenosis actively absorbs carbon dioxide. After the ice melts, these organisms take part in the formation of spring phytoplankton outbursts [57]. In the work by [58], the inflow of anthropogenic CO₂ (CO₂-antr) is calculated for 1993–2000, which made up 10 g C/m² a year. Accumulation of CO₂-antr for the same period for the whole Sea of Okhotsk is estimated be 0.78 Tg (1 Tg = 10¹² g) C yr⁻¹ over the entire Sea of Okhotsk. A few studies are devoted to the general analysis of hydrochemical parameters in other parts of the Sea of Okhotsk [59,60,61,62], but they do not consider Aniva Bay. Moreover, a significant climate effect was established from the Sea of Okhotsk on the northern Pacific region from the Sea of Okhotsk [63]. For example, the annual discharge of dissolved organic carbon from the Sea of Okhotsk into intermediate layer waters of the northern Pacific Ocean made up 68–72 Tg C yr⁻¹ and dissolved organic nitrogen of 5.4 Tg N yr⁻¹ [64].

Of course, volcanic eruptions can have a noticeable effect on the chemical composition of the atmosphere and ocean [65,66]. It is also necessary to pay attention to the work related to investigations about the coupled CO₂ diffusion–reaction process influencing CO₂ transport and the evolution of carbonate precipitation in diffusion-limited zones [67].

Thus, a detailed analysis of hydrochemical parameters in Aniva Bay is lacking and we are trying here to fill this gap. In this work, we supplement and reconsider the previously published preliminary results devoted to the HCh regime in Aniva Bay waters [41,42,43].

The principle motivation for this research is the need to have an understanding of spatial and temporal variability of HCh parameters for the upcoming work on modeling the carbon cycle and gas emissions. The aim of this study is to review all the available data and publications devoted to the analysis of spatial and temporal variability of oxygen, pH, inorganic phosphorus, inorganic nitrogen, and silica in Aniva Bay. Based on the available dataset of in situ measurements over 1975–1994, we try to reconstruct the annual cycle of HCh parameters and render maps of their spatial variability.

The next important purpose of this article is to provide a wide review of studies on hydrochemistry and modeling of the Sea of Okhotsk. Unfortunately, most of the studies have been published in highly specialized journals and books in the Russian language and access is therefore limited to Russian speakers.

This article consists of 5 sections. The Introduction focuses on the urgency of the topic, our motivation and aims. Section 2 gives a review of previously published research papers about HCh parameters in the Sea of Okhotsk and Aniva Bay. Section 3 describes the initial data and methods of their processing. Section 4 presents the results of the analysis of HCh parameters seasonal variability in Aniva Bay. Section 5 is devoted to discussion, and Section 6 provides our conclusions.

2. Review of Hydrochemical Regime of the Sea of Okhotsk and Aniva Bay Waters

In [38], initial data for the hydrochemical analysis of the Sea of Okhotsk waters were taken from all the available research vessel deep water measurements obtained from 51,607 stations during 1930–1988. It is worth noting that water temperature was recorded at all oceanographic stations, while dissolved oxygen was measured only at 9236 stations, with inorganic forms of P, N, and Si measured at 1478, 1637, and 2239 stations, respectively. All the available data, after clearing, were sorted in 1-degree squares and then these datasets were divided into months. For every square, statistics were calculated (mean, maximum, minimum values, and mean-squared deviations for all elements at all standard sample depths), which were related to its center. From the obtained numerical data, maps were plotted for the spatial and temporal distribution of parameters across the entire Sea of Okhotsk water territory. Vertical distribution graphs were also plotted for separate squares.

The conducted studies permitted the formulation of the basic features of the HCh regime of the Sea of Okhotsk in general. Thus, we confirmed the consistency of the main salt composition of the water in different parts of the Sea of Okhotsk on the surface and in the water column. A similarity was noted between the main salt composition of seawater and the main salt composition of waters in the adjacent part of the Pacific Ocean, which has the overall tendency of vertical distribution for nearly all elements of the salt composition—their concentrations increase with depth. It was also shown that it makes sense to divide and analyze the Sea of Okhotsk water column into several layers depending on the intensity of thermodynamical, chemical, and biological processes taking place there. These layers are surface, intermediate, and deepwater ones.

The surface layer clearly demonstrates seasonal changes in HCh parameters (apart from dynamically active sea zones) down to the depths of 150–250 m. In spring and summer, most of the water territory of the sea exhibits subsurface oxygen (up to 14.5 mg/L) and phosphorus (P-PO₄) maximums, as well as subsurface minimums of P-PO₄ and silica, which coincide with the oxygen maximum. This layer also has the highest values of hydrogen parameters (pH), up to 8.1–8.4. The depth of oxygen compensation varies from 30 to 70 m. The concentration of nitrite-nitrogen (N-NO₂) does not exceed 0.001–0.0015 Nmg/L (in dynamically active zones, it rises up to 0.002–0.004 Nmg/L), P-PO₄—0.09 Pmg/L, and silica—1.7 Simg/L.

It was found that there are no seasonal changes in values of HCh parameters in the intermediate layer. The lower boundary of this layer is the oxygen minimum layer where its concentration drops until 1.1–1.8 mg/L (saturation 12–18%). The thickness of this oxygen minimum layer varies from 200–500 m in most parts of the sea to 1200–1500 m in the Kuril Straits in zones with active mixing. Remarkable characteristics are lower oxygen concentration, minimum pH (down to 7.5–7.8) and nitrite-nitrogen values, and maximum alkaline values (up to 2.50–2.70 mmol/L). Concentrations of inorganic phosphorus and silica reach their maximums at the lower boundary of the layer (up to 0.08–0.14 Pmg/L and 4.0–6.0 Simg/L, respectively).

Characteristic features of HCh variability in the deepwater layer are increased oxygen concentrations up to 2.5–3.0 mg/L across the entire basin (saturation up to 22–27%), decreased silica and stable P-PO₄, pH, and alkaline. Concentrations of nitrite-nitrogen drop until analytical zero here [38].

A series of studies of comprehensive oceanological research of the Sea of Okhotsk ecosystem was carried out by the laboratory of marine ecology of the Russian Research Institute of Fishery and Oceanography (VNIRO) in 1990–1994 [39]. A joint monograph presents the results of full-scale expeditions in the Sea of Okhotsk in 1990–1994 carried out by the laboratory of Marine Ecology VNIRO, together with researchers from the Pacific branch of the Federal State Budget Scientific Institution, the Russian Federal Research Institute of Fisheries and Oceanography (TINRO) and the Shirshov Institute of Oceanology of Russian Academy of Sciences. The indisputable value of these studies is the systematic approach to the analysis of the results obtained during five expeditions with the use of modern equipment. The conclusions in this work regarding hydrochemistry, based on high-accuracy data from flow automatic analyzers, permitted us to check and verify currency fields with respect to HCh parameters and to assess phytoplankton production with respect to a decrease in dissolved inorganic nitrogen (inorganic nutrients) in the eutrophic layer over the period from the beginning of its vegetation to the moment of research.

New methods of registration of organic and mineral compounds of phosphorus, nitrogen, and silica made it possible to finalize the balance of various inorganic nutrient forms, single out the “new” primary production based on the inorganic nutrients recycling, which does not reach higher trophic levels, i.e., does not create a biological ground for a fishery. A stable maximum of chlorophyll-a concentration above the thermocline was established during a detailed study of the vertical structure of phytopigments. In most cases, this maximum is much higher than the upper extremum at depths of 10–20 m. It is exactly here where the “new” primary production is created in summer when the main inorganic nutrients supply above the thermocline has already been assimilated. In the second half of summer, there is a change in phytocoenosis around the lower maximum of chlorophyll. Diatomic plankton slowly decreases and, in some cases, it is substituted for peridinials.

Instrumental measurements confirmed the fact that inorganic nutrients content rises due to autumn–winter convective mixing. The concentration of inorganic nutrients is quite high across the entire Sea of Okhotsk: 0.0558–0.0744 mgP/L for inorganic phosphorus, 0.210–0.420 mgN/L in total for all forms of inorganic nitrogen, and 0.840–1.680 mgSi/L for silica. Maximum concentrations are registered in the areas of continental slopes of the Sakhalin and the Kamchatka, where they match the anti-cyclonic eddies. The deepwater layer of the sea demonstrates the least variety of inorganic nutrients, although shallow-water shelf areas also have lower concentrations. In summer, the content of inorganic nutrients drops drastically compared to winter and, unlike the pre-vegetational period with its lack of nitrogen, shows a deficiency of silica. The lack of inorganic forms of nitrogen and phosphorus is mostly compensated by their regeneration during the decomposition of organic matter and metabolic processes of marine hydrobionts, as well as the ability of phytoplankton to consume nitrogen and phosphorus organic compounds. Meanwhile, silica is regenerated rather slowly; hence, it is relatively deficient compared to nitrogen and phosphorus [39]. Unfortunately, the water territory of Aniva Bay was not included in such research.

In [42,43], one can find the most comprehensive analysis of HCh regime of the La Perouse Strait waters, conducted on the basis of mean average calculations of O₂, pH, P-PO₄, N-NO₂, and Si with the “ATLAS” software for 1975–1994 [40,68]. However, these data provide only partial information on the HCh regime in Aniva Bay. The characteristics of Aniva Bay obtained from this application were used as input parameters for the modeling of hydrological, HCh parameters, and transformation of petroleum hydrocarbons in the bay [43,69,70,71,72]..

Seasonal research of HCh parameters (salinity, contents of oxygen, ammonium ions, nitrate nitrogen, inorganic phosphorus, silica, and chlorophyll-a) was conducted in the southern part on Aniva Bay and the adjacent La Perouse Strait water territory in 1996–1997 [73]. The analysis of rather patchy information on the seasonal inorganic nutrients variation demonstrates low values of inorganic phosphorus—down to 0.0803 mgP/L. For ammonium nitrogen (N-NH₄) the concentrations varied between 0 to 0.0676 mgN/L, for nitrate nitrogen (N-NO₃)—from 0 to 0.2366 mgN/L, for silica—from 0.0224 to 0.9436 Si/L. Maximum concentrations of inorganic forms of N and P were recorded in June, and Si in March.

Complex seasonal research of hydrochemical and hydrobiological parameters in Aniva Bay was carried out on standard cross-sections by SakhNIRO in 2001–2013 (five sets of measurements a year in 2001, 2002, 2009, and three sets in 2013) [44,45,47,49,73]. In hydrological spring 2001–2002 (April–June), the concentrations of all inorganic nutrients above the thermocline were low due to high amounts of phytoplankton. However, below the thermocline these numbers were quite high, reaching 0.0088 mgN/L for nitrite nitrogen, 0.5418 mgN/L for nitrate nitrogen, 0.0701 mgP/L for inorganic phosphorus, and 1.1312 mgSi/L for silica. The same inorganic nutrients distribution (minimum above the thermocline and maximum below it) was observed in summer (August). In autumn (October–November) low values of inorganic nutrients in the surface water layer remained because of the phytoplankton bloom, while the inorganic nutrients reserve in deepwater layers remained unchanged. In 2009 and 2013, the peculiarities of vertical inorganic nutrients distribution in spring demonstrated poorer upper layer. The concentrations were as follows: N-NO₂—0–0.0103 mgN/L; N-NO₃—0–0.4551 mgN/L; P-PO₄—0–0.1734 mgP/L; and Si—0.0111–1.6382 mgSi/L. In summer (July–first third of September) the inorganic nutrients vertical distribution was similar. As for the values of concentrations, inorganic phosphorus dropped significantly, while silica increased because of the shift in the type of phytoplankton from diatomic to peridinial. At the same time, the values for inorganic nitrogen were similar to those in spring: N-NO₂—0–0.0101 mgN/L; N-NO₃—0–0.3792 mgN/L; P-PO₄—0–0.1002 mgP/L; and Si –0.0437–1.8161 mgSi/L. In autumn (October–November) the maximum concentrations of nitrite and nitrate nitrogen remained almost unchanged (maximum 9.8 mgr N-NO₂/L and 382.2 mgr N-NO₃/L), while the values for inorganic phosphorus and silica rose considerably (up to 0.1967 mgP/L and 2.5875 mgSi/L).

Thus, we can conclude that the most complete hydrochemical research of Aniva Bay was conducted between 1975 and 1994. In the last 30 years, the studies have been inconsistent, which makes it impossible to assess carbon balance in the water basin correctly based on field measurements.

5. Discussion

As was shown in [72], the Aniva Bay basin is occupied by a subarctic water mass. In the warm season, it is represented by a subsurface Sea of Okhotsk water mass of summer modification and the very Sea of Okhotsk water mass (SOWM). In the cold season, it is only SOWM. The La Perouse Strait basin, which includes Aniva Bay, contains the surface Pacific water mass in the form of the Soya Current all year round (part of the subtropical water mass of the Sea of Japan). The subarctic water structure occupies from 60% (in May) to 90% (in November) of the strait. Between the mentioned structures, there is a distinct border throughout the year, but it does not stop the interpenetration of water masses under certain hydrological conditions, which makes the hydrological and HCh regime of both La Perouse Strait and Aniva Bay much more complicated.

Seemingly, there is enough data for the analysis of the HCh regime in Aniva Bay (we used data from more than 600 stations for oxygen and about 300 for biogenic elements); however, these measurements are spaced in time very inhomogeneously. Oceanographic expeditions were in different years on different dates, so it does not seem possible to reconstruct the annual cycle of parameters correctly with statistical methods. The number of measurements also changes greatly from year to year (from 1981 to 1987, there were fewer observations than in other years) (Figure 17). Hence, it is not possible to obtain trustworthy climatic norms for the studied parameters for the time intervals recommended by the World Meteorological Organization [83], since an uneven number of observations in separate years influences the mean values of calculated parameters.

The question of statistical analysis of raw data, especially if they are distributed unevenly in space and time, is always controversial. However, the method chosen here of objective interpolation with a running average proves its consistency in the example of correct results of the combined hydroecological CNPSi model and the oceanographic model of Bergen University [84,85,86].

In the future, it is assumed that simulations should be carried out within the works on carbon polygon, with reviewed input data and three models applied: a new version of “ATLAS”-2018 [68], a new version of CNPSi-model [87], and a more actual, compared to the oceanographic model of Princeton University [88].

Still, the lack of field observations in the cold season does not permit us to evaluate the “sea breath” during the ice cover settling and melting, which has been seriously underestimated just recently [56,57]. The analysis of all-year-round variations of inorganic nutrients and NF concentrations is crucial for the understanding of biohydrochemical patterns of their transformation, specific features of spatial distribution and the formation of plankton communities that define the CO₂ absorption. That is why it is essential to assess the hydrochemical water conditions as a basis for the experimental carbon site and the evaluation of carbon emission and absorption. Inorganic nutrient concentrations can even be used to identify water masses, which has not been completed yet for the studied water basin.

Also important is the issue of the increasing anthropogenic impact that occurs as a result of the industrial development of the Sakhalin region. The construction of a liquefied natural gas terminal, the development of related infrastructure, and a large number of tourists—all these factors will lead to a change in the balance of the ecosystem.

It is also necessary to mention important works related to the use of various nanorods for the conversion of CO₂, which subsequently leads to the restoration of the environment [89,90].

Our results showed that the availability of nutrients is sufficient for the mass placement of mariculture marine farms in the coastal zone of Aniva Bay. It is especially important we are shown that the areas with the maximum nutrient run-off (mineral phosphorus and nitrogen) are Salmon Bay and the east coast. It is a better site for sanitary seaweed plantations for wastewater treatment or rapid absorption of CO₂ from the atmosphere [91]. Our calculations show that when growing brown algae Saccharina japonica in Aniva Bay for the purpose of decarbonization, it is possible to additionally sequestrate up to 49.5 thousand tons of C [86]. Successful cultivation of algae, their timely removal from the environment (before the destruction of the thallus) and the production of biochar from them with further industrial use [92,93,94] will increase the absorption of CO₂ from the atmosphere and have a positive impact on the carbon balance of the entire region. Experimental measurements of CO₂ fluxes in the area of mariculture farms have shown a noticeable photosynthetic absorption of CO₂ from the environment [95,96,97,98].

The authors are not aware of any works with analysis of long-run seasonal variations of concentrations of organic and inorganic nutrient fractions in the Sakhalin shelf water basins. When hydrochemical research of the World Ocean was on the rise in the 1970s due to a new instrumental base, Aniva Bay was an internal body of water. Then, the introduction of a 200-mile economic zone in 1979 made it practically inaccessible for international researchers. All studies conducted by the manages of the Sakhalin shelf projects were strictly formalized and the research results are their property. So, in our calculations, we used only (mostly) open access national observation materials. Thus, the analysis of instrumental observations data and numerical simulation results presented here can be recommended for the baseline evaluation and calculation of carbon emissions on the carbon site in Aniva Bay.

The authors understand that field data from the analysis of hydrochemical parameters have limitations and uncertainties in assessing greenhouse gas emissions from marine ecosystems.

The greatest uncertainty in estimating fluxes GHGs by field research data comes from sampling conditions and spatial–temporal variability of hydrochemical parameters: sampling time, depth, surface insolation, waves, tidal phase, time after rain or storm, and other natural parameters. Sampling difficulty from large areas in a short time period and sampling inability during the ice period also limit the use of direct field data to calculate the annual balance of GHGs. Some greenhouse gases, like methane, can be challenging to quantify accurately using only hydrochemical data (bacterioplankton investigation is needed).

Climate change has made it impossible to use a long historical series of direct hydrochemical observations to assess the current state of the carbon balance of marine ecosystems.

We would like to show the several alternative approaches to flux GHG estimation used by scientific groups:

–: In situ gas sensors and continuous monitoring systems (including eddy covariance techniques) on ocean platforms, vessels, and ferries can provide real-time data on gas emissions more effectively than discrete sampling;
–: Biomolecular marker techniques like DNA-based analyses or genetic markers can identify specific microbial communities responsible for greenhouse gas production, offering insights into potential emission sources;
–: Isotope tracing methods can trace the origin and transformation of greenhouse gases, providing insights into their production pathways and helping quantify emissions more accurately;
–: Remote sensing and imaging spectral investigations by unmanned aerial vehicles and satellite-based remote sensing technologies offer broad coverage and continuous monitoring of large marine areas;
–: Integrated modeling approaches with machine learning techniques.

6. Conclusions

This study is a comprehensive analysis of a great number of published materials regarding hydrochemistry in Aniva Bay and the Sea of Okhotsk. Apart from the literature review, we analyzed the diagrams of the annual cycle of long-run mean hydrochemical values in Aniva Bay from averaged data of r/v observations with a time step of 5 days. It allowed us to not only confirm the main features of the hydrochemical regime described earlier [43], but to add significant details to this picture.

The annual cycle of dissolved oxygen (especially in surface layers) is negatively related to the water temperature and, basically, so are their mirror reflections. The main maximum is observed in spring (>9.3 mg/L or >114% at the 10 m depth) after ice melting in the bay, and the main minimum (<5.6 mg/L or <97% on the surface) in the beginning of September. A more detailed analysis in the upper 30 m layer reveals a local minimum in mid-July and a local maximum in the first half of August. In the central part of the bay, a maximum of oxygen is observed at the horizon of 20 m, probably caused by eddy dynamics.

The following features are observed in the annual pH cycle for all depths: maximums 8.25–8.26 in April–May, September, and November; minimum in August. The local pH maximum in the euphotic layer (8.22) is formed in July and is a consequence of additional inorganic nutrient inflow with precipitation and increased rate of CO₂ absorption during short-term active photosynthesis.

The annual cycle of phosphates in Aniva Bay is characterized by winter maximum (>19 µg/L on the surface, >47 µg/L on the 50 m depth) and summer minimum (<10 µg/L on the surface, <40 µg/L on the 50 m depth). A more detailed analysis demonstrates two local maximums: in the second half of June in the 0–50 m layer (comparable to winter maximum in absolute values) and in mid-August (comparable to winter maximum in absolute values on the 0–20 m depth and exceeding it by 20–30% in the 30–50 m depth).

There are two maximums of silica concentrations at the end of June and the beginning of September. Mean values are >340 µg/L on the surface and >800 µg/L on the 50 m depth. The main minimum is in mid-August with <70 µg/L on the surface and ~640 µg/L on the 50 m depth.

The summer maximum nitrites concentrations in the 0–50 m layer have average values of >2 µg/L in mid-June on the surface and >3.5 µg/L on the 50 m depth in mid-September. Its temporal scale matches the pH and silica minimums, as well as the maximum in the year’s cycle of phosphates and a local oxygen maximum.

Hence, this study presents the analysis of the data from instrumental observations of hydrochemical parameters for the previous climatic period. It is an important step towards the assessment of the carbon balance in Aniva Bay, chosen as the carbon site in Sakhalin Oblast. The performed data unification minimizes errors of further numerical calculations of baseline absorption and emission of carbon in Aniva Bay.

The obtained results can be used to study the ecological state of coastal waters, the carbon cycle and the GHG analysis in other water areas. As recommendations for politicians, businesspersons and scientists, we can offer the following basic approaches: complex approach including simultaneously measuring meteorological, hydrological and biohydrochemical parameters; using the autonomous stations for continuous monitoring; uniform measurement standards; and involvement of scientists and high-level specialists. These approaches will make it possible to implement both environmental and business projects related to the environmental sustainability of the coastal zone, the carbon cycle, mariculture, transport and recreation.

The obtained results systematize and significantly complement previously performed studies. Based on our results, expeditionary research in Aniva Bay can be planned, and an important practical significance is the use of our results for the research and validation of biogeochemical models.

Knowledge of seasonal and spatial variability of hydrochemical parameters in the water area of Aniva Bay, jointly with data on hydrological factors (water circulation, ice conditions, etc.), is a very important factor when choosing a location for the construction of a carbon landfill, carbon farms and other marine facilities, as well as for calculating various risks associated with carbon emissions.

[ad_2]