Leaf anatomy plasticity of Alchornea triplinervia (Euphorbiaceae) under distinct light regimes in a Brazilian montane Atlantic rain forest

 Alchornea triplinervia trees occur in a montane Atlantic rainforest at the Macaé de Cima State Ecological Reserve (Nova Friburgo, RJ, Brazil). They are found in two adjacent secondary forests at distinct successional stages: a late-secondary (shaded and humid) and an early-secondary forest (drier, higher light intensities reaching the understory). Leaf samples collected in these environments were compared in regard to various anatomic parameters. Histochemical tests, anatomic measurements and counting indicated no significant variation for the basic leaf anatomy. Nevertheless, as compared to the late-secondary forest, the leaves of the early-secondary forest individuals showed an increased bulk of sclerenchyma at the main nervation and petiole, gelatinous fibers with hygroscopic walls all along the central vascular system, thicker cuticle, a higher abaxial stomatal density, abundance of tannin in the mesophyll, a significantly thicker palisade and spongy parenchyma, and compactation of the spongy parenchyma. This anatomical variation indicates an increased xeromorphism of leaves under the drier conditions of the early-secondary forest. Leaf anatomy plasticity may contribute to the high ecological plasticity of this widespread neotropical species. Received: 11 November 1996 / Accepted: 5 February 1997     

Trait variation and trait stability in common marmosets (Callithrix jacchus) inhabiting ecologically distinct habitats in northeastern Brazil

Understanding the set of factors that promote and constrain a species’ ability to exploit ecologically distinct habitats is central for addressing questions of intraspecific variability in behavior and morphology. In this study, we compared newly collected data with published data on body measurements, group size and composition, daily path length, home range, and reproductive output in wild common marmosets naturally inhabiting two contrasting environments in northeastern Brazil: the Atlantic Forest (AF), which is characterized by high biodiversity and reduced seasonality in food availability and the Caatinga (CAT), which is characterized by a severe hot and dry season lasting from 5 to 11 months, drought‐resistant plant species, and reduced primary productivity. Despite marked differences in ecological conditions, CAT marmosets and AF marmosets differed minimally in daily path length, home range, reproductive output, and infant survivorship. CAT marmosets were found to live in smaller groups containing fewer adult females than AF marmosets, and also were characterized by a greater surface area to body mass ratio, a trait that may represent an adaptation to the hot and dry conditions of the Caatinga. We propose that in conjunction with body mass reduction, minor adjustments in behavior, the exploitation of cacti as a source of water and nutrients, and access to exudates as a dependable year‐round food resource, common marmosets successfully used the same adaptive pattern to maintain high reproductive output and infant survivorship in exploiting these two ecologically distinct environments.     

Intra- and interspecific variation in genome size in Lathyrus (Leguminosae)

2C nuclear DNA amounts for 24 species of Lathyrus (section Lathyrus ) were determined using flow cytometry. A greater than two-fold variation was observed, ranging from 10.2 pg in L. basalticus to 24.2 pg in L. latifolius. In general, the perennial species had greater DNA amounts than the annuals. Significant intraspecific variation was observed in five species of Lathyrus (from 10.1% in L. annuus to 28% in L. tingitanus). A positive correlation was observed between DNA values obtained by flow cytometry and those previously determined by microdensitometry. Finally, the distribution of DNA amounts in species within section Lathyrus appears to be continuous.     

Phenology and fecundity in 11 sympatric pioneer species of Macaranga (Euphorbiaceae)in Borneo

Reproductive traits of tropical tree species vary predictably in relation to successional stage, but this variation may be due to the species' phylogenetic histories rather than selective pressures imposed by regeneration requirements. Reproductive phenology, tree size at the onset of reproduction, and fecundity of 11 sympatric, closely related Macaranga species were studied to investigate within-species variation in reproductive traits in relation to resource availability, and among-species variation in relation to other life-history traits (shade tolerance, seed size and maximum tree size, H(max)) and consequently the requirements for forest-gap colonization. Nine species reproduced in synchronous episodes, and two species reproduced continuously over 32 mo. Episodic reproduction was most intense in 1992 following a severe drought. For several species, reproductive trees had greater light availability, lower fecundity in lower light levels, and lower growth rates than nonreproductive trees, reflecting resource-limited reproduction. Among species, H(max) was negatively correlated with shade tolerance and seed size. Tree size at the onset of reproduction and fecundity was strongly linked to this axis of life-history variation, but phenological pattern was not. Absolute tree size at the onset of reproduction was positively correlated with H(max) and negatively correlated with shade tolerance. Relative size at reproductive onset was not correlated with shade tolerance or H(max). Fecundity ranged four orders of magnitude among species and was correlated positively with H(max) and negatively with seed size and shade tolerance. The interrelationships among these reproductive and other life-history traits are strongly correlated with the species' requirements for gap colonization.     

Extensive intraspecific chloroplast DNA (cpDNA) variation in the alpine Draba aizoides L. (Brassicaceae): haplotype relationships and population structure

Chloroplast DNA (cpDNA) sequence variation is currently the most widely used tool for the inference of phylogenetic relationships among plants at all taxonomic levels. Generally, noncoding regions tend to evolve faster than coding sequences and have recently been applied to the study of phylogenetic relationships among closely related taxa. An implicit assumption of many of these studies is that intraspecific cpDNA variation is either absent or low and therefore will not interfere with the reconstruction of interspecific relationships. A survey of cpDNA sequence variation in the common alpine plant species Draba aizoides L. was undertaken to assess levels of intraspecific cpDNA sequence variation. These levels were compared to levels of interspecific sequence divergence between D. aizoides and related alpine Draba species. Intraspecific cpDNA sequence divergence was extensive in D. aizoides, and intraspecific differences were often larger than interspecific differences. cpDNA haplotype relationships were explored using a maximum parsimony approach and minimum-spanning networks. Results from both methods were largely congruent but comparisons provided interesting insights into the presumed evolutionary history of cpDNA haplotypes. A combined effect of cpDNA introgression and complex lineage sorting was inferred to explain the pattern of cpDNA variation found in D. aizoides. Our results suggest that intraspecific cpDNA variation can be extensive and that intraspecific variation needs to be taken into account when inferring phylogenetic relationships among closely related taxa.     

Examining variation in the leaf mass per area of dominant species across two contrasting tropical gradients in light of community assembly

Understanding variation in key functional traits across gradients in high diversity systems and the ecology of community changes along gradients in these systems is crucial in light of conservation and climate change. We examined inter‐ and intraspecific variation in leaf mass per area (LMA) of sun and shade leaves along a 3330‐m elevation gradient in Peru, and in sun leaves across a forest–savanna vegetation gradient in Brazil. We also compared LMA variance ratios (T‐statistics metrics) to null models to explore internal (i.e., abiotic) and environmental filtering on community structure along the gradients. Community‐weighted LMA increased with decreasing forest cover in Brazil, likely due to increased light availability and water stress, and increased with elevation in Peru, consistent with the leaf economic spectrum strategy expected in colder, less productive environments. A very high species turnover was observed along both environmental gradients, and consequently, the first source of variation in LMA was species turnover. Variation in LMA at the genus or family levels was greater in Peru than in Brazil. Using dominant trees to examine possible filters on community assembly, we found that in Brazil, internal filtering was strongest in the forest, while environmental filtering was observed in the dry savanna. In Peru, internal filtering was observed along 80% of the gradient, perhaps due to variation in taxa or interspecific competition. Environmental filtering was observed at cloud zone edges and in lowlands, possibly due to water and nutrient availability, respectively. These results related to variation in LMA indicate that biodiversity in species rich tropical assemblages may be structured by differential niche‐based processes. In the future, specific mechanisms generating these patterns of variation in leaf functional traits across tropical environmental gradients should be explored.     

Development and plasticity of endangered shrub Lindera melissifolia (Lauraceae) seedlings under contrasting light regimes

Lindera melissifolia (Walt.) Blume seedlings were raised in a growth chamber to determine the effects of light availability on shoot growth pattern, and basic leaf and stem growth. Lindera melissifolia seedlings exhibited a sympodial shoot growth pattern for 3 months following emergence from the soil medium, but this pattern was characterized by a reduction in leaf blade area approximately 30 days after emergence, followed by increases in leaf blade area. Seedlings receiving low light were 76% taller than seedlings receiving high light. Seedlings receiving low light also had larger leaf blade dimensions, blade area, seedling leaf area, and greater mass. Seedlings raised in high light had a greater proportional distribution of biomass in the roots, suggesting possible water stress from greater vapor pressure deficits. Furthermore, these seedlings displayed sharp angles of blade inclination and blade folding – acclimation that reduces exposure to light and subsequent higher leaf temperatures in open environments. These differences in morphological response to light resulted in high phenotypic variability in L. melissifolia seedlings. Lindera melissifolia seedling development showed a brief period of phenotypic plasticity, followed by ontogenetic plasticity. The short period of phenotypic plasticity may, however, have profound ecological implications for the conservation and recovery of this federally endangered shrub. Further experimentation should take into account the development of ontogenetic standards for comparisons of plant traits in addition to temporal standards.     

Three species of Pagodia (Corynexochida,Trilobita) from the Furongian (Cambrian) of northern Anhui,China and their intraspecific variation

More than one thousand cranidia, librigenae, and pygidia of the trilobite Pagodia were collected from the Cambrian (upper Jiangshanian through Stage 10) at Fenghuangshan in northern Anhui, China. They were identified as Pagodia lotos Walcott, 1905, Pagodia depressa (Walcott, 1905), and Pagodia bia Walcott, 1905. Based on one rare nearly complete exoskeleton of P. depressa and comparisons among the three species of Pagodia, the redefinition for those species is attempted as follows. Pagodia lotos: Glabella rectangular, slightly convergent anteriorly; anterior border ridge-like, anterior border furrow deep; pygidium semicircular, axis tapered rearward, with 4 axial rings and a terminal piece. P. depressa: Cranidium subquadrate, highly convex; glabella cylindrical, convergent anteriorly; pygidium semicircular. P. bia: Glabella rectangular, constricted in the midlength; fixigena wide (tr.); anterior border very narrow (sag.) and ridge-like; pygidium semicircular with wide borders; pygidial axis ended closer to the pygidial border furrow. We also discuss their intraspecific variations, commonly seen on the exfoliated specimens: palpebral ridge evident or obliterated, the length (sag.) of anterior area variable, the decoration of exoskeleton surface (with granules variable in size, convex or pitted).     

Mitochondrial DNA variation in Chinese and Indian rhesus macaques (Macaca mulatta)

DNA was extracted from the buffy coats or serum of 212 rhesus macaques (Macaca mulatta) sampled throughout the species' geographic range. An 835 base pair (bp) fragment of mitochondrial DNA (mtDNA) was amplified from each sample, sequenced, aligned, and used to estimate genetic distances from which phylogenetic trees were constructed. A tree that included sequences from rhesus macaques whose exact origins in China are known was used to determine the regional origin of clusters of haplotypes, or haplogroups, defined by the trees. Indian rhesus sequences formed one large homogeneous haplogroup with very low levels of nucleotide diversity and no geographic structure, and a second much smaller haplogroup apparently derived from Burma. The sequences from Burma and eastern and western China were quite divergent from those in the major haplogroup of India. Each of these sequences formed separate clusters of haplotypes that exhibited far greater nucleotide diversity and/or population structure. Correspondingly, sequences from Indian rhesus macaques that are considered to represent different subspecies (based on morphological differences) were intermingled in the tree, while those from China reflected some, but not all, aspects of subspecific taxonomy. Regional variation contributed 72% toward the paired differences between sequences in an analysis of molecular variance (AMOVA), and the average differences between the populations of eastern and western China were also statistically significant. These results suggest that Indian and Chinese rhesus macaques were reproductively isolated during most, if not all, of the Pleistocene, during which time Indian rhesus macaques experienced a severe genetic bottleneck, and that some gene flow westward into India was subsequently reestablished. Samples from breeding centers in three different provinces of China included sequences from rhesus macaques that originated in both eastern (or southern) and western China, confirming anecdotal reports that regional breeding centers in China exchange breeding stock. Genetic differences among rhesus macaques (even those acquired from the same regional breeding center) that originate in different geographic regions and are employed as subjects in biomedical experiments can contribute to phenotypic differences in the traits under study.     

Leaf traits and performance vary with plant age and water availability in Artemisia californica

Background and AimsLeaf functional traits are strongly tied to growth strategies and ecological processes across species, but few efforts have linked intraspecific trait variation to performance across ontogenetic and environmental gradients. Plants are believed to shift towards more resource-conservative traits in stressful environments and as they age. However, uncertainty as to how intraspecific trait variation aligns with plant age and performance in the context of environmental variation may limit our ability to use traits to infer ecological processes at larger scales.MethodsWe measured leaf physiological and morphological traits, canopy volume and flowering effort for Artemisia californica (California sagebrush), a dominant shrub species in the coastal sage scrub community, under conditions of 50, 100 and 150 % ambient precipitation for 3 years.Key ResultsPlant age was a stronger driver of variation in traits and performance than water availability. Older plants demonstrated trait values consistent with a more conservative resource-use strategy, and trait values were less sensitive to drought. Several trait correlations were consistent across years and treatments; for example, plants with high photosynthetic rates tended to have high stomatal conductance, leaf nitrogen concentration and light-use efficiency. However, the trade-off between leaf construction and leaf nitrogen evident in older plants was absent for first-year plants. While few traits correlated with plant growth and flowering effort, we observed a positive correlation between leaf mass per area and performance in some groups of older plants.ConclusionsOverall, our results suggest that trait sensitivity to the environment is most visible during earlier stages of development, after which intraspecific trait variation and relationships may stabilize. While plant age plays a major role in intraspecific trait variation and sensitivity (and thus trait-based inferences), the direct influence of environment on growth and fecundity is just as critical to predicting plant performance in a changing environment.     

