Effective Safety Assessment of Aged Concrete Gravity Dam based on the Reliability Index in a Seismically Induced Site

: The seismic performance of the aged-concrete gravity dam (aged-CGD) by safety assessment based on the reliability-index is the main focal point of this study. Determination of reliability-index has been handled by the site seismic hazard analysis with the help of PSHRisk-tool (developed by the authors) and risk assessment. Incorporated with the uncertainties, the failure probability has been carried out by the IDA and fragility analysis. For the nonlinear ﬁnite element model of the CGD, the concrete damage plasticity (CDP) model is adopted. To investigate the aging effect, the hygro-chemo-mechanical model has been taken for different years consideration. Through the failure risk assessment of serviceability and safety level, the target reliability index has been determined here for an existing CGD in Korea. Despite several types of research on the CGD safety assessment, the main novelty of this proposed approach will help the dam operator to check the safety barrier for the aged-CGD. A safety index is investigated by comparing the target reliability-index of the age-CGD with the reliability-index for two potential earthquake levels. However, the proposed approach can implement to check the safety range of any seismic site for any set of earthquakes.


Introduction
Safety evaluation is intimately involved with the word reliability. For the structures, reliability analysis depicts the probabilistic failure consequences by several natural hazardous conditions. Therefore, the reliability-based safety index will help someone understand the safety determination of a structure taking the probabilistic meaning of the traditional concept. The definition of reliability analysis is implied through a function containing all random variables and their stochastic distribution [1]. Structural risk analysis is the expansion of the reliability analysis, which leads to the failure function as a structural response. The failure probability comes from the fragility analysis of the structure, where the vulnerability is defined as the combination of the failure function and the seismic hazard analysis of the site. Consequently, this chapter will briefly discuss the background which will express the related literature review and the scope of this study.

Background and Scope
The approach conducts here the determination of safety assessment of aged-concrete gravity dam (aged-CGD) by a key term reliability index, which is usually expressed either by Damage Probability Matrixes (DPM) or probabilistic failure consequences. Mohammadi et al. [1] illustrated the relationship between the structural redundancy and reliability index of RC frames where the ultimate result shows a result increase of reliability index with increasing the structural redundancy. In this study for determining

Theoretical Background
This section will explain the methodology related to express the proposed approach for determining a reliability-based safety index of the aged-concrete gravity dam in a seismically hazardous site. The proposed approach will show an investigated methodology adopting a case study on Bohyeonsan CGD in the Korean Peninsula. Figure 1 presents the pictorial view of the steps and also the required methodology area unit delineated as follows:

Classical Reliability Assessment
The classical term of "reliability" deals with the complement of the probability of failure (= 1 − P f ) for the given uncertainties of the system [11], but more properly, it is the probability of safety of the structure over a given period of time [12]. The limit-state function G(X) (sometimes also referred to as performance function) can be expressed by two domains such as the safe portion (D s ) and failure portion (D f ) [13]: The time-independent (Figure 2a) failure probability is described by a combined probability density function (PDF) counting the random vector of the safe variables as: 3 of 23 where F R and F S represents the Cumulative Density Functions (CDFs), and the randomness of R and S is expressed by f R and f S [12].  With the consideration of the random function of time, the probabilistic analysis of structural capacity will change as shown in Figure 2b. The structural capacity decreases with time due to the material degradation, while the reliability index for the safety check will also be decreased. A time-varying structural load can be idealized as an order of randomly occurring pulses like S 1 , S 2, . . . S n at times t 1 , t 2 ,...t n with τ (small duration when the degradation due to load is zero) as shown in Figure 2b. Jeppsson [14] explained the Appl. Sci. 2021, 11,1987 4 of 23 several methods of the time-integrated and discrete approaches for the reliability analysis, where the time-variant reliability analysis is defined as the following equation: (4) where F R,t and f S,t are the impromptu CDF of R and the impromptu PDF of S at time t respectively explained by Li et al. [15].
To determine the failure probability P f , several methods have been discussed previously [16][17][18], where three strategies are available to evaluate it such as approximation, simulation and adaptive surrogate-modeling-based methods. The FORM, i.e., first-order reliability method and SORM, i.e., second-order reliability method which is the type of approximation method and the Monte Carlo simulation (MCS), Importance Sampling (IS) [19], Latin Hyperbolic Sampling (LHS) [20] and Subset Simulation (SS) [21] methods are the type of simulation method. Alongside, Adaptive Kriging Monte Carlo Sampling (AK-MCS) can be discussed under adaptive surrogate-modeling-based methods that are based on building a Kriging (aka Gaussian process regression) surrogate model from a small initial sampling of the input vector X available in Marelli, Schöbi and Sudret [13].

