Leaf chlorophyll estimates of temperate deciduous shrubs during autumn senescence using a SPAD-502 meter and calibration with extracted chlorophyll

Abstract

Key message

The soil plant analysis development (SPAD) meter proved effective at estimating leaf Chl content of native and non-native temperate deciduous shrubs. However, due to the change in the relationship between SPAD values and Chl content during the course of autumn senescence, it may be necessary to establish separate calibration equations to cover a range of leaf Chl concentrations.

Context

Non-destructive estimates of leaf chlorophyll (Chl) content by hand-held chlorophyll meters such as the Minolta SPAD-502 provide an effective time-efficient method of collecting field data on senescence.

Aims

To establish useable species-specific calibration equations of relationships between chlorophyll meter readings and leaf Chl content.

Methods

During autumn, we collected SPAD values and quantum yield of photosystem-II (Φ_PSII) over a 10-week period from 1044 leaf samples representing native and non-native temperate deciduous shrub species growing in the wild. Subsequently, Chl was extracted from leaf discs (1 cm diameter), incubated at 65 °C for 2 h in dimethyl sulfoxide, and quantified by spectroscopy. Relationships between extracted Chl and SPAD values were established using linear, quadratic, and exponential equations.

Results

Quadratic functions proved the most reliable fit for the data. Interestingly, non-native species tended to have higher leaf Chl content and SPAD values but also exhibited higher variability than native species. The strength (r²) of the SPAD-Chl relationships was weaker than those reported for agricultural or greenhouse-grown plants, but strengthened over time as Chl, Φ_PSII, and SPAD values declined.

Conclusion

The SPAD meter proved effective at estimating leaf Chl content of native and non-native temperate deciduous shrubs. These estimates may be used in future to validate remote sensing (satellite and near-surface) vegetation indices in late autumn when trees are leafless.

1 Introduction

Despite the pivotal role that shrubs play in temperate deciduous forests by, for example, contributing to carbon and water exchange, nutrient cycling and providing food and habitat for a range of organisms, they have tended to be under-researched relative to trees (Donnelly and Yu 2017, 2019). Shrubs generally leaf out earlier than trees (Harrington et al. 1989; Gill et al. 1998) and non-native shrubs, in particular, often remain photosynthetically active after trees become leafless (Resasco et al. 2007; Wilfong et al. 2009), thus extending the C uptake period of forests. Therefore, the contribution of shrubs to forest productivity and C exchange is particularly important at the extremes of the growing season. Tracking leaf coloration in autumn is often achieved by direct in situ observations (Menzel et al. 2006; Yu et al. 2016; Donnelly et al. 2018) typically of individual forest trees or more recently over much larger geographical scales using satellite-derived phenological metrics extracted from NDVI (normalized difference vegetation index) and EVI (enhanced vegetation index) time series (Zhang and Goldberg 2011; Zhang et al. 2012; Liu et al. 2017). However, an increasing number of reports of discrepancies in the timing of coloration has emerged with the start of satellite-derived leaf color tending to be earlier than direct in situ observations in urban (Donnelly et al. 2018) and forested (Zhao et al. 2020) locations. There are potentially many explanations for these differences, including geographical scale, species composition, and sampling technique. However, another possibility may be that during autumn senescence, chlorophyll (Chl) begins to degrade before visible symptoms are evident but which may be detectable by satellite sensors which are sensitive to small changes in leaf spectral properties. In order to understand the autumn phenology of shrubs under natural conditions, it is necessary to consider different species and in particular to quantify declines in chlorophyll content and variation in autumn senescence between native and non-native shrubs. One way to examine in situ autumn leaf senescence of shrubs is by using a non-destructive chlorophyll meter.

The use of a hand-held chlorophyll meter to estimate foliar Chl content is a cost effective and convenient method of tracking leaf senescence during the autumn season and different meters have been used across a broad range of species (reviewed by Parry et al. 2014 and Table 1). A large number of readings may be achieved in a relatively short period of time without having to remove or destroy leaf tissue for Chl extraction and laboratory analysis. However, in order to ensure accurate and reliable relationships between chlorophyll meter readings and actual leaf Chl content, it is necessary to establish species-specific calibration equations (Markwell et al. 1995; Jifon et al. 2005; Netto et al. 2005; Pinkart et al. 2006; Hawkins et al. 2009; Coste et al. 2010).

Table 1 Sample publications reporting use of non-destructive methods for chlorophyll content determination from a range of agricultural, horticultural, forestry, and other plant species.

