Am. J. Bot. Join the BSA
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (35)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Wiemann, M. C.
Right arrow Articles by Wheeler, E. A.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Wiemann, M. C.
Right arrow Articles by Wheeler, E. A.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Wiemann, M. C.
Right arrow Articles by Wheeler, E. A.
(American Journal of Botany. 1998;85:1796-1802.)
© 1998 Botanical Society of America, Inc.


Estimation of temperature and precipitation from morphological characters of dicotyledonous leaves1

Michael C. Wiemanna, Steven R. Manchester5,a, David L. Dilchera, Luis Felipe Hinojosab and Elisabeth A. Wheelerc

a Florida Museum of Natural History, University of Florida, Gainesville, Florida 32611; b Laboratorio de Sistemática y Ecología Vegetal, Universidad de Chile, Santiago, Chile; and c Department of Wood and Paper Science, North Carolina State University, Raleigh, North Carolina 27695


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 RECOMMENDATIONS FOR FUTURE MODEL...
 REFERENCES
 
The utility of regression and correspondence models for deducing climate from leaf physiognomy was evaluated by the comparative application of different predictive models to the same three leaf assemblages. Mean annual temperature (MAT), mean annual precipitation (MAP), and growing season precipitation (GSP) were estimated from the morphological characteristics of samples of living leaves from two extant forests and an assemblage of fossil leaves. The extant forests are located near Gainesville, Florida, and in the Florida Keys; the fossils were collected from the Eocene Clarno Nut Beds, Oregon. Simple linear regression (SLR), multiple linear regression (MLR), and canonical correspondence analysis (CCA) were used to estimate temperature and precipitation. The SLR models used only the percentage of species having entire leaf margins as a predictor for MAT and leaf size as a predictor for MAP. The MLR models used from two to six leaf characters as predictors, and the CCA used 31 characters. In comparisons between actual and predicted values for the extant forests, errors in prediction of MAT were 0.6°–5.7°C, and errors in prediction of precipitation were 6–89 cm (=6–66%). At the Gainesville site, seven models underestimated MAT and only one overestimated it, whereas at the Keys site, all eight models overestimated MAT. Precipitation was overestimated by all four models at Gainesville, and by three of them at the Keys. The MAT estimates from the Clarno leaf assemblage ranged from 14.3° to 18.8°C, and the precipitation estimates from 227 to 363 cm for MAP and from 195 to 295 cm for GSP.

Key Words: canonical correspondence analysis • leaf morphology • leaf physiognomy • multiple linear regression • paleoclimate • simple linear regression


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 RECOMMENDATIONS FOR FUTURE MODEL...
 REFERENCES
 
The relationship between leaf morphology and climate was first documented by Sinnott and Bailey (1915), and since the publication of their findings leaf characters have been used to estimate temperature and precipitation. Bailey and Sinnott (1915) compared woody dicotyledonous species from four temperature regimes. They found that species that have entire-margined leaves predominate in frigid, subtropical, and tropical environments, whereas leaves with toothed margins are more frequent in temperate regions. Because of their small size, short life cycles, and ability to pass unfavorable periods underground or as seeds, herbs are more adaptable and variable than trees, and as a result the correlations between leaf margin and climate are weaker in herbs than in trees and large shrubs (Bailey and Sinnott, 1916).

Wolfe (1971) presented a comparison of mean annual temperature (MAT) and percentage of species with entire-margined leaves; for 19 modern floras, an increase from 10 to 86% of entire-margined species corresponded to an increase from 4° to 28°C in temperature. Wolfe's (1979) plot of MAT as a function of the percentage of entire-margined species in eastern Asian forests shows a remarkably strong correlation; similar plots presented by Greenwood (1992) for Australia show the same trend but with weaker correlations. Leaf length is also related to MAT; Carpenter et al. (1994) found a high correlation (r = 0.88) between MAT and the mean length of leaves collected from Australian leaf litter.

Studies in Australia (Webb, 1968), Ghana (Hall and Swaine, 1981), and worldwide (Dilcher, 1973) show that leaf size decreases with decreasing site rainfall. Givnish (1984) found that, among lowland tropical sites, leaves are narrower in drier environments. Wilf et al. (1998) found a strong (r = 0.760) relationship between the logarithm of mean annual precipitation (MAP) and the logarithm of average leaf area. Dolph and Dilcher (1980a, b) found that leaf size also varies with temperature, with larger leaves found at warmer sites.

The relationships among environment and leaf anatomy and morphology can be used to predict climate from an assemblage of extant or fossil leaves. This physiognomic method is based on the similarity in morphology and anatomy of the constituent species in areas with similar climates, even if the floras have very different taxonomic compositions. Bivariate analyses may be used to determine climate measures from single leaf morphological characters. The relationships may be illustrated by plots or may be expressed as simple linear regression (SLR) equations. SLR assumes that the relationship between a response and a single predictor is a straight line. This assumption is good for site MAT vs. percentage of woody dicotyledonous species having entire-margined leaves, as can be seen by Wolfe's (1979, fig. 8) plot, which has a coefficient of determination of r = 0.983 (Wing and Greenwood, 1993). Other relationships may be linear but much weaker; the plots of site MAP vs. percentage of dicotyledonous species with large leaves given by Wilf et al. (1998, fig. 3), which have coefficients of determination of r = 0.439 and r = 0.554, are examples.

Multivariate analyses, which simultaneously analyze multiple measurements, might be expected to improve the predictive ability of models. One commonly applied multivariate technique is multiple linear regression (MLR), in which a response is the function of more than one predictor variable. MLR models must be linear in the parameters, but may be linear (first-order) or curvilinear (higher than first-order) with respect to the predictors (Neter et al., 1996). Even though the bivariate relationships between the response and the individual predictors may not be linear, the relationship between the response and two or more (n) predictors taken together may be linear in (n + 1)-dimensional space. Thus an assumption of a linear bivariate relationship between response and individual predictors is not necessary to derive a valid multiple linear regression model.

Useful multivariate models may also be developed through ordination methods, which find weighted combinations that maximize the dispersion of the characters included in the ordination. Ordination methods include correspondence analysis (CA), which maximizes dispersion of the characters themselves, and canonical correspondence analysis (CCA), which integrates regression and ordination methods. As used in community ecology, CCA relates species abundances to explanatory variables, with no assumptions of linearity in the data. In CCA, ordination axes, which are linear combinations of the explanatory variables, are extracted from the data so as to maximize the dispersion among the species, with the constraint that the ordination axes must be uncorrelated with each other (i.e., they are orthogonal). Multiple regression can then be used to quantify the relationships between species abundances and scores on the ordination axes (ter Braak, 1986, 1987, 1988, 1995; Herman et al., 1996).

Wolfe (1993, 1995) developed ordination methods to estimate a number of temperature and precipitation variables based on the simultaneous use of multiple leaf characters. His original CLAMP (Climate-Leaf Analysis Multivariate Program) method used CA to predict climate based on 29 foliar characters related to margin, size, apex, base, and shape. Additional leaf characters, such as venation pattern (Dilcher, 1974), stomatal area (Davis and Taylor, 1980), or blade thickness (Roth, 1984), may also be related to environmental variables, but these are not included in Wolfe's database. Herman and Spicer (1996, 1997) used Wolfe's CLAMP database, with an additional two leaf size characters, to estimate paleotemperature and paleoprecipitation for four fossil assemblages using CCA. Kovach and Spicer (1996) also used Wolfe's data to calculate climate by CCA; they found that the method worked well for MAT in the range of 10°–20°C, but that temperatures above and below this range, as well as measures of annual, growing season, and dry season precipitation, could not be accurately estimated. They suggest that the inaccuracies at MAT < 10°C may be caused by an interaction of factors that affect leaf physiognomy only at low temperatures and that precipitation might not be in a form immediately available for plant use.

Wing and Greenwood (1993), Gregory (1994), and Gregory and McIntosh (1996) applied regression analysis to Wolfe's data and obtained regression equations that used from two to six characters as predictor variables. Jacobs and Deino (1996) used the equations derived by Wing and Greenwood (1993) to estimate modern climate at two sites in Uganda and paleoclimate at two sites in Kenya. Their estimates of MAT at the modern sites were within 1.1°C of the actual values, but MAP was overestimated by 96% at one site and 148% at the other.

The purpose of this study was to assess the ability of three types of models (SLR, MLR, and CCA) to predict climate and to estimate the paleoclimate of the Eocene Clarno Nut Beds. To do this we applied the relationships between leaf morphological characters and climate measures (temperature and precipitation) to assemblages of leaves from two extant forests in Florida and the extinct Clarno Nut Beds forest of Oregon. Leaf morphology was used to estimate MAT, MAP, and growing season precipitation (GSP) from models developed by other investigators, and to estimate MAT and GSP from regression models derived by us using Wolfe's data for 144 sites (CLAMP 3B database). We have used Wolfe's (1993) criterion for the calculation of GSP, which is the amount of precipitation that falls during those months in which the monthly mean temperature is at least 10°C. In the case of our two Florida sites, MAP is equal to GSP.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 RECOMMENDATIONS FOR FUTURE MODEL...
 REFERENCES
 
We computed correlation coefficients and their significance levels among leaf morphological characters, MAT, and GSP, using data from Wolfe's CLAMP 3B database. CLAMP 3B is a tabulation of data on 31 leaf characters and ten meteorological parameters from 144 sites located in temperate and tropical America, and Japan. The correlations were run using SAS Release 6.04. Although previous versions of CLAMP included MAP, the latest version does not. This is because winter precipitation at sites with less than a 12-mo growing season is not significant to the plants (J. Wolfe, personal communication).

In order to test the ability of various models to predict temperature and precipitation, two validation sites with 30-year climate records were chosen to provide leaf samples for analysis. One site, "Dilcher's Florida Woods," is a 25-ha woodland located 3 km east of Gainesville, Florida (29.6° N latitude, 82.2° W longitude, 20 m elevation). The MAT, MAP, and GSP of the site are, respectively, 21.0°C, 134 cm, and 134 cm, measured over the period 1951–1980 (NOAA, 1985) at Gainesville Airport. The species from the other site, located in the southern Florida Keys (24.7° N latitude, 81.4° W longitude, <2 m elevation), were collected on No Name Key, Big Pine Key, and Little Torch Key. The MAT, MAP, and GSP of the closest weather station, on Key West (24.5° N, 81.7° W), are, respectively, 25.4°C, 100 cm, and 100 cm, also for the period 1951–1980 (NOAA, 1985). Samples of living leaves, representing both sun and shade conditions, were collected from trees on the sites and preserved in a herbarium press. When a species was represented by more than one tree, the leaves of that species were combined. Only leaves of woody dicots were included in the sample; herbaceous species were excluded because their leaves do not fall off and contribute to leaf litter (Gregory, 1994) and because the correlations between their leaf characters and climate are weak (Bailey and Sinnott, 1916). Whether or not a species contributes to the leaf litter is important in paleoclimate determination by foliar physiognomy because litter is most likely the source of fossil leaves. The leaf assemblages, which contained 24 species for the Gainesville site and 30 species for the Florida Keys site, were used to compare the results of various methods of climate estimation to actual climate.

The climate models used to compare predicted with actual MAT, MAP, and GSP at the validation sites were also used to determine paleoclimate at a fossil site. This site, the Nut Beds of the middle Eocene Clarno Formation, is known for its well-preserved silicified fruits, seeds, leaves, and wood (Scott, 1954; Scott and Wheeler, 1982; Manchester, 1994) and is located 3 km east of Clarno, Oregon, at 44.9° N latitude, 120.4° W longitude and ~550 m elevation. An assemblage collected by S. R. Manchester of >500 fossil leaf specimens, representing 74 distinct leaf types, was used to estimate paleoclimate. The leaf assemblage from the Clarno Nut Beds has not been fully described (Manchester, 1981), but is available for study at the Florida Museum of Natural History. The Nut Beds site is ~44 million years old, based on fission-track (43.6 and 43.7 ± 10%; Vance, 1988) and Ar/Ar (43.76 ± 0.29; Turrin, in Manchester, 1994, p. 13) dating of pumice fragments collected from the fossil-bearing tuff.

The leaves in the living and fossil assemblages were examined for the presence or absence of each of 31 morphological characters, described by Wolfe (1993) and Herman et al. (1996), and listed in Table 1. The character "no teeth" is the same as "margin entire" of Bailey and Sinnott (1916) except that spinose leaves are included in the "no teeth" category. The percentage occurrences of the characters, shown in Table 1 for each of the assemblages, were used to calculate climate using SLR, MLR, and CCA models developed by other researchers and using MLR models derived by us using Wolfe's CLAMP 3B database. The SLR models used absence of teeth on leaf margins as the predictor variable for MAT and the presence of large leaves as the predictor variable for MAP. The models were derived from plots of MAT vs. percentage entire-margined species published by Wolfe (1979, fig. 8) for species of eastern Asia, Greenwood (1992, fig. 12) for species of Australia, and Wilf (1997, fig. 1) for species of temperate and tropical America, and from the plot of MAP vs. percentage large-leaved species published by Wilf et al. (1998, fig. 3) based on Wolfe's CLAMP database. Converting the plots into linear equations gives


View this table:
[in this window]
[in a new window]
 
Table 1. Leaf physiognomic characters, abbreviations used in MLR equations (in parentheses), and correlation coefficients of each character with MAT and GSP. Correlation coefficients were computed from Wolfe's 144 sites (CLAMP 3B database) and are significant at the 0.01 level unless marked ns. Percentage of leaf types given for assemblages from the Clarno Nut Beds (74 fossil leaf types), Gainesville (24 extant species), and the southern Florida Keys (30 extant species).

 
MAT = 1.14 + 0.306 x (% entire)   (Wolfe, 1979; Wing andGreenwood, 1993)

MAT = 4.4 + 0.22 x (% entire)   (Greenwood, 1992)

MAT = 2.24 + 0.286 x (% entire)   (Wilf, 1997)

MAP = 47.5 + 6.18 x (% large leaves)    (Wilf et al., 1998)

In these equations, (% entire) is the percentage of leaves in an assemblage that have entire (untoothed) margins, and (% large leaves) is the percentage of leaves in an assemblage that are size mesophyll I or larger (i.e., area >=33 cm)

We also calculated climate from the MLR equations published by Wing and Greenwood (1993), Gregory (1994), and Gregory and McIntosh (1996), which are based on Wolfe's data or subsets of Wolfe's data. They used angular transformations in their regression analyses of leaf characters to estimate climate variables.

In order to reduce multicollinearity, Wing and Greenwood (1993) eliminated some of the potential regressors, selecting from the characters in Table 1 one character from each leaf form category that had low correlation with the other characters and reasonably high correlation with the climate variable of interest. They used only 74 sites, eliminating those with very low (less than -2.0°C) cold-month mean temperatures. Their multiple regression models for MAT and MAP are as follows.

MAT = 2.536 + 17.372(none) + 2.896(emarg) - 8.592(<1:1)   r = 0.863

MAP = 11.489 + 167.948(atten) + 377.735(mesoII)   r = 0.497

Gregory (1994) used 29 of the 31 characters listed in Table 1; she deleted the characters "teeth regular" and "teeth close" because of their high collinearity with "no teeth." She used 84 of the 86 sites that were included in Wolfe's database at that time, eliminating the two coldest sites because of their anomalous residuals. Her stepwise regression included quadratic functions of the character frequencies. Her leaf size leptophyll III corresponds to leptophyll II as defined by Herman et al. (1996) (K. Gregory, Lamont-Doherty Earth Observatory of Columbia University, personal communication). We use leptophyll II in this paper, and have made the substitution in her equation, given below.

MAT = 15.32 + 10.34(none) + 5.48(emarg) - 15.32(<1:1) - 15.29(leptII) - 5.79(micrII)   r = 0.938

