Evaluation of historical and future thermal conditions for almond trees in north-eastern Portugal

The Trás-os-Montes agrarian region (TM; north-eastern Portugal) is a traditionally growing area of almond trees. However, climate change may significantly alter the current growing conditions and threaten sustainability. Chilling and forcing conditions in TM are assessed herein, also considering different varietal phenological timings. The dynamic model/chill portions (CP) and chilling hour (CH) models are used to assess the chilling phase. For the forcing phase, growing degree hours (GDH) and growing degree days (GDD) models are selected, hinting at differences between simpler (CH and GDD) and more complex models (CP and GDH). Furthermore, the climate change projections for these models are assessed. The models are computed using daily temperatures for baseline (Iberia01, 1971–2015) and two future periods (EURO-CORDEX: medium-term, 2041–2060, and long-term, 2061–2080), following two anthropogenic forcing scenarios (RCP4.5 and RCP8.5). For the future, a five-member ensemble of regional-global climate model chains is used. Projections for CP and CH show decreases, mostly for RCP8.5, whereas projections for GDH and GDD reveal an increase in heat conditions until the beginning of summer. GDH shows that potentially damaging high temperatures during summer, effectively lead to a decrease in heat accumulation for almond trees. CP and GDH, more biologically effective, seem to be more accurate in capturing climate change impacts. For the three varietal groups, the late cultivar may experience a greater reduction of CP and GDH compared to early and mid-season varieties. Nonetheless, the increase in heat forcing should be more detrimental than the decrease in chilling conditions. Although almond production may be compromised in the future, adequate adaptation options are suggested to mitigate future losses of yield/quality and warrant its sustainability.


