Front. Public Health Frontiers in Public Health Front. Public Health 2296-2565 Frontiers Media S.A. 10.3389/fpubh.2024.1400680 Public Health Original Research Model prediction of radioactivity levels in the environment and food around the world’s first AP 1000 nuclear power unit Wang Peng 1 † Huang Wengzhe 2 † Zou Hua 1 Lou Xiaoming 1 Ren Hong 1 Yu Shunfei 1 Guo Jiadi 1 Zhou Lei 1 Lai Zhongjun 1 Zhang Dongxia 1 Xuan Zhiqiang 1 Cao Yiyao 1 * 1Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, Zhejiang, China 2Yangming College, Ningbo University, Ningbo, Zhejiang, China

Edited by: Antoaneta Ene, Dunarea de Jos University, Romania

Reviewed by: Hesham M. H. Zakaly, Ural Federal University, Russia

Hamdy Awad, Al-Azhar University, Assiut Branch, Egypt

 *Correspondence: Yiyao Cao, yycao@cdc.zj.cn

†These authors have contributed equally to this work

 15 05 2024 2024 12 1400680 14 03 2024 29 04 2024 Copyright © 2024 Wang, Huang, Zou, Lou, Ren, Yu, Guo, Zhou, Lai, Zhang, Xuan and Cao. 2024 Wang, Huang, Zou, Lou, Ren, Yu, Guo, Zhou, Lai, Zhang, Xuan and Cao

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

 Objectives

Model prediction of radioactivity levels around nuclear facilities is a useful tool for assessing human health risks and environmental impacts. We aim to develop a model for forecasting radioactivity levels in the environment and food around the world’s first AP 1000 nuclear power unit.

 Methods

In this work, we report a pilot study using time-series radioactivity monitoring data to establish Autoregressive Integrated Moving Average (ARIMA) models for predicting radioactivity levels. The models were screened by Bayesian Information Criterion (BIC), and the model accuracy was evaluated by mean absolute percentage error (MAPE).

 Results

The optimal models, ARIMA (0, 0, 0) × (0, 1, 1)4, and ARIMA (4, 0, 1) were used to predict activity concentrations of 90Sr in food and cumulative ambient dose (CAD), respectively. From the first quarter (Q1) to the fourth quarter (Q4) of 2023, the predicted values of 90Sr in food and CAD were 0.067–0.77 Bq/kg, and 0.055–0.133 mSv, respectively. The model prediction results were in good agreement with the observation values, with MAPEs of 21.4 and 22.4%, respectively. From Q1 to Q4 of 2024, the predicted values of 90Sr in food and CAD were 0.067–0.77 Bq/kg and 0.067–0.129 mSv, respectively, which were comparable to values reported elsewhere.

 Conclusion

The ARIMA models developed in this study showed good short-term predictability, and can be used for dynamic analysis and prediction of radioactivity levels in environment and food around Sanmen Nuclear Power Plant.

 radioactivity levels AP1000 time-series analysis ARIMA food environment 2024C03229 2022RC120 2022KY130 2022KY132 2023KY643 2024KY903 2024KY892 ZFY-2022-K-001 ZFY-2022-K-006 ZKF-2023-K-001 Zhejiang Provincial Key Research and Development Zhejiang Health Science and Technology Plan South Zhejiang Institute of Radiation Medicine and Nuclear Technology section-at-acceptance Radiation and Health Introduction

As a clean energy with near zero emissions, nuclear power has been vigorously developed in China in recent years as one of the important measures to achieve carbon neutrality by 2060 (1). Currently, there are 55 nuclear reactors in operation, and 18 under construction in China (2). During the operation of nuclear power plants, radioactive debris or effluents are inevitably discharged into the environment through air and water. As the continuous expansion of nuclear energy, the levels of radioactivity in environment and food have become a major concern for residents around the nuclear power plant (3–5), especially after the Fukushima Daiichi Nuclear Power Plan accident in 2011.

Radioactive substances present in the environment could potentially induce radiation exposure to human via external radiation and internal radiation through absorption, inhalation, and ingestion. As per the findings of the United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR 2000), an estimated 8% of natural human radiation exposure can be attributed to the ingestion of water and food (6). In order to safeguard the well-being of inhabitants, the World Health Organization (WHO) has established thresholds for gross alpha (0.5 Bq/L) and gross beta (1.0 Bq/L) in drinking water as a means of ensuring radiation safety (7). 90Sr is a high-yield byproduct of nuclear fission (8, 9), possessing a half-life of 28.8 years. Its primary route of entry into human body is through the food chain, where it can accumulate in teeth, bones, and muscle tissues. Given its relatively long half-life and high radiotoxicity, 90Sr is recognized as a crucial artificial radionuclide for evaluating radiation risks to both the environment and human health (10, 11).

