Spatial Variability of Ozigo Wood Beams under Long-Term Loadings in Various Environmental Exposures

: Timber is a renewable material that should be more used for sustainable construction. While the mechanical behavior and durability of some species have been widely studied in the past, few studies are available for the Ozigo ( Dacryodes buettneri ) specie. This paper deals with the spatial variability of Ozigo beams subjected to long-term loadings and different environmental conditions. These beams were previously subjected to long-term creep in three environments (air-conditioned, unsheltered, and sheltered) at Masuku in the south-east of Gabon. Various specimens were extracted from these beams to determine its moisture content and subjected to three-point bending tests to obtain the modulus of elasticity and failure stress at various points in the space. The results obtained showed that, after long-term loadings, environmental exposure combined with mechanical loading, play a key role in the mechanical properties of the timber beams. A reduction of strength was found for the specimens extracted from the unsheltered and sheltered outdoor exposures in comparison with those extracted from the air-conditioned exposure. Concerning the spatial variability, statistical tests conﬁrm that there is signiﬁcant spatial correlation. It was also found that the spatial variation of properties in the beam is not stationary because it was affected by loading and support conditions.


Introduction
In tropical countries such as Gabon, the forest covers more than 80% of the territory with about 480 various types of species.To diversify its economy, the Gabonese government decided to ban the export of logs from January 2010, thus paving the way for local processing into finished products and promotion of this sustainable resource in the construction of timber structures that are characterized by a low environmental footprint [1][2][3][4].Predicting structural responses and estimating the durability of tropical timber structures subjected to thermo-hygro-mechanical loading [5] appears as an important challenge for promoting its use in the construction sector, especially for the Ozigo species (Dacryodes buettneri) which is an endemic species after Okume (Aucoumea klaineana Pierre).The Ozigo specie is a leafy tree of tropical origin from western Central Africa (Gabon and other equatorial countries).This specie is mainly used for plywood construction, in veneer, and finished or semi-finished products [6].
The mechanical behavior of timber structures depends strongly on climate variations [7].With over 1900 mm of annual precipitation, around 85% relative humidity and a mean temperature near 27 • C, the environmental conditions in Gabon are very arduous and impact the durability of timber structures during its service life [8].Therefore, these environmental conditions should be considered in the design of wooden structures to ensure adequate levels of serviceability and safety, especially for structures exposed in tropical environments.In addition, there is a variability of mechanical parameters due to the heterogeneous character of wood material that depends on the age, location of wood in the tree, loading time [9,10], climatic conditions [11,12], nodes, growth of annual rings [13,14], etc.This variability is both random and spatially distributed and is commonly called 'spatial variability' [15,16].
In the literature background, several authors investigated the spatial variability of materials [17].Some of them developed statistical and mathematical models, based on random field tools, to characterize, to represent or to propagate it.For example, Clerc et al. [18] and Rakotovao Ravahatra et al. [19] proposed procedures for characterizing the spatial variability from one-dimensional discrete data.Both studies were applied to reinforced concrete structures.The outcomes of these studies could be used to represent the spatial variability using several methods such as the Karhunen-Lòeve expansion, the optimal linear estimation method and/or Expansion Optimal Linear Estimation [20,21].A detailed overview of these methods is provided by Ghanem and Spanos [22] and Sudret and Der Kiureghian [23].The application of the outcomes and methods mentioned previously allow estimating the influence of spatial variability on safety of structures.Stewart [24] discussed about the spatial variability of damage and evaluated the related maintenance costs for the corroding reinforced concrete structures.This study highlighted the importance of considering the spatial variability to accurately estimate the probability of deteriorating structural components.Zhu et al. [25] also found that neglecting the spatial variability of reinforced concrete properties leads to an overestimation of the probability of compression failure.
Few works have been focused on the spatial variability of timber structures.Moshtaghin et al. [13,26] performed an experimental characterization of the longitudinal mechanical properties (elastic modulus and compressive strength) of light wood at different scales.They found that the dimension of the specimens should be carefully selected to avoid scale effects.A specific procedure to correct these scale effects was proposed by the authors; this fact should be accounted for when the size of the specimen is very small in comparison to the correlation length.Tran et al. [27] proposed a framework based on Bayesian networks to study the effect of the spatial variability of decay on the durability of timber beams.This framework could be used to evaluate the durability condition of the beam using partial observations; however, the spatial variability of mechanical properties was not considered in this work.
The main objective of this study is to analyze the variability and the spatial variability of mechanical properties (modulus of elasticity (MOE) and failure stress (FS)) of Ozigo beams exposed to various long-term environmental conditions and mechanical loading.Towards this aim, the beams were previously subjected to long-term creep (9 years) under three environments (air-conditioned environment, unsheltered outdoor environment, and sheltered outdoor environment) at the south-east of Gabon [28].In addition, this study also evaluates the effects of climatic and hydric parameters on the mechanical properties of wood subjected to different environmental conditions.This long-term work will address the crucial need to promote the sustainable development and use of wood-based materials and structures in Gabon and in the region.The study also responds to the practical need to adapt the Eurocode 5 for the design and verification of timber constructions in a tropical environment.In the mid-term, the main findings of this paper will contribute to a wider use of this renewable and sustainable material in the construction sector.