Gregory and McIntosh (1996) used the same 29 characters to derive models that did not use quadratic functions. They derived two MLR models, one based on all 106 sites included in Wolfe's CLAMP database at that time and another based on only the sites with a cold month mean temperature of at least -2°C. We have used the model based on the larger database for two reasons: we have already chosen one model derived from the reduced database, the Wing and Greenwood (1993) model discussed above, and we wanted to examine models that were as inclusive as possible so that if a flora comes from a cold site it can be legitimately evaluated. The Gregory and McIntosh (1996) model for MAT, based on all 106 sites, follows.

MAT = -11.262 + 10.282(lobed) + 23.258(none) - 12.211(leptII) + 7.022(Bacut) - 16.099(<1:1) + 11.484(1–2:1)   r = 0.893

For the models derived by Wing and Greenwood (1993) and Gregory (1994), the quantities indicated in parentheses represent the arcsines (in radians) of the square roots of the proportions of leaves in the assemblage having the indicated characters, which are summarized in Table 1. For the model derived by Gregory and McIntosh (1996), 0.005 was added to each proportion before taking the arcsine of the square root.

In addition to the models developed by other researchers, we used Wolfe's CLAMP 3B data to develop our own multiple linear regression models using the STEPWISE function in SAS Release 6.04. The untransformed data from all 144 sites and 31 morphological characters were used to derive MLR models for MAT and GSP at significance levels of {alpha} = 0.001 for the F statistics.

The CCA model used the CANOCO Version 3.12 program developed by C. J. F. ter Braak and available from Microcomputer Power; we used Wolfe's data from 144 sites (CLAMP 3B) to establish the ordinations. The method for doing this using leaf physiognomic data is described by Herman et al. (1996). The leaf character percentages for the Gainesville, Florida Keys, and Clarno assemblages were entered as passive samples (i.e., without climate data). The Axis 1 and 2 scores from the CANOCO output were used to calculate MAT, and the Axis 2 and 5 scores were used to calculate GSP, from power functions furnished by J. Wolfe.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 RECOMMENDATIONS FOR FUTURE MODEL...
 REFERENCES
 
Examination of the correlation coefficients between the climate variables and the leaf characters (Table 1) shows that leaf morphology has a stronger relationship with temperature than with precipitation. Whereas nine of the characters have correlation coefficients >0.7 with MAT, none of the correlation coefficients with GSP are this high. The five characters most highly correlated with MAT all relate to leaf margin. In contrast, leaf margin characters correlate poorly with GSP. The leaf characters most strongly correlated with GSP include L:W ratio, leaf shape, presence of an attenuate apex, and leaf size.

The STEPWISE procedure, using the untransformed data from CLAMP 3B, gave the following MLR models.

MAT = 9.865 + 0.207(none) - 0.058(Brnd) - 0.202(<1:1)   r = 0.898

GSP = 31.6 - 3.393(leptII) + 2.400(atten) - 2.671(Bcord) + 2.360(2–3:1) + 3.122(3–4:1)   r = 0.796

The quantities indicated in parentheses represent the percentages of leaves in an assemblage having the characters summarized in Table 1.

The character most strongly correlated with MAT, leaf margin with no teeth, is less abundant in the fossil assemblage than in the extant assemblages (Table 1), suggesting that Eocene Clarno was cooler than present-day Florida. Using this character alone, estimation of MAT by SLR gave values ranging from 14.3° to 15.1°C for Clarno, compared with values ranging from 17.2° to 18.9°C for Gainesville (actual = 21.0°C) and from 26.0° to 31.1°C for the Keys (actual = 25.4°C) (Table 2).


View this table:
[in this window]
[in a new window]
 
Table 2. MAT, MAP, and GSP for Clarno Nut Beds, Gainesville, and the southern Florida Keys, calculated from simple linear regression (SLR), multiple linear regression (MLR), and canonical correspondence analysis (CCA) models.

 
The leaves in the fossil assemblage are larger than those in the extant assemblages: 51% of them are in the mesophyll categories, whereas only 28% of the species in the Gainesville assemblage and 10% in the Keys assemblage have leaves this large (Table 1). This suggests that Clarno was probably wetter than present-day Florida. Another character that is significantly correlated with precipitation, leaf apex attenuate (r = 0.37), is more abundant in the fossil assemblage, also indicating that Clarno was wetter than Florida. Estimation of MAP by the single character "large leaves," using the model of Wilf et al. (1998), gave a value of 363 cm for Clarno, 221 cm for Gainesville (actual = 134 cm), and 109 cm for the Keys (actual = 100 cm) (Table 2).

Use of several characters and MLR gave MAT estimates of 14.7°–18.8°C for Clarno, 15.4°–21.8°C for Gainesville, and 27.7°–30.6°C for the Keys. The estimates obtained using all 31 characters and CCA are 17.0°C for Clarno, 18.9°C for Gainesville, and 28.7°C for the Keys (Table 2).

For Clarno, MAP was estimated as 227 cm using Wing and Greenwood's (1993) MLR equation, whereas GSP was estimated as 195 cm using our MLR equation and 295 cm using CCA. For Gainesville, the corresponding values are 159, 172, and 223 cm, respectively, and for the Keys they are 94, 112, and 146 cm, respectively (Table 2).

The SLR, CCA, and three of the four MLR models underestimated MAT at Gainesville—the SLR models by 2.1°–3.8°C, the CCA model by 2.1°C, and the MLR models by 0.6°–4.9°C; the other MLR model overestimated it by 0.8°C. The opposite trend was found for the Keys, where MAT was overestimated by every model; the SLR models by 0.6°–5.7°C, the CCA model by 3.3°C, and the MLR models by 2.3°–5.2°C. All four precipitation models overestimated MAP in Gainesville (by 19–66%). Three of them overestimated it in the Keys (by 9–46%), and one of them underestimated it, but by only 6%.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 RECOMMENDATIONS FOR FUTURE MODEL...
 REFERENCES
 
Because MAT was underestimated by seven of the eight models at Gainesville and overestimated by all eight models in the Keys, the errors in estimation are probably due, at least in part, to characteristics of the floras. Examination of the estimates of precipitation show that it was overestimated by all four models at Gainesville and by three of the four models in the Keys. Only one model underestimated precipitation, and that error was only 6%. Thus, the errors in estimation of precipitation may be due to characteristics of the floras, or it may be that the models are biased. In their investigation of the usefulness of CCA for paleoclimate determination, Kovach and Spicer (1996) found that precipitation variables were skewed, and they used logarithmic transformations to compensate. Wilf et al. (1998) showed an approximately linear relationship between the logarithm of leaf area and the logarithm of MAP. Therefore, more valid models might be obtainable by using logarithmic functions for precipitation models.

The estimates of MAT at Clarno ranged from 14.3° to 18.8°C, a difference of 4.5°C among models. The models were derived from data based on measurements of leaves sampled directly from living plants (Wolfe, 1993), which may not accurately represent leaf litter or, consequently, fossil leaves. The composition of assemblages of fossil leaves is also affected by the processes of transport, degradation, and burial, and these may introduce additional biases. Therefore, the use of these models to deduce paleoclimate may be in error at least in part due to differences in the characteristics of fossil and living leaves. Greenwood (1992) found that litter is dominated by leaves from the canopy trees, and that leaves were smaller in the litter than in the canopy. Roth and Dilcher (1978) found that leaves deposited in an Indiana lake bottom were smaller than leaves from the surrounding forest and that fewer entire-margined species were recovered from the lake. Untransported, undecayed leaf litter reflects the forest tree composition in temperate forests (Burnham, 1993; Burnham, Wing, and Parker, 1992), but not in subtropical and tropical forests, probably due to the high species richness of these forests (Burnham 1989, 1993, 1994). Burnham (1997) used the Wing and Greenwood (1993) MLR model to calculate MAT of Santa Rosa National Park in Costa Rica based on leaf litter collections. She found that it underestimated MAT by more than 7°C; however, the same model also gave an estimate of MAT for Clarno that was 2.1°–4.1°C cooler than those obtained by the other multivariate models (Table 2), indicating that the low estimates may be due at least in part to bias in the model. Greenwood and Wing (1995) state that calculation of paleoclimate using their MLR approach yields MAT estimates 2°–3°C cooler than estimates based on leaf margins or nearest living relatives because it is sensitive to the low proportions of species in Eocene floras with emarginate apices, and because of deciduousness induced by light rather than temperature seasonality. In deciduous forests, leaves with L:W = 1:1 predominate, whereas in evergreen forest the predominant ratio is L:W > 3:1 (Greenwood and Basinger, 1994). However, for the Clarno assemblage, the result from using Wing and Greenwood's (1993) model (MAT = 14.7°C), which was based on leaf margin, apex, and length-to-width ratio, was not much different from the SLR results (MAT = 14.3°–15.1°C) (Table 2), which are based only on leaf margin; furthermore, Gregory's (1994) model, which also included emarginate apex and length-to-width ratio, gave the highest estimate of any model (MAT = 18.8°C).

It is not clear which, if any, of the models shown in Table 2 give more accurate climate determinations. The only model to estimate MAT to within 3°C at both validation sites was the MLR model that we derived based on untransformed leaf physiognomy data. Our precipitation model also gave reasonably close estimates (within 30%) of GSP at the two sites. However, because both sites were in Florida, possible biases in the models might not have been evident; validation sites that are considerably colder and wetter need to be sampled.

At the two validation sites, SLR seemed to be as good at predicting MAT as some of the multivariate models, although the errors exceeded 3°C in at least one site for each SLR model. Wilf (1997) calculated MAT at tropical and temperate sites using SLR and MLR and also found that SLR was at least as precise as MLR; he maintains that using additional characters beyond leaf margin contributes little information and may increase error because scoring procedures for many of the characters are subjective, making it difficult to score them reliably. Although we concur with his contention that difficult-to-score characters may increase error and have experienced the same scoring problems that he described in determining the presence or absence of an acute base (Wilf, 1997, p. 379), we are inclined to agree with Wolfe (1995) that multiple characters might provide more information than a single character alone, especially considering the interactions between temperature and precipitation. Furthermore, use of additional characters might act as a "buffer" to reduce the effect of an anomalously high or low frequency of a single character in an assemblage of leaves. A comparison of the number of characters used in the models with the accuracy of prediction of MAT is illustrative of this point. As can be seen in Table 3, the largest error (5.7°C) was found in a model that used only one character. In contrast, the model that used all 31 characters had a maximum error of only 3.3°C, even though many of the characters may have poor correlation with MAT.


View this table:
[in this window]
[in a new window]
 
