The phenological response of spring-flowering woody species to temperature has been widely approved and correlates very well with the air temperature (Spano et al., 1999; Menzel, 2002). After flower bud break, morphological development requires a certain thermal contribution to continue to develop after their latent phase (Orlandi et al., 2002; Anderson et al., 2005; Foley et al., 2009). In general, our results confirm a clear relationship between the beginning of pollination and air temperature.
It is commonly accepted that heat is only effective after sufficient chill has been met (Cannell and Smith, 1983; Fuchigami and Nee, 1987; Cesaraccio et al., 2004; Rea and Eccel, 2006). The models to estimate chill accumulation are under disagreement due to their variability among climates and the lack of physiological knowledge involved in their development (Luedeling, 2012). The choice of which type of model to use is often made a priori in phenology modeling studies, and models are then fitted to available data (Guo et al., 2013).
The chilling requirement based on “chill portions” has been considered as one of the best models for temperate nuts and fruits (Luedeling et al., 2009; Zhang and Taylor, 2011; Guo et al., 2013). This Dynamic Model has been found to be effective for other woody species (Xu et al. 2020), naming apricot (Zhuang et al., 2016; Benmoussa et al., 2017), olive (Aguilera et al., 2013), apple (Parkes et al., 2020), pistachio (Zhang and Taylor, 2011; Benmoussa et al., 2017), walnut (Luedeling et al., 2009), almond (Benmoussa et al., 2017) and cherry (Luedeling et al., 2013a). As a result, we kept the same parameter values proposed by Fishman et al. (1987) to calculate the chill portions of four taxa in this study. The Dynamic Model has often been found to be the most suitable among common chill models, especially in warm regions (Benmoussa et al., 2017).
The biothermic requirements that govern the different phenological development of the woody species in Mediterranean areas and estimated by different methodologies and threshold temperatures were the objective of numerous studies (Alba and de la Guardia, 1998; Galán et al., 2005; Bonofiglio et al., 2008; Aguilera et al., 2014, 2015; Achmakh et al., 2015). Although numerous methods for calculating GDD have been reported, agreement in estimating the best threshold temperatures has not received widespread acceptance (Ruml et al., 2010; Achmakh et al., 2015).The threshold temperature resulting effective in this research differs from those recorded for woody species in previous studies. The base temperatures reported ranged between 5 and 12.5°C for Olea in Andalusia region and Italy (Galán et al., 2005; Orlandi et al., 2006), between 0 and 25°C in Tunis (Aguilera et al., 2014) and between 4 and 12°C for Quercus in different climatic zones in Spain (Garcia-Mozo et al., 2000). Previously in the area of Tétouan, threshold temperatures of 10 and 7°C performed well at forecasting the initiation of pollen season and full anthesis of O. europaea, respectively (Achmakh et al., 2015). This fact can be considered to depend mainly on methodology used both to begin the heat accumulation and to statistically evaluate the best GDD threshold temperatures. Indeed, in the previous study undertaken in Tétouan, Achmakh et al. (2015) tested a range of base temperatures from 0 to 10°C. Aguilera et al. (2014) reported that the thresholds values decrease proportionally as the latitude increases, different threshold were obtained and ranged from 16°C in Perugia to 23°C in Zarzis. These results showed that the lower threshold or base temperature may vary between geographical areas. A relationship may exist between specific climatic conditions and taxa (Benmoussa et al., 2017). However, we have observed that the threshold temperatures are the same for the four taxa studied.
The peak of chilling, i.e. the day of the year with the maximum accumulation of positive chilling portions, has not been considered as the day to begin the heat accumulation in any taxa studied. This later fitted as the day after the end of chilling. This result suggests, despite arriving at the peak chilling day, the chilling requirements have not been met and the endodormancy phase release would need more days of chilling (Lang, 1987).
According to several studies, chilling and heat accumulation appeared discontinuous (Guo et al., 2013; Benmoussa et al., 2017) due to different physiological processes inside the bud that happen during different stages within the chilling and forcing periods. Dormancy initiation and release and plants sensitivity to chill are governed by genetical, genomic and physiological mechanisms (Ríos et al., 2014). Other factors such as endogenous signals, temperature, photoperiod and water availability regulate the flowering transition (Benmoussa et al., 2017).
Moreover, when buds chilling is met and conditions are favorable for budbreak, the active phenological development starts (De la Rosa et al., 2000). In this regard, the results of the present study demonstrate a good approximation. Here, for all studied taxa, the cumulative chilling decline after the peak of chilling provided the best start date to begin the heat accumulation period.
On these terms, a period of approximatively one month between the peak of chilling and the heat accumulation start date was detected for all taxa, accumulating the same amount of chilling. However, different GDD were required for their flowering onset. This finding might indicate that in woody species flowering in early spring, such as Morus, Pistacia and Quercus species, endodormancy release occurs earlier which is consequently related to their adaptative capacity to the environment concerned. However, for Olea, the endodormancy release seems to occur later. Indeed, Aguilera et al. (2014) suggested that a delay in the peak flowering date of olive observed in the Italian and Spanish areas might be provoked by a long transition between the end of endodormancy and the initiation of ecodormancy. Several studies have reported that, for the olive, low temperatures are needed only to break dormancy in previously initiated buds, as occurs in other fruit-tree species (Pinney and Polito, 1990; Rallo and Martin, 1991; Fernández-Escobar et al., 1992).
Summer temperature, winter chilling units and heat accumulation in woody species have been studied by numerous authors (Arnold, 1959; Frenguelli and Bricchi, 1998; Chuine et al., 1999; Garcia-Mozo et al., 2000; Galán et al., 2001) in order to clarify their involvement in the settlement and intensity of dormancy as well as cold acclimatisation. In the Tingitana peninsula, the coldest temperature is recorded during December-February (Aboulaich et al., 2013). Chilling requirements were not completely met under present climate conditions and insufficient in 2009 and 2010 for all populations, due to a warm winter. This result suggests that woody species in the southern of the Mediterranean sea need more heat than chilling showing their natural plasticity. However, the physiological effects of suboptimal chilling/forcing combinations should be investigated by horticultural researchers, which may have implications on crop yield or quality (Luedeling et al., 2013b).
The different results observed in chilling and heat requirements between flower and vegetative buds may be due to different environmental factors during the growing plant ‘life cycle (Rahemi and Pakkish, 2009). Thus chilling and heat requirement amounts vary between years to year depending on photoperiod, plant nutrition (Powell, 1987), temperature (Jacobs et al., 2002) and other environmental factors, may change the (Rahemi and Pakkish, 2009).
In our analysis, chill accumulation trend showed a significant increase for most of taxa studied over the past 12 years. This could be due to the significant decreasing trend of annual and monthly minimum temperatures observed in Tétouan in our previous study (Boullayali et al., 2021). Heat accumulation increased significantly with rising mean temperatures during the forcing period. However, during 2008–2019, heat accumulation trend of Morus, Olea, Pistacia and Quercus decreased by -0.89 to -0.19 GDD per year. These tendencies in heat requirements have been related to abundant chill accumulated, this conclusion was also confirmed by Alburquerque et al. (2008).
Corresponding with many other studies, heat accumulation is the main driver for flowering dates of woody taxa. Similarly, to this study, Guo et al. (2013) demonstrated in the cold winter climate of Beijing that forcing effects were much stronger than those of chilling on chestnut and jujube. Alonso et al. (2005) showed a similar result on almond in the cool climate of Zaragoza in northeast Spain. Sparks (1993) indicated that bud break of pecan could occur even when no chill is accumulated, and heat requirements are satisfied. In contrast, Egea et al. (2003) showed that in Spain the initiation of flowering of almond cultivars was influenced mostly by chilling requirements, with a little heat requirements contributing. This result was corroborated by other Studies concerning namely apricot, almond, pistachio and sweet cherry cultivars (Ruiz et al., 2007; Alburquerque et al., 2008; Rahemi and Pakkish, 2009; Campoy et al., 2012).
Regarding our results, it seems that warming due to increasing mean temperature during some parts of the dormancy season decreased chill accumulation rates, whereas the chill trend increased during the whole period studied. Since chilling requirements of woody taxa are met in all winters under present climate conditions of Tétouan, changes in winter chill may cause delays in spring phases in the near future. In contrast, impacts due to climate change are thus likely to arise from increased rates of heat accumulation rather than from changes in chill (Guo et al., 2013). These two figures may translate into delayed bud break and partial flowering compromising fruit production and resulting in low yield production.