Materials
Ozigo is the tropical specie considered in this work.It has a light brown color and is a light wood with a density of 0.59 tons/m 3 ± 0.05 [29,30].The specimens were extracted from three wood beams (89 × 176 × 3090 mm 3 ) previously subjected to a creep test in tropical atmospheres from three environments during 5 years, Figure 1 [28].After the 5 years of loading, the beams were unloaded and remained exposed to their same environmental conditions over 3 years (recovery period) before being stored.Any treatment against decay was applied to study the natural durability of Ozigo beams under realistic exposure conditions; however, any presence of fungi or insects was observed after the exposure time.For the beams exposed outdoors, they were subjected to the rain and therefore to shrinkage swelling; however, this exposure has not generated cracks of the beams.
ical atmospheres from three environments during 5 years, Figure 1 [28].After the 5 years of loading, the beams were unloaded and remained exposed to their same environmental conditions over 3 years (recovery period) before being stored.Any treatment against decay was applied to study the natural durability of Ozigo beams under realistic exposure conditions; however, any presence of fungi or insects was observed after the exposure time.For the beams exposed outdoors, they were subjected to the rain and therefore to shrinkage swelling; however, this exposure has not generated cracks of the beams.
The beams were sawn into nine sections 340 mm long.Figure 2 shows the section of the beam and the exact location of each specimen extracted in the beams.These sections were subdivided into three zones in height (upper zone (UZ), middle zone (MZ), and lower (LZ)), Figure 2a,b.From each section were extracted: two specimens for bending test, two for tensile tests and four for compression tests.Note that this paper only presents and analyses the results from the three-point bending tests (81 samples of section 15 × 15 mm 2 and 300 mm long).The specimens concerned by bending tests are shown in Figure 2b and the real samples are presented in Figure 2c.Before loading the specimens, the mechanical behavior of wood inside each beam cannot vary.However, after long-term loading, the upper and tensile zones are in compressive and tensile areas, respectively.Consequently, specimens in each zone will be submitted to similar mechanical conditions during the tests.The beams were sawn into nine sections 340 mm long.Figure 2 shows the section of the beam and the exact location of each specimen extracted in the beams.These sections were subdivided into three zones in height (upper zone (UZ), middle zone (MZ), and lower (LZ)), Figure 2a,b.From each section were extracted: two specimens for bending test, two for tensile tests and four for compression tests.Note that this paper only presents and analyses the results from the three-point bending tests (81 samples of section 15 × 15 mm 2 and 300 mm long).The specimens concerned by bending tests are shown in Figure 2b and the real samples are presented in Figure 2c.Before loading the specimens, the mechanical behavior of wood inside each beam cannot vary.However, after long-term loading, the upper and tensile zones are in compressive and tensile areas, respectively.Consequently, specimens in each zone will be submitted to similar mechanical conditions during the tests.Figure 3 presents the experimental devices used to measure the geometry, moisture content and temperature of each specimen.These devices are: (i) a caliper (Figure 3a) for measuring the dimensions; (ii) a Romus Mini 93250 brand spiked moisture meter (Figure 3b) for determining the moisture content; and (iii) a thermo-hygrometer brand Testo 510 (Figure 3c) for estimating the temperature and relative humidity during the tests.All measurements were carried out before each bending test.
Sustainability 2021, 13, 5356 5 of 17 Figure 3 presents the experimental devices used to measure the geometry, moisture content and temperature of each specimen.These devices are: (i) a caliper (Figure 3a) for measuring the dimensions; (ii) a Romus Mini 93250 brand spiked moisture meter (Figure 3b) for determining the moisture content; and (iii) a thermo-hygrometer brand Testo 510 (Figure 3c) for estimating the temperature and relative humidity during the tests.All measurements were carried out before each bending test.Table 1 provides the number of specimens, and the average density and moisture content per zone.At least 8 specimens per zone and 26 specimens per beam are used herein.The average values of density and moisture content provided in Table 1 illustrate the variability of these parameters for the considered beam zones and environmental exposure conditions.However, since the coefficient of variation of the density is very small (about 3%), the statistical and spatial variability analyses will only focus on the moisture content and mechanical parameters.Table 2 provides the average exposure conditions (temperature T, relative humidity RH) and the moisture content at the beginning of the exposure [28].Based on the moisture content, we note that the original beams were loaded in a wet state (semi-dry and commercially dry wood), and thus dried under loading.By comparing the average values and standard deviation of the three exposures, it is observed that the temperature and relative Table 1 provides the number of specimens, and the average density and moisture content per zone.At least 8 specimens per zone and 26 specimens per beam are used herein.The average values of density and moisture content provided in Table 1 illustrate the variability of these parameters for the considered beam zones and environmental exposure conditions.However, since the coefficient of variation of the density is very small (about 3%), the statistical and spatial variability analyses will only focus on the moisture content and mechanical parameters.Table 2 provides the average exposure conditions (temperature T, relative humidity RH) and the moisture content at the beginning of the exposure [28].Based on the moisture content, we note that the original beams were loaded in a wet state (semi-dry and commercially dry wood), and thus dried under loading.By comparing the average values and standard deviation of the three exposures, it is observed that the temperature and relative humidity are lower for the air-conditioned exposure; but these parameters are very close for the unsheltered and sheltered conditions.Therefore, the main difference between unsheltered and sheltered conditions will be related to the exposure to rain and sun radiation for the unsheltered condition.Indeed, the exposure of the beams to rain and sun increase the risk of fungi attack, the shrinkage-swelling phenomenon, and creation of cracks and microcracks, which can in turn influence measured mechanical properties.
Table 2. Climatic characteristics of the three environments during the 9 years [28] and moisture content of beam before loading.

