3.1 Variation in ambient climatic parameters
Ambient temperature and photosynthetically active radiation (PAR) varied significantly by month and phytogeographical zone (P < 0.05) (Fig. 2a). The highest temperatures (41.48°C ± 0.22) and instantaneous PAR values (369.10 ± 70.84 µmol photons m-2s−1) were recorded during the dry months (April, May October, November, and December) in both sites (Fig. 2). In contrast, from June to September, ambient temperature dropped and PAR showed considerable variability. The lowest temperature (32.63°C ± 0.21) was recorded in July, while the lowest PAR (96.28 ± 26.37 µmol photons m− 2s− 1) occurred in August, coinciding with peak of the rainy season. Cloud cover during this period strongly influenced PAR levels. A significant interaction between months and phytogeographic zone was also observed (P < 0.001), indicating that the temporal patterns of temperature and PAR differed between sites.
Statistical analysis revealed a significant difference in VPD according to month, phytogeographical zone and the interaction between month and phytogeographical zone (P < 0.05) (Fig. 2b). From June to September, the presence of moisture in the air and the drop in temperatures led to a decrease in the VPD. The VPD remained low from June to September (Fig. 2b), indicating a very low capacity of the air to absorb water during this period. After September, the VPD increased and reached 5.7 ± 0.07 kPa in December in both phytogeographical zones.
3.2. Variation of Shea leaf temperature and thickness
Leaf temperature values revealed a significant difference depending on the month and phytogeographic zone (P < 0.001) (Fig. 3a). However, no significant differences were found depending on tree diameter class (P = 0.15). Conversely, significant interactions were revealed between diameter class and month (P < 0.001), and between diameter class and phytogeographic zone (P < 0.001). The highest leaf temperatures were observed in both phytogeographic zones in April, and the lowest temperatures in July, namely 41.45 ± 0.27°C and 32.73 ± 0.15°C in the NSZ, and 41.10 ± 0.25°C and 31.55 ± 0.23°C in the SSZ, respectively. Low leaf temperatures were recorded in both phytogeographical zones during the months of June, July, August and September, while high temperatures were recorded in April (Fig. 3a). The [15;30[ diameter class recorded the lowest leaf temperature (35.78 ± 0.23°C) in the South Sudanese phytogeographical zone, compared to the highest leaf temperature (37.27 ± 0.24°C) recorded in the North Sudanese phytogeographical zone.
Leaf thickness differed significantly according to phytogeographical zone (P < 0.001) and month of the year (P < 0.001). However, no significant difference in leaf thickness was observed according to tree diameter class (P = 0.06). Conversely, significant interactions of diameter class*month of the year and phytogeographical zone, as well as month of the year*diameter class*phytogeographical zone were revealed (P < 0.001). Leaves from the northern Sudan zone were thicker (0.652 mm ± 0.015) than those from the southern Sudan zone (0.600 mm ± 0.015). A difference in leaf thickness was observed at the start of leaf expansion in early April; however, this difference was no longer evident from August until leaf fall in December (Fig. 3b). A decrease in leaf thickness was observed from June to October in both phytogeographical zones. Similarly, the diameter class [30; 45] recorded the lowest values of leaf thickness (0.53mm ± 0.028) in the southern Sudan phytogeographic zone, while the highest(0.68 mm ± 0.027) was recorded in the northern Sudan phytogeographic zone in the diameter class ≥ 45.
3.3. Variation of shea photosynthetic parameters
Values of relative chlorophyll content differed significantly for month and phytogeographical zone (P < 0.001), but not for diameter class (P = 0.98). Furthermore, significant interactions were found between diameter class and month, and between month, diameter class and phytogeographical zone (P < 0.001). At the beginning of bud break in April, the relative chlorophyll content was similar in the two phytogeographical zones. Over time, however, a difference emerged. Specifically, the mean chlorophyll content of shea leaves from the SSZ was higher (42.9 ± 0.4 SPAD) than that of leaves from the NSZ (40.8 ± 0.37 SPAD). The maximum chlorophyll content was recorded in July and August with values 47.28 ± 0.47 and 48.46 ± 0.47 SPAD, respectively. Relatively low chlorophyll content was recorded in April (33.53 ± 0.77 SPAD) and December (34.54 ± 0.94 SPAD). In both phytogeographic zones, maximum chlorophyll content was recorded in August: 47.17 ± 0.68 SPAD in the Northern Sudan phytogeographic zone, and 49.78 ± 0.56 SPAD in the Southern Sudan (Fig. 4a). The lowest relative chlorophyll content (39.07 ± 0.70 SPAD) was recorded in the North Sudanian phytogeographic zone, while the highest value (44.89 ± 0.75 SPAD) in the South Sudanian phytogeographic zone both for the ≥ 45 diameter class (Fig. 4b).
There was a significant difference in photosystem II yield components as a function of month (P < 0.05). Similarly, there was a significant difference between phytogeographic zones for φNO (P = 0.04) and φNPQ (P = 0.001), but not for φII (P = 0.20) (Fig. 5). Tree diameter class did not show significant effect for any of these components (P > 0.05). Significant interaction between diameter class and photogeographic zone was observed for only φII (P = 0.002) as well as between month, diameter class and photogeographic zone for all components (P < 0.001).
The φ II varied from 0.57 ± 0.02 to 0.65 ± 0.01 and from 0.54 ± 0.01 to 0.65 ± 0.001, respectively, between April and December in the SSZ and NSZ (Fig. 5a). Conversely, the non-photochemical quenching efficiency (φNPQ) increased in the opposite direction to the photosystem II effective quantum yield (φII) (Fig. 5b). The lowest values of (φNO) were recorded in April in both phytogeographic zones (Fig. 5c). The φNO and φNPQ were 0.25 ± 0.002 and 0.13 ± 0.005 respectively in the southern Sudanian and northern Sudanian phytogeographic zones, compared with 0.24 ± 0.006 and 0.13 ± 0.005 respectively in the southern Sudanian and northern Sudanian phytogeographic zones.
Significant monthly differences were revealed in non-photochemical quenching (NPQt) and linear electron flux (LEF) (P < 0.001). Similarly, the test revealed significant differences in NPQt between phytogeographic zone (P < 0.001), while no significant differences were observed for LEF (P = 0.42) as well as in linear electron flow (P = 0.91) or NPQt (P = 0.31) according to Shea tree diameter class.. No significant interactions were found between diameter class and phytogeographic zone for LEF (P = 0.11) and for NPQt (P = 0.68). However, significant interactions were revealed between diameter class and month, and between month, diameter class and phytogeographic zone (P < 0.001) for NPQt and LEF.
Non-photochemical quenching was higher in the North Sudanese phytogeographic zone (0.78 ± 0.08) than in the South Sudanese zone (0.73 ± 0.06). The highest NPQt value was recorded in April (1.42 ± 0.20), compared to 0.35 ± 0.03 in July. Low NPQt values were recorded in June (0.37 ± 0.04) and July (0.031 ± 0.03), while high values were recorded in December (1.68 ± 0.48) and April (1.70 ± 0.35) in the South and North Sudanian phytogeographic zones respectively (Fig. 6a). Linear electron flow followed the same trend, with the highest value observed in April (70.11 ± 6.42 µmol electrons m⁻² s⁻¹) compared to July (28.26 ± 3.73 µmol electrons m⁻² s⁻¹). Similarly, low LEF values were recorded in August (23.23 ± 4.10 µmol electrons m⁻² s⁻¹) and July (27.23 ± 4.65 µmol electrons m⁻² s⁻¹), while high values were recorded in April (66.83 ± 7.71 µmol electrons m⁻² s⁻¹) and (73.46 ± 10.27 µmol electrons m⁻² s⁻¹) for the southern and northern Sudanian phytogeographic zones, respectively (Fig. 6b)
3.4. Relationship between climatic variables and shea photosynthetic parameters
The regression between leaf temperature and PAR was significant (P < 0.001) and positive across phytogeographic zones. Regression coefficients remained low, with values of 0.08 and 0.06 in the northern and southern Sudanian phytogeographic zones, respectively .
A significant (P < 0.001) and positive regression was observed in the relationship between ambient temperature and leaf temperature, depending on the phytogeographic zone. The coefficient of determination for leaf temperature and ambient temperature indicates that these two parameters are closely related. Indeed, 84% of leaf temperature is related to variations in ambient temperature in the northern Sudanese phytogeographic zone and 78% in the southern Sudanese phytogeographic zone (Fig. 7a). An increase in ambient temperature leads to an increase in leaf temperature (Fig. 7a).
The linear regression between leaf temperature and vapour pressure deficit was significant (P < 0.001) and positive depending on the phytogeographic zone. The regression coefficients were 0.75 and 0.59 respectively in the North Sudanese and South Sudanese phytogeographic zones.
The coefficient of determination of the linear regression was 0.89 between LEF and PAR in NSZ and 0,86 in SSZ (Fig. 8). Although the coefficient of determination of LEF was identical in the two phytogeographical zones, the half saturation point was different. The half-saturation point was 137.54 µmol photons m− 2s− 1 and 175.60 µmol photons m− 2s− 1 for PAR in the southern Sudanian and northern Sudanian phytogeographical zones, respectively. A significant (P < 0.01) and positive regression was observed between LEF and VPD on the one hand, and between LEF and ambient temperature according to phytogeographic zone on the other hand. However, the regression coefficients remained low and varied between 0.02 and 0.05 depending on the phytogeographic zone.
The regression was positive and significant (P < 0.001) between NPQt and PAR according to phytogeographic zone. The regression coefficient was 0.25 for NPQt in the Southern Sudanian phytogeographical zone and 0.45 in the Northern Sudan phytogeographical zone (Fig. 9). A significant (P < 0.001) and positive regression was observed between NPQt and VPD, as well as between NPQt and ambient temperature, depending on the phytogeographic zone. However, regression coefficients remained low, varying between 0.04 and 0.06 depending on the phytogeographic zone
The regression coefficients of the photosystem II yield components show that the photochemical efficiency of photosystem II (φII ) and the quantum yield of unregulated non-photochemical energy loss in photosystem II (φNO) decrease with increasing photosynthetically active radiation (PAR) (Fig. 10). The linear regression was negative and significant (P < 0.001) for PAR, φII and φNO, depending on the phytogeographic zone (Fig. 10a, 10c). The respective regression coefficients are 0.49 and 0.57 for φII and 0.31 and 0.07 for φNO in the northern and southern phytogeographic zones of Sudan. The quantum yield of regulated thermal dissipation in photosystem II (φNPQ) showed a positive, significant regression (P < 0.001) with PAR: φNPQ increased with increasing PAR (Fig. 10b), with R² values of 0.78 and 0.55 for the northern and southern Sudan phytogeographic zones, respectively. Almost 50% of the variation in φII and φNPQ can be explained by PAR (Fig. 10b).
Significant negative relationships (P < 0.001) were observed between φII and VPD, and between φNO and VPD, depending on the phytogeographic zone. However, regression coefficients remained low, ranging from 0.04 to 0.13 depending on the phytogeographic zone. As for φNPQ, the regression was positive and significant (P < 0.001) with respect to vapour pressure deficit (VPD), depending on the phytogeographic zone. Similarly, regression coefficients were low, at 0.06 for the NSZ and 0.05 for the SSZ.
Significant negative relationships (P < 0.001) were observed between φII and ambient temperature, and between φNO and ambient temperature, depending on the phytogeographic zone. However, regression coefficients remained low, ranging from 0.01 to 0.18 depending on the phytogeographic zone. The regression of ambient temperature and φNO had R² values of 0.09 and 0.18 in the NSZ and SSZ, respectively. As for φNPQ, the regression was positive and significant (P < 0.001) with respect to vapour pressure deficit (VPD), depending on the phytogeographic zone. Similarly, the R² regression coefficients were low, at 0.08 for the NSZ and 0.07 for the SSZ.
The regression was positive and significant (P < 0.001) between stomatal conductance and VPD according to phytogeographic zone. The coefficient of determination R² of the linear regression of stomatal conductance and VPD shows that 26% of stomatal conductance is explained by VPD in the northern Sudanian phytogeographical zone, whereas it is explained by 44% in the southern Sudanian phytogeographical zone. A higher sensitivity of stomatal conductance to VPD was observed in the southern Sudanian phytogeographical zone (slope of the equation = 190 mol m− 2s− 1) compared to the northern Sudanian phytogeographical zone (slope of the equation = 140 mol m− 2s− 1) (Fig. 11).
The linear regression between PAR and stomatal conductance was positive and significant (P = 0.015) in the North Sudanese phytogeographic zone, but not significant in the South Sudanese phytogeographic zone. (R² = 0.03).
The regression was positive and significant (P < 0.01) between ambient temperature and stomatal conductance according to phytogeographic zone with low regression coefficient (R²=0.03) in both zones.