Non-destructive optical methods of estimating leaf Chl content using instruments such as either Minolta soil plant analysis development (SPAD) or Opti-Sciences Chlorophyll Content Meter (CCM) have been widely used in agriculture (wheat, rice, cotton, sorghum, soybean, maize, pea, potato, etc.), horticulture (tomato, apple, variety of citrus species, coffee, turfgrass), and viticulture (grapevine) (Table 1). In addition, similar studies have been carried out on (i) temperate tree species either growing in the wild or in plantations (oak, sycamore, maple, birch, beech, eucalyptus), (ii) deciduous shrubs and lianas grown in a common garden, and (iii) a range of tropical and sub-tropical fruit tree species (Table 1). The underlying principle of this method is based on absorption of red light by Chl. The SPAD instrument compares absorbance in the red (~ 650 nm, peak Chl absorbance) and near-infrared (~ 940 nm non-Chl absorbance) parts of the spectrum with the difference being proportional to foliar Chl content (Markwell et al. 1995). The main application of this approach has been in assessing physiological condition, photosynthetic potential, health/nutrient status, senescence rates, and stress levels of target plants for both research and management purposes. Together, these studies clearly demonstrate the reliability of hand-held meters at estimating foliar Chl content for a wide variety of species and applications over a range of environmental conditions. While other leaf pigments such as carotenoids play an important role in autumn senescence, hand-held chlorophyll meters such as SPAD and CCM target wavelengths other than those absorbed by carotenoids and are thus not suitable for examining leaf carotenoids (Peñuelas and Filella 1998). Carotenoid content may show strong relationships with SPAD but only because leaf Chl and carotenoid content are typically highly correlated (e.g., Keskitalo et al. 2005).

Despite the wide range of species for which relationships between non-destructive Chl readings and leaf Chl content have been established (Table 1), including for a number of forest trees, we found only two reports relating to a wetland shrub species (Hawkins et al. 2009) and none relating specifically to temperate deciduous shrubs growing in the wild. In fact, Coste et al. (2010) noted that portable chlorophyll meters were underutilized in forest science while Donnelly and Yu (2017, 2019) reported that shrubs tended to be overlooked in forest research. Therefore, to address this knowledge gap, the aim of the current research was to construct calibration equations to quantify the relationship between chlorophyll meter readings (Minolta SPAD-502) and extracted leaf Chl content for a range of native and non-native temperate deciduous shrub species growing in the wild. We aimed to test two hypotheses (i) the SPAD chlorophyll meter is a robust method to estimate leaf Chl content of temperate deciduous shrubs and (ii) the relationship between SPAD and Chl changes during the autumn senescence process.

2 Materials and methods

2.1 Shrub species, SPAD meter readings, and fluorescence

Leaf color change of eight shrub species representative of the understory layer of a small (11 ha) urban woodland on the campus of the University of Wisconsin-Milwaukee, USA, was monitored during autumn 2018. Shrubs were categorized into native (Ribes americanum Mill., wild currant; Prunus virginiana L., chokecherry; Viburnum acerifolium L., maple leaf viburnum; Viburnum lentago L., nannyberry; Cornus alternifolia L., pagoda dogwood) and non-native (Lonicera morrowii A. Gray, honeysuckle; Rhamnus cathartica L., common buckthorn; Ligustrum vulgare L., European privet) species.

SPAD values (Konica-Minolta, Japan, SPAD-502) and Chl a fluorescence measurements were recorded twice weekly from DOY 250 to 296 on all species and subsequently on two occasions (DOY 302 and 312) when leaves were present only on non-native species. Quantum yield of photosystem II, Φ_PSII (Maxwell and Johnson 2000), was measured with a Diving-PAM fluorometer (Walz GmbH, Effeltrich, Germany) following Bott et al. (2008). On each occasion SPAD and Φ_PSII readings were recorded at the same location on each of three leaves from three individual plants of each species (nine readings per method per species). All shrubs were shaded by a tall canopy and readings were made between 9:00 and 10:30 am each day, so leaves were exposed only to low light levels and no direct sunlight prior to recording Chl a fluorescence and SPAD values. Diving-PAM excitation light intensity and gain settings were kept constant throughout the whole study to ensure directly comparable Chl a fluorescence values. The same leaves (or leaf part) were subsequently removed and placed in individual labelled plastic bags prior to freezing at − 28 °C for Chl extraction at a later date. All raw data are publically available (Donnelly et al. 2019).

2.2 Chlorophyll extraction

