Temporal Variations in Temperature and Moisture Soil Proﬁles in a Mediterranean Maquis Forest in Greece

: Soil moisture (SM) and temperature (ST) are critical factors in forest eco-hydrological research. In this study, we investigated the inter-and intra-annual changes in SM and ST proﬁles in a mixed Mediterranean maquis forest stand together with soil and meteorological parameters. Hourly data from three ﬁeld measurements points at four depths ( − 5, − 20, − 40 and − 70 cm) for 6 years were interpolated using the kriging method to produce annual SM and ST proﬁles. The results indicate that air temperature highly affects the upper 5 cm of the mineral soil. In general, it increases with depth in winter at an average rate of 0.036 ◦ C/cm and decreases in summer (0.035 ◦ C/cm), presenting higher values compared to air temperature from April to August and lower ones during the rest of the period. Precipitation is the main factor driving SM variations up to a superﬁcial soil depth of 40 cm. The upper soil layer (0–40 cm) inﬁltrates water faster and presents high SM variability, especially in monthly and seasonal (year to year) time steps. The maquis forest stands are likely to be strongly affected by climate change, therefore the results of this study could be useful in hydrological and climate change studies focused on maquis vegetation water management.


Introduction
One of the most important variables in the hydrological cycle is water storage in the unsaturated zone of soils.Many rainfall-runoff and climatic models applying different spatial scales require soil moisture (SM) as an input parameter [1,2].The study of temporal variations in SM in the unsaturated zone of forest soils is very complex [3,4].The factors affecting water transport into the deeper layers have been well documented by many researchers in previous studies.SM is critically affected by the type and structure of forest vegetation [5][6][7], which determine the patterns of intercepted precipitation [8][9][10][11].Vegetation also provides a diverse hydrologic response at the catchment scale [12].The temporal variations of SM in different soil layers are important for assessing the hydrological response of the catchment [13] since they are directly related to water storage capacity and evaporation rates [14].The need for SM measurements varies depending on the scope of the hydrological assessment required by each study.Point-specific field measurements (in situ) using sensors or ground networks are characterised by high accuracy and allow continuous monitoring [15,16].On the other hand, remote sensing monitoring products have a lower accuracy but a high spatial resolution of more than 250 m [17]; for this reason they can provide important information for large areas.
Studying only SM dynamics may limit our understanding of the field hydrological processes [16], as these are influenced by several factors.According to Kumar et al. [18], monitoring of both SM and soil temperature (ST) is an essential input for hydrological simulations and predictions.The most important factors affecting the values of ST are solar radiation and wind speed [19,20].Together with soil thermal parameters which are responsible for heat transfer into the soil, such as heat capacity and thermal conductivity, they usually form heterogeneous conditions for heat transfer at different positions in the soil profile [21].This leads to spatial variability in ST.ST in the upper layers (0-10 cm) can explain up to 95% of the temporal changes in soil respiration [22], while the variations in ST are influenced by many different factors, such as solar radiation, precipitation, vegetation, surface air temperature [19,20] and microbial activity in the soil [23].Therefore, studying SM together with ST and other hydrophysical properties of soils (texture, hydraulic conductivity, root density, etc.) and assessing their complex interactions can enhance our understanding of soil-water-plant relationships [12,24,25].
Air temperature is important in controlling vegetation respiration rates [26,27], soil water and nutrient availability, seed germination, plant growth and nutrient enrichment [28].It is associated with ST through complex processes and several factors [29,30].Precipitation greatly affects soil water availability and SM in rainfed ecosystems.Its volumes and distribution differ greatly, both spatially and temporally, especially in the Mediterranean [24,[31][32][33][34][35][36][37].Summer water deficits and high temperatures as well as winter low temperatures are identified as the main factors affecting growth in forests [38,39].Considering that the water stress period for vegetation is increasing due to climate change [32,40,41], forest adaptation measures are imperative [42,43]; thus, the realistic representation of water availability in soils is critical [44] for the management of natural ecosystems.
The long-term investigation of SM and ST temporal variations, especially on an annual cycle detecting the highs and lows during specific seasons or periods, can produce important results concerning the water stress of plants and its impacts on vegetation's survival, growth and phenology, especially at local scales.Vereecken et al. [16], reported that spatio-temporal SM field observations would be useful to determine soil water balances if investigated in conjunction with other hydrological processes.Dinca et al. [45] identified high variability in time and depth for soil water content under different climatic and vegetation type conditions, reporting that SM values at different depths can be explained by four elements, i.e., vegetation type, soil type, precipitation and temperature.
The critical role of vegetation was also addressed by Vergarechea et al. [4], who suggested the use of mixed stands as more resilient to face climate change.The aim of this study is to investigate the temporal variations in SM and ST profiles at 0-70 cm depths in a mixed sclerophyllous evergreen broadleaved maquis forest experimental stand in Greece.Understanding the temporal and vertical changes in these abiotic factors and their interactions will be useful for eco-hydrological modelling in similar Mediterranean environments, especially in the case of the complex canopy architecture of a maquis forest consisting of different vegetation floors (dominant, sub-floor, herbaceous and litter) and for a long (6-year) period.Considering that the Symposium of International Association of Hydrological Sciences (IAHS) in 2018 identified 23 "unsolved problems" in hydrology, ref.
[3] and among them was the necessity to determine the hydrologic laws that change with scale, this study may contribute to this direction since it applies field measurements in small temporal and spatial (forest stand) scales.

Study Area
The study was carried out in a forest plot (45 m × 33 m) located in western Greece (38 • 50 46 N, 21 • 18 16 E, alt: 340 m a.s.l.) as shown in Figure 1.The dominant tree species of the plot is Quercus ilex, and the sub-floor hosts mainly Arbutus unedo, Phillyrea latifolia and Erica arborea shrub species.The herbaceous floor and its composition with hygrophyte and mesohygrophyte species (Smilax aspera, Asplenium trichomanes, Asplenium onopteris, Asplenium adiantum nigrum, Carex distachya, Selaginella denticulate) and moss (Scleropodium purum) confirm the exceedance of 1000 mm annual rainfall in the area.This floor is not particularly dense but becoming denser in the gaps of shrubs and on their edges; its coverage can seasonally (mainly in spring) exceed 90% of the ground.In general, the evergreen broadleaved species are dense, presenting a LAI = 5.8 and height ranging from 6 to 15 m.The parental soil material is flysch and the soil type is Eutric Cambisol [46].
particularly dense but becoming denser in the gaps of shrubs and on their edges; its coverage can seasonally (mainly in spring) exceed 90% of the ground.In general, the evergreen broadleaved species are dense, presenting a LAI = 5.8 and height ranging from 6 to 15 m.The parental soil material is flysch and the soil type is Eutric Cambisol [46].The climate of the area is warm temperate with hot and dry summers, classified in type Csa according to the modified Köppen-Geiger world map classification [47].Based on UNEP s [48] aridity climate classification, the climate of the broader area is humid [49] presenting a Thornthwaite s [50] aridity index value of 0.78 during the recent climatic period .This is significantly decreased compared to the previous period  when its value was 0.85, indicating a rapid local climate transition to more arid conditions [41].According to the data from the meteorological station installed in a forest opening near the plot (420 m) covering the period from 1973 to 2020, the average annual air temperature is 15.4 ± 0.64 °C, varying seasonally from 7.4 ± 1.37 °C in winter to 23.9 ± 1.38 °C in summer.The precipitation has an average annual value of 1049 ± 225.5 mm and presents what is expected for the Mediterranean climate seasonal variation: unevenly distributed between seasons with higher amounts falling in winter (40.2%) and autumn (33.3%).The pluviothermic diagram in Figure 2 confirms a water deficit period of about of 3 months starting from June until the beginning of September.The climate of the area is warm temperate with hot and dry summers, classified in type Csa according to the modified Köppen-Geiger world map classification [47].Based on UNEP's [48] aridity climate classification, the climate of the broader area is humid [49] presenting a Thornthwaite's [50] aridity index value of 0.78 during the recent climatic period .This is significantly decreased compared to the previous period  when its value was 0.85, indicating a rapid local climate transition to more arid conditions [41].According to the data from the meteorological station installed in a forest opening near the plot (420 m) covering the period from 1973 to 2020, the average annual air temperature is 15.4 ± 0.64 • C, varying seasonally from 7.4 ± 1.

Air Temperature and Precipitation Data
Air temperature was recorded with a time step of 15 min for a 6-year period from 1 January 2013 to 31 December 2018 inside the meteorological station using an MP101A-T7-WAW sensor (Rotronic Measurement Solutions).The monitoring of the air temperature was continuous, and there were gaps only for 68 days due to malfunction of the sensor (3.1%).There was no gap filling, and the gaps were excluded from the analysis.
Bulk precipitation (PR) was measured by a Belford rain gauge installed in the abovementioned forest opening with its orifice at a height of 1.5 m above ground.Its recording tapes can continuously measure PR up to a maximum of 300 mm.The intensity of each rainfall event was analysed in an hourly time step from the readings of the recording tapes.Individual rainfall events were defined assuming that there was at least a six-hour interval with zero rain between them.The volume of measured PR (mm) was also compared with the average volume collected from three volumetric gauges installed next to the Belford rain gauge at the same height.
Two different flows are essential for the investigation of SM temporal variation: throughfall precipitation, which is the amount of PR water that penetrates the forest canopy, and stemflow, which is the flow from PR that drains down the trunks of trees.The sum of these flows constitutes the effective rainfall (PRef) that reaches the ground and is responsible for the formation of SM.
For the measurement of throughfall, a 25-volumetric collector network installed under the forest canopy (Figure 1), following the reports of Rodrigo and Avila [51] and Llorens and Domingo [52], was used.The collectors orifices were at 1.30 m above soil surface, and they were vertically adjusted regardless of the topographical slope.Effective rainfall (PRef) daily values were estimated as the sum of throughfall and stemflow, assuming that throughfall is produced in rainfall events with PR > 4 mm, except for the cases where the small rainfall events occurred on consecutive days, as suggested by Baloutsos et al. [53,54] for the same site.The authors [53,54]  Air temperature was recorded with a time step of 15 min for a 6-year period from 1 January 2013 to 31 December 2018 inside the meteorological station using an MP101A-T7-WAW sensor (Rotronic Measurement Solutions).The monitoring of the air temperature was continuous, and there were gaps only for 68 days due to malfunction of the sensor (3.1%).There was no gap filling, and the gaps were excluded from the analysis.
Bulk precipitation (PR) was measured by a Belford rain gauge installed in the abovementioned forest opening with its orifice at a height of 1.5 m above ground.Its recording tapes can continuously measure PR up to a maximum of 300 mm.The intensity of each rainfall event was analysed in an hourly time step from the readings of the recording tapes.Individual rainfall events were defined assuming that there was at least a six-hour interval with zero rain between them.The volume of measured PR (mm) was also compared with the average volume collected from three volumetric gauges installed next to the Belford rain gauge at the same height.
Two different flows are essential for the investigation of SM temporal variation: throughfall precipitation, which is the amount of PR water that penetrates the forest canopy, and stemflow, which is the flow from PR that drains down the trunks of trees.The sum of these flows constitutes the effective rainfall (PR ef ) that reaches the ground and is responsible for the formation of SM.
For the measurement of throughfall, a 25-volumetric collector network installed under the forest canopy (Figure 1), following the reports of Rodrigo and Avila [51] and Llorens and Domingo [52], was used.The collectors' orifices were at 1.30 m above soil surface, and they were vertically adjusted regardless of the topographical slope.Effective rainfall (PR ef ) daily values were estimated as the sum of throughfall and stemflow, assuming that throughfall is produced in rainfall events with PR > 4 mm, except for the cases where the small rainfall events occurred on consecutive days, as suggested by Baloutsos et al. [53,54] for the same site.The authors [53,54] also proposed the following equation for the estimation of daily stemflow for the rainfall events with PR > 10 mm: Stemflow = 0.0934 PR − 0.7299 (March to August) 0.1064 PR − 1.2416 (September to February) (1) 2.2.2.Soil Moisture and Soil Temperature Data SM and ST data were recorded every hour in three representative plot replicates (Figure 1) at four soil depths (−5, −20, −40, −70 cm) using ECH2O EC-5TM water content and temperature sensors (Decagon Devices Inc., Pullman, WA, USA), with an accuracy of ±1 • C for ST and ±0.03 m 3 /m 3 or ±3% VWC for SM, using the Topp equation which is typical in mineral soils with solution electrical conductivity lower than 10 dS/m.For this study, no additional sensor calibration was performed considering adequate data quality, since for all soil samples the electrical conductivity was lower than the above mentioned threshold.However, it should be mentioned that the vicinity of the sensors to roots and stones, as well as the carbon content of the soil, may have an impact on the quality of measurements, though this impact was considered to be systematic.In this work, SM and ST are used as indicative parameters to assess vegetation water stress.The triplicate design for the network of sensors adopted in this study is considered appropriate for monitoring the variability of SM and ST at different depths and seasons [55].After the installation of the sensors, we allowed a four-month wet period (September to December 2012) for the stabilization of the measurements and to identify any quality issues.The daily values of SM and ST were calculated as averages from the hourly data.In this work, the annual SM and ST profile patterns were also employed, and they were produced by the application of SURFER 13.6.618software (Golden Software, LLC, 2016) using the kriging interpolating technique which estimates the unknown values at all points of a defined spatial domain using a weighted average of all known values around the unknown point [56].

Soil Sampling and Analysis
The soil sample collection was carried out by means of systematic sampling, following the guidelines suggested by the soil expert panel of the ICP Forests (http://icp-forests.net) (accessed on 13 February 2023) into which our experimental site is integrated and in the framework of project "Biosoil".Inside the plot, along three 25 m long transects 5 m apart, six soil pits were excavated 5 m apart from each other on each transect.The organic horizons, comprising litter (L) and the fragmented horizon (FH) layers, were collected from each pit by means of a metallic frame having an area of 15 × 15 cm.Mineral soil layers were collected with a shovel from depths 0-10 cm, 10-20 cm, 20-40 cm and 40-80 cm.From a total of 18 soil samples, six (approximately 2 kg each) per horizon and soil depth were mixed to have three pooled samples per horizon and depth.The bulk density of the mineral soils was measured with a cylinder of known volume (129 cm 3 ).The value was converted to bulk density of fine earth (<2 mm) after subtracting the volume of coarse material (stones and gravel).The percentage of large stones was estimated visually in the field.All soil samples were transferred to the laboratory and were air-dried.Apart from the L layer, the other samples passed through a 2 mm sieve.After the analysis, the samples were dried at 105 • C for 24 h, and the results were expressed in dry weight of soil at that temperature.
The soil-water retention curve (pF) and hydraulic conductivity function K(θ) were produced for the assessment of the soil hydraulic characteristics.Soil samples were taken during the installation of the SM and ST probes at the three measurement points (Figure 1) at depths of −20, −40 and −70 cm.Two samples were extracted per depth, resulting in a total number of 18 samples that were analysed in the laboratory.In total, 7 points of the pF curve were estimated from laboratory measurements using the method of the sand-kaolin box (following the protocol ISO 11274) for pressure heads down to −330 hPa and 5 more points for even smaller pressure heads using a pressure membrane apparatus (down to −5000 hPa).The 12 experimental pressure points were calibrated and verified by using the RETC software version 6.02 (©2005-2009) for the pF curve fitting.The remaining 6 samples were used for the results cross-validation.
To determine the hydraulic conductivity in saturation (K s ), the double-ring infiltrometer method was employed and 9 in situ experiments (3 on each measurement point) were conducted.The K s was found to have an average value of 4.5 cm/h which is considered moderately high [57].This value was imported in RETC for the prediction of the K(θ) curve along with the Van Genucthen-Mualem [58] parameters.
The texture analysis was completed using the hydrometer method [59].The pH was measured electrometrically (pH meter ORION 410-A) using a glass electrode in water and 0.01 M CaCl 2 (1:5 soil weight per volume of water and CaCl 2 solution).The concentrations of total N and organic C in soils were measured with a CN analyser (Vario MAX) [60].

Soil Texture, Nutrient Status and Soil Hydraulic Properties
The soil is silty loam, and its reaction is slightly acidic (Table 1).The coefficients of variation (CV-defined as percentage of standard deviation over the mean) given in the parentheses of Table 1 show a small spatial variation (<15%) for all soil properties [61].Moderate spatial variation exists only in 20-40 cm (21.56%) and 40-80 cm (21.64%) layers for the sand texture fraction, which does not affect the texture class (silty loam) at any measurement point [62].According to B.A.I. [63], the nitrogen and organic carbon concentrations are satisfactory for plant demands.The ratio C/N is rather high in the litter (L) layer, but it goes down abruptly in the FH and mineral layers.This means that the decomposition rates are satisfactory.Michopoulos et al. [64] found that the soil organic matter in the maquis plot under consideration had a higher decomposition rate than that in a mountain fir forest in central Greece despite the fact that the fir forest had higher SM but lower ST.This suggests that ST is more important than SM regarding microbial activity when comparing different ecosystems.Soil hydraulic behaviour is important for the assessment of soil-plant-atmosphere relations.The soil water availability and the water retention in soil profile are critical for plants and for the design and implementation of adaptation strategies to mitigate the impacts of climate change, especially in natural ecosystems.The soil water retention pF (fitted and extrapolated), the hydraulic conductivity K(θ) (predicted) curves derived from RETC software (R 2 = 0.995) and the main hydraulic characteristics are presented in Figure 3. Important values for the pF curve are SMs in field capacity (θ FC = 0.258 cm 3 /cm 3 , pF = 2.525), in wilting point (θ WP = 0.111 cm 3 /cm 3 , pF = 4.18) and in saturation (θ s = 0.482 cm 3 /cm 3 ).According to these values, the soil profile down to −70 cm can potentially hold 337.0 mm of water in saturation conditions, which is much greater than the amount of water available for vegetation (102.9 mm).It should also be stated that the fitted residual water content (θ r ) is in line with the findings of Dexter and Richard [65] who reported that the amount of adsorbed water on soil particle surfaces is zero at pF = 6.6.The resulting Van Genuchten parameters from the application of RETC, for the pF curve and their standard error, are a = 0.0260 ± 0.0023, n = 1.2777 ± 0.0105 and m = 0.2174.The results of the prediction of the K(θ) curve (green line) show, as expected, that K is decreasing very abruptly from K s (4.5 cm/h) to the value of K = 0.05 cm/h (θ = 0.400 cm 3 /cm 3 ) which results in extreme small velocities.

Temporal Variations of Air Temperature
The annual average air temperature (Tav) for the 6-year study (15.Results indicate that all monthly air temperatures of the years 2013 and 2015 are within the limit value chosen by the authors of "mean ± 1SD" from the climatic means, suggesting normal meteorological conditions in terms of temperature.Similarly, the year 2014 had the expected monthly temperatures for all months, except for February, when slightly warmer conditions (by +2.2 °C) prevailed.In 2016, February and April were also

Temporal Variations of Air Temperature
The annual average air temperature (T av ) for the 6-year study (15.

Temporal Variations of Air Temperature
The annual average air temperature (Tav) for the 6-year study (15.Results indicate that all monthly air temperatures of the years 2013 and 2015 are within the limit value chosen by the authors of "mean ± 1SD" from the climatic means, suggesting normal meteorological conditions in terms of temperature.Similarly, the year 2014 had the expected monthly temperatures for all months, except for February, when slightly warmer conditions (by +2.2 °C) prevailed.In 2016, February and April were also Results indicate that all monthly air temperatures of the years 2013 and 2015 are within the limit value chosen by the authors of "mean ± 1SD" from the climatic means, suggesting normal meteorological conditions in terms of temperature.Similarly, the year 2014 had the expected monthly temperatures for all months, except for February, when slightly warmer conditions (by +2.2 • C) prevailed.In 2016, February and April were also warmer (by +3.3 and +2.8 • C, respectively), whereas in 2017, January temperatures were cooler (by −2.2 • C).
Finally, in 2018, April, May and November temperatures have exceeded the climatic normal by +4.3, +2.4 and 2.0 • C, respectively.The above pattern indicates that during our study period the temperatures of the environment were rather normal for the area and in some cases warmer.
The general intra-annual pattern of temperature is coherent with the typical distribution of the Mediterranean climate presenting seasonal variability: higher values in summer (+24.1 ± 0.84 • C), lower in winter (+7.9 ± 1.23 • C) and intermediate values of +14.3 ± 1.42 • C and +16.6 ± 0.93 • C for the transitional seasons of spring and autumn, respectively (Table 2).The warmest month is August (25.7 ± 0.66 • C) and the coolest is January (6.9 ± 1.18 • C).

Precipitation Variations and Characteristics
The total rainfall events that occurred in the 6-year period of the study were 531.Annual PR ranged from 1090.4 mm in 2017 to 1525.0 mm in 2014, suggesting that all years of the study are wetter compared to the long-term annual PR average of the 1973-2020 period (1049 mm) with percentages ranging from +4% in 2017 to +45% in 2014.The monthly PR values during the study period (2013-2018) showed that November was the wettest month of the year (221 mm) and August the driest one (6 mm).In general, the annual PR yield mainly occurred in winter (549 mm) and autumn (403 mm), and to a lower degree, in spring (285 mm) and summer (63 mm).
The monthly PR deviations from the respective long-term means are shown in Figure 5, along with the SD of the monthly PR of the 1973-2020 climatic period.In 2013, January (+128.3%),February (+95.0%) and November (+82.9%) received higher and December (−65.7%)lower amounts of PR since their respective PR values are out of the threshold of "mean PR ± 1SD".In 2014, only wetter months were recorded, and these were April (+113.1%),May (+344.6%)and October (+82.7%),whereas in 2015, September (+112.5%)and November (+105.5%)were wetter and April (−75.2%) and December (no rain) drier.In 2016, January (+87.2%),September (+231.1%)and March (+109.1%)had a higher PR and December (no rain) a lower one.Finally, December 2017 (+107.9%)and February (+186.3%),March (+94.0%)and June (+272.3%)2018 were wetter, while March (−67.7%)2017 and April (−82.2) 2018 were drier.An increase in rainfall variability is one of the key findings of all anthropogenic scenarios in global climate model projections [66,67].Our results indicate a high monthly variability in PR which is consistent with the findings of other studies implemented in the Mediterranean [35,68].The above mentioned percentages refer to the percentage differences from the respective monthly climatic PR of the period 1973-2020.To investigate the rainfall characteristics, the 531 rainfall events that occurred during the study period were grouped into categories based on their mean intensity and duration (Figure 6).The great majority (87.8%) of the incidents was of low intensity (<4 mm/h).These are associated with 74.4% of the total PR, and their low intensity enhances the ability of the soil to store water in the root zone and/or to infiltrate water deeper.The extreme intense rainfall events are rather scarce, considering that in our 6-year study only 13 incidents were recorded with an average intensity greater than 8 mm/h and only one of these had an intensity higher than 20 mm/h.The 13 relatively high-intensity rainfall events accumulated only 5.2% of the incoming total PR in the ecosystem, and 10 of them had small durations (<4 h) that could result in risks to the forest soil matrix (erosion or increased runoff).Regarding event duration, in most cases (62%) rainfall lasts more than 4 h.For the rest of the events (duration with less than 4 h), it should be mentioned that, in general, they provided 15.4% of the total PR to the ecosystem and in most cases are characterized by low mean intensities (83.5% on average).The effective part of precipitation (PRef) also presents similar variability.
The average monthly PRef for the study period (Table 2) ranges from 3.5 ± 3.9 mm (49.3% of PR) in July to 176.0 ± 76.2 mm (79.8% of PR) in November.It is also perceived that the wettest month is November followed by January, while summer months August  To investigate the rainfall characteristics, the 531 rainfall events that occurred during the study period were grouped into categories based on their mean intensity and duration (Figure 6).The great majority (87.8%) of the incidents was of low intensity (<4 mm/h).These are associated with 74.4% of the total PR, and their low intensity enhances the ability of the soil to store water in the root zone and/or to infiltrate water deeper.The extreme intense rainfall events are rather scarce, considering that in our 6-year study only 13 incidents were recorded with an average intensity greater than 8 mm/h and only one of these had an intensity higher than 20 mm/h.The 13 relatively high-intensity rainfall events accumulated only 5.2% of the incoming total PR in the ecosystem, and 10 of them had small durations (<4 h) that could result in risks to the forest soil matrix (erosion or increased runoff).To investigate the rainfall characteristics, the 531 rainfall events that occurred during the study period were grouped into categories based on their mean intensity and duration (Figure 6).The great majority (87.8%) of the incidents was of low intensity (<4 mm/h).These are associated with 74.4% of the total PR, and their low intensity enhances the ability of the soil to store water in the root zone and/or to infiltrate water deeper.The extreme intense rainfall events are rather scarce, considering that in our 6-year study only 13 incidents were recorded with an average intensity greater than 8 mm/h and only one of these had an intensity higher than 20 mm/h.The 13 relatively high-intensity rainfall events accumulated only 5.2% of the incoming total PR in the ecosystem, and 10 of them had small durations (<4 h) that could result in risks to the forest soil matrix (erosion or increased runoff).Regarding event duration, in most cases (62%) rainfall lasts more than 4 h.For the rest of the events (duration with less than 4 h), it should be mentioned that, in general, they provided 15.4% of the total PR to the ecosystem and in most cases are characterized by low mean intensities (83.5% on average).The effective part of precipitation (PRef) also presents similar variability.
The average monthly PRef for the study period (Table 2) ranges from 3.5 ± 3.9 mm (49.3% of PR) in July to 176.0 ± 76.2 mm (79.8% of PR) in November.It is also perceived that the wettest month is November followed by January, while summer months August  Regarding event duration, in most cases (62%) rainfall lasts more than 4 h.For the rest of the events (duration with less than 4 h), it should be mentioned that, in general, they provided 15.4% of the total PR to the ecosystem and in most cases are characterized by low mean intensities (83.5% on average).The effective part of precipitation (PR ef ) also presents similar variability.
The average monthly PR ef for the study period (Table 2) ranges from 3.5 ± 3.9 mm (49.3% of PR) in July to 176.0 ± 76.2 mm (79.8% of PR) in November.It is also perceived that the wettest month is November followed by January, while summer months August (6.0mm) and July (7.1 mm) are the driest ones.Both monthly SD and coefficients of variation (CV) presented similar variations during the year in both PR and PR ef .It is worth noting that extreme values of CV for August are due to the occurrence of only one rainfall event in 2017 (3.8 mm PR) and four rainfall events in 2018 (32.4 mm PR), while the other four years of observation no rainfall events occurred.In general, the annual PR ef ranged in our study from 843.8 mm in 2017 to 1046.8 mm in 2014, representing 77.4% and 68.6% of PR, respectively.Similar percentages were also presented by other studies.Baloutsos et al. [69] measured PR ef values of 68.0% ± 3.2% in the same forest plot for a seven-year period.Llorens et al. [10] measured a bulk interception of 24% of PR in a Pinus sylvestris forest patch in a Mediterranean mountainous area, which indicates that 76% is attributed to PR ef .Avila and Rodrigo [70] estimated in two Quercus ilex stands in Spain PR ef percentages of 79.1% and 77.5%, respectively.Vicente et al. [36] simulated PR ef values in four mixed Mediterranean forests in Spain which ranged from 56.8% to 75.9% of PR.
In the forest environment, the complex interactions between the canopy and rainfall characteristics determine the throughfall inputs [71] and therefore the interception rates [8].As rainwater drips through the forest canopy, the droplets' diameter, velocity and kinetic energy change [8], thus introducing higher variability of the water distribution on the forest soil surface and highly affecting the soil's moisture and temperature profiles as well as the soil's infiltration properties, both spatially and temporally.In our study site, the kinetic energy of PR ef (mainly THR) is expected to be small due to the rather dense canopy (LAI = 5.84) in conjunction with the rainfall events that are mainly of high duration and low intensity.This enhances the ecosystem's ability to infiltrate water with higher rates and limits the possibility for surface runoff formation, even if hydraulic conductivity is relatively low.

Temporal Variations of Soil Temperature
Annually, ST showed a mean value in −5 cm of 14.5 ± 1.64 • C, in −20 cm of 14.4 ± 1.30 • C, in −40 cm of 14.3 ± 0.99 • C and in −70 cm of 14.7 ± 0.87 • C. The general pattern of ST at all depths depicts the same distribution with air temperature (see Table 2, T av ), also presenting a seasonal variability with higher values in summer, lower in winter and intermediate values for the transitional seasons of spring and autumn (Table 3).Results from previous studies showed that ST memory increases with soil depth and varies seasonally [30,72].Soil temperature memory varies with season, and deep soil layers have more significant seasonal differences than subsurface soil.The CV index showed a downward trend with depth but with a steeper decline up to the −40 cm.The ST varies inside the soil profile, presenting values with differences among the various depths that range from 0.4 ± 0.96 • C in September to 3.2 ± 0.87 • C in December.With respect to CV and depth, results show larger values in winter months (18.6% in January, 14.1% in December and 13.6% in February).Both SD and CV decrease their values with depth during all months; the warmest month is August with an ST range with depth 2.5 ± 0.60 • C, while the coolest one is January (range 2.6 ± 1.52 • C).
Figure 7a illustrates the mean daily time of ST at four different depths (−5 cm, −20 cm, −40 cm and −70 cm) along with the mean daily T av for the study period.Compared to T av , the soil profile is warmer in the lower depths at the beginning of the year (January, February).It is notable that the diurnal range of ST is highly variable with depth, presenting high diurnal changes at the topsoil zone that rapidly decrease with depth, as is clearly depicted in Figure 7b.The rather stable ST throughout the 24 h day period is mainly observed for the −40 cm and −70 cm soil zones regardless of the season, whereas at the upper soil layers (−5 cm and −20 cm) the diurnal ST ranges maximize in summer.
in September to 3.2 ± 0.87 °C in December.With respect to CV and depth, results show larger values in winter months (18.6% in January, 14.1% in December and 13.6% in February).Both SD and CV decrease their values with depth during all months; the warmest month is August with an ST range with depth 2.5 ± 0.60 °C, while the coolest one is January (range 2.6 ± 1.52 °C).
Figure 7a illustrates the mean daily time of ST at four different depths (−5 cm, −20 cm, −40 cm and −70 cm) along with the mean daily Tav for the study period.Compared to Tav, the soil profile is warmer in the lower depths at the beginning of the year (January, February).It is notable that the diurnal range of ST is highly variable with depth, presenting high diurnal changes at the topsoil zone that rapidly decrease with depth, as is clearly depicted in Figure 7b.The rather stable ST throughout the 24 h day period is mainly observed for the −40 cm and −70 cm soil zones regardless of the season, whereas at the upper soil layers (−5 cm and −20 cm) the diurnal ST ranges maximize in summer.In March, the temperature differences between Tav and the four sensors (at four depths) decrease and tend to equalize; hence, the soil profile shows the same ST throughout its depth which is similar with Tav.However, for a few days (during April), the temperature profile reverses, presenting higher values in the surface that decrease with depth.During this period, the differences between Tav and STs become positive and maximized in August, decreasing thereafter until the end of September when all STs and Tav become equal.By October, the ST profile reverses again and for the next months until the end of February the soil is warmer when compared with the air above, and ST increases with depth in the soil profile.In general, the decreasing rate of ST with depth in winter is about 0.036 °C/cm, whereas in summer the respective rate is of the same magnitude (0.035 °C/cm) but decreasing.
August is the month when the smallest SD and CV values of ST in −70 cm appeared.This strengthens the conclusion that the lack of precipitation and the high temperatures for this month stabilize ST values in the soil profile [72].Thus, until the beginning of September, ST shows its greatest values, as expected, at −5 cm and then at −20 cm, −40 cm and −70 cm.In September and until the end of the year, the shape of ST in the soil profile shows the shape it had in the first two months of the year.Regarding the fluctuations of ST with soil depth, results show high values at the depth of −5 cm and gradually smoothing out  In March, the temperature differences between T av and the four sensors (at four depths) decrease and tend to equalize; hence, the soil profile shows the same ST throughout its depth which is similar with T av .However, for a few days (during April), the temperature profile reverses, presenting higher values in the surface that decrease with depth.During this period, the differences between T av and STs become positive and maximized in August, decreasing thereafter until the end of September when all STs and T av become equal.By October, the ST profile reverses again and for the next months until the end of February the soil is warmer when compared with the air above, and ST increases with depth in the soil profile.In general, the decreasing rate of ST with depth in winter is about 0.036 • C/cm, whereas in summer the respective rate is of the same magnitude (0.035 • C/cm) but decreasing.
August is the month when the smallest SD and CV values of ST in −70 cm appeared.This strengthens the conclusion that the lack of precipitation and the high temperatures for this month stabilize ST values in the soil profile [72].Thus, until the beginning of September, ST shows its greatest values, as expected, at −5 cm and then at −20 cm, −40 cm and −70 cm.In September and until the end of the year, the shape of ST in the soil profile shows the shape it had in the first two months of the year.Regarding the fluctuations of ST with soil depth, results show high values at the depth of −5 cm and gradually smoothing out from −20 cm to −70 cm where ST fluctuations are very small.The results clearly show that ST is affected by T av at different depths, but less so as depth increases.This is due to the stronger persistence effects that prevail in the deeper soil layers originating from previous conditions.According to Bilgili [73], these ST reversals are very important, especially during the winter period, as they generate soil biological activities.Zhang et al. [30] found a significant impact of air temperature on the ST values of the shallower soil layers.Finally, compared to daily T av , the whole soil profile has a higher temperature from the beginning of the year until March and from October until the end of the year.
The simulation in Figure 8 confirms the temporal variation in ST.The values of ST at least at the depth of −5 cm were clearly affected more intensely by air temperature.However, the rates of change both on a daily and annual basis are noticeably smaller compared to T av .The simulation also confirms the finding that the largest ST variation occurs up to the depth of about 40 cm.According to Terradas and Save [39], transpiration rates showed a linear significant decline with decreasing soil temperature, and Lloyd and Taylor [27] found that the effective activation energy for respiration varies inversely with Tav.Combining these findings with our results we can conclude that respiration rates are higher in the upper layers, independently from the root mass.Nevertheless, this is a case for further The simulation also confirms the finding that the largest ST variation occurs up to the depth of about 40 cm.According to Terradas and Save [39], transpiration rates showed a linear significant decline with decreasing soil temperature, and Lloyd and Taylor [27] found that the effective activation energy for respiration varies inversely with T av .Combining these findings with our results we can conclude that respiration rates are higher in the upper layers, independently from the root mass.Nevertheless, this is a case for further investigation.Even though the shading of the soil due to the presence of vegetation is strong (LAI = 5.84), at least in the first 20 cm the temperature response is significantly affected.
Until the end of March, the temperature of the soil profile does not exceed 12.0 • C. The transitional period of the quarter April-June presents many fluctuations in the soil temperature values, both in relation to depth and over time.Indicatively, year 2013 shows from the beginning of April until the middle of May a higher temperature by at least 1.0 • C in the first 25-30 cm compared to the lowest depths.Meanwhile, the year 2016 shows from the beginning of April until the middle of June a gradual increase in ST with approximately parallel lines for the entire soil profile.Despite all the different daily fluctuations, the temperature of the soil profile does not exceed, at the end of June in any year, 20.0 • C, except in the upper layers (up to −20 cm) in 2013, 2014, 2016 and 2017, presenting within this quarter an increase of approximately 8.0 • C.
Within the next quarter (July-September), for all years, the temperature rises to 24.0-25.0• C, showing again the pattern of the first months of the year, i.e., only up to the depth of the first 40 cm.ST in the last 30 cm depth exceeded 20 • C in all years only for the very hot period from the beginning of August until the first days of September.By the end of September, the ST of the whole profile decreases to 18.0 • C with different rates at different depths for all years.
In the last quarter of each year, a gradual decrease in ST is observed throughout the profile.The strong variation in all years that appears from mid-November to mid-March, in the first 30-40 cm is mainly due to the occurrence of many PR rainfall events and the significant and rapid temporal storage of water in the soil.The same strong variation appears during the summer period when the high values of air temperature influence the ST values.

Temporal Variations in Soil Moisture
The results for the whole study period in monthly values are presented in Table 4. Overall, annual SM at −5 cm had a mean value 0.185 ± 0.203 cm 3 /cm 3 ; at −20 cm 0.174 ± 0.205 cm 3 /cm 3 ; at −40 cm 0.208 ± 0.114 cm 3 /cm 3 ; and at −70 cm 0.221± 0.119 cm 3 /cm 3 .Seasonal variability was also observed with higher values in winter (0.247 ± 0.020), lower ones in summer (0.154 ± 0.024) and intermediate values in spring (0.229 ± 0.026) and autumn (0.174 ± 0.049) (Table 4).Spring is the crucial season during which water needs are high for tree growth and determines the duration and the severity of water stress [37].According to Barbeta et al. [31], maquis vegetation with dominant species such as Quercus ilex, Arbutus unedo and Phillyrea latifolia form a dimorphic root system, which has the ability to differentiate soil water uptake, depending on the soil moisture content.Therefore, vegetation adapts to the prevailing meteorological conditions, thus shaping water absorbing conditions at different depths.This could be the reason that explains decreasing SM values (due to hydraulic decrease) during the transitional spring season.Additionally, it could be attributed to the fact that quick drainage may promote aeration in surface soil, hence optimizing the occurred conditions in the upper layers' roots [74].
The monthly CV displays, according to depth (0-70 cm), a range from 48.4% in September and 31.3% in October to 4.7% in December and 4.8% in February.This is because the dry soil profile is moistened with the water of autumn's rainfall events but at rates inextricably linked to their intensities which in our study are bigger during the autumn months (2.6 and 2.3 mm/h, respectively) than in the winter (both 2.0 mm/h) ones.Overall, the smallest values of CV at all depths are recorded in December for −5 cm (10.9%), in March for −20 cm (6.6%) and −40 cm (6.2%), and in January for −70 cm (6.2%).To the contrary, the highest CV values are recorded in the autumn months.This behaviour could be attributed to the fact that during the winter months the soil moisture has higher values due to many rainfall events, while during the autumn months, even though it rains, the infiltrated water mainly moistens the soil profile because this season follows the dry summer period [30].Additionally, it has also been established from previous studies that the existence of more roots in the upper layers of the soil increase and decrease the soil's moisture more intensely, while the lower layers are usually wetter, showing smaller variations in their soil moisture values [43,75].
During the period of June to early of September, vegetation is in water stress due to low PR values (<20 mm), higher ST values (>18 • C) and solar radiation, and therefore, gradually increased values of evapotranspiration.This situation leads of course to an increase in respiration rates, mainly due to an increase in decomposers' respiration activity depending on their composition, litterfall amounts and quality [26,76].This is very important because it is strongly establishing the relationship between carbon and nitrogen contents [77,78].Hence, transpiration needs for tree vitality are served by water availability at lower depths, and therefore, SM values continuously decrease at all depths.
The seasonal SM variation indicated that, in general, drier conditions persisted in winter 2015 and 2018 and wetter ones in 2013, whereas spring of 2017 was the driest season and 2014 the wettest.During the critical summer periods of 2015, 2016, 2017 and 2018, dry conditions persisted, but the wettest summer was in 2014.In addition, autumn 2018 was the driest and the wettest in 2016.However, the above patterns varied highly with respect to soil depth.
The daily temporal variations in depth-wise soil moisture are presented in Figure 9 along with rainfall events, suggesting that SM in the upper soil zone (0-20 cm) presents high fluctuations and quick responses to PR [12,43], while in the −70 cm depth the range of the soil moisture is smaller and characterized by temporal stability with depth [79,80].In general, the process of increasing moisture in the soil profile follows the distribution of PR.The prevailing low temperatures during the winter months and the occurrence of many rainfall events increase SM values.Nevertheless, the inclination of vertical isothermic curves (Figure 9) in spring varies, which leads us to the conclusion that the response in ST is very rapid down to 70 cm, while in SM, the response is slower mainly under 20 cm and could last until the beginning of the summer [45,74].The SM and ST patterns are inversed during summer compared to winter, and the respective curves display similar but opposite inclinations.In autumn, following the dry summer period when SM minimizes, rainfall events impact highly on SM and ST, which present rapid responses.Their isolines become vertical, with similar SM and ST for the total soil profile (regardless of the soil depth), indicating that the soil profile is uniformly wetting and cooling.Such distribution is generally observed when the soil is rather dry and PR ef amounts exceed 51.9 mm for 3-6 consecutive rainy days (e.g., 252-257 days of the year in 2016 and 321-327 days of the year in 2018).Additionally, from the start of the growing season (April), a gradual decrease in SM values has been observed.This downward trend was interrupted only with the occurrence of precipitation [45,55].It is worth noting that the lack of PR during the growing season (14 PRef rainfall events, 71.9 mm, 2.3 mm/h mean intensity) in 2017 formed a special waterstressed soil profile for vegetation s needs from mid-April until late October.SM values at all depths for this period did not exceed 0.200 cm 3 /cm 3 for at least 5.5 months (late May until late of November) during each monitoring year.Only in 2014, a very wet growing Additionally, from the start of the growing season (April), a gradual decrease in SM values has been observed.This downward trend was interrupted only with the occurrence of precipitation [45,55].It is worth noting that the lack of PR during the growing season (14 PR ef rainfall events, 71.9 mm, 2.3 mm/h mean intensity) in 2017 formed a special water-stressed soil profile for vegetation's needs from mid-April until late October.SM values at all depths for this period did not exceed 0.200 cm 3 /cm 3 for at least 5.5 months (late May until late of November) during each monitoring year.Only in 2014, a very wet growing season (34 PR ef rainfall events, 384.7 mm, 2.0 mm/h mean intensity) limited the water-stressed period to only 3 months (mid-June to early September) and with lower intensity (higher SM values).This situation indicates that the occurrence of rainfall events during the growing season is very important for the eco-hydrology of the study area since it increases water availability, and therefore, the capability of vegetation to overcome the dry summer season or an even larger drought period.
In terms of the water availability for the vegetation, daily SM values revealed that the whole soil profile is in general inside the range between field capacity (θ FC ) and wilting point (θ WP ) for all years.The results showed that during the wet semester, from October to March of the next year, daily SM values at the depth of −70 cm are above θ FC which indicates the occurrence of soil water percolation.The percentages of these days are ranging from 71.4% (October 2013-March 2014) to 65.4% (October 2015-March 2016) but only 10.4% for the dry semester October 2017-March 2018).To the contrary, SM values at the depth of 5 cm from June to late September are below θ WP , which indicates the important role of air temperature increase especially in the top layer of the soil profile.
The combination of phenological observations with the SM and ST patterns described above can enhance our knowledge about the vegetation's behaviour.For example, in our study we detected that the magnitudes of SM were minimized (0.106 on average of all depths) in October 2018 compared to all other months of all years, indicating that the natural vegetation was facing strong water stress.Similarly, a low SM value was also recorded in October 2017 (0.115).The water stress and the warmer conditions, especially in 2018, have probably affected the flowering of Arbutus unedo that generally occurs in our site during October and November.Based on our phenological observations in 2017 and 2018, the beginning of flowering was delayed and started in December.

Conclusions
This study focuses on understanding the dynamics of soil moisture (SM) and soil temperature (ST) on a temporal basis in an experimental evergreen broadleaved mixed maquis plot in Greece by employing six years (2013-2018) of data for the continuous monitoring of the soil profile (0-70 cm) along with aerial meteorological attributes.
Air temperature influences more intensely the ST at the superficial soil depth (0-40 cm).Depth-wise ST changes are more pronounced in winter, when the average ST increases with depth, presenting a maximum difference of +2.5 • C at −70 cm compared to the superficial soil layer (0-5 cm).Of similar magnitude, though negative (−2.4 • C), is the respective ST difference in summer, indicating cooler soil conditions at greater depths.In contrast, the soil profile appears to be uniformly heated at all soil depths during spring, presenting negligible ST differences (indicatively −0.3 • C between −70 cm to 0 cm), whereas in autumn the soil profile pattern for ST is similar to winter but with smaller differences (maximum +1.2 • C between −70 cm to 0 cm).
The increase in SM is closely related to PR ef .During the wet half of the year (October to March), SM at the upper soil layers (−5 cm, −20 cm) shows greater fluctuations compared to the rest of the year.In contrast, SM decreases continuously at all depths during the dry semester (April to September); this is due to the increased water demand of the vegetation and the correspondingly developed dimorphic root system.In spring, which is associated with the beginning of the growing season, the amount of precipitation could lead to limiting the length of the water stress period by up to 2.5 months, considering that summer precipitation in the region is generally scarce or not adequate to significantly increase the amounts of stored soil water in the root zone.This is very important for vegetation to be viable and sustain high growth rates.Before the beginning of the rainy season (in autumn), the soil profile presents minimum SM values (0.065 in August) quite lower than the wilting point (0.111).Thereafter, the impact of precipitation on the formation of the SM profile is critical.It was found that a period of three to six consecutive rainy days, with cumulative PR ef greater than 51.9 mm, is sufficient to uniformly moisten the entire profile.
The above patterns of SM and ST are even more variable monthly, imposing the need for continuous monitoring of the soil profile.They also underline the importance of each of the time steps adopted for the analysis to produce reliable assessments of the effects of SM and ST on natural vegetation growth development and forest soil micrometeorological properties.The mean soil moisture content, computed for 6 years of the study period in the forest stand, was 0.201 cm 3 /cm 3 .This value and the value of 337.0 mm of water that the soil holds under saturated conditions indicate a sufficient potential for soil water storage to sustain local evergreen vegetation.In addition, the current study suggests that the SM and ST temporal variations may be useful tools for identifying the changes in phenology in rainfed ecosystems, especially in the Mediterranean region, but further research is required.In future work we will investigate the role of vegetation on the formation of SM and ST profiles in terms of their phenology and water requirements.

Figure 1 .
Figure 1.Experimental plot and installed monitoring equipment.The red dot of the map in the right corner indicates the study site.

Figure 1 .
Figure 1.Experimental plot and installed monitoring equipment.The red dot of the map in the right corner indicates the study site.
37 • C in winter to 23.9 ± 1.38 • C in summer.The precipitation has an average annual value of 1049 ± 225.5 mm and presents what is expected for the Mediterranean climate seasonal variation: unevenly distributed between seasons with higher amounts falling in winter (40.2%) and autumn (33.3%).The pluviothermic diagram in Figure 2 confirms a water deficit period of about of 3 months starting from June until the beginning of September.

Figure 2 .
Figure 2. The pluviothermic diagram of the 1973-2020 period indicating air temperature Tav (left axis) and bulk precipitation PR (right axis) average monthly values.

Figure 3 .
Figure 3. Soil moisture relation curves against pressure head, h (red line), and hydraulic conductivity, K (green line), as derived from laboratory results and RETC software.Shaded area depicts the soil layer that includes available water for the plants.The continuous red line indicates the fitted pF curve produced by using measured retention data for pressures up to 5 bars, whereas the dashed red line is constructed after extrapolation and refers to higher pressures (θWP was estimated for 15 bars pressure).
7 °C) is quite similar with its long-term respective value for the climatic period 1973-2020 (15.4 °C), ranging on an inter-annual basis from 15.3 °C in 2014 to 16.4 °C in 2018.The air temperature variability is even more sound monthly, and specific months were much warmer or cooler from what was expected based on the climatic profile of the region.The deviations from the local long-term climatic mean for every month are depicted in Figure 4 along with the standard deviations (SD) of the monthly Tav of the 1973-2020 climatic period.

Figure 4 .
Figure 4. Deviations of mean monthly Tav (°C) from the monthly means of the climatic period 1973-2020 (black line in value 0 on the y axis).Shaded areas depict the ± standard deviation (SD) of the same period (confidence interval ±95%).

Figure 3 .
Figure 3. Soil moisture relation curves against pressure head, h (red line), and hydraulic conductivity, K (green line), as derived from laboratory results and RETC software.Shaded area depicts the soil layer that includes available water for the plants.The continuous red line indicates the fitted pF curve produced by using measured retention data for pressures up to 5 bars, whereas the dashed red line is constructed after extrapolation and refers to higher pressures (θ WP was estimated for 15 bars pressure).

21 Figure 3 .
Figure 3. Soil moisture relation curves against pressure head, h (red line), and hydraulic conductivity, K (green line), as derived from laboratory results and RETC software.Shaded area depicts the soil layer that includes available water for the plants.The continuous red line indicates the fitted pF curve produced by using measured retention data for pressures up to 5 bars, whereas the dashed red line is constructed after extrapolation and refers to higher pressures (θWP was estimated for 15 bars pressure).
7 °C) is quite similar with its long-term respective value for the climatic period 1973-2020 (15.4 °C), ranging on an inter-annual basis from 15.3 °C in 2014 to 16.4 °C in 2018.The air temperature variability is even more sound monthly, and specific months were much warmer or cooler from what was expected based on the climatic profile of the region.The deviations from the local long-term climatic mean for every month are depicted in Figure 4 along with the standard deviations (SD) of the monthly Tav of the 1973-2020 climatic period.

Figure 4 .
Figure 4. Deviations of mean monthly Tav (°C) from the monthly means of the climatic period 1973-2020 (black line in value 0 on the y axis).Shaded areas depict the ± standard deviation (SD) of the same period (confidence interval ±95%).

Figure 4 .
Figure 4. Deviations of mean monthly T av ( • C) from the monthly means of the climatic period 1973-2020 (black line in value 0 on the y axis).Shaded areas depict the ± standard deviation (SD) of the same period (confidence interval ±95%).

Figure 5 .
Figure 5. Deviations of total monthly precipitation (mm) from the monthly means of the climatic period 1973-2020.(Shaded areas depict the ± standard deviation (SD) of the same period (confidence interval ±95%)).

Figure 6 .
Figure 6.Distribution frequencies of the total 531 rainfall events in bin classes of (a) mean intensity and (b) duration.The lines are the cumulative frequencies of the classes.

Figure 5 .
Figure 5. Deviations of total monthly precipitation (mm) from the monthly means of the climatic period 1973-2020.(Shaded areas depict the ± standard deviation (SD) of the same period (confidence interval ±95%)).

Figure 6 .
Figure 6.Distribution frequencies of the total 531 rainfall events in bin classes of (a) mean intensity and (b) duration.The lines are the cumulative frequencies of the classes.

Figure 6 .
Figure 6.Distribution frequencies of the total 531 rainfall events in bin classes of (a) mean intensity and (b) duration.The lines are the cumulative frequencies of the classes.

Figure 7 .
Figure 7. Average (a) daily values of soil and air temperatures Tav and (b) monthly diurnal temperature range, at the four depths for an annual period.

Figure 7 .
Figure 7. Average (a) daily values of soil and air temperatures T av and (b) monthly diurnal temperature range, at the four depths for an annual period.

Hydrology 2023 , 21 Figure 8 .
Figure 8. Soil profile temperatures (ST) in °C for the years 2013-2018 based on daily values.

Figure 8 .
Figure 8. Soil profile temperatures (ST) in • C for the years 2013-2018 based on daily values.

Hydrology 2023 , 21 Figure 9 .
Figure 9. Changes in soil moisture SM profile in cm 3 /cm 3 based on daily values (left axis) and daily effective precipitation PRef (right axis), for the years 2013-2018.

Figure 9 .
Figure 9. Changes in soil moisture SM profile in cm 3 /cm 3 based on daily values (left axis) and daily effective precipitation PR ef (right axis), for the years 2013-2018.

Table 1 .
Average values of texture fractions (%), pH and concentrations of C and N (g/kg).Values of coefficient of variations (%) are shown inside parenthesis.

Table 2 .
Monthly T av , average monthly PR and PR ef (mm) for the study period 2013-2018.

Table 3 .
Soil temperature (ST) variations with depth for the period 2013-2018.

Table 4 .
Soil moisture (SM) variations with depth for the period 2013-2018.