Introduction
Almond species [Prunus dulcis Mill. (D.A Webb) syn. P. amygdalus Batsch] belongs to the family Rosaceae, sub-family Prunoideae and subgenre Amygdalus (L.) (AJAP 2017;Bento et al. 2020). Almond trees are native to warm and arid regions of southwest Asia and lower mountain slope regions of central Asia (Benmoussa et al. 2017;Oliveira et al. 2018). Over the centuries, the commercial interest in almonds increased, contributing to the growth of trade routes (Cabo and Matos 2017). Therefore, the almond expansion occurred in 3 main areas: Asia, America (especially in California) and the Mediterranean region. In 2020, the worldwide almond production was 3.9 × 10 6 t (FAOSTAT 2022). America is the continent with the largest almond production (2.2 × 10 6 t), followed by Asia (6.48 × 10 5 t), Europe (5.61 × 10 5 t), Africa (2.97 × 10 5 t) and Oceania (2.21 × 10 5 t) (Cabo and Matos 2017;FAOSTAT 2022). In America, California has been the largest producing region due to its suitable pedoclimatic conditions and structured commercial organization, supporting the demands of both national and international markets (Aguiar et al. 2017;Jin et al. 2020). In some Southern European regions, the prevailing Mediterranean-type climatic characteristics, with mild rainy winters and warm and dry summers, also benefit the almond tree development (NPR 2021). Spain (4.17 × 10 5 t), Italy (8.1 × 10 4 t) and Portugal (3.2 × 10 4 t) are indeed the largest producers in Europe.
In 2020, Portugal ranked 15th in the world in terms of production, with an almond cultivation area of 5.2 × 10 4 ha (FAOSTAT 2022). Besides the fact mainland Portugal has favourable climatic conditions for almond tree growing, the country has also benefited from the implementation of modern irrigation systems and the change from traditional to intensive cropping systems, particularly in the south of the country. Nonetheless, in the northern part of the country, the Trás-os-Montes (TM) agrarian region has a long tradition of almond production, which has led to a great diversity of varieties (Santos and Teixeira 2020;Doll et al. 2021). TM was the region that suffered fewer changes in traditional practices, remaining mostly rainfed, with emphasis on the preservation of native varieties (Doll et al. 2021). These varieties show different development and growth characteristics (Alonso et al. 2005;Sakar et al. 2019), resulting in differences in terms of productivity. The traditional varieties that are representative of this region are "José Dias", "Duro da Estrada Grado" and "Duro Amarelo Grado" (Oliveira et al. 2018;Santos and Teixeira 2020). Over the decades, other varieties have been cultivated, focusing on a higher resistance to climate change and pests/diseases, performances and productivity, such as the "Ferragnès", "Penta" and "Tardona" (Cordeiro et al. 2005).
Successful flowering and yield production require the fulfilment of dormancy, which is directly related to the air temperature (Díez-Palet et al. 2019). Dormancy can be divided into two main stages: endodormancy and ecodormancy (Prudencio et al. 2018b). The chill accumulation is associated with the endodormancy phase, when cool temperatures guarantee that plants retain the nutrients incorporated over the growing season, to promote floral differentiation and allow the branches to harden which makes them more resistant to frost (Díez-Palet et al. 2019;Guillamón et al. 2020;Djaman et al. 2021). Once the chill requirements are fulfilled, budbreak may occur (Paper et al. 2010;Prudencio et al. 2018a), and the endodormancy phase gives place to the ecodormancy phase, where heat accumulation is essential. Warm temperatures contribute to the sprouting and flowering buds, fruit ripening and vegetative growth also known as forcing phase (Tabuenca et al. 1972;Gaeta et al. 2018). Chill and heat requirements are species-and cultivar-specific (Prudencio et al. 2018b).
For the next decades, a significant regional temperature rise is predicted (Benmoussa et al. 2017;Salama et al. 2021;Freitas et al. 2022). These changes may have an impact Page 3 of 21 89 on the dormancy stages and endanger growth and physiological development (Santos et al. 2017). In case the chilling requirements are not fulfilled, it could be manifested with delayed budburst and crop growth, partial anthesis and fruit development may be threatened (Egea et al. 2003;Prudencio et al. 2018b;Djaman et al. 2021). For the forcing phase, warmer temperatures can lead to phenological advancements, unbalanced fruit maturation, enhancing stress conditions, causing disturbances in fruit production and bloom-advancing impacts Lorite et al. 2020). Moreover, higher temperatures throughout the growing season may also increase the frequency of occurrence of late spring frost, which can damage the flowers (Lorite et al. 2020). According to Guillamón et al. (2022), late frosts can affect 31-40% of the flowering of both extra-early and early varieties.
Several mathematical models have been used to calculate chill and heat accumulation. In general, these models are applied to different species and climatic conditions (Luedeling and Brown 2011). Moreover, these are useful for predicting the impacts of climate change and fundamental for outlining mitigation and adaptation measures (Salama et al. 2021). To determine chill conditions, the chilling hour (CH) and the dynamic (DM, henceforth CP) models are two of the most used models (Fishman et al. 1987a, b;Weinberger 1950). The former is a simpler model, widely used, while the latter is a complex model, which more accurately simulates the plant's physiological response to cold temperatures (more biologically effective). To calculate heat accumulation, the growing degree day (GDD) and the growing degree hour (GDH) models are also widely used, also representing a simple model and a more elaborated model, respectively (Anderson et al. 1986). The CP and the GDH use optimum temperature thresholds and compensation mechanisms to consider the effects of unfavourable extreme or unusually warm temperatures (Gaeta et al. 2018;Díez-Palet et al. 2019;Salama et al. 2021). Due to these factors, the CP and the GDH do not respond linearly to temperature, which is particularly important under future warmer climates.
In Portugal, studies considering the almond tree chill and heat conditions are still incipient, considering the high importance of almonds for the country's economy. As such, the present study aims to fill this knowledge gap. The objectives are four-fold: (1) to describe the climatic characteristics of the study area; (2) to understand the evolution of chill and heat accumulation under climate change; (3) to compare simple and complex models for chill and heat estimation; (4) to propose adaptative measures towards climate change risk reduction.