Leaf samples were removed from the freezer and discs, of 1 cm in diameter, were cut from each using a standard borer. Each leaf disc was placed in a 15-ml polyethylene centrifuge tube prior to adding 10 ml dimethyl sulfoxide (DMSO) and the lid secured. The tubes were placed in a water bath and incubated at 65 °C (Hiscox and Israelstam 1979) for 2 h. Preliminary tests indicated that 30 min (Hiscox and Israelstam 1979) or 1 h was insufficient time to render the leaf discs, for all species, colorless. Furthermore, use of DMSO in this study was justified as extracts did not exhibit a brown color like some tree leaf samples reported by Minocha et al. (2009). Samples were removed from the water bath and a 3 ml extract of each supernatant was pipetted individually into a quartz cuvette for spectrophotometer (Beckman-Coulter DU-640) analysis, using a DMSO blank prior to recording absorbance of all leaf samples at 750, 663, 645, and 480 nm.

Total Chl, Chl a, and Chl b concentrations were calculated following the widely (e.g., Richardson et al. 2002; Netto et al. 2005) used equations reported in Arnon (1949), while carotenoid concentration was determined following Wellburn (1994):

$$ \mathrm{Total}\ \mathrm{Chl}\ \left(\mathrm{g}\ {\mathrm{l}}^{-1}\right)=0.0202\ {A}_{645}+0.00802\ {A}_{663} $$

$$ \mathrm{Chl}\ a\ \left(\mathrm{g}\ {\mathrm{l}}^{-1}\right)=0.0127\ {A}_{663}+0.00269\ {A}_{645} $$

$$ \mathrm{Chl}\ b\ \left(\mathrm{g}\ {\mathrm{l}}^{-1}\right)=0.0229\ {A}_{645}+0.00468\ {A}_{663} $$

$$ \mathrm{Total}\ \mathrm{carotenoids}\ \left(\mathrm{mg}\ {\mathrm{l}}^{-1}\right)=\left(1000\ {A}_{480}-2.14\ \mathrm{Chl}\ a-70.16\ \mathrm{Chl}\ b\right)/220 $$

Absorbance values were corrected for turbidity using absorbance at 750 nm. Total Chl, Chl a, Chl b, and carotenoid content were calculated on a leaf area basis (mg cm⁻² leaf area). All raw data are publically available (Donnelly et al. 2019).

2.3 Temperature data

Temperature (°C) was recorded and logged every 10 min by HOBO sensors at two locations within the study area. Data from the two sensors were averaged and plotted for the duration of the study.

2.4 Outliers, relationship equations, and statistical analysis

On visual inspection of the scatter plots of SPAD vs. total Chl, it was evident that some of the points may have been outliers. In order to statistically determine if any of the points were true outliers, histograms of the residuals of the regression for SPAD and Chl content for each species were constructed and any point that lay outside the normal distribution was considered an outlier and omitted from further analysis (ESMs 1 and 2).

IBM SPSS version 25.0 was used to fit linear, quadratic, and exponential regression equations to the data to explore the relationship between SPAD values and extracted leaf Chl content and to determine which model was the best fit for the data based on r² and root mean square error (RMSE). MATLAB version 2019a was used to (i) calculate and plot a 15-day moving window of the r² values for the relationship between SPAD and Chl content to determine the robustness of the relationship over time and (ii) compare changes in parameters over time between native and non-native groups using ANCOVA (analysis of covariance).

3 Results

3.1 Treatment of outliers

Prior to construction of the statistical models, we identified a small number of outliers within the data set. Based on the distribution plots of the residuals, six out of 1044 data points were identified as outliers, and omitted from further analysis as they lay outside the normal distribution of the data set (ESMs 1 and 2). The outliers were for two native (wild currant and dogwood) and one non-native (privet) species and all showed a higher than normal Chl content for a given SPAD value.

3.2 Temperature, leaf fluorescence, and chlorophyll content over time

Temperature (°C) dropped significantly (p < 0.0001) from an average of 15.9 °C (DOY 250–285) to 6.3 °C (DOY 285–312) (Fig. 1a) which was reflected in a marked decline in Φ_PSII values for both native and non-native shrubs (Fig. 1b). Interestingly, the temperature drop did not appear to have a detectable impact on the gradual decline in total leaf Chl content (mg cm⁻²) or SPAD values (Fig. 1c, e). However, the Chl:carotenoid ratio declined over time suggesting that leaf Chl degraded at a faster rate than carotenoids (Fig. 1d).

Fig. 1

