Uncertainty and Sensitivity Analysis at Low Value of Determination Coe ﬃ cient of Regression Analysis: Case of I-129 Release from RBMK-1500 SNF under Disposal Conditions

: As in other nuclear countries, the operation of the Ignalina nuclear power plant in Lithuania has led to the accumulation of around 22 thousand assemblies of spent nuclear fuel (SNF). The development of geological disposal program involves an iterative assessment of the system safety supported by scientiﬁc research on radionuclides migration and related processes. This study focused on the application of Contribution to the Sample Mean (CSM) and Contribution to Sample Variance (CSV) methods to complement the uncertainty and sensitivity analyses of the time-dependent ﬂux of I-129 from the engineered barriers of a conceptual disposal facility for RBMK-1500 SNF (RBMK is abbreviation of “High Power Channel-type Reactor” (in Russian)). The analysis was performed using a MATLAB platform (8.0.0.783 (R2012b), MathWorks, MA, USA). The mean and variance ratios derived from CSM and CSV plots were applied to estimate the e ﬀ ect of reduced uncertainty range on mean ﬂux and its variance, and the uncertainty analysis was also complimented. Increasing the lower bounding value of defect size enlargement time range to 4.6 × 10 4 years would lead to a lower mean ﬂux until 5 × 10 4 years after repository closure. Later on (up to 1 million years after repository closure), the only reduction of the upper bounding value of the SNF dissolution rate range would a ﬀ ect a decreased mean ﬂux. Tlarge—defectsizeenlargementtime, SNFDR—SNFmatrixdissolutionrate, around the canister, IRF—instant De—e di ﬀ usivity of I-129 in bentonite, CSM—contribution to the sample mean, CSV—contribution to the sample variance.


Introduction
As in other nuclear countries, the operation of the Ignalina nuclear power plant in Lithuania has led to the accumulation of around 22 thousand assemblies of spent nuclear fuel (SNF). Each country's responsibility for the safe management of its SNF is acknowledged worldwide. Within the European Union (EU), directive 2011/70/EURATOM [1] contains the provision for every member state (country) to be responsible for the implementation of a safe and sustainable solution for SNF and radioactive waste management and disposal.
In the case of SNF disposal, one of the key issues is the demonstration of the long-term safety of the disposal system. Safety assessment covers the understanding and analysis of a broad range of processes in relation to disposal system evolution, radionuclide distribution, and the impact on human health and the environment. Lack of reliable data leads to overestimating environmental impact and modeling results with large uncertainties.
The role of sensitivity analysis (SA) in line with uncertainty analysis (UA) is widely recognized as an effective tool to understand the influence and to rank the most important parameters for reducing criteria. However, there are several other indicators, which would be assessed and discussed within the safety case as a complement to dose/risk. These includes performance indicators, such as radionuclide transfer times, radionuclides concentrations in the near field (EBS and part of surrounding host rock), radionuclide fluxes in the near field, characteristics that control "dilution" in time and space (e.g., waste-form dissolution or release rates, canister failure rate), etc. [20].
The most recent results of probabilistic uncertainty and sensitivity analysis of I-129 release through the engineered barriers of a conceptual disposal facility for RBMK-1500 SNF were reported in [21]. In the study [21], the peak flux and the time-dependent flux of I-129 were analyzed as variables of interest (model output), and the importance ranking was performed based on regression/correlation coefficients. Sensitivity measures in terms of standardized regression coefficients (SRC) appeared to be sufficient for importance ranking for peak flux. Meanwhile, for time-dependent flux, it was concluded that regression analysis might not be reliably ranking the parameter importance as the coefficient of determination (R 2 ) became very low (<0.3) 10 4 -10 5 years after repository closure, indicating strong nonlinearity between the flux and input parameters. For these reasons, this study was devoted to the analysis that would determine the most important parameter for the period where the regression analysis was unreliable.
This study estimated the contribution of each parameter to the mean and variance of model output. Based on derived CSM and CSV plots, the ranking of parameter importance was established for the time-dependent flux of I-129 through the engineered barriers of a conceptual geological repository for RBMK-1500 SNF disposal. It also applied CSM-and CSV-based measures to explore the effect of parameter uncertainty reduction. The effect of reducing the uncertainty of input parameters was assessed in terms of the mean and variance. The advantage of CSM and CSV tools is no necessity of additional model runs. Also, the derivation of complemental results regarding importance ranking was done, where the regression-based method was unreliable, and indications about the effect of newly justified input parameter uncertainty range on the model output were available at no additional computational cost.

