Abstract
The penetration of chlorides in concrete is the main reason for the beginning of corrosion in reinforced concrete structures exposed to a marine environment, reducing their service life. This article proposes correlations between the surface chloride content (Cs) of Fick’s model and the ambient factor (kCl) of a resistivity model. Literature data for three types of cement were used to simulate chloride diffusivity: (1) CEM I 52.5 R; (2) CEM II/A-P 42.5 R (with 6%–20% of pozzolans); and (3) CEM type II (with 12% of pozzolans and 8% of silica fume). Concretes containing these cements were analyzed for different environmental conditions: atmosphere zone, seawater immersion, and tidal/splash zones. The surface chloride contents were determined by a combination of Fick’s and resistivity models. In this article, it was noticed that the estimated Cs value varies with the exposure class from 0.484% to 0.644% (total chloride concentration by weight of cement). An equation has been proposed to correlate Cs with the ambient factor kCl of the resistivity model.
Introduction
Chloride penetration in concrete
Reinforcement corrosion is the degradation process that has the most impact on the durability of reinforced concrete structures. The natural conditions of high alkalinity in concrete normally protect the reinforcement from corrosion reactions. However, some agents may induce corrosion even in an alkaline environment, such as chloride, which destroys the passivation layer in the reinforcement–concrete interface.
Chloride-induced corrosion is an important subject studied in the literature. Balestra et al. (2019a) studied the reinforcement corrosion risk of marine concrete structures evaluated through electrical resistivity. Farzad et al. (2019) studied the effect of moisture on macrocell current and evaluated the suitability of ultra-high performance concrete as a repair material. Rossi et al. (2020) evaluated the effect of defects at the steel–concrete interface for chloride-induced pitting corrosion using the X-ray computed tomography. Thus, the effects of different variables on corrosion kinetics have been studied over the years.
Modeling the chloride penetration in concrete is a recurrent subject in researches (Andrade et al., 2013; Mazer et al., 2018; Petcherdchoo, 2013). Especially in the last 30 years, these models progressed considerably in terms of the complexity of the input variables and the number of factors taken into consideration (Meira et al., 2014; Torres-Luque et al., 2014). Some chloride penetration models are presented below.
Saetta et al. (1993) used a numerical procedure based on the finite element method to model the chloride ingress in concrete (equation (1))
where Dc1 and Dc2 (cm2/s) are the reference and adjusted chloride diffusion coefficients, respectively; T1 is the reference temperature (296 K) and T2 (K) is the current temperature; Ea is the activation energy of diffusion (kJ/mol), which is related to the water/cement ratio of the concrete; R is the universal gas constant (8.314 J/mol K); ζ measures how much the diffusivity decreases with time and varies from 0 to 1; te (days) is the actual time of exposure to chloride; RHc is the relative humidity at which Dc1 drops halfway between its maximum and minimum values; and RH is the current relative humidity considered (%).
Bastidas-Arteaga et al. (2010) used a stochastic approach to determine the influence of weather and global warming in chloride ingress into concrete. Mazer et al. (2018) proposed a mathematical model to estimate the chloride diffusion coefficient using fuzzy logic, considering parameters of the concrete including water/cement ratio, compressive strength, and concrete curing temperature. Wang et al. (2019a) proposed a semi-empirical prediction model of chloride-induced corrosion for application in uncracked reinforced concrete.
Balestra et al. (2019b) recently presented the modified Holliday equation in order to model chloride profiles (equation (2)) of field naturally degraded structures that present a sharp peak of maximum concentration, contemplating the convection and diffusion zones
where C(x,t) is the chloride concentration (% by weight of concrete);
Several other models can be found in the literature to predict chloride penetration. However, the application of most of these models presents difficulties mainly due to the following three factors: (1) boundary conditions that reproduce the reality of chloride penetration; (2) variables that change with the exposure time and the location; and (3) equations too complex to be applied under practice conditions.
Despite the difficulties discussed, service life models are extremely useful for the durability design of concrete structures. These models provide the opportunity to pre-arrange a maintenance program and, therefore, to guarantee the minimum service life specified in standards (above 50 or 100 years). Another benefit is the possibility of considering the local environmental variables acting on different parts of the structure in a microclimate concept.
Electrical resistivity model
Electrical resistivity indicates the amount and connection of pores in the concrete and can be used as a parameter to model the transport of different agents in this material (Medeiros-Junior et al., 2019). The resistivity and diffusivity relation can be established according to Einstein’s law (equation (3)), which relates to the movement of electric charges and the conductivity of the medium (Andrade, 1993). Factor A in equation (3) is dependent on the external ionic concentration. In the case of chloride ions, Andrade (2004) named A as kCl (cm3·Ω/year)
where Def is the effective diffusion coefficient (cm2/year), which is also known as the steady-state diffusion coefficient; and ρef is the effective resistivity at 28 days on saturated conditions (Ω·cm).
Andrade (2004) proposed a chloride diffusivity model based on concrete electrical resistivity. The reaction factor (rCl) is the factor responsible for transforming the Def (steady state) diffusion coefficient into the Dap (non-steady state) diffusion coefficient. Def considers only the ionic transport while Dap also takes into account the binding of chlorides with cement matrix. According to Castellote and Andrade (2006), the ratio between Def and Dap (=rCl) addresses the binding ability of the cement matrix. This reaction or binding factor is responsible for taking into account the delay in chloride penetration due to chloride binding and depends mainly on the cement composition. Thus, by definition, rCl is equal to Def/Dap (Andrade 2004, 2014). Def and Dap parameters can be obtained through a multi-regime test (UNE 83987, 2014) and then applied to equation (4) to determine rCl. Porosity (ε, % by volume) multiplied by Dap is used in order to homogenize the units in which Def and Dap are expressed
Dividing both sides of equation (3) by rCl factor, the apparent diffusion coefficient (in non-steady state) is determined using equation (5) (Andrade, 2004). Therefore, when the first-half part of equation (5) is divided by rCl, Def becomes Dap. The concepts of effective and apparent for electrical resistivity are similar to concepts for the diffusion coefficient mentioned before. The rCl factor represents the number of times the effective resistivity is apparently increased when the chloride combination is considered (Andrade, 2004; Andrade et al., 2014). That is, as part of the initially free chloride will become chemically bonded with the cementitious matrix, the resistivity of the concrete tends to increase because less free chloride will be available to interfere with the electrical conductivity. Therefore, the apparent resistivity (ρap = ρef·rCl) was proposed by Andrade (2004), which takes into account the slower chloride advance
where Dap is the diffusion coefficient in non-steady state (cm2/year); rCl is the chloride reaction factor; ρap is the concrete apparent resistivity (Ω·cm); and kCl is the ambient factor (cm3·Ω/year).
Finally, using the well-known “square root law” (equation (6)), the corrosion initiation period (ti) is estimated by equation (7) (Andrade, 2004), considering the moment of depassivation of the reinforcement as a threshold condition
where ti is the corrosion initiation period (years); x is the penetration depth of the chloride threshold concentration, which is equal to concrete cover (cm) when ti is considered; and ρt is the resistivity corrected over time.
In addition, the concrete age factor (q) should be used in equation (7) to consider that the concrete electrical resistivity increases over time, as shown in equation (8) (Andrade, 2004). ρt is the effective resistivity corrected over time due to the effect of concrete aging. Over time, concrete tends to increase electrical resistivity due to the evolution of cement hydration reactions and the consequent pore refinement (Mendes et al., 2018). Thus, the effect of increased resistivity over time due to concrete aging is taken into account in the Andrade (2004) model for chloride penetration by the factor q.
where t is the time considered for application of the age factor; t0 is the age of first resistivity measurement, which is 28 days; and q is the age factor of resistivity.
Andrade and D’Andrea (2010) calculated kCl values to different exposure classes defined by the European Standard EN 206-1 (2007). These kCl values were proposed based on rough estimations of experimental results of real structures exposed for a long period to the aggressive attack of chloride ions.
According to Andrade and D’Andrea (2010), the environment factor (kCl) is equal to 5000, 10,000, 17,000, and 25,000 cm3·Ω/year to the exposure classes XS1a, XS1b, XS2, and XS3, respectively. These classes are classified as follows (EN 206-1, 2007):
XS1: parts of structures in contact with marine aerosol, but without direct contact with seawater. It represents cases of structures in the coastal zone or close to the coastal zone. This class can be subdivided according to the distance to the shoreline as XS1a (distance > 500 m) and XS1b (distance < 500 m).
XS2: parts of marine structures subjected to permanent seawater immersion.
XS3: parts of marine structures located in tidal or splash zones.
Due to the difficulty of finding concrete structures with all variables required and available for a complete durability study, the number of samples analyzed to estimate the environmental coefficient for this model was relatively scarce. Therefore, it is assumed that the results are close to the desired conditions, although probably different from the overall average value (Andrade and D’Andrea, 2010).
Model based on Fick’s law
The chloride diffusion model (equation (9)) based on Fick’s second law (Collepardi et al., 1972; Crank, 1975) is well known and applied in the studies of the service life of the reinforced concrete structure (Medeiros-Junior et al., 2015; Wang et al., 2014)
where Cx is the chloride concentration (%) at depth x and time t; Cs is the surface chloride content (%); x is the depth (cm); Dap(t) is the apparent diffusion coefficient in time (cm2/s); t is the exposure time (s); and erf is the mathematical error function.
According to Mangat and Molloy (1994) and Maage et al. (1996), a mathematical model expressed by a power function (equation (10)) may estimate the effect of concrete maturity on diffusion
where t1 is the maturity age of the concrete—a time considered for application of the maturity factor; and t0 is a reference time for the first diffusion coefficient measurements generally assumed to be 28 days. The effect of diffusion coefficient reduction over time is represented by n in equation (10). According to Andrade et al. (2011), n can be related to the age factor q from the resistivity model using equation (11)
Discussion concerning resistivity and Fick’s models
Some similarities can be identified by observing the resistivity (equation (7)) and Fick’s (equation (9)) models, since both equations model the chloride transport in concrete. Table 1 shows some similarities of each variable used in these models.
The similarity between variables from resistivity and Fick’s models.
According to Table 1, the environment factor (kCl) from the resistivity model is equivalent to the Cs on Fick’s model. Considering the time needed for the depassivation of the reinforcement, Cx is considered as the threshold chloride concentration (Cthr) at the depth equal to the concrete cover (x = xc), and t is the service life of concrete structures (t = ti).
Thus, the aim of this article is to verify which Cs of Fick’s model is equivalent to several environment factors (kCl) from the resistivity model for different cement types by numerical simulation based on data already available in the literature. This relationship was identified for each of the exposure classes proposed in the EN 206-1 (2007) standard for chloride attack (XS1a, XS1b, XS2, and XS3). Experimental data obtained from the literature were used for simulations of chloride diffusivity and for studying the correlations between both environmental factors (Cs and kCl).
The advantages and significance of this study are to stimulate the use of durability models in predicting chloride penetration. The lack of defined variables for different types of environment and materials inhibits the use of these mathematical equations. Therefore, although good models for predicting chloride diffusivity exist, their practical application is hampered by the lack of data. This article attempts to narrow this gap by defining an equation to obtain environment factor data of Fick’s model from the resistivity model and vice versa.
In practice, determining the concentration of chlorides on the concrete surface (Cs) is time-consuming and expensive. However, this variable is essential for the model based on Fick’s second law (equation (9)). Even though some sensors (Meira et al., 2010) can be used for this determination, the application in real structures is still rare. Due to this difficulty, sometimes the concentration of chlorides in the first millimeters is assumed as being Cs. Some studies have also determined Cs by fitting the profiles of chlorides on Fick’s law (Medeiros-Junior et al., 2015; Wang et al., 2014). However, this practice has been questioned in the literature, because it results in an extrapolation of Cs and ignores the peak of chlorides (Andrade et al., 1997; Ann et al., 2009; Balestra et al., 2019b). The main reason for this article is trying to relate Cs with kCl values that were previously determined for different marine environments. This methodology highlights the importance of this study for the advancement of concrete durability studies in the marine environment.
In addition, electrical resistivity is a non-destructive technique that allows the continuous evaluation of processes occurring in a concrete structure by characterizing the electrical resistance inside the material. Electrical resistivity has the advantage of characterizing the corrosion kinetics, and it has also been used as a parameter for concrete design (Mendes et al., 2018). Therefore, this article defines environmental factor values for some common types of marine environments (exposure classes: XS1a, XS1b, XS2, and XS3) through an easily applied equation that allows the conversion of environmental factors from two different models.
Methods
The chloride reaction factor (rCl), age factor (q), resistivity at 28 days (ρef), and apparent diffusion coefficient at 28 days (Dap) for three usual types of cement were taken from previous researches (Andrade et al. (2014); Andrade (2014)). Three different cement types were selected for this study: CEM I 52.5 R and CEM II/A-P 42.5 R (Portland-composite cement with 6%–20% of pozzolans (PN)), according to the EN 206-1 (2007) standard, and CEM type II moderate sulfate resistance cement with 12% of PN and 8% of silica fume (SF). Andrade et al. (2014) used seven test specimens for CEM I 52.5 R, three test specimens for CEM II/A-P 42.5 R, and four test specimens for CEM type II with 12% PN + 8% SF. The amount of cement (kg/m3) ranges from 300 to 450, and the water/cement ratio ranges from 0.28 to 0.60. The details regarding material specifications can be found in Andrade et al. (2014).
According to the data available in Andrade et al. (2014), the chloride reaction factor (rCl) ranges from 0.7 to 2.6 (CEM I 52.5 R), 3.5 to 6.2 (CEM II/A-P 42.5 R), and 2.5 to 3.7 (CEM type II with 12% PN + 8% SF). The resistivity at 28 days (ρef) ranges from 6591 to 12,742 Ω·cm (CEM I 52.5 R), 4533 to 5072 Ω·cm (CEM II/A-P 42.5 R), and 20,365 to 45,299 Ω·cm (CEM type II with 12% PN + 8% SF). The apparent diffusion coefficient at 28 days (Dap) ranges from 6.24 to 43.38 × 10−8 cm2/s (CEM I 52.5 R), 3.85 to 11.5 × 10−8 cm2/s (CEM II/A-P 42.5 R), and 2.62 to 4.93 × 10−8 cm2/s (CEM type II with 12% PN + 8% SF). The age factor (q) is 0.22 for CEM I 52.5 R and 0.37 (recommended value for up to 20% PN) for CEM II/A-P 42.5 R and CEM type II with 12% PN + 8% SF.
Since CEM II has additions (PN and/or SF) in its composition, this explains the lower diffusion coefficient of CEM II when compared to CEM I, due to the pore refinement caused by these additions. The reaction factor is higher due to the presence of aluminates, which chemically fix a greater amount of chlorides in the cementitious matrix. The age factor for CEM II is also higher since the pozzolanic reactions are slower and progressive (Andrade, 2014; Medeiros-Junior and Lima, 2016).
In this article, these cement types were analyzed under different environmental conditions of chloride attack: XS1a, XS1b, XS2, and XS3, proposed by the EN 206-1 (2007) standard. The kCl values used for each class were suggested by Andrade and D’Andrea (2010), which are already discussed in the introduction of this article.
The first age for the measure of resistivity and apparent diffusion coefficient (t0, from equations (8) and (10)) was 28 days. The exposure period (t1, from equations (8) and (10)) considered was 91 days due to the time needed for the cement hydration, since at this age, over 80% of the cement is already hydrated (Taylor, 1997). The chloride threshold concentration assumed was 0.4% by weight of cement, considering total chloride concentration. The breakdown of the passivating film and subsequent corrosion process begins when a certain concentration of chloride ions reaches the reinforcement surface. The threshold value is highly relevant in chloride penetration studies in concrete structures and can be expressed in different ways. Some common examples are the percentage by weight of cement or concrete or from the Cl−/OH ratio of the solution present in the pores of the cementitious matrix. In a comprehensive literature review of the critical chloride ion content, Angst et al. (2009) commented that due to the existence of several factors involved (reinforcement characteristics, cement type, water/cement ratio, and temperature), the establishment of a single critical content for application in all real conditions is not reasonable. However, the value of 0.4% by weight of cement is one of the most used values in the literature (Meira et al., 2014; Mohammed and Hamada, 2006; Pacheco and Polder, 2016; Wang et al., 2019b; Wu et al., 2019). Although considered conservative, this value is adopted by standards in some countries—such as Brazil, Portugal, and England—in order to meet most exposure situations of a reinforced concrete structure to chlorides.
The combination of the models (resistivity and Fick) results in equation (12), which was used to determine Cs for each selected condition. Equation (12) comes from the substitution of t from equation (7) in equation (9). The concrete cover (xc) does not interfere in the analyses since this value is present in the numerator and the denominator of equation (12). A constant (=31,536,000) must be used in equation (12) to transform time t1 from years to seconds. The rearranged equation is presented in equation (13)
Results and discussion
Figures 1 to 4 show the surface chloride content (Cs) for the exposure classes XS1a, XS1b, XS2, and XS3, respectively. These values are based on a threshold concentration, considering the total chloride concentration in concrete (Cthr = 0.4% by weight of cement).