Site temperature (a), quantum yield of photosystem II ΦPSII (b), total leaf chlorophyll content (c), total leaf chlorophyll:carotenoid ratio (d), and SPAD values (e) during the study period (DOY 250–312 September–November) for 5 native (N = 608 leaves) and 3 non-native (N = 430 leaves) shrub species. Note: scales on x axes differs slightly between (a) and the boxplots

3.3 Leaf chlorophyll content

Total leaf Chl content across the 8 shrub species ranged from 0.0018 to 0.0638 mg cm⁻² with native species exhibiting a slightly lower (p < 0.001) content than non-native species (0.0018 to 0.048 mg cm⁻² vs. 0.002 to 0.0638 mg cm⁻²) (Figs. 1c and 2). Similar trends were observed for Chl a and b but a decline in Chl a:Chl b ratio over time indicated faster decline in Chl a than Chl b (ESM 3). The Chl a:Chl b ratio decline was faster in native than non-native species (p < 0.002, ANCOVA). Total leaf Chl content also decreased during autumn senescence though at a similar rate for native and non-native shrub species (Fig. 1c; p = 0.7076 ANCOVA). Carotenoid content also declined as autumn progressed but at a faster rate in native than non-native species (ESM 4; p < 0.002 ANCOVA).

Fig. 2

Community (all shrub species (N = 1038) (a), native (N = 608) (b), and non-native (N = 430) (c)) linear (solid line), quadratic (dotted line), and exponential (dashed line) relationships between total leaf chlorophyll content and SPAD meter readings for 5 native and 3 non-native temperate deciduous shrub species. Table 2 presents corresponding conversion equations

3.4 Relationship between SPAD values and leaf chlorophyll content

The scatter plots showed strong positive correlations (p < 0.001) between SPAD values and extracted leaf Chl concentration for (i) the shrub community as a whole, (ii) native and non-native species groupings, (iii) individual native species, and (iv) individual non-native species (Figs. 2, 3, and 4). The corresponding linear, quadratic, and exponential calibration equations are presented in Tables 2 and 3. Overall, as indicated by relatively high r² and low RMSE values, the exponential model performed the poorest while the quadratic model performed the best but only marginally better than the linear model. Non-native shrub species showed both the strongest (honeysuckle) and weakest (buckthorn) species-specific correlations, whereas natives exhibited consistently strong relationships (Table 3).

Fig. 3

Scatter plots of relationships between total leaf chlorophyll content and SPAD meter readings for 5 native temperate deciduous shrub species. Table 3 presents corresponding linear, quadratic, and exponential calibration equations

Fig. 4

Scatter plots of relationships between total leaf chlorophyll content and SPAD meter readings for 3 non-native temperate deciduous shrub species. Table 3 presents corresponding linear, quadratic, and exponential calibration equations

Table 2 Calibration equations to convert SPAD (Minolta 502) values to total leaf chlorophyll content. Equations derived from 1038 (community), 608 (native), and 430 (non-native) randomly selected leaves from a range of native and non-native temperate deciduous shrub species respectively. All model fits were statistically significant (p < 0.0001)

Table 3 Calibration equations to convert SPAD (Minolta 502) values to total leaf chlorophyll content. Equations derived from randomly selected leaves from a range of native and non-native temperate deciduous shrub species. All model fits were statistically significant (p < 0.0001)

In general, the relationship between SPAD and extracted Chl was weaker at higher values showing higher variability, and this pattern was particularly evident in non-native shrub species (Figs. 2 and 4). Individual species such as buckthorn and privet showed a high degree of scatter throughout the range of Chl contents (Fig. 4).

In order to explore how the Chl-SPAD relationship changed as the autumn season progressed, we examined the relationship using a 15-day moving window of the r² values between SPAD and leaf Chl content (Fig. 5). The r² values indicated a weak relationship between SPAD and Chl in early-mid autumn which strengthened later in the season when both SPAD and Chl content were low. This pattern was more pronounced in non-native than native species (Fig. 5).

Fig. 5

15-day moving window of the r² values for the linear relationship between SPAD readings and total leaf chlorophyll content over the autumn season for 5 native (a) and 3 non-native (b) shrub species

4 Discussion

