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 (20)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kouwenberg, L. L. R.
Right arrow Articles by Visscher, H.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Kouwenberg, L. L. R.
Right arrow Articles by Visscher, H.
Agricola
Right arrow Articles by Kouwenberg, L. L. R.
Right arrow Articles by Visscher, H.
(American Journal of Botany. 2003;90:610-619.)
© 2003 Botanical Society of America, Inc.


Ecology

Stomatal frequency adjustment of four conifer species to historical changes in atmospheric CO21

Lenny L. R. Kouwenberg2,6, Jennifer C. McElwain3, Wolfram M. Kürschner2, Friederike Wagner2, David J. Beerling4, Francis E. Mayle5 and Henk Visscher2

2Laboratory of Palaeobotany and Palynology, Utrecht University, 3584 CD Utrecht, Netherlands; 3Department of Geology, Field Museum of Natural History, 1400 S. Lake Shore Drive, Chicago, Illinois 60605 USA; 4Department of Animal and Plant Sciences, University of Sheffield, Sheffield S10 2TN UK; 5Department of Geography, University of Leicester, Leicester LE1 7RH UK

Received for publication July 30, 2002. Accepted for publication November 12, 2002.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
The species-specific inverse relation between atmospheric CO2 concentration and stomatal frequency for many woody angiosperm species is being used increasingly with fossil leaves to reconstruct past atmospheric CO2 levels. To extend our limited knowledge of the responsiveness of conifer needles to CO2 fluctuations, the stomatal frequency response of four native North American conifer species (Tsuga heterophylla, Picea glauca, Picea mariana, and Larix laricina) to a range of historical CO2 mixing ratios (290 to 370 ppmV) was analyzed. Because of the specific mode of leaf development and the subsequent stomatal patterning in conifer needles, the stomatal index of these species was not affected by CO2. In contrast, a new measure of stomatal frequency, based on the number of stomata per millimeter of needle length, decreased significantly with increasing CO2. For Tsuga heterophylla, the stomatal frequency response to CO2 changes in the last century is validated through assessment of the influence of other biological and environmental variables. Because of their sensitive response to CO2, combined with a high preservation capacity, fossil needles of Tsuga heterophylla, Picea glauca, P. mariana, and Larix laricina have great potential for detecting and quantifying past atmospheric CO2 fluctuations.

Key Words: CO2 • conifers • leaf morphology • North America • quantification methods • stomata • stomatal density


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
The current rise in atmospheric CO2 levels has stimulated research on physiological responses of land plants to changes in the concentration of this greenhouse gas. In growth experiments with elevated CO2 levels, the number of leaf stomata decreased for 40% of the species studied (Woodward and Bazzaz, 1988 ; Woodward and Kelly, 1995 ; Royer, 2001 ). This inverse relation between stomatal frequency and atmospheric CO2 is even more apparent under lower-than-present experimental CO2 levels (Woodward, 1987 ; Woodward and Bazzaz, 1988 ; Malone et al., 1993 ; Royer, 2001 ). A genetic basis for the response is provided by the identification of a CO2-sensitive gene in Arabidopsis thaliana that is involved in the control of stomatal development (Gray et al., 2000 ).

The relation between stomatal frequency and atmospheric CO2 concentration (mixing ratio; measured as parts per million by volume [ppmV]) is confirmed by studies of leaves that were grown under the lower CO2 levels of the past 100 yr. Such leaves can be obtained from herbarium collections or subfossil leaf assemblages. The main advantage of historical leaf material over experimental material is that the stomatal frequency responses to CO2 have occurred under natural conditions. Effects of experimental single-step CO2-enrichment (typically CO2-doubling) on seedlings may not be representative of long-term responses to incremental increases of 1 to 2 ppmV CO2 per year or per growing season.

Notably leaves of woody angiosperm taxa show a negative correlation between mean stomatal frequency and the historical CO2 rise from 280 to 360 ppmV (Woodward, 1987 ; Peñuelas and Matamala, 1990 ; Paoletti and Gellini, 1993 ; Kürschner et al., 1996 ; Wagner et al., 1996 ; Wagner, 1998 ). Modelled response curves for species of Betula and Quercus predict different response limits to a CO2 increase (~400 and ~340 ppmV, respectively), indicating that nonlinear stomatal frequency responses vary from one tree species to another (Kürschner et al., 1997 ). The models also suggest that the maximum effect of the current CO2 increase on stomatal frequency has already been reached, thereby explaining the common lack of stomatal frequency responses in CO2-doubling experiments.

Conifers and broad-leaved angiosperms have basic differences in stomatal formation and leaf development (Esau, 1977 ). In angiosperms, stomata and epidermal cells are initiated at multiple points on the developing leaf surface, resulting in a relatively random distribution of stomata on the entire leaf surface. Epidermal cells and stomata in conifer needles are all initiated at the base of the needle, developing henceforth in longitudinal files during needle growth (Croxdale, 2000 ). Conifers have a less-pronounced decrease in stomatal conductance than angiosperms in response to long-term, elevated CO2 exposure (Saxe et al., 1998 ; Medlyn et al., 2001 ). Because stomatal conductance is in part determined by the number of stomata on the leaf surface (Jones, 1992 ), in theory different responses could be related to differences in the adjustment of stomatal frequency to CO2 levels.

To date, only six experiments have been performed to test whether the typical arrangement of stomata in conifer needles affects stomatal frequency responses to atmospheric CO2 (Royer, 2001 ). A frequency decrease in Pinus sylvestris under elevated CO2 (Beerling, 1997 ; Lin et al., 2001 ) suggests that conifers also may have the capacity to adjust their stomatal numbers. However, stomatal frequency in other species (Pinus pinaster, P. banksiana, P. palustris, Pseudotsuga menziesii, and Picea abies) is not reduced at elevated CO2 levels (Stewart and Hoddinott, 1993 ; Guehl et al., 1994 ; Dixon et al., 1995 ; Pritchard et al., 1998 ; Apple et al., 2000 ). This lack of reduction could indicate that in many conifer species, stomatal frequency is not very sensitive to CO2 changes. On the other hand, the controlled-environment experiments may have been performed under CO2 levels that are well above the response limit to a CO2 increase.

In addition to experimental studies, the CO2 responsiveness of two conifer species has been confirmed by analyses of fossil and herbarium materials. Van de Water et al. (1994) have reported elevated stomatal frequencies in fossil needles of Pinus flexilis corresponding to a CO2 range from 180 to 290 ppmV. Recently, in an herbarium study on needles of Metasequoia glyptostroboides, stomatal frequency decreased linearly over a CO2 range from 310 to 340 ppmV (Royer et al., 2001 ). Additional experiments on Metasequoia under CO2 concentrations of 430 and 790 ppmV showed a nonlinear response, similar to the response curves of angiosperm species, levelling off at about 500 ppmV. Obtaining these response curves is hampered by the fact that most experiments involve a single-step doubling of CO2. To compare the CO2 response curves of conifers to those of angiosperms, the stomatal frequency must be measured over a gradual increase in CO2.

Historical data sets enable the construction of accurate response curves that can be used as training sets for the quantification of past atmospheric CO2 levels on the basis of stomatal frequency analysis of fossil leaves (Van der Burgh et al., 1993 ; Beerling et al., 1995 ; Kürschner et al., 1996 ; Rundgren and Beerling, 1999 ; Wagner et al., 1999 ; McElwain et al., 2002 ). Considering the dominance of conifers in temperate and boreal forest ecosystems, the ability to use fossil conifer needles for quantifying past CO2 levels would greatly extend the spatial and temporal coverage of such reconstructions. In North America, well-preserved needles of Tsuga heterophylla, Picea glauca, P. mariana, and Larix laricina are often present in Quaternary lake and peat deposits (Dunwiddie, 1986 ; Cwynar, 1987 ; Mayle and Cwynar, 1995 ).

To assess the CO2 responsiveness of conifer needles, as well as their potential as biosensors of paleo-atmospheric CO2 fluctuations, the present study focuses upon the stomatal frequency response of these species to a range of historical CO2 values (290–370 ppmV). Customized counting strategies are introduced to account for the typical stomatal patterning in conifers needles.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Materials
Stomatal frequency was analyzed for four native North American conifer species: Tsuga heterophylla, Picea mariana, P. glauca, and Larix laricina. The measurements were performed at two institutes: (1) on Tsuga heterophylla needles at the Laboratory of Palaeobotany and Palynology (Utrecht University) and (2) on the latter three species at the Department of Animal and Plant Sciences (University of Sheffield).

Some of the needles of Tsuga heterophylla used in the study were obtained from the collection of the National Herbarium of the Netherlands, complemented with needles from living trees collected in the field in 1998 and 2000. The material originates from 13 localities in the Northwest Pacific area of the USA and Canada; 8–11 needles from each locality were processed (Fig. 1). Additionally, assemblages of subfossil Tsuga heterophylla needles were derived from an 8-cm-long peat core taken at the margin of Jay Bath, a pond in Mount Rainier National Park (Washington, USA; 46°46' N, 121°46' W). Eleven of the annual layers in this core, spanning a period from 1980 to 1998, yielded 3–5 Tsuga heterophylla needles each. A total of 160 Tsuga heterophylla needles from herbaria, living trees, and the peat core were analyzed for stomatal frequency.



View larger version (29K):
[in this window]
[in a new window]
 
Fig. 1. Source localities of the Tsuga heterophylla needles in the Northwest Pacific region. BC = British Columbia (Canada), WA = Washington (USA), and OR = Oregon (USA)

 
The material of Picea glauca, P. mariana, and Larix laricina originates from herbarium sheets of the Royal Botanical Gardens Kew, cultivated trees in the United Kingdom, and living trees at various locations in the USA and Canada. Each locality provided 1–4 needles, for a total of 69 Picea glauca/mariana needles and 31 Larix laricina needles used in this study.

Methods
The needles of Tsuga heterophylla were bleached with a 4% natriumhypochloride solution to remove the mesophyll. The remaining cuticle was then stained with safranin and mounted in glycerin jelly on a microscopic slide. Computer-aided measuring of epidermal cell parameters on needle cuticles of Tsuga heterophylla was performed on a Leica Quantimet 500C/500+ Image Analysis system (Wetzlar, Germany). Stomata on Picea and Larix needles were counted directly from whole unprepared needles using a Leica epifluorescence microscope. Regression analysis and Student's t test were performed using SPSS 8.0 for Windows statistical software (Chicago, Illinois, USA).

The presence of tephra from the 1980 eruption of Mt. St. Helens (Mullineaux, 1974 , 1996 ) enabled a precise age assessment of the annual layers apparent in the peat core from Jay Bath.