Geographic variation in otolith morphology among freshwater populations of Aphanius dispar (Teleostei,Cyprinodontiformes) from the southeastern Arabian Peninsula

Aphanius Nardo is a large genus of teleost fishes in the Old World, with 19 described species. Several of these species have only recently been recognized and additional species can be expected from isolated populations in remote areas. We show here that otolith morphology and statistical analyses of otolith variables can contribute to the detection of genetic differentiation in Aphanius. We studied samples of eight Aphanius dispar populations from the southeastern part of the Arabian Peninsula. Two populations originate from freshwater habitats far inland and probably have been isolated since the late Holocene some 4,000 years ago, three populations come from freshwater habitats with occasional connections to the sea, one population originates from a coastal site, and two were artificially introduced populations. The coastal population is interpreted as possessing the basic otolith type of A. dispar. The basic otolith type also occurs in fishes from the freshwater habitat which is located closest to the coast. Otoliths from the two other freshwater populations with occasional connections to the coast differ slightly from the basic type. However, the two populations from the long‐term isolated freshwater habitats far inland show distinct morphological changes. Our results are consistent with the hypotheses that i) otolith morphology is primarily genetically determined and is little influenced by physical parameters of the habitat, and ii) isolated A. dispar populations may be capable of evolving into new species within short periods of time. J. Morphol., 2009. © 2008 Wiley‐Liss, Inc.     

A preliminary study of genetic variation in populations of Monstera adansonii var. klotzschiana (Araceae) from North-East Brazil, estimated with AFLP molecular markers

BACKGROUND AND AIMS: This study sought genetic evidence of long-term isolation in populations of Monstera adansonii var. klotzschiana (Araceae), a herbaceous, probably outbreeding, humid forest hemi-epiphyte, in the brejo forests of Ceará (north-east Brazil), and clarification of their relationships with populations in Amazonia and the Atlantic forest of Brazil. METHODS: Within-population genetic diversity and between-population dissimilarity were estimated using AFLP molecular markers in 75 individuals from eight populations located in Ceará, the Brazilian Atlantic Forest and Amazonia. KEY RESULTS: The populations showed a clinal pattern of weak genetic differentiation over a large geographical region (F(ST) = 0.1896). A strong correlation between genetic and geographical distance (Mantel test: r = 0.6903, P = 0.002) suggests a historical pattern of isolation by distance. Genetic structure analysis revealed at least two distinct gene pools in the data. The two isolated Ceará populations are significantly different from each other (pairwise Phi(PT) = 0.137, P = 0.003) and as diverse (Nei's gene diversity, average H(e) = 0.1832, 0.1706) as those in the Atlantic and Amazon forest regions. The population in southern Brazil is less diverse (Nei's gene diversity, average H(e) = 0.127) than the rest. The Ceará populations are related to those of the Atlantic forest rather than those from Amazonia (AMOVA, among-groups variation = 11.95 %, P = 0.037). CONCLUSIONS: The gene pools detected within an overall pattern of clinal variation suggest distinct episodes of gene flow, possibly correlated with past humid forest expansions. The Ceará populations show no evidence of erosion of genetic diversity, although this was expected because of their isolation. Their genetic differentiation and relatively high diversity reinforce the importance of conserving the endangered brejo forests.     

Inferring the patterns and causes of geographic variation in Ephippiger ephippiger (Orthoptera, Tettigoniidae) using geographical information systems (GIS)