Radionuclide Transport Model for RBMK-1500 SNF Disposal
For studying the long-term performance of a disposal system, commonly involves using numerical models for the analysis of events and processes that take place in (or could affect) the system. There are several processes (thermal, chemical, physical, etc.) influencing the disposal facility (repository) and its surrounding environment necessitating a complex analysis of the processes and their interactions. Thus, the analyzed model output differs from study to study (i.e., peak flux [22], peak time [23], maximum concentration [5], peak dose [24], the amount of radionuclides being released at particular times after repository closure [25], dose over time [26], dose conversion factors [6]). The uncertainty and sensitivity of physical parameters that impact contaminant transport to the environment are considered as a model output to be analyzed (e.g., pressure in the borehole [27]). Nevertheless, the main aim of such studies is to help to assess radionuclide release and migration within the disposal system and its surrounding environment to evaluate the DGR system's safety.

Disposal System
There is an international consensus that geological disposal facility located at significant depth in stable rock formations offer a reliable and sustainable solution for the long-term management of high-level radioactive waste. The long-term safety of DGRs is ensured by applying the multi-barrier concept. The multi-barrier concept involves the consideration, justification, and implementation of a series of natural (host rocks) and engineered (waste form, package, backfill, and seal materials) barriers [28].
For the disposal of RBMK-1500, SNF in Lithuania DGR constructed in the crystalline rock could be considered. The radionuclide migration is assessed in consideration of engineered barriers, including SNF matrix itself, a copper canister, and surrounding bentonite ( Figure 1). More details on the repository concept and radionuclide release model are provided in [21]. The main function of the engineered barrier system of DGR is to prevent and/or limit the release of radionuclides into the surrounding geological environment. As indicated by [29], it is impossible to ensure complete containment for hundreds of thousands of years after repository closure; thus, the eventual release of radionuclides must be minimized after the loss of integrity of the engineered barriers surrounding the SNF matrix. Once the contaminants are released from the engineered barriers, they could be dispersed in the surrounding geological environment and transferred into the biosphere. The main assumptions considered during the assessment of I-129 transfer in the near field region (EBS) were as follows [21]:

