Issue 3/2023

First data on rainfall interception in an Atlas cedar forest (Cedrus atlantica Manetti) in the Aurès (eastern Algeria)

Article by Toufik Benhizia, Salim Lebbal and Abdelghafour Abaidia | 28.09.2023 - 10:03

Schlüsselbegriffe: semi-aride Wälder, Stammabfluss, Interzeption, Kronen­durchlass, Klimawandel

Available at https://doi.org/10.53203/fs.2303.2

Download as PDF

See below the issue 3/2023 as E-Paper or have a look at our E-Paper archive dating back to 1955.

Abstract

With climate change, data on throughfall, stemflow and rainfall interception in forests are much needed by hydrologists throughout the world, including the Mediterranean region. Such data are rare for Algeria, especially for cedar forests. Thus, the present study is interested in analysis of throughfall, stemflow and rainfall interception in a cedar forest dominated by Cedrus atlantica Manetti in the province of Batna (Eastern Algeria) and establishing models for forecasting the different rainfall fractions in this forest type. Specific containers (plastic bottles and collars) were installed at stems of 30 trees, under tree crowns and in an open area not covered by tree crowns between January 2018 and December 2019 to collect the different studied fractions of water (rainfall, throughfall and stemflow). Our results reveal that throughfall was with 65.41% of the total rainfall the greatest fraction in the studied cedar forest. Stemflow represented only a very small fraction of total rainfall with 3.86%. The mean rainfall interception was 30.72 % of the total rainfall. The relationships between rainfall, throughfall and stemflow were assessed by linear regressions. We found strong correlations between throughfall, stemflow and rainfall with coefficient of determination R2 ranging between 0.92 and 0.99. The equations developed by this study revealed that the studied cedar forest must receive at least 1.70 and 4.46 mm of rainfall, to trigger throughfall and stemflow, respectively. 

Zusammenfassung

Angesichts des Klimawandels sind Daten über Kronendurchlass, Stammabfluss und Interzeptionin Wälder für Hydrologen auf der ganzen Welt, einschließlich der Mittelmeerregion, sehr gefragt. Dennoch sind solche Daten in Algerien, insbesondere für Zedernwälder, nicht verfügbar. Diese Studie beschäftigt sich mit der Quantifizierung des Kronendurchlass, des Stammabflusses und der Interzeption von Niederschlag in einem Zedernwald dominiert durch Cedrus atlantica Manetti in der Provinz Batna (Ostalgerien) und mit der Entwicklung von Modellen zur Vorhersage der verschiedenen Niederschlagsfraktionen in diesem Waldtyp. Zwischen Januar 2018 und Dezember 2019 wurden spezielle Behälter (Plastikflaschen und -manschetten) installiert an 30 Bäumen, unter Baumkronen und auf Freiflächen ohne Kronenbedeckung um die verschiedenen untersuchten Niederschlagsfraktionenzu sammeln.

Die Ergebnisse zeigten, dass der Kronendurchlass 65.41 % des Niederschlags im untersuchten Zedernwald betrug. Der Stammabfluss stellte nur einen sehr kleinen Teil des Gesamtniederschlags dar mit 3.86 %. Die mittlere Interzeption betrug 30.72 % des Gesamtniederschlags. Die Beziehungen zwischen Niederschlag, Kronendurchlass und Stammabfluss wurden mit linearen Regressionen analysiert. Wir beobachteten starke Korrelationen zwischen Kronendurchlass, Stammabfluss und Niederschlag mit Bestimmtheitsmass R2 zwischen 0.92 und 0.99. Die in dieser Studie entwickelten Gleichungen ergaben, dass der untersuchte Zedernwald mindestens 1.70 bzw. 4.46 mm Niederschlag erhalten muss, damit Kronendurchlass bzw. Stammabfluss messbar ist.

1 Introduction

Rainfall is one of the most important sources of water in forest ecosystems and it is responsible for forest productivity, especially in seasonally dry forests (Tague et al. 2019; Missaoui et al. 2020; Yang et al. 2022; 2023; Zhao et al. 2023). Rainfall is the main input to the soil water balance and plays an important role in the water cycle and the plant-soil-atmosphere system (Biniak-Pieróg 2017; Klamkowski et al. 2011; Levia et al. 2017; Schwingshackl et al. 2017; Han et al. 2018; Magliano et al. 2019; Sun et al. 2022; Treder et al. 2022). The distribution of rainfall by trees is an important hydrological process, particularly in the context of climate change (Llorens & Domingo, 2007; Sadeghi et al. 2018; Yang et al. 2022; Zhao et al. 2023). The two main processes, evapotranspiration and interception, are strongly influenced by vegetation and its composition (Vicente et al. 2018). With climate change, any alterations in the precipitation or temperature regimes will directly affect the dynamics of the interception of rains by forests. Previous studies pointed out that with climatic variations, rainfall may be concentrated in a shorter time period with more intense rainfall events and in consequence higher drought risks (Benhamiche et al. 2014; Zhao et al. 2023). It is increasingly recognized that a better understanding of water use patterns and their response to climate change depends on our ability to better quantify the dynamics of rainfall fractions (Yang et al. 2023).

Therefore, an understanding of the redistribution of rainfall by the action of the forest cover is very useful to understand the water balance and the nutrients of the ecosystem, but also to develop the silvicultural practices, that mitigate the effects of climate change (Jeong et al. 2019; Tague et al. 2019; Shi et al. 2022; Zhao et al. 2023). During a rainfall event, the forest cover causes, on the one hand, a reduction in the rainfall amount reaching the forest floor, and on the other hand, a redistribution of the latter into three fractions, throughfall, stemflow and interception (Dohnal et al. 2014; Sadeghi et al. 2020; Jančo et al. 2021; Zhao et al. 2023). Throughfall penetrates the tree crowns, while stemflow is water running along the surface area of tree trunks (Lévia et al. 2017). The amounts of rainfall, that is captured by forest canopy and then returned to the atmosphere by evaporation, are often referred to as interception (Šraj et al. 2008; Dohnal et al. 2014; Jančo et al. 2021).

The Mediterranean area is more strongly affected by climate change than other regions globally (GIEC 2014; Lionello et al. 2018; IPCC 2021). The expected changes include hotter and drier summers as well as a general decrease in rainfall as large as minus 40% in certain regions (Barros et al. 2014; Peñuelas et al. 2017; Arar et al. 2020; Tuel et al. 2020).

The Atlas cedar, Cedrus atlantica Manetti, is endemic to North African mountains and constitutes the southwesternmost species of the genus Cedrus (Bouahmed et al. 2019). In Algeria, it constitutes a much values forest asset of remarkable ecological, economical and cultural importance (Médail & Quézel 2003). It is also classified as an endangered species and carries the status of a rare plant (Yahi et al. 2007; Touati et al. 2021). The Algerian cedar forests occupy an area of about 2089 km², where the stands located in the Aurès regions are the most xeric in Algeria (Arar et al. 2020) and they have suffered from tree mortality linked to drought and climate change for several years (Allen 2009; Yahi & Djellouli 2010; Kherchouche et al. 2012; Bezzih et al. 2021). Bouahmed et al. (2019) indicated that the decline of Algerian cedar forests is driven by a shift of climate towards drier conditions. This is problematic as Atlas cedar is characterized by strong height growth, that is highly sensitive to drought (Ladjal et al. 2007).

Growth-climate relationships for Atlas cedar showed a significant positive, cumulative and stable temporal effect of precipitation variability (Slimani et al. 2014). A negative growth shift was triggered by a climate shift towards drier conditions in the 1980s (Navarro-Cerrillo 2019). A very high decrease in suitable areas of Atlas cedar was forecasted based on the predictions of Arar et al. (2019) for 2070; which indicates that the species has little chance to survive under the conditions of future climate change scenarios. More arid conditions would raise Atlas cedar mortality rates and trigger compositional shifts toward forests dominated by species better able to tolerate drought, such as Quercus ilex (Navarro-Cerrillo 2019).

Studies of rainfall interception have long received considerable attention worldwide (Gash et al. 1980; Llorens & Domingo 2007; Muzylo et al. 2009; Zhao et al. 2023) and remains an important topic for future environmental studies (Zhang et al. 2023). Many researchers have already recommended continuing studies on the redistribution of rainfall by forests, particularly for semi-arid and arid regions, to better understand the interception process in the context of climate change and to balance the water demands of humans and the ecosystem (Friesen & Van Stan 2019; Jeong et al. 2019; Ma et al. 2019; Levia et al. 2019; Yang et al. 2022; Zhao et al. 2023). It is noticed that the majority of available studies on rainfall interception have focused on cumulative annual amounts and not on individual values of interception by rainfall events (Muzylo et al. 2009). Although semi-arid environments represent a third of the world’s land area (Li et al. 2020; Mirzabaev et al. 2022), little is known from a point of view of rainfall interception compared to other climates (Llorens & Domingo et al. 2007; Muzylo et al. 2009; Magliano et al. 2019).

Despite the obvious significance of this topic, there is currently no data on the interception of rainfall in Algeria by the action of the Atlas cedar. Our study aims fill this gap and to provide models to predict the interception of rainfall by traa canopies for the first time for the Atlas cedar. A novel aspect of our study is exploring rainfall interception for single events for an endemic species in a semi-arid climate. 

2 Materials and methods

Ced_Fig1.jpg

Figure 1: The location of the study site and distribution of Cedrus atlantica (Caudullo et al. 2020).

Abbildung 1: Der Standort des Untersuchungsortes und die Verbreitung von Cedrus atlantica (Giovanni et al. 2020).

2.1 Measurements of rainfall, throughfall, stemflow and interception

The study was conducted on the interception of rainfall by the action of Atlas cedar, in an area located in the northeastern part of Algeria in the province of Batna (Fig. 1). More precisely, the Atlas cedar forest was located in the Larbaa mountains (35°39'50''N, 6°12'31''E) at an altitude of 1680 m with average tree density of 600 trees per hectare and an estimated leaf area index (LAI) between 1.7 and 2.5. We estimated LAI based literature (Breda 1999; Breda et al. 2002) and the observation that the canopy of our studied forests is relatively sparse and suffered from drought effects. The tree height ranged between 14 and 18 m and their age was approximated 90 years. Tree diameter at breast height ranged between 26 and 32 cm and the basal area of the studied cedar forest was 32 m2/ha.  

The data provided by the meteorological station of Ain Skhouna (35°75'79''N, 6°30'57''E), about 43 km in the northeast of the study site, revealed that the average rainfall in the Batna region over 25 years (1989-2013) was 332.83 mm with an average air temperature of 15.58 °C. 

Ced_Fig2.jpg

Figure 2: Devices for collecting stemflow (A) and throughfall (B). 

Abbildung 2: Geräte zum Sammeln von Stammabfluss (A) und Kronendurchlass (B).

The collectors used to measure rainfall, throughfall and stemflow were deduced from previous studies (Gash et al. 1980; Molina & del Campo 2012; Llorens & Domingo 2007; Attarod et al. 2015). Because of its adaptability to forecast rainfall interception, particularly for coniferous forests in semi-arid climates, the Gash model was chosen (Llorens & Domingo 2007; Muzylo et al. 2009). Water volumes (rainfall, throughfall and stemflow) were measured continuously (after each rainfall event) for two years from January 2018 until December 2019. 

The different fractions of water were collected, just after a rainfall event finished, using plastic bottles with the upper part cut and turned into a funnel (Fig. 2). The receiving area for each bottle was 0.0227 m². Throughfall was quantified by four bottles installed under the central and peripheral parts of each tree. Previous works have stated, that rain gauges should be placed not only in the central area of the crown, near the trunk of the trees, but also in the periphery to ensure more reliable results (Bartík et al. 2016; Dohnal et al. 2014; Minďáš et al. 2018). The bottles used to measure the throughfall were installed under crowns of 30 Atlas cedar trees, for a total of 120 bottles. Our overall surface for receiving throughfall water (2.72 m²) is considered sufficient to correctly estimate this fraction of rainfall (Rodrigo & Avila 2001; Llorens & Domingo 2007). On a bare plot further 20 bottles were placed identical to those used to measure the throughfall and in approximately eight meters distance from the Atlas cedar trees, to measure rainfall undisturbed by tree canopies.

Stemflow was measured using plastic collars fixed around the trunks of 16 trees and sealed with mastic. These collars were connected with a plastic container using a pipe. The collected rainfall volumes are converted into mm using the formula (Vialard-Goudou & Richard 1956):

with P rainfall (mm), S reception surface (cm²) and V volume of water collected (cm3).

Stemflow was calculated by dividing the collected volumes by the crown surface (Livesley et al. 2014). The interception in mm were estimated according to (Llorens & Domingo 2007; Moreno-Pérez et al. 2018): 

with I interception (mm), P rainfall (mm), Tf throughfall (mm) and Sf is stemflow (mm).

 

2.2 Data analysis

The results of rainfall, throughfall, stemflow and interception were subject to analysis of variance (ANOVA) at the 5% error threshold, using SPSS software version 10.0.5 (SPSS Inc.). The relationships between the different fractions studied (P, Tf, Sf, I) were estimated by linear regressions using Microsoft Excel 2013. We calculated the amount of rainfall needed to trigger Tf and Sf by reformulating the linear functions for Sf and Tf equal 0.1 mm. 

3 Results

3.1 Rainfall

The results of two years of measurements showed that the cumulative daily rainfall was 773.6 mm, created by 118 days with rainfall. The year 2018 received more rainfall with 436.6 mm, compared to year 2019 with 337 mm.

Daily rainfall fluctuated between 0.9 and 56 mm with an average of 4.28 mm per day of rainfall. The distribution by rainfall class ranging from 0 to 60 mm and their frequencies are available in Table 1.

Ced_Tab1.jpg

Table 1: Classes, frequencies and accumulations of daily rainfall.

Tabelle 1: Klassen, Häufigkeiten und Anhäufungen der täglichen Niederschläge.

Table 1 demonstrated the dominance of rainfall less than or equal to 5 mm per day, throughout the study period. In 2018, most rainfall was created by rainfall events with 30-60 mm per day, while in 2019 most rainfall originated from events with 10-20 mm per day. It should be noted that only about 16% of days in the measurement period received rainfall greater 0 mm. From these 118 days with rainfall, 73 days has rainfall between 0 and 5 mm and the number of days with rainfall greater 30 mm was very low with 4 days out of 118. Surprisingly, the cumulative rainfall of these four days was 162 mm (19.66 % of the total rainfall).

 

3.2 Throughfall

Our results revealed that among 118 rainy days, only 92 events could trigger throughfall, penetrating the crowns of the Atlas cedar forest. Accumulation throughfall over two years was 509.32 mm, which represents an average of 65.41% of total annual rainfall. We note that the highest throughfall rate (68.69%) was recorded in year 2018 (Fig. 3), which was also received more rainfall compared to the year 2019 (436.6 and 337 mm, respectively).

Ced_Fig3.jpg

Figure 3: Annual percentages for throughfall, stemflow and interception during the study period.

Abbildung 3: Jährliche Prozentsätze für Kronendurchlass, Stammabfluss und Interzeption während des Untersuchungszeitraums.

The average daily throughfall for 2018 (0.82 mm/day) was much greater compared to the throughfall for the year 2019 (0.54 mm/day). It should also be noted that during the year 2019, the throughfall for the two months of June and July was nil. Throughout the study period, the largest amounts of throughfall were recorded between March and May. Figure 4 presents the monthly variations for the two years 2018 and 2019.

Ced_Fig4.jpg

Figure 4: Monthly variation of throughfall during the study period.

Abbildung 4: Monatliche Variation des Kronendurchlasswährend des Untersuchungszeitraums.

3.3 Stemflow

The number of rainfall days that could trigger stemflows, through the cedar tree trunks, was only 37 days, that is a rate of 31.35% of the total number of rainfall days during our study period. The results also showed that the average stemflow rate was around 3.86%. The stemflow for 2019 was slightly higher than for 2018 (4.01 and 3.71%, respectively) (Fig. 3).

The average volume per rainfall event of water brought to the soils of the cedar forest through stemflow fluctuates between 0.03 and 0.05 mm for 2019 and 2018 respectively.

On an interannual scale, the two months of June and July were characterized by zero stemflow, while, the largest stemflows were recorded between March and May (2.2 - 3.4 mm) (Fig.5).

Ced_Fig5.jpg

Figure 5: Monthly variation of stemflow volumes throughout the study period.

Abbildung 5: Monatliche Variation der Stammabfluss während des Untersuchungszeitraums.

3.4 Interception

The average interception rate found for this study is 30.72%. ANOVA analysis confirmed a significant difference between years (p= 0.010) (Table 2). The interception was higher in the least rainy year. It was 34.15% in 2019 (where the rainfall was 337 mm) and 27.30% in 2018 (where the rainfall was 436.6 mm) (Fig. 3). The mean interception varies between 0.32 and 0.33 mm per rainfall day throughout the study period.

As for the monthly variation, the ANOVA analysis revealed a significant difference (p= 0.002) (Table 2). The interception was much greater in March 2019 (20.2 mm) and May 2018 (24 mm) (Fig.6).

Ced_Tab2.jpg

Table 2: ANOVA analysis of the studied parameters at an intra-annual and inter-annual scale. Asterisk indicate significant difference at p < 0.05.

Tabelle 2: ANOVA-Analyse der untersuchten Parameter auf einer inter- und intra-annualen Skala. Stern zeigt signifikante Unterscheide mit p < 0.05.

Ced_Fig6.jpg

Figure 6: Monthly variation of rainfall interception throughout the study period.

Abbildung 6: Monatliche Variation der Niederschlagsinterzeption während des gesamten Untersuchungszeitraums.

3.5 Relationship between rainfall and throughfall

The results relating to the relationship between rainfall and throughfall during both years, under the Atlas cedar, are presented in Figure 7. The latter shows a positive linear relationship between throughfall (Tf) and rainfall (P) (R2 = 0.99). There is a proportional relation between these two parameters.

The equations presented in Figure 7 show that Atlas cedar trees must receive at least 1.70 mm of rain to allow the rain to pass to the forest soils and subsequently feed the groundwater.

Ced_Fig7.jpg

Figure 7: Relationship between rainfall and throughfall in the studied cedar forest.

Abbildung 7: Zusammenhang zwischen Niederschlag und Kronendurchlass im untersuchten Zedernwald.

3.6 Relationship between rainfall and stemflow

Figure 8 shows that stemflow is positively correlated with rainfall amounts (R2 =0.92). The relationships between the rainfall and stemflow using a linear trend function are presented in Figure 8. By rearranging equation 8, we can conclude that the stemflow is triggered once 4.46 mm of rain is exceeded.

