AJEV
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am. J. Enol. Vitic. 60:1:87-92 (2009)
Copyright © 2009 by the American Society for Enology and Viticulture.
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
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 Web of Science (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Steele, M. R.
Right arrow Articles by Merzlyak, M. N.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Steele, M. R.
Right arrow Articles by Merzlyak, M. N.
Agricola
Right arrow Articles by Steele, M. R.
Right arrow Articles by Merzlyak, M. N.

Research Note

Nondestructive Estimation of Anthocyanin Content in Grapevine Leaves

Mark R. Steele1,2, Anatoly A. Gitelson1,*, Donald C. Rundquist1 and Mark N. Merzlyak3

1 Center for Advanced Land Management Information Technologies, School of Natural Resources, University of Nebraska–Lincoln, Lincoln, NE 68583; 2 Agricultural Research and Development Center, University of Nebraska–Lincoln, Ithaca, NE 68033; and 3 Department of Physiology of Microorganisms, Faculty of Biology, Moscow State University, 119991, GSP-1 Moscow, Russia.

* Corresponding author (email: agitelson2{at}unl.edu)


    Abstract
 Top
 Abstract
 Materials and Methods
 Results and Discussion
 Conclusions
 Literature Cited
 
The anthocyanin (Anth) content in leaves provides valuable information about the physiological status of plants. Thus, there is a need for accurate, efficient, practical methodologies to estimate this biochemical parameter. Reflectance measurement is a means of quickly and nondestructively assessing leaf Anth content in situ. The objective of this study was to test the overall performance and accuracy of nondestructive techniques for estimating Anth content in grapevine leaves. Relationships were established between Anth content and four vegetation indices: NIR (near-infrared)/green, red/green, anthocyanin reflectance index (ARI, based on reflectances in bands within the green and the red-edge regions), and a modified anthocyanin reflectance index (MARI, based on reflectances in green, red edge, and NIR). The algorithms for Anth retrieval were calibrated. The accuracy of Anth prediction was evaluated using an independent data set containing sampled leaves from two field-grown grape cultivars (Saint Croix and Saint Pepin) with no adjustment of the coefficients after initial calibration. Although Anth in the validation data set was widely variable, from 3 to 45 nmol cm-2, the ARI and MARI algorithms were capable of accurately predicting Anth content in grapevine leaves with a root mean square error below 3 nmol cm-2 and 2.3 nmol cm-2, respectively. Such an approach has potential for developing simple hand-held field instrumentation for accurate nondestructive Anth estimation and for analyzing digital airborne or satellite imagery to assist in making informed decisions regarding vineyard management.

Key words: anthocyanins, grapes, leaves, nondestructive, reflectance


Recent interest in red anthocyanin (Anth) pigments is partially due to their pharmacological and nutritional effects and possible role in prevention of common chronic diseases, including atherogenesis, coronary heart diseases, thrombosis, and cancer. The polyphenols in red wine are believed to be the main contributors to the beneficial health effects of this beverage. The well-known French paradox is related to the low incidence of cardiovascular problems among the populace of France, despite their consumption of rich foods. The broad range of biological activities of anthocyanins is often described as their anti-oxidant properties (see review, Borkowski et al. 2005).

Anthocyanin biosynthesis is genetically determined, and plant leaves differ widely, according to species and cultivar, in their ability to synthesize Anth. In many plants, the bright red Anth coloration is abundant in young and senescing leaves (Lee 2007). Anth biosynthesis is often initiated due to drought, insect pests, potassium deficiency, extreme temperature, and excessive light. This behavior may allow the pigment to serve as an indicator of plant stress (e.g., Neill and Gould 1999). It is generally accepted that one key physiological function of Anth in higher plants is its photo-protective role. The pigment, which is localized in vacuoles of epidermal cells, serves as a filter: an internal light trap for excessive solar radiation (Chalker-Scott 1999, Close and Beadle 2003, Steyn et al. 2002). In particular, the photo-protective role of epidermal Anth in young grapevine leaves has been reported (Liakopoulos et al. 2006).

Anthocyanins have traditionally been determined through wet-chemical methods, including pigment extraction in a solvent, spectrophotometric determination of absorbance by the Anth solution, and conversion from measured absorbance to Anth content (e.g., Gitelson et al. 2001). However, laboratory procedures are laborious, time-consuming, and destructive to leaves. Spectral reflectance may provide viticulturists with an efficient, nondestructive method of monitoring Anth content and the associated physiological status of the fruit crop. Various models (called vegetation indices) have been developed for retrieving pigment content in leaves via spectral information (Richardson et al. 2002). Such indices, when applied to measurements acquired by hand-held field instrumentation or aerial imagery, should enhance the ability to make informed decisions regarding vineyard management. Recently, vegetation indices were designed to estimate Anth content in leaves nondestructively (Gamon and Surfus 1999, Gitelson et al. 2001, 2006, van den Berg and Perkins 2005). These indices are based on relationships involving reflectance at several specific wavelengths with varying sensitivity and response to changes in Anth content.

Gamon and Surfus (1999) used a ratio of reflectances at red ({rho}red) and green ({rho}green) wavelengths as a proxy of Anth content with the following:


Formula 1(1)

Anth absorbs in situ at ~550 nm (Gitelson et al. 2001), and the red peak of chlorophyll (Chl) absorption in situ is ~670 nm, but absorption by Chl also occurs in the green range of the spectrum (Gitelson et al. 2003). Thus, the {rho}red/{rho}green estimates Anth content by comparing reflectance in the red band of Chl absorption to reflectance in the green band where both Chl and Anth absorb.

Van den Berg and Perkins (2005) developed and tested a modified chlorophyll content meter (CCM-200) for measuring Anth in sugar maple leaves displaying fall colors. The device was optimized for Anth detection by replacing the light emitting diode centered at 655 nm with one centered at 530 nm. The authors used the ratio of absorbance ({alpha}) in the near-infrared band (940 nm) and the green (530 nm) as a proxy of Anth. An Anth content index (ACI) was suggested in the form: ACI = {alpha}green/{alpha}NIR.

An anthocyanin reflectance index (ARI) for estimating Anth content is based on hemispherical reflectance measurements of various types of leaves, compiled using an integrating sphere (Gitelson et al. 2001):


Formula 2(2)

ARI compares the reciprocal of reflectance in the green with that in the red-edge region and attributes the difference to absorption by Anth. The ARI was later modified (labeled MARI) to include NIR reflectance, {rho}NIR, as a means of accounting for variability in leaf thickness (Gitelson et al 2001, 2006):


Formula 3(3)

The goal of the current study was to investigate the performance of several reflectance-based nondestructive techniques for estimating Anth content in grapevine leaves, specifically, (1) to establish the relationship between the Anth content measured analytically and four selected vegetation indices (ACI, red/green, ARI, and MARI) and to calibrate the algorithms for Anth estimation; and (2) to validate the algorithms using an independent data set for two grape cultivars (Saint Croix and Saint Pepin), thus determining the accuracy of Anth-content prediction without adjustment of the coefficients after initial calibration. The reflectance measurements described in this work were carried out by means of a fiber optic attached to a leaf clip (i.e., no expensive integrating sphere). Our interest was in developing and testing a simple method for making pigment measurements in a routine, efficient manner in a vineyard setting.


    Materials and Methods
 Top
 Abstract
 Materials and Methods
 Results and Discussion
 Conclusions
 Literature Cited
 
Forty-two leaves were sampled during four field campaigns in the summer of 2006. Leaves were selected based on visual characteristics: selected leaves ranged from dark green with little red to completely red with little green. Field campaigns took place on: (1) 1 June 2006 (10 Saint Croix leaves); (2) 9 June 2006 (10 Saint Pepin leaves); (3) 19 Sept 2006 (12 Saint Croix leaves); and (4) 20 Sept 2006 (10 Saint Croix leaves). Both Saint Croix and Saint Pepin are French-hybrid wine cultivars grown in the Midwest where vinifera vines are not well-suited to climatic conditions. The selected leaves were detached from the vine, immediately sealed in a plastic bag with a small amount of water, and placed in a cooler with ice.

After field sampling was completed, the leaves were transported to the lab. When the coloration of the entire leaf was not uniform (as often occurred during and after veraison), areas of homogeneous color on each leaf were identified and delineated with a marker. Six thickness measurements per leaf were made with an Absolute ID-S 1012 Digimatic Indicator (Mitutoyo, Aurora, IL) and mean thickness was determined for each leaf. Care was taken to restrict the measurements to areas between veins.

Reflectance measurements.   Spectral-reflectance measurements for the leaves were collected for each of the two grape cultivars noted above using a clip with a 2.3- mm diam bifurcated fiber optic attached to an Ocean Optics USB2000 radiometer (Dunedin, FL) and an Ocean Optics LS-1 tungsten halogen light source. The USB2000 radiometer measures radiance with a spectral resolution of ~1.5 nm in a wavelength range from 350 to 1000 nm. The LS-1 light source uses a regulated power supply and a tungsten halogen filament bulb. The light source was turned on 15 minutes prior to scanning to allow the bulb and filament to warm and stabilize.

The plastic leaf clip, which positioned the fiber against individual grapevine leaves, consisted of a black polyvinyl chloride (PVC) attachment and the bifurcated glass fiber optic (transmissive between 400 and 1000 nm). The black PVC clip held each leaf at a 60° angle relative to the fiber to reduce specular reflectance from the surface of the leaf. The clip also held a black foam background, with a nominal reflectance of 3% within the spectral range of the instrument, upon which each leaf was placed during spectral sampling. The low reflectance properties of the background minimized extraneous reflectance from the reflected light being transmitted through the leaf.

The radiometer was calibrated before each data-collection session using a Labsphere Spectralon reference panel (North Sutton, NH) with a nominal reflectance of 99%. The reference panel was held tightly against the fiber optic, and a spectral scan was recorded. The sensor was operated by the CALMIT Data Acquisition Program (University of Nebraska-Lincoln), which uses a single calibration scan collected at the time nearest to the acquisition of the target scans to compute reflectance. The reflectance spectra were calculated as a ratio of leaf radiance to the radiance of the reference panel.

For accurate representation of reflectance of the marked area, six reflectance measurements were acquired for each leaf. The locations of spectral measurements were carefully distributed throughout the entire marked area superimposed on the leaf. The average of the six scans per sample was calculated to establish a single representative reflectance spectrum per leaf from which index values were calculated.

Pigment extraction.   After reflectance measurements were acquired, two or three discs (1 cm diam) were cut from the marked area on each leaf. Pigment contents were calculated in a two-step procedure. The punches were weighed and ground in 100% methanol using a mortar and pestle until the pulp turned white in color and all pigments were extracted. The resulting homogenate was centrifuged in test tubes for 6 min. Absorption spectra of the supernatants were recorded, and Chl-a, Chl-b, and carotenoid contents were calculated using absorption coefficients (Lichtenthaler 1987).

Once chlorophyll-absorption spectra were collected, HCI (final concentration ~0.1 %) was added to the solution, resulting in the appearance of the characteristic Anth coloration of the extract. Next, absorption spectra of the acidified extracts were taken and absorbance at 530 nm was corrected for the contribution of pheophytins (Lichtenthaler 1987). Anth concentrations were quantified using an absorption coefficient of 30 mM cm-1 at 530 nm (Strack and Wray 1989). Pigment content was expressed as a function of leaf area.

Calibration and validation.   Reflectance spectra and the corresponding measured pigment contents were combined and split into two data sets: calibration and validation. The data sets were arranged in ascending order, with odd samples assigned to calibration groups and even samples to validation groups.

ACI was tested using reflectances in green {rho}green (530 nm) and NIR {rho}NIR (940 nm) spectral bands. Thus, the modified ACI had the following form:


Formula 4(4)

Red/green, ARI, and MARI were calculated using average reflectance values for each leaf in following bands (Gamon and Surfus 1999; Gitelson 2006): {lambda}green = 540–560 nm; {lambda}red = 660–680 nm; {lambda}red edge = 690–710 nm; and {lambda}NIR = 760–800 nm. The index values from the calibration data set were regressed against the corresponding measured Anth contents to calibrate the algorithms. The calibrated algorithms were applied to predict Anth content using the validation data set. The predicted Anth contents were compared to Anth values determined analytically as part of the validation data set, and the root mean square error (RMSE) of Anth prediction was calculated.


    Results and Discussion
 Top
 Abstract
 Materials and Methods
 Results and Discussion
 Conclusions
 Literature Cited
 
Analytical laboratory measurements yielded a wide variety of pigment contents, with Anth ranging from 3.90 to 44.43 nmol cm-2 and total Chl ranging from 5.1 to 47.6 nmol cm-2, similar to ranges observed in other studies (Gamon and Surfus 1999, Gitelson et al. 2001, 2006). Data for the two sampled cultivars at various growth stages and for both the calibration and validation data sets are shown (Figure 1Go). There was a wide variation in content and overall composition of both pigments in the sampled leaves. A Chl to Anth ratio for each data set was calculated to underscore that variation: 0.11 to 6.08 for the data set acquired on 1 June; 0.11 to 8.7 for 9 June; 0.07 to 8.7 for 19 Sept; and 0.21 to 10.8 for 20 Sept. The best fit functions for the relationship between Chl and Anth for each data set and all data taken together are power functions of the form Chl = kxAnth-n. The determination coefficients for these relationships were all high: 0.73 for the 1 June data set; 0.82 for 9 June; 0.86 for 19 Sept; and 0.84 for 20 Sept. When all data were combined, the relationship was Chl = 90.5Anth-0.67 with R2 = 0.67. For individual and combined data sets, moderate-to-high Chl content (>25 nmol cm-2) corresponded to low Anth (<10 nmol cm-2) while leaves with high Anth had much lower Chl (<20 nmol cm-2).


Figure 1
View larger version (18K):
[in this window]
[in a new window]

 
Figure 1 Relationship between anthocyanin and total chlorophyll content in grapevine leaves. Data correspond to leaves from the two cultivars at different growth stages. Note the wide variation in pigment content and composition for leaves comprising each data set. The best fit functions of the relationship Chl vs. Anth for each data set (dashed lines) and all data combined (solid line) are shown.

 
The leaves studied were between 10 and 90 days old with leaf thickness ranging between 0.25 and 0.6 mm. Leaves with moderate-to-high Anth content were younger and thinner (between 0.25 and 0.35 mm) than mature leaves with high Chl content (thicker than 0.4 mm). NIR reflectance, a proxy of leaf thickness and density (Slaton et al. 2001), was much higher in green mature leaves with low Anth content and lower in red-tinted to reddish leaves with higher Anth content (Figure 2Go).


Figure 2
View larger version (23K):
[in this window]
[in a new window]

 
Figure 2 Reflectance spectra of leaves with different anthocyanin contents.

 
While reflectance in the red-edge region ~700 nm remains insensitive to changes in Anth content, reflectance in the green range (~550 nm) notably decreases (Figure 2Go) with increase in Anth because that pigment absorbs in-situ in the green range (Gitelson et al. 2001). However, the green reflectance is also affected by Chl content as it decreases with an increase in Chl. This overlapping of Chl and Anth absorption in the green spectral region was the main challenge in developing algorithms for nondestructive Anth estimation.

The modified ACI (Eq. 4) was related to Anth content in neither the calibration (Figure 3AGo, R2 = 0.06) nor validation (R2 = 0.0005, not shown) data sets. Relationships between red/green, ARI, and MARI vegetation indices (Eqs. 1, 2, and 3) and the analytically determined Anth content were established; all three co-varied in a statistically significant manner. These relationships yielded calibration of three algorithms for nondestructive Anth estimation.


Figure 3
View larger version (9K):
[in this window]
[in a new window]

 
Figure 3 Relationships between reflectance indices and anthocyanin content measured analytically in the lab for leaves comprising the calibration data set.

 
Red/green ratio (R2 = 0.54) (Figure 3BGo):


Formula 5(5)

Anthocyanin reflectance index (R2 = 0.91) (Figure 3CGo):


Formula 6(6)

Modified anthocyanin reflectance index (R2 = 0.93) (Figure 3DGo):


Formula 7(7)

Algorithms (Eqs. 1, 2, 3) were also calibrated separately for Saint Pepin leaves. The calibrations yielded relationships in which slopes and offsets were not statistically significant from Eq. 5 to 7 (not shown).

Because of the absence of a reliable close relationship between the modified ACI and Anth, we did not validate this index. The results of validating the other three algorithms (Figure 4Go) were consistent with the findings in the calibration phase. Anth prediction by the red/green was the poorest with the lowest determination coefficient (R2 = 0.72) and highest RMSE (7.63 nmol cm-2) (Figure 4AGo). The derived model was:


Figure 4
View larger version (13K):
[in this window]
[in a new window]

 
Figure 4 Anthocyanin prediction by vegetation indices.

 

Formula 8(8)

Anth prediction by ARI was much more accurate than the red/green, with an increase in accuracy more than two-fold (R2 = 0.96, RMSE < 2.93 nmol cm-2) (Figure 4BGo). The model was:


Formula 9(9)

MARI was the most accurate in Anth prediction (R2 = 0.97, RMSE < 2.23 nmol cm-2) (Figure 4CGo), and a model of the form:


Formula 10(10)

Importantly, the slope of the best-fit function of the relationship between Anth content predicted by MARI and the analytically measured was almost equal to one and the intercept was close to the origin (Figure 4CGo).

The poor performance of the modified ACI is understandable. In the modified ACI = {rho}NIR/{rho}530, the reflectance {rho}530 is the term sensitive to content of both pigments, Anth and Chl, and it decreases as their contents increase. In addition, NIR reflectance, {rho}NIR, is affected by leaf thickness (Slaton et al. 2001). Thus, the modified ACI depends upon three variables (Chl, Anth, and leaf thickness) and when they vary independently, it becomes insensitive to Anth. The close relationship between ACI and Anth, as found by others (van den Berg and Perkins 2005), may be explained by the very close relationship between Chl and Anth in the autumn sugar maple leaves studied. Analyzing table 1 in van den Berg and Perkins (2005), we found that Chl and Anth contents were related with R2 > 0.96. Thus, ACI was an accurate proxy of Anth as well as Chl (the determination coefficient of the polynomial relationship ACI vs. Chl was greater than 0.95). It appears that ACI is an effective tool for Anth estimation only under specific conditions: when Chl and Anth are closely related and leaf thickness does not vary.

The red/green ratio has the same limitation as ACI. When Chl and Anth contents are closely and inversely related, the ratio can be used to measure Anth content. However, when Chl and Anth contents are slightly related or not related at all, the red/green ratio depends upon pigment composition rather than on Anth content alone (Gitelson et al. 2001).


    Conclusions
 Top
 Abstract
 Materials and Methods
 Results and Discussion
 Conclusions
 Literature Cited
 
ARI and MARI were effective in the nondestructive prediction of anthocyanin in the leaves of two French-hybrid grapevine cultivars despite differences in pigment composition, leaf thickness, age, and pubescence. Spectral bands selected for this study were 10 nm in width and both algorithms produced very accurate prediction of Anth, ranging from 3.90 to 45 nmol cm-2, with root mean square error below 3 nmol cm-2 by ARI and below 2.25 nmol cm-2 by MARI. Thus, a radiometer with only two spectral bands (540–560 nm and 690–710 nm) is required for nondestructive Anth estimation using ARI. For precise Anth estimation by MARI, one additional spectral band in the NIR range (760–800 nm) is required. ARI and MARI were validated with two hybrid cultivars of grapes, but it does not appear that the algorithms are cultivar-specific among the two varieties examined. Thus, the algorithms are likely to allow accurate Anth determination in leaves of the pure V. vinifera vines used for wine production worldwide. More work should be done to test the algorithms in various plant species.

Manuscript submitted June 2008; revised September 2008

Accepted for publication October 2008


    Literature Cited
 Top
 Abstract
 Materials and Methods
 Results and Discussion
 Conclusions
 Literature Cited
 
Borkowski, T., H. Szymusiak, A. Gliszczynaka-Swiglo, I.M.C.M. Rietjens, and B. Tyrakowska. 2005. Radical scavenging capacity of wine anthocyanins is strongly pH-dependent. J. Agric. Food Chem. 53:5526–5534.[Web of Science][Medline]

Chalker-Scott, L. 1999. Environmental significance of anthocyanins in plant stress responses. Photochem. Photobiol. 70:1–9.[Web of Science]

Close, D.C., and C.L. Beadle. 2003. The ecophysiology of foliar anthocyanin. Bot. Rev. 69:149–161.

Gamon, J., and J. Surfus. 1999. Assessing leaf pigment content and activity with a reflectometer. New Phytol. 143:105–117.[Web of Science]

Gitelson, A.A., M.N. Merzlyak, and O.B. Chivkunova. 2001. Optical properties and nondestructive estimation of anthocyanin content in plant leaves. Photochem. Photobiol. 74:38–45.[Web of Science][Medline]

Gitelson, A.A., U. Gritz, and M.N. Merzlyak. 2003. Relationships between leaf chlorophyll content and spectral reflectance and algorithms for non-destructive chlorophyll assessment in higher plant leaves. J. Plant Physiol. 160:271–282.[Web of Science][Medline]

Gitelson, A.A., G.P. Keydan, and M.N. Merzlyak. 2006. Three-band model for noninvasive estimation of chlorophyll, carotenoids, and anthocyanin contents in higher plant leaves. Geophys. Res. Lett. 33:L11402.

Lee, D. 2007. Nature’s Palette: The Science of Plant Color. University of Chicago Press, Chicago.

Liakopoulos, G., D. Nikolopoulos, A. Klouvatou, K.A. Vekkos, Y. Manetas, and G. Karabourniotis. 2006. The photoprotective role of epidermal anthocyanins and surface pubescence in young leaves of grapevine (Vitis vinifera). Ann. Bot. 98:257–265.[Abstract/Free Full Text]

Lichtenthaler, H.K. 1987. Chlorophyll and carotenoids: Pigments of photosynthetic biomembranes. Meth. Enzymol. 148:331–382.

Neill, S., and K. Gould. 1999. Optical properties of leaves in relation to anthocyanin concentration and distribution. Can. J. Bot. 77:1777–1782.

Richardson, A.D., S. Duigan, and G. Berlyn. 2002. An evaluation of noninvasive methods to estimate foliar chlorophyll content. New Phytol. 153:185–194.[Web of Science]

Slaton, M., E. Hunt Jr., and W. Smith. 2001. Estimating near-infrared leaf reflectance from leaf structural characteristics. Am. J. Bot. 88:278–284.[Abstract/Free Full Text]

Steyn, W.J., S.J.E. Wand, D.M. Holcroft, and G. Jacobs. 2002. Anthocyanins in vegetative tissues: A proposed unified function in photoprotection. New Phytol. 155:349–361.[Web of Science]

Strack, D., and V. Wray. 1989. Anthocyanins. In Methods in Plant Biology. Vol. 1. Plant Phenolics. J.B. Harborne (ed.), pp. 325–356. Academic Press/Harcourt Brace Jovanovich, London.

van den Berg, A.K., and T.D. Perkins. 2005. Nondestructive estimation of anthocyanin content in autumn sugar maple leaves. HortScience 40:685–686.[Abstract/Free Full Text]




This article has been cited by other articles:


Home page
Am. J. Bot.Home page
A. A. Gitelson, O. B. Chivkunova, and M. N. Merzlyak
Nondestructive estimation of anthocyanins and chlorophylls in anthocyanic leaves
Am. J. Botany, October 1, 2009; 96(10): 1861 - 1868.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
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 Web of Science (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Steele, M. R.
Right arrow Articles by Merzlyak, M. N.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Steele, M. R.
Right arrow Articles by Merzlyak, M. N.
Agricola
Right arrow Articles by Steele, M. R.
Right arrow Articles by Merzlyak, M. N.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS