Inﬂuence of Soil Moisture vs. Climatic Factors in Pinus Halepensis Growth Variability in Spain: A Study with Remote Sensing and Modeled Data

: The inﬂuence of soil water content on Aleppo pine growth variability is analyzed against climatic variables, using satellite and modeled soil moisture databases. The study was made with a dendrochronological series of 22 forest sites in Spain with different environmental conditions. From the results of the correlation analysis, at both daily and monthly scales, it was observed that soil moisture was the variable that correlated the most with tree growth and the one that better identiﬁed the critical periods for this growth. The maximum correlation coefﬁcients obtained with the rest of the variables were less than half of that obtained for soil moisture. Multiple linear regression analysis with all combinations of variables indicated that soil moisture was the most important variable, showing the lowest p -values in all cases. While identifying the role of soil moisture, it was noted that there was appreciable variability between the sites, and that this variability is mainly modulated by water availability, rather than thermal conditions. These results can contribute to new insights into the ecohydrological dynamics of Aleppo pine and a methodological approach to the study of many other species.


Introduction
It is generally appreciated that forests, especially those located in water limited areas, are one of the environments most affected by climate change [1]. At the same time, it is well known that ecosystems play a crucial role in the interactions between natural systems that face climate change and the atmospheric system [2]. Forests are directly involved in critical issues for natural systems such as energy balance [3], carbon dioxide budget [4], or freshwater availability [5], and forest evolution is sharply conditioned by climate change.
A better understanding of how environmental factors shape plant function is crucial for predicting the consequences of environmental change for ecosystem functioning [6]. Precise knowledge of the factors controlling vegetation dynamics is also essential for improving ecosystem monitoring [7]. Climatic factors are identified as the main determinants in vegetation distribution, and climate dynamics are considered to drive vegetation shifts in many areas of the world [8,9]. Nevertheless, other factors, such as soil water content or topography and its derived variables, may be of great importance in the interaction of environmental factors affecting forest dynamics [10].
Global warming and the subsequent water scarcity may negatively impact forest productivity in extensive areas in many regions [11]. Forests of the Mediterranean region will be among the most affected because climate change exacerbates Mediterranean stresses [12]. A better understanding of tree growth responses to climate and to other related drawbacks could provide valuable insight into the vulnerability of those forests.
Aleppo pine is the most widespread pine species in the Mediterranean basin [13,14], and therefore is very representative of the bioclimatic conditions in this water-limited environment [15,16]. For these reasons, Aleppo pine has captured the interest of many researchers in recent decades from different points of view and all around the Mediterranean [17][18][19][20]. The study of the influence of environmental factors on the growth of this species has been approached from diverse perspectives, but in most cases has been focused on climatic and terrain variables [21].
For example, in one location in southern Italy, [22] studied Aleppo pine growth compared with weekly precipitation and the integrated temperature excess above 8 • C. Ref. [23] analyzed radial growth and wood density profiles with precipitation and temperature variables in one location in southern France. Ref. [24] studied the Aleppo pine forest growth response to climate by analyzing radial growth series and their relationship to precipitation and temperature conditions at different spatial scales. Ref. [25] studied the wood formation process and the intra-annual density fluctuations in Aleppo pines correlated with monthly temperature and precipitation. Ref. [26] analyzed tree-ring widths and wood anatomical features in tree rings under contrasting environmental conditions in Spain and Slovenia. They applied principal component analysis to elucidate the relationship between those tree parameters and climate, based on the monthly mean temperature and the monthly sums of precipitation. Ref. [21] analyzed climatic variables and local site conditions, and concluded that the local water balance was the main factor driving Aleppo pine productivity in their study area in the South of France.
Usually, the few studies in which soil moisture has been used have been conducted using in situ measurements or estimations at the forest stand scale, and with very shortlived series, of a few growing seasons or years. Accordingly, very few studies have been conducted using soil water content as a variable to analyze tree growth dynamics and variability in large areas and with suitable long series. This fact is particularly relevant considering that this variable refers to the volume of water supplied to plants [27], and for that reason, it is generally considered a relevant factor. Annual growth of trees is governed by numerous biotic and abiotic factors, with available soil moisture being one of the most important [28], as has been known for a long time [29]. For example, [30] studied the relationships between climate, soil moisture, and phenology in two woody areas in West Africa using in situ soil moisture monitoring stations. Ref. [31] studied the role of soil moisture, estimated by an ecohydrological model, in Aleppo pine growth-climate relationships in four locations in Spain. Ref. [32] studied the relationship between soil moisture and forest productivity in a small forest catchment in USA using a forest growth model. The need to pay more attention to soil moisture is in line with the finding that climate change will intensify the soil water content deficit in the coming years in many regions of the world [33].
The lack of studies associating these topics with applied soil moisture is mainly due to the difficulty in obtaining soil moisture databases with long series and over large areas. Soil moisture monitoring requires costly sensors, is time consuming and, most of the time, requires direct surveying. In the past, it has been very difficult to obtain long-term observations because monitoring campaigns are usually linked to specific experiments for one or several years, at best. In many cases, only point measurements are available, and they are hardly applicable to areas of interest.
A suitable strategy to overcome these constraints could be the use of soil moisture databases generated by remote sensing, hydrological modeling, or reanalysis. Several satellites are now providing global soil moisture retrieval using both active [34] and passive [35,36] microwave technologies. These satellite missions provide soil moisture products for both the surface and root-zone layer, with suitable revisit times and at a variety of spatial resolutions. Many of these databases still have a limited timespan for use in long-term analysis. However, other satellite databases, such as the Climate Change Initiative (CCI) Soil Moisture [37] series ongoing since 1978, allow long-term analysis for multiple applications. Although it is known that the retrieval of soil moisture with Remote Sens. 2021, 13, 757 3 of 20 microwave remote sensing is complicated in areas with dense vegetation [38], there are increasing numbers of examples where satisfactory results have been obtained [39,40]. Ref. [41] found that soil moisture estimates from microwave-based active and passive satellites outperform several model-based products over vegetated areas analyzed at global scale. Ref. [42] validated the CCI database under different environmental conditions in the Iberian Peninsula and obtained good results over forest areas (correlation coefficients around 0.9, and very low error and bias values).
Another way to apply soil moisture data in forest studies is to use databases obtained from hydrological modeling or reanalysis. There are many options and modeling approaches that generate soil moisture as a model output and provide databases at a variety of spatial resolutions, with long-term series and for many regions, or even globally [43][44][45][46]. In some cases, a hydrological model has been applied to a specific site [31] to generate soil moisture outputs at forest stand scale. The main limitation is that the required information at the desired scale and location is not always available. A suitable approach to overcome that constraint is the use of available soil moisture databases obtained through modeling [47], and previously validated in the area of interest.
Recently, [48] demonstrated that satellite soil moisture is sensitive enough to track the Aleppo pine growth pattern along the year. The satellite soil moisture was able to detect a bimodal pattern of tree growth with a maximum in May and a secondary peak in autumn. This temporal pattern was much clearer than that obtained using the modeled soil moisture or the precipitation. The study showed that, despite the limited spatial resolution of microwave satellite databases, the products that are available, such as CCI, can be useful for ecohydrological analysis and monitoring of this species. The opportunity to access global and daily root-zone soil moisture via satellite databases provides a huge advantage over other approaches that have been used so far.
The main goal of the present study was to investigate the role of soil moisture to explain the variability of the Aleppo pine growth in 22 locations in Spain, with different environmental conditions, from the north to the south of the country. With the aim of analyzing the influence of soil moisture on tree growth, several statistical approaches were applied, and two soil moisture databases, one from satellites and one from modeling, were used for the first time. Accordingly, another goal of the work was to study the influence of soil moisture compared with several climatic variables commonly used in this kind of studies. The climatic variables and indicators analyzed were minimum and maximum temperature, precipitation, standardized precipitation index (SPI), and potential evapotranspiration (PET). Finally, the role of soil moisture on tree growth was analyzed based on the climatic characteristics of each location, to study whether there is environmental modulation of the influence of the soil water content on Aleppo pine phenology.