Table 3. Number of leaf characters in MAT model and range in errors in predicted MAT at the two validation sites (absolute values) for each model.

 
Angular transformations and use of quadratic functions did not improve the ability of a model to predict climate at the Florida sites. The best overall models were our MLR, which predicted MAT within 2.8°C, and CCA, which predicted MAT within 3.3°C (Table 3); neither of these used a transformation or power function of the raw data. The unjustified use of power functions may be especially serious. If the relationships between MAT and the frequencies of occurrence of leaf characters are in fact approximately linear, then relatively small biases in the floras can result in relatively large errors in estimation due to the effects of exponentiation. Fitting a curve to data can result in a high coefficient of determination even if the true fit is a straight line, but the resulting model will not represent the phenomenon well and its application to new data may give large errors. An example is Gregory's (1994) model, which uses quadratic functions and has a coefficient of determination of 0.938, but which underestimated MAT by 4.9°C at the Gainesville site and overestimated it by 4.6°C at the Keys site.

Including too many independent variables in a model can have the same effect as fitting a curve when the actual phenomenon is linear. Inclusion of variables that do not belong decreases a model's descriptive abilities and increases the problem of roundoff errors; on the other hand, omitting important variables can lead to serious bias (Neter et al., 1996). Keeping in mind the difficulty of accurately and consistently evaluating leaf base shape (Wilf, 1997), the contribution of the character "base round" to our MLR model for MAT should be examined to rule out the possibility that its inclusion is just an artifact of the data and the STEPWISE procedure, particularly because its correlation with MAT is non-significant. The other two characters in this model, "leaf with no teeth" and "length less than width," in contrast, can be rigorously defined and easy to evaluate, and they have high correlations with MAT. Our MLR model for GSP includes two characters that measure leaf shape: "length 2–3 times greater than width" and "length 3–4 times greater than width." Their highly significant correlation coefficient is +0.49, and it would seem that one or the other of them might profitably be eliminated from the analysis.


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 RECOMMENDATIONS FOR FUTURE MODEL...
 REFERENCES
 
1) Some of the leaf morphological characters are difficult to score reliably, and more accurate models might be developed if these characters are not used.

2) Some characters are unimportant in climate prediction, and others duplicate information. These might profitably be deleted to simplify scoring and to reduce "noise."

3) Multiple character models as they are presently formulated may not be any better than single character models.

4) Angular transformations of character data do not seem to improve models.

5) The precipitation models gave biased estimates; a logarithmic transformation of precipitation data may be appropriate to derive improved models.

6) Leaf size is one of the most important predictors of precipitation. It may, however, be worth exploring measurements of size that are not based on categories, since the ones presently defined for the CLAMP databases may not be optimum for climate determination, and it is often difficult to decide which category a leaf falls into, especially if it has a shape that does not correspond to one of the template choices. On the other hand, an advantage of using size categories determined by use of a template is ease of use, especially for the measurement of fossils.


    RECOMMENDATIONS FOR FUTURE MODEL DERIVATION AND TESTING
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 RECOMMENDATIONS FOR FUTURE MODEL...
 REFERENCES
 
1) Have several individuals score the same assemblages of leaves to determine which leaf characters can be evaluated consistently, and eliminate those that cannot.

2) Compare the correlation coefficients among leaf characters to determine which ones may be superfluous due to collinearity.

3) Verify whether data manipulation, such as the use of angular or logarithmic transformations, yields models with increased predictive capacity.

4) Use STEPWISE regression procedures at various F levels to compile a list of leaf characters to be used as potential predictor variables.

5) Derive regression and correspondence analysis models based on a more limited number of well-defined, easy-to-evaluate leaf characters.

6) Test models using validation sites from a wide range of climates to determine which models give the best predictions and to evaluate the precision of leaf physiognomy as a climate estimator.


    FOOTNOTES
 
1 This research was supported in part by D. L. Dilcher Research Funds from the University of Florida. Terry Lott and Anthony Davanzo assisted with the collection and identification of the leaves from Gainesville, and Daniel Uriarte, Edwin DuPont, and Luis Contreras helped to score them. Garry Eckerson and Dianna Stevenson helped with the collection and identification of the leaves from the Florida Keys. The authors thank Jack Wolfe for making available his extensive climate and leaf physiognomic data. Robyn J. Burnham, Kathryn M. Gregory, and Peter Wilf offered valuable suggestions and criticisms that greatly improved the manuscript. This is Report Number 494 of the series University of Florida Contributions to Paleobiology . Back

5 Author for correspondence. Back


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 RECOMMENDATIONS FOR FUTURE MODEL...
 REFERENCES
 
Bailey, I. W., and E. W. Sinnott.1915A botanical index of Cretaceous and Tertiary climates. Science 41:831 –834.

———, and ———.1916The climatic distribution of certain types of angiosperm leaves. American Journal of Botany 3:24 –39.

Burnham, R. J.1989Relationships between standing vegetation and leaf litter in a paratropical forest: implications for paleobotany. Review of Palaeobotany and Palynology58 :5 –32.

———.1993Reconstructing richness in the plant fossil record. Palaios8 :376 –384.[CrossRef][ISI]

———.1994Patterns in tropical leaf litter and implications for angiosperm paleobotany. Review of Palaeobotany and Palynology81 :99 –113.[CrossRef][ISI]

———.1997Stand characteristics and leaf litter composition of a dry forest hectare in Santa Rosa National Park, Costa Rica. Biotropica29 :384 –395.[CrossRef][ISI]