The atmospheric CO2 mixing ratios for the past 100 yr that were used in the response curves were derived from instrumental measurements at Mauna Loa since 1957 (Keeling and Whorf, 2002 ) and shallow Antarctic ice cores (Etheridge et al., 1996 ; Indermühle et al., 1999 ). Records of temperature and precipitation at Mount Rainier were obtained from the IRI/LDEO Climate Data Library (website: http://ingrid.ldeo.columbia.edu).

Expression of stomatal frequency
Stomatal frequency is conventionally expressed as stomatal density (the number of stomata per unit leaf area) or as stomatal index (the proportion of stomata expressed as a percentage of total epidermal cells). Stomatal index is generally favored because it accounts for the influence of lateral cell expansion related to water availability on stomatal density (Salisbury, 1927 ; Kürschner et al., 1996 ; Kürschner, 1997 ).

The stomata on conifer needles are arranged in rows, occurring in either single files (Larix laricina) or grouped in bands (Tsuga heterophylla, Picea glauca/mariana; Fig. 2). Therefore, in contrast to broad-leaved species, stomatal placement and occurrence on conifer needles is not well characterized through the measurement of either stomatal density or index in small counting fields.



View larger version (67K):
[in this window]
[in a new window]
 
Fig. 2. Abaxial surface of a Tsuga heterophylla needle viewed with a dissecting stereomicroscope. Band = band of stomatal rows on each side of the midvein. Dotted lines: approximate positions of the transects used for measuring the number of stomatal rows in both bands (number of rows). White squares: approximate positions of the counting fields used to measure stomatal density (SD)

 
In Larix laricina, the arrangement of stomata in single files across the entire needle surface hampers the conventional use of small counting fields to determine stomatal density because variation in the number of rows might not be fully captured. Therefore, stomatal frequency in Larix laricina is expressed as the total number of stomata per millimeter of needle length. Also, the typical stomatal pattern precludes the application of stomatal index.

The lower surface of the hypostomatal Tsuga heterophylla needles displays two broad bands of stomatal rows on each side of the central vein (Fig. 2), while four stomatal bands are present on the surface of the four-sided amphistomatal Picea glauca/mariana needles. In this way, the leaf surface is divided into stomate-free (midvein, leaf margins) and stomatal regions (the bands). Consequently, the stomatal frequency on the needles of these species not only varies in terms of stomatal density/index within the bands, but also by the extent of the stomatal regions on the needle surface. Thus, for Tsuga heterophylla and Picea glauca/mariana stomatal index and density within the bands were measured, as well as the width of the bands and the total number of stomata per millimeter of needle length.

To account for the varied stomatal patterning in conifer needles, the following customized stomatal quantification methods were applied in this study.7

Stomatal density (SD) was measured as the number of stomata per square millimeter of leaf area. On Tsuga heterophylla needles, SD was measured in 10 counting fields (0.057 mm2) within the stomatal bands at a magnification of 500x. For Picea glauca/mariana needles, SD was measured on the abaxial surface in 2–10 counting fields spanning the entire width of a band along the length of a calibrated 0.800-mm scale bar at a magnification of 200x.

Stomatal index (SI) was measured as (the number of stomata divided by the number of stomata plus epidermal cells) x100. On Tsuga heterophylla needles, SI was calculated based on stomatal and epidermal cell counts in 12 counting fields spanning the width of a band along a 0.400-mm needle length. On Picea glauca/mariana needles, SI was calculated based on stomatal and epidermal cell counts in 2–10 counting fields spanning the width of a band along the length of a calibrated 0.800-mm scale bar.

Stomatal number per Length (SNL) was measured as (the number of abaxial stomata plus the number of adaxial stomata) divided by needle length in millimeters. On Larix laricina needles, SNL was measured in 2–10 counting fields along the length of a calibrated 0.800-mm scale bar at a magnification of 200x.

Stomatal rows (SR) was measured as the number of stomatal rows in both stomatal bands. This method quantifies the extent of the stomatal regions on the Tsuga and Picea needles. Band width is highly dependent of the number of rows in a band, as measured in 35 perfectly preserved modern Tsuga heterophylla needles (n = 700; Pearson correlation, 0.953). Band width is expressed as number of rows rather than absolute width (in millimeters) because the former is easier to measure (especially in less well-preserved [fossil] needles). For Tsuga heterophylla, the number of rows in both bands was counted along 10 transects perpendicular to the needle axis and restricted to the middle third part of the needle because the number of rows decreases strongly at the tip and the base of the needles (Fig. 2). For the amphistomatal Picea glauca/mariana, mean number of rows per needle (on both sides) was determined.

Stomatal density per length (SDL) was determined using the equation SDL = SD x SR. Stomatal density per length combines stomatal density within the bands and the width of the bands, providing a measure of the total number of stomata per millimeter of needle length in Tsuga heterophylla and Picea glauca/mariana needles. Because SR is the number of rows instead of the true band width (in millimeters), SDL represents a general indication of the number of stomata per millimeters needle length, capable of reflecting changes in stomatal frequency. The expression of band width in number of rows rather than absolute distance also explains the apparent discrepancy of expressing SDL per square millimeters instead of per millimeter.

True stomatal density per length (TSDL) was determined using the equation TSDL = SD x band width (in millimeters). The number of stomata per millimeter of needle length (TSDL) can be approximated more accurately when the band width is expressed in millimeters instead of number of rows. The TSDL was calculated only for Tsuga heterophylla, using the linear relation between SR and band width (in millimeters) as measured in 35 perfectly preserved modern Tsuga heterophylla needles (n = 700; Pearson correlation, 0.953).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
The SI calculated for needles of Tsuga heterophylla from four samples of herbarium material (Fig. 3A) did not change significantly over a CO2 range from 290 to 367 ppmV (R2 = 0.22; P = 0.531). By contrast, three other measured stomatal parameters decreased in relation to the CO2 rise from 290 to 367 ppmV: SD (Fig. 3B), SR (Fig. 3C), and TDSL (Fig. 3D). The SD decreased from 205 to 177 stomata/mm2 (R2 = 0.5077; P < 0.001), a response rate of 1.99% decrease per 10 ppmV of CO2 increase. The SR shows a response rate of –2.85% per 10 ppmV of CO2 increase, based on the decrease in number of rows from 17.5 to 14.0 (R2 = 0.3389; P = 0.002). The TSDL decreases from 237 to 164 stomata/mm; a response rate of –4.49% per 10 ppmV of CO2 increase (R2 = 0.5873; P < 0.001).



View larger version (16K):
[in this window]
[in a new window]
 
Fig. 3. Response of stomatal parameters of Tsuga heterophylla to a CO2 increase from 290 to 370 ppmV. Black diamonds represent subfossil and modern needles from Jay Bath (Mount Rainier, Washington), open diamonds modern and herbarium needles from other localities. Error bars indicate ±1 SE. Solid lines indicate best fit in classical regression analysis. (A) Response of stomatal index (SI) of Tsuga heterophylla to a CO2 increase from 290 to 370 ppmV. Mean SIs do not change significantly (SI = 0.0038 x [CO2] + 16.19; r2 = 0.22; P = 0.531). (B) SD: stomatal density (SD = –0.3786 x [CO2] + 316.25; r2 = 0.4892; P < 0.001). (C) SR: number of stomatal rows in both bands (SR = –0.0514 x [CO2] + 32.875; r2 = 0.3588; P = 0.002). (D) TSDL: true stomatal density per millimeter of needle length (TSDL = –1.0501 x [CO2] + 549.67; r2 = 0.5873; P < 0.001)

 
The highly similar stomatal frequencies of Picea glauca and P. mariana allows their treatment as a single group. No significant relationship was observed between SI and CO2 mixing ratio for these species (data not shown). Mean SDL in P. glauca and P. mariana decreased from 106 to 50 stomata/mm2 (Fig. 4) in response to the CO2 rise from 288 to 360 ppmV (R2 = 0.5048; P = 0.001), resulting in a response rate of –7.34% per 10 ppmV of CO2 rise for these species.



View larger version (11K):
[in this window]
[in a new window]
 
Fig. 4. Response of stomatal density per millimeter of needle length (SDL) on Picea glauca/mariana needles to a CO2 increase from 283 to 360 ppmV based on both adaxial and abaxial measurements. Error bars indicate ±1 SE. Solid line indicates best fit in classical regression analysis (SDL = –0.7949 x [CO2] + 335.41; r2 = 0.5048; P = 0.001)

 
Larix laricina responds to the CO2 increase from 283 to 319 ppmV with a decrease in mean SNL from 24.4 to 16.5 stomata/mm2 (Fig. 5), a response rate of –8.09% per 10 ppmV of CO2 increase (R2 = 0.3884; P = 0.023).



View larger version (9K):
[in this window]
[in a new window]
 
Fig. 5. Response of stomatal number per millimeter of needle length (SNL) on Larix laricina needles to a CO2 increase from 283 to 319 ppmV. Error bars indicate ±1 SE. Solid line indicates best fit in classical regression analysis (SNL = –0.1983 x [CO2] + 80.413; r2 = 0.3884; P = 0.023)

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Counting strategies
In Tsuga heterophylla and Picea glauca/mariana, stomatal frequency is significantly reduced as a response to increasing CO2. The SD decreases as does the stomatal region. The SI, on the other hand, although commonly used as the most sensitive indicator for CO2-related change in stomatal frequency in broad-leaved angiosperm species, does not decrease as atmospheric CO2 concentrations increase. Thus, the observed decrease in SD on the needles must be accompanied by a proportionally equal reduction in epidermal cells, keeping SI constant. This discrepancy in response between SD and SI can be explained by the mode of stomatal formation in conifers.

Conifer leaves develop linearly from a single growth center at the needle base, similar to monocots, instead of from multiple point sources on the leaf as in broad-leaved angiosperms (Croxdale, 2000 ). Epidermal cells appear in longitudinal files, with the cells at the tip of the leaf maturing first. Stomatal formation in conifers takes place during a late stage of leaf growth. Guard cells of Pseudotsuga menziesii, for example, do not become apparent until the needle has reached two-thirds of its final size (Owens, 1968 ). When the stomatal precursor cells form in specific rows, the epidermal cells in the adjacent files of cells divide asymmetrically to produce subsidiary cells lying directly next to the stomatal precursor cells and epidermal cells in the adjacent cell file (Fig. 6; Croxdale et al., 1992 ; Larkin et al., 1997 ). In many conifer species and monocots, the latter epidermal cells can continue to divide a specific number of times creating a species-specific fixed number of epidermal cells in the "stomatal-epidermal complex" (Tomlinson, 1974 ). These complexes can contain between four and 12 epidermal cells each, depending on the species (Florin, 1931 ). When stomata are densely packed, as in Tsuga and Picea, nearly all epidermal cells encountered in the stomatal bands belong to the stomatal-epidermal complexes, except for a few very elongated pavement cells in rows of epidermal cells (Fig. 7). In these cases of densely packed stomata, the number of epidermal cells in the bands will mainly depend on the number of stomata present, keeping SI constant. The SD reflects stomatal initiation rate much better than the SI. Because epidermal cells from the stomatal complexes should be excluded from the ratio, the application of SI as an expression of stomatal frequency in conifers should be restricted to species with a significant number of pavement epidermal cells. This phenomenon may also explain in part why little or no stomatal response has been observed in grasses growing in elevated CO2. The number of subsidiary cells can be affected by CO2 levels (Boetsch et al., 1996 ), but pavement and stomatal complex epidermal cells in mature conifer leaves cannot always be easily distinguished (Florin, 1931 ). Therefore, the size of the stomatal complexes and the number of pavement epidermal cells in a conifer species should be taken into account before using a SD- or SI-based counting strategy.



View larger version (19K):
[in this window]
[in a new window]
 
Fig. 6. Stomatal patterning during leaf development in grasses (note: stomatal patterning in conifers and monocots is highly comparable [Croxdale, 2000]). (i) Adjacent cell files on a developing leaf. In the middle cell file asymmetric divisions occur with the stomatal initial cell forming towards the leaf tip and the larger cell product forming closer to the leaf base, resulting in an alternation of stomatal initials, which develop directly into a guard mother cell, and future epidermal neighbor cells. New transverse cell walls are offset from walls in adjacent files. (ii) Subsidiary cells form in asymmetric divisions in adjacent cell files on both sides of the guard mother cell. (iii) The guard mother cell divides symmetrically to form two guard cells that surround a pore. The lower neighbor cell originated in the asymmetric division in (i); the darker line outlines the two cells produced by that division. (Reproduced with permission from the American Society of Plant Biologists: Larkin et al., 1997 )

 


View larger version (76K):
[in this window]
[in a new window]
 
Fig. 7. Abaxial cuticle of Tsuga heterophylla. Stomatal complexes are outlined in black. (P) indicates "pavement" epidermal cells (those not originating as part of a stomatal complex)

 
When the other stomatal parameters are considered, TSDL/SDL appears to be the quantification method most sensitive to CO2. The SD in Tsuga heterophylla decreases with increasing CO2, but at a lower response rate than TSDL, while the band width changes too, only with much greater variability than TSDL. In P. glauca and P. mariana, the same picture emerges; SR has a strong response to CO2, but with high variability. The significance increases when SD and SR are combined in SDL (J. C. McElwain, personal observation). The TSDL/SDL is thus proposed as the most appropriate method to use in cases where the stomatal rows in a hypo- or amphistomatal conifer species are grouped in bands, as in Tsuga heterophylla and Picea glauca/mariana. As a further validation for the use of the number of rows as a proxy for band width, TSDL (based on actual measurements of band width in millimeters) was calculated for 35 modern Tsuga heterophylla needles. Measured TSDL and SDL correlate extremely well (Pearson correlation coefficient of 0.951), which indicates that the number of rows can be applied as a very useful proxy for band width. Actually, the number of rows might be an even more consistent proxy of stomatal region than the measured band width itself. When a few stomata lie away from the main band or large areas are devoid of stomata, the measurement of band width in stomatal rows rather than of absolute distance will minimize the overestimation of the stomatal region.

Concerning conifer taxa with single stomatal rows (such as Pinus and Larix), care should also be taken to capture the variability in the number of rows as well as the variability in the number of stomata per row. This can be achieved by measuring SNL, the total (abaxial and adaxial) number of stomata per millimeter of needle length. If small counting fields are used (less than half of the needle width) for conventional SD measurements, minor variations in the number of rows will not be reflected, and stomatal row numbers should be included in the quantification of stomatal frequency.

Validation of the stomatal frequency response to CO2
In the present study SD-based counting strategies are applied instead of SI-based methods. The SD depends not only on stomatal initiation rate, but also on lateral cell expansion after stomatal formation, which can be strongly influenced by environmental factors other than atmospheric CO2. To determine if the observed changes in stomatal frequency in the conifers can be ascribed to the rise in CO2 over the last century, the potential influence of environmental and biological variations (geography, altitude, light regime, humidity, temperature, and ontogeny) on stomatal frequency has to be evaluated.

Geography
Because the conifer needles in this study were collected in varying geographical regions, interpretation of the observed responses could be complicated by differences in climate at the source localities and genotypic dissimilarities between different populations. Contrasting sampling strategies for Tsuga heterophylla vs. Picea glauca/mariana and Larix laricina were used to deal with the potential influence of variation in climate and populations on the stomatal frequency.

All Tsuga heterophylla needles (from living trees, herbarium sheets, and subfossil peat) were collected in the Northwest Pacific region of the USA/Canada and have thus grown under similar temperate conditions with high precipitation, minimizing the potential influence of large climatic variations. The modern samples (Fig. 3: 363 and 367 ppmV) were collected from living trees at several sites in British Columbia and Washington. The considerable variation in stomatal frequency in these modern samples may reflect intrinsic variability between trees, phenotypical responses to micro-climatic differences between these sites, and differences in local CO2 mixing ratios present at forested sites (Tarnawski et al., 1994 ). Despite the variation in stomatal frequency, the TSDL (the most sensitive quantification method) of the modern samples is consistently low and the decreasing trend in mean stomatal frequency in the data set is substantial and significant.

Subfossil and modern needles, which were derived from a single population of trees around Jay Bath over the last 20 yr, are indicated by black diamonds in Fig. 3D. The TSDL changes in this restricted data set (within a single population) do not differ significantly from the TSDL decrease observed in the total data set, which includes data points from other localities (and populations). Thus, there are no indications that the observed stomatal frequency changes in Tsuga heterophylla reflect variations in stomatal frequency between different populations.

A reverse sampling strategy was employed for Picea glauca/mariana and Larix laricina. Herbarium material from these species was chosen from a wide range of localities in northern USA and Canada, within the species' natural climatic ranges and tolerances. Herbarium leaves of cultivated trees from Scotland, North America, and England were also included in the calibration data sets because no significant differences in stomatal frequencies were observed between trees growing in cultivation and in nature under the same CO2 levels. By maximizing variation in climate and populations in the data set, the potential effect of biological and environmental factors other than CO2 on stomatal frequency of Picea and Larix is expressed within the statistical confidence intervals of the regression data set. Because the decrease in stomatal frequency appears both within a single population and in a data pool of maximized environmental variation, the stomatal response to CO2 of the studied conifer species is unlikely to be an artefact caused by differences in populations and climate associated with the geographical variation in the data sets.

Altitude
Although the CO2 mixing ratio in air remains constant over altitudinal gradients, the CO2 partial pressure (as measured in pascals) is lower at higher elevation because of the decrease in air pressure. Stomatal frequency responds to altitude-controlled changes in CO2 partial pressure when CO2 mixing ratio is unaltered (Woodward and Bazzaz, 1988 ). To check whether the observed stomatal frequency decrease in Tsuga heterophylla could be related to altitudinal differences in the data set (sea level to 1600 m), an alternative response curve was plotted using only needles grown at altitudes in the range of Jay Bath (1200–1600 m; Fig. 8). The stomatal frequency response in Tsuga heterophylla needles from this restricted altitudinal range is not significantly different from the response in the original data set with needles from a broader altitudinal range (Fig. 3D). For Picea glauca, P. mariana, and Larix laricina, leaf material from herbarium sheets collected from high elevations was not included in the data sets to avoid the effects of decreased CO2 partial pressure with increasing elevation on stomatal frequency. Thus, the observed stomatal frequency response in the conifer species is not significantly affected by the altitudinal variation in the data set.



View larger version (13K):
[in this window]
[in a new window]
 
Fig. 8. Response of true stomatal density per millimeter of needle length (TSDL) of Tsuga heterophylla to a CO2 increase from 290 to 370 ppmV. The material depicted in this figure are those needles from Fig. 3D that grew at 1300–1600 m altitude. Black diamonds represent subfossil and modern needles from Jay Bath (Mount Rainier, Washington, USA); open diamonds represent herbarium needles from other localities. Error bars indicate ±1 SE. Solid line indicates best fit in classical regression analysis (TSDL = –1.0256 x [CO2] + 545.95; r2 = 0.6223; P < 0.001)

 
Light intensity
Because light intensity is known to strongly influence epidermal morphology and SD (Nordhausen, 1903 ; Lichtenthaler, 1985 ; Kürschner, 1997 ), additional stomatal frequencies were measured for Tsuga heterophylla needles that were grown under sunny and shady light regimes at the botanical gardens in Utrecht (Netherlands). No significant difference in either SD, SR, or TSDL could be detected in the needles grown under contrasting light regimes (shade: 177 ± 36 stomata/mm, sun: 174 ± 18 stomata/mm; P = 0.883 for TSDL; L. L. R. Kouwenberg, unpublished data). Hence, the observed relation between CO2 and TSDL in Tsuga heterophylla is not significantly affected by the varying light intensity under which the herbarium and fossil needles may have developed. The observed TSDL values of these needles are highly comparable to the modern material from the Northwest Pacific region of North America.

Water availability
Water availability determines cell expansion and thus stomatal density on leaf surfaces (Salisbury, 1927 ; Tichá, 1982 ). Therefore it must be considered whether the observed SD decrease in the data set could be a result of local changes in soil water availability during the last century. Annual and spring precipitation records measured at Longmire (Mount Rainier) from 1930 to 2002 were compared to the stomatal frequency records of Tsuga heterophylla (Fig. 9A). These records do not show any long-term unidirectional increase in precipitation that might have triggered the observed reduction in stomatal frequency.



View larger version (21K):
[in this window]
[in a new window]
 
Fig. 9. (A) Precipitation at Longmire (Mount Rainier, Washington, USA), 1930–2002. The upper curve represents annual precipitation; the thick black line is the 19-yr running mean to emphasize the long-term trends. The lower curve represents spring precipitation (February–May); the thick gray line indicates the 19-yr running mean. (B) Mean annual temperature at Longmire, 1915–2002. Thick black line represents the 19-yr running mean to visualize the long-term trend. Data in (A) and (B) from the IRI/LDEO Climate Data Library (website: http://ingrid.ldeo.columbia.edu). (C) The atmospheric CO2 mixing ratios from instrumental measurements at Mauna Loa, Hawaii, USA, since 1957 (black line; Keeling and Whorf, 2002 ) and Antarctic ice cores 1900–1990 (white squares are measurements from Law Dome [Etheridge et al., 1996 ]; black squares are measurements from Taylor Dome [Indermühle, 1999 ]). (D) True stomatal density per millimeter of needle length (TSDL) for Tsuga heterophylla. Black diamonds represent subfossil and modern needles from Jay Bath (Mount Rainier); open diamonds represent modern and herbarium needles from other localities. Error bars indicate ±1 SE. Thick line represents the long-term trend (drawn by visual match)

 
Temperature
Experiments with Betula pubescens indicate that large temperature changes (from 12 to 30°C) can cause a significant increase in stomatal formation as reflected in SI (Wagner, 1998 ). Thus, monthly temperature records from Longmire were also compared with the stomatal frequency records of Tsuga heterophylla (Fig. 9B). Both spring and annual temperatures increased 0.8°C between 1915 and 1940 but showed no further long-term unidirectional change. Thus, differences in growing temperature are highly unlikely to have influenced the decrease in stomatal frequency.

Ontogeny
Marked variations in SD are reported for angiosperm leaves during leaf ontogeny (Tichá, 1982 ). Considering the specific leaf development and stomatal initiation in conifers, ontogenetical differences in stomatal frequency on conifer needles could be expected to be less pronounced. Measurements of SD, SR, and TSDL on developing Tsuga heterophylla needles from the same branch were indeed not significantly different in these stomatal frequency parameters for young and fully grown needles (mean needle length of 7.2 and 12.0 mm, respectively; L. L. R. Kouwenberg, unpublished data).

Stomatal frequency response rates
The four conifer species analyzed in this study show a significant reduction in stomatal frequency as a response to a CO2 rise of 80 ppmV over the last century. The response rates, a 4.49% decrease in TSDL per 10 ppmV of CO2 increase for Tsuga heterophylla and a 7.34% and 8.09% decrease in SDL for Picea mariana/glauca and Larix laricina, respectively, are comparable with the SI response rates for angiosperm tree taxa commonly used for CO2 analysis of fossil leaves (Betula: –5.47% of the SI per 10 ppmV of CO2, Quercus –3.82%). Thus, over a CO2 range between 290 and 370 ppmV, stomatal frequency adjustment in conifers can occur at a similar rate as in woody angiosperms. Furthermore, earlier suggestions that the response of conifers in general would level off at CO2 values of 280 ppmV (Van de Water et al., 1994 ) are refuted by the observation that both Tsuga heterophylla and Picea glauca/mariana have not reached their response limit yet at the current CO2 level of 370 ppmV. Thus, similar to angiosperms, conifers (Pinus, Metasequoia, Tsuga, Picea, Larix) may have species-specific response ranges. Several conifer species still respond to increases in CO2 above the response limits for angiosperm species used for stomatal frequency analysis, making them very suitable for the reconstruction of paleo-atmospheric CO2 concentrations well above present levels (Royer et al., 2001 ).

Conclusions
The capability of conifers to adjust their stomatal frequency to changes in atmospheric CO2 mixing ratios has now been confirmed in four more conifer species. Tsuga heterophylla, Picea glauca, P. mariana, and Larix laricina show a strong reduction in stomatal frequency in accordance with a change in historical CO2 levels from 290 to 370 ppmV. Because of their sensitive response over a broad range of CO2 levels, combined with a high preservation capacity, fossil needles of the studied conifer species show great potential for paleo-atmospheric CO2 reconstructions. Considering the dominance of conifers in temperate and boreal forest ecosystems (in both the Northern and Southern Hemisphere), the spatial and temporal coverage of such reconstructions may thus be markedly extended. This is corroborated already by results of ongoing research in the Northwest Pacific region of the USA (L. L. R. Kouwenberg et al., unpublished manuscript and Atlantic Canada (McElwain et al., 2002 ). Furthermore, the observation that all four conifer species have not yet reached their response limits to CO2 suggest that paleo-CO2 estimates derived from conifers rather than angiosperms may be preferable during those times in the Earth's history when greenhouse gas concentrations were much higher than present (such as during the Cretaceous).

The specific stomatal patterning on conifer needles has precluded the application of stomatal index measurements as a CO2-sensitive expression of stomatal frequency. New quantification methods, based on the number of stomata per millimeter of needle length, are proposed as the most accurate measure of stomatal frequency in conifer species. In general, stomatal quantification methods should be tailored to specific leaf development and subsequent stomatal patterning of the species studied.


    FOOTNOTES
 
1 The authors thank the curators of the Herbarium at the Royal Botanic Gardens Kew and the National Herbarium of the Netherlands for loan of herbarium material. We thank the Mount Rainier National Park for permission to sample at Mount Rainier and Peter Dunwiddie, Johan van der Burgh, and Ton van Druten for assistance in the field. We acknowledge the constructive remarks of two anonymous reviewers. This research was supported by The Leverhulme Trust (F/118/AP), the Royal Society and the Netherlands Organisation for Scientific Research (NWO project 750.198.07). This is the Netherlands School for Sedimentary Geology (NSG) publication no. 20030201. Back

6 Author for reprint requests (L.L.R.Kouwenberg{at}bio.uu.nl ) Back

7 Methods have been developed independently in Utrecht and Sheffield. Back


    LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Apple M. E. D. M. Olszyk D. P. Ormrod J. Lewis D. Southworth D. T. Tingey 2000 Morphology and stomatal function of Douglas fir needles exposed to climate change: elevated CO2 and temperature. International Journal of Plant Sciences 161: 127-132[CrossRef][ISI][Medline]

Beerling D. J. 1997 Carbon isotope discrimination and stomatal responses of mature Pinus sylvestris L. trees exposed in situ for three years to elevated CO2 and temperature. Acta Oecologica 18: 697-712[CrossRef]

Beerling D. J. H. H. Birks F. I. Woodward 1995 Rapid late-glacial atmospheric CO2 changes reconstructed from the stomatal density record of fossil leaves. Journal of Quaternary Science 10: 379-384

Boetsch J. J. Chin M. Ling J. Croxdale 1996 Elevated carbon dioxide affects the patterning of subsidiary cells in Tradescantia stomatal complexes. Journal of Experimental Botany 47: 925-931

Croxdale J. L. 2000 Stomatal patterning in angiosperms. American Journal of Botany 87: 1069-1080[Abstract/Free Full Text]

Croxdale J. B. Johnson J. Smith B. Yandell 1992 Stomatal patterning in Tradescantia: an evaluation of the cell lineage theory. Developmental Biology 93: 1078-1082

Cwynar L. 1987 Fire and forest history of the North Cascade Range. Ecology 68: 791-802[CrossRef][ISI]

Dixon M. D. Le Thiec J. P. Garrec 1995 The growth and gas exchange responses of soil-planted Norway spruce (Picea abies (L.) Karst) and red oak (Quercus rubra L.) exposed to a naturally occurring drought. New Phytologist 129: 265-273[CrossRef][ISI]

Dunwiddie P. 1986 A 6000-year record of forest history on Mount Rainier, Washington. Ecology 67: 58-68

Esau K. 1977 Anatomy of seed plants. John Wiley and Sons, New York, New York, USA

Etheridge D. M. L. P. Steele R. L. Langenfels R. J. Francey J. M. Barnola V. I. Morgan 1996 Natural and anthropogenic changes in atmospheric CO2 over the last 1000 years from air in Antarctic ice and firn. Journal of Geophysical Research 101: 4115-4128[CrossRef]

Florin R. 1931 Untersuchungen zur Stammengeschichte der Coniferales und Cordaitales. Kungliga Svenska Vetenskapsakademiens Handlingar 10: 1-588

Gray J. E. G. H. Holroyd F. Van der Lee P. C. Sijmons F. I. Woodward W. Schuch A. M. Hetherington 2000 HIC; a gene involved in controlling stomatal development in responses to changes in atmospheric CO2 concentration. Nature 408: 713-715[CrossRef][Medline]

Guehl J. M. C. Picon G. Aussenac P. Gross 1994 Interactive effects of elevated CO2 and soil drought on growth and transpiration efficiency and its determinants in two European tree species. Tree Physiology 14: 707-724[ISI][Medline]

Indermühle A. T. F. Stocker F. Joos H. Fischer H. J. Smith M. Wahlen B. Deck D. Mastroianni J. Tschumi T. Blunier R. Meyer B. Stauffer 1999 Holocene carbon-cycle dynamics based on CO2 trapped in ice at Taylor Dome, Antarctica. Nature 398: 121-126[CrossRef]

Jones H. G. 1992 Plants and microclimate. Cambridge University Press, Cambridge, UK

Keeling C. D. T. P. Whorf 2002 Atmospheric CO2 concentrations–Mauna Loa Observatory, Hawaii, 1958–2001. Website: http://cdiac.esd.ornl.gov/ndps/ndp001.html

Kürschner W. M. 1997 The anatomical diversity of recent and fossil leaves of the durmast oak (Quercus petraea Lieblein/Quercus pseudocastanea Goeppert): implications for their use as biosensors of paleoatmospheric CO2 levels. Review of Palaeobotany and Palynology 96: 1-30

Kürschner W. M. J. van der Burgh H. Visscher D. L. Dilcher 1996 Oak leaves as biosensors of late Neogene and early Pleistocene paleoatmospheric CO2 concentrations. In R. Poore and L. C. Sloan [eds.], Pliocene climates Marine Micropaleontology 27: 299-312

Kürschner W. M. F. Wagner E. H. Visscher H. Visscher 1997 Predicting the response of leaf stomatal frequency to a future CO2-enriched atmosphere: constraints from historical observations. Geologische Rundschau 86: 512-517[CrossRef][ISI]

Larkin J. C. D. M. Marks J. Nadeau F. Sack 1997 Epidermal cell fate and patterning in leaves. Plant Cell 9: 1109-1120[CrossRef][ISI][Medline]

Lichtenthaler H. K. 1985 Differences in morphology and chemical composition of leaves grown at different light intensities and qualities. In N. R. Baker, W. J. Davies, and C. K. Ong [eds.], Control of leaf growth, 201–221. Cambridge University Press, Cambridge, UK

Lin J. M. E. Jach R. Ceulemans 2001 Stomatal density and needle anatomy of Scots pine (Pinus sylvestris) are affected by elevated CO2. New Phytologist 150: 665-674[CrossRef][ISI]

Malone S. R. H. S. Mayeux H. B. Johnson H. W. Polley 1993 Stomatal density and aperture length in four plant species grown across a subambient CO2 gradient. American Journal of Botany 80: 1413-1418[CrossRef][ISI]

Mayle F. E. L. Cwynar 1995 Impact of Younger Dryas cooling event upon lowland vegetation of maritime Canada. Ecological Monographs 65: 129-154[CrossRef]

McElwain J. C. D. J. Beerling F. E. Mayle 2002 Stomatal evidence for a decline in atmospheric CO2 concentration during the Younger Dryas stadial: a comparison with Antarctic ice core records. Journal of Quaternary Science 17: 21-29[CrossRef][ISI]

Medlyn B. E. et al 2001 Stomatal conductance of forest species after long-term exposure to elevated CO2 concentration: a synthesis. New Phytologist 149: 247-264[CrossRef][ISI]

Mullineaux D. R. 1974 Pumice and other pyroclastic deposits in Mount Rainier National Park, Washington. Geological Survey Bulletin 1326: 1-83

Mullineaux D. R. 1996 Pre-1980 Tephra-fall deposits erupted from Mount St. Helens, Washington. U.S. Geological Survey Professional Paper 1563

Nordhausen M. 1903 Über Sonnen- und Schattenblätter. Berichte der Deutschen Botanischen Gesellschaft 21: 30-45

Owens J. N. 1968 Initiation and development of leaves in Douglas fir. Canadian Journal of Botany 46: 271-278

Paoletti E. R. Gellini 1993 Stomatal density in beech and holm oak leaves collected over the last 200 years. Acta Ecologica 14: 173-178

Peñuelas J. R. Matamala 1990 Changes in N and S leaf content, stomatal density and specific leaf area of 14 plant species during the last three centuries of CO2 increase. Journal of Experimental Botany 41: 1119-1124[Abstract/Free Full Text]

Pritchard S. G. C. Mosjidis C. M. Peterson G. B. Runion H. H. Rogers 1998 Anatomical and morphological alterations in longleaf pine needles resulting from growth in elevated CO2: interactions with soil resource availability. International Journal of Plant Sciences 159: 1002-1009[CrossRef][ISI]

Royer D. L. 2001 Stomatal density and stomatal index as indicators of paleoatmospheric CO2 concentration. Review of Palaeobotany and Palynology 114: 1-28[CrossRef][ISI][Medline]

Royer D. L. S. L. Wing D. J. Beerling D. W. Jolley P. L. Koch L. J. Hickey R. A. Berner 2001 Paleobotanical evidence for near-present-day levels of atmospheric CO2 during part of the Tertiary. Science 292: 2310-2313[Abstract/Free Full Text]

Rundgren M. D. J. Beerling 1999 A Holocene CO2 record from the stomatal index of subfossil Salix herbacea L. leaves from northern Sweden. Holocene 9: 509-513[Abstract/Free Full Text]

Salisbury E. J. 1927 On the causes and ecological significance of stomatal frequency, with special reference to the woodland flora. Philosophical Transactions of the Royal Society of London B 216: 1-65

Saxe H. D. S. Ellsworth J. Heath 1998 Tree and forest functioning in an enriched CO2 atmosphere. New Phytologist 139: 395-436[CrossRef][ISI]

Stewart J. D. J. Hoddinott 1993 Photosynthetic acclimation to elevated atmospheric CO2 and UV radiation in Pinus banksiana. Physiologia Plantarum 88: 493-500[CrossRef]

Tarnawski M. G. T. G. A. Green B. Buedel A. Meyer H. Zellner O. L. Lange 1994 Diel changes of atmospheric CO2 concentration within, and above, cryptogam stands in a New Zealand temperate rain forest. New Zealand Journal of Botany 32: 329-336[ISI]

Tichá I. 1982 Photosynthetic characteristics during ontogenesis of leaves. 7. Stomata density and sizes. Photosynthetica 16: 375-471[ISI]

Tomlinson P. B. 1974 Development of the stomatal complex as a taxonomic character in the monocotyledons. Taxon 23: 109-128

Van der Burgh J. H. Visscher D. L. Dilcher W. M. Kürschner 1993 Paleoatmospheric signatures in Neogene fossil leaves. Science 260: 1788-1790[Abstract/Free Full Text]

Van de Water P. K. S. W. Leavitt J. L. Betancourt 1994 Trends in stomatal density and 13C/12C ratios of Pinus flexilis during last glacial-interglacial cycle. Science 264: 239-242[Abstract/Free Full Text]

Wagner F. 1998 The influence of environment on the stomatal frequency in Betula. Ph.D. thesis, Laboratory of Palaeobotany and Palynology, Utrecht University, Utrecht, The Netherlands. LPP Contributions Series 9: 1-102

Wagner F. R. Below P. de Klerk D. L. Dilcher H. Joosten W. M. Kürschner H. Visscher 1996 A natural experiment on plant acclimation: lifetime stomatal frequency response of an individual tree to annual atmospheric CO2 increase. Proceedings of the National Academy of Science USA 93: 11705-11708[Abstract/Free Full Text]

Wagner F. S. J. P. Bohnke D. L. Dilcher W. M. Kürschner B. van Geel H. Visscher 1999 Century-scale shifts in early Holocene atmospheric CO2 concentration. Science 284: 1971-1973[Abstract/Free Full Text]

Woodward F. I. 1987 Stomatal numbers are sensitive to increases in CO2 concentration from pre-industrial levels. Nature 327: 617-618[CrossRef][ISI]

Woodward F. I. F. A. Bazzaz 1988 The response of stomatal density to CO2 partial pressure. Journal of Experimental Botany 39: 1771-1781[Abstract/Free Full Text]

Woodward F. I. C. K. Kelly 1995 The influence of CO2 concentration on stomatal density. New Phytologist 131: 311-327[CrossRef][ISI]




This article has been cited by other articles:


Home page
Reviews in Mineralogy and GeochemistryHome page
L. L.R. Kouwenberg, W. M. Kurschner, and J. C. McElwain
Stomatal Frequency Change Over Altitudinal Gradients: Prospects for Paleoaltimetry
Reviews in Mineralogy and Geochemistry, October 1, 2007; 66(1): 215 - 241.
[Abstract] [Full Text] [PDF]


Home page
GeologyHome page
M. Haworth, S. P. Hesselbo, J. C. McElwain, S. A. Robinson, and J. W. Brunt
Mid-Cretaceous pCO2 based on stomata of the extinct conifer Pseudofrenelopsis (Cheirolepidiaceae)
Geology, September 1, 2005; 33(9): 749 - 752.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Bot.Home page
F. Wagner, D. L. Dilcher, and H. Visscher
Stomatal frequency responses in hardwood-swamp vegetation from Florida during a 60-year continuous CO2 increase
Am. J. Botany, April 1, 2005; 92(4): 690 - 695.
[Abstract] [Full Text] [PDF]


Home page
GeologyHome page
Atmospheric CO2 fluctuations during the last millennium reconstructed by stomatal frequency analysis of Tsuga heterophylla needles
Geology, January 1, 2005; 33(1): 33 - 36.



Home page
ANN BOT (LOND)Home page
L. L. R. KOUWENBERG, W. M. KURSCHNER, and H. VISSCHER
Changes in Stomatal Frequency and Size During Elongation of Tsuga heterophylla Needles
Ann. Bot., October 1, 2004; 94(4): 561 - 569.
[Abstract] [Full Text] [PDF]


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