Ephippiger ephippiger is a tcltigoniid cricket that is fairly widespread in Europe. In southern France, between the Alps and the Pyrenees, it is geographically highly variable lor a number of different traits including morphology, allozymcs, behaviour and genetic markers. This variation has resulted in considerable confusion over the taxonomic and evolutionary status of forms of the organism. Through the use of geographical information systems (GIS) it has been possible to integrate both trait and environmental data sets from a variety of sources. The GIS is used to interpolate trait and environmental surfaces that are then examined for coincidence. Four general spatial patterns are identified. These can be explained through a combination of secondary contact, environmental adaptation and drift. Covariance between trait and environmental surfaces point to patterns which are primary in origin, whereas other dines probably reflect secondary contact following divergence in rcfugiae. Multivariate statistics support the validity of a major subdivision identified by GIS. This study is an example of the considerable potential for GIS as an investigative tool in evolutionary studies.     

The geographic structure of morphological variation in eight species of fiddler crabs (Ocypodidae: genus Uca) from the eastern United States and Mexico

Species with larger geographic distributions are more likely to encounter a greater variety of environmental conditions and barriers to gene flow than geographically‐restricted species. Thus, even closely‐related species with similar life‐history strategies might vary in degree and geographic structure of variation if they differ in geographic range size. In the present study, we investigated this using samples collected across the geographic ranges of eight species of fiddler crabs (Crustacea: Uca) from the Atlantic and Gulf coasts of North America. Morphological variation in the carapace was assessed using geometric morphometric analysis of 945 specimens. Although the eight Uca species exhibit different degrees of intraspecific variation, widespread species do not necessarily exhibit more intraspecific or geographic variation in carapace morphology. Instead, species with more intraspecific variation show stronger morphological divergence among populations. This morphological divergence is partly a result of allometric growth coupled with differences in maximum body size among populations. On average, 10% of total within‐species variation is attributable to allometry. Possible drivers of the remaining morphological differences among populations include gene flow mediated by ocean currents and plastic responses to various environmental stimuli, with isolation‐by‐distance playing a less important role. The results obtained indicate that morphological divergence among populations can occur over shorter distances than expected based on dispersal potential. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 248–270.     

Determination of ploidy and nuclear DNA content in populations of Atriplex halimus (Chenopodiaceae)

Nuclear DNA contents were determined by flow cytometry for 20 populations of the perennial C₄ shrub Atriplex halimus L. (Chenopodiaceae) originating from the Mediterranean basin and Fuerteventura (Canary Islands). Two populations were also analysed for chromosome number: one (from Ibiza, Spain), with a 2C nuclear DNA content of 2.40 pg, was shown to be diploid (2 n  = 2 x  = 18), whilst the other (from Sicily, Italy), with 5.11 pg, was tetraploid (2 n  = 2 x  = 36). With respect to nuclear DNA content, two groups of populations were detected, diploids with 2.40–2.44 pg and tetraploids with 4.77–5.13 pg. The diploid populations were mainly from the western Mediterranean (Spain and France) and Fuerteventura, whereas tetraploids were generally, but not exclusively, from more arid areas in North Africa and the eastern Mediterranean. In general, the diploid and tetraploid populations corresponded to the subspecies halimus and schweinfurthii , respectively. For certain populations having morphologies intermediate between those considered typical of these two subspecies, nuclear DNA contents showed them to be tetraploid. There was significant variation in nuclear DNA content among the tetraploid populations, with greater values in the more easterly populations.  © 2005 The Linnean Society of London, Botanical Journal of the Linnean Society , 147 , 441−448.     

Functional connectivity and genetic diversity of Eulaema atleticana (Apidae,Euglossina) in the Brazilian Atlantic Forest Corridor: assessment of gene flow

The mobility and dispersal of organisms affect population genetics and dynamics, and consequently affect persistence and the risk of extinction. Thus, it is important to understand how organisms move in the fragmented landscapes in order to manage populations and predict the effects of habitat changes on species persistence. This study evaluated the functional connectivity of an orchid bee (Eulaema atleticana Nemésio, 2009) with a high fidelity to forest habitats in the Brazilian Atlantic Forest Corridor by analyzing genetic diversity, spatial genetic structure, and gene flow estimated from microsatellite and mitochondrial markers. Genetic diversity was not correlated with area of the forest fragments, or with forest isolation. At the mosaic scale, Eulaema atleticana showed no significant or low genetic differentiation, indicating genetic homogeneity among forest fragments. A previous field study indicated that Eulaema atleticana was one of the most sensitive Euglossina bees to forest fragmentation but the present molecular analyses demonstrates that current gene flow is sufficient to maintain genetic variability at the mosaic scale.     

Wide variation in spatial genetic structure between natural populations of the European beech (Fagus sylvatica) and its implications for SGS comparability

Identification and quantification of spatial genetic structure (SGS) within populations remains a central element of understanding population structure at the local scale. Understanding such structure can inform on aspects of the species' biology, such as establishment patterns and gene dispersal distance, in addition to sampling design for genetic resource management and conservation. However, recent work has identified that variation in factors such as sampling methodology, population characteristics and marker system can all lead to significant variation in SGS estimates. Consequently, the extent to which estimates of SGS can be relied on to inform on the biology of a species or differentiate between experimental treatments is open to doubt. Following on from a recent report of unusually extensive SGS when assessed using amplified fragment length polymorphisms in the tree Fagus sylvatica, we explored whether this marker system led to similarly high estimates of SGS extent in other apparently similar populations of this species. In the three populations assessed, SGS extent was even stronger than this previously reported maximum, extending up to 360 m, an increase in up to 800% in comparison with the generally accepted maximum of 30-40 m based on the literature. Within this species, wide variation in SGS estimates exists, whether quantified as SGS intensity, extent or the Sp parameter. Consequently, we argue that greater standardization should be applied in sample design and SGS estimation and highlight five steps that can be taken to maximize the comparability between SGS estimates.     

A pre‐adaptive approach for tropical forest restoration during climate change using naturally occurring genetic variation in response to water limitation

Effective reforestation of degraded tropical forests depends on selecting planting material suited to the stressful environments typical at restoration sites that can be exacerbated by increased duration and intensity of dry spells expected with climate change. While reforestation efforts in nontropical systems are incorporating drought‐adapted genotypes into restoration programs to cope with drier conditions, such approaches have not been tested or implemented in tropical forests. As the first effort to examine genetic variation in plant response to drought in a tropical wet forest, we established a watering experiment using five replicated maternal lines (i.e. seedlings from different maternal trees) of five dipterocarp species native to Borneo. Apart from the expected species level variation in growth and herbivory (3‐fold variation in both cases), we also found intraspecific variation so that growth in some cases varied 2‐fold, and herbivory 3‐fold, among genetically different maternal lines. In two species we found that among‐maternal line variation in growth rate was negatively correlated with tolerance to water limitation, that is, the maternal lines that performed the best in the high water treatment lost proportionally more of their growth during water limitation. We argue that selection for tolerance to future drier conditions is not only likely to impact population genetics of entire forests, but likely extends from forest trees to the communities of canopy arthropods associated with these trees. In tropical reforestation efforts where increased drought is predicted from climate change, including plant material resilient to drier conditions may improve restoration effectiveness.     

Spatial and temporal variation of soil solution chemistry and ion fluxes through the soil in a mature Norway Spruce (Picea abies (L.) Karst.) stand

In this study we investigated the spatial and temporal variation in soil solution chemistry and of water and ion fluxes through the soil in a forest ecosystem. Our aim was to evaluate the relevance of these variations for the accuracy of average areal soil solution concentrations and ion fluxes with seepage at 90 cm depth.Twenty spatially distinct subcompartments of approximately 1 m² were established within a mature stand of Norway spruce and ceramic suction lysimeters were installed at depths of 20, 35 and 90 cm. A tensiometer was placed close to each suction lysimeter, and one throughfall sampler was established for each subcompartment.Soil solution samples were analysed for major ions (H⁺, Na⁺, K⁺, Ca²⁺, Mg²⁺, Mn²⁺, Fe³⁺, Al³⁺, Cl^-, NO ₃ ^- , SO ₄ ^2- . We calculated water fluxes for each subcompartment separately by a numeric simulation of the soil water flux close to the lysimeters. The ion fluxes at each lysimeter were calculated by multiplying the simulated water fluxes with the ion concentrations on a fortnightly base. Averaging these 20 independent ion fluxes gave the areal average flux and an estimate of its statistical accuracy. The spatial variation of ion concentrations in the soil solution was high with coefficients of variance ranging from 5% to 128%. Part of the spatial variation was related to stem distance. Temporal variation of the concentrations was less than spatial for most ions. The spatial variation of water and ion fluxes with seepage was also substantial; for example the fluxes of SO ₄ ^2- -S calculated for each subcompartment ranged from 21 to 119 kg ha^-1 yr^-1, with an arithmetic average of 47 kg ha^-1 yr^-1. For H₂O, Mg²⁺, Cl^-, and SO ₄ ^2- , the spatial heterogeneity of seepage fluxes was largely explained by the heterogeneity of throughfall fluxes. No such relationship was found for nitrogen.Despite using 20 replicates, the 95% confidence intervals of the average annual areal fluxes with seepage were found to be 20–30% for most ions.