———, S. L. Wing, and G. G. Parker.1992 The reflection of deciduous forest communities in leaf litter: implications for autochthonous litter assemblages from the fossil record. Paleobiology18 :30 –49.

Carpenter, R. J., R. S. Hill, and G. J. Jordan.1994Cenozoic vegetation in Tasmania: macrofossil evidence. In R. S. Hill [ed.], History of the Australian vegetation: Cretaceous to Recent, 276–298. Cambridge University Press, Cambridge.

Davis, J. M., and S. E. Taylor.1980Leaf physiognomy and climate: a multivariate analysis. Quaternary Research14 :337 –348.

Dilcher, D. L.1973A paleoclimatic interpretation of the Eocene floras of southeastern North America. In A. Graham [ed.], Vegetation and vegetational history of northern Latin America, 39–59. Elsevier Scientific Publishing Company, Amsterdam.[CrossRef][ISI]

———.1974Approaches to the identification of angiosperm leaf remains. Botanical Review 40:1 –157.[Abstract]

Dolph, G. E., and D. L. Dilcher.1980aVariation in leaf size with respect to climate in Costa Rica. Biotropica12 :91 –99.[Abstract/Free Full Text]

———, and ———.1980bVariation in leaf size with respect to climate in the tropics of the Western Hemisphere. Bulletin of the Torrey Botanical Club107 :154 –162.[CrossRef]

Givnish, T. J.1984Leaf and canopy adaptations in tropical forests. In E. Medina, H. A. Mooney, and C. Vázquez-Yánes [eds.], Physiological ecology of plants of the wet tropics, 51–84. Dr. W. Junk, The Hague.

Greenwood, D. R.1992Taphonomic constraints on foliar physiognomic interpretations of Late Cretaceous and Tertiary palaeoclimates. Review of Palaeobotany and Palynology71 :149 –190.

———, and J. F. Basinger.1994The paleoecology of high-latitude Eocene swamp forests from Axel Heiberg Island, Canadian High Arctic. Review of Palaeobotany and Palynology81 :83 –97.

———, and S. L. Wing.1995Eocene continental climates and latitudinal temperature gradients. Geology23 :1044 –1048.

Gregory, K. M.1994Palaeoclimate and palaeoelevation of the 35 ma Florissant flora, Front Range, Colorado. Palaeoclimates 1:23 –57.

———, and W. C. McIntosh.1996Paleoclimate and paleoelevation of the Oligocene Pitch-Pinnacle flora, Sawatch Range, Colorado. Geological Society of America Bulletin108 :545 –561.

Hall, J. B., and M. D. Swaine.1981Distribution and ecology of vascular plants in a tropical rain forest. Dr. W. Junk, The Hague.

Herman, A. B., and R. A. Spicer.1996Palaeobotanical evidence for a warm Cretaceous Arctic Ocean. Nature 380:330 –333.

———, and ———.1997New quantitative palaeoclimate data for the Late Cretaceous Arctic: evidence for a warm polar ocean. Palaeogeography, Palaeoclimatology, Palaeoecology 128:227 –251.

———, ———, J. A. Wolfe, and M. E. Collinson.1996A guide to obtaining a CLAMP analysis. Unpublished manuscript.

Jacobs, B. F., and A. L. Deino.1996Test of climate-leaf physiognomy regression models, their application to two Miocene floras from Kenya, and Ar/Ar dating of the Late Miocene Kapturo site. Palaeogeography, Palaeoclimatology, Palaeoecology 123:259 –271.

Kovach, W. L., and R. A. Spicer.1996Canonical correspondence analysis of leaf physiognomy: a contribution to the development of a new palaeoclimatological tool. Palaeoclimates2 :125 –138.

Manchester, S. R.1981Fossil plants of the Eocene Clarno Nut Beds. Oregon Geology 43:75 –81.

———.1994Fruits and seeds of the Middle Eocene Nut Beds Flora, Clarno Formation, Oregon, Palaeontographica Americana58 :1 –205.

Neter, J., M. H. Kutner, C. J. Nachtsheim, and W. Wasserman.1996Applied linear statistical models, 4th ed. Richard D. Irwin, Inc., Chicago.

NOAA.1985Climates of the states, 3rd ed. Gale Research, Detroit, MI.

Roth, I.1984Stratification of tropical forests as seen in leaf structure. Dr. W. Junk, The Hague.

Roth, J. L., and D. L. Dilcher.1978Some considerations in leaf size and leaf margin analysis of fossil leaves. Courier Forschungsinstitut Senckenberg 30:165 –171.

Scott, R. A.1954Fossil fruits and seeds from the Eocene Clarno Formation of Oregon. Palaeontographica 96:66 –97.

———, and E. A. Wheeler.1982Fossil woods from the Eocene Clarno Formation of Oregon. IAWA Bulletin n.s. 3: 135–154.

Sinnott, E. W., and I. W. Bailey.1915Foliar evidence as to the ancestry and early climatic environment of the angiosperms. American Journal of Botany 2:1 –22.

ter Braak, C. J. F.1986Canonical correspondence analysis: a new eigenvector technique for multivariate direct gradient analysis. Ecology 67:1167 –1179.

———.1987The analysis of vegetation-environment relationships by canonical correspondence analysis. Vegetatio69 :69 –77.

———.1988CANOCO—a FORTRAN program for canonical community ordination. Microcomputer Power, Ithaca, NY.