Study area
The TM agrarian region, located in north-eastern Portugal (Fig. 1a), is a territorial division with traditional cultural practices regulated by the Portuguese government. It integrates the agricultural subregions of "Alto Tâmega-Barroso", "Marão-Padrela", "Terra Fria", "Planalto Mirandês" and "Terra Quente" (Lima 1999). The almond tree is one of the species that stands out in the region, traditionally rainfed and mostly located in the "Terra Quente" (translates to "warm land") subregion. Herein, production is mostly limited by weather conditions and water availability (Doll et al. 2021 Office (INE 2022), almond production has been systematically recorded since 1986, reaching 1.7 × 10 4 t in 2021 (Fig. 1c). The characterization of the almond area was based on the digital inventory provided by the Geographical Institute of Portugal (Carta de Uso e Ocupação do Solo de Portugal Continental-COS 2007) (DGT 2018). Areas with the designation "Pomares de amendoeira", "Pomares de amendoeira com vinha" and "Pomares de amendoeira com olival" were selected as data source (Fig. 1b).
For future periods, a 5-member ensemble of the global circulation model-regional climate model (GCM-RCM) chain experiment was used (Table 1). This multi-model ensemble covers inter-model variability and reduces single-model uncertainties, explicitly in model setup, initializations and parameterizations (Santos et al. 2017;Freitas et al. 2022). Additionally, two future anthropogenic radiative forcing scenarios were taken into account: RCP4.5 (intermediate emission pathway) and RCP8.5 (high-emissions pathway) (IPCC 2019). These scenarios cover a wide range of uncertainties regarding future greenhouse gas emissions, atmospheric concentrations and climate change magnitudes, with global mean near-surface temperature increases until the end of the twenty-first century ranging from ~ 2 to ~ 4 °C, in RCP4.5 and RCP8.5, respectively (IPCC 2019). RCP4.5 corresponds to an additional (concerning pre-industrial levels) radiative forcing of 4.5 Wm −2 with stabilization after the mid-twenty-first century, while an additional radiative forcing of 8.5 Wm −2 is foreseen until the end of the twenty-first century for RCP8.5 (Thomson et al. 2011;Barredo et al. 2017).
A latitude-longitude sector was defined to optimize the study area, which contains vast continuous almond orchards (Fig. 1b). Gridded daily minimum and maximum air temperatures were extracted from the abovementioned datasets, for both baseline and future periods. Subsequently, the Iberia01 gridded temperatures were extracted for the recent-past period to calibrate the climate model variables with the "Empirical Quantile Mapping" bias-correction method (Cofiño et al. 2018). This methodology was critical to ensure that the bias between observed and modelled variables was correctly adjusted (Fraga et al. 2020;Fonseca et al. 2023). Subsequently, this methodology was then applied to future climate model data.

Climatic characterization
To characterize the climate of TM, the Iberia01 database (variables: daily mean temperature, in °C, and daily precipitation, in mm) was used to calculate specific parameters. The mean temperature for January (JanMT, in °C) and July (JulMT, in °C) were calculated, as these months correspond to the extremes of the mean seasonal cycle of temperature in the target region. The annual accumulated precipitation (AAP, in mm) and the ratio between the precipitation in the winter half and the summer half of the year (QWS, dimensionless) were also computed. QWS was a simple measure of precipitation seasonality. In TM, on average, from October (previous year) to March (next year) (winter half) there was much higher precipitation than from April to September (summer half), which was a typical feature of northern hemisphere Mediterranean-type climates. The regional orographic characteristics (elevation, in meters above mean sea level) were determined by a digital elevation model (DGT 2018). Table 1 Ensemble of global circulation model-regional climate model (GCM-RCM) chains selected for the present study (Kjellström et al. 2016;Kiesel et al. 2021;Larsson Ivanov et al. 2022) Global circulation model (GCM) Regional climate model (RCM) 89 Page 6 of 21

Estimating chill and heat accumulation
The chill and heat requirements of almond trees have been determined using several models. To analyse the differences between them, simple and more complex models were selected herein. In the present study, the chilling during winter dormancy was evaluated using the chilling hours (CH) model (Weinberger 1950) and the dynamic model (DM, henceforth CP) (Fishman et al. 1987a, b). The heat forcing during the growing season was calculated according to the growing degree days (GDD) and growing degree hours (GDH) models (Miranda et al. 2021;Miranda 2020;Anderson et al. 1986). According to Paper et al. (2010) and Weinberger (1950), the CH model is a simple and extensive system to calculate the chill conditions. The model quantifies hours between 0.0 and 7.2 °C (Luedeling and Brown 2011). However, CH does not recognize the adverse effect of warm temperatures on the chill accumulation period (Salama et al. 2021). The chillR package (version 0.72.7; Luedeling 2022) was used to apply the DM in R Core Team (2022), expressed in chill portions (CP) (Fishman et al. 1987a, b). The DM provided accurate estimates at various sites and conditions and takes into consideration the disruptive effect of intermittent higher temperatures in chill accumulation, which is particularly important in warm climate conditions (Díez-Palet et al. 2019). Its response curve is a function of hourly temperatures, with an optimum temperature of 4 °C . In the case of high temperatures, the model applies a neutralizing effect (Luedeling and Brown 2011). The hourly temperatures were estimated from daily data using the built-in package functions. The fruclimadapt package (version 0.4.4) was used to estimate GDD and GDH values (Miranda et al. 2021;Miranda 2020;Anderson et al. 1986). GDD and GDH were determined from the daily mean minimum and maximum temperatures. The GDD model considers temperatures above the base temperature (Tbase temperature = 4.5 °C), which was calculated as a daily temperature summation (Blanco-Ward et al. 2019). The GDH model was based on the assumption that heat accumulates when the hourly temperature is between a base temperature (4 °C) and the optimum temperature (25 °C), leading to a maximum heat accumulation (Benmoussa et al. 2017;Díez-Palet et al. 2019). Moreover, the model design takes into account warm temperatures with a critical temperature (36 °C). Temperatures above 36 °C cause slow and premature growth and may bring negative effects on the plant stomatal closure (Gu 2015;Santos et al. 2017).
A simulation was performed for the study area, considering the parameters highlighted above. Given the large spatial variability, some assumptions were introduced. For this reason, the chilling phase was estimated to start accumulating on October 1st (of the previous year) and end on February 28th (of the following year). The chilling requirement is often fulfilled in January, though it may shift to February under warm conditions (Santos et al. 2017). For heat accumulation, the period between January 1st and September 30th was considered, from flowering to fruit maturity (Sakar et al. 2019). These relatively long periods were adjusted to the thermal needs of several varieties, according to their different maturation periods, being slightly overlapped. The CP, GDH, CH and GDD were computed separately for each climatic dataset. Subsequently, for future climates, multi-model means were used, to take into account model uncertainty. Finally, the results of the present period were compared to the two future periods: mid-and long-term. For the sake of succinctness, the CH and GDD model results can be found in Supplementary Material.

Complementary analysis
As already explained, in the heat accumulation phase, temperatures above 36 °C can lead to heat stress, leading the GDH/day values to decrease. To understand these effects, the percentage of days where GDH/day values decrease after June 1st was calculated. In addition, the difference between historical and future scenarios/periods was computed to evaluate to what extent CP and GDH accumulated values will change in the future (increase/ decrease). These parameters are useful to growers because they indicate the possibility of insufficient chill and heat conditions during the crop lifetime (Fraga and Santos 2021).

Chill and heat conditions for almond varieties
The almond tree species has a large diversity of varieties and clones, with a distinct duration of phenological phases. Herein, we chose three representative groups of almond varieties, selected based on their flowering time classification (Table 2): early (represented by "Desmayo", "Duro da Estrada Grado" and "Duro Amarelo Grado"), mid-season ("Casa Nova", "Nonpareil" and "Ferraduel") and late ("Tarraco", "Penta" and "Tardona"). Therefore, we computed the amount of CP and GDH required by each group, based on their average phenological timings (Table 2), and compared the baseline with the future (midand long-term). In this study, we assume that the phenological timings (chilling and forcing phase and flowering) are the same in the present and future scenarios/periods.

Regional characterization
January is usually the coldest month of the year. According to the results, the JanMT (mean temperature for January) values range from < 2 to 9 °C (Fig. 2a). The region shows 40% of the areas that correlate with JanMT between 4.1 and 5 °C, standing out on a large area in the northeast of TM and 28% and 25% of areas corresponding to 3.1-4 °C and 5.1-6 °C classes, respectively. In the southern part of TM, temperatures are higher (6.1-9 °C), due to the influence of the Atlantic Ocean. Considering the grid box, previously defined for the almond tree distribution (Fig. 1b), the orchard is especially distributed in the southeast of TM. This means that the JanMT values corresponding to the almond distribution range between 4 and 7 °C.
In contrast, July temperatures are usually the highest of the year, and JulMT (mean temperature for July) values range from 16 to 26 °C for TM (Fig. 2b). According to JulMT, the warmer temperatures are observed in the southeast of the TM region (22.1-25 °C), which coincides with a vast area of almond trees. Overall, 47% of the total area is represented by 22.1-24 °C and 14% of the area corresponds to temperatures above 24 °C, in the southeast. Temperatures below 21 °C are seen in the west of the region. The physiological activity of almond trees can be strongly reduced at temperatures below 15 °C and above 35 °C (Bento et al. 2020).
Regarding AAP (annual accumulated precipitation) (Fig. 2c), for the region, the values range between 250 and 2000 mm. Most of the region (70%) corresponds to values between 501 and 1000 mm. The AAP values tend to decrease from the west (2000 mm) to the east (250 mm). According to the results, 38% and 34% of the area correspond to 501-750 mm and 751-1000 mm classes, respectively. The orchard area is characterized by AAP of 251-750 mm. Almond trees are very resistant to drought periods, being native to regions with 300-600 mm (Aguiar et al. 2017). The QWS (ratio between the precipitation in the winter half and the summer half of the year) values range from < 1.7 to 2.7 (Fig. 2d), which represents the precipitation seasonality. For this parameter, 34% of the TM corresponds to 1.91-2.1 QWS, in the southeast, whereas the highest values, from 2.61 to 2.7 QWS, stand out to the west. Apropos the almond orchard, QWS values vary between < 1.7 and 2.3, which means that the difference in precipitation between winter and summer is not as high as elsewhere in the region.
Elevation ranges roughly from < 250 to 2000 m (Fig. 3e). In total, 23% of the area corresponds to the elevation class 251-500 m, 44% corresponds to 501-750 m, 24% corresponds to 751-1000 m and 4% corresponds to < 250 m and 1001-1250 m, the other class are residuals. Lower elevation regions are found in the southwest. The mountainous areas  ; 1971-2015; dimensionless); e elevation (m). For representation purposes, a downscaling methodology was applied, following the work by Freitas et al. (2022) are in the north, depicting higher elevations. The almond tree areas correspond to elevations from 64 to 971 m (southeast).