The majority of studies (Table 1) quantifying the relationship between leaf Chl content and SPAD values focused on a small number of typically agricultural crop species used to inform fertilizer management practices. The only studies reporting multi-species calibration equations were for tropical and subtropical tree species (Schaper and Chacko 1991; Marenco et al. 2009; Coste et al. 2010) while all the temperate deciduous trees for which equations were reported were either grown in a greenhouse (birch Betula papyrifera Richardson et al. (2002)) or under managed conditions (birch Betula pendula Uddling et al. (2007); sycamore Acer pseucoplatanus, oak Quercus robur, beech Fagus sylvatica Percival et al. (2008)). This study and review supports the comment by Coste et al. (2010) that SPAD measurements are underutilized in forest science and particularly in temperate deciduous forests. The current study is the first to report calibration equations (linear, quadratic, and exponential) for leaf Chl content and SPAD values for a range of native and non-native temperate deciduous shrub species growing under natural conditions. As reported in previous studies, albeit on different species, the SPAD meter proved to be a robust method of estimating leaf Chl content of deciduous shrubs in support of our first hypothesis.

A significant shift towards cooler temperatures during the study period (after DOY 285) coincided with a marked decline in Φ_PSII coupled with higher variability particularly for native shrubs. Interestingly, the temperature drop had no immediately detectable effect on the declining pigment (Chl or carotenoids) content or SPAD values. This suggests that the photosynthetic apparatus in photosystem II is more sensitive to abrupt decreases in temperature than the total pigment content of a leaf. A similar rapid drop in photosynthetic quantum yield, together with a gradual decline in Chl following a cold shock, was also reported for Aspen (Populus tremula) leaves (Keskitalo et al. 2005). In the current study, as the season progressed, decline in Chl:carotenoid ratios also indicated a slower decline for leaf carotenoid content than for Chl. Chl content degraded at a similar rate for both native and non-native species but carotenoid content decreased at a faster rate in native than non-native shrubs. This retention of carotenoids during loss of Chl has been reported for several species (Lichtenthaler 1987; Keskitalo et al. 2005) and may contribute to extended autumn productivity in non-native shrub species (Harrington et al. 1989; Knight et al. 2007). The range of Chl a:Chl b ratios (~ 1–4.5) in the shrubs was wider than in birch (Richardson et al. 2002), but similar to other deciduous shrubs (Parry et al. 2014). The faster decline in Chl a than Chl b (declining Chl a: Chl b ratio) suggests selective retention of chlorophyll antennae pigments during loss of photosystem chlorophylls as reported for Parthenocissus vine during autumn senescence (Lichtenthaler 1987) and a range of other species (Parry et al. 2014).

In this study we examined the SPAD-Chl relationship over a broad range of values, representing similar ranges reported for tree and other shrub species (e.g., Richardson et al. 2002; Netto et al. 2005; Uddling et al. 2007). Over this broad range, we found that a quadratic fit best described the relationship between total foliar Chl content and SPAD values for a range of temperate deciduous shrubs. Despite finding only marginal differences in the strength of the fit between linear and quadratic relationships, other authors reported that polynomial functions were better at describing the relationship between SPAD values and leaf Chl than linear functions for a range of species including wheat (Monje and Bugbee 1992; Finnan et al. 1997), paper birch (Richardson et al. 2002), coffee (Netto et al. 2005), citrus species (Jifon et al. 2005), spicebush (Hawkins et al. 2009), and Arabidopsis thaliana (Ling et al. 2011). Jifon et al. (2005) proposed that quadratic equations were better than linear models at describing the SPAD-Chl relationship in species with a wide range of chlorophyll concentrations. Furthermore, other authors reported exponential (soybean and maize, Markwell et al. 1995; birch, wheat potato Uddling et al. 2007), linear (grape, Steele et al. 2008), and homographic (neotropical trees, Coste et al. 2010) functions best described the SPAD-Chl relationship. In all the above studies, the reported coefficient of determination (r²) values were consistently above 0.90 compared with the lower r² values (0.32–0.82) we obtained suggesting that their models were able to explain a greater amount of the variability within the data. Likely reasons for the weaker relationships obtained for our plant species include (i) the longer period over which leaves were sampled, and (ii) the fact that our plants were growing in the wild as opposed to either in a greenhouse or a managed agricultural system.