In view of the radiation exposure risks of and high level of public concern, prediction of radioactivity levels in the environment and food around the nuclear power plant are crucial to ensure radiation safety for the public and the environment.

Time series analysis is widely recognized as a valuable predictive model, and the Autoregressive Integrated Moving Average model (ARIMA) being one of the most important models. The ARIMA model, a hybrid of autoregressive and moving average models, was initially proposed in 1976 by Box and Jenkins (12). So far, it has been widely used in economic management, meteorological prediction, environmental prediction, and disease prediction (13–16). In addition, ARIMA model also has been applied in the field of radioactive monitoring and evaluation. In four districts of Istanbul, the model was used to predict concentrations of 226Ra, 232Th, and 40K (17). After Fukushima Daiichi Nuclear Power Plan accident, Hemn Salh et al. reported the use of ARIMA models for prediction of air radiation dose rates around (18). Some scholars have already used it to predict radon concentrations and thus to predict earthquakes (19, 20).

Sanmen Nuclear Power Plant (SNPP), located in Sanmen, Zhejiang, China, adopts the world’s most advanced third-generation pressurized water reactor (AP1000) technology, which is one of the achievements of China’s efforts to develop nuclear power. With the operation of SNPP, its impacts on the environment and residents’ health have become a growing concern. To assess the impacts, we have continuously monitored radioactivity levels in the environment and food around SNPP since 2011. In this study, we analyzed the historical monitoring data and used ARIMA models, which is a time-series analysis technology, to fit and predict radioactive levels around SNPP for the first time. The predicted data can provide a basis for environmental impact assessment and human health risk assessment around SNPP.

 Data and methods The study region

All monitoring stations were located in Sanmen County. There were four monitoring stations for gross α and gross β in water representing surface water, factory water, tap water, and well water, respectively. Water samples were collected quarterly. There were three monitoring stations for 90Sr in collected quarterly food including mullet, crucian carp, cabbage, and rice. There were 30 monitoring stations for ambient radiation, which were monitored quarterly. The images of the collected samples are shown in Figure 1. All data were obtained from Zhejiang Provincial Center for Disease Control and Prevention and Taizhou City Center for Disease Control and Prevention. In this study, all these monitoring data were averaged across the corresponding monitoring stations as shown in Figures 2, 3. The dataset of gross α in water was excluded because the activity concentrations were mostly below the detection limits.

Samples collected in this study.

