Estimation of Radioactivity Release Activity Using Non-Linear Kalman Filter-Based Estimation Techniques

: The estimation of radioactivity release following an accident in a nuclear power plant is crucial due to its short and long-term impacts on the surrounding population and the environment. In the case of any accidental release, the activity needs to be estimated quickly and reliably to effectively plan a rapid emergency response and design an appropriate evacuation strategy. The accurate prediction of incurred dose rate during normal or accident scenario is another important aspect. In this article, three different non-linear estimation techniques, extended Kalman ﬁlter, unscented Kalman ﬁlter, and cubature Kalman ﬁlter are proposed in order to estimate release activity and to improve the prediction of dose rates. Radionuclide release rate, average wind speed, and height of release are estimated using the dose rate monitors data collected in proximity of the release point. Further, the estimates are employed to improve the prediction of dose rates. The atmospheric dispersion phenomenon of radioactivity release is modelled using the Gaussian plume model. The Gaussian plume model is then employed for the calculation of dose rates. A variety of atmospheric and accident related scenarios for single source and multiple sources are studied in order to assess the efﬁcacy of the proposed ﬁlters. Statistical measures have been used in order to validate the performance of the proposed approaches.


Introduction
Nuclear power plants (NPP) and installations are potential sources of release of radionuclides into the atmosphere.Some radionuclides may release routinely and others may arise out of accidental situations.Radionuclides release estimation in the vicinity of a nuclear power plant is thus a crucial and regular activity.It plays a deciding role while performing various important tasks, such as emergency planning, sheltering, evacuation management, and in measuring environmental impact during normal or accidental scenarios [1].Particularly, at the initial release stage, a good estimate of the radionuclide spread is of paramount importance.These requirements make radionuclide release estimation a key regulatory problem.
The early release of radionuclides from nuclear installations forms radioactive plume in the atmosphere.These plumes are transported and dispersed by the winds and turbulence in the surroundings.The release of radioactivity can cause significant doses in the surroundings [2].For instance, these radionuclides may get deposited on vegetation and find their way to human metabolism [3].Environmental radiation monitoring detectors or dose rate monitors are placed in the vicinity of an NPP to measure the dose rates.These monitors provide on-line data, which can be used to estimate the source term and other relevant parameters [4].During radiological emergency, it is sometimes difficult to exactly know the amount of release or plume rise.Moreover, the data received from such measurements usually contain noise and affect their use in estimating the exact source of release [5].The accurate and timely prediction of radioactivity release is essential to evaluate the adequacy of emergency response measures, and to know its short-term or long-term environmental impact [6].
The environmental monitors placed around a nuclear installation employ radionuclide release information to measure dose rates [7].The calculation of dose rate is dependent on the plume shape and the distribution of radionuclide concentration.The dispersion and transport of radionuclides in the atmosphere is governed by the atmospheric dispersion modelling [8][9][10].At the early release stage, a rapid evaluation of the imminent radionuclide spread is required [11].Thus, the dispersion model should be fast, reliable, and easy-to-use, such as Gaussian models.The Gaussian plume dispersion (GPD) model has been widely used in the literature with different simplifications [12].The GPD model is suitable for modelling radionuclide release from point sources.It only requires a limited amount of input information and are usually time efficient during setting up and execution.The rendered results are trustworthy for near-field dispersion cases.Based on this rationale, the GPD model has been considered in this work to model radionuclide release in atmosphere and for the calculation of dose rates.
In the literature, different estimation techniques have been applied in order to solve problems of dispersion parameter estimation, such as predication improvement or data assimilation [13][14][15][16][17][18][19][20][21][22][23] and source term estimation or source inversion [24][25][26][27][28][29][30][31][32][33][34].A data assimilation technique based on Kalman filter has been proposed by Rojas-Palma et al. to improve the prediction capability of a real-time on-line decision support system [13].The work has been extended by Astrup et al. for off-site nuclear emergency management for the mesoscale atmospheric dispersion program [14].An Extended Kalman Filter (EKF) based approach has been suggested by Drews et al. for source term estimation in the case of short-range radionuclide release utilizing the offline recorded radiation monitoring data [24].However, their work did not consider improvement of gamma dose rates.Drews et al. further extended the application of the Kalman filter while using experimental online radiation monitoring data from 41 Ar releases [25].Santhosh et al. applied the Kalman filter technique to design part of a diagnostic system for accident management during radioactivity release [26].Ensemble Kalman filters (EnKF) [15,16] have been proposed by Zheng et al. to develop data assimilation techniques combining model predictions and measurements for the design of emergency response system using the Monte Carlo atmospheric dispersion model.EnKF [27] is further employed by for source estimation using multi-direction wind tunnel experiments.Modified EnKF have been proposed by Zhang et al. to improve predictions during nuclear accidents [17,18].However, the measurements considered are tracer concentration, which is different from the actual nuclear accident scenario.In some recent studies, particle filter based techniques have been proposed by Hiemstra et al. for model prediction and in forecasting the contamination [19,20].A preliminary study based on Unscented Kalman Filter (UKF) has been done in [28].An adaptive neuro-fuzzy inference system and linear regression filter are proposed for data assimilation [22].Application of other source estimation methods in real nuclear accident scenarios, such as the Chernobyl and Fukushima accidents [29][30][31][32], have been demonstrated.
Most of the work reported in literature employs estimation techniques for radionuclide release problem while a few studies attempted nonlinear estimation techniques.Particle filters generally encounter difficulties to generate meaningful samples in high-dimensional spaces.The EnKF belongs to a broader category of particle filters.In the case of EnKF, the accuracy of the state estimates is dependent on ensemble size which is heuristically chosen.The EnKF requires several equally weighted ensembles, which increases the computational time.On the contrary, the unscented and cubature Kalman filters (CKF) are computationally efficient filters in which the sample points are chosen deterministically and the number of sample points required is of the same order as the dimension of the system.The performance of unscented and cubature Kalman filters can be enhanced by careful tuning of their parameters.In this work, different nonlinear Kalman filter-based approaches are proposed in order to estimate radionuclide release activity.Particularly, EKF, UKF, and CKF techniques are proposed to estimate concentration of release, height of release, and wind speed from the measured dose rates obtained from radiation monitors.In addition, the estimation algorithms work to improve prediction of dose rates.The proposed algorithms are tested in a variety of atmospheric and beyond design basis scenarios.The efficacy of the proposed techniques has been validated using simulations in Matlab/Simulink.The performance of proposed algorithms are compared while using statistical measures of root mean squared error and sample correlation coefficients.The main contributions of the paper are listed below: The rest of the paper is organized, as follows: Section 2 describes dynamical models of atmospheric dispersion and dose rate calculation.Section 3 presents different nonlinear estimation techniques.Section 4 demonstrates the application of estimation techniques in radionuclide release estimation.Finally, conclusions are given in Section 5.