Chill accumulation
The annual mean evolution of daily CP and CP accumulated for historical and future periods are demonstrated in Fig. 3  stabilizes at ~ 110 CP until the end of the season. For the future scenarios and periods (midand long-term), the CP data are analysed according to mean, minimum and maximum values attained from the 5 climate models (CP-model range henceforth). For the RCP4.5 scenario and the mid-term period (Fig. 3b), the end of the stabilization period of the daily CP is anticipated (December to January). During the chilling phase (October to February), the accumulated CP is reduced to 76 CP, a decrease of − 5 CP compared to the baseline (81 CP), which is only reached after 9 days (March 9th). For 2061-2080 (Fig. 3d), the accumulated value is 75 CP (− 6 CP) at the end of February, requiring 11 additional days (March 11th) to obtain the same CP as in the past. The CP reduction is more severe under RCP8.5 (stronger anthropogenic forcing) and more evident in the long-term period. In the chilling phase, the accumulated values of 75 CP and 69 CP correspond to 2041-2060 (Fig. 3c) and 2061-2080 (Fig. 3e), respectively. The CP/day stabilization periods are also reduced. For 2041-2060, to obtain the same 81 CP of the recent-past, 17 additional days are required (− 8 CP), while for 2061-2080, 33 days (April 2nd) are needed. According to the climate model means, in all future scenarios/periods the baseline value of 81 CP will not be met until the end of February. Nonetheless, regarding the climate model range for the accumulated CP, these values may be reached, taking into account the maximum CP values, while the model minimum can have very low values of CP. These very low CP values are particularly noticeable under RCP8.5 and more pronounced in the long-term period.
The results of the CH model (SM 1) indicate that the chill accumulation will be lower in the future compared to the recent-past period. SM 1a is following the baseline, the total number of chilling hours is 1375 CH. The CH/day values vary from 0 and 15 CH. For RCP4.5, 2041-2060 (SM 1b), during the chilling phase 1082 CH, a lower CH value compared recent-past period (− 293 CH). In the 2061-2080 period (SM 1d), the accumulated value obtained at the end of the chilling phase is 1079 CH (− 296 CH). For RCP8.5, the values are 1022 CH (− 353 CH) and 894 CH (− 481 CH), corresponding to the mid-(SM 1c) and long-term periods (SM 1e), respectively. Contrarily to the CP, the CH model does not consider the disruptive effect of a warm temperature period on the chill accumulation. Therefore, the results may not be as realistic as in the DM. Hence, the CH/day continues to increase from December to February.

