Service Life Design of Concrete Structures Made of High-Volume Limestone Powder Concrete—Case of the Carbonation-Induced Corrosion

: One of the paths to CO 2 emissions reduction in the concrete industry is to use low-clinker cements, providing at the same time the performance of concrete that is adequate for application in concrete structures. This paper explores the impact of the clinker replacement with high amounts of limestone powder (21–70% in the powder phase) on concrete carbonation resistance. To quantify this impact, the empirical relationship between the carbonation resistance and the compressive strength of the high-volume limestone powder concrete (HVLPC) was established. For that purpose, the regression analysis was applied on the experimental results collected from the published research. The service life analysis based on the full probabilistic approach was performed using the fib Model Code 2010 prediction model and proposed empirical relationship. The first-order reliability method (FORM) was applied to solve the limit state function of reinforcement depassivation with a reliability index equal to 1.3. The obtained minimum concrete cover depths were 40–110% higher compared to those prescribed in the current European standard EN 1992-1-1:2004 for indicative strength classes. Based on the full probabilistic analysis, recommended cover depths are given for all carbonation exposure classes, commonly applied concrete strength classes, and service lives of 50 and 100 years.


Introduction
A key challenge of this century is the need to mitigate and adapt to climate change, and avoid negative impacts on the environment, the economy, and social well-being.Across the globe, the negative impacts of climate change are already recorded and set to continue: higher risks from natural hazards, global sea rise, significant species loss, ocean acidification, etc. [1,2].To avoid the worst, it is imperative to "hold the increase in the global average temperature to well below 2 • C above pre-industrial levels and to pursue efforts to limit the temperature increase to 1.5 • C above pre-industrial levels" [3]-the main objectives of the flagship Paris Climate Agreement (COP 21).Achieving results with the urgency needed to meet COP 21 commitments requires targeting high-emitting industries and activities.The construction industry largely contributes to Europe's environmental footprint: 50% of natural raw materials use, 40% of total energy consumption [4] (as the single largest consumer), and 36% of all greenhouse gas (GHG) emissions [5], while construction and demolition waste (CDW) makes 46% of the total waste generated [6].
Concrete and other cement-based materials are responsible for a large share of the construction industry's environmental impacts.Concrete is the most widely used construction material in the world, with 33 billion tons produced annually [7]; it is composed of cement, aggregates, water, mineral, and chemical admixtures.Due to the cement (clinker) production process and the chemical reactions involved, it is responsible for 8-9% of annual anthropogenic CO 2 emissions [8].Therefore, decarburization of the concrete/cement sector plays an important role in achieving a carbon neutral future.
Decreasing the clinker content by replacing it with so-called supplementary cementitious materials (SCMs) with low embodied CO 2 presents one of the viable solutions to this problem.Today we already use commercial cements with part of the clinker replaced with SCMs: CEM II/A with up to 20% fly ash (FA), ground granulated blast furnace slag (GGBS) and limestone (LS), CEM II/B with up to 35% FA, and GGBS, LS, and CEM III with up to 95% GGBS.Bearing in mind the limited global availability of commonly used SCMs (FA and GGBS) [9], alternative SCMs, abundant in nature, should be looked for.Limestone, if finely ground, can serve this purpose [10].However, compared to FA and GGBS, it is less effective and has a higher impact on the mechanical and durability-related concrete properties due to its very low reactivity [11][12][13][14][15]. Therefore, higher clinker replacements with limestone powder (>15-20%) present a problem for the structural concrete at all performance aspects.
One of the important durability-related properties of the concrete is its carbonation resistance.This depends mainly on the porosity and pore size distribution on the one hand, and the amount of carbonatable constituents, i.e., alkaline reserve, on the other.On both levels, the limestone powder concrete is inferior to Portland cement concrete due to the dilution effect of the limestone.Several experimental campaigns reported so far in the literature [11,[16][17][18][19] as well as review papers [20][21][22][23] showed that the replacement of the clinker with up to 15% limestone powder practically does not affect the carbonation resistance, while higher replacement ratios reduce it significantly, the reduction being larger with a larger replacement ratio.However, it was shown in the previous research that it was possible to significantly improve the high-volume limestone powder concrete (HVLPC) performance by optimizing the paste composition, lowering the water-to-cement ratio using a larger superplasticizer content, improving the particle packing, and incorporating the nano-calcium carbonate [22,[24][25][26][27][28][29][30][31].
Carbonation resistance of the concrete in the concrete cover of reinforced concrete (RC) structures is not an only, but it is a decisive factor for the RC structures' resistance to carbonation-induced corrosion, and therefore their durability.The carbonationinduced corrosion of steel reinforcement is one of the major sources of the RC structures' deterioration-according to Jones et al. [32], two thirds of all structural concrete is exposed to this type of deterioration mechanism.
Modern standards in the construction area tend to rely on the performance-based design for all limit states including durability, as is the case for the fib Model Code 2010 [33], for instance.For a performance-based design, it is necessary to have a reliable prediction model for the carbonation-induced reinforcement corrosion.The carbonation process, as mostly a diffusion-controlled process, is usually described with the Fick's first law in the form: where x c is the carbonation depth, k is the carbonation rate coefficient, and t is time.Unlike that for the concrete with FA and GGBS blended cements [22], very few mathematical prediction models for the limestone powder concrete have been developed so far [11,34,35].These models are complex and require knowledge on many parameters that influence the carbonation process.For everyday engineering practice, however, simpler models are needed.One such simpler model is given in the previously mentioned fib Model Code 2010 [33], originating from the fib Model Code for service life design [36].Although simple, this prediction model requires an experimental determination of the carbonation depth under defined conditions.Therefore, the property that reflects the concrete carbonation resistance has to be experimentally determined.Current European standards EN 12390-10:2018 [37] and EN 12390-12:2020 [38] prescribe the test conditions for determining the carbonation resistance under natural and accelerated conditions, respectively.Tests under natural conditions should last for one year, while for the tests under accelerated conditions, which are shorter (70 days), non-standard laboratory equipment-carbonation chamber-is required.To avoid this experimental test, finding an empirical relationship with concrete properties that must be tested in any case, such as concrete compressive strength, would significantly improve the feasibility of the model.For the sake of reliability, this relationship would preferably be established for carbonation resistance under natural conditions.Several authors have already tried to establish such relationships: based on the measured carbonation depths on existing RC structures without information on the exact concrete type [39], or for high volume fly ash and recycled aggregate concrete [40].Gonçalves et al. [41] reported two different power relationships between carbonation resistance and compressive strength: one for concrete with CEM I and CEM II/A (including CEM II/A-L) and the other for concrete with CEM II/B (but excluding CEM II/B-L(LL)), CEM III, CEM IV, and CEM V. Shah et al. [42] found that a meaningful correlation between the carbonation rate coefficient and compressive strength could not be established due to high scattering.This conclusion was drawn on the basis of their own experimental results and results reported in Rathnarajan et al. [43] and Vu et al. [44], which included concrete mixes with various SCMs (GGBS, FA, calcined clay, limestone, and combinations).On the contrary, Lollini et al. [13] concluded that an exponential relationship between the carbonation rate coefficient and compressive strength existed for concretes with different SCMs, although only based on their own test results.They did not propose specific relationships since consensus on this matter was still lacking according to authors.As for high-volume limestone powder concrete (>20% in the powder phase), one proposal was found in the published research, but for accelerated conditions.Marques et al. [45] proposed a linear relationship between the carbonation rate coefficient and compressive strength for CEM II/B-L concrete (up to 35% limestone powder) based on their own test data with a relatively small range of compressive strengths (between 40 and 50 MPa).
Service life modeling based on physico-chemical prediction models and a full probabilistic approach should give estimates on RC structures' durability that are more realistic compared to the deemed-to-satisfy approach mostly based on the empirical evidence.Several researchers have performed a probabilistic carbonation service life analysis but their results were limited to limestone powder concrete mixes tested in their own investigation [13,45,46].Neves et al. [47] and Belgacem et al. [48] established an empirical relationship between the accelerated carbonation rate coefficient and air permeability coefficient and performed a probabilistic service life analysis, giving recommendations for the practical design.However, these recommendations were developed on the basis of their own test results on CEM II/A-L concrete mixes and therefore cannot be generalized.
Therefore, a reliable simplified prediction model for carbonation resistance of HVLPC currently does not exist.Models proposed so far in the literature are based either on test results from one laboratory with a small range of involved parameters or include mixes with various SCMs in lower contents.Consequently, it is not possible to design the HVLPC reinforced structure for a prescribed service life, i.e., it is not possible to calculate the required concrete cover depths.Current standards are valid for some of the cements available in the market and these present a serious obstacle to the application of HVLPC in RC structures.In this paper, a simplified prediction model based on test results from 14 different laboratories available in the literature was proposed.Using this empirical model and full probabilistic service life analysis, the HVLPC concrete cover depths for exposure classes XC1-XC4 were also proposed.

Objectives
The objectives of this work were: • To establish the empirical relationship between the carbonation resistance and the compressive strength of the HVLPC (21-70% limestone powder content in the concrete powder phase).The relationship was obtained using the regression analysis applied to the experimental results collected from the published research.
• To propose the concrete cover depths for the HVLPC, depending on the limestone powder content and carbonation exposure class.The service life analysis based on the full probabilistic approach was performed using the fib Model Code 2010 [33] prediction model and proposed empirical relationship.

Methodology
For describing the carbonation process development with time t, the fib Model Code 2010 [33] prediction model was used: where: x c is the carbonation depth (mm); k is the carbonation rate coefficient (mm/year 0.5 ) defined as: is the execution transfer parameter (-); C s is the CO 2 concentration in the air (kg/m 3 ); W(t) is the weather function (-); R -1 NAC is the inverse effective carbonation resistance of concrete under natural conditions ((mm 2 /years)/(kg/m 3 )).
The inverse effective carbonation resistance R −1 NAC reflects the carbonation resistance of concrete depending mostly on the water-to-cement ratio and binder type, while other parameters in Equations ( 2) and (3) take into account the influence of the environmental and execution conditions.It is recommended in [33,36] ACC , which is experimentally determined under defined accelerated conditions (ACC): where: k t is the regression parameter for the test effect of the ACC test (-); ε t is the error term for inaccuracies that can occur when using the ACC test method ((mm 2 /years)/(kg/m 3 )).
It was assumed that model (Equations ( 2) and ( 3)) is valid for HVLPC as well, except for the inverse effective carbonation resistance that is a material property.In order to establish the empirical relationship between R −1 NAC and compressive strength of the limestone powder concrete, experimental results were collected from the previous research in the 1995-2023 period.Collected test results were then filtered to satisfy the following conditions: limestone powder was the only SCM in the concrete mix with more than 20% participation in the powder phase, whether it was contained in the commercial CEM II/B-L(LL) cement or added to CEM I concrete mix in a certain amount; minimal mean compressive strength of concrete f cm on the cylinder was equal to 20 MPa; measured carbonation depth was not less than 1 mm; minimum natural exposure duration was 140 days; all necessary curing and environmental conditions were reported.Since the curing conditions have a large impact on the carbonation front propagation, a minimal curing period of 7 days was adopted.Applying these filters led to 75 experimental results of carbonation depth measured under natural indoor and sheltered outdoor conditions [12,13,16,18,25,26,[49][50][51][52][53][54][55][56], on which the regression analysis was performed.Those data originate from 14 different laboratories; however, from some experimental campaigns, only 2 or 3 measured data satisfied the mentioned conditions and were included.The data range is shown in Table 1.After establishing a relationship between R −1 NAC and f cm , the service life was predicted using Equation (2).By comparing the probability of the carbonation depth reaching the nominal concrete cover depth (c nom ), the limit state function of reinforcement depassivation can be written in the following form: The first order reliability method (FORM) was applied to calculate the limit function for the defined probability of failure.However, for the calculation of the limit state function, it is necessary to define the distribution parameters for each variable in Equation (5).
Given that concrete cover is a stochastic variable, it can be described with a lognormal distribution, bearing in mind that it cannot have negative values.The nominal concrete cover (c nom ), defined in EN 1992-1-1:2004 [57], represents the mean value of the concrete cover depths (see Table 2).It consists of the minimum cover necessary to protect the reinforcement from corrosion (c min,dur ) increased by an absolute value of the accepted negative deviation (∆c dev ), which depends on the quality of execution works.For cast-insitu structures, the prescribed ∆c dev is 10 mm.It is reasonable to assume that the concrete cover standard deviation (σ c ) is controlled by the execution requirements, which are often interpreted as 5% quantile of concrete cover values [58]: The coefficient k e represents the environmental conditions: where RH real is the environmental relative humidity (%), RH ref is the referent relative humidity (65%), f c is an exponent (5.0), and g c is another exponent (2.5).
Distribution parameters for the environmental relative humidity were adopted according to previous studies [59,60], for each exposure class defined in EN 1992-1-1:2004 [57].The execution transfer parameter (k c ) represents the curing conditions and is defined as: where t c is the curing period (days), and b c is the regression exponent (-).
The curing period (t c ) of 7 days was adopted assuming that this is the standard on-site curing period.It is also the time defined in the fib Model Code for service life design [36] for the accelerated carbonation test.The distribution parameters of the exponent b c were adopted according to the recommendation from the same document [36].In addition, the distribution parameters for the CO 2 concentration (C s ) were adopted according to values defined in [36], taking into account the constant increase in concentration in the past century.For the weather function W(t), the maximum value (1.0) was taken, given that only indoor and sheltered-from-rain samples were analyzed in this study.Finally, the distribution of predicted R −1 NAC was calculated according to the procedure explained in Section 4.2.The summary of all input parameters is shown in Section 4.3.
To solve the limit state function defined in Equation ( 5), it is necessary to define the required minimum reliability index (β).Here, the value of 1.3 (P f ≤ 0.10) was adopted as per the fib Model code for the service life design [36] requirement for the depassivation limit state.A service life of 50 and 100 years was used in the calculations.By solving the limit state function defined in Equation ( 5), it is possible to obtain the required concrete cover depth for the defined parameters, or to determine the service life for the defined concrete cover.

Relationship between R −1
NAC and Compressive Strength of the Limestone Powder Concrete In the analyzed literature, experimental results were reported either as the measured carbonation depth x c or carbonation rate coefficient k; in the latter case, carbonation depth was calculated using Equation (1).
The inverse effective carbonation resistance of concrete, determined under natural conditions (natural carbonation test-NAC) was calculated from Equations ( 2) and (3): using the reported curing and exposure conditions and measured carbonation depths.
Values of R −1 NAC calculated in that way and the reported compressive strengths were considered as mean values.
Concrete mixes were divided into two groups according to limestone percentage: (1) between 21% and 35%, and (2) between 36% and 70% limestone powder in the concrete mix.The limestone percentage was calculated as the ratio between the limestone amount and the powder phase amount (Portland cement + limestone powder).Such division was adopted in order to differentiate between the limestone content available in commercial cement CEM II/B-L(LL) and higher limestone content.
The relationship between R −1 NAC and f cm was best described with the power regression model.However, reducing the power regression to the linear regression model of the log-transformed variables offers a simpler mathematical tool for further application.The linear model of the log-transformed variables is: where b 0 and b 1 are the linear regression coefficients representing intercept and slope, respectively.
The results are shown in Figures 1 and 2 for concrete mixes with between 21% and 35%, and between 36% and 70% limestone powder, respectively.Since the models did not differ substantially, an additional model for the concrete mix with limestone content between 21% and 70% was also estimated and shown in Figure 3.The coefficient of determination R 2 for the linear regression models of log-transformed variables ranges between 0.76 and 0.80 depending on the mix (Table 3).The coefficient estimates, their standard errors, and significance are also given in Table 3.The null hypotheses that b 0 = 0 and b 1 = 0 can both be rejected for all mixes at virtually any significance level.Each regression model was also shown to have independent, homoscedastic, and normally distributed residuals.The normality of the residuals was tested by applying the Kolmogorov-Smirnov goodness-of-fit test at the 5% significance level.The relationship between R −1 NAC and fcm was best described with the power regression model.However, reducing the power regression to the linear regression model of the logtransformed variables offers a simpler mathematical tool for further application.The linear model of the log-transformed variables is: where b0 and b1 are the linear regression coefficients representing intercept and slope, respectively.
The results are shown in Figures 1 and 2 for concrete mixes with between 21% and 35%, and between 36% and 70% limestone powder, respectively.Since the models did not differ substantially, an additional model for the concrete mix with limestone content between 21% and 70% was also estimated and shown in Figure 3.The coefficient of determination R 2 for the linear regression models of log-transformed variables ranges between 0.76 and 0.80 depending on the mix (Table 3).The coefficient estimates, their standard errors, and significance are also given in Table 3.The null hypotheses that b0 = 0 and b1 = 0 can both be rejected for all mixes at virtually any significance level.Each regression model was also shown to have independent, homoscedastic, and normally distributed residuals.The normality of the residuals was tested by applying the Kolmogorov-Smirnov goodness-of-fit test at the 5% significance level.From this analysis, the relationship between the natural inverse carbonation resistance and mean compressive strength of the limestone powder concrete can be proposed in the following way (with strong correlation, R 2 between 0.76 and 0.8): for 21% to 35% limestone powder concrete, including commercial CEM II/B-L(LL) cement, and: for 36% to 70% limestone powder concrete.
As already mentioned, since the relationships for mixes with 21-35% and 36-70% limestone powder are very similar, a single relationship can be proposed for the whole range of 21% to 70% limestone powder, Figure 3: for 21% to 70% limestone powder concrete, including commercial CEM II/B-L(LL) cement, with statistical parameters of the linear regression model given in Table 3.
The carbonation resistance depends mainly on the porosity and amount of carbonatable constituents.As already mentioned, HVLPC has lower carbonation resistance compared to pure-clinker concrete due to the low reactivity of limestone powder on the one hand, and low clinker content on the other.The first problem (porosity) can be resolved with a properly designed mix even with high limestone powder content, i.e., very good compressive strength can be achieved.For that, it is necessary to optimize the HVLPC mix design, which can be completed in several ways.The most effective way is to use a low content of water and in that way keep the water-to-clinker ratio in limits that do not significantly decrease the strength [24][25][26].Since the water content required for adequate workability is usually above the content needed for full clinker hydration, relatively high amounts of plasticizer are needed to keep this excess water at a minimum.Secondly, use of the finely ground limestone improves particle packing density and density of the ITZ (filler effect).Mixing limestone powders of different fineness can also contribute to compressive strength for the same reasons [24].Finally, adding the limestone powder in contents larger than that for the simple replacement of clinker increases the total powder paste volume, and improves filling of the voids and workability at the same time [27].Most of the mixes in the database with limestone powder content higher than 35% were designed in that way and achieved compressive strength higher than 30 MPa, even at clinker contents as low as 150-180 kg/m 3 .The second problem, a low amount of carbonatable constituents, depends mainly on the clinker content: the lower the clinker content, the lower the alkaline reserve.Due to low Ca(OH) 2 content, the carbonation rate is faster, i.e., carbonation resistance is lower and the mix optimization methods mentioned here cannot help.For that reason, HVLPC has lower carbonation resistance compared to pure-clinker concrete for the same compressive strength.The obtained results are: the similar relationships (11) and ( 12) for 21-35% and 36-70% limestone powder groups, respectively, mean that the impact of porosity, which dictates CO 2 diffusivity, prevails over the impact of the available alkaline reserve.
Keeping in mind the broad range of limestone powder contents (21-70%), the correlation (13) described with the R 2 equal to 0.78 can be considered strong enough for the simplified prediction model.It is important to note that the proposed equations are valid for concrete mixes that are wet-cured for at least 7 days.If shorter curing periods were included (1 and 3 days), a single relationship with a strong correlation between inverse natural carbonation resistance and compressive strength could not be obtained due to the significant impact of short periods on the carbonation resistance.This was however considered acceptable for the simplified prediction model.Moreover, curing periods shorter than 7 days are not likely in practice-they are rather an exception to the rule in special cases.

Distribution of Predicted R −1
NAC for Selected Values of Mean Compressive Strength Logarithmic transformation of the natural inverse carbonation resistance and the mean compressive strength as the regression variables was used to formulate the linear regression model given in Equation (10).This allowed us to establish the regression line as the expected response of ln R −1 NAC for a given ln f cm , as well as to compute the prediction (confidence) interval for the values of ln R −1 NAC for a given ln f cm .While this procedure is straightforward for linear regression, interpretation of the results with log-transformed variables needs careful consideration.For the sake of clarity, we use the following notation: y = ln R −1 NAC and x = ln f cm .In simple linear regression, the prediction interval is symmetric around the regression line (see Figure 1b, for example).This follows from the regression theory, which postulates that the residuals are normally (symmetrically) distributed with zero mean and variance σ 2 .Consequently, predicted response variable y is also normally distributed with conditional mean being the regression line and with the same variance σ 2 .Furthermore, the leastsquares estimators of regression coefficients b0 and b1 are random variables by their nature.The variability of b0 and b1 means that the estimated regression line (10) may deviate from a true one in terms of both slope and intercept.For a fixed value of the explanatory variable x = x 0 , the estimated conditional mean response µ Y|X=x 0 = b1 •x 0 + b0 is normally distributed as: where x is the mean value of the explanatory variable, S xx is the sum of squares of x, b 1 and b 0 are true regression coefficients, and n is the number of observations.The variance of the estimated prediction of the response variable ŷ0 for the given x 0 is obtained by combining the variances of the estimated regression model and the variance of the residuals ε: The above variance is used to compute the 100(1 − α)% confidence intervals for the predicted value y 0 of the response variable: where t n−2,1−α/2 is the Student's t variate with n − 2 degrees of freedom and probability 1 − α/2.The confidence interval of log-transformed R −1 NAC that is symmetric due to the normality assumption becomes asymmetric after exponentiation: (17) It should be noted, however, that the inverse transformation (exponentiation) of the linear regression equation does not yield the conditional mean response in the space of the original variable.Since the predicted ln R −1 NAC for the given x 0 is normally distributed, the predicted R −1 NAC is log-normally distributed.When exponentiated, the mean value of the log-normal distribution is not equal to the mean value of the normal distribution, although it is equal to the median value of normal distribution [61].Therefore, the power curves shown in Figures 1a, 2a and 3a do not represent the mean predicted response of R −1 NAC , but the median predicted response.The lower and upper confidence limits of the asymmetric 90% confidence interval in Figures 1a, 2a and 3a represent the same probabilities (5% and 95%) of the conditional distribution as in Figures 1b, 2b and 3b.
With log-normally distributed values of the predicted R −1 NAC for the given f cm , the relationship between the moments of the original and log-transformed normal variables was used to define the mean and variance of estimated R −1 NAC prediction for the given f cm as: where R denotes the estimated R −1 NAC prediction and ŷ0 denotes the estimated ln R −1 NAC prediction for the given x 0 , i.e., for the given f cm .This allows for computing the moments of the estimated R −1 NAC prediction for selected values of f cm and using them in the original space of R −1 NAC for the service life analysis.Distribution parameters for R −1 NAC have been determined for commonly applied concrete classes defined in EN 1992-1-1:2004 [57].The distribution parameters, concrete classes, as well as the corresponding f cm are shown in Table 4 for the 21-70% limestone powder concrete since service life analysis was performed for this group of concrete mixes.

Service Life Analysis
A summary of the input parameters for the limit state function of carbonation-induced depassivation is shown in Table 5. Obtained results of the FORM analysis are presented as relationships between the reliability index and period of reinforcement depassivation.Comments are also given regarding the requirements of the currently valid standard EN 1992-1-1:2004 [57].Figure 4 shows this relationship for a 50-year service life (t SL ) concrete cover (c nom ) for all exposure classes (XC1-XC4).For exposure class XC1, HVLPC mixes showed that the period of carbonation-induced reinforcement depassivation was 16 years for concrete class C20/25, which is defined as an indicative minimum concrete strength class [57].As the strength class increased, so did the service life.Although exposure class XC1 was the least aggressive environment, concretes up to class C30/37 did not meet the service life of 50 years (the highest depassivation period was 40 years).Concrete class C35/45 was the first to achieve a service life of 50 years with the prescribed concrete cover (t SL = 57 years).With a further increase in strength, the service life increased up to 79 years for the C40/50 concrete class.By solving the limit state function for the depassivation period of 50 years, the obtained concrete cover depth for strength class C20/25 was 42 mm, instead of the prescribed 25 mm [57].For exposure class XC1, HVLPC mixes showed that the period of carbonation-in duced reinforcement depassivation was 16 years for concrete class C20/25, which is de fined as an indicative minimum concrete strength class [57].As the strength class in creased, so did the service life.Although exposure class XC1 was the least aggressive en vironment, concretes up to class C30/37 did not meet the service life of 50 years (the high est depassivation period was 40 years).Concrete class C35/45 was the first to achieve service life of 50 years with the prescribed concrete cover (tSL = 57 years).With a furthe increase in strength, the service life increased up to 79 years for the C40/50 concrete class By solving the limit state function for the depassivation period of 50 years, the obtained concrete cover depth for strength class C20/25 was 42 mm, instead of the prescribed 2 mm [57].
A similar situation was observed for exposure class XC2, with the indicative strength class C25/30 and the minimum concrete cover of 35 mm proposed in [57].For this concret A similar situation was observed for exposure class XC2, with the indicative strength class C25/30 and the minimum concrete cover of 35 mm proposed in [57].For this concrete class, the service life was 18 years.The required service life would be achieved with a concrete cover of 56 mm (21 mm higher than the prescribed value).The only concrete class that satisfied the service life of 50 years was the class C40/50 (54 years).
A higher indicative minimum concrete strength class is prescribed for higher exposure classes.Concrete class C30/37 is prescribed as minimum indicative [57] for classes XC3 and XC4.For exposure class XC4, the values of the service life were similar to those for class XC2, for all concrete classes.The differences in service life duration were up to 1 year.This is not surprising because a 5-mm-higher cover depth is prescribed for the exposure class XC4 as being more prone to carbonation compared to class XC2.The required service life for class C30/37 would be achieved with a concrete cover of 52 mm (12 mm higher than the prescribed value).
The exposure class XC3 had the lowest service life, because moderate humidity represents the most aggressive environment for carbonation development.Although wetting and drying cycles, as defined for exposure class XC4, slow down the carbonation process compared to moderate humidity [59,60], a 5-mm-higher concrete cover depth is also prescribed compared to class XC3 [57].The higher prescribed concrete cover depth for class XC4 is a consequence of taking into account the possibility of rapid corrosion progress after depassivation has been reached [60], which was not considered as a part of the service life in this study.For exposure class XC3, not a single concrete strength class met the service life of 50 years with the prescribed concrete cover depth of 35 mm.The required concrete cover depth to satisfy the service life of 50 years, for the prescribed concrete class C30/37, was 56 mm (21 mm more than the prescribed value).
In total, the concrete cover depths for HVLPC mixes were larger by 12 to 21 mm (depending on the exposure class) compared to the prescribed covers in EN 1992-1-1:2004 [57].The reason for this is their lower carbonation resistance due to the binder dilution caused by replacing the clinker with limestone [23].
Bearing in mind that in EN 1992-1-1:2004 [57] a service life of 100 years is achieved by increasing the concrete cover depth to only 10 mm (Table 2), the depassivation period for this service life duration was also analyzed.Figure 5 shows the relationships between the reliability index and period of reinforcement depassivation for a prescribed 100-year service life of a concrete cover (c nom ) for all carbonation exposure classes.HVLPC mixes had the smallest depassivation period for exposure class XC3.Not single concrete class met the service life of 100 years.The longest service life was obtaine for the highest strength class C40/50 and it was 65 years.The depassivation period for th indicative minimum strength class C30/37 was only 32 years.The required concrete cove depth to satisfy the service life of 100 years was 77 mm (32 mm more than the prescribe value and 21 mm higher than the required 50-year concrete cover).
The analysis showed that in order to ensure a service life of 100 years, it is necessar to increase the 50-year concrete cover depth by approximately 20 mm in the case o HVLPC.The concrete cover depths were 23 mm larger for classes XC1 and XC4, i.e., 3 mm for classes XC2 and XC3, compared to the prescribed covers in EN 1992-1-1:2004 [57 Calculated values of the minimum concrete cover depths (cmin,dur) from the durability poin of view for commonly applied concrete classes are shown in Table 6.Given that the cal culation was made with the nominal cover depth (cnom), as explained in the Methodology For exposure class XC1, the period of depassivation was only 34 years for the indicative minimum concrete strength class C20/25.The required concrete cover depth for this strength class and for a depassivation period of 100 years was 58 mm, 23 mm higher than the prescribed value [57].With increased concrete strength, the service life also increased.Although exposure class XC1 was the least aggressive environment, as in the case of the service life of 50 years, only the concrete classes C35/45 and C40/50 met the service life of 100 years.
A similar service life was obtained for both XC2 and XC4 exposure classes.Only the C40/50 concrete class had a depassivation period of nearly 100 years.In the case of exposure class XC2, the depassivation period for indicative minimum concrete strength class C25/30 was 31 years.The required service life for this strength class would be achieved with a concrete cover of 78 mm (33 mm higher than the prescribed value and 22 mm higher than the required 50-year concrete cover).The depassivation period for concrete class C30/37 in the case of exposure class XC4 was 45 years.The required concrete cover depth to satisfy the service life of 100 years was 73 mm (32 mm more than the prescribed value and 21 mm higher than the required 50-year concrete cover).
HVLPC mixes had the smallest depassivation period for exposure class XC3.Not a single concrete class met the service life of 100 years.The longest service life was obtained for the highest strength class C40/50 and it was 65 years.The depassivation period for the indicative minimum strength class C30/37 was only 32 years.The required concrete cover depth to satisfy the service life of 100 years was 77 mm (32 mm more than the prescribed value and 21 mm higher than the required 50-year concrete cover).
The analysis showed that in order to ensure a service life of 100 years, it is necessary to increase the 50-year concrete cover depth by approximately 20 mm in the case of HVLPC.The concrete cover depths were 23 mm larger for classes XC1 and XC4, i.e., 32 mm for classes XC2 and XC3, compared to the prescribed covers in EN 1992-1-1:2004 [57].Calculated values of the minimum concrete cover depths (c min,dur ) from the durability point of view for commonly applied concrete classes are shown in Table 6.Given that the calculation was made with the nominal cover depth (c nom ), as explained in the Methodology, the minimum cover depth was calculated by reducing the nominal value by the value of the accepted negative deviation (∆c dev ).As expected, increasing the concrete class led to a decrease in the cover depth for all exposure classes.The biggest difference was for exposure classes XC1, XC2, and XC3, where it was necessary to increase the minimum cover depth for the indicative minimum class by two times compared to cement concretes of the same strength.However, the values shown in Table 6 enable a simpler application of high-volume limestone powder concretes in terms of durability (carbonation-induced corrosion).

Conclusions
Due to the dilution effect of the limestone, the carbonation resistance of the highvolume limestone powder concrete is lower than that of the pure-clinker concrete.This effect can present a limiting factor for the application of the HVLPC in the reinforced concrete structures despite the environmental benefit from significantly reducing the clinker content.Furthermore, the current standards [33,36,57] are not adapted to concrete mixes with such a type of binder.
The property that reflects the concrete carbonation resistance (whether this be the inverse carbonation resistance or carbonation rate coefficient) is determined on the basis of the measured carbonation depths under defined environmental conditions.Therefore, it has to be experimentally determined, which either is a long-lasting procedure or requires specialized equipment not available in every concrete laboratory.Hence, the attempt was made to establish an empirical relationship between the natural inverse carbonation resistances and mean compressive strength of the HVLPC that is to use the concrete property that has to be tested anyway.To do that, a regression analysis was performed on the experimental results collected from the previous research in the 1995-2023 periods.Applying certain filters on collected data led to 75 experimental results of the carbonation depth measured under natural indoor and sheltered outdoor conditions.Selected data originated from 14 different laboratories.For this analysis, concrete mixes were divided into two groups according to limestone percentage: (1) between 21% and 35%, and (2) between 36% and 70% limestone powder in the concrete powder phase.Based on the linear regression model of the log-transformed variables, empirical relationships were proposed, together with the statistical description of the models.Since the obtained equations for both concrete groups were similar, finally one single correlation was proposed for the whole 21-70% limestone powder content range.The coefficient of determination in this case was equal to 0.78.Keeping in mind the broad range of various parameters in the collected test data, the correlation described with such R 2 can be considered strong enough for the simplified prediction model.It is important to note that the proposed relationships are valid for concrete mixes wet cured for at least 7 days and indoor and outdoor sheltered natural conditions.
A full probabilistic service life analysis regarding the carbonation-induced depassivation of reinforcement was performed using the fib Model Code 2010 prediction model and proposed empirical relationship.The analysis was performed for all carbonation exposure classes, commonly applied concrete strength classes (C20/25-C40/50), and service lives of 50 and 100 years.Minimum concrete cover depths for HVLPC were recommended.The obtained results showed that for all exposure classes and both service life durations, the required minimum concrete cover depths c min,dur for HVLPC were significantly larger compared to those prescribed in EN 1992-1-1:2004 for indicative strength classes.For a service life of 50 years, the calculated minimum cover depths were 32 mm for XC1, 46 mm for XC2 and XC3, and 42 mm for XC4.For a service life of 100 years, the calculated cover depths were 48 mm for XC1, 68 mm for XC2, 67 mm for XC3, and 63 mm for XC4.Altogether, the increase varied between 40 and 110%, or between 12 and 33 mm, depending on the exposure class and service life duration.These results clearly indicate that HVLPC should be used for the structures where a service life of 50 years is acceptable.Adding a safety factor of 10 mm for cast-in-situ structures according to EN 1992-1-1:2004 leads to nominal covers of 75-80 mm for a service life of 100 years, which may not be acceptable in practice.
The basic limitations of this study are inherent to every empirical study: its reliability depends on the number and representativeness of the input data.Although analyzed test data originate from a number of different laboratories and cover a wide range of involved parameters, the total number was 75.Increasing the number of experimental results would certainly improve the quality of prognosis.Moreover, as already mentioned, reported results are valid for concrete mixes wet cured for at least 7 days and indoor and outdoor sheltered natural conditions.This is a consequence of the lack of test results.Therefore, more experimental testing of HVLPC carbonation resistance for different curing periods and under different natural conditions is recommended in future research.
The analysis in this study is performed under the assumption that the end-of-service life is the end of the initiation period, when the conditions for the start of corrosion are reached.However, a part of the propagation period, i.e., a certain amount of reinforcement corrosion, can be included in the service life definition, in addition to the initiation (depassivation) period.For instance, that part can be taken as the time until first cracks due to reinforcement corrosion appear.Such a design approach could enable broader HVLPC application in RC structures and should be a subject of the future research.
Figures 1-3 also show the 90% confidence intervals for predicted values of R −1 NAC .
Figures 1 through 3 also show the 90% confidence intervals for predicted values of R −1 NAC.

Figure 1 .
Figure 1.Relationship between the natural inverse carbonation resistance and mean compressive strength of concrete for mixes with 21-35% limestone powder: (a) power regression and (b) linear regression of log-transformed variables.

Figure 1 .
Figure 1.Relationship between the natural inverse carbonation resistance and mean compressive strength of concrete for mixes with 21-35% limestone powder: (a) power regression and (b) linear regression of log-transformed variables.Buildings 2023, 13, x FOR PEER REVIEW 8 of 18

Figure 2 .Figure 2 .
Figure 2. Relationships between the natural inverse carbonation resistance and mean compressive strength of concrete for mixes with 36% to 70% limestone powder: (a) power regression and (b) linear regression of log-transformed variables.

Figure 2 .
Figure 2. Relationships between the natural inverse carbonation resistance and mean compressive strength of concrete for mixes with 36% to 70% limestone powder: (a) power regression and (b) linear regression of log-transformed variables.

Figure 3 .
Figure 3. Relationships between the natural inverse carbonation resistance and mean compressive strength of concrete for mixes with 21% to 70% limestone powder: (a) power regression and (b) linear regression of log-transformed variables.

Figure 3 .
Figure 3. Relationships between the natural inverse carbonation resistance and mean compressive strength of concrete for mixes with 21% to 70% limestone powder: (a) power regression and (b) linear regression of log-transformed variables.

Figure 4 .
Figure 4. Reliability index (β) values during time of exposure for 50-year service life concrete cover (cnom) and different exposure classes.

Figure 4 .
Figure 4. Reliability index (β) values during time of exposure for 50-year service life concrete covers (c nom ) and different exposure classes.

Buildings 2023 , 1 Figure 5 .
Figure 5. Reliability index (β) during time of exposure for a 100-year service life concrete cover (cnom and different exposure classes.

Figure 5 .
Figure 5. Reliability index (β) during time of exposure for a 100-year service life concrete cover (c nom ) and different exposure classes.

Table 1 .
Range of concrete and environmental data in analyzed experimental results.

Table 2 .
[57]nominal concrete cover depths for service lives of 50 and 100 years according to[57].

Table 3 .
Linear regression of R −1 NAC on ln fcm: model adequacy as measured by coefficient of determination, estimated coefficients, and their properties.

Table 3 .
Linear regression of R −1 NAC on ln f cm : model adequacy as measured by coefficient of determination, estimated coefficients, and their properties.

Table 4 .
Moments of conditional log-normal distribution of predicted natural inverse carbonation resistance R −1NAC for the selected values of the mean compressive strength f cm .

Table 5 .
Input parameters for the limit state function of reinforcement depassivation.
* Lower and upper limit of the beta distribution.

Table 6 .
Recommended values of minimum concrete cover depths (c min,dur ) for 21% to 70% limestone powder concrete and different exposure classes.