•
The bentonite material surrounding copper canister will be fully saturated by the groundwater by the time I-129 release from the SNF assemblies; • As soon as a small initial canister defect becomes large, the void space within the canister (void between SNF assemblies and the channel in the canister insert) will be filled with the groundwater (app. 0.5 m 3 ); • Mechanisms (instant release of a part of the inventory and congruent release of the rest part of the inventory from degrading SNF matrix) take place in the canister and contribute to the radionuclide flux from the canister; • Radionuclides released from SNF assemblies interact with a limited amount of water inside the canister, and dissolved radionuclide are transported from the canister in liquid form (mainly by diffusion); • I-129 is released through the bentonite barrier and diffuses into the water flowing in a (conceptual) fracture intersecting the disposal tunnel.
The selected uncertain input parameters and the properties of the parameter distributions for the probabilistic transport analysis are provided in Table 1. The model realized using the computer code AMBER [34] assessed time-dependent radionuclide release from the SNF matrix, dissolution, radioactive decay, and contaminant transport by diffusion through engineered barriers. The uncertainty of the main transport-related parameters, including defect size enlargement time, was characterized by probability density function for each parameter. The model output was a large number of the time-dependent flux values over an extended period (10 3 -10 6 years after repository closure), which allowed for the evaluation of the mean flux, the quantiles, and the distribution of the peak flux, as illustrated in Figure 2. The parameter values were sampled randomly (the Monte Carlo method). In total, 1000 simulations were performed with sampled parameter values. The results of regression analysis (Standardized Rank Regression Coefficient (SRRC) as sensitivity measures) reported in [21] for time-dependent radionuclide flux were the main motivators for the current study and are presented in Figure 3.
Regression model for the analysis of model output dependence on the input parameters is expressed as follows: y m , x i,m are the average values of model output Y = (y 1 , . . . , y N ), input parameter X i = (x 1i , . . . , x Ni ), S y , S x i are respective standard deviations of model output () and input parameter ( are standardized regression coefficients, and N-number of model runs. SRRC is obtained through regression analysis on rank-transformed data. The rank-transformation is replacement of the original values of the input variables and the model output by their rankings (ranking 1 for the smallest value). This technique is used to linearize monotonic nonlinear relations so that linear regression analysis can be applied to the rank-transformed data [35]. . Standardized rank regression coefficients (SRRC) as parameter importance measures and the coefficient of determination R 2 over the analyzed period for I-129 flux as presented in [21] (R 2 -coefficient of determination, Tlarge-defect size enlargement time, SNF DR-SNF (spent nuclear fuel) matrix dissolution rate, Qeq-equivalent groundwater flow rate around the canister, IRF-instant release fraction, De-effective diffusivity of I-129 in bentonite).
As previously mentioned, the SA of radionuclide migration through the engineered barriers of a conceptual geological repository for RBMK-1500 SNF disposal was performed [21]. The maximal (peak) flux and the time-dependent flux of I-129 were analyzed, and the SA was done based on regression/correlation coefficients. The regression analysis of I-129 peak flux yielded a determination coefficient (R 2 ) of 0.77, showing the relationship between the peak flux of this nuclide and the main transport-related parameters was not strongly linear. However, the standardized regression coefficients could still be used to rank parameter importance. For the peak flux of I-129, the effective diffusion coefficient in bentonite (De) and the instant release fraction (IRF) were found to be the most important input parameters.
Meanwhile, the regression analysis results for the time-dependent flux of iodine indicated that given parameters impact varied with different time steps. Effective diffusivity of I-129 in bentonite (De) and the instant release fraction (IRF) was reported as the most important parameters from 10 3 -10 4 years after repository closure. Within 10 4 -10 5 years after closure, the importance of De and IRF was significantly decreased, and the time when initial defect increases became important. Based on the regression analysis, defect size enlargement time (Tlarge) became the most important parameter after approx. 2 × 10 4 years after repository closure.
However, R 2 became very low (<0.3) during this period, indicating strong non-linearity between the flux and the parameters; thus, the regression analysis of the rank-transformed data might not provide reliable results for determining parameter importance.
For these reasons, this study was devoted to the analysis that would determine the most important parameter for the period where the regression analysis was unreliable. The contribution of each parameter uncertainty to the mean and variance of model output was estimated through CSM and CSV plots, and conclusions of parameter importance could be derived then.

Contribution to Sample Mean (CSM)
The CSM plot was introduced by [7] and further developed by [14]. To create a CSM plot, a random (or quasi-random) sample of size N for inputs Xi and corresponding model outputs Y need to be generated. The dependence (function) for a given input parameter Xi is then derived with the application of the following steps as provided in [14]: The randomly (quasi-randomly) generated values of X i from its probability distribution function (PDF) are sorted in ascending order, generating the series of values The output values y i corresponding to sorted input X i are obtained too y i,1 , y i,2 , . . . , y i,k , . . . , y i,N , k = 1, . . . , N; 3.
An ancillary variable M i is calculated, whose values are calculated from y i,q as i.e.: The function CSM X i is obtained by normalization of M i , i.e., dividing the values m i,k by the sample mean of model output y m :

5.
Cumulative relative frequency (cumulative fraction) of the sorted input parameter X i lies in the interval [0,1] and could be calculated: 6.
Then, the function CSM Each data point of q (k) X i corresponds to a particular value of input parameter value in series x i,1 , x i,2 , . . . , x i,k , . . . , x i,N (values sorted in ascending order). Each data point (q ) represents the fraction of the output mean due to any given fraction of the values of the input lower than that of corresponding to q (k) X i [14]. For example, the fraction of output mean at point (q X i = 0.25) is determined by the 25% of the lowest values of the parameter X i . CSM plot allows the contribution to the sample mean of the analyzed output Y to be assessed over a certain range of input parameter X i values. For example, the impact of 20% of the lowest X i values on the mean of Y could be analyzed over the range [0,0.2] on the x-axis. The impact on the mean is then assessed from the y-axis. The contribution to the output's mean could vary along with the range of the input parameter's values, which is reflected by the CSM curve shape.
The ranking of the most important input parameters with regard to the output mean is based on the greatest distance (D m ) of the CSM plot from the diagonal. The CSM curve could cross the diagonal several times. According to [14], the ranking of parameter importance could be based on the sum of maximum distances (in absolute values) of the CSM curve from the diagonal (D m ).

Contribution to Sample Variance (CSV)
Contribution to sample variance is an extension of CSM first introduced by [15]. For CSV plot, a random sample of size N for inputs X i and corresponding model outputs Y must again be generated. The dependency (function) for a given input parameter X i is obtained using the following steps according to [15]: The output mean y m is computed, and each value of output Y is transformed by subtracting the mean value y m. The transformed output Yt has zero mean value.

2.
The randomly (quasi-randomly) generated values of input X i are sorted in ascending order, generating the series of values . . , N, and the corresponding set

3.
Function CSV is obtained as follows

4.
Cumulative relative frequency (cumulative fraction) of the sorted input parameter X i lies in the interval [0,1] and could be calculated by Equation (4).
Then, the function CSV X i is plotted versus the cumulative relative frequency of q X i . CSV also lies in the interval [0,1]). The plot consists of pairs (q CSM or CSV curves that plot close to the diagonal indicate that the contribution to the mean or the variance is equal with each value of the input parameter being analyzed. For example, the CSV value at quantile q = 0.1 provides an estimate of the model output variance due to 10% of the smallest values of the input parameter. According to [14], if at some regions of the plot, the CSM (CSV) curve is increasi ng quickly, the contribution to the output mean (variance) is large. Where the CSM (or CSV) curve is flat, this means that the contribution to the output mean (or variance, respectively) will be small.
Besides identifying the most important input, CSM and CSV methods provide insight into which part of the input parameter uncertainty range affects uncertainty in model outputs (its mean and variance).

Revised Mean and Variance Ratio Functions
To complement the UA and SA of I-129 release from the conceptual RBMK-1500 SNF disposal facility, a part of the current study was devoted to getting insights on the effect of parameter uncertainty reduction with application of revised mean and variance ratio functions. Authors of study [36] revised the variance ratio proposed by [15] and introduced revised mean and variance ratio functions (HM and HV, respectively), to be derived as follows: The mean y m of model output Y is computed.

2.
The generated values of the input X i are sorted in ascending order, generating the series of . . , N, and the corresponding set of model output values

3.
Cumulative fractions of the input parameter q 1 , q 2 lying in the interval [0,1] are defined: 4. The mean ratio function HM i and variance ratio function HV i for q 1 ∈ [0, 1], q 2 ∈ [0, 1] can be estimated by the following expressions These ratio functions allow direct observation of change (for example, decrease) of mean and variance of model output due to reduced parameter uncertainty. The effect on model output could be explored on the graph for any pair of selected new q 1 and q 2 .

Results and Discussion
The results are arranged in three subsections. The extension of SA with the means of derived CSM and CSV plots for I-129 flux is presented first. The second subsection presents the results and discussion about the effect of parameter uncertainty range reduction on the uncertainty of the peak flux in terms of statistical uncertainty measures, such as mean and variance. Finally, the analysis of the effect of parameter uncertainty reduction on the uncertainty of flux over time (in terms of the mean) is presented.
At time t 1 and times t 6 -t 9 , the most important parameters for the mean flux were clear (the CSM of these parameters had the largest distance from the diagonal D m ): effective diffusivity of I-129 in bentonite (De) was the most important parameter at the very beginning of the release, and the SNF dissolution rate dominated in the long-term. This correlates well with the results of the regression analysis presented in Section 1.1.2.
The CSM plots showed that at certain points in time during between 5 × 10 3 -5 × 10 4 years after closure (t 2 , t 3 , t 4 , t 5 ), the defect size enlargement time was important (see Figure 4). The importance of the defect size enlargement time was observed at earlier points in time (at 5 × 10 3 years) in comparison to what was observed from the regression analysis (at 2 × 10 4 years after the closure). As previously discussed (Section 1.1.2), the coefficient of determination (R 2 ) decreased significantly during this period, and the regression-based ranking was considered as not reliable enough.
The CSM plots where the curve is close to the diagonal (where equivalent groundwater flow rate around the canister (Qeq) and all remaining parameters except for the SNF dissolution rate at 200 thousand years post repository closure) revealed a similar contribution to the mean flux over any quantile range of parameter values. Thus, these parameters could be fixed at any value within the range of the uncertainty without any effect on the mean flux at the corresponding time (Figure 4).   The actual change of the mean (variance), which could be achieved by reducing the parameter value range to X i (q 1 , q 2 ) (where q 1 < q 2 < 1), has been presented and discussed in the next sections using revised mean ratio functions. Figure 5 presents the CSV plots that indicate the influence of each parameter on the flux variance at different time points. The results showed the complexity of the model and varying contribution (of different parameters) to the model output variance over a simulated period (up to 1 million years after repository closure). value range to i X (q1, q2) (where q1 < q2 < 1), has been presented and discussed in the next sections using revised mean ratio functions. Figure 5 presents the CSV plots that indicate the influence of each parameter on the flux variance at different time points. The results showed the complexity of the model and varying contribution (of different parameters) to the model output variance over a simulated period (up to 1 million years after repository closure). At the beginning of radionuclide release from the engineered barriers (t1 = 1 × 10 3 years), the largest contribution to the flux variance came from De (effective diffusivity of I-129 in bentonite), as observed from its curve (Figure 5a). This is similar to what was observed from the CSM plots (regarding the contribution to the mean flux). At later time points, the impact of this parameter was less significant, and the time of defect size enlargement (Tlarge) became the most significant contributor (t2-t6, ~5 × 10 3 -10 5 years after the repository closure). Based on the CSV plots, parameter Qeq had a more significant influence on variance than on the mean flux up to 5 × 10 3 years after repository closure.
Uneven contribution to the variance was observed for the SNF dissolution rate (SNF DR). As for the Tlarge parameter in the CSM plots, the SNF dissolution rate on the CSV plot crossed the diagonal, indicating that the smallest and the largest values of this parameter contributed more to the output variance than the values from the middle part of the quantile range. Table 2 provides the most significant parameters for time-dependent radionuclide flux into geosphere based on the CSM and CSV plots. In the majority of cases, the ranking for both methods was the same, except for t = 1 × 10 5 years after closure when the most significant contribution to the mean flux came from SNF DR, and the most significant contribution to flux variance came from At the beginning of radionuclide release from the engineered barriers (t 1 = 1 × 10 3 years), the largest contribution to the flux variance came from De (effective diffusivity of I-129 in bentonite), as observed from its curve (Figure 5a). This is similar to what was observed from the CSM plots (regarding the contribution to the mean flux). At later time points, the impact of this parameter was less significant, and the time of defect size enlargement (Tlarge) became the most significant contributor (t 2 -t 6 ,~5 × 10 3 -10 5 years after the repository closure). Based on the CSV plots, parameter Qeq had a more significant influence on variance than on the mean flux up to 5 × 10 3 years after repository closure.
Uneven contribution to the variance was observed for the SNF dissolution rate (SNF DR). As for the Tlarge parameter in the CSM plots, the SNF dissolution rate on the CSV plot crossed the diagonal, indicating that the smallest and the largest values of this parameter contributed more to the output variance than the values from the middle part of the quantile range. Table 2 provides the most significant parameters for time-dependent radionuclide flux into geosphere based on the CSM and CSV plots. In the majority of cases, the ranking for both methods was the same, except for t = 1 × 10 5 years after closure when the most significant contribution to the mean flux came from SNF DR, and the most significant contribution to flux variance came from Tlarge. The effect of SNF DR on the variance was also observed to be uneven, and the smallest and the largest values of SNF DR contributed most to variance; however, this was not the case for the CSM plot.

Effect of Parameter Range Reduction on Mean and Variance of I-129 flux
As mentioned above, the derivation of CSM and CVM plots provided an opportunity to assess the range of the input parameter values to focus on to affect the output mean or variance. Several cases were analyzed in terms of analysis of the effect of uncertainty reduction. Cases were devoted to exploring the effect of five parameters' range reduction: 3D mean (variance) ratio HM (HV) plots enable the assessment of any effect on the mean (variance) for a range of parameter values X i (at q 1 , at q 2 ), satisfying the condition q 1 ≤ q 2 (see Section 2.3 for definition).
Scheme for the exploration of the effect of the input parameters uncertainty reduction on model output with the means of ratio HM for each parameter (while keeping the other parameters at their original ranges) is presented in Figure 6. The same procedure is valid in the case of HV. Figure 7a-e present the 3D plots of the revised mean ratio (HM) for all parameters at t 5 = 5 × 10 4 years corresponding to the time of maximal mean flux. As the plots demonstrate, for defining reduced ranges for parameter X i values (with assignment of the new lowest value of the range to be equal to value at selected q 1 and the highest value to be equal to value at selected q 2 ), the change of new mean (variance) could be estimated from the ratio HM (HV). If the ratio is >1, then the new estimate would be higher than that observed with the original parameter uncertainty range. If the ratio is <1, then the mean (variance) of the model output with a reduced parameter uncertainty range would be lower in comparison to the original one. As follows from Figure 7, increasing the minimal value of parameter Tlarge to be equal to q 1 > 0.7 and the decrease of the maximal value to be equal to q 2 ≥ q 1 would result in significant mean reduction (HM < 1). But assigning the Tlarge values at q 1~q2~0 .65 would result in increasing the mean by a factor of~6. The 3D plot of HM(q 1 , q 2 ) over q 1 and q 2 for the parameter SNF DR (Figure 7b) showed a trend leading to an increased mean by a factor of~3, in the case of a significant increase in the lower bounding value (assigning it equal to value at q 1 between 0.7 and 1). For the same parameter, decreasing the upper bound and increasing of the lower bound, a decrease in the mean flux would be expected. A similar trend was observed for instant release fraction (IRF). The effect of decreased mean due to the reduction of parameter uncertainty range was very limited for parameters Qeq and De (HM around 1). Figure 8 presents the 3D plots of the revised variance ratio (HV) for all parameters at t 5 = 5 × 10 4 years (corresponding to the time of maximal mean flux).
From the plots in Figure 8, the effect of changing the parameter uncertainty range on the flux variance is observable. The reduction of uncertainty of Tlarge values to a certain range would lead to a significant decrease in the flux variance. But the plot also showed the presence of a parameter range that would lead to the increased flux variance.
For ease of interpretation, the change in the parameter value range was explored in 2D plots derived from the 3D plots above. a significant decrease in the flux variance. But the plot also showed the presence of a parameter range that would lead to the increased flux variance. For ease of interpretation, the change in the parameter value range was explored in 2D plots derived from the 3D plots above.

Mean Ratio Function
2D plots of the mean ratio function were derived considering the following aspects:

Mean Ratio Function
2D plots of the mean ratio function were derived considering the following aspects: (a) a plot of HM keeping the lower bounding value of input parameter X i at q 1 = 0 and varying the upper bounding value (from q 2 = 0 to q 2 = 1); (b) a plot of HM keeping the upper bounding value of input parameter X i fixed at q 2 = 1 and varying the lower bounding value (from q 1 = 0 to q 1 = 1). Tlarge range to q2 > 0.65 would lead to an increased mean flux; assigning the lower value (q2 < 0.65) would yield a reduced mean flux at the time being analyzed (t5 = 5 × 10 4 years). The value of Tlarge at q2 = 0.65 for the sorted series of randomly sampled values (size N = 1000) was equal to 4.5 × 10 4 years. Similarly, increasing the lower bound up to value at q1 = ~0.65 while keeping q2 constant and at its original value would result in an increased mean flux (ratio HM > 1). However, an increase to greater values (q1 > 0.65) would lead to a significant reduction of the mean flux (see Figure 9).
The parameters SNF DR and IRF led to the lower mean flux by decreasing the upper bound of their value range. Parameters Qeq and De had a small impact on mean flux at time t5 = 5 × 10 4 years after closure.  The plots in Figure 9 indicated where an effort could be made to improve the model output at t 5 = 5 × 10 4 years after the closure.
Revision and justification of new uncertainty range of parameter Tlarge and decreasing its uncertainty would have a variable effect on the mean flux. Decreasing the upper bounding value of Tlarge range to q 2 > 0.65 would lead to an increased mean flux; assigning the lower value (q 2 < 0.65) would yield a reduced mean flux at the time being analyzed (t 5 = 5 × 10 4 years). The value of Tlarge at q 2 = 0.65 for the sorted series of randomly sampled values (size N = 1000) was equal to 4.5 × 10 4 years. Similarly, increasing the lower bound up to value at q 1 =~0.65 while keeping q 2 constant and at its original value would result in an increased mean flux (ratio HM > 1). However, an increase to greater values (q 1 > 0.65) would lead to a significant reduction of the mean flux (see Figure 9).
The parameters SNF DR and IRF led to the lower mean flux by decreasing the upper bound of their value range. Parameters Qeq and De had a small impact on mean flux at time t 5 = 5 × 10 4 years after closure.

Variance Ratio Function
Similar to the effect on the mean flux, the effect of reducing the input parameter uncertainty on radionuclide flux (in terms of the variance) was estimated. Figure 10 presents the variance ratio functions due to:

Variance Ratio Function
Similar to the effect on the mean flux, the effect of reducing the input parameter uncertainty on radionuclide flux (in terms of the variance) was estimated. Figure 10   As Figure 10 shows, increasing the lower bound of the input parameter Tlarge would result in increased radionuclide flux variance at the analyzed time (5 × 10 4 years). This effect was similar to the estimated impact on the mean flux. Decreasing the upper bound of the parameter to value at q2 = ~0.7 or less would result in the reduction of the flux variance only, while a reduction to value at q2 > 0.7 would lead to slightly increased flux variance. As Figure 10 shows, increasing the lower bound of the input parameter Tlarge would result in increased radionuclide flux variance at the analyzed time (5 × 10 4 years). This effect was similar to the estimated impact on the mean flux. Decreasing the upper bound of the parameter to value at q 2 =~0.7 or less would result in the reduction of the flux variance only, while a reduction to value at q 2 > 0.7 would lead to slightly increased flux variance.
The impact of parameter SNF DR on the I-129 flux variance was similar to the effect on the mean. The increased lower bound of parameter range would lead to the higher flux variance, while the lower upper bound would affect decreased flux variance by a factor of~2.
Justified uncertainty reduction by increasing the lower bound of parameter Qeq had some potential to result in decreased variance, which would be lower in comparison to some parameters, such as Tlarge, SNF DR, and IRF. For IRF, only a lower value of the upper bound would lead to a significant decrease in the flux variance. The effect of uncertainty reduction of parameter De would decrease the flux variance by up to 2 times only if the parameter uncertainty range is reduced significantly (to q 2 < 0.2).

Effect of Parameter Uncertainty Reduction on Mean Flux over Time
The time of the maximal mean flux could change with reduced parameter uncertainty. Thus, by fixing the parameter bound at certain quantile, the possible effect on the mean flux should be explored over the other time points.
For example, while fixing the upper bound of Tlarge at q 2 = 0.6 would cause the mean flux at time point t 5 to be reduced, the mean flux would be higher at earlier time points (t 3 , t 4 ) by a factor of 1.5 (see Figure 11c,d).
As Figure 11 demonstrates, reduced mean radionuclide flux could be observed in case of decreased upper bound of parameters SNF DR (ratio HM is less than 1) especially at a longer time scale (after 1 × 10 5 years after closure). The figure also showed that for parameter SNF DR, decreasing the upper bound to the parameter value at, for example, q 2 = 0.8 (for the current sorted set of sampled values (N = 1000) it corresponded to 5.7 × 10 −7 1/year) would lead to a decrease in the mean flux by at least a factor of 1.14. The impact of parameter SNF DR on the I-129 flux variance was similar to the effect on the mean. The increased lower bound of parameter range would lead to the higher flux variance, while the lower upper bound would affect decreased flux variance by a factor of ~2.
Justified uncertainty reduction by increasing the lower bound of parameter Qeq had some potential to result in decreased variance, which would be lower in comparison to some parameters, such as Tlarge, SNF DR, and IRF. For IRF, only a lower value of the upper bound would lead to a significant decrease in the flux variance. The effect of uncertainty reduction of parameter De would decrease the flux variance by up to 2 times only if the parameter uncertainty range is reduced significantly (to q2 < 0.2).

Effect of Parameter Uncertainty Reduction on Mean Flux over Time
The time of the maximal mean flux could change with reduced parameter uncertainty. Thus, by fixing the parameter bound at certain quantile, the possible effect on the mean flux should be explored over the other time points.
For example, while fixing the upper bound of Tlarge at q2 = 0.6 would cause the mean flux at time point t5 to be reduced, the mean flux would be higher at earlier time points (t3, t4) by a factor of ~1.5 (see Figure 11c,d).
As Figure 11 demonstrates, reduced mean radionuclide flux could be observed in case of decreased upper bound of parameters SNF DR (ratio HM is less than 1) especially at a longer time scale (after 1 × 10 5 years after closure). The figure also showed that for parameter SNF DR, decreasing the upper bound to the parameter value at, for example, q2 = 0.8 (for the current sorted set of sampled values (N = 1000) it corresponded to 5.7 × 10 −7 1/year) would lead to a decrease in the mean flux by at least a factor of 1.14.
Reducing of the upper bound of IRF (at earlier times t1-t5) and De (soon after the start of radionuclide release) would also affect decreased mean flux as the HM ratio was less than 1.
On the other hand, increasing the lower bound of parameters SNF DR and IRF would lead to the increased mean as the ratio HM > 1 (see Figure 12). Meanwhile, increasing the lower bound of parameter Tlarge to q1 = 0.65 would result in a lower mean flux at time points t2, t3, t4 and almost no change at time point t5 ( Figure 11). If such reduction cannot be justified, then changing q1 to less than 0.65 would be expected to result in a lower mean at earlier times, but this would lead to an increased mean to some extent at time t5. For the current set of sampled values (N = 1000), q1 = 0.65 for parameter Tlarge corresponded to 4.6 × 10 4 years.
As previously mentioned, increasing the lower bound of parameters SNF DR and IRF would lead to an effect on the increased mean as the ratio HM > 1 (Figure 12). Reducing of the upper bound of IRF (at earlier times t 1 -t 5 ) and D e (soon after the start of radionuclide release) would also affect decreased mean flux as the HM ratio was less than 1.
On the other hand, increasing the lower bound of parameters SNF DR and IRF would lead to the increased mean as the ratio HM > 1 (see Figure 12).
Meanwhile, increasing the lower bound of parameter Tlarge to q 1 = 0.65 would result in a lower mean flux at time points t 2 , t 3 , t 4 and almost no change at time point t 5 (Figure 11). If such reduction cannot be justified, then changing q 1 to less than 0.65 would be expected to result in a lower mean at earlier times, but this would lead to an increased mean to some extent at time t 5

Summary of Discussion
It is clear that for the time-dependent variable of interest (model output), the manner of reducing input parameter uncertainty should only be done after an analysis of HM (HV) ratio evolution over time.
As presented, the performed analysis showed that the parameter ranking for I-129 time-dependent flux differed for different time points after repository closure. In the majority of cases, the ranking based on CSM and CSV methods was the same: soon after radionuclide release, the most significant parameter to contribute to mean flux and variance was the effective diffusivity of I-129 in bentonite (De); later, the defect size enlargement time (Tlarge) became dominant (up to 10 5 years after the closure). At t = 1 × 10 5 years after closure, the most significant contribution to the mean flux came from SNF dissolution rate (SNF DR), and the most significant contributor to the flux variance came from Tlarge. For longer periods (up to 1 million years after closure), the SNF DR was the main contributor to the mean flux and its variance.
The smallest and the largest values of SNF DR were the most significant contributors to the flux variance, which was not observed from the CSM plot.
Importance De (soon after the start of release) and SNF DR (in a long-term perspective) correlated with the regression analysis results; however, the greater importance of the defect size enlargement time was observed for the earlier time points (before 1 × 10 5 years after the closure).
The analysis also showed that the effect of parameter uncertainty range reduction on the mean flux and its variance could be quantified without additional model runs by using the mean and variance ratio functions.
This study focused on the analysis of I-129 migration through engineered barriers. The observations were nuclide specific, and also depended on repository evolution conditions; therefore, in the case of other radionuclides or other disposal conditions, the analysis and parameter ranking exercise should be performed.
As CSM and CSV plots can be developed using the same set of sampled parameter values and do not require any additional model runs that are already available from probabilistic analyses, these tools provide additional valuable information for further research and have the potential to be applied more widely in the context of radioactive waste disposal. As previously mentioned, increasing the lower bound of parameters SNF DR and IRF would lead to an effect on the increased mean as the ratio HM > 1 (Figure 12).

Summary of Discussion
It is clear that for the time-dependent variable of interest (model output), the manner of reducing input parameter uncertainty should only be done after an analysis of HM (HV) ratio evolution over time.
As presented, the performed analysis showed that the parameter ranking for I-129 time-dependent flux differed for different time points after repository closure. In the majority of cases, the ranking based on CSM and CSV methods was the same: soon after radionuclide release, the most significant parameter to contribute to mean flux and variance was the effective diffusivity of I-129 in bentonite (De); later, the defect size enlargement time (Tlarge) became dominant (up to 10 5 years after the closure). At t = 1 × 10 5 years after closure, the most significant contribution to the mean flux came from SNF dissolution rate (SNF DR), and the most significant contributor to the flux variance came from Tlarge. For longer periods (up to 1 million years after closure), the SNF DR was the main contributor to the mean flux and its variance.
The smallest and the largest values of SNF DR were the most significant contributors to the flux variance, which was not observed from the CSM plot.
Importance De (soon after the start of release) and SNF DR (in a long-term perspective) correlated with the regression analysis results; however, the greater importance of the defect size enlargement time was observed for the earlier time points (before 1 × 10 5 years after the closure).
The analysis also showed that the effect of parameter uncertainty range reduction on the mean flux and its variance could be quantified without additional model runs by using the mean and variance ratio functions.
This study focused on the analysis of I-129 migration through engineered barriers. The observations were nuclide specific, and also depended on repository evolution conditions; therefore, in the case of other radionuclides or other disposal conditions, the analysis and parameter ranking exercise should be performed.
As CSM and CSV plots can be developed using the same set of sampled parameter values and do not require any additional model runs that are already available from probabilistic analyses, these tools provide additional valuable information for further research and have the potential to be applied more widely in the context of radioactive waste disposal.

Conclusions
With uncertainty and sensitivity analysis at a low value of determination coefficient of regression analysis, additional conclusions on parameter importance ranking have been drawn in case of I-129 flux from EBS of a conceptual disposal facility for RBMK-1500 SNF: 1.
The CSM identified defect size enlargement time on the I-129 time-dependent flux to have greater importance relative to the effective diffusivity in bentonite and instant release fraction of I-129 (identified in regression analysis) at earlier time points (for a period of 5 × 10 3 -5 × 10 4 years after the repository closure). This importance ranking overcame the results from the regression analysis.

2.
The importance of defect size enlargement time was confirmed with the CSV method; its largest contribution to the variance of I-129 flux into the geosphere occurred at a time from 5 × 10 3 -10 4 years after repository closure.

3.
Soon after the onset of radionuclide release, the most significant contributing parameter to the mean flux and its variance was effective diffusivity of I-129 in bentonite. For longer periods (up to 1 Million years after repository closure), the SNF dissolution rate was observed as the most significant contributor to the mean flux and its variance. These observations were in line with the results of the regression analysis. 4.
At t = 5 × 10 4 years after repository closure (time of maximal mean flux), the most significant contributor to the mean flux was the SNF dissolution rate; however, the most significant contributor to flux variance was the defect size enlargement time.

5.
The effect of decreased mean flux before t = 5 × 10 4 years after repository closure was observed from reduced defect size enlargement time uncertainty: if the lowest bounding value of the defect size enlargement time is increased to be 4.6 × 10 4 years (parameter value at q 1 = 0.65), then the mean flux at t = 5 × 10 3 -2 × 10 4 years would be lower. Almost no effect on the mean flux would occur at t = 5 × 10 4 years after repository closure. If such input parameter range is not justified, then the increase to q 1 < 0.65 would lead to a lower mean at earlier time points and an increased mean (to some extent) at t = 5 × 10 4 years after repository closure. 6.
If justification of the upper bounding value of the defect size enlargement time would lead to a value less than 4.6 × 10 4 years (parameter value at q 2 less than 0.65), this would result in a decreased mean flux at t = 5 × 10 4 years after the repository closure would be expected; however, this would cause increase in the mean flux at earlier time points. Fixing the upper bounding value of Tlarge to a value at q 2 = 0.6 would lead to the lower mean flux at t = 5 × 10 4 years, but at earlier time t = 1 × 10 4 years, t = 2 × 10 4 years, the mean of flux would be greater by a factor of~1.5. 7.
The effect of decreased mean flux in the long-term (up to 1 million years after closure) was observed in case of a justified reduction of the upper bound of the SNF dissolution rate only. Reducing parameter SNF DR uncertainty range from [10 −8 , 10 −6 ] 1/year to [10 −8 , 5.7 × 10 −7 ] 1/year (q 2 = 0.8) would lead to a decrease of mean flux by at least a factor of 1.14. 8.
CSM, CSV plots, and derived mean (variance) ratios have the potential to be applied more widely in the context of radioactive waste disposal as a means of complementing regression-based sensitivity analyses.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.