Heat accumulation
The forcing phase follows the chilling phase and it is estimated between January 1st and September 30th, in agreement with the GDH model (Fig. 4), from flowering to fruit maturity. The evolution of GDH/day and accumulated GDH for the recent-past and future periods under both scenarios are illustrated in Fig. 4. For the baseline (Fig. 4a), the reference value is 5.8 × 10 4 GDH during the period from January to September. Overall, from January to June there is a continuous increase in daily GDH values (384 GDH), followed by a subtle negative slope until August, decreasing constantly afterward. The GDH/day curve shows two peaks, in June and September. These two peaks may correspond to the upper-temperature limit that is considered for GDH accumulation. In the RCP4.5 scenario, the accumulated GDH is 6.1 × 10 4 GDH for the two periods ( Fig. 4b and d). For RCP8.5, the increase in GDH is more pronounced (Fig. 4c and e) for the 2061-2080 period, 6.2 × 10 4 (Fig. 4e). Generally, in the future scenarios/periods, the accumulated GDH tends to increase, though differences appear to be small (compared to recent-past values), possibly due to the GDH temperature thresholds. From June to September, there is a negative curve when the GDH value reduces due to high temperatures (> 36 °C).
This effect is particularly visible in 2061-2080, RCP8.5. Nonetheless, all future scenarios/periods indicate that the reference value 5.8 × 10 4 GDH will be reached before the end of the forcing phase.
Concerning the results of the GDD model (SM 2), heat conditions show an increase in the future compared to the baseline. It should be noted that the GDD model does not have an upper-temperature threshold, such as the GDH. For the baseline, the daily GDD values range between 0 and 22 GDD, and the accumulated GDD reached 2665 GDD (SM 2a). For RCP4.5, in 2041-2060, during the forcing phase a value of 3010 GDD is achieved, higher when compared to recent-past (344 GDD).