Monitoring stations of water and food around SNPP. The map was produced by software of Lantu Draw (URL link: https://www.ldmap.net/).

Monitoring stations for ambient radiation exposure of ambient environmental around SNPP. The map was produced by software of Lantu Draw (URL link: https://www.ldmap.net/).

 Sample preparation and analysis Gross <italic>α</italic> and <italic>β</italic> in water

Figure 4 illustrates the schematic flow for determination of gross α and β in water, according to the Chinese national standard (21).

Procedure for determination of gross α and β in water.

 <sup>90</sup>Sr in food

Figure 5 illustrates the schematic flow for determination of 90Sr in food, according to the Chinese national standard (22).

Procedure for determination of 90Sr determination in food.

 Ambient radiation

The cumulative ambient dose (CAD) was monitored by thermoluminescent dosimeter (TLD) (23). The LiF (Mg, Cu, and P) powder was placed into a hermetically sealed container in order to fabricate TLD. Each monitoring point was equipped with two TLDs installed at a height of 2 m. The TLDs were measured by TLD reader.

 Data

The model training dataset (Supplementary Tables S1–S3) included 90Sr activity concentrations in food from the second quarter (Q2) of 2011 to the fourth quarter (Q4) of 2022, gross β activity concentrations in water from the first quarter (Q1) of 2016 to Q4 of 2022 and CAD from Q1 of 2011 to Q4 of 2022. The monitoring data from Q1 to Q4 in 2023 were used as a test dataset to assess the predictability of the models using the mean absolute percentage error (MAPE) between the forecasted and observed values. Finally, the best model was applied to forecast the radioactivity levels from Q1 to Q4 of 2024.

 Methods Auto-regressive integrated moving average [ARIMA (<italic>p</italic>,<italic>d</italic>,<italic>q</italic>)]

An ARIMA model is defined by three parameters: p, d, and q, where p is the order of Auto-Regressive (AR) term, d is the order of differencing required to make the time-series stationary, and q is the order of Moving Average (MA) term.

AR (p) usually explains the present value Xt, unidirectionally it terms of its previous values Xt − 1,Xt − 2,, Xt − p, and the current residuals εt. It can be expressed as Equation (1):

X t = ϕ 1 X t − 1 + ϕ 2 X t − 2 + ⋯ + ϕ p X t − p + ε t

The model illustrates a linear association between the current observed value of the sequence at time t and the past observed values at the preceding p time points. This type of regression was known as autoregression because it is based on its own historical data. The regression with the observed values from the previous p time points is referred to as p-order autoregression. φi (i = 1, 2…p) is its partial regression coefficient.

MA (q) refers to the current value of the time series Xt in terms of its current and previous residuals εt − 1, εt − 2,…, εt − q. It can be expressed as Equation (2):

X t = ε t − θ 1 ε t − 1 − θ 2 ε t − 2 − ⋯ − θ q ε t − q

The model indicates that the value of the sequence at time t is independent of the past observed values at the preceding q time points, however, it exhibits a linear relationship with the preceding q stochastic disturbances. Therefore, the model is referred to as the q-order moving average model. The εt is indicative of a stochastic disturbance sequence, also referred to as a white noise sequence, that is characterized by independence and adherence to a normal distribution. θi (i = 1, 2…q) is its partial regression coefficient.

The ARIMA model is the combination of AR model and MA model algorithms. I in the ARIMA (p,d,q) refers to Integrated. When time-series is stationary, the ARIMA (p,d,q) model is ARMA (p,q), can be expressed as Equation (3):

X t = ϕ 1 X t − 1 + ϕ 2 X t − 2 + ⋯ + ϕ p X t − p + ε t − θ 1 ε t − 1 − θ 2 ε t − 2 − ⋯ − θ q ε t − q

Transform the Equation (3) to Equation (4).

X t − ϕ 1 X t − 1 − ϕ 2 X t − 2 − ⋯ − ϕ p X t − p = ε t − θ 1 ε t − 1 − θ 2 ε t − 2 − ⋯ − θ q ε t − q

In order to facilitate ease of expression, a backshift operator is employed, which is akin to a time pointer. The multiplication of the present sequence value by a backshift operator is tantamount to shifting the temporal position of said sequence value one moment into the past. If denoted as B, the backshift operator yields Equation (5) and (6):

B m X t = X t − m B m ε t = ε t − m

Upon utilization of the backshift operator, Equation (4) is transformed to Equation (7).

X t = θ B ϕ B ε t

Where φ(B) is non-seasonal autoregressive polynomial, θ(B) is non-seasonal movingaverage polynomial, which are expressed as Equation (8) and (9):

φ B = 1 − φ 1 B − φ 2 B 2 − .. − φ p B p θ B = 1 − θ 1 B − θ 2 B 2 − .. − θ q B q

If the time-series is non-stationary, which should be converted to stationary by differencing, the model used is ARIMA (p,d,q). Assume ΔXt is the sequence which obtained after first-order differencing of Xt, which is expressed as Equation (10):

Δ X t = X t − X t − 1

Therefore, the equation of ARIMA (p,d,q) as follows Equation (11), which transforms from Equation (7).

Δ d X t = θ B ϕ B ε t Multiple seasonal auto-regressive integrated moving average [MSARIMA <italic>(p, d, q) (P, D, Q)</italic>]

The Seasonal-ARIMA [SARIMA(p, d, q)(P, D, Q)] model incorporates a non-seasonal ARIMA(p, d, q) model along with supplementary seasonal terms (P, D, Q)s, which accounts for the seasonality inherent in the time-series data over S time steps, corresponding to a singular seasonal period. P is the order of seasonal autoregressive term, D is the order of seasonal differencing, Q is the order of seasonal moving average term and S is the length of the seasonal cycle. The complete expression of the SARIMA model can be written as Equation (12):

Δ d Δ S D X t = θ B V B s ϕ B U B s ε t

Where εt was white noise; U(BS) is seasonal autoregressive polynomial, and V(BS) is seasonal moving average polynomial, which are expressed as Equation (13) and (14):

U B S = 1 − u 1 B S − u 2 B 2 S − .. − u P B P S V B S = 1 − v 1 B S − v 2 B 2 S − .. − v Q B Q S

When P = Q = D = 0, it means that there is no seasonal in the model. In this case, the model reduces to the standard ARIMA model.

 Results and discussion Model establishment

Figures 4, 5 show the time-series diagrams and autocorrelation function (ACF) plots of the original sequences of the CAD, 90Sr activity concentrations in food, and gross β activity concentrations in water. The time-series diagram for the CAD (Figure 6A) shows fluctuations within a certain range and most ACFs (Figure 7A) fall into the confidence interval and tend to zero rapidly. This indicates that the sequence is stable and does not need differencing. In contrast, both the time-series diagrams (Figure 6B) and ACF plots (Figure 7B) for the 90Sr activity concentrations in food exhibit clear periodic changes, indicating that the sequences require seasonal differencing. For gross β activity concentrations in water, the time series diagram (Figure 6C) appears relatively stable without any obvious trend, but all ACFs (Figure 7C) fall within the confidence interval. Therefore, it was determined to be white noise, and further analysis was halted.

Time-series diagram of monitoring data. (A) CAD; (B) 90Sr activity concentrations in food; (C) Gross β activity concentrations in water.

ACF plot of original sequence. (A) CAD; (B) 90Sr activity concentrations in food; (C) Gross β activity concentrations in water.

After performing first-order seasonal differencing on the original time-series of 90Sr activity concentrations in food, the ACF plot of the resulting time-series is shown in Figure 8. The majority of ACFs fall within the confidence interval and quickly reach zero, indicating that the new sequence is stable and does not exhibit any significant seasonal fluctuations.

ACF plot of the sequence of 90Sr activity concentrations in food after performing first-order seasonal differencing.

The determination of p, q, P and Q is a crucial aspect for establishing the ARIMA model. Some researchers suggested that the most effective method to determine these values involves analyzing ACF and partial autocorrelation function (PACF) of time-series after differencing and seasonal differencing. The lag number of the peak entering the confidence level in ACF and PACF plots is used to determine the values of p, q, P and Q (24–28). However, some scientists criticized this estimation method for not being sufficient every time (29). Besides the ACF and PACF plots, there are other methods that can be used to determine the optimal ARIMA model, such as the Bayesian Information Criterion (BIC). The model with the lowest BIC value is considered as the best fit.

In this study, we tested various parameter values for p, q, P and Q (with a maximum value of 4) in ascending order, and selected the models with the lowest BIC values. After model identification and parameter estimation, we selected the models that met all standards, as shown in Table 1.

Optimal models.

Name Optimal model BIC Ljung-Box Q-test
p values
90Sr in food ARIMA(0, 0, 0) × (0, 1, 1)4 −3.024 0.388
CAD ARIMA(4, 0, 1) −5.532 0.927

Another key point of model construction is the diagnosis of residual sequence and model parameters. If the residual sequence of the model follows a normal distribution and appears random like white noise, it suggests that the model has extracted most of the information from the original sequence (30). The residual sequences that are considered as white noise could be verified by Ljung-Box Q-test (31, 32) with p value greater than 0.05. In addition, the parameters of the model should also be significantly different from zero (p < 0.05) to ensure the validity of the model. The models shown in Table 1 all passed both residual sequences Ljung-box Q-test and model parameters diagnosis, which were eventually selected as the optimal models.

 Model fitting and prediction

Figure 9 shows the model forecasting results covering the training, testing and prediction periods. Each diagram consists of the observed, fitted and forecasted data, as well as the upper and lower confidence limits. The fitted values in the diagram follow the same trend as the observed values and are within the confidence interval, indicating that the fitting effect of the model is good. To further the accuracy of the model, MAPEs were calculated between the observed and the forecast values for 90Sr activity concentrations in food, CAD from Q1 to Q4 in 2023, which were 21.4 and 22.4%, respectively (Table 2). These low error values indicate that the model has a high level of accuracy in predicting the future values of these variables.

Time-series diagram for the model fit and forecast. (A) 90Sr activity concentrations in food; (B) CAD.

Observed and predicted radioactivity levels around SNPP from Q1 to Q4 in 2023.

Quarter Observed value Forecast value 95% confidence interval Absolute percentage error
90Sr activity concentration in food (Bg/kg) CAD (mSv) 90Sr activity concentration in food (Bg/kg) CAD (mSv) 90Sr activity concentration in food (Bg/kg) CDA (mSv) 90Sr activity concentration in food (%) CAD (%)
1 0.40 0.183 0.43 0.133 0.150–0.984 0.044–0.222 7.5 27.3
2 0.92 0.111 0.77 0.130 0.270–1.764 0.036–0.223 16.3 17.1
3 0.18 0.108 0.22 0.109 0.076–0.497 0.015–0.204 22.2 0.9
4 0.048 0.099 0.067 0.055 0.024–0.154 −0.040-0.150 39.6 44.4
Mean absolute percentage error (MAPE) 21.4 22.4

In this study, the established optimal models were applied to forecast the radioactivity levels around SNPP in 2024. As shown in Table 3, the activity concentrations of 90Sr in food from Q1 to Q4 were predicted to be 0.067–0.77 Bq/kg, and the CAD to be 0.067–0.129 mSv. The forecasted values of 90Sr in food are lower than the concentration limit recommended by the Chinese standard (33) and comparable to the activity concentrations of 90Sr in food from other regions of the world (Table 4) (34–40). The forecasted values of CAD are comparable to the level of Qinshan Nuclear Power Plant (0.073–0.093 mSv) (41).

Forecast values of radioactivity levels around SNPP from Q1 to Q4 in 2024.

Quarter Forecast value
90Sr activity concentration in food (Bq/kg) CAD (mSv)
1 0.43 0.129
2 0.77 0.123
3 0.22 0.112
4 0.067 0.067

Comparison of activity concentrations of 90Sr in foods from different regions of the world.

Location Sample type(sample time) Activity concentration (Bq/kg)
Vicinity of Ningde Nuclear Power Plant, China (34) Rice, Vegetables, Marine fish, Freshwater fish (2013–2017) <0.017–0.677
Qinshan Nuclera Power Plant, China (35) Rice, Salsola, Mullet, Crucian carp (2012–2019) 0.04–1.3
Tianwan Nuclear Power Plant, China (36) Rice, Chinese cabbage, wheat (March 2000 to April 2002) 0.024–0.23
Cuban (37) Vegetables, Fish (/) 0.0034–0.28
Austrian (38) Cereals, Cabbage, Freshwater fish (1997) 0.09–0.12
Mayak Industrial Association (39) Wheat, Cabbage, Carrots (2008–2010) 0.029–0.24
Niigata, Japan (40) Undaria pinnatifida (1999–2007) <0.016–0.036

Based on radioactivity monitoring data around SNPP, ARIMA models were established to fit and predict the short-term changes. However, there were instance where the forecast values and actual values of individual quarters differ greatly, such as 90Sr in food in Q1 and Q2 of 2019, CAD in Q4 of 2017 and 2020, and Q2 of 2021 (Figure 7). These differences may be related to large fluctuations of external conditions, such as climate (42), which were accounted by the ARIMA model. While ARIMA models focus on the role of time factors in fitting and forecasting, they do not analyze and discuss the relationship between the prediction object and the influencing factors (43).

When the prediction time becomes longer, these external influencing factors bring greater changes. Therefore, despite ARIMA models can maintain high accuracy when using historical data for short- and medium-term prediction, their accuracy may decrease when predicting further into the future. To improve the models accuracy in predicting radioactivity levels around SNPP for a longer period, continuous radioactivity monitoring data collection should be carried out to adjust and refine the model over time (44). This will ensure that the models reflect the variability and trends of radioactivity levels, leading to the best possible predictions.

In this study, the monitoring data from 2011 to 2022 as the training dataset. In June 2018, after the operation of SNPP, significant changes have taken place in external conditions. In our previous research, we found that the contributions of radioactive substances released to the environment after SNPP operation was negligible (45, 46), we think this change does not affect the establishment of the model. In addition, if a nuclear accident or other unknown external input occurs, the model cannot predict. When the observed data are significantly higher than predicted, excluding external input, it indicates that SNPP significantly releases radioactive substances.

The ARIMA model has been successfully applied to SNPP, but its applicability to other nuclear power plants with different technologies, operating practices or environmental backgrounds, such as Qinshan Nuclear Power Plant, will depend on the regularity and randomness of the time series of their historical radioactivity data. If the time series has a certain regularity and norandom, the ARIMA model has good potential to be applied.

 Conclusion

This study attempts to develop ARIMA models for medium and short-term prediction of radioactive levels in the environment and food around SNPP, utilizing historical time series data spanning from 2011 to 2023. The data from 2011 to 2022 were used as the training set, while the data from 2023 served as the testing set. The time series of gross β in water was a white noise series, which has no value in establishing a model. In contrast, the time series of 90Sr in food and CAD were stationary, non random, and have short-term correlation, making them suitable for establishing ARIMA models. In this study, our established models showed good consistency between the fitted values and observed values, coupled with the relatively small MAPEs, suggesting satisfactory fitting effect and accuracy in the ARIMA models.

The CAD and 90Sr activity concentrations in food in 2024 were predicted using the established model. The forecasted values of 90Sr in food are lower than the recommended threshold by the Chinese standard, and comparable to the active concentrations of 90Sr in food from other regions of the world. The forecasted values of CAD are comparable to the level of Qinshan Nuclear Power Plant.

However, we acknowledge the limitation of the ARIMA models, such as some external factors, e.g., climate, were not taken into account. This limitation might be compensated by continuous feeding of monitoring data thus to adjust and improve the models performance over time. In short, ARIMA models can be used as an additional tool for environmental radioactivity monitoring and human health risk assessment around SNPP.

 Data availability statement

The original contributions presented in the study are included in the article/Supplementary material; further inquiries can be directed to the corresponding author.

 Author contributions

PW: Conceptualization, Methodology, Project administration, Software, Writing – original draft, Writing – review & editing. WH: Data curation, Investigation, Software, Writing – original draft. HZ: Supervision, Writing – review & editing. XL: Supervision, Writing – review & editing. HR: Data curation, Investigation, Methodology, Writing – original draft. SY: Data curation, Investigation, Writing – original draft. JG: Data curation, Investigation, Writing – original draft. LZ: Data curation, Investigation, Writing – original draft. ZL: Data curation, Writing – original draft. DZ: Data curation, Writing – original draft. ZX: Writing – original draft, Formal analysis. YC: Supervision, Writing – review & editing.

 Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by Zhejiang Provincial Key Research and Development Plan (No. 2024C03229), Zhejiang Health Science and Technology Plan (Nos. 2022RC120, 2022KY130, 2022KY132, 2023KY643, 2024KY903, and 2024KY892), and Project of South Zhejiang Institute of Radiation Medicine and Nuclear Technology (Nos. ZFY-2022-K-001, ZFY-2022-K-006, and ZKF-2023-K-001).

 Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2024.1400680/full#supplementary-material

 References The State Council of the People’s Republic of China (2021). The government work report of China in 2021. Available at: http://www.gov.cn/zhuanti/2021lhzfgzbg/index.htm (Accessed March 13, 2024). International Atomic Energy Agency (2024). The database on nuclear power reactors of China. Available at: https://pris.iaea.org/PRIS/CountryStatistics/CountryDetails.aspx?current=CN (Accessed March 13, 2024). Pidgeon NF Lorenzoni I Poortinga W. Climate change or nuclear power—no thanks! A quantitative study of public perceptions and risk framing in Britain. Glob Environ Chang. (2008) 18:69–85. doi: 10.1016/j.gloenvcha.2007.09.005 Cardis E Hatch M. The Chernobyl accident--an epidemiological perspective. Clin Oncol (R Coll Radiol). (2011) 23:251–60. doi: 10.1016/j.clon.2011.01.510, PMID: 21396807 Tuo F Xu C Zhang J Zhou Q Li W Zhao L . Radioactivity analysis following the Fukushima Dai-ichi nuclear accident. Appl Radiat Isot. (2013) 78:77–81. doi: 10.1016/j.apradiso.2013.04.002 UNSCEAR. Sources and effects of ionizing radiation. United Nations Scientific Committee on the Effects of Atomic Radiation, New York. Available at: https://www.unscear.org/unscear/en/publications/2000_1.html WHO. Guidelines for Drinking-Water Quality. 4th ed. Geneva: World Health Organization (2011). Miki S Fujimoto K Shigenobu Y Ambe D Kaeriyama H Takagi K . Concentrations of 90Sr and 137Cs/90Sr activity ratios in marine fishes after the Fukushima Dai-ichi nuclear power plant accident. Fish Oceanogr. (2017) 26:221–33. doi: 10.1111/fog.12182 Tazoe H Yamagata T Tsujita K Nagai H Obata H Tsumune D . Observation of dispersion in the Japanese coastal area of released 90Sr, 134Cs, and 137Cs from the Fukushima Daiichi nuclear power plant to the sea in 2013. Int J Environ Res Public Health. (2019) 16:4094. doi: 10.3390/ijerph16214094, PMID: 31652941 Merz S Shozugawa K Steinhauser G. Analysis of Japanese radionuclide monitoring data of food before and after the Fukushima nuclear accident. Environ Sci Technol. (2015) 49:2875–85. doi: 10.1021/es5057648, PMID: 25621976 Wang P Wang R Zou H Lou XM Ren H Zhou L . Investigation of environmental radioactivity levels and human health risk around San'ao nuclear power plant (China) before operation. Environ Technol Innov. (2023) 32:103451. doi: 10.1016/j.eti.2023.103451 Box GEP Jenkins GM. Time Series Analysis: Forecasting and Control. San-Franisco: Holden-Day (1976). Sandhir V Kumar V Kumar V. Prognosticating the spread of Covid-19 pandemic based on optimal Arima estimators. Endocr Metab Immune Disord Drug Targets. (2021) 21:586–91. doi: 10.2174/1871530320666201029143122, PMID: 33121426 Yu L Wu CX Xiong NN. An intelligent data analysis system combining ARIMA and LSTM for persistent organic pollutants concentration prediction. Electronics. (2022) 11:652. doi: 10.3390/ELECTRONICS11040652 Selvaraj JJ Arunachalam V Coronado-Franco KV Romero-Orjuela LV Ramírez-Yara YN. Time-series modeling of fishery landings in the Colombian Pacific Ocean using an ARIMA model. Reg Stud Mar Sci. (2020) 39:101477. doi: 10.1016/j.rsma.2020.101477 Dimri T Ahmad S Sharif M. Time series analysis of climate variables using seasonal ARIMA approach. J Earth Syst Sci. (2020) 129:149. doi: 10.1007/s12040-020-01408-x Külahcı F Aközcan S Günay O. Monte Carlo simulations and forecasting of Radium-226, Thorium-232, and Potassium-40 radioactivity concentrations. J Radioanal Nucl Chem. (2020) 324:55–70. doi: 10.1007/s10967-020-07059-y Salh H Külahcı F Aközcan S. A mobile simulation and ARIMA modeling for prediction of air radiation dose rates. J Radioanal Nucl Chem. (2021) 328:889–901. doi: 10.1007/s10967-021-07726-8 Ambrosino F Thinová L Briestenský M Šebela S Sabbarese C. Detecting time series anomalies using hybrid methods applied to radon signals recorded in caves for possible correlation with earthquakes. Acta Geod Geophys. (2020) 55:405–20. doi: 10.1007/s40328-020-00298-1 Stránský V Thinová L. Radon concentration time series modeling and application discussion. Radiat Prot Dosim. (2017) 177:155–9. doi: 10.1093/rpd/ncx207, PMID: 29045758 Ministry of Health of the People’s Republic of China. Standard Examination Methods for Drinking Water-Radiological Parameters: GB/T 5750.13–2006. Beijing: Standards Press of China (2006). Ministry of Environmental Protection of the People's Republic of China. Radiochemical Analysis of Strontium-90 in Water and Ash of Biological Samples: HJ 815–2016. Beijing: Standards Press of China (2016). General Administration of Quality Supervision. Thermoluminescence Dosimetry Systems for Personal and Environmental Monitoring: GB/T10264-2014. Beijin: Standards Press of China (2014). Taneja K Ahmad S Ahmad K Attri SD. Time series analysis of aerosol optical depth over New Delhi using Box–Jenkins ARIMA modeling approach. Atmos Pollut Res. (2016) 7:585–96. doi: 10.1016/j.apr.2016.02.004 Imai C Armstrong B Chalabi Z Mangtani P Hashizume M. Time series regression model for infectious disease and weather. Environ Res. (2015) 142:319–27. doi: 10.1016/j.envres.2015.06.040 Lin YL Chen M Chen GW Wu XQ Lin TQ. Application of an autoregressive integrated moving average model for predicting injury mortality in Xiamen, China. BMJ Open. (2015) 5:e8491. doi: 10.1136/bmjopen-2015-008491, PMID: 26656013 Alsharif MH Younes MK Kim J. Time series ARIMA model for prediction of daily and monthly average global solar radiation: the case study of Seoul, South Korea. Symmetry. (2019) 11:240. doi: 10.3390/sym11020240 Lugendo P Tadesse T Wortmann C Neale C. Climate variability implications for maize yield food security and rural poverty in Tanzania. Annu Agric Policy Conf. (2018). Gairaa K Khellaf A Messlem Y Chellali F. Estimation of the daily global solar radiation based on Box–Jenkins and ANN models: a combined approach. Renew Sust Energ Rev. (2016) 57:238–49. doi: 10.1016/j.rser.2015.12.111 Katimon A Shahid S Mohsenipour M. Modeling water quality and hydrological variables using ARIMA: a case study of Johor River, Malaysia. Sustain Water Resour Manag. (2018) 4:991–8. doi: 10.1007/s40899-017-0202-8 Bouznad I Guastaldi E Zirulia A Brancale M Barbagli A Bengusmia D. Trend analysis and spatiotemporal prediction of precipitation, temperature, and evapotranspiration values using the ARIMA models: case of the Algerian highlands. Arab J Geosci. (2021) 13:1281. doi: 10.1007/s12517-020-06330-6 Hassan J. ARIMA and regression models for prediction of daily and monthly clearness index. Renew Energy. (2014) 68:421–7. doi: 10.1016/j.renene.2014.02.016 Ministry of Health of the People's Republic of China. Limited Concentrations of Radioactive Materials in Foods: GB14882-94. Beijing: Standards Press of China (1994). Zheng QS Yan Z Hua HL Xing ZS Min ZS Jia L . Investigation on 90Sr concentrations in foods near Ningde nuclear power plant in 2013-2017. Chin J Radiol Med Prot. (2019) 39:931–5. doi: 10.3760/cma.j.issn.0254-5098.2019.12.010 Cao YY Zhao ZX Wang P Yu SF Lai ZJ Zhang MB . Long-term variation of 90Sr and 137Cs in environmental and food samples around Qinshan nuclear power plant, China. Sci Rep. (2021) 11:20903. doi: 10.1038/s41598-021-00114-y, PMID: 34686687 Lu JG Huang YJ Li F Wang LH Li SZ Hsia YF. The investigation of 137Cs and 90Sr background radiation levels in soil and plant around Tianwan NPP, China. J Environ Radioact. (2006) 90:89–99. doi: 10.1016/j.jenvrad.2006.06.002, PMID: 16887243 Tomás Zerquera J Prendes Alonso M Fernández Gómez IM Rodríguez Castro GV Martínez Ricardo N López Bejerano G . Studies on internal exposure doses received by the Cuban population due to the intake of radionuclides from the environmental sources. Radiat Prot Dosim. (2006) 121:168–74. doi: 10.1093/rpd/ncl002, PMID: 16513820 Mück K Sinojmeri M Whilidal H Steger F. The long-term decrease of 90Sr in the environment and its transfer to man after a nuclear fallout. Radiat Prot Dosim. (2001) 94:251–9. doi: 10.1093/oxfordjournals.rpd.a006497 Kazachonok NN Kostyuchenko VA Popova IY Polyanchikova GV. Research into the factors affecting the uptake of 90Sr and 137Cs by plants from radioactively contaminated soils of the Southern Urals. Health Phys. (2012) 103:56–7. doi: 10.1097/HP.0b013e318249ba1a, PMID: 22647915 Morita T Fujimoto K Kasai H Yamada H Nishiuchi K. Temporal variations of 90Sr and 137Cs concentrations and the 137Cs/90Sr activity ratio in marine brown algae, Undaria pinnatifida and Laminaria longissima, collected in coastal areas of Japan. J Environ Monit. (2010) 12:1179–86. doi: 10.1039/B920173D, PMID: 21491686 Cao YY Lin JP Zhai KL Jiang W Zou H Ren H . Long-term investigation of environmental radioactivity levels and public health around the Qinshan nuclear power plant, China. Sci Rep. (2022) 12:4945. doi: 10.1038/s41598-022-09091-2, PMID: 35322138 Dahlgaard H Eriksson M Nielsen SP Joensen HP. Levels and trends of radioactive contaminants in the Greenland environment. Sci Total Environ. (2004) 331:53–67. doi: 10.1016/j.scitotenv.2004.03.023, PMID: 15325141 Gao SC Wu HM. Prediction of infection incidence of carbapenem-resistant pseudomonas aeruginosa by ARIMA time series model. Chin J Hosp Stat. (2021) 28:401–4. doi: 10.3969/j.issn.1006-5253.2021.05.004 Contreras J Espinola R Nogales JF Conejo JA. ARIMA models to predict next-day electricity prices. IEEE Trans Power Syst. (2003) 18:1014–20. doi: 10.1109/TPWRS.2002.804943 Ren H Yu SF Wang ZY Zheng TT Zou H Lou XM . Assessment of radiation exposure and public health before and after the operation of Sanmen nuclear power plant. Front Public Health. (2023) 11:1131739. doi: 10.3389/fpubh.2023.1131739, PMID: 36815151 Wang P Yu SF Zou H Lou XM Ren H Zhou L . Levels, sources, variations, and human health risk assessment of 90Sr and 137Cs in water and food around Sanmen nuclear power plant (China) from 2011 to 2020. Front Public Health. (2023) 11:1136623. doi: 10.3389/fpubh.2023.1136623, PMID: 36908432