Ced_Fig8.jpg

Figure 8: Relationship between rainfall and stemflow in the studied cedar forest.

Abbildung 8: Zusammenhang zwischen Niederschlag und Stammabfluss im untersuchten Zedernwald.

3.7 Relationship between rainfall and interception

Figure 9 shows that the volumes of interception are negatively correlated with the volumes of rainfall with coefficient of determination R2 varies between 0.18 and 0.63.

Ced_Fig9.jpg

Figure 9: Relationship between rainfall and interception.

Abbildung 9: Zusammenhang zwischen Niederschlag und Interzeption.

4.1 Rainfall

The average annual rainfall of 386.8 mm is higher than the average rainfall (332.83 mm) of the Batna region for a period of 25 years (1989-2013). Nevertheless, hydrologists are interested not only in annual rainfall, but also in rainfall events, as the latter are more important for throughfall, stemflow and interception than annual rainfall sums. In the present investigation, rainfall events were dominated by rainfall below 5 mm with an average of 4.28 mm per event.

Similarly, Benhizia et al. (2020) found, forthe Aleppo pine in a semi-arid region, that the daily rainfall was low, as daily rainfall higher than 30 mm was rare, but when it occurs, it can constitute a very significant proportion of the annual rainfall (33.52%). Moreover, other researchers stated, that in semi-arid and Mediterranean climates, daily rainfall of less than 10 mm constitutes almost all of the rainfall (Pérez-Suárez et al. 2013; Llorens et al. 2018). Our results suggest that the rainy days that can pass through the canopy of the Atlas cedar and feed the soil and the water table are much reduced. Indeed, Andreu et al. (2001) indicated that only daily rainfall greater than 15 mm was capable of producing groundwater recharge.

Vicente et al. (2018) have already asserted that rainfall events will become less frequent, which directly affects the deep percolation to groundwater. Similarly, Barros et al. (2014) and Vicente et al. (2018) reported that drought conditions are also expected to intensify, which could trigger plant mortality and species replacement events and consequently the forest water balances will be affected.

 

4.2 Throughfall

Studies on throughfall under the action of the Atlas cedar in a forest environment are rare or even nonexistent. Sensoy and Tanyel (2022) published the only study on Cedrus libani in an urban environment. Unfortunately, these researchers worked on rainfall greater than 10 mm and they found a throughfall rate of 51%. This rate is considered low compared to the present study (65.41%).

The scarcity of work on the redistribution of rainfall under the action of the Atlas cedar in forest environments led us to compare our results with studies that have worked on other conifers. The throughfall rate of the present study remains high compared to the results of Benhizia et al. (2020) which found for an Aleppo pine forest, a throughfall that varies between 51.36 and 56.95%. However, rainfall volumes and their characteristics as well as forest cover are not the same. The rainfall volumes of the present study are higher than the results of the latter researchers (386.8 vs. 312 mm). This difference may be, on the one hand, responsible for the increase in the rate of throughfall in our cedar forest and, on the other hand, dead branches and cones are abundant in the pine forest studied by the aforementioned researchers; and they are not too abundant in our cedar forest. Thus, Jeong et al. (2019) have claimed that for conifers, the drainage is low for stands characterized by dead branches because of the additional gain of rainwater from dead branches.

Besides, several previous studies stated that the throughfall in coniferous forests varied between 48.5 and 85% (Rapp & Romane 1968; Cape et al. 1991; Barbier et al. 2009; Zhou et al. 2013; Fan et al. 2014; Sun et al. 2014; Aydın et al. 2018; Jeong et al. 2019; Magliano et al. 2019; Dong et al. 2020; Yan et al. 2022).

Because the intra-annual variation was considerable, it must be pointed out here that the low throughfall, between June and August in the studied cedar forest, can negatively affect the germination of the Atlas cedar. Furthermore, Bentouati and Oujehih (1999) mentioned that the growth of Atlas cedar forests in the Aurès is conditioned by the substrate and the rainfall. They indicated that the total production of cedar forests in the Aurès is between 31 and 1007 m3/ha depending on rainfall, density and age of the forest trees. In our case, it seems that cedar forests that rest on deep soils will benefit more from throughfall compared to stands that grow on soils not deep.

 

4.3 Stemflow

Studies on stemflow are not very abundant, particularly for conifers. Among these rare studies, some revealed a stemflow rate that fluctuates between 0.6 and 1.74% (Llorens et al. 1997; Llorens & Domingo 2007; Xiao et al. 2007; Ma et al. 2019; Dong et al. 2020). The stemflow rate found in the present study (3.86%) is relatively high compared to the results of the aforementioned researchers. The difference can be attributed to daily rainfall patterns and stand characteristics. Therefore, Fan et al. (2014) claimed that stemflow rate is closely related to daily rainfall amounts. Besides, the rate found by the present study can be very beneficial for cedar forests. Stemflow represents a minor proportion compared to throughfall, but it can concentrate 17.9 to 56.6 times the amount of precipitation in the soil profile adjacent to the root zone (Zhao et al. 2023). 

Furthermore, other studies indicated that the stemflow represents only a very small proportion of rainfall and its rate ranges between 0.8 and 9.1% depending on precipitation, density, diameter and height of trees (Barbier et al. 2009; Saito et al. 2013; Pérez-Suárez et al. 2014; Aydın et al. 2018; Dong et al. 2020; Yang et al. 2022).

 

4.4 Interception

The interception rate of the present study (30.72%) is relatively low compared to the results of Benhizia et al. (2020) who found an average interception rate of 42.5% in an Aleppo pine forest, although Keim et al. (2006) found for Pinus spp an interception rate of 28.5%. Nevertheless, the average rainfall found in the present study (386.8 mm) is higher than that found by Benhizia et al. (2020) (311.95 mm). In addition, the cumulative rainfall between 30 and 60 mm in the present study, is greater than that in the study of Benhizia et al. (2020). Consequently, the low interception found in the present study can be attributed to rainfall concentrated on four days among 118 rainy events.

The other factor that probably reduced the interception rate in the present study, is the reduced number of dead branches and cones on the cedar trees compared to the pine trees where they were abundant. Indeed, some researchers have already reported that cones and dead branches can significantly modify the redistribution of rainfall by the action of conifers due to the additional water gain of rain by the latter (Shinohara et al. 2010; Jeong et al. 2019).