Complementary analysis
In this analysis, we compare the percentage of days when the value is lower than June 1st (Fig. 5), which usually corresponds to the highest GDH value (1st peak date) (Fig. 4), for recent-past and future periods/scenarios. These results indicate that, from June to October, the temperature will rise significantly, resulting in lower GDH values, as temperatures rise above the critical threshold (36 °C) (Miranda 2020;Fraga et al. 2019). In the baseline, the maximum percentage of 64% is reached. For RCP4.5, the maximum is 81% and 89% for 2041-2060 and 2061-2080, respectively. For RCP8.5, the maximum is 79% and 99% for 2041-2060 and 2061-2080, correspondingly. As previously mentioned, this future percentage increase is associated with the limits set by the model when temperatures reach 36 °C.
The CP (Fig. 6a) and GDH (Fig. 6b) accumulated and the respective differences between each future period/RCP and historical are illustrated in Fig. 6. From October 1st (previous year) to September 30th (next year), the CP accumulated is calculated. From October until January, a sharp decrease is visible, with a stabilization from January until the end of February, followed by a strong decrease in CP until June, when a stabilization occurs. At the end of the chilling phase (end of February), anomalies are already observed, especially for the 2061-2080 period, with − 8 and − 12 for RCP4.5 and RCP8.5, respectively. Conversely, in the GDH case (Fig. 6b), future values tend to increase until June, decrease until September and then increase again. The decrease in the difference values, from June to September, is associated with the high summer temperature and the increase projected for the future. At the end of December, the anomalies reach their maximum value between 4.8 × 10 3 and 6.2 × 10 3 , for RCP4.5, 2041-2060and 2061-2080, and RCP8.5, 2041-2060. In RCP8.5, 2061-2080, the difference corresponds to 8.5 × 10 3 .