Exposure
Parameters T ( Due to the lack of a climate chamber, the specimens were not conditioned at 12% of moisture content.They were kept in an air-conditioned room where average temperature and relative humidity were 24.3 • C and 43.8, respectively.Since the moisture content is different for each zone (i.e., Table 1), the results of the mechanical tests were adjusted to a moisture content of 12% for comparison purposes by following the procedure indicated in Section 2.2.2.

Experimental Protocol
The experimental device used for the 3 points bending tests is an 'United' hydraulic press testing system machine of the series N • 0314523 (Figure 4).The maximum load capacity is 100 kN.The device is connected to a data acquisition system managed by the Datum 5.0 software.The bench test consists of two cylindrical horizontal supports (114.44 mm long and 12.80 mm of diameter).Two small aluminum plates 28.11 × 8.50 × 0.41 mm 3 are interposed between the test tube and supports.The test bench is equipped also of a displacement sensor, positioned at the mid-span of the specimen, to measure the arrow.The force increases gradually until the failure of specimen.
Sustainability 2021, 13, 5356 6 of 17 humidity are lower for the air-conditioned exposure; but these parameters are very close for the unsheltered and sheltered conditions.Therefore, the main difference between unsheltered and sheltered conditions will be related to the exposure to rain and sun radiation for the unsheltered condition.Indeed, the exposure of the beams to rain and sun increase the risk of fungi attack, the shrinkage-swelling phenomenon, and creation of cracks and microcracks, which can in turn influence measured mechanical properties.Due to the lack of a climate chamber, the specimens were not conditioned at 12% of moisture content.They were kept in an air-conditioned room where average temperature and relative humidity were 24.3 °C and 43.8, respectively.Since the moisture content is different for each zone (i.e., Table 1), the results of the mechanical tests were adjusted to a moisture content of 12% for comparison purposes by following the procedure indicated in Section 2.2.2.

Experimental Protocol
The experimental device used for the 3 points bending tests is an 'United' hydraulic press testing system machine of the series N ° 0314523 (Figure 4).The maximum load capacity is 100 kN.The device is connected to a data acquisition system managed by the Datum 5.0 software.The bench test consists of two cylindrical horizontal supports (114.44 mm long and 12.80 mm of diameter).Two small aluminum plates 28.11 × 8.50 × 0.41 mm 3 are interposed between the test tube and supports.The test bench is equipped also of a displacement sensor, positioned at the mid-span of the specimen, to measure the arrow.The force increases gradually until the failure of specimen.

Assessment of Mechanical Properties
The mechanical tests allow to determine two parameters: failure stress (FS) and Modulus of Elasticity (MOE).The determination of FS is given by Equation (1) [31]: where F max is the maximum force applied during the bending test, l is the length, b is the thickness and h is the height of the specimen (Figure 4).The values of MOE (Equation ( 2)) are obtained according to the NF EN 408 [31] standard: where ∆F is the increase in force between 0.1 to 0.4F max or 0.2 to 0.3F max , and ∆w is the increase in the arrow.
For a comparative study with the results obtained by Manfoumbi [28], the MOE is adjusted to equivalent values for a 12% of moisture content according to NF EN 338 [32] standard by using Equation ( 3): where E 12% represents the modulus of elasticity at 12% of the moisture content, and E w is the modulus of elasticity corresponding to the internal moisture content w of each specimen.

Statistical and Spatial Variability Analyses
The data obtained from experimental tests is statistically analyzed.The statistical analyses are useful to study the effects of environmental conditions on the studied parameters (MC, MOE, and FS) as well as to estimate their correlations.The first stage consists of gathering the data per beam (or exposure type) and to compare the mean, standard deviation, maximum and minimum values, and histograms for each studied parameter.In a second stage, a similar analysis is carried out by gathering data per zone (upper, middle, and lower).In the last stage, the Pearson correlation is computed to estimate the potential correlations between the three studied parameters per environmental condition and zone.
The main goal of the spatial variability analysis is to determine if the database is spatially correlated and stationary.The Ljung-Box (LB) test [33] is performed to examine if autocorrelation exists in this spatial data series.The null hypothesis of the LB test considers that the residuals are independently distributed.Once the existence of the autocorrelation is confirmed by the LB test, the Augmented Dicky Fuller (ADF) test [34] is carried out to determine the stationarity of the spatial data series.The null hypothesis for the ADF test implies that the data series is not stationary.

Statistical Analysis of Measured Data
Figure 5 presents the values of moisture content, modulus of elasticity and failure stress obtained for the three environments for all specimens extracted from the beams (Figure 2).The MOE and FS presented in this figure were adjusted by Equation (3) to consider an equivalent moisture content of 12%.In general, it is observed that all parameters vary for each beam and exposure conditions.Therefore, a statistical analysis will be first carried out to study the influence of the exposure conditions.The data is first gathered per beam/zone and analyzed for each kind of exposure.The histograms for each considered parameter and each beam are summarized in Figure 6.These histograms show different patterns of average values and dispersion of data for each exposure.The mean, standard deviation, and coefficient of variation (CoV) for each exposure/beam/zone are given in Table 3. Concerning moisture content, it is observed that the lower mean value was obtained for the beams under unsheltered exposure.This behavior is due to the fact that the beam was exposed to a dried period before the tests.For the modulus of elasticity and the failure stress, it is noted that higher mean values for beams are estimated for the air-conditioned exposure.This means that the strength of the beam was preserved when it was subjected to controlled environmental conditions.For MOE and FS, the dispersion of the values for the sheltered beams covers almost all the range of the values for air-conditioned and unsheltered exposure conditions (Figure 6).The mean values of the module of elasticity and failure stress are compared to those obtained by Manfoumbi [28] that are also presented in Table 3.The average MOE values were estimated from the instantaneous deflection of the beams measured after loading.The average FS values were computed from samples extracted from the beams but without any differentiation concerning its exposure.For this reason, Table 3 provides the same values for all exposures.It is noted that the mean values of MOE for all beams are lower than that the indicated by Manfoumbi [28].This decrease of the wood rigidity can be explained by the shrinkage-swelling, creep and mechano-sorptive phenomena, which are more accentuated in environments exposed to the weather than in the air-conditioned environment.Indeed, loss or absence of water in the wood leads to the disappearance of some links [35] such as low-energy (van Der Waals) bonds between water molecules and wood polymers (cellulose, hemicellulose, and lignin); therefore, timber unlinked fibers could be cracked due to water loss.The results obtained partly explain this observation since we observed a very low loss of rigidity (2%) in an air-conditioned environment, compared to sheltered outdoor environments (21%) and unsheltered (35%) in which the surrounding conditions change during the year.Further analysis of the timber microstructure will be very useful to better understand the mechanisms leading to the wood rigidity reduction.Concerning the FS, our values are in average 12 to 16% higher than the reported by Manfoumbi [28].Nevertheless, the values reported by Manfoumbi [28] are in the range of those found in this work.
Comparing the average values of the parameters in the zones, it is observed in Table 3 that most part of minimum values for each exposure are in the upper zone.Conversely, most maximum values are in the middle zone.The lower zone mainly contains medium average values.Any pattern was identified for the CoV because the range of variation for each zone is not very large.The difference between the zones could be related to the long-term behavior where tensile/compressive stresses reduced the resistance in the upper/lower zones, respectively.The middle zone, that is normally less loaded, preserving the maximum resistance.The reduction of MOE and FS for the unsheltered and sheltered exposure is due to timber aging, loading duration, and exposure conditions [36][37][38][39][40][41][42].All these parameters considerably influence the resistance of wood due to its nature and anisotropic behavior [43,44].More particularly, the variations in the moisture content, associated with the exposure and loading, play a key role in the longevity of the wooden structure, because they influence its aging process.The loss of secondary metabolites [45] such as tannins (polyphenols) would also be the cause of decrease in rigidity of the test samples since the tannins act as defender of wood structure by limiting the growth of wood biological agents [46] and the influence of climatic variations (modification of the conformation of molecular structures) which can alter the main polymers of wood [47].However, other chemical mechanisms such as the strong crystallization of celluloses [48] in wood can greatly reduce the stiffness of the material.Further analysis of chemical changes will provide additional information allowing to better understand their role on the loss of rigidity.
Determining the potential correlation between the parameters is important when the information will be after used for probabilistic modeling purposes.The Pearson correlations and p-values for the studied parameters (MC, MOE, and FS) are given in Table 4.It is noted that the absolute values of the correlations r MC−MOE and r MC-FS are lower than 0.56 with p-values larger than 5% for the most part of the cases.This indicates the moisture content is not linearly correlated with the mechanical parameters (MOE and FS).This is explained by the fact that the moisture contents depend on the surrounding weather conditions at the time when the tests were realized.Contrarily, it is possible to state that there is a positive correlation between the mechanical parameters MOE and FS because the values of r MOE-FS are close to 1 with p-values lower than 5% for most part of cases.This positive correlation should be considered for propagating uncertainties or spatial variability for these parameters.

Spatial Variability Analysis
The spatial distribution of the measured parameters given in Figure 5 shows that there exists significant spatial variability for all collected data.The LB test was carried out to confirm the existence of a spatially correlated database.The statistics and p-values for the LB test applied to the data obtained in this work are given in Table 5.By considering a significance level of 0.05, it is observed that in all cases the p-values are larger than this value confirming that autocorrelation is statistically significant for all spatial data series.Figure 5 also shows that there are particular values in the areas of support and application of loads.To better illustrate these particular behaviors, Figure 7 provides the spatial variability of each parameter estimated by doing an average of the three values per location x estimated for upper, middle, and lower zones for each exposure.An average value gathering all measurements (nine values per location x) is also shown in this figure (green line).Lower values are observed close to the support zones for the moisture content.This behavior is due to the fact that the beams are more exposed to the environment in the support areas.Some trends with low values of MOE and FS are also observed in the support zone close to 2 m.These trends could be related with the loading conditions that induced more concentrated stresses in these support zones.Since these trends could imply that the spatial variability of the studied parameters is not stationary, statistical tests will be carried out to evaluate the stationarity of these random fields.The analysis of autocorrelation of a spatial database allows to quantify the spatial regularity of a phenomenon and to determine the extent of the spatial dependence [49].Figure 8 presents the autocorrelation of the residuals for all beams and exposure conditions.It is observed in all cases that the autocorrelation decreases when the measurements are distant (larger lag #).This suggests that exists a spatial dependency structure [50,51]; however, in several cases, the autocorrelation values are significantly larger/lower than zero indicating a complex autocorrelation structure.Therefore, ADF tests were carried out to confirm the existence of the autocorrelation and to determine if the random fields are stationary.The statistics and p-values for the ADF tests are also provided in Table 5.For the same significance level (0.05) the estimated p-values indicate that the data series is not stationary.This result is not surprising considering the trends observed in Figure 7 for the support or loading areas.
tions.It is observed in all cases that the autocorrelation decreases when the measurements are distant (larger lag #).This suggests that exists a spatial dependency structure [50,51]; however, in several cases, the autocorrelation values are significantly larger/lower than zero indicating a complex autocorrelation structure.Therefore, ADF tests were carried out to confirm the existence of the autocorrelation and to determine if the random fields are stationary.The statistics and p-values for the ADF tests are also provided in Table 5.For the same significance level (0.05) the estimated p-values indicate that the data series is not stationary.This result is not surprising considering the trends observed in Figure 7 for the support or loading areas.

Conclusions
In this paper, the influence of long-term loading and three exposure conditions on the spatial variability of the moisture content, the modulus of elasticity and the failure of three Ozigo beams are investigated.The results show that the strength of the timber decreases for the specimens exposed to an unsheltered outdoor environment.The strength loss is also larger for the upper and lower zones of the beam than to the middle zone because of the long-term loading.These effects could significantly reduce the service life of timber structures exposed to sheltered and unsheltered outdoor environments.The experimental database allows estimating the correlations between the studied parameters.The correlation analysis showed that the moisture content is not linearly correlated with the modulus of elasticity and the failure stress.However, a positive correlation was estimated between the modulus of elasticity and the failure stress.
In addition, the spatially distributed data was used to study the spatial variability of the studied parameters.First, a LB test statistically confirmed the existence of spatial correlation.Then the spatial distribution of the data and the shapes of the autocorrelation

Conclusions
In this paper, the influence of long-term loading and three exposure conditions on the spatial variability of the moisture content, the modulus of elasticity and the failure of three Ozigo beams are investigated.The results show that the strength of the timber decreases for the specimens exposed to an unsheltered outdoor environment.The strength loss is also larger for the upper and lower zones of the beam than to the middle zone because of the long-term loading.These effects could significantly reduce the service life of timber structures exposed to sheltered and unsheltered outdoor environments.The experimental database allows estimating the correlations between the studied parameters.The correlation analysis showed that the moisture content is not linearly correlated with the modulus of elasticity and the failure stress.However, a positive correlation was estimated between the modulus of elasticity and the failure stress.
In addition, the spatially distributed data was used to study the spatial variability of the studied parameters.First, a LB test statistically confirmed the existence of spatial correlation.Then the spatial distribution of the data and the shapes of the autocorrelation functions suggested that there are some non-stationary trends related to support and loading conditions.This non-stationarity was confirmed by performing the ADF test.These results highlight the importance of the loading and exposure conditions on the spatial variability of the studied materials that induced non-stationarity in the spatial variability.Additional tests on unloaded beams will be very useful to confirm that this nonstationarity is mainly due to loading and support conditions.Future work will also focus on the characterization and then modeling of random fields considering exposure conditions and loading.In addition, the proposed methodology will be applied to other tropical and tempered species.These data and models will ultimately be used for the evaluation and updating of reliability of timber structures by considering spatial variability [27].

Sustainability 2021, 13 , 5356 4 of 17 Figure 2 .
Figure 2. Spatial representation of the beam: (a) Beam cutting plane.(b) Arrangement of specimens in a section of the beam.(c) Cross-section of the specimens from the three environments.Figure 2. Spatial representation of the beam: (a) Beam cutting plane.(b) Arrangement of specimens in a section of the beam.(c) Cross-section of the specimens from the three environments.

Figure 2 .
Figure 2. Spatial representation of the beam: (a) Beam cutting plane.(b) Arrangement of specimens in a section of the beam.(c) Cross-section of the specimens from the three environments.Figure 2. Spatial representation of the beam: (a) Beam cutting plane.(b) Arrangement of specimens in a section of the beam.(c) Cross-section of the specimens from the three environments.

Figure 4 .
Figure 4. (a) Three-point bending test bench.(b) Experimental device for three-point bending tests.Figure 4. (a) Three-point bending test bench.(b) Experimental device for three-point bending tests.

Figure 4 .
Figure 4. (a) Three-point bending test bench.(b) Experimental device for three-point bending tests.Figure 4. (a) Three-point bending test bench.(b) Experimental device for three-point bending tests.

Figure 5 .
Figure 5. Spatial distributed data for the considered parameters and zones.

Figure 6 .
Figure 6.Probability densities for each parameter and kind of exposure (a) moisture content, (b) modulus of elasticity, (c) failure stress.

Figure 7 .
Figure 7. Spatial variability of average values of the studied parameters (a) moisture content, (b) modulus of elasticity, (c) failure stress.

Figure 8 .
Figure 8. Autocorrelation of the studied parameters.

Figure 8 .
Figure 8. Autocorrelation of the studied parameters.

Table 1 .
Number of specimens per beam and average density and moisture content per zones.

Table 1 .
Number of specimens per beam and average density and moisture content per zones.

Table 2 .
[28]atic characteristics of the three environments during the 9 years[28]and moisture content of beam before loading.

Table 3 .
Mean, standard deviation, maximum and minimum values for each parameter, zone, and kind of exposure.

Table 4 .
Pearson's correlations and p-values for the studied parameters and exposures.

Table 5 .
Statistics and p-values (in brackets) for the LB and ADF tests carried out for studied parameters and exposures.