Intensity and Damage Measure
For the engineering structures, the IDA-based fragility analysis can be expressed by the intensity measure (IM) as cumulative absolute velocity (CAV), cumulative absolute displacement (CAD), arias intensity (AI), spectral acceleration (S a ), spectral displacement (S d ), peak ground acceleration (PGA), etc. In this study, the PGA has been being used to express the IM in IDA. PGA is used today in most cases because of its connection with the inertial forces and for the specific types of structure like concrete gravity dam (stiff structure) to evaluate the maximum dynamic force [22]. As IDA is involved in performing a series of dynamic analyses of a structural model under a set of GMs, each GM is scaled up to a range of IM. The EQs are scaled to such a level of force to excite the structure through the IDA method so that the structure experiences an ultimate collapse after showing the elastic to inelastic behavior. The structural response can be presented by the IDA with scaler IMs and damage measure (DM) used as an engineering demand parameter (EDP). Generally, the maximum inter-story drift, the separate peak story drift, peak floor acceleration and damage index are used as EDP. The relative displacement of the CGD is used in several types of research previously for seismic assessment such as Nahar et al. [23], Sen [24] and Tekie and Ellingwood [25]. Therefore, this study also has been followed the crest displacement as an EDP for the IDA curve, from where the fragility analysis has been carried out. For the seismic design of the hydraulic structure, the specification combining with the international and Chinese two level of intensity approach has been adopted here as (1) a maximum design earthquake (MDE) and (2) a maximum credible earthquake (MCE). The limit states are taken as the serviceability limit state (with a slope of 0.8Ke in IDA) and the safety limit state (with a slope of 0.2Ke in IDA), which is determined here from the IDA curve [26,27]. The Ke is the slope of each IDA curve that is determined according to the two different levels of structural response.

Fragility Formulation
The seismically exciting structure is subjected to damage to such an extent as during ground shaking; the seismic analysis can be carried out through fragility analysis. With a lot of computational efforts, the fragility curve shows a relationship between the probability Appl. Sci. 2021, 11,1987 5 of 23 of exceeding taken limits of EDPs and the exceedance of damage [28]. Mathematically the fragility function can be expressed as follows: where at I M = x, P denotes the GM probability which motivates the structural collapse, ϕ( ) is the function of standard normal cumulative distribution, θ is the median and β is the standard deviation of the IM and lastly, C is the damage state [8]. The fragility parameter followed in this study is the maximum likelihood estimation (MLE) [8,29], where the determination of the parameters can be made by the following: where the number of IM intensities are notified by m, p = 1 (when the case exceeds the limit state (LS)) or p = 0 (not exceed the LS) and q = 1 − p. The fragility curve is exposed from elastic to the collapse level for a set of GM by IDA. Besides, the higher threshold value of the IM is determined by adopting the high confidence of low probability of failure (HCLPF) methodology. HCLPF is the identification of that GM, where the 5% failure confidence remains at 95%.