Evaluation of chill and heat conditions for flowering
The chill and heat conditions are now assessed taking into account the 3 varietal groups defined for the flowering stage: early, mid-season and late (Fig. 7). For the baseline, for budbreak the early varietal group shows 42 CP, the mid-season 63 CP and the late season 81 CP. For flowering, during the forcing phase, early varieties accumulate up to 2.5 × 10 3 GDH, mid 3.0 × 10 3 GDH and late 4.6 × 10 3 GDH. Concerning the RCP4.5 and regarding budbreak, a reduction in CP values occurs for all varietal groups and is more significant in the long-term. In general, compared to the baseline, in 2061-2080 a reduction of around − 8 CP can be observed. Conversely, for flowering, heat accumulation tends to increase, with Fig. 5 Percentage of days when the GDH value is lower than the GDH value on the June 1st (1st GDH peak) higher intensity for the late varieties, as they endure the warming later in the season. For RCP8.5, the most pertinent CP reduction is found in 2061-2080, where compared to the baseline, there is a more pronounced decrease in late (− 12 CP) than in early varieties (− 9 CP). In comparison with recent-past, the GDH increase is around 1.5 × 10 3 at the flowering period, for all varieties.  Fig. 7 Chill (CP) and heat (GDH) accumulation throughout the almond tree annual cycle, for the recent-past and RCP4.5 and RCP8.5 future scenarios (annual means for each period: 1971-2015, 2041-2060 and 2061-2080). The flowering periods for early, mid-season and late varieties are also represented according to Table 2 4 Discussion The constant search, by producers and scientific teams, for tools to identify the challenges and present solutions is increasing, and the evolution of technology has contributed to the effectiveness of this task. Climatic and phenological descriptions help the producers select the cultivar and the most effective agricultural practices for each situation (Valentini et al. 2001). As is known, the success of the chill and forcing phase is essential for the well-being of the plant and productivity. There is an inter-dependent relationship between the chill and forcing phase, e.g. the later the chill phase ends, the later the forcing phase begins (Díez-Palet et al. 2019). Adequate development of floral buds and fruits, efficient pollination process and fruit set depend on the thermal requirements being fulfilled (Prudencio et al. 2018b). In the future, the predicted climate change may have detrimental impacts on the chill and heat conditions, threatening the species' stability, and the almond tree sector as a whole.
The TM agrarian region, where winters are relatively cool and rainy while summers are hot and dry, is considered suitable for almond orchards, particularly in the southeast ("Terra Quente") (Bento et al. 2020). Chilling and forcing conditions were analysed for the almond trees in this region, for recent-past and future scenarios, also taking into account different varietal requirements. Given the results for the recent-past, the chill and heat conditions in the TM region may exhibit different behaviours in the future, chill accumulation tends to decrease and heat accumulation tends to increase, which may influence the prevalence of this crop.
Currently, this region shows CP values ranging from 75 to 95 CP, which is in agreement with Santos et al. (2017). For the future, the CP reduction is almost − 12 CP, which may have negative impacts. The chill reduction also applies to the results of the CH model. In the baseline, the value is 1375 CH and, in the future, the reduction is − 482 CH. It should be noted that the CH (− 35%) reduction is much stronger than the CP (− 15%). Even though the results express a chill reduction over the TM region under future climates, almond trees should still be able to fulfil their chilling requirements. However, if a stronger chill reduction indeed occurs, this could lead to complications, such as delayed phenological stages, asynchronous flowering and loss of flower buds and flowers (Egea et al. 2003;Prudencio et al. 2018b;Djaman et al. 2021;Tominaga et al. 2022). Moreover, it is fundamental to highlight the main difference between these chill models, the DM takes into account that in the case of high temperatures, the model applies a compensating effect, which is not the case with the CH model (Luedeling and Brown 2011). For this reason, it is believed that the DM is more accurate than the CH model, particularly regarding climate change projections.
Regarding forcing conditions, heat accumulation is projected to rise due to higher temperatures. According to the baseline, heat accumulation is 5.8 × 10 4 GDH. Considering the most severe scenario, the GDH will rise 2.8 × 10 3 and 3.5 × 10 3 , for mid-and long-term, respectively. In baseline, the GDD value is 2665 and the results show a rise of 430 GDD and 626 GDD, for periods mid-and long-term, in the RCP8.5. The increase in the GDD (+ 24%) is much more prominent than the GDH (+ 7%). This is due to a reduction in GDH during the June to September period (RCP8.5; long-term). Effectively, higher temperatures during this period (temperatures above 36 °C) affect the plant's stomatal closure (Gu 2015) and negatively influence fruit development, vegetative growth and flower differentiation which will influence the following year's yield (Cantin et al. 2020;Fraga and Santos 2021), and thus are discarded by the GDH model. Conversely, the GDD model does not establish a limit for maximum temperatures in summer. As the results suggest, regarding climate change projections, the use of the CP and GDH models, compared to CH and GDD, is recommended, as these models use statistical tools to reduce the unfavourable effect of warm temperatures (Gaeta et al. 2018;Salama et al. 2021). The most noticeable differences, between CP and GDH models, occur between the baseline and the RCP8.5 (longterm). According to the results and analysis, future simulations for the DM and GDH models could be useful tool in crop management and in forecasting climate change impacts, and growers should take this into account.
Regarding the linkage between the chilling and forcing models and almond phenology, 3 varietal groups, from the TM region, were analysed: early, mid-season and late varieties (Table 2). Almond trees require relatively low chilling to budbreak, already in late winter, well before other fruit trees (Aguiar et al. 2017;Prudencio et al. 2018b). For chill accumulation, in the baseline period, CP ranges from 42 to 81 (early to late seasons). In other studies, the CP values range from 8 to 55 CP (early to extra-late-41 to 106 days) (Benmoussa et al. 2017;Prudencio et al. 2018a, b;Díez-Palet et al. 2019;El Yaacoubi et al. 2019;Djaman et al. 2021). In future conditions (RPC8.5; long-term), a reduction in CP values is observed, being more significant for late varieties (− 12 CP). This reduction should not compromise the almond tree chilling requirements in the TM region, even under the most severe climate scenario. Nonetheless, a stronger reduction in chill may still occur. If this is the case, late varieties may experience greater impacts, possibly leading to more pronounced changes in phenology, such as delayed flowering, and other physiological disorders (Prudencio et al. 2018a, b). As an example, in south-eastern Spain, studies found that late varieties, such as Tardona and Penta, can have yield losses associated with chill requirements not being met (Prudencio et al. 2018a, b). Conversely, early and mid-season varieties should be less affected by a chill reduction, due to the lower chilling requirements.
From budbreak onwards, heat accumulation will increase under both future scenarios, with increased severity under RCP8.5. Regarding the varietal groups, late varieties will tend to experience this enhanced warming more than early or mid-season. This fast and intensive increase in heat accumulation may contribute to the advancement of flowering and maturity dates (Prudencio et al. 2018b). The occurrence of earlier flowering dates can increase the exposure of flowers and fruits to the frost damage, especially early flowering cultivars. Other effects of the temperature increase in the forcing phase are photoassimilates inhibition, disturbances in fruit formation with a smaller kernel, anticipated harvest and yield, and quality attributes reduction (Luedeling et al. 2013;Lorite et al. 2020). Studies report that for the heat accumulation up to flowering, GDH values can vary between 5.5 × 10 3 (e.g. Desmayo) and 9.0 × 10 3 (Penta) GDH (Egea et al. 2003;Prudencio et al. 2018b;Díez-Palet et al. 2019). Overall, the increase in temperatures during the season, which contribute to the decreased chilling and increased forcing will tend to shorten the duration of the phenological stages, which has already been reported by several studies .
It should be mentioned that the current study does not consider other factors that may also influence almond development and growth (Wolkovich et al. 2022). The impact of the increase in pest and disease pressure (Aguiar et al. 2017). The impact of spring frost, which can also significantly affect almond productivity, especially in cold regions, such as TM (Segara and Ossama 2007;Lorite et al. 2020). Early-bloom cultivars are more vulnerable to frost impacts compared to late-bloom varieties (Guillamón et al. 2022). Moreover, other climatic variables could affect the species' growth, such as the decreased annual precipitation that may increase drought stress (Doll et al. 2021). Additionally, this current study provides a simulation based on climatic factors that are prone to several uncertainty sources, such as model and scenario uncertainty (Bedia et al. 2013). Furthermore, we assume that the phenological timings (chilling and forcing phase and flowering) are the same in the present and future scenarios/periods, which may not occur.
Adaption measures may be used to reduce potential yield/quality losses and ensure the sector's long-term viability. To address the lack of chill, growers can take advantage of the application of bio-stimulants or even adopt more suitable scion-rootstock combinations (Prudencio et al. 2018a, b). To decrease heat accumulation, the implementation of sunscreens or shading systems is recommended . In the case of almond species, growers focus on developing varieties that are more resistant to prevent the effects of spring frosts. Earliest varieties are more sensitive to frost events, since the flowers are already in development, and can be easily damaged by these extreme events (Guillamón et al. 2022). Contrarily, late varieties present good adaptability to different climatic conditions, which show naturally more resistance to this event, owing to the later flowering period (Martínez-Gómez et al. 2017). To reduce the decreased precipitation impacts, in a sustainable manner, the implementation of irrigation systems can be an option that helps during the hottest and driest months, or crop relocation (Bento et al. 2020).
Heatwaves can also be harmful to the species, and to reduce their impact select more heat-resistant varieties and apply for fertilizer and bio-stimulating plant protection under environmental stress are proposed (Freitas et al. 2021). To combat the other climate impacts, it is recommended to cultivate varieties that are most suited to the local climate and are more resilient to climate change, as well as to employ modern breeding techniques, like mapping the genetic basis of chilling requirements, to improve the viability and productivity of almond trees (Luedeling 2012). Additionally, the support for the breeding of new enhanced varieties for predetermined chilling and heat condition values is essential for climate change risk reduction (Guillamón et al. 2022). To address these issues, additional studies are required, such as calculating the occurrence of extreme weather events (e.g. spring late frost and heatwaves), precipitation irregularity over the year, water supply in the future and the respective impacts they have on the almond trees.

Conclusions
Climate challenges are being experienced with each passing day. According to the climatic projections, the TM agrarian region, defined by traditional agricultural practices and almond varieties, may be more susceptible to significant climate change impacts. As the temperature is expected to increase, chill and heat conditions tend to change. As seen in this study, increased temperature during the forcing phase should be more problematic than during the chilling phase. As can be seen, the use and development of mathematical models are fundamental, since they are auxiliary tools for crop prediction and management. These models can also be considered as a decision-support tool for farmers and agricultural consultants, namely, concerning the selection of varieties and cultural and management practices. Adaptation measures must be timely planned and implemented so that the sector does not suffer economic and biodiversity losses. Therefore, investment should be made in more resilient and sustainable varieties and more sophisticated and adapted agricultural practices to adapt to climate change, reducing the negative impacts. Lastly, the present study contributes to a better understanding of the evaluation of thermal conditions for almond species in the TM agrarian region. Furthermore, since the orographic and climatic characteristics of the TM region may be similar to other regions of the world, this study may be used as a reference. It is hoped that this study will be useful in almond management and inspire further work.