Surface chloride contents to XS1a.

Surface chloride contents to XS1b.

Surface chloride contents to XS2.

Surface chloride contents to XS3.
According to Figures 1 to 4, Cs did not present great variations between the types of cement. The standard deviations ranged from 0.023 to 0.085. This indicates that the surface concentration of chlorides is a measure that depends more on the exposure environment—the proximity of the coast, temperature, wind, and relative humidity—than on the cement type. The coefficients of variation ranged from 4.79% to 13.17%. Figure 5 shows the mean value of Cs for each exposure environment.

Mean surface chloride contents for each exposure environment.
According to Figure 5, Cs grows with the aggressiveness of chlorides of the environment. This behavior seems reasonable since the surface chloride content tends to be reduced once the concrete structure is further away from the coast. The Cs value of marine structures located in tidal or splash zones is 33.2% higher when compared to that of structures at 500 m from the shoreline. Meira et al. (2006, 2007) observed this behavior by monitoring the chloride concentration through equipment called “wet candle.”
The results of Figure 5 were statistically analyzed, in order to verify if the differences are due to the effective changes in the control variables. The one-way analysis of variance (ANOVA) method was applied, with a significance level of 5%, followed by a multiple-variance analysis to verify the differences between the exposure classes. Tables 2 and 3 present the results of the statistical analysis.
One-way ANOVA.
ANOVA: analysis of variance; SQ: sum of squares; FD: freedom degree; AS: average squares.
Multiple-variance analysis.
According to Tables 2 and 3, all classes generated concentrations that are statistically different from each other. According to Figure 5, Cs varies with the exposure class from 0.484% to 0.644%. These results seem reasonable, according to the data of surface chloride content measured and estimated in the field by other studies in the literature. Table 4 shows several surface chloride concentrations found in the literature. According to Table 4, there is no consensus on the Cs value on previous studies related to chloride ingress in concrete. Different authors have found Cs ranging from 0.4% to more than 4.0%, either through chloride profile fitting or laboratory experiments. However, in general, some tendencies can be noticed: the results are higher for tidal > spray zone > atmosphere. Furthermore, Cs grows with time, with the water/cement ratio, and with the proximity to the seashore. As shown in Table 4, the following studies found Cs from 0.4% to 0.7% by weight of cement: Mustafa and Yusof (1994), Costa and Appleton (1999), Yoon (2007), Meira et al. (2010), O’Connor and Kenshel (2013), Safehian and Ramezanianpour (2013), Valipour et al. (2013), and Kim et al. (2016). These values are very close to the results in Figure 5, confirming and validating the results obtained in this article.
Surface chloride concentration found in the literature.
W/B: water/cement ratio; OPC: ordinary portland cement; FA: fly ash; SLAG: blast furnace slag; PFA: pulverised fuel ash; MK: metakaolin; SF: silica fume.
An explanation for the variation of Cs found in Table 4 is the chloride threshold concentration (Cthr). The Cthr assumed in this article was 0.4%; however, there is also no consensus in the literature about the limit chloride concentration to initiate the corrosion process. Values from 0.2% up to 1.5% have already been documented and summarized by Meira et al. (2014). Figure 6 shows an example for a class XS1a of how Cs depends on the Cthr parameter.

