Abstract
The discharge of cooling water into surface waters containing residual heat energy is a pressing issue in the environmental impact of nuclear power. Cooling water flow can cause negative changes in the water and ecological processes of natural waters, in particular an imbalance of organic matter, in combination with a complex of natural factors. The purpose of the article was to perform a statistical analysis of the dynamics of organic matter content in the water of the Styr River in the Rivne Nuclear Power Plant (Rivne NPP) impact zone to assess the manifestation of the environmental impact of the discharge of cooling water with residual heat. The monitoring data series consisted of the results of monthly monitoring of the rate of chemical oxygen demand, the rate of biochemical oxygen demand in 5 days, and total organic carbon in the Styr River sections before and after the discharge of Rivne NPP cooling water in 2018–2022.
HIGHLIGHT
The purpose of the article was to perform a statistical analysis of the dynamics of organic matter content in the water of the Styr River in the Rivne Nuclear Power Plant impact zone to assess the manifestation of the environmental impact of the discharge of cooling water with residual heat.
INTRODUCTION
The thermal regime of a water body is a crucial factor for the quality of the ecosystem, and the thermal impact on water bodies caused by nuclear power plant (NPP) cooling water discharges can significantly affect the aquatic environment and its biota (Barescut et al. 2009). Because water treatment technologies can be used to minimize water consumption (Kuznietsov et al. 2024), recirculating cooling water systems (RCWS) are the most commonly used cooling systems in the power sector. Chemical and temperature effects occur when RCWS effluent is discharged into a body of water (Kuznietsov et al. 2023a). In the RCWS of NPPs, water-soluble components of the cooling water are heated, evaporated, and concentrated. National regulations and permits regulate the maximum allowable temperature increase, and chemical content of NPP cooling water discharges (Zhang et al. 2023). Water discharge into natural bodies of water can lead to changes in the chemical equilibrium of their constituents and pose a potential risk to human health (Zak et al. 2021).
The temperature effects of water discharges are determined by direct and indirect effects: direct effects include increased activity with accelerated digestion, increased food demand, reproductive disorders, and destruction of sensitive tissues of the nervous system of aquatic organisms; indirect effects cause negative changes in aquatic ecological processes, in particular, disturbance of the nutrient balance (Sandstrom 1997). There have been studies on the impact of the temperature effects of cooling water discharge on aquatic ecosystems, but the studies that have been conducted at real sites have only looked at the spatial and temporal interpretation of temperature (Jiang et al. 2018). In addition to thermal discharges, aquatic ecosystems face the impact of chemicals discharged from NPP with cooling water return flows, which should also be considered (Vries et al. 2008). Changing chemical balances in surface waters can be considered important indicators of aquatic ecological responses to NPP water discharges (Kuznietsov & Biedunkova 2023a).
Rivers are important sites for the transport of organic matter and are critical components of the global carbon (C) cycle. However, little is known about the longitudinal and temporal changes in organic matter in rivers (Soria-Reinoso et al. 2022). The implementation of measures to protect water resources from pollution and their rational use is an urgent task, which is also essential for the sustainable development of the entire energy sector (Kuznietsov & Biedunkova 2023b). Studies usually consider the temperature and chemical effects of cooling water discharge separately (Kuznietsov et al. 2023b). Our study links chemical to temperature effects, and the established correlations allow predicting the content of chemical indicators after cooling water discharge and focus on the potential impact and relationship of temperature with the content of organic matter in terms of chemical oxygen demand (COD), biochemical oxygen demand (BOD5), and total organic carbon (TOC).
BOD5 represents the oxygen demand needed microorganisms for the decomposition of organic matter, COD reflects the oxygen consumption during the chemical decomposition of organic matter, and TOC represents the concentration of carbon dioxide produced during the catalytic combustion of organic matter. Aguilar-Torrejón et al. (2023) show that the values of BOD5, COD, and TOC vary and depend on the type and parameters of the water analysed. COD and BOD5 are the most accurate measures of the increase in the amount of organic matter that can be oxidized by chemical or biological processes (Costa et al. 2018), and recommend using them together to assess organic pollution in water bodies, with COD values generally higher than BOD5 values and depending on the type of water analysed (Recoules et al. 2019). BOD5 is directly related to the amount of microbial contamination, and TOC is the most comprehensive analysis if the aim is to detect all types of organic matter present in water (Si et al. 2019). Despite the use of toxic chemicals for COD and very long analysis times (5 days) for BOD5, these methods are widely used and have not been completely replaced by other measurements, including TOC (Aguilar-Torrejón et al. 2023).
COD, TOC, and BOD5 can be considered complementary methods for monitoring organic matter content (Dubber & Gray 2010) Considering the relationship and different natures of COD and BOD5, the ratio of BOD5:COD is used to indicate biodegradability. If the ratio of COD:BOD5 is less than 0.1, it indicates the presence of organic matter that is difficult to biodegrade. The lower limit of the BOD5:COD ratio for organic matter to be biodegradable is 0.4 (Farraji et al. 2015). The BOD5:TOC ratio is also proposed as a measure of biodegradability (Ledakowicz 1998). Thus, the degree of biodegradation calculated from the BOD5:COD or BOD5:TOC ratio is commonly used for wastewater but is much more universal (Kowalski 1987). In our study, the BOD5:COD and BOD5:TOC ratios were used to assess the biodegradability of organic matter at different sites before and after the Rivne NPP discharge to determine the possible change in the nature of organic matter due to the possible impact of the cooling water discharge.
The purpose of our study is to investigate the spatial and temporal changes in the contents of BOD5, COD, and TOC and to establish correlation dependencies based on subject-specific approaches to the distribution of their content data in surface waters. The study was carried out on the example of the Styr River in the zone of influence of the water discharges of the Rivne NPP. The novelty of our study lies in the application of statistical regression analysis to the assessment and impact of cooling water discharges on the natural waters of the river. The practical value of this research lies in the possibility of applying this model to other RCWS power plants with cooling water discharges into water bodies.
METHODS
Characterization of methods for measuring the concentrations of BOD5, COD, and TOC used in the study (CI is the measurement range)
Indicator . | CI . | Relative measurement error δ (%) . | Method of measurement (standard in Ukraine) . |
---|---|---|---|
TOC (mgC/dm3) | 0.3–100* | From 0.3 to 10: δ = ±10; more than 10: δ = ±5 | DSTU ISO 1484:2003 |
BOD5, (mgO2/dm3) | 0.5–15 | From 0.5 to 2 inclusive: δ = ±(90–27); over 2 to 5 inclusive: δ = ±(27–11); over 5 to 15 inclusive: δ = ±(11–5) | KND 211.1.4.024 − 95 |
COD (mgO/dm3) | 5–100 | From 5 to 10 inclusive: δ = ±(65–34); over 10 to 30 inclusive: δ = ±(34–14); over 30 to 100: δ = ±(14–9) | KND 211.1.4.021 – 95 |
Indicator . | CI . | Relative measurement error δ (%) . | Method of measurement (standard in Ukraine) . |
---|---|---|---|
TOC (mgC/dm3) | 0.3–100* | From 0.3 to 10: δ = ±10; more than 10: δ = ±5 | DSTU ISO 1484:2003 |
BOD5, (mgO2/dm3) | 0.5–15 | From 0.5 to 2 inclusive: δ = ±(90–27); over 2 to 5 inclusive: δ = ±(27–11); over 5 to 15 inclusive: δ = ±(11–5) | KND 211.1.4.024 − 95 |
COD (mgO/dm3) | 5–100 | From 5 to 10 inclusive: δ = ±(65–34); over 10 to 30 inclusive: δ = ±(34–14); over 30 to 100: δ = ±(14–9) | KND 211.1.4.021 – 95 |
Note: According to DSTU EN 1484:2003, sample preparation for TOC determination includes pre-filtering the sample through a 0.45 ηm pore size filter.
Location of water sampling and monitoring sites of the water of the Styr River before (А−) and after (В−) water discharge from the Rivne NPP (Ukraine).
Location of water sampling and monitoring sites of the water of the Styr River before (А−) and after (В−) water discharge from the Rivne NPP (Ukraine).
Statistical processing of the study results included determination of the range of data series (min-max), arithmetic mean (M), standard deviation (±SD), coefficient of variation (CV), Pearson coefficient (r), significance of the connection (p), and coefficient of determination (R2) of the respective sample and factor analysis of data using the Minitab software package (Version 21.4.1, Minitab, LLC). Dataset analysis was performed using Pearson correlation analysis to evaluate the relationship between water quality variables (Barakat et al. 2016; Kuznietsov & Biedunkova 2024). Тhe principal component analysis (PCA) was applied efficiently to our data to identify the underlying interrelationship among the parameters (Hajigholizadeh & Melesse 2017). Partial least-squares regression (PLSR), as a multivariate regression method, was used to determine the linear relationship between a set of dependent response variables (X) and a set of predictor variables (Y) (Kahaer & Tashpolat 2019).
RESULTS AND DISCUSSION
Changes in the concentration of BOD5, COD, TOC, and temperature in the water of the Styr River in the area of influence of the Rivne NPP water discharges show a wide range of fluctuations of these indicators (Table 2). The statistical analysis of the laboratory control data shows that the average values of the water content of the Styr River differ slightly in the sections of the river before and after the Rivne NPP water discharges. The difference in the temperature of the Styr River water before and after the discharge of cooling water from the Rivne NPP did not exceed the established temperature increase effect under the conditions of the Permit of the State Agency for Water Resources of Ukraine (SAWRU) (Permit SAWRU 2020) of 3 °C and had an average value of 0.9 °C, with m = ± 0.54 °C.
Statistical descriptive the concentrations of BOD5, COD, and TOC and water temperature in different monitoring sites of the Styr River for 2018–2022
Indicator . | M . | ± SD . | min . | max . | CV . |
---|---|---|---|---|---|
The water intake of the Rivne NPP | |||||
BOD5 (mgO2/dm3) | 1.18 | 0.19 | 0.86 | 3.86 | 25.6 |
COD (mgO/dm3) | 42.6 | 20.3 | 17.6 | 83.2 | 43.9 |
TOC (mgC/dm3) | 11.6 | 4.2 | 5.23 | 20.0 | 56.5 |
Temperature (°C) | 11.3 | 10.6 | 0.8 | 24.3 | 92.5 |
After the water discharge of the Rivne NPP | |||||
BOD5 (mgO2/dm3) | 1.22 | 0.22 | 0.80 | 3.96 | 27.6 |
COD (mgO/dm3) | 40.6 | 22.3 | 15.1 | 82.4 | 43.8 |
TOC (mgC/dm3) | 12.2 | 3.9 | 6.22 | 25.4 | 42.3 |
Temperature (°C) | 12.2 | 11.2 | 1.7 | 25.6 | 94.5 |
Indicator . | M . | ± SD . | min . | max . | CV . |
---|---|---|---|---|---|
The water intake of the Rivne NPP | |||||
BOD5 (mgO2/dm3) | 1.18 | 0.19 | 0.86 | 3.86 | 25.6 |
COD (mgO/dm3) | 42.6 | 20.3 | 17.6 | 83.2 | 43.9 |
TOC (mgC/dm3) | 11.6 | 4.2 | 5.23 | 20.0 | 56.5 |
Temperature (°C) | 11.3 | 10.6 | 0.8 | 24.3 | 92.5 |
After the water discharge of the Rivne NPP | |||||
BOD5 (mgO2/dm3) | 1.22 | 0.22 | 0.80 | 3.96 | 27.6 |
COD (mgO/dm3) | 40.6 | 22.3 | 15.1 | 82.4 | 43.8 |
TOC (mgC/dm3) | 12.2 | 3.9 | 6.22 | 25.4 | 42.3 |
Temperature (°C) | 12.2 | 11.2 | 1.7 | 25.6 | 94.5 |
Dynamics of temporal change in the concentrations of BOD5, COD, TOC, and temperature in the water of the Styr River before (a) and after (b) water discharge from the Rivne NPP for 2018–2022.
Dynamics of temporal change in the concentrations of BOD5, COD, TOC, and temperature in the water of the Styr River before (a) and after (b) water discharge from the Rivne NPP for 2018–2022.
Seasonal variability of the temperature (a) and the concentrations of BOD5 (b), COD (c), and TOC (d) in the water of the Styr River for 2018–2022.
Seasonal variability of the temperature (a) and the concentrations of BOD5 (b), COD (c), and TOC (d) in the water of the Styr River for 2018–2022.
Relative variable importance components of BOD5, COD, ТОС, and temperature in the water of the Styr River before (а) and after (b) water discharge from the Rivne NPP for 2018–2022.
Relative variable importance components of BOD5, COD, ТОС, and temperature in the water of the Styr River before (а) and after (b) water discharge from the Rivne NPP for 2018–2022.
Dendrogram data of ВOD5, COD, and ТОС and temperature in the water of the Styr River before (А−) and after (В−) water discharge from the Rivne NPP for 2018–2022.
Dendrogram data of ВOD5, COD, and ТОС and temperature in the water of the Styr River before (А−) and after (В−) water discharge from the Rivne NPP for 2018–2022.
Dendrogram data of ВOD5, COD, ТОС, and temperature in the water of the Styr River before (А−) and after (В-) water discharge from the Rivne NPP for 2018–2022.
Dendrogram data of ВOD5, COD, ТОС, and temperature in the water of the Styr River before (А−) and after (В-) water discharge from the Rivne NPP for 2018–2022.
PCA рlot data of BOD5, COD, and TOC contents of the group in the water of the Styr River before (А−) and after (В-) water discharge from the Rivne NPP for 2018–2022.
PCA рlot data of BOD5, COD, and TOC contents of the group in the water of the Styr River before (А−) and after (В-) water discharge from the Rivne NPP for 2018–2022.
The normalized object scores and variable loadings on each PC were scaled proportionally to the root of the variance accounted for by that PC, as shown in Figure 7. The PC consists of moderate positive loadings for BOD5, COD, and TOC. The variables (Styr River before water intake and after discharge of the Rivne NPP) were found to be highly correlated and their contribution was significant in grouping the samples collected from the Styr River upstream and downstream sampling sites. It is important to note that although PCA can provide insight into the data on the distribution of water discharge components, it is not necessarily an optimal method for feature extraction, and other methods such as correlation analysis with subject-specific data partitioning approaches may be to complement PCA (Aguilar-Torrejón et al. 2023).
Contour plot the concentration of BOD5, COD, and ТОС in water of the Styr River before (а) and after (b) water discharge from the Rivne NPP for 2018–2022.
Contour plot the concentration of BOD5, COD, and ТОС in water of the Styr River before (а) and after (b) water discharge from the Rivne NPP for 2018–2022.
The relationship correlations of BOD5, COD, and TOC content in the water of the Styr River before (A−) and after (B−) water discharge from the Rivne NPP for 2018–2022: (a) COD, (b) BOD5, and (c) TOC.
The relationship correlations of BOD5, COD, and TOC content in the water of the Styr River before (A−) and after (B−) water discharge from the Rivne NPP for 2018–2022: (a) COD, (b) BOD5, and (c) TOC.
Summary statistics of correlations of BOD5, COD, and TOC concentrations and temperature (t, °C) in the water of the Styr River before (X) and after (Y)
Indicator (Y) . | Correlation between (X) . | R2 . | p . | Equation of correlation . |
---|---|---|---|---|
COD (0–80 mgО/dm3) | TOC (0–20 mgС/dm3) | 70.10 | 0.0003 | Y = 14.53 + 0.12*X + 2.438*t |
BOD5 (0–6 mgО/dm3) | 79.18 | 0.0004 | Y = 6.82 + 9.31*X + 1.243*t | |
TOC (0–20 mgС/dm3) | BOD5 (0–6 mgО2/dm3) | 72.27 | 0.0002 | Y = 6.82 + 4.31*X + 1.251*t |
COD (0–80 mgО/dm3) | 82.23 | 0.0002 | Y = 6.45 + 0.64*X + 0.184*t | |
BOD5 (0–6 mgО2/dm3) | COD (0–80 mgО/dm3) | 71.30 | 0.0005 | Y = 9.71–0.367*X + 0.152*t |
TOC (0–20 mgС/dm3) | 85.40 | 0.0009 | Y = −5.37 + 0.86*X + 0.135*t |
Indicator (Y) . | Correlation between (X) . | R2 . | p . | Equation of correlation . |
---|---|---|---|---|
COD (0–80 mgО/dm3) | TOC (0–20 mgС/dm3) | 70.10 | 0.0003 | Y = 14.53 + 0.12*X + 2.438*t |
BOD5 (0–6 mgО/dm3) | 79.18 | 0.0004 | Y = 6.82 + 9.31*X + 1.243*t | |
TOC (0–20 mgС/dm3) | BOD5 (0–6 mgО2/dm3) | 72.27 | 0.0002 | Y = 6.82 + 4.31*X + 1.251*t |
COD (0–80 mgО/dm3) | 82.23 | 0.0002 | Y = 6.45 + 0.64*X + 0.184*t | |
BOD5 (0–6 mgО2/dm3) | COD (0–80 mgО/dm3) | 71.30 | 0.0005 | Y = 9.71–0.367*X + 0.152*t |
TOC (0–20 mgС/dm3) | 85.40 | 0.0009 | Y = −5.37 + 0.86*X + 0.135*t |
Matrix plot for indicators (temperature (Т) in °C, COD in mgO/dm3, BOD5 in mgO2/dm3, and TOC in mgC/dm3) in water of the Styr River before (А−) and after (В-) water discharge from the Rivne NPP for 2018–2022.
Matrix plot for indicators (temperature (Т) in °C, COD in mgO/dm3, BOD5 in mgO2/dm3, and TOC in mgC/dm3) in water of the Styr River before (А−) and after (В-) water discharge from the Rivne NPP for 2018–2022.
Correlation coefficient r-Pearson for indicators (temperature (Т) in °C, COD in mgO/dm3, BOD5 in mgO2/dm3, and TOC in mgC/dm3) in water of the Styr River before (А−) and after (В-) water discharge from the Rivne NPP for 2018–2022.
Correlation coefficient r-Pearson for indicators (temperature (Т) in °C, COD in mgO/dm3, BOD5 in mgO2/dm3, and TOC in mgC/dm3) in water of the Styr River before (А−) and after (В-) water discharge from the Rivne NPP for 2018–2022.
Changes in the BOD5:COD, BOD5:TOC ratio (a) and correlation relationships between the ratio of BOD5:COD (b), BOD5:TOC (c) indicates in water the Styr River before (А− and after (В-) water discharge from the Rivne NPP for 2018–2022.
Changes in the BOD5:COD, BOD5:TOC ratio (a) and correlation relationships between the ratio of BOD5:COD (b), BOD5:TOC (c) indicates in water the Styr River before (А− and after (В-) water discharge from the Rivne NPP for 2018–2022.
In addition, the identified seasonal variations of COD, BOD5, and TOC (Figures 2(a) and 3) in the Styr River section up to the Rivne NPP outfall may be the initial data for identifying natural influencing factors. The data analysis methodology used in this study can separate the influencing anthropogenic factor is to assess the difference (increase) in COD, BOD5, TOC concentrations in the water of the Styr River before (A−) and after (B−) water discharge from the Rivne NPP (Figure 9 and Equations (1)–(3)). Moreover, taking into account the temperature increase in the Styr River water before (A−) and after (B−) due to the discharge of heated water from the Rivne NPP, the anthropogenic factor can be identified by water temperature (Table 3) when only one COD or BOD5 or TOC indicator and the corresponding temperature at the site before (A−) and after (B−) water discharge are known. Thus, it is proposed to use the geometric sum of the measurement error according to the applied measurement methodology as a criterion for the presence of anthropogenic influence on the formation of COD, BOD5, and TOC concentrations as a result of water discharge. Identification and separate influencing factors may be the subject of further studies. It is known (Kuznetsov & Tichomirov 2017) that considerable attention is being paid to the problems of NPP operation, so further research will be important in the future.
CONCLUSIONS
Based on the results of long-term measurements of BOD5, COD, and TOC concentrations in the Styr River water in the zone of influence of Rivne NPP wastewater, the dynamics of changes in the content of organic substances, the factors of variability, and correlations were determined. Currently, the results of the study are the initial data for further monitoring of possible abnormal changes in BOD5, COD, and TOC content and trends in the C cycle in the Styr River water, including changes due to the impact of anthropogenic factors of Rivne NPP effluent.
It should be noted that the seasonal variability of BOD5, COD, and TOC concentrations is characterized by minimum values in the cold season and maximum values in the warm season. The values of water content in the Styr River differ slightly in the sections of the watercourse before and after the Rivne NPP discharges. Statistically significant correlations were found between the indicators at different sites and between BOD5, COD, TOC, and temperature. The regression analysis revealed a statistically significant (p = 0.0005) direct relationship between the concentrations of COD and TOC in the water before (A− and after (B−) water discharge (r = 0.89, R2 = 81.45%). The calculated equations allow predicting the organic matter input with return cooling water based on the initial BOD5, COD, and TOC content in the Styr River water before water discharge from the Rivne NPP, taking into account temperature effects.
The calculated BOD5:TOC and BOD5:TОС ratios can provide useful information on the organic matter content of the water, but this is only one of many factors that need to be considered when assessing water quality and organic pollution potential, which will be the subject of further research. The obtained correlation for the BOD5:TOC ratio in the Styr River before and after water discharge from an NPP is very strong (r = 0.91, R2 = 95.12%) and statistically significant (p = 0.0001).
In general, the results of the study indicate that there is no negative impact of the Rivne NPP water discharge on the organic matter balance and therefore no aquatic ecological reactions to the thermal effects of the Rivne NPP cooling water discharge. Methodology of this study can be used to assess the aquatic content of carbon forms in water bodies and to evaluate the water discharges of any NPP with RCWS.
DATA AVAILABILITY STATEMENT
All relevant data are included in the paper or its Supplementary Information.
CONFLICT OF INTEREST
The authors declare there is no conflict.