SPAD values are typically recorded over the course of the senescence period as leaves age but leaves used for calibration determination are generally sampled over a much shorter time period typically on the same day or over a few days to weeks, suggesting that the full range of leaf ages and phenophases may not be attained. However, we took leaf samples for chlorophyll extraction throughout the entire autumn period (10 weeks). Therefore, this study captured a wider range of leaf ages and phenophases which better represents the seasonal change component reflective of a typical SPAD monitoring campaign, especially if applied to phenological or ecological studies. Furthermore, our sample plants were growing in a natural habitat compared with either in controlled environments (e.g., Richardson et al. 2002; Hawkins et al. 2009; Coste et al. 2010; Ling et al. 2011; Xiong et al. 2015) or under managed agricultural conditions (e.g., Finnan et al. 1997; Jifon et al. 2005; Netto et al. 2005; Steele et al. 2008) both of which would be expected to result in reduced intraspecific variation. In fact, Jifon et al. (2005) reported that differences in leaf thickness between field- and greenhouse-grown citrus leaves significantly influenced SPAD readings, with thicker field-grown leaves exhibiting higher values for a given Chl content. Therefore, changes in leaf thickness with age, such as leaf thinning as senescence progresses (Ito et al. 2006), may account for some of the variability observed in the relationship between SPAD and leaf Chl content reported here. Furthermore, the senescence period for agricultural crops typically occurs over a shorter period (several weeks) than for temperate deciduous shrubs and trees (several months). Consequently, even though the Chl-SPAD calibration equations reported here were reasonably robust over the study period, it may be useful to establish predictive equations for different phenophases as plants undergo dynamic changes in pigment content and photosynthetic processes over the growth period (Ito et al. 2006).

In agreement with previous reports for birch (Richardson et al. 2002; Uddling et al. 2007), citrus species (orange, lemon, grapefruit) (Jifon et al. 2005), macrolichens (Liu et al. 2019), and agricultural crops (wheat and potato) (Uddling et al. 2007), we found a weaker relationship between SPAD and Chl content at high concentrations of leaf Chl which supports our second hypothesis that the relationship between SPAD and Chl changes over time during autumn senescence and pigment breakdown in leaves. A possible explanation for this change may be due to an uneven distribution of Chl in leaves with high Chl content (Monje and Bugbee 1992; Richardson et al. 2002). This so-called sieve effect (Parry et al. 2014) was more pronounced in non-native than native shrub species. The greater variation in Chl content of non-native species may indicate a greater degree of plasticity for foliar Chl concentration (Harrington et al. 1989). Indeed, Martinez and Fridley (2018) suggested that non-natives may be more plastic than natives in relation to seasonal light changes which may result in more carbon uptake under the same environmental conditions but only if physiological function is maintained (Harrington et al. 1989; Caplan et al. 2018). Even though leaves may be green in late autumn, their photosynthetic capacity is uncertain as components of the photosynthetic apparatus may degrade before total Chl content shows signs of decline (Lichtenthaler 1987; Augspurger et al. 2005; Caplan et al. 2018). These factors may influence interpretation of satellite-derived end of season phenology. For example, commonly used satellite-based phenological metrics, such as NDVI/EVI time series, tend to exhibit an increase in late autumn which is likely related to the extended growing season of shrubs and understory plants after deciduous trees have lost their leaves.

As suggested by previous studies, inter-specific variation in leaf Chl content for a given SPAD value could be as a result of (i) the uneven structure and arrangement of Chl molecules throughout the leaf (Fukshansky et al. 1993; Uddling et al. 2007), (ii) the pattern and density of veins across the leaf surface (McClendon and Fukshansky 1990), (iii) the thickness of the leaf which may change with age (Jifon et al. 2005; Marenco et al. 2009), and (iv) a high rate of chloroplast movement leading to high leaf absorbance changes observed in shade leaves (Davis et al. 2011) all of which may lead to an uneven distribution of Chl within a leaf and may account for at least some of the species-specific differences observed in the current study. Even though we made every effort to avoid veins, especially prominent and mid-rib veins, it is highly likely that at least some leaf samples and SPAD values included a small portion of vein, thus contributing to this observed inter-specific variation.

5 Conclusions

The hand-held SPAD-502 meter produced reliable estimates of leaf Chl content for a range of, native and non-native, temperate deciduous shrubs. The progressively stronger relationship between SPAD values and Chl content during the course of autumn senescence as Chl degraded suggests the need for establishing calibration equations to cover a range of leaf Chl concentrations if absolute Chl content is required. Overall, coefficients of determination were higher for quadratic functions and for the native shrub community compared with the non-native community indicating a stronger relationship between SPAD and leaf Chl content. Given the interspecific variation observed in this study, we recommend deriving species-specific calibration equations for different shrub species to determine the relationship between SPAD and Chl content. Nonetheless, SPAD can provide a useful means of comparing relative Chl content if required. In future, we anticipate exploring the potential use of SPAD values to examine the, often observed, peak in remote sensing vegetation indices in late autumn when upper canopy trees are leafless.