Tree Samples and Study Area
We used a previously developed dataset, selecting the annual tree-ring width series of 22 Aleppo pine forest stands covering a wide range of environmental conditions across the Iberian Peninsula ( Figure 1 and Table 1). Site selection was based on the homogeneity of the forestland cover and soil conditions (marls and limestone) at the sampled sites [48]. In each forest site, a total number ranging from 15 to 40 dominant or co-dominant Aleppo pine trees, separated by at least 10 m, with no obvious signs of damage or disease, were sampled. Trees were cored at 1.3 m height, at two perpendicular directions, using Pressler increment borers.
Remote Sens. 2021, 13, 757 4 of 20 using a negative exponential function followed by a cubic smoothing spline with a 50% cutoff frequency and a 30-year response period [15,49]. The autocorrelation filter was then applied to the detrended series to remove correlations between consecutive measurements and to obtain residual series containing only the high-frequency variations in yearto-year growth series. Finally, a bi-weight robust estimation of the mean was applied to construct each local residual chronology [50].    Thus, using the annual tree-ring width series of 22 Aleppo pine forest stands, the dendrochronological database was derived using standard dendrochronological protocols [49]. Then, for each forest, a ring-width index (TRI) chronology was constructed by detrending tree-ring width measurements to eliminate the biological age trend in the radial growth. To preserve the inter-annual scale variability and eliminate the age-related trend in the radial growth, we standardized the individual tree-ring width series data were standardized using the ARSTAN program [49]. The long-term trends were removed using a negative exponential function followed by a cubic smoothing spline with a 50% cutoff frequency Remote Sens. 2021, 13, 757 5 of 20 and a 30-year response period [15,49]. The autocorrelation filter was then applied to the detrended series to remove correlations between consecutive measurements and to obtain residual series containing only the high-frequency variations in year-to-year growth series. Finally, a bi-weight robust estimation of the mean was applied to construct each local residual chronology [50].