———.1995Ordination. In R. H. G. Jongman, C. J. F. ter Braak, and O. F. R. van Tongeren [eds.], Data analysis in community and landscape ecology, 91–173. Cambridge University Press, Cambridge.

Vance, J. A.1988New fission track and K-Ar ages from the Clarno formation, Challis age volcanic rocks in north central Oregon. Geological Society of America, Rocky Mountain Section, Abstracts with Programs20 : 473 (Abstract).

Webb, L. J.1968Environmental relationships of the structural types of Australian rain forest vegetation. Ecology 49:296 –311.

Wilf, P.1997When are leaves good thermometers? a new case for leaf margin analysis. Paleobiology 23:373 –390.

———, S. L. Wing, D. R. Greenwood, and C. L. Greenwood.1998Using fossil leaves as paleoprecipitation indicators: an Eocene example. Geology 26:203 –206.

Wing, S. L., and D. R. Greenwood.1993Fossils and fossil climate: the case for equable continental interiors in the Eocene. Philosophical Transactions of the Royal Society of London, Series B341 :243 –252.

Wolfe, J. A.1971Tertiary climatic fluctuations and methods of analysis of Tertiary floras. Palaeogeography, Palaeoclimatology, Palaeoecology 9:27 –57.

———.1979Temperature parameters of humid to mesic forests of eastern Asia and relation to forests of other regions of the northern hemisphere and Australasia. U.S. Geological Survey Professional Paper 1106, Washington, DC.

———.1993A method of obtaining climatic parameters from leaf assemblages. U.S. Geological Survey Bulletin 2040, Washington, DC.

———.1995Paleoclimatic estimates from tertiary leaf assemblages. Annual Review of Earth and Planetary Science 23:119 –142.




This article has been cited by other articles:


Home page
PALAIOSHome page
J. D. KRIEGER, R. P. GURALNICK, and D. M. SMITH
GENERATING EMPIRICALLY DETERMINED, CONTINUOUS MEASURES OF LEAF SHAPE FOR PALEOCLIMATE RECONSTRUCTION
Palaios, March 1, 2007; 22(2): 212 - 219.
[Abstract] [Full Text] [PDF]


Home page
PaleobiologyHome page
N. S. Gupta, M. E. Collinson, D. E. G. Briggs, R. P. Evershed, and R. D. Pancost
Reinvestigation of the occurrence of cutan in plants: implications for the leaf fossil record
Paleobiology, September 1, 2006; 32(3): 432 - 449.
[Abstract] [Full Text] [PDF]


Home page
Geological Society, London, Special PublicationsHome page
P. A. Hayes, J. E. Francis, D. J. Cantrill, and J. A. Crame
Palaeoclimate analysis of Late Cretaceous angiosperm leaf floras, James Ross Island, Antarctica
Geological Society, London, Special Publications, January 1, 2006; 258(1): 49 - 62.
[Abstract] [PDF]


Home page
PALAIOSHome page
R. J. BURNHAM, B. ELLIS, and K. R. JOHNSON
Modern Tropical Forest Taphonomy: Does High Biodiversity Affect Paleoclimatic Interpretations?
Palaios, October 1, 2005; 20(5): 439 - 451.
[Abstract] [Full Text] [PDF]


Home page
PALAIOSHome page
D. R. GREENWOOD
Leaf Margin Analysis: Taphonomic Constraints
Palaios, October 1, 2005; 20(5): 498 - 505.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Bot.Home page
D. L. Royer, P. Wilf, D. A. Janesko, E. A. Kowalski, and D. L. Dilcher
Correlations of climate and plant ecology to leaf size and shape: potential proxies for the fossil record
Am. J. Botany, July 1, 2005; 92(7): 1141 - 1151.
[Abstract] [Full Text] [PDF]


Home page
PALAIOSHome page
Paleotemperature Estimation Using Leaf-Margin Analysis: Is Australia Different?
Palaios, April 1, 2004; 19(2): 129 - 142.



Home page
Proc. Natl. Acad. Sci. USAHome page
P. Wilf, K. R. Johnson, and B. T. Huber
Correlated terrestrial and marine evidence for global climate changes before mass extinction at the Cretaceous-Paleogene boundary
PNAS, January 21, 2003; 100(2): 599 - 604.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
E. A. Kowalski and D. L. Dilcher
Warmer paleotemperatures for terrestrial ecosystems
PNAS, January 7, 2003; 100(1): 167 - 170.
[Abstract] [Full Text] [PDF]


Home page
PaleobiologyHome page
Estimation of low-latitude paleoclimates using fossil angiosperm leaves: examples from the Miocene Tugen Hills, Kenya
Paleobiology, September 1, 2002; 28(3): 399 - 421.



Home page
Am. J. Bot.Home page
R. J. Burnham, N. C. A. Pitman, K. R. Johnson, and P. Wilf
Habitat-related error in estimating temperatures from leaf margins in a humid tropical forest
Am. J. Botany, June 1, 2001; 88(6): 1096 - 1102.
[Abstract] [Full Text] [PDF]


Home page
PaleobiologyHome page
Relationships between leaf morphology and climate, Bolivia: implications for estimating paleoclimate from fossil floras
Paleobiology, December 1, 2000; 26(4): 668 - 688.



This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (35)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Wiemann, M. C.
Right arrow Articles by Wheeler, E. A.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Wiemann, M. C.
Right arrow Articles by Wheeler, E. A.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Wiemann, M. C.
Right arrow Articles by Wheeler, E. A.


HOME HELP <