Modelling of Radioactivity Release
Radioactivity releases travel and disperse in the atmosphere and they are diffused with the wind in its direction.The movement of release is governed by the physical properties of radionuclide and the ambient atmosphere.Generally, the initial values of speed, momentum, and temperature of release are different from the ambient conditions.Due to the combination of these factors and buoyancy, the radionuclides release rise in the vertical direction until their effect get dissipated by mixing with the ambient atmosphere.This vertical rise of release is called plume rise and it can change the effective height of release, which denotes the summation of plume rise and stack height.Once the release enters the atmosphere, it is transported by wind and diffused by turbulent eddies, both during the plume rise and subsequently.Moreover, the turbulent motion of the atmosphere causes random movement of the release resulting in its progressive lateral and vertical spread by mixing with air.The combination of this transport and diffusion phenomena is called atmospheric dispersion [3,[8][9][10].

Atmospheric Dispersion Model
Accident management of an NPP is decided by the atmospheric dispersion model, which predicts the spatio-temporal diffusion of a radionuclide containment in the atmosphere.Atmospheric dispersion phenomenon in the case of release can be represented using different models.The GPD model [8] is adopted in this work due to its suitability for point sources in the case of plain and homogeneous terrain conditions, and due to its feasibility in modelling both short and long term impacts effectively.The GPD model is based on the assumptions that meteorological conditions are homogeneous and stationary with constant source and the radionuclide concentration is normally dispersed along the centreline [9].
The concentration of radionuclide activity is represented, as where x, y, and z denote downwind distance from the release point, crosswind distance from centreline of plume, and height above the ground, respectively.q, u, and h are release rate, wind speed, and height of radionuclide release respectively.σ y (x) and σ z (x) are plume dispersion parameters in y and z direction, respectively.They can be derived from observed wind speed fluctuations through similarity theory, or they can be inferred from Pasquill-Gifford type atmospheric stability classification [9].The Pasquill-Gifford stability categorisation is done while using wind speed together with the incoming solar radiation during daytime and cloudiness during night time.The stability of the atmosphere is a property that depends upon temperature profile and determines the vertical acceleration of atmospheric motions [10].The Pasquill-Gifford stability are derived using the Eimutis-Konicek relations of the form where the parameters a y , a z , b y , b z , and c z depend on the Pasquill-Gifford stability class which classify the amount of atmospheric turbulence present [8].

Dose Rate Model
The evaluation of dose rate for an atmospheric release of radioactivity is carried out at various stages of nuclear fuel cycle.It starts from the design phase of a nuclear facility and it continues to its operational phase for ascertaining compliance during normal operation and for emergency planning during accidents [3].
Dose rate computation from radionuclide release is related to meteorological parameters and the release itself [10].The dose rate absorbed at a receptor position is given by: where E γ is photon energy, ρ is mass density of air, and ζ (x, y, z) is radionuclide activity concentration at point (x, y, z) in the plume.The linear mass attenuation factor and energy absorption factor for photons in air are represented by µ and µ a , respectively.The photon yield or the number of photons per disintegration is assumed to be equal to one.The distance between plume point and receptor point is given by The dose build-up factor for Compton scattered radiation is

Nonlinear Estimation Techniques
A general non-linear discrete-time system is given by where u ∈ R n u , x ∈ R n x , and y ∈ R n y denote input, state, and output vector, respectively; f (.) and h (.) are known nonlinear vector functions; The covariance of process noise w ∈ R n w and measurement noise v ∈ R n v are assumed to be uncorrelated white with zero mean and covariances

Extended Kalman Filter
EKF is a well-known approach for nonlinear stochastic state estimation.It works by performing the first-order Taylor series expansion of the model.The implementation of the EKF is based on the computation of Jacobian matrices at each time step in order to determine the local linearized model of the system for which the Kalman filter equations can be employed [35].
The main steps of the EKF algorithm are described, as follows: 1.
The EKF is initialized using an initial estimate of state vector and error covariance matrix as For k ∈ 1, 2, . . ., N 2. Time Update (a) For the given nonlinear system (6), the Jacobian matrices are calculated as (b) Prediction Transformation: The a priori state and covariance are computed as (c) Observation Transformation: The mean and covariance of predicted output and cross-covariance between state and output are given as The Kalman gain is calculated while using the covariance matrices as (b) The state vector and error covariance matrix are updated as

Unscented Kalman Filter
The UKF is a derivative-free nonlinear Kalman Filter algorithm to avoid the analytic linearisation that is required by the EKF [36].The implementation of UKF is based on unscented transformation (UT), which is used to handle the nonlinearity in the transformation.The UT provides a statistical alternative to the analytical linearisation approach using Jacobian matrices, which is used in the EKF.The UT uses a small set of deterministically selected points, called sigma-points, in order to capture the probability distribution and generalizes to nonlinear system without the burdensome analytic derivation as in the EKF.The confidence level of these points from the assumed prior distribution is determined based on the selected scaling parameters for the UT.The scaling parameters affect the spread of the sigma points as well as the weight vectors that are used in reconstructing the a posteriori statistics.
The UT-based UKF algorithm is described, as follows: 1.
The UKF is initialized using the initial estimate of state vector and error covariance matrix, as For k ∈ 1, 2, . . ., N 2. Time Update (a) The UT transformation employs primary, secondary, and tertiary scaling parameters, in which the primary scaling parameter α determines the spread of the sigma-points around the mean, the secondary scaling parameter β includes information regarding the prior distribution, and the tertiary parameter κ minimizes the higher order errors of mean and covariance approximations.Using these three scaling parameters, an additional scaling parameter λ and weight vectors, η m and η c are defined as and (b) The Sigma points are computed using the initial state estimate ( xk−1 ) and covariance (P k−1 ) as where Γ is a 1 × n x vector of ones.(c) Prediction Transformation: Sigma points are fed to the non-linear state prediction as The a priori state and covariance are computed based on weighted averages of the transformed sigma-points as (d) Observation Transformation: Sigma-points are transformed while using observation function as The mean and covariance of predicted output and cross-covariance between state and output are given as

Measurement Update
(a) The Kalman gain is calculated using the covariance matrices as (b) The state vector and error covariance matrix are updated as (c) Sigma points are updated while using the a prior state estimate ( xk ) and covariance (P k ) as

Cubature Kalman Filter
The CKF is a recently developed derivative-free nonlinear Kalman filtering algorithm.The performance of the CKF gives nonlinear optimal approximation, filtering accuracy, and numerical stability.The CKF builds on the numerical-integration perspective of Gaussian filters [37].It uses spherical-radial cubature rules of third degree in order to evaluate Gaussian-weighted integrals encountered in the nonlinear Bayesian filter.The CKF calculates the means and covariances of state variables at a set of cubature points.
The CKF algorithm is described, as follows: 1.

Initialization
The CKF is initialized using an estimate of state vector and error covariance matrix, as For k ∈ 1, 2, . . ., N 2. Time Update (a) Cubature points are computed using the initial state estimate ( xk−1 ) and covariance (P k−1 ) as (b) Prediction Transformation: Cubature points are fed to the non-linear state prediction, as The a priori state and covariance are computed based on weighted averages of the transformed cubature-points as (c) Observation Transformation: Cubature-points are transformed using observation function as The mean and covariance of predicted output and cross-covariance between state and output are given as The Kalman gain is calculated while using the covariance matrices, as (b) The state vector and error covariance matrix are updated as (c) Cubature points are updated using the a prior state estimate ( xk ) and covariance (P k ) as

Application to Radionuclide Release Estimation
The nonlinear estimation techniques proposed in Section 3 are applied to estimate radionuclides release rate, wind speed, and effective height of release.Eight radiation detectors (D1 to D8), placed 1 km away from the source at the ground level, as shown in Figure 1.The source is placed at the origin indicated with a triangle.The mean wind direction is assumed to be constant at 30 degrees to x-axis in a counter-clockwise direction.The measurement data comprise of a set of dose rates from D1 to D8 detectors is generated using the GPD model discussed in Section II.The initial value of release rate, mean wind speed, and release height are assumed as 1 × 10 7 Bq/s, 4 m/s, and 100 m, respectively [28].Once the release has occurred and the dose rate have been recorded, the model can be used in its inverse mode to estimate the system states and input parameters.This task is performed by coupling the model to the measurements.For the purpose of estimation, the GPD model can be represented in state space form using a random walk model.It is suitable if the unknown parameters are time-varying.The random walk model is given by, where s k and o k denote system states and measurement, respectively.The vectors w and v represent process noise and measurement noise, respectively.q, u, and h denote the release rate, mean wind speed, and effective height of release, respectively.D 1 to D 8 respectively represent dose rates measured at detectors.The proposed algorithms are then applied to estimate different states using environmental radiation monitoring data collected close to the release point.Countermeasures are applied at 300 s, which gradually brings back the release rate to its initial value in the next 500 s.The release is then kept constant afterwards.The dose rate detectors measure the amount of release in terms of doses and they are assumed to be added with white Gaussian noise.These noisy measurements are then used in designing estimation algorithms.The estimated dose rates using EKF, UKF, and CKF are shown with measurements at eight different detector locations in Figure 2. The estimates of radionuclides release rate, mean wind speed, and effective height of release with true values are shown in Figure 3.All of the estimation algorithms able sufficiently to estimate the true values of the dose rates.The UKF and CKF clearly outperform the EKF during state estimation.The EKF shows perturbations during tracking the state variation.The UKF and CKF algorithms are able to track the step and ramp variations change in release rate in addition to other states and measurements.The estimated states from UKF and CKF are found to be in good agreement with the true values.

Variation in Mean Wind Speed
The efficacy of the proposed approaches is validated by testing another beyond design basis scenario of a wind storm [28].The wind is regarded to be time-dependent, in which the speed varies with time while the direction remains fixed.The wind speed is considered to be 4 m/s for the first 200 s.It is then linearly changed to 24 m/s in 500 s and maintained there for the rest of the duration.The noisy measurements of dose rates at different detectors are collected and then used in estimation algorithms.The dose rates estimated using EKF, UKF, and CKF are shown with measurements at eight different detector locations in Figure 4. Estimates of release rate, wind speed, and height of release are shown in Figure 5.All three techniques are able to estimate the dose rates close to their true values.
The UKF and CKF perform better than the EKF during state estimation.The EKF shows comparatively more perturbations in the estimation of mean speed and height of release than the UKF and CKF.All of the algorithms are able to track the ramp variation in wind speed.The estimates are in good agreement with the true measurements.

Variation in Effective Height of Release
The effective height of release is a sum of the physical stack height and the plume rise.The physical stack height remains constant, while the plume rise can vary, which can alter the effective height of release [3].During radionuclides release generally the plume is of non-buoyant nature, in which the radionuclide discharge temperature is near ambient.During an accident, the discharge temperature can increase and it will affect the plume rise, which, in turn, change the effective height of release.Here, a scenario is considered, in which the effective release height is step changed by 25 m from its initial value at 100 s.It then linearly decreased to its initial value during 300 to 800 s and maintained at its initial value for rest of the duration [17].The noisy measurements of dose rates at different detectors are collected and used in estimation algorithms.The estimated dose rates using EKF, UKF, and CKF are shown with measurements at eight different detector locations in Figure 6. Figure 7 shows the estimates of radionuclides release rate, mean wind speed, and height of release with true values.
The UKF and CKF are effectively able to predict the dose rates at D1, D2, D3, and D5 detectors, while the predictions at detectors placed away from plume dispersion contain perturbations.The UKF and CKF show initial variation in the estimation and takes some time before tracking the true values.The state estimates from EKF contain perturbations.It can be observed that all the three techniques are able to track states, except some initial variation.The algorithms are able to handle the step and ramp variation in release height in while tracking other states and measurements.

Variation in Atmospheric Condition
This simulation studies the effectiveness of estimation techniques during variation in atmospheric condition.The atmospheric condition changes due to variation in incoming solar radiation, night time cloud cover, and surface wind speed [10].For instance, solar radiation increases atmospheric instability through warming of the Earth's surface, whereas clear nights push conditions toward stability [3].Here, the atmospheric condition are considered to change due to atmospheric turbulence in random fashion.Figure 8 shows the variation in atmospheric turbulence parameters.Because of the variation in atmospheric condition, the wind speed is also assumed to vary in a random fashion from its initial value.The noisy measurements of dose rates at different detectors are collected and then used in estimation algorithms.Dose rates estimated using EKF, UKF, and CKF are shown with measurements at eight different detector locations in Figure 9. Estimates of release rate, wind speed, and height of release are shown in Figure 10.The UKF and CKF techniques perform better than the EKF during prediction as well as during estimation.The EKF shows comparatively more perturbations in the estimation of mean speed and height of release than the UKF and CKF.Further, the dose rates estimated by UKF and CKF are in close agreement with the true value than the EKF.

Multiple Release Points
To analyse the effectiveness of the estimation techniques when there are multiple sources of radionuclide release estimation, a case study that comprised of three release sources has been considered.These three sources represent three different units with locations (x s , y s , z s ) and release concentrations ζ s (x s , y s , z s ) for s = 1, 2, 3 and their placement is shown in Figure 11 by the points labelled S1 to S3.The initial values of release rates from sources S1 to S3 are considered as 1 × 10 7 Bq/s, 1.2 × 10 7 Bq/s, and 1.3 × 10 7 Bq/s, respectively.In a worst case scenario, where an earthquake affects all three units simultaneously, there would be radionuclides release from all three sources.The total concentration that results from all three sources is given by the sum of individual concentration.For simulation purposes, here, the total release rate variation is considered, as follows.The total release rate is varied to two times of its initial value from 50 s to 150 s and maintained there for the next 200 s.Countermeasures are applied at 350 s, which gradually bring back the total release rate to its initial value in the next 500 s.The release is then kept constant afterwards.The noisy measurements of dose rates at different detectors are collected and then used in estimation algorithms.Estimates of dose rates using EKF, UKF, and CKF are shown with corresponding measurements in Figure 12. Figure 13 shows estimates of release rate, wind speed, and height of release.The estimation algorithms estimate the dose rates and states with variations in tracking.It is to be noted that the measurements comprised of multiple sources while the estimation algorithms modelled only a single source.All three techniques predict the dose rate with some overshoot and take some time in predicting the measurements.The UKF and CKF are found to be perform better than the EKF during state estimation.

Statistical Performance Assessment
The performance of the estimation algorithms is quantitatively assessed using different statistical measures.The root mean squared error (RMSE) is calculated between the reference value and estimated value.The RSME is given by where N is the length of data.r k and pk , respectively denote the reference and estimated values.
The RMSE for different scenarios is calculated for measurements and states and are listed in Tables 1 and 2, respectively.The bold values represent the best estimate among all.It can be noted that the RMSE values by the UKF and CKF techniques are much lower than that obtained by the EKF approach.From the Tables 1 and 2, it can be noticed that the estimation accuracy for the algorithms is similar for cases 4.1.2,4.1.3,and 4.1.4.The non-linear estimation algorithms that are based on statistical linearisation track the variation in release rate, effective release height, and mean wind speed better than the analytical linearisation technique.The estimation and prediction capabilities can be assessed by calculating the fraction bias (FB).It is given by where rk and pk denote the mean values of r k and p k , respectively.The FB for different scenarios is calculated for states and measurements and they are listed in Tables 3 and 4, respectively.The positive and negative values of FB denote underestimation and overestimation, respectively, whereas the zero value denotes correct estimation.It can be noted that the value of FB is lowest in the case of UKF during state estimation, while all of the algorithms perform similarly during prediction.Further, the prediction capability of the estimation algorithms in predicting the true value from the noisy measurements can be assessed by estimating the sample correlation coefficient.It is given by The sample correlation coefficient is calculated between the reference value and predicted values.Table 5 shows the value of correlation coefficient for different estimation algorithms in predicting dose rates.It can be noted that all of the estimation algorithms predict the measurement accurately.The sample correlation using UKF and CKF are closer to one as compared to that obtained by EKF.

Tuning of Covariance Matrices
The process noise covariance matrix Q k and the measurement noise covariance matrix R k are constants and they do not depend on k.In the literature, various noise covariances matrix estimation approaches, such as the maximum-likelihood methods, correlation methods, covariance matching methods, and Bayesian methods are proposed.The various methods differ in assumptions related to the system model, underlying ideas and principles, properties of the estimates, and number and essence of the design parameters.Any of the above-mentioned methods can be employed to estimate the covariance matrices Q k and R k [38].Selecting the optimum value of Q k and R k is crucial for the successful application of estimation algorithms however, the optimal selection is arduous.If Q k is chosen less than the actual value, the estimation algorithm will weigh the model more over the measurements.This can affect the performance during the correction of system states.In contrast, if Q k is chosen more than the actual value, then the increased values of the state covariance matrix will give noisy and uncertain estimates.A tuning approach that is based on computing the variance of innovations has been adopted here.The values of Q K are varied for a fixed value of R k and the variance of innovation is calculated.The values of Q k and R k are selected, such that there is no significant variation in the variance of innovations on further reducing Q k for a fixed value of R k .

Discussion
The proposed nonlinear Kalman filter techniques simultaneously estimate different parameters related to radioactivity release and improve the forecast of dose rate measurement.The unscented and cubature Kalman filters are found to have higher accuracy and robustness for non-linear models than the extended Kalman filter.Instead of linearising the functions as is done by the extended Kalman filter, the unscented and the cubature Kalman filters use a set of points and propagate these points through the actual non-linear function.These points are chosen, such that the lower and possibly the higher-order moments match those of the Gaussian random variable.The mean and the covariance can be recalculated from the propagated points, yielding more accurate results when compared to the ordinary function linearisation.Different scenarios have been studied to analyse the performance of the proposed estimation techniques.In case of increment in release rate, the dose rate monitors suddenly detect large dose rates.This sudden increment in dose can also be due to a reduction in wind speed.Thus, the estimation algorithm initially estimate a reduction in wind speed.However, due to the prediction and correction nature of Kalman filters, they are able to figure out the exact reason of increased dose rate and, thus, are able to estimate estimate the correct states.Then, due to countermeasure the release rate decreases, which leads to a decrement in the dose rates.This decrement in dose can also be due to an increment in wind speed which take away the plume.After prediction and correction, the estimation algorithms are able to estimate the correct states which leads to increment in the estimate of wind speed.In the case of increment in wind speed, the dose rate monitors suddenly detect low dose rates.This low dose rates can also be due to a reduction in release.Thus, the estimation algorithm initially estimate a reduction in release rate.Now, when the wind is strong, the released plume will quickly leave the monitored area, and the dose rate monitors fail to capture the plume.Moreover, the estimate of release heights are slightly different than the steady state condition due to poor measurements.In case of increment in plume rise, it becomes difficult for the plume to reach the ground, which leads to the insufficiency of the measurements.This will lead to sudden increment in wind speed.However, after prediction and correction steps with more dose rate data, the Kalman filter will converge to its steady state value.Thus, the accuracy of the nonlinear Kalman filter-based techniques is dependent on the accuracy of the model and the measurements.The employed Gaussian plume dispersion model is suitable in the range of tens of km for plain and homogeneous terrain conditions.If the measurements are sufficiently clean then the accuracy of the proposed technique does not differ much in between these characteristic lengths.However, if the measurements are heavily noisy then the performance of the algorithms will be affected.The selection of process and measurement noise covariance matrices further affect the performance.Thus, a careful selection of Q k and R k matrices is important.In addition, the considered Gaussian plume dispersion model implies some limitations on model's applicability for detailed studies.The presented model does not try to represent wind direction, which is assumed as constant.Any variation in the wind direction will act as modelling uncertainty and it will affect the accuracy of the state estimation.For a better prediction, the wind direction can be considered as another state of the model.For the detection of source location, the estimation algorithm can be combined with any backtracking algorithm to locate the centre of the source by extrapolating the states backwards [39].This will further increase the accuracy of the release rate.For locating multiple release sources, the nonlinear estimation algorithm can be combined with other techniques, such as the optimization methods that are used in the field of collective robotics [40,41].

Conclusions
Radioactivity release estimation is an important task during normal or accident scenarios.In this work, nonlinear Kalman filter-based approaches have been proposed for radionuclide release activity estimation.Specifically, Extended Kalman Filter, Unscented Kalman Filter, and Cubature Kalman Filter algorithms have been proposed in order to estimate radionuclides release rate, mean wind speed, and effective height of release.The algorithms are further employed to predict dose rates at different detectors.The estimation algorithms are found to reliably estimate the source information.It has been found that the UKF and CKF give better estimates than the EKF algorithm.The technique employed the Gaussian plume model of atmospheric dispersion for the computation of doses.The effectiveness of the estimation algorithms has been validated for different simulated scenarios.It has been observed from the simulation results that the estimation algorithms perform satisfactorily.Different statistical measures have been computed for analysing the performance of the proposed approaches.In the future, the nonlinear estimation algorithms will be designed in order to incorporate mobile sensors placed on autonomous vehicles in addition to the fixed sensor arrangement to further enhance the estimates.

Figure 1 .
Figure 1.Placement of different dose rate monitors (distance in km).

4. 1 .
Case Studies of Different Scenarios 4.1.1.Variation in Radionuclides Release Rate A simulation analysis is performed in order to study the effectiveness of proposed techniques during an explosion scenario [17].To represent an explosion, it is considered that the release rate is step changed to two times of its initial value at 100 s and maintained there for the next 200 s.

Figure 2 .Figure 3 .
Figure 2. True, measured, and estimated dose rates at detectors D1 to D8 for variation in release rate.

Figure 4 .Figure 5 .
Figure 4. True, measured, and estimated dose rates at detectors D1 to D8 for variation in wind speed.

Figure 6 .Figure 7 .
Figure 6.True, measured, and estimated dose rates at detectors D1 to D8 for variation in release height.

Figure 9 .Figure 10 .
Figure 9. True, measured, and estimated dose rates at detectors D1 to D8 for variation in atmospheric condition.

Figure 11 .
Figure 11.Placement of multiple release sources in a triangular arrangement (distance in km).

Figure 12 .Figure 13 .
Figure 12.True, measured, and estimated dose rates at detectors D1 to D8 for multiple release sources.

Table 3 .
Fractional Bias between true and estimated states for EKF, UKF, and CKF.

Table 4 .
Fractional bias between true and estimated dose rates for EKF, UKF, and CKF.

Table 5 .
Sample correlation coefficients between true and estimated dose rates for EKF, UKF, and CKF.