Besides, the time between two successive rainfalls can also influence the interception rate. Rodrigo and Avila (2001) have set the canopy drying time at 4 hours, whereas Cape et al. (1991) estimate it at 12 hours. Moreover, the rainfall regime of the Mediterranean is probably the first to be responsible for the interception rates (Bellot et al. 1999; Llorens et al. 2018).

Nevertheless, the interception rate found in our study remains, in general, close to the results recorded in coniferous forests. Thus, previous studies have reported interception rates between 8 and 81% under coniferous species (Llorens et al. 1997; Xiao et al. 2007; Barbier et al. 2009; Zhou et al. 2013; Fan et al. 2014; Aydın et al. 2018; Dong et al. 2020; Jančo et al. 2021; Yang et al. 2022). 

Regarding the significant interannual variation, Moreno-Pérez et al. (2018) stated that the interception rate is higher during dry years.

 

4.5 Relationships with rainfall

The equations found in the present study, that elicit the relationship between rainfall, throughfall and stemflow are consistent with previous studies. Indeed, several researchers (Magliano et al. 2019; Dong et al. 2020; Zhao et al. 2023) have already reported these positive correlations between rainfall, throughfall and stemflow. They indicated that a major rainfall event tends to favor a higher percentage of throughfall and stemflow.

The volume necessary to trigger the throughfall in the Atlas cedar (1.70 mm), is close to that found by Benhizia et al. (2020). The latter found a volume of 1.9 mm, for the Aleppo pine in the Aurès (Eastern Algeria). Similarly, Dong et al. (2020) claimed that the rainfall capable of generating throughfalls in the Pinus tabuliformis stand is between 1.5 and 1.9 mm depending on the age.

On the other hand, the volume necessary for the initiation of stemflow (4.46 mm) is less compared to that found by Benhizia et al. (2020), which found a volume of 6 mm is required to start stemflow in Aleppo pine. Moreover, the volumes capable of triggering throughfall and stemflow in the present investigation are higher compared to those stated by Zhang et al. (2023) who claimed that the average rainfall required for the initiation of throughfall and stemflow, in the forest canopies, begins on average at 1.2 and 3.3 mm, respectively. This difference can be attributed to the bark properties of cedar but and to the hydrophobous dust accumulating on the needles and trunks of trees, especially after long drought periods. 

Furthermore, the difference in rainfall volumes triggering the stemflow in cedar and pine forests is most likely due to the presence of dead branches on Aleppo pine trunks, especially for the first few meters from the level of the soil. Dead branches attached to the stem trunks will likely hinder the stemflow in Aleppo pine and Cedar forests.

Applying the developed equations suggests that rainfall exceeding 5 mm will contribute to the water supply of Atlas cedar forests.

5 Conclusions

Our study presents the first data on the redistribution of rainfall in an Algerian Atlas cedar forest. The results of the two years of measurements recorded a cumulative rainfall of 773.6 mm distributed over 118 events of rain, including only four events with rainfall above 30 mm. Moreover, the rate of rainfall that crossed the crowns of the Atlas cedar trees by throughfall and stemflow was 65.41% and 3.86% respectively, while the interception rate was 30.72%. Thus, the interception was considerable and varied between 0.32 and 0.33 mm per rainfall day throughout the study period.

The aforementioned data are important and they can help foresters for effective management, in particular to balance the water demand of a forest ecosystem and the demand for water by humans. Thinning and removal of dead branches in dense Atlas cedar forests may be beneficial by increasing the amount of rainfall reaching forest soils and increasing the groundwater recharge with positive effects on water supply to societies and ecosystems. This study of rainfall interception has theoretical and practical importance in providing insight into the roles of certain climatic variables and cedar stand traits in the redistribution of rainfall, which can help determine of the impacts of reforestation with cedar on the local water balance.

As a perspective and to bring new knowledge relating to rainfall interception under the action of cedar forests in Aures (eastern Algeria), it is recommended to consider the following rarely considered meteorological variables: rainfall intensity, wind speed and wind direction, vapor pressure deficit and the number of wetting/drying cycles within a given rainfall event.

References

Allen R. M., 2009. Warming caused by cumulative carbon emissions towards the trillionth tone. Nature, 458: 1163–1166.

Andreu J.M., Delgado J., García-Sánchez E., Pulido-Bosch A., Bellot J., Chirino E., Ortíz de Urbina J.M., 2001.Caracterización del funcionamiento y la recargadelacuíferodelVentós-Castellar (Alicante). Revista de la Sociedad Geológica de España, 14 : 247–254.

Arar A., Nouidjem Y., Bounar R., Tabet S., Kouba Y., 2020. Modeling of the current and future potential distribution of Atlas cedar (Cedrus atlantica) forests revealed shifts in the latitudinal, longitudinal and altitudinal range towards more humid conditions. Ecological Questions, 31(3): 49–61.

Arar A., Tabet S., Nouidjem Y., Bounar R., Chenchouni H., 2019. Projected small-scale range reductions of Cedrus atlantica forests due to climate change at the Belezma National Park (Algeria). In Exploring the Nexus of Geoecology, Geography, Geoarcheology and Geotourism: Advances and Applications for Sustainable Development in Environmental Sciences and Agroforestry Research: Proceedings of the 1st Springer Conference of the Arabian Journal of Geosciences (CAJG-1), Tunisia 2018 (pp. 15-19). Cham: Springer International Publishing.

Attarod P., Sadeghi S.M.M., Pypker T.G., Bagheri H., Bagheri M., Bayramzadeh V., 2015. Needleleaved trees impacts on rainfall interception and canopy storage capacity in an arid environment. New Forests, 46: 339-355.

Aydın M., Sen S.G., Celik S., 2018. Throughfall, stemflow, and interception characteristics of coniferous forest ecosystems in the western black sea region of Turkey (Daday example). Environmental monitoring and assessment,190: 1-10.

Barbier S., Balandier P., Gosselin F., 2009. Influence of several tree traits on rainfall partitioning in temperate and boreal forests: a review. Annals of Forest Science, 66(6): 1-11.