Soil Moisture
For this study, two soil moisture databases from satellites and modeling have been used. The first one corresponds to version 4.4 of the CCI Soil Moisture Combined product satellite database [51]. From the three different CCI soil moisture products available, the combined product was selected for this study because it is the longest satellite database, and it has given the best results in validations conducted worldwide, including the Iberian Peninsula [42,52]. This database provides daily surface soil moisture from seven passives and four actives satellites that have soil moisture measurements from November 1978 to June 2018 over a regular 25 km × 25 km grid. A detailed description of the algorithm used to obtain this product from the merge of the L2 products of these eleven satellites can be found in [37,53]. In this research, the pixels corresponding to the 22 tree samples were selected, and the surface soil moisture data was filtered applying the soil moisture flag ('no data inconsistency'), in order to use the values with highest quality. Once the surface soil moisture time series from the CCI database were selected and filtered, the Soil Water Index (SWI) model [54] was used to obtain the soil moisture in the root zone. This model only uses as input the surface soil moisture and a parameter T related to the water time travel through the soil profile. To obtain this parameter, the SMAP L4 Surface and Root Zone product [55], which provides daily surface and root zone (0-100 cm depth) soil moisture over a regular grid of 9 km × 9 km from March 2015 to December 2019, was used. From this product, the surface and root zone soil moisture time series corresponding to the 22 tree samples were selected. Then, the surface soil moisture was used as input in the SWI model, and the T parameter was varied from 1 to 100 days to obtain 100 root zone soil moisture time series for each tree sample. These 100 root zone time series from the SWI model were compared with the root zone soil moisture from the SMAP L4 product to select the best T parameter for each tree sample based on the highest correlation, according to the study of [56]. Once the best value of T was selected for each tree sample, this parameter, together with the surface soil moisture time series from the CCI combined product, was used in the SWI model to obtain the CCI root zone soil moisture time series of the first 100 cm depth.
The second soil moisture database used in the present study is the set of soil moisture estimations from the Lisflood (LF) model [47]. LF is a spatially distributed hydrological rainfall-runoff model developed by the floods group of the Natural Hazards Project of the Joint Research Centre (JRC) of the European Commission [57]. LF was validated by [58] along Europe, with satisfactory results over the Iberian Peninsula. The soil moisture provided by this model is used by the European Drought Observatory (EDO) monitoring system. This hydrological model provides daily soil moisture data at three different layers with a spatial resolution of 5 km × 5 km and is available from January 1991 to December 2018. Only the time series of the first two layers, named the topsoil [59], were selected for the pixels corresponding to the 22 tree samples. The average of these two layers was considered the root zone soil moisture, assuming a soil volume similar to that of the CCI root zone, i.e., 0-100 cm depth, and in accordance with the root zone depth of the Aleppo pine [60].
These two databases were compared in a previous study [48], and have been widely used in studies with different soil moisture applications [61][62][63][64].

Precipitation and Temperature
Precipitation and temperature are the climatic variables most commonly used in studies of the impact of environmental factors on tree growth and phenology [65,66] in Remote Sens. 2021, 13, 757 6 of 20 general and of Aleppo pine growth specifically [20,26,67]. The precipitation and temperature data used in the present study were obtained from the Spanish PREcipitation At Daily scale (SPREAD) and Spanish TEmperature At Daily scale (STEAD) databases, respectively [68,69]. Both databases are provided on a 5 km × 5 km resolution grid over the Iberian Peninsula, with daily data for precipitation, maximum temperature, and minimum temperature. Both databases were created with measurements from more than 12,000 meteorological stations with reference values through generalized linear models. These were based on the 10 closest measurement stations using latitude, longitude, and altitude as covariates. A detailed description of this computational process to reconstruct the original daily precipitation time series can be found in [68]. These two datasets have been widely validated and used in many works since they are the spatially distributed databases for these two variables with the greatest temporal coverage (1950-2012 for precipitation, 1901-2014 for temperature) over the Iberian Peninsula [69][70][71][72].
As for the soil moisture databases, only the time series of precipitation and maximum and minimum temperature pixels corresponding to the 22 locations of the tree samples were used in this research.

Standardized Precipitation Index (SPI)
The SPI is a meteorological drought index [73] and is calculated from the historical record of the precipitation collected at a location, where the accumulated precipitation during a specific study period is compared with the same period of time throughout the long-term record history in the same location. This index is considered by the World Meteorological Organization as the universal reference meteorological index [74].
The SPI generator application from the National Drought Mitigation Center (NDMC) (https://drought.unl.edu/droughtmonitoring/SPI/SPIProgram.aspx, accessed on 22 December 2020) of the University of Nebraska was used to calculate the SPI-1. The SPI obtained in this study is based on a representation of the historical record of precipitation with a gamma distribution, where the positive values of SPI correspond to wet conditions, while the negative values correspond to drought conditions. The long-term time series, from 1950 to 2012, of SPREAD precipitation data from the 22 locations of the tree samples were used as input to obtain the monthly value of SPI-1 in those locations. Only SPI-1 was analyzed, since only daily and monthly temporal scales were considered for the rest of the climatic variables used in this study. SPI was used even though precipitation itself was also included because SPI is a suitable indicator of precipitation anomalies.

Potential Evapotranspiration
Evapotranspiration is an essential variable in the water balance and clearly identifies the interaction of vegetation in the hydrological system. It is also well known that evapotranspiration is decisively involved in the growth of trees in different ecosystems [65]. In the present work, PET was correlated with Aleppo pine growth with the aim of studying whether a temporal pattern or a critical period exists.
The Climatic Research Unit gridded Time Series (CRU TS) v4.03 is a database that provides monthly data of a set of climatic variables over a regular 0.5 • × 0.5 • grid from 1901 to 2018 [75]. Due to this long-term time series, this dataset has been widely used in different climatic studies worldwide [76][77][78]. For this study, PET time series from the CRU TS were obtained from the pixels corresponding to each tree sample. This variable is calculated through the Penman-Monteith equation [79,80].
Prior to analysis, and to test the suitability of using PET with such a spatial resolution, CRU PET data were validated in all sampling sites. As sufficient climatic information was not available in these areas for using the Penman-Monteith method, PET was calculated with the STEAD database and the Thornthwaite method. The result of the comparison between CRU PET and PET-Thornthwaite was very satisfactory, with an average R 2 of 0.93. Finally, in this study, CRU PET was used considering that the Penman-Monteith method has more physical basis.

Correlation Analysis
In order to study the role of soil moisture and climatic factors in Aleppo pine growth and identify any temporal patterns, a correlation analysis was performed between the annual TRI series and the different variables considered in this study. For this, Pearson's correlation coefficient (R) was calculated, and the statistical significance was settled for p-values < 0.05. This analysis was carried out at daily and monthly temporal resolutions and applied to a study period spanning from January 1979 (January 1991, in the case of LF) to December 2012.
For the analysis at daily resolution, the soil moisture from both LF and CCI, precipitation, and maximum and minimum temperature series were analyzed. The daily time series of these variables were obtained by averaging the original daily time series with a 30-day moving window. These values were only computed for those days when at least half of the data in the window were available. By doing so, the variability of the original series was smoothed, and the time patterns could be identified more clearly. Then, the annual TRI series were correlated with the series of each day from the beginning of the hydrological year, i.e., from 1st October of the previous year, to 31st December of the corresponding year, as in [48]. Thus, 457 R values were obtained for each tree-ring chronology. In this way, it is possible to identify any temporary pattern for each of the variables. Finally, the monthly median of the 22 R values and the percentage of significant results for each day were also calculated with the aim of simplifying and improving the representation of the results.
The same analysis was carried out for the monthly series. For the variables with daily resolution series, monthly values were obtained by calculating the average in those months where at least half of the daily data were available. The variables used for this case were soil moisture from both LF and CCI, precipitation, maximum and minimum temperature, SPI-1 and PET. The annual TRI series were correlated with the series of each month from October of the previous year to December of the corresponding year. Thus, 15 R values were obtained for each study site. In the same way, the median of the 22 resulting R data sets and the percentage of significant results for each month were also calculated, similar to the approach for the daily scale. In addition, for this analysis, the standard deviation between samples for each month was calculated.
With the aim of investigating whether the relationship between soil moisture and Aleppo pine growth dynamics is modulated by environmental conditions or whether that relationship follows a geographical pattern, the R (LF vs. TRI) data obtained were correlated with the mean annual temperature, maximum temperature, minimum temperature, precipitation (P), potential evapotranspiration (PET), P/PET ratio (aridity index), and altitude, latitude, and longitude of every forest site.

Multiple Linear Regression Analysis
With the aim of analyzing the influence of soil moisture on tree growth compared with climatic variables, a multiple linear regression (MLR) analysis was performed. Monthly series were used; the independent variables that were incorporated into the MLR model were the soil moisture, precipitation, SPI, PET, and maximum and minimum temperature; and the dependent variable was the TRI. Due to the data gaps in the CCI database during the first half of the series, and therefore the difficulty of obtaining a sufficiently long monthly series, this analysis was only performed with the LF database. Multicollinearity between independent variables is identified as a problem when trying to study influences on a dependent variable [81]. The variance influence factor (VIF) can detect the presence of multicollinearity. Therefore, in order to avoid this issue, different MLR analyses combining the independent variables were performed with the most correlated variables omitted and ensuring a VIF < 5 for the remaining variables [82]. This resulted in six different MLR analyses in which the soil moisture was maintained with the precipitation or SPI and the PET, maximum temperature or minimum temperature. For each analysis, the p-value was calculated for each independent variable, thus obtaining its importance in the MLR model. Therefore, a monthly p-value was obtained for each variable and each sample from October of the previous year to December of the corresponding year. Finally, the median of the 22 resulting p-values and the percentage of p-values < 0.05 for each variable for each month were also calculated. Figure 2 shows the temporal pattern of the relationship of each of the climatic variables studied with Aleppo pine growth, with soil moisture, from both CCI and LF, being most highly correlated along the seasonal cycle. For the LF database, the period with highest correlation spans from April to August, with correlations ranging from 0.4 to 0.6, reaching 80% of significant cases. This period, identified as critical for the Aleppo pine growth, starts in spring, as shown for the CCI database, but lasts until mid-summer. However, the secondary peak was not detected by LF. This might be because the LF model was not able to adequately reproduce the soil water recharge that takes place after the summer under Mediterranean conditions. In fact, [31], which also used modeled soil moisture, found both the main spring period extending to mid-summer and the secondary autumn peak seen with the CCI database. Regarding the number of significant samples, the difference obtained between the two soil moisture databases might be due to the smaller number of data points available in the CCI database along the first half of the series. Another difference observed between the two sets of soil moisture database results is the variability between samples. The 25th and 75th percentiles of the correlation coefficient with the LF database differ from the median by between 0.2 For the satellite soil moisture database, the maximum correlation was obtained during spring, i.e., from April to June, along with the increment in significant cases. The maximum of both median R value (0.78) and number of significant cases (50%) was reached in mid-May. In addition, all correlation coefficients exceed the significance level during these months. Ref. [31] obtained similar R values for Aleppo pine samples located in semiarid areas of the southeast Iberian Peninsula. A secondary peak was also observed in October of the previous year, with R values at approximately 0.5; although it does not reach the level of significance, it coincides with an increase in the number of significant samples, as seen in [48]. These results detecting this secondary peak are in agreement with those obtained Remote Sens. 2021, 13, 757 9 of 20 in the study carried out by [83]. These authors found two growth phases in Aleppo pine, coinciding with the critical soil moisture periods detected in this study. The first one was observed in spring, the period of maximum correlation of soil moisture. The second one was observed in autumn, coinciding with the time when the soil water recharge occurs under this kind of climate condition.

Analysis of the Influence of Soil Moisture and Climatic Factors
For the LF database, the period with highest correlation spans from April to August, with correlations ranging from 0.4 to 0.6, reaching 80% of significant cases. This period, identified as critical for the Aleppo pine growth, starts in spring, as shown for the CCI database, but lasts until mid-summer. However, the secondary peak was not detected by LF. This might be because the LF model was not able to adequately reproduce the soil water recharge that takes place after the summer under Mediterranean conditions. In fact, ref. [31], which also used modeled soil moisture, found both the main spring period extending to mid-summer and the secondary autumn peak seen with the CCI database. Regarding the number of significant samples, the difference obtained between the two soil moisture databases might be due to the smaller number of data points available in the CCI database along the first half of the series. Another difference observed between the two sets of soil moisture database results is the variability between samples. The 25th and 75th percentiles of the correlation coefficient with the LF database differ from the median by between 0.2 and 0.3 for the whole period, while with the CCI database, this difference reached 0.7 in mid-July days. This might be related to the different spatial resolutions of the two products used.
When the climatic variables were considered in the daily temporal scale analysis, the R values obtained did not exceed 0.4 in either case (Figure 2c,d,e). Some temporal patterns can be observed in the relationship of these factors with the growth of Aleppo pine when an increase in the correlation coefficients coincides with an increase in the number of significant samples. Nevertheless, these correlations are weaker than those obtained for soil moisture, and the variability found between the different samples was also smaller than that obtained for this variable. The precipitation showed a peak of correlation between May and June (R = 0.37), exceeding the significance threshold on a few days, but was much lower than that obtained during the critical period identified with soil moisture. This correlation vanishes in summer, coinciding with the period when the cambial activity of the pine ceases due to dry conditions [25,83]. In the study carried out by [31], it was observed that the precipitation also showed a high correlation with tree growth in spring months. Ref. [84] determined that this variable was the main factor influencing the Aleppo pine growth with correlation values of approximately 0.4 in spring. In addition, [85,86] obtained similar results in the south of the Iberian Peninsula, as did [87] in the northeast. A similar temporal pattern was obtained by [88] for Aleppo pine trees located in Greece. This period, when precipitation seems to have more influence on tree growth, coincides with the period observed to coincide with soil moisture in those months when it seems to be more determinant. However, all the R values obtained for precipitation in the previously cited works, and in this study, are lower than those obtained for the soil moisture, both with CCI and LF.
In the case of temperature, a temporal pattern of correlation can be seen, with two peaks of relative importance, but moderate R values (Figure 2c,d). The first one, with positive R values (0.36), is observed with the minimum temperature in February, and the second peak, with negative R values (−0.34), is observed with the maximum temperature in July. Similar results were obtained in other works [31,86,[89][90][91][92]. This implies that high temperatures have a negative influence on the growth of Aleppo pine in summer, while warm minimum temperatures have a positive influence in winter [87,93]. This could be in agreement with the well-described plasticity in the annual rhythms of cambial activity of the Aleppo pine in response to seasonal climatic variations [83,94,95]. Thus, the growth of the species during the growing season is usually subject to "double stress", characterized by two stops or slowdowns in cambial activity, one during winter, caused by low temperatures, and one during the summer, triggered by high temperatures and a lack of precipitation [96,97].
Temperature has been an environmental factor commonly used to explain the variability of Aleppo pine growth. Correlation values similar to those obtained in the present study were obtained in other works [89,90,92]. However, the results obtained both in this study and in those cited above, indicated less correlation than that obtained with the soil moisture.
The results obtained at monthly temporal resolution (Figure 3) showed similar temporal patterns to those obtained with the daily analysis. Once again, soil moisture was the variable that was the most highly correlated with Aleppo pine growth, and this was seen with both CCI and LF. The maximum correlation period in this analysis was identified from April to June with the CCI database, while with LF, it extended to August. Moreover, the R median values were greater for the CCI database during the period when soil moisture seems to be critical for Aleppo pine growth. However, the secondary peak observed in autumn with the satellite database in the previous analysis was not clearly observed here. For LF soil moisture results, both correlation coefficients and significant cases lightly decreased for the entire period at this time resolution. This would indicate that daily scale would be preferable as the temporal scale for this kind of analysis, as critical periods could be shorter than one month. Ref. [31], using the same time scale for their analysis and modeled soil moisture, obtained stronger correlation results than the LF results and similar results to those of the CCI, but using different soil depths (from 30 to 50 cm). When the variability of the results between samples was studied, for CCI, the standard deviation ranged between 0.4 to 0.2, while in LF, these values ranged between 0.3 and 0.2 throughout the period, similar to that obtained in the daily scale analysis. In addition, when the CCI soil moisture reached its maximum correlation values, a decrease in the standard deviation was observed, indicating that soil moisture has a more homogeneous influence on Aleppo pine growth between samples during this period.
As in the analysis carried out at the daily scale, this analysis also showed a lower relationship between climatic variables and TRI series than those obtained with the soil moisture ( Figure 3). The R values obtained were always lower than 0.3, and barely 40% of significant samples were found in those months where R reaches its maximum value. Regarding the variability, an almost constant standard deviation ranging between 0.1 and 0.2 was observed for the climatic variables, exceeding in several months the mean R value, similar to the variability results obtained in the study by [24]. Precipitation and maximum and minimum temperature reached the maximum correlation coefficients in the same months as in the daily analysis, although with lower values: 0.25, -0.29, and 0.26, respectively. PET reached its maximum correlation in February (R = 0.23) and July (R = −0.26), coinciding with the periods of maximum correlation for both maximum and minimum temperatures. The correlation values obtained were similar to those obtained with precipitation and temperatures; however, these maximum correlations were not even half of that obtained for soil moisture. In the case of SPI, the variable reaches its maximum correlation values in January and May (R = 0.28 in both cases). Recently, several drought indices have been used to study the impact of extreme climatic situations on tree growth, and in some cases, a temporal pattern was detected [92,[98][99][100][101]. In the present study, no temporal pattern was clearly identified with this variable, but the results obtained were consistent with those obtained by [102]. The Aleppo pine has been shown to be more resistant to more extreme drought conditions than other pine species [92]. However, it is interesting to highlight that an increase in this variable indicates better conditions in terms of water availability and positively affects tree growth [100]. The results obtained in the present study show a slight influence of this precipitation anomaly in winter and spring periods, although its relevance is much lower than that observed for soil moisture itself. CCI soil moisture reached its maximum correlation values, a decrease in the standard d viation was observed, indicating that soil moisture has a more homogeneous influence o Aleppo pine growth between samples during this period. As in the analysis carried out at the daily scale, this analysis also showed a low relationship between climatic variables and TRI series than those obtained with the so moisture ( Figure 3). The R values obtained were always lower than 0.3, and barely 40% significant samples were found in those months where R reaches its maximum value. R garding the variability, an almost constant standard deviation ranging between 0.1 an 0.2 was observed for the climatic variables, exceeding in several months the mean R valu similar to the variability results obtained in the study by [24]. Precipitation and maximu and minimum temperature reached the maximum correlation coefficients in the sam months as in the daily analysis, although with lower values: 0.25, -0.29, and 0.26, respe tively. PET reached its maximum correlation in February (R = 0.23) and July (R = −0.26 coinciding with the periods of maximum correlation for both maximum and minimu temperatures. The correlation values obtained were similar to those obtained with preci itation and temperatures; however, these maximum correlations were not even half of th obtained for soil moisture. In the case of SPI, the variable reaches its maximum correlatio values in January and May (R = 0.28 in both cases). Recently, several drought indices ha been used to study the impact of extreme climatic situations on tree growth, and in som cases, a temporal pattern was detected [92,[98][99][100][101]. In the present study, no temporal pa tern was clearly identified with this variable, but the results obtained were consistent wi those obtained by [102]. The Aleppo pine has been shown to be more resistant to mo extreme drought conditions than other pine species [92]. However, it is interesting to hig light that an increase in this variable indicates better conditions in terms of water avail bility and positively affects tree growth [100]. The results obtained in the present stud Therefore, the critical period for soil moisture coincided with the period when the maximum correlation values of precipitation and SPI were reached. However, the correlation of the climatic variables at any temporal scale was even less than a half of that obtained with both soil moisture databases. The results of the correlation analysis seem to indicate that the fundamental variable for Aleppo pine growth is water availability, and the key variable seems to be the soil water content [100]. This becomes more evident for the spring and mid-summer, coinciding with the growth phases [83].

Combined Influence of Soil Moisture Together with Climatic Factors
From the results obtained with the different combinations of variables in the MLR analysis (Figures 4 and 5), it was observed that soil moisture was the most important variable, showing the lowest p-values in all cases. In the spring and summer months, median p-values are lower than the significance threshold, and the percentage of significant samples exceeds 50%. In all cases, the minimum p-value was reached in June, ranging from 0.02 to 0.04. These results are in agreement with that obtained in the correlation analysis, and the period when the amount of water available in the soil reaches the most critical relevance for the growth of trees was again identified. of soil moisture in the MLR models. The results obtained from the correlation analysis showed that soil moisture plays a fundamental role in tree growth, and this stood out when considering all environmental factors studied together. These results confirmed the hypothesis raised by [100], that soil water availability better estimates Aleppo pine growth than other variables such as precipitation.  Other studies have incorporated different climatic variables into models, commonly precipitation and temperature, in order to simulate tree growth [65,103]. However, soil moisture has not been used for this purpose. It would be a good option to use soil moisture as an adequate variable to investigate tree growth, as it has been shown from the results obtained with the MLR.

Environmental Modulation of the Role of Soil Moisture in Aleppo Pine Growth
It has been seen so far that soil moisture plays a predominant role in relation to the other environmental variables studied in the present work and those that are commonly used in studies investigating the effects of climate factors on Aleppo pine growth [66]. It is also of interest to investigate whether this relationship follows a spatial or environmental pattern. From the results analyzed, it has been seen that there is variability between Conversely, the climatic variables used did not show a clear influence on Aleppo pine growth in this analysis (Figures 4 and 5). The percentage of significant samples never exceeded 40%, and the median p-values were always higher than the significance threshold. This reveals a strong influence of soil moisture on Aleppo pine growth with respect to the climatic variables, which becomes clearer in spring and summer months. In addition, the p-values obtained with the climatic variables with different combinations in the MLR analyses were very similar. This could be due to the high correlation between precipitation and SPI, on the one hand, and between maximum temperature, minimum temperature, and PET on the other. Although climatic variables exert influence on the growth of Aleppo pine when they are analyzed individually [26,89], when all factors analyzed are considered together, the effect of the climatic variables was masked because of the higher weight of soil moisture in the MLR models. The results obtained from the correlation analysis showed that soil moisture plays a fundamental role in tree growth, and this stood out when considering all environmental factors studied together. These results confirmed the hypothesis raised by [100], that soil water availability better estimates Aleppo pine growth than other variables such as precipitation.
Other studies have incorporated different climatic variables into models, commonly precipitation and temperature, in order to simulate tree growth [65,103]. However, soil moisture has not been used for this purpose. It would be a good option to use soil moisture as an adequate variable to investigate tree growth, as it has been shown from the results obtained with the MLR.

Environmental Modulation of the Role of Soil Moisture in Aleppo Pine Growth
It has been seen so far that soil moisture plays a predominant role in relation to the other environmental variables studied in the present work and those that are commonly used in studies investigating the effects of climate factors on Aleppo pine growth [66]. It is also of interest to investigate whether this relationship follows a spatial or environmental pattern. From the results analyzed, it has been seen that there is variability between samples and that not all places show the same degree of relationship. It is therefore interesting to know whether this relationship between soil moisture and Aleppo pine growth is modulated by the environmental or geographical conditions of the places where it has been analyzed. It is also important to determine whether there is a predominant factor that accentuates or attenuates the role of soil moisture. To investigate this issue, the R (LF vs. TRI) data obtained were correlated with the annual mean temperature, maximum temperature, minimum temperature, precipitation, potential evapotranspiration, P/PET ratio, and altitude, latitude, and longitude, of every location. P/PET is a common aridity index, usually used to characterize the degree of water limitation of an ecosystem or region [104].
The results of this analysis show that almost no relationship exists with the thermal conditions of the areas (Figure 6a-c). Only mean and minimum temperature show a significant and direct relationship with the R values of October of the previous year. This could be related to the secondary peak of TRI in autumn [48], indicating that such a pattern in the relationship between soil moisture and tree growth is more evident where the conditions are warmer in that season. However, the R values between LF and TRI are clearly correlated with environmental factors that express the water availability of the sampling sites. Precipitation and P/PET show an inverse relationship, whereas in the case of the PET ratio, the relationship is direct (Figure 6d-f). This is consistent, as it indicates that the relationship between soil water content and tree growth is closer as water availability is lower at a given location. The relationship is very clearly expressed, as it is statistically significant (p < 0.01) in the spring and summer months, the time of maximum biological activity and aridity under Mediterranean conditions. A similar pattern was observed by [31], who found that the relationship between soil moisture and Aleppo pine growth was stronger in more arid sites.

Conclusions
From the two approaches applied in this study, univariate and multiple linear regression analysis, it was seen that soil moisture is the most important variable in Aleppo pine growth, compared with the climatic factors analyzed. This work found that the influence of those variables, commonly used for this kind of study, is lower than that observed for soil moisture itself. It was also seen how the variables related to water dynamics (precipitation or SPI) showed a clearer relationship than other climate factors with the growth of the trees. Previous studies showed that factors such as temperature conditions had a considerable influence, while this study showed a much less clear relationship than that obtained for soil water content.
When soil moisture was taken into account together with the different climatic variables through multiple linear regression analysis, the results pointed in the same direction. It was clearly observed how soil moisture is the variable that plays the most important role in the growth of Aleppo pine. In fact, soil moisture outperforms the other factors considered, and no significant influence was observed when they were used in the analysis.
From the analysis of the relationship between soil moisture and tree growth, remarkable variability was observed between the studied sites. For that reason, it was determined whether the link between both variables was modulated by the environmental conditions. From the correlation of that relationship with the climatic characteristics and the geographical location of the forests where the samples were obtained, it was determined that thermal conditions have no influence. However, the relationship with the variables that define the availability of water, mainly the amount of precipitation and the aridity index, is very clear, especially at the time of year when, according to the characteristics of the It is noticeable that no relationship was found with altitude ( Figure 6h,i), considering that the variability of cases is high (Table 1), ranging between 170 and 1676 m.a.s.l. It is well known that altitude has great importance in plant phenology [105], especially because it controls and modulates thermal conditions. It is also known that altitude interferes with tree growth, but usually this is because this terrain factor conditions aspects such as competition for light or nutrients [106]. In the present study, the results show that altitude has no relevance in relation with the influence of soil water content on Aleppo pine growth.
However, geographical location seems to be clearly related. Both the latitude and longitude of the sampled sites show a clear, inverse relationship with R ( Figure 6). As the latitude decreases, the value of R is higher, and the link between soil moisture and tree growth is stronger. This is clearer in the months of April to June, where the correlation is statistically significant, and it is precisely in the months when the correlation between soil moisture and growth is higher (Figures 2 and 3). The longitude results are similar, since in this case, both latitude and longitude are related. Owing the configuration of the eastern part of the Iberian Peninsula (Figure 1), as the latitude decreases (to the south), the longitude does as well (to the west). Therefore, the geographical location of the samples becomes a modulating factor for the role of soil moisture. This is in relation to the fact that, in this region, aridity increases to the south and, especially, to the southeast [107]. This result is clearly consistent with what we observed in relation to indicators such as precipitation or P/PET ratio (aridity index). Therefore, again, it is observed that the role of soil moisture in the growth of Aleppo pine is modulated by the aridity or water availability conditions.

Conclusions
From the two approaches applied in this study, univariate and multiple linear regression analysis, it was seen that soil moisture is the most important variable in Aleppo pine growth, compared with the climatic factors analyzed. This work found that the influence of those variables, commonly used for this kind of study, is lower than that observed for soil moisture itself. It was also seen how the variables related to water dynamics (precipitation or SPI) showed a clearer relationship than other climate factors with the growth of the trees. Previous studies showed that factors such as temperature conditions had a considerable influence, while this study showed a much less clear relationship than that obtained for soil water content.
When soil moisture was taken into account together with the different climatic variables through multiple linear regression analysis, the results pointed in the same direction. It was clearly observed how soil moisture is the variable that plays the most important role in the growth of Aleppo pine. In fact, soil moisture outperforms the other factors considered, and no significant influence was observed when they were used in the analysis.
From the analysis of the relationship between soil moisture and tree growth, remarkable variability was observed between the studied sites. For that reason, it was determined whether the link between both variables was modulated by the environmental conditions. From the correlation of that relationship with the climatic characteristics and the geographical location of the forests where the samples were obtained, it was determined that thermal conditions have no influence. However, the relationship with the variables that define the availability of water, mainly the amount of precipitation and the aridity index, is very clear, especially at the time of year when, according to the characteristics of the Mediterranean climate, water scarcity is greater. Therefore, the soil water content is decisive for the growth of Aleppo pine, and this variable is much more determinant as the environmental conditions are more arid. This result is very relevant considering that most climate change projections show an increase in water deficits in regions such as the Mediterranean, where water-limited ecosystems are predominant.
The results obtained in this study have shown that the role of soil moisture is decisive for this tree species, which is characteristic of the vegetation in water-limited environments, and that it is possible to undertake its study using approaches and methodologies, such as satellite or modeled soil moisture, that are fully available for almost any territory and at appropriate temporal and spatial scales. The use of these new approaches allows knowledge of the inherent environmental conditions of forest areas to be advanced. These tools can also be very useful to better understand the functioning of ecosystems as sensitive as forests and to face the challenges of climate change. In addition, recent studies highlighted the potential of Aleppo pine as a suitable climate proxy able to reconstruct past precipitation conditions on semiarid areas of its distribution area [16]. Despite its importance, information on past soil moisture conditions is extremely scarce, and our study demonstrates the high potential of dendroclimatological studies to be used to improve our knowledge of soil moisture evolution over the last few centuries in highly sensitive areas.