Cs and Cthr relationship for XS1a.
Another relevant consideration is the variation of Cs in time. However, according to Ann et al. (2009), the time-dependent characteristic of the surface chloride content has been rarely considered in chloride ingress. Therefore, the values proposed by this article refer to the final concentration of Cs, after their stabilization in time. For further study on this topic, Petcherdchoo (2013) showed a literature review on how to consider the increasing of Cs in time. Figure 7 shows the relationship found between Cs of Fick’s model and kCl of resistivity model.

Cs and kCl relationship according to the conditions used in this study.
The proposed correlations in Figure 7 show a high R2 (0.98). This equation allows determining the surface chloride content in the concrete as a function of the exposure classes (kCl) of the concrete resistivity model and vice versa. These proposed correlations are based on the application of calculation models conducted in this article and require validation in real concrete structures for safe use in service life models.
Conclusion
This article proposes correlations between the surface chloride content (Cs) of Fick’s model and the ambient factor (kCl) of a resistivity model. All environmental classes have shown Cs values that are statistically different from each other. Cs did not show great variations between the types of cement. This indicates that the surface concentration of chlorides is highly dependent on the exposure environment. Cs grows with the aggressiveness of chlorides of the environment since the surface chloride content tends to increase when the concrete structure is closer to seawater. The estimated Cs values vary with the exposure class from 0.484% to 0.644% (by weight of cement). These values are very close to some results found in the literature. Although surface concentration may be time-dependent, the values proposed by this article refer to the Cs after their stabilization in time. An equation has been proposed to correlate Cs with the ambient factor kCl of the resistivity model. The results contribute to facilitating the application of both models for the conditions used in this article. An experimental validation would be essential to reinforce the values obtained in this article and the limits of application of the proposed equation. Therefore, although probably different from the overall average value, the results are assumed to be close to some real conditions.
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