Seismic Hazard Analysis
The IDA-based fragility curves for different performance levels need to combine with the site seismic hazard analysis to determine risk assessment for any structural system [30,31]. Probabilistic seismic hazard analysis (PSHA) is a well-known basic theory to determine probabilistically the earthquake shaking for a seismically induced site. Presented in Kramer [32], Baker [33] and Pailoplee and Palasri [34], accounting for all the aleatory uncertainties, which will cover the surrounding i earthquake sources from the generation point and all moment magnitudes, PSHA expresses the exceedance rate [34] for a given IM as follows: where λ(M i > m min ) is the probability of earthquakes happening which is greater than m min , and λ(I M > x) is the rate of I M > x. f M i m j and f R i (r k ) are the probability density functions for the earthquake size and the source to site distance, respectively. P(I M > x m j r k ) is the probability of I M > x for an earthquake of size m occurred at source-to-site distance r. λ is the inverse of the return period T R , from where the structural yearly failure probability of any specific earthquake level can be determined. The PDF of a taken seismic source magnitude-frequency follows the Gutenberg-Richter relationship [35] as expressed by the following: where the annual rate of exceedance of m is denoted by λ m , a and b values are known as the Gutenberg-Richter recurrence parameters, β = 2.303b and α = 2.303a. All process in the hazard analysis has been carried out by the own developed software "PSHRisk-tool" (available in https://www.kim2kie.com/3_ach/PSHRisk-tool/PSHRisk_ tool.php, accessed on 25 October 2020) [36] has been used, where can be found the detailed about the theory and application.

Reliability Index
This section focuses on the safety check that has been adopted by the reliability index. The concept of the reliability index is adjuvant determinative to compute the failure probability [37]. The failure condition of the structural measurement can be defined by the failure probability in its design period [1]. The reliability index for the respective failure probability can be determined as follows: To calculate the failure probability for reliability analysis, the measurement of risk can help to fulfill the destination. According to the report presented by ICOLD and CIGB [38], generally, the risk calculation is made by the following expression: Figure 3 shows the overall process to calculate the reliability index with the help of Equations (5), (8), (10) and (12). The determination of the failure probability in risk is also similar, which was first applied in Zion power plant. This is called the Zion method, and it can be estimated by using the two procedures. One is estimated by the product of fragility and differential of seismic hazard curve slope, and the other one is determined by the product of seismic hazard and differential of fragility curve slope. Here, the former one has been used to determine the reliability index as shown in Figure 3.

Dam Specification
For a numerical example, the Bohyeonsan concrete gravity dam has been taken, where Figure 4a shows the real view and Figure 4b shows the sensor locations of the CGD. The details of this CGD are available in Nahar et al. [23], Cao et al. [39] and Rahman et al. [40]. This Dam is located in the upper stream of Gohyeoncheon, which is the second tributary of the Kumho River in South Korea. During the Pohang earthquake, the base acceleration and the system output response were recorded using the two sensors located at the bottom and top of the dam, respectively, and another sensor was used at the surrounding of the dam on the ground to collect ground motion (free filed) shown in Figure 4b. The earthquake epicenter was at 36.12 • N and 129.36 • E, 37 km from the Bohyeonsan dam. Table 1 shows the detailed physical feature of the Bohyeonsan CGD, which is taken from the k-water academy situated in South Korea.

FEM Configuration
The 3D finite element model (FEM) is configured with accounting for the topographic and geological characteristics of the real CGD. The dimensional FEM is modeled with the help of software ABAQUS considering fluid-soil-structure interaction (FSSI). The fluid (reservoir) and foundation (rock) is extended up to 1.5 times [24] of the dam height along each side. The dam body is modeled as a plastic 3D solid element, and the foundation is a modeled homogeneous, elastic and isotropic 3D solid element. The elements for both dam and foundation domains are acted here as the quadratic tetrahedral elements C3D10. Zero pressure is applied at the headless boundary and at the top of the reservoir to consider the damping effect arising from the propagation of pressure waves. Zero displacements are imposed on horizontal translation degrees of freedom at the bellow of the dam and reservoir. Figure 5a shows all dimensions of the dam section, and Figure 5b shows the numerical FEM meshing of the fluid-foundation-dam interaction system. The mesh numbers counted for the FEM in the dam's inside, outside material, foundation, and water were 6776, 41575, 14,012 and 37,260, respectively. Besides showing the whole model, Figure 5b also indicates the axis direction, where the x-axis indicates the direction parallel with the fluid (water) flow, the y-axis denotes the direction perpendicular to the fluid flow, and the z-axis was upward vertically. To illustrate the dam condition for the construction period, the material property is taken from Table 1 following the dimensions and meshing techniques as shown in Figure 5b. After the FEM optimization through RSM (detained in Section 4), the material properties will show the present condition of the CGD. The optimization has been carried out by modal verification with the FEM result and the record data. The FEM was excited by the free field data of the prescribed earthquake and applied to the CGD foundation. Neglecting the less intensive earthquake component, only the horizontal was considered except the vertical component of GM. Assigning the static condition gravity load for the self-weight, the implicit integration was taken along with the non-linear dynamic analysis. The Rayleigh method [41] is adopted here for considering the damping matrix, where 5% damping ratio has been adopted for both the inside and outside dam section (more detailed in Nahar et al. [23]). The epistemic uncertainties as shown in Figure 6 such as aging effect, concrete damage plasticity (CDP) model, parametric randomization and water height effect have been incorporated in this study, where the water has been modeled considering an isotropic 3D solid element as an acoustic model with element type quadratic tetrahedral AC3D10. The aging effect and the CDP model have been discussed in detail in the next. The formulation of the FEM system in a two-dimensional discretized dynamic equation is followed here as explained in Løkke [42] and Sen [24].
The material property was assigned in FEM as the concrete damage plasticity (CDP) model. The damage condition in both tension and compression can be evaluated by the inelastic potential behavior namely CDP. By including the loading combinations both static and dynamic [43], CDP containing four constitutive parameters [44] to define the yield surface in ABAQUS. Here, for assigning these parameters, the data from Nahar et al. [23] were used. The illustration of the tension and compression behavior can be expressed either by the stress-strain curve for compression and also in crack length investigation in tension softening. The stress-strain curve is drawn by two adopted expressions as followed by EN1992-1-1 up to nominal ultimate strain (Figure 7a), and the rest can be defined by the sinusoidal expression (Figure 7b). According to EN1992-1-1, the mathematical expression up to the ultimate strain for the relationship between the compressive stress, σ c , and shortening strain, ε c , is determined by the following equation: where the concrete compressive stress is expressed by , the concrete compressive strain is denoted by ε c , ε c1 is the peak stress compressive strain in the concrete is defined by f c and k = 1.05E cm | ε c1 | f cm . The nominal ultimate strain ε cu1 = 0.0035 were adopted from EN1992-1-1 for all aged concrete (0-75years). Therefore, the extension of the stress-strain curve was made with the sinusoidal part between points using the following equation: Pavlović et al. [45] well explained the parameters for the calculation of this sinusoidal part. The stress-strain behavior of our selected CGD is determined using the above expressions and pictured in Figure 8a for the four (0, 25, 50 and 75 years) different aging conditions. Figure 8b shows the tension stiffening behavior considered the crack opening ratio by adopting the method of concrete exponential tension softening model [46]. The maximum tensile crack starts from σ t , which follows the equation explained in Mirza et al. [47] and is also used in Nahar et al. [23]. The equation for defining the tensile softening behavior is expressed by the following exponential function: where the displacement for crack opening is denoted by w, when the stress cannot be transferred, then, the notation is defined by w c , and c 1 and c 2 are material constants for normal concrete. Figure 8 shows that by increasing the age, the compression and tensile performance of the concrete are decreased. As the hygro-chemo-mechanical model affects the modulus of elasticity with time, the tensile strength is reduced, and this tensile stress for each year has been adopted for the limit state of this study.

Aged-CGD Modeling
The durability of aged concrete is affected because of the joint mechanisms of the external load and environmental incidents such as chemical extension reaction, thermal effect due to heat transfer, freeze-thaw actions, etc. [48]. These affect the concrete strength by making the porosity in it, which is phenomena by hygro-chemo-mechanical cause induced porosity. Gogoi and Maity [49] and Gogoi and Maity [50] used the damage index considering this phenomenon to determine the degraded concrete strength for modeling the aged-CGD. Keeping solidarity with that concept, this study applied the hygro-chemomechanical effect on the existing Bohyeonsan CGD in Korea. As explained in Gogoi and Maity [49] and Gogoi and Maity [50], the hygro-chemo-mechanical affect the concrete modulus of elasticity. By adopting that parameter, the same as in Nahar et al. [23], the damaged modulus of elasticity (E e ) for the prescribed CGD is compared with the sound modulus of elasticity (E 0 ). Figure 9a shows the graphical representation of the contact surface of the dam and water, from which, in Figure 9b, it is shown that the modulus of elasticity for the CGD is reduced with age in both inside and outside materials.

Optimization of FEM
The optimization of the FEM model has been carried out by the RSM to know the present condition of the Bohyeonsan CGD. This will give the correct estimation of the degraded modulus of elasticity for showing the aging effect on safety evaluation. RSM is a statistical and mathematical technique for system identification initially developed and explained by Box and Wilson [51]. In generating an empirical model building, RSM plays an important role [52]. The objective of this method is to optimize the output (y) variables generally affected by the input variables (x), which can be expressed by Equation (17) [53].
where the responses of the structure are indicated by f (x 1 , x 2 , . . . x k ), and error of that response is described by ε. In order to skip the unnecessary true function, for many reasons, a flexible mark off polynomial function can be adopted. The two most common forms of a polynomial function are as follows: The coefficients β 0 , β 1 · · · · · · are parameters that need to be estimated from the data. Several mathematical tools are available for the DoE, and CCD developed by Box and Wilson [51] is a very useful tool to reduce the experimental phase. This equation is also used recently for evaluating the seismic capacity and risk prediction of CGD in Nahar et al. [23]; Cao et al. [39] and Cao et al. [54]. The CCD is classified into three types according to the axial points (an indication of the curvature) as central composite circumscribed, central composite inscribed and central composite face-centered. In this study, using CCD, the total number of samples for running the experiment is calculated by the following Equation (20).
where p is the factor number, and n c denotes the center-point number of the sample, which is generally considered as 5 or 6. Using Equation (20) the number of numerical simulations is calculated as nine with two input variables such as CoE (coefficient of modulus of elasticity) and ρ (mass density). To get the optimized modulus of elasticity, the CoE will help by multiplying with the present variables. The recorded acceleration data was collected for the Pohang Earthquake, which was the strongest earthquake in the Korean Peninsula [55]. The simulation had been done with the nine combinations of two material properties (input variables), and after optimization by RSM, the present value of those two parameters has been come out as 0.787 and 2.32 (tone/m 3 ), respectively. In addition, the final inside and outside modulus of elasticity is 10,835 MPa and 13,269 MPa, respectively. The response spectra comparison between the before and after optimization has been carried out and showed in Figure 10. Figure 10 indicates that the similarities of the top RS after optimization is acceptable for the structural system identification. The FEM with the present parameters then verified by comparing the fundamental frequencies by the frequency domain decomposition (FDD) of the experimental data. According to Brincker et al. [56] the FDD can decompose the system response from recorded data by following a simple technique that is related to the decomposition of each of the estimated spectral density matrices explained in Ko et al. [57].
The FDD from the sensor database is visualized in Figure 11a, where the first peak depicts the fundamental frequency of the CGD. Analysis of the mode shape of 3D FEM considering the optimized material properties is also illustrated in Figure 11b. From the modal analysis, the most dominating mode shape shows the fundamental frequency of the FEM. From Figure 11, the coefficient of variance (CoV) between the FDD result and FEM result is almost the same. The acceptable value for the CoV should be less than 15% [58] for the FEM verification, where this study shows that the CoV is 2%. Therefore, it can be said that the above FEM for the seismic safety assessment of the aged-CGD is properly validated and verified.

Ground Motion Selection
Non-linear analysis of the structure required the right selection of GM, which can provide the real scenario during an earthquake. This selection depends on several uncertainties that need to be characterized as aleatory uncertainties by some necessary steps. The steps elapse to determine the most vulnerable seismic contribution to the 2% of 50 years shaking through site seismic model source-to-site distribution, uncertainty and PDF of earthquake magnitude, seismic hazard analysis and deaggregation plot. Then a set of GM can be selected either from PEER or artificially generated [59] and scaled with the target spectra [60,61]. According to Choun et al. [62], the seismic record of Korea is divided into four seismic source models. For selecting the most effective earthquake range in the CGD site, seismic hazard analysis has been carried out by "PSHRisk-tool" [36]; Figure 12 shows the deaggregation plot for 2% of 50 years shaking PGA. From Figure 12, the peaks correspond to magnitude-distance contributions that contribute more to the hazard than has been considered. In correspondence with the analysis, around a range of moment magnitude (5.0 to 6.5) and the closet ruptured distance ranges (50 km to 204.7 km) is adopted to select the GMs for reducing the computational efforts as shown in Figure 13a. The uncertainties of the GM are reduced by matching with the taken design response spectra (DRS) (RG 1.60) [63,64]. By avoiding the vertical components because of the lesser acceleration, Figure 13a shows the scatter plot and Figure 13b shows the matched 30 GMs with the DRS. The scatter plot represents a specific PGA of GM distribution through the required earthquake magnitude and rupture range.

Seismic Failure Assessment
This section explains the seismic analysis by the validated and verified FEM, which is subjected to 30 selected GMs. As explained in the previous section, the FEM was only subjected to the horizontal two-directional GM components. One was applied to the x-axis, and another one was applied to the y-axis. From the modal analysis, the most vulnerable direction has been notified in the x-directional response, which leads to the result in evaluation and analysis in this manuscript to the x-direction only.
The analysis was also performed considering four different years for showing the response of the aging effect during an earthquake. By adopting the calibrated FEM through system identification, at the initial stage (0 years), the young's modulus of elasticity is considered as non-damaged E 0 [23]. The other cases are adopted by taking aging effects as 25, 50, and 75 years, respectively. The model was at first subjected to the gravity load in all cases and simultaneously the non-linear time history analysis had been performed at the foundation base and foundation side. As for CGD one of the most key limit states is the relative displacement against the fluid flow, the crest displacements acted here as the DM in IDA. Following the IM-DM points in the IDA from Jalayer and Cornell [65], in each increment of GM, the mean and standard deviation (SD) were calculated to show the 50% fractile curve (mean), 16% fractile curve (mean-SD) and 84% fractile curve (mean+SD) were calculated. Figure 14 shows the IDA curves for the 30 earthquakes for each specified year along with the mean, curves of the dam crest serviceability. The dynamic analysis is allowed here in IDA from 0.1g to 1g for each earthquake scaling. To identify the two-limit states from the IDA, the mean, 16% and 84% fractile curves of the dam crest are visualized in Figure 15. The serviceability and safety points have been calculated from the slope of the curve as discussed in Section 2.2.1.  The seismicity parameters for the failure probability have been tabulated for two performance levels in each year in Table 2. For the mean IDA, the crest displacement in the serviceability stage is observed from 0.259-0.159 (g) by increasing the age of the CGD. Besides the mean, the 16% and 84% limits these values from 0.269-0.175 (g) and 0.227-0.180 (g), respectively. Therefore, for serviceability level, these earthquakes will not obviously infect the dam. The controlling IM for the 0-year aged-CGD is from 0.227-0.269, corresponding which the displacement limit state for the fragility analysis. Similarly, for safety level, the controlling IM for 0 years aged-CGD is identified from 0.405-0.438. The case study of the considered dam is then seismically analyzed by these limit values by considering the key point as the crest displacement. The 3D of the existing multipurpose dam has been analyzed probabilistically through the time history, nonlinear analyses considering the effect of the hygro-chemo-mechanical effect on the structural vulnerability analysis. Crest displacement was carried out from the IDA curve to illustrate the fragility curves. By using Equation (5), fragility curves have been shown in Figure 16 to summarize the failure pattern for the mean, 16 and 84% fractile in serviceability and safety level. In case of a maximum credible earthquake, from any year consideration in Figure 16, the serviceability and safety show almost 1. If it is considered for the maximum credible earthquake, these serviceability points are observed as 0.90-0.95, wherein for safety point, it is almost also 1. From the analysis of the IDA, the corresponding fractile PGA for considering the crest displacement shows here the safe range in either the MDE or the MCE level. In the global case, it may cross the design earthquake level. To the effect of site seismic hazard uncertainties, this study adopted a further reliability-based safety assessment. Therefore, the seismic hazard analysis will show the site-specific safety evaluation in the next section as well.

Site Seismic Hazard Analysis
The site of the taken dam is located in South Korea, which is situated in a moderate seismic zone. The site-specific magnitude is appropriate for use in seismic risk assessment because there are only a few strong motion data and seismological information for the dam sites in Korea. As a result, according to the Choi et al. [10] for Korea, there are four seismic source models as shown in Figure 17 along with the dam site. More than 2000 earthquake records were used for making the seismic source models [66], and the magnitude-frequency relationship was presumed the same as the Gutenberg-Richter study [35]. The site-specific magnitude range is important to predict potential issues. Which GM intensity will contribute more to the hazard may cause the vulnerability to risk assessment of the CGD. Therefore, all the sources in each seismic model have been used to illustrate the hazard curve as shown in Figure 18a-d. The total hazard curve for each source model can then be extracted using the ground motion prediction equation (GMPE) from Cornell et al. [67], and all the total hazards along with the average curve have been visualized in Figure 18e.  The idealized reason for failure in the aged-CGD was characterized by the function of time in the previous study of the author of this manuscript [23]. The connection between the seismic capacity and the modulus of elasticity was correlated with the analysis of some specifically selected GMs. In this case, the aging effect has been more characterized by the real scenario of the site seismic risk analysis. To determine the safety, the dam reliability index is strongly supported by the mean annual rate of earthquake exceedance. The failure probability of the reliability estimation covers the probable failure pattern by considering all the epistemic and aleatory uncertainties.

Safety Assessment of Aged-CGD
The safety check has been carried out here through the design IM for different earthquake levels. To give a guideline for the dam operator, Operational (OP) and Life Safety (LS) [68] is selected to determine the reliability index for earthquake levels.
In the safety check, the procedure discussed in Section 2.2.4 is applied to determine the target reliability index of the CGD. The fragility curve from the failure probability for different year consideration has been multiplied by the integration of the slope of the site seismic hazard curve. Figure 19a-d shows the reliability index curve for the example CGD for each considered year, showing the time-variant seismic safety evaluation. Figure 19 shows the reliability curve for mean, 16% and 84% fractile for both the serviceability and safety stage. For both levels of the target PGA (serviceability and safety), the HCLPF point from the fragility curve has been determined and indicated as target (serv.) and target (safety) as shown in Figure 19. The target reliability index then is calculated by the identification of the structural capacity. Table 3 shows the list of HCLPF points from the fragility curve for each year.
(safety) as shown in Figure 19. The target reliability index then is calculated by the identification of the structural capacity. Table 3 shows the list of HCLPF points from the fragility curve for each year.  Figure 20, the regression analysis is summarized paired with their annual frequency of exceedance for the 73 years (OP) and 975 years (LS) in hazard levels. To compare with the target value, the reliability index for different earthquake levels has also been taken for two of these performance levels (OP and LS). From the average curve, as shown in Figure 20, the design earthquake PGA is 0.075g and 0.221g, respectively, for the taken CGD site. Using a peak ground acceleration value of 0.302g, which is equal to that estimated for the records having a 2% in 50 years hazard level, it can be extracted directly from the mean fractile of the fragility curve. From Figure 19, the design PGA, as shown in Figure 20, two, 0.075g and 0.221g, were taken to measure the reliability index.  From Figure 20, the regression analysis is summarized paired with their annual frequency of exceedance for the T R = 73 years (OP) and T R = 975 years (LS) in hazard levels. To compare with the target value, the reliability index for different earthquake levels has also been taken for two of these performance levels (OP and LS). From the average curve, as shown in Figure 20, the design earthquake PGA is 0.075 g and 0.221 g, respectively, for the taken CGD site. Using a peak ground acceleration value of 0.302 g, which is equal to that estimated for the records having a 2% in 50 years hazard level, it can be extracted directly from the mean fractile of the fragility curve. From Figure 19, the design PGA, as shown in Figure 20, two, 0.075 g and 0.221 g, were taken to measure the reliability index. To visualize the trend of reliability index of the prescribed aged-CGD target (capacity), response with time can be compared with demand (EQ level) as shown in Figure 21. Figure 21a,b shows the comparison between the capacity and demand reliability index for all fractile cases in serviceability and safety level, respectively. With increasing time this CGD shows a safe index in serviceability level even each year, wherein safety level the CGD looks unsafe after 45 years. Therefore, adopting the proposed methodology described in this study is useful to check the safety index in a specific location. The reliability index for the performance level of serviceability (check for OP earthquake level) and safety (check for LS earthquake level) for each fractile (16%, 50% and 84%) is listed in Table 4. The target reliability index was determined from the HCLPF capacity of the dam and compared with each earthquake level to check the safety assurance. Moreover, to check the safety, a 16% to 84% fractile range has been taken in this study. Table 4 shows that the dam for serviceability performance level is safe for the whole year and for all fractile cases. Nevertheless, for safety performance level, the dam fails to satisfy for 50 and 75 years with 84% fractile cases. This process can make it obvious the need to take a more confidence level in the fragility analysis, on which the reliability index is dependent. More prediction of the earthquake performance levels will also be a great solution to check the safety evaluation.

Conclusions
The study presented a practical approach to checking whether a massive structure is safe under seismic excitation. To do so, the classical reliability methodology has been used to give the safety index of an aged-CGD. To investigate the site-specific hazard effect of the aged-CGD, the PSHA is used extensively in this study, where PSHRisk-tool has been used to determine the hazard curve and deaggregation plot. IDA-based seismic performance assessment is carried out to set the limit state in two-levels. The failure probability for the structural response has been adopted by fragility analysis in each performance level. A proper 3D FEM of required CGD in Korea is reflected in its dynamic responses through the time history analysis. To excite the FEM, a set of 30 EQs was selected in that site for analyzing the maximum contribution of IM from the hazard curve.
Seismic risk evaluation combining the failure probability with the integration of the slope from the hazard curve is considered here for illustrating the reliability index curve. The hygro-chemo-mechanical effect is taken for defining the aging condition. The target reliability index is expressed by two different performance levels, which encompasses damage levels of the CGD. It is found that the GM characteristics showed a great effect on the dispersion of the fragility curves.
Moreover, it is notified at the 50 and 75 aged-CGD, the material property is become deteriorated due to strength reduction. The reliability index curve showed the safety pattern was for the same year, which indicates it is less safe than the other year's consideration. For the two EQs level carried from the previous study, the damage index for the corresponding EQs level has been listed to compare the reliability index with the target one (from the capacity of the structure). This safety evaluation process will give information about the safety level of a CGD through the reliability index check. The safety range will be satisfied if the reliability index for any earthquake level is greater than the target reliability index of the aged-CGD. The life safety has been taken in this study to show the maximum reliability index, where the target reliability index divides the safe and unsafe zone for the CGD for a considered earthquake. The safety measurement of the seismic site in Korea depicts the aged-CGD as safe in all taken earthquake levels but unsafe for some fractile case of the different fractile level.
Lastly, it can be said that the main advantage of this proposed approach is that it provides a practical guideline to the user for investigating the safety range of potential earthquakes. The time-dependent aging consideration is a new aid to the reliability index for predicting the safety check of the deteriorated structure. This method can be implemented in another structure or can be updated suitably with another set of earthquakes.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.