Barros V.R., Field C.B., Dokken D.J., Mastrandrea M.D., Mach K.J., Bilir T.E., Chatterjee M., Ebi K.L., Estrada Y.O., Genova R.C., 2014. Climate Change Impacts, Adaptation, and Vulnerability – Part B: Regional Aspects – Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, New York, p. 688.

Bartík M., Jančo M., Střelcová K., Škvareninová J., Škvarenina J., Mikloš M., Vido J., Waldhauserová P.D., 2016. Rainfall interception in a disturbed montane spruce (Picea abies) stand in the West Tatra Mountains. Biologia, 71: 1002–1008.

Bellot J., Sanchez J.R., Chirino E., Hernandez N., Abdelli F., Martinez J.M., 1999. Effect of different vegetation type cover on the soil water balance in semi-arid areas of south eastern Spain. Physics and Chemistry of the Earth, Part B: Hydrology, Oceans and Atmosphere, 24(4): 353–357. 

Bellot J., Avila A., Rodrogo A., 1999. Throughfall and Stemflow. In: Ecological Studies (Rodà F., Retana J., Gracia C.A., Bellot J., eds). Springer-verlag: Heidelberg, pp. 209-222.

Benhamiche N., Madani K., Laignel B., 2014. Impact of climate changes on water resources in Algeria. Vulnerability of agriculture, water and fisheries to climate change: toward sustainable adaptation strategies, pp.193-205. Springer Science+Business Media: Dordrecht.

Benhizia T., Lebbal S., Abdenbi Z. A., Zeraib A., 2020. Action du couvert forestier du pin d’Alep (Pinus halepensis Mill.) sur la redistribution des pluies sous un bioclimat Semi-aride de l’Est Algérien. Annales de l’INRGREF, 21: 17-34.

Bentouati A., Oudjehih B., 1999. Première étude de la croissance et de la productivité du Cèdre de l'Atlas (Cedrus atlantica Manetti) dans le massif de Bélezma (Aurés Algérie). Forêt méditerranéenne, 20(2): 115-119.

Bezzih H., Malki H., Aissi A., 2021. Characterising of plant cover evolution in the aurès region (eastern Algeria): using landsat imagery. Geoadria, 26(2):111-124.

Biniak-Pieróg M., 2017. Monitoring of atmospheric precipitation and soil moisture as basis for the estimation of effective supply of soil profile with water. Monografie, 207.

Bréda N., 1999. L'Indice foliaire des couverts forestiers: mesure, variabilité et rôle fonctionnel. Revue forestière française, 51(2): 135-150.

Bréda N., Soudani K., Bergonzini J.C., 2002. Mesure de l'indice foliaire en forêt. Ecofor: Clamecy. 155p.

Carrera-Hernandez J.J., Gaskin S.J., 2008. Spatio-temporal analysis of potential aquifer recharge: Application to the basin of Mexico. Journal of Hydrology, 353: 228–246. 

Cape J.N., Brown A.H.F., Robertson S.M.C., Howson G., Paterson I.S., 1991. Interspecies comparisons of throughfall and stemflow at three sites in northern Britain. Forest Ecology and Management, 46 (3-4): 165-177. 

Caudullo G., Welk E., San-Miguel-Ayanz J. 2020. Cedrus atlantica chorology. figshare. Dataset. https://doi.org/10.6084/m9.figshare.12287441.v1  

Dohnal M., Černý T., Votrubová J., Tesař M., 2014. Rainfall interception and spatial variability of throughfall in spruce stand. Journal of Hydrology and Hydromechanics, 62, 277–284.

Dong L., Han H., Kang F., Cheng X., Zhao J., Song X., 2020. Rainfall partitioning in Chinese pine (Pinus tabuliformis Carr.) stands at three different ages. Forests, 11(2): 243. 

Fan J., Oestergaard K.T., Guyot A., Jensen D.G., Lockington D.A., 2014. Spatial variability of throughfall and stemflow in an exotic pine plantation of subtropical coastal Australia. Hydrological Processes, 25(5): 793-804. 

Friesen J., Van Stan J.T., 2019. Early European observations of precipitation partitioning by vegetation (2019). A Synthesis and evaluation of 19th century findings. Geosciences, 9(10): 423.

Gash J.H.C., Wright I.R., Lloyd C.R., 1980. Comparative estimates of interception loss from three coniferous forests in Great Britain. Journal of Hydrology, 48(1-2): 89-105. 

GIEC, 2014. Changements climatiques 2014. Rapport de synthèse du GIEC, Suisse, 161 p.

Han X., Gao G., Chang R., Li Z., Ma Y., Wang S., Wang C., Lü Y., Fu B., 2018. Changes in soil organic and inorganic carbon stocks in deep profiles following cropland abandonment along a precipitation gradient across the Loess Plateau of China. Agriculture, Ecosystems & Environment, 258: 1-13.

IPCC, 2021.  Climate Change 2021. The physical science basis working group i contribution to the sixth assessment report of the intergovernmental panel on climate change. 40 P.

Jančo M., Mezei P., Kvas A., Danko M., Sleziak P., Minďáš J., Škvarenina J., 2021. Effect of mature spruce forest on canopy interception in subalpine conditions during three growing seasons. Journal of Hydrology and Hydromechanics, 69(4): 436-446.

Jeong S., Otsuki K., Farahnak M., 2019. Relationship between stand structures and rainfall partitioning in dense unmanaged Japanese cypress plantations. Journal of Agricultural Meteorology, 75(2): 92-102.

Keim R.F., Skaugs A.E., Weiler M., 2006. Storage of water on vegetation under simulated rainfall of varying intensity. Advances Water Resources, 29(7): 974-986. 

Kherchouche D., Kalla M., Gutiérrez E.M., Attalah S., Bouzghaia M., 2012. Impact of droughts on Cedrus atlantica forests dieback in the Aurès (Algeria). Journal of Life Sciences, 6: 1262-1269.

Klamkowski K., Treder W., Tryngiel-Gać A., 2011. Influence of mulch on rainfall interception in apple orchard. Sodininkystėir Daržininkystė, 30(2): 15-21.

Ladjal M., Deloche N., Huc R., Ducrey M., 2007. Effects of soil and air drought on growth, plant water status and leaf gas exchange in three Mediterranean cedar species: Cedrus atlantica, C. brevifolia and C. libani. Trees - Structure and Function, 21(2): 201–213.

Levia D.F., Hudson S.A., Llorens P., Nanko K., 2017. Throughfall drop size distributions: A review and prospectus for future research. Wiley Interdisciplinary Reviews: Water, 4(4): e1225.

Levia D.F., Nanko K., Amasaki H., Giambelluca T.W., Hotta N., Iida S., Mudd R.G., Nullet M.A., Sakai N., Shinohara Y., Sun X., Suzuki M., Tanaka N., Tantasirin C., Yamada K., 2019. Throughfall partitioning by trees. Hydrological Processes, 33:1698–1708.

Li Y., Li J., Are K.S., Huang Z., Yu H., Zhang Q., 2019. Livestock grazing significantly accelerates soil erosion more than climate change in Qinghai-Tibet Plateau: Evidenced from 137Cs and 210Pbex measurements. Agriculture, Ecosystems & Environment, 285: 106643.

Lionello P., Scarascia, L., 2018. The relation between climate change in the Mediterranean region and global warming. Regional Environmental Change, 18: 1481–1493. 

Livesley S.J., Baudinette B., Glover D., 2014. Rainfall interception and stem flow by eucalypt street trees: The impacts of canopy density and bark type. Urban Forestry & Urban Greening, 13(1): 192-197. 

Llorens P., Domingo F., 2007. Rainfall partitioning by vegetation under Mediterranean conditions. A review of studies in Europe. Journal of Hydrology, 335, 37-54.

Llorens P., Poch R., Latron J., Gallart F., 1997. Rainfall interception by a Pinus sylvestris forest patch overgrown in a Mediterranean mountainous abandoned area. I. Monitoring design and results down to the event scale. Journal of Hydrology, 199: 331-345.

Llorens P., Domingo F., 2007. Rainfall partitioning by vegetation under Mediterranean conditions: A review of studies in Europe. Journal of Hydrology, 335(1-2): 37-54. 

Llorens P., Gallart F., Cayuela C., Planasdemunùt M.R., Casellas E., Molina A.J., Moreno De Las Heras M., Bertran G., Sánchez-costa E., Latron J., 2018. What have we learnt about Mediterranean catchment hydrology? 30 years observing hydrological processes in the Vallcebre research catchments. Geographical Research Letters, 44(2): 475-502. 

Ma C., Lib X., Luoa Y., Shaoa M., Jiaa X., 2019. Modelling of rainfall interception in growing and dormant seasons for a pine plantation and a black locust plantation in semi-arid Northwest China. Journal of Hydrology, 577: 123849.

Magliano P.N., Whitworth-Hulse J.I., BaldiG., 2019. Interception, throughfall and stemflow partition in drylands: Global synthesis and metaanalysis. Journal of Hydrology, 568: 638-645. 

Medail F., Quezel P., 2003. Conséquences écologiques possibles des changements climatiques sur la flore et la végétation du bassin méditerranéen. Bocconea, 16(1): 397-422.

Minďáš J., Bartík M., Škvareninová J., Repiský R., 2018. Functional effects of forest ecosystems on water cycle – Slovakia case study. Journal of Forest Science, 64: 331–339.

Mirzabaev A., Stringer L.C., Benjaminsen T.A., Gonzalez P., Harris R., Jafari M., Stevens N., Tirado C.M., Zakieldeen S., 2022. Cross-Chapter Paper 3: Deserts, Semiarid Areas and Desertification. In: Climate Change 2022: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change [H.-O. Pörtner, D.C. Roberts, M. Tignor, E.S. Poloczanska, K. Mintenbeck, A. Alegría, M. Craig, S. Langsdorf, S. Löschke, V. Möller, A. Okem, B. Rama (eds.)]. Cambridge University Press, Cambridge, UK and New York, NY, USA, pp. 2195–2231.

Missaoui K., Gharzouli R., Djellouli Y., Messner F., 2020. Phenological behavior of Atlas cedar (Cedrus atlantica) forest to snow and precipitation variability in Boutaleb and Babors Mountains, Algeria. Biodiversitas Journal of Biological Diversity,21(1): 239-245.

Molina A.J., del Campo A.D., 2012. The effects of experimental thinning on throughfall and stemflow: A contribution towards hydrology-oriented silviculture in Aleppo pine plantations. Forest Ecology and Management, 269: 206-213. 

Moreno-Pérez M.F., Pérez-Arellano R., Roldán-Cañas J., 2018. Influence of interannual rainfall variability on the interception process in a continental mediterranean climate. Revista de la Facultad de Ciencias Agrarias UNCuyo, 50(2): 139-154. 

Muzylo A., Llorens P., Valent F., Keizer J.J., Domingo F., Gash J.H.C., 2009. A review of rainfall interception modeling. Journal of Hydrology, 370(1-4): 191-206.

Navarro-Cerrillo R.M., Sarmoum M., Gazol A., Abdoun F., Camarero, J.J., 2019. The decline of Algerian Cedrus atlantica forests is driven by a climate shift towards drier conditions. Dendrochronologia, 55: 60–70.

Pérez-Suárez M., Arredondo-Moreno J.T., Huber-Sannwald E., Serna-Pérez A., 2013. Forest structure, species traits and rain characteristics influences on horizontal and vertical rainfall partitioning in a semiarid pine–oak forest from Central Mexico. Ecohydrology, 7(2): 532-543. 

Peñuelas J., Sardans J., Filella I., Estiarte M., Llusià J., Ogaya R., Carnicer J., Bartrons M., Rivas-Ubach A., Grau O., 2017. Impacts of global change on Mediterranean forests and their services. Forests, 8(12): 463. 

PriceA.G., Carlyle-Moses D.E., 2003. Measurement and modelling of growing-season canopy water fluxes in a mature mixed deciduous forest stand, southern Ontario, Canada. Agricultural and Forest Meteorology, 119(1-2): 69-85.

Rapp M., Romane F., 1968. Contribution à l'étude du bilan de l'eau dans les écosystèmes méditerranéens: Egouttement des précipitations sous des peuplements de Quercus ilex L. et de Pinus halepenis Mill. Oecologia Plantarum, 3: 271-284.

Rodrigo A., Avila A., 2001. Influence of sampling size in the estimation of mean throughfall in two Mediterranean holms oak forest. Journal of Hydrology, 243: 216‐227.

Sadeghi S.M.M., Van Stan J.T., Pypker T.G., Tamjidi J., Freisen J., Farahnaklangroudi M., 2018. Importance of transitional leaf states in canopy rainfall partitioning dynamics. European Journal of Forest Research, 137: 1-10.

Sadeghi S.M.M., Gordon D.A., Van Stan J.T., 2020. A global synthesis of throughfall and stemflow hydrometeorology. In: Van Stan J.T., Gutmann E., Friesen J. (Eds.): Precipitation Partitioning by Vegetation: A Global Synthesis. Springer Nature, Cham, pp. 49–70.

 Saito T., Matsuda H., Komatsu M., Xiang Y., Takahashi A., Shinohara Y., Otsuki K. 2013. Forest canopy interception loss exceeds wet canopy evaporation in Japanese cypress (Hinoki) and Japanese cedar (Sugi) plantations. Journal of Hydrology, 507: 287–299.

Schwingshackl C., Hirschi M., Seneviratne S.I., 2017. Quantifying spatiotemporal variations of soil moisture control on surface energy balance and near-surface air temperature. Journal of Climate, 30(18): 7105-7124.

Sensoy H., Tanyel M., 2022. Effect of heavy rain conditions on throughfall in evergreens and conifers in urban settings. Polish Journal of Environmental Studies, 31(1): 271-279.

Shi X., Bao Y., Liang Q., 2022. Applying the revised gash analytical model of precipitation interception to Pinus sylvestris var. mongolica in horqin sandy land, northeastern China. Fresenius Environmental Bulletin, 31(1): 572-582.

Shinohara Y., Ide J., Higashi N., Komatsu H., Kume T., Chiwa M., Otsuki K., 2010. Observation of canopy interception loss in an abandoned coniferous plantation. Nihon Ringakkai Shi/Journal of the Japanese Forestry Society, 92(1):54-59.

Slimani S., Derridj A., Gutierrez E., 2014. Ecological response of Cedrus atlantica to climate variability in the Massif of Guetiane (Algeria). Forest Systems, 23(3): 448-460.

Slimani S., Kherchouche D., Bekdouche F., Gutierrez E., 2021. Tree-ring reconstruction of March-June precipitation from the Atlas cedar forest of Mount Takoucht, Béjaïa (northern Algeria). Forest Systems, 30(3): e011–e011.

Sraj M., Brilly M., Mikos M., 2008. Rainfall interception by tow deciduous Mediterranean forests contrasting stature in Slovenia. Agricultural and Forest Meteorology, 148: 121-134.

Sun S., Xiang W., Ouyang S., Hu Y., Zhao Z., Xu X., Yue K., Chen L., Zeng Y., Lei P., Deng X., 2022. Higher canopy interception capacity of forests restored to the climax stage in subtropical China. Hydrological Processes, 36(3):p.e14538.

Sun X., Onda Y., Kato H., 2014. Incident rainfall partitioning and canopy interception modelling for an abandoned Japanese cypress stand. Journal of Forest Research, 19: 317-328.

Tague C., Moritz M., Hanan E., 2019.The changing water cycle: The eco-hydrologic impacts of forest density reduction in Mediterranean (seasonally dry) regions. Wiley Interdisciplinary Reviews: Water, 6(4):e1350.    

Touati L., HamelT., Meddad-Hamza A., de Bélair G., 2021. Analysis of rare and endemic flora in northeastern Algeria: the case of the wilaya of Souk Ahras. Bulletin de la Société Royale des Sciences de Liège, 90: 213–240.

Treder W., Klamkowski K., Tryngiel-Gać A.,Wójcik K., 2022. Assessment of rainfall efficiency in an apple orchard. Journal of Water and Land Development, 53: 51-57.

Tuel A., Eltahir E.A., 2020. Why is the Mediterranean a climate change hot spot? Journal of Climate, 33(14): 5829–5843. 

Vialard-Goudou A., Richard C., 1956. Étude pluviométrique, physico-chimique et économique des eaux de pluie à Saigon (1950-1954). Agronomie Tropicale, 11(1): 78-92.

Vicente E., Vilagrosa A., Ruiz-Yanetti S., Manrique-Alba A., González-Sanchís M., Moutahir H., Chirino E., del Campo A., Bellot J., 2018. Water balance of mediterranean Quercus ilex L. and Pinus halepensis Mill. forests in semiarid climates: A review in a climate change context. Forests, 9(7): 426.

Xiao Y., Chen L.H., Yu X.X., Yang X.B., Sun Q.Y., 2007. Influence on precipitation distribution of Pinus tabuleaefomis forest in Miyun Reservoir. Journal of soil and water conservation, 21(3): 154-157.

Yahi N., 2007. Les cédraies algériennes : phytoécologie et phytosociologie dynamique et conservation des peuplements. Thèse de Doctorat, USTHB, 218 pages.

Yahi N., Djellouli Y., 2010.Groupements forestiers et préforestiers à Cedrus atlantica Manetti d’Algérie: État des connaissances et dynamique des syntaxons. Revue Forestière Française, LXII (3-4):309-319.

Yang J., He Z., Feng J., Lin P., Du J., Guo L., Liu Y., Yan, J., 2023. Rainfall interception measurements and modeling in a semiarid evergreen spruce (Picea crassifolia) forest. Agricultural and Forest Meteorology, 328: 109257.

Yang T., Li Y., Ouyang X., Wang B., Ge X., Tang L., 2022. Rainfall physical partitioning and chemical characteristics in evergreen coniferous and deciduous broadleaved forest stands in a high nitrogen deposition region, China. Forests, 13(10): 1644.

Zhang Y., Yuan C., Chen N., Levia D.F., 2023. Rainfall partitioning by vegetation in China: A quantitative synthesis. Journal of Hydrology, 617: 128946.

Zhao W., Ji X., Jin B., Du Z., Zhang J., Jiao D., Zhao L., 2023. Experimental partitioning of rainfall into throughfall, stemflow and interception loss by Haloxylonammodendron, a dominant sand-stabilizing shrub in northwestern China. Science of the Total Environment, 858: 159928.

Zhou B., Han H.R., Kang F.F., Cheng X.Q., Song Y.L., Liu K., Li Y., 2013. Characteristics of precipitation distribution in Pinus tabulaeformis plantations under different canopy coverage in Taiyue Mountain. Acta Ecologica Sinica, 33:  1645-1653.