Integrating high and moderate spatial resolution image data to estimate forest age structure

The collection of spatial information to quantify changes to the state and condition of the environment is a fundamental component of conservation or sustainable utilization of tropical and subtropical forests, Age is an important structural attribute of old-growth forests influencing biological diversity in Australia eucalypt forests. Aerial photograph interpretation has traditionally been used for mapping the age and structure of forest stands. However this method is subjective and is not able to accurately capture fine to landscape scale variation necessary for ecological studies. Identification and mapping of fine to landscape scale vegetative structural attributes will allow the compilation of information associated with Montreal Process indicators lb and ld, which seek to determine linkages between age structure and the diversity and abundance of forest fauna populations. This project integrated measurements of structural attributes derived from a canopy-height elevation model with results from a geometrical-optical/spectral mixture analysis model to map forest age structure at a landscape scale. The availability of multiple-scale data allows the transfer of high-resolution attributes to landscape scale monitoring. Multispectral image data were obtained from a DMSV (Digital Multi-Spectral Video) sensor over St Mary's State Forest in Southeast Queensland, Australia. Local scene variance levels for different forest tapes calculated from the DMSV data were used to optimize the tree density and canopy size output in a geometric-optical model applied to a Landsat Thematic Mapper (TU) data set. Airborne laser scanner data obtained over the project area were used to calibrate a digital filter to extract tree heights from a digital elevation model that was derived from scanned colour stereopairs. The modelled estimates of tree height, crown size, and tree density were used to produce a decision-tree classification of forest successional stage at a landscape scale. The results obtained (72% accuracy), were limited in validation, but demonstrate potential for using the multi-scale methodology to provide spatial information for forestry policy objectives (ie., monitoring forest age structure).

Habitat selection and breeding success in a forest-nesting Alcid, the marbled murrelet, in two landscapes with different degrees of forest fragmentation

We studied habitat selection and breeding success in marked populations of a protected seabird (family Alcidae), the marbled murrelet (Brachyramphus marmoratus), in a relatively intact and a heavily logged old-growth forest landscape in south-western Canada. Murrelets used old-growth fragments either proportionately to their size frequency distribution (intact) or they tended to nest in disproportionately smaller fragments (logged). Multiple regression modelling showed that murrelet distribution could be explained by proximity of nests to landscape features producing biotic and abiotic edge effects. Streams, steeper slopes and lower elevations were selected in both landscapes, probably due to good nesting habitat conditions and easier access to nest sites. In the logged landscape, the murrelets nested closer to recent clearcuts than would be expected. Proximity to the ocean was favoured in the intact area. The models of habitat selection had satisfactory discriminatory ability in both landscapes. Breeding success (probability of nest survival to the middle of the chick rearing period), inferred from nest attendance patterns by radio-tagged parents, was modelled in the logged landscape. Survivorship was greater in areas with recent clearcuts and lower in areas with much regrowth, i.e. it was positively correlated with recent habitat fragmentation. We conclude that marbled murrelets can successfully breed in old-growth forests fragmented by logging.

The effect of nutrient enrichment on growth, photosynthesis and hydraulic conductance of dwarf mangroves in Panama

1. Dwarf stands of the mangrove Rhizophora mangle L. are extensive in the Caribbean. We fertilized dwarf trees in Almirante Bay, Bocas del Toro Province, north-eastern Panama with nitrogen (N) and phosphorus (P) to determine (1) if growth limitations are due to nutrient deficiency; and (2) what morphological and/or physiological factors underlie nutrient limitations to growth. 2. Shoot growth was 10-fold when fertilized with P and twofold with N fertilization, indicating that stunted growth of these mangroves is partially due to nutrient deficiency. 3. Growth enhancements caused by N or P enrichment could not be attributed to increases in photosynthesis on a leaf area basis, although photosynthetic nutrient-use efficiency was improved. The most dramatic effect was on stem hydraulic conductance, which was increased sixfold by P and 2.5-fold with N enrichment. Fertilization with P enhanced leaf and stem P concentrations and reduced C : N ratio, but did not alter leaf damage by herbivores. 4. Our findings indicate that addition of N and P significantly alter tree growth and internal nutrient dynamics of mangroves at Bocas del Toro, but also that the magnitude, pattern and mechanisms of change will be differentially affected by each nutrient.

Canopy carbon and oxygen isotope composition of 9-year-old hoop pine families in relation to seedling carbon isotope composition, growth, field growth performance, and canopy nitrogen concentration

Carbon isotope composition (delta C-13), oxygen isotope composition (delta O-18), and nitrogen concentration (N-mass) of branchlet tissue at two canopy positions were assessed for glasshouse seedlings and 9-year-old hoop pine (Araucaria cunninghamii Ait. ex D. Don) trees from 22 open-pollinated families grown in 5 blocks of a progeny test at a water-limited and nitrogen-deficient site in southeastern Queensland, Australia. Significant variations in canopy delta C-13, delta O-18, and N-mass existed among the 9-year-old hoop pine families, with a heritability estimate of 0.72 for branchlet delta C-13 from the upper inner canopy position. There was significant variation in canopy delta C-13 of glasshouse seedlings between canopy positions and among the families, with a heritability estimate of 0.66. The canopy delta C-13 was positively related to canopy N-mass only for the upper outer crown in the field (R = 0.62, p < 0.001). Phenotypic correlations existed between tree height and canopy delta C-13 (R = 0.37-0.41, p < 0.001). Strong correlations were found between family canopy delta C-13 at this site and those at a wetter site and between field canopy delta C-13 and glasshouse seedling delta C-13. The mechanisms of the variation in canopy delta C-13 are discussed in relation to canopy photosynthetic capacity as reflected in the N-mass and stomatal conductance as indexed by canopy delta O-18.

Insulin-like growth factor 2 cDNA cloning and ontogeny of gene expression in the liver of the marsupial brushtail possum (Trichosurus vulpecula)

The cDNA sequence for insulin-like growth factor 2 (IGF-2) was determined from the liver of the marsupial brushtail possum (Trichosurus vulpecula) using reverse transcription followed by polymerase chain reaction (RT-PCR) with gene-specific primers. The 359 bp of possum sequence encompassed the mature peptide, 27 bp of the signal peptide, and 125 bp of the E-peptide. Alignment of the deduced amino acid sequence with those from other species indicated that the mature peptide was 71 amino acids in length, 4 amino acids longer than most other mammals. At both the nucleotide and amino acid levels there was a high degree of sequence identity with IGF-2 from other mammalian and nonmammalian species. Amino acid identity ranged from 94.4% with a variant form of human IGF-2 to 80.3% with zebrafinch IGF-2. Northern analysis revealed that radiolabeled possum IGF-2, cDNA hybridized to multiple transcripts in the liver of both adult possums and 150-day-old pouch young and that the overall level of expression was greater in pouch young. Semiquantitative RT-PCR with total RNA from liver samples of pouch young aged 12 to 150 days postpartum and adults confirmed that IGF-2 gene expression was two to three times more abundant in pouch young than in adults but there was no significant change in the level of expression during pouch life. Unlike other mammalian species, in which there is a decline in levels of liver IGF-2 gene expression around the time of birth, levels in the marsupial brushtail possum remain elevated for at least 150 days after birth. This suggests that the decline in liver IGF-2 expression in marsupials and eutherians occurs at a similar stage of development and may reflect a role for this growth factor during the postnatal growth and development of the marsupial, (C) 2001 Academic Press.

Growth hormone (GH)/GH receptor expression and GH-mediated effects during early bovine embryogenesis

Pituitary growth hormone (GH) stimulates postnatal growth and metabolism. The role of CH and its receptor (GHR) during prenatal development, however, is still controversial. As shown by reverse transcription polymerase chain reaction (RT-PCR), bovine in vitro fertilization embryos synthesized the transcript of GHR from Day 2 of embryonic life onwards. Real time RT-PCR revealed that synthesis of GHR mRNA was increased 5.9-fold in 6-day-old embryos compared with 2-day-old embryos. Using in situ hybridization, the mRNA encoding GHR was predominantly localized to the inner cell mass of blastocysts. The GHR protein was first visualized 3 days after fertilization. GH-specific transcripts were first detected in embryos on Day 8 of in vitro culture. As shown by transmission electron microscopy, GH treatment resulted in elimination of glycogen storage in 6- to 8-day-old embryos and an increase in exocytosis of lipid vesicles. These results suggest that a functional GHR able to modulate carbohydrate and lipid metabolism is synthesized during preimplantation development of the bovine embryo and that this GHR may be subject to activation by embryonic GH after Day 8.

Denitrification, leaching and immobilisation of N-15-labelled nitrate in winter under windrowed harvesting residues in hoop pine plantations of 1-3 years old in subtropical Australia

A field study was carried out to investigate the impacts of windrowed harvesting residues on denitrification, immobilisation and leaching of N-15-labelled nitrate applied at 20 kg N ha(-1) to microplots in second-rotation hoop pine (Araucaria cunninghamii) plantations of 1-3 years old in southeast Queensland, Australia. The PVC microplots were 235 mm in diameter and 150 mm. long, and driven into the 100 mm soil. There were three replications of such microplots for each of the six treatments which were areas just under and between 1-, 2- and 3-year-old windrows of harvesting residues. Based on gaseous N losses estimated by the difference between the recoveries of bromide (Br) applied at 100 kg Br ha(-1) and N-15-labelled nitrate, denitrification was highest (23% based on N-15 loss) in the areas just under the 1-year-old windrows 25 days after a simulated 75 mm rainfall and following several natural rainfall events. There was no significant difference in N-15 losses (14-17%) among the other treatments. The N-15 immobilisation rate was highest for microplots in the areas between the 1-year-old windrows and generally higher for microplots in the areas just under the windrows (30-39%) than that (26-30%) between the windrows. Direct measurement of N-15 gas emissions (N-15(2) + (N2O)-N-15) confirmed that the highest denitrification rate occurred in the microplots under the 1-year-old windrows although the gaseous N-15 loss calculated by gas emission was only about one-quarter that estimated by the N-15 mass balance method. A significant, positive linear relationship (P < 0.05) existed between the gaseous N-15 losses measured by the two methods used. The research indicates that considerable mineral N could be lost via denitrification during the critical inter-rotation period and early phase of the second rotation. However, the impacts of windrowed harvesting residues on N losses via denitrification might only last for a period of about 2 years. Published by Elsevier Science B.V.

Growth of subtropical ECM fungi with different nitrogen sources using a new floating culture technique

Eight species of ectomycorrhizal (ECM) fungi in the genera Amanita. Gymnoboletus, Lactarius, and Russula were isolated from subtropical plant communities in eastem Australia. Two species were isolated from each of rainforest, Nothofagus forest, Eucalyptus forest, and Eucalyptus dominated wallum (heath) forest. These communities differ strongly in their soluble soil nitrogen (N) composition. The ability of the fungi to use inorganic (nitrate, ammonium) and organic (amide, peptide, protein) nitrogen sources was determined. As the fungi did not grow in liquid culture, a 'floating culture' technique was devised that allows hyphal growth on a screen floating on liquid medium. With some exceptions, fungal biomass production in floating culture closely reflected fungal growth on solid media assessed by total colony glucosamine content. Most isolates grown in floating culture had similar glucosamine concentrations on all N sources, with isolate specific concentrations ranging from 6 to 12 mug glucosamine g(-1) DW. However, Russula spp. had up to 1.7-fold higher glucosamine concentrations when growing with glutamine or ammonium compared to nitrate, glutathione or protein. Floating cultures supplied with 0.5, 1.5. 4.5, or 10 mm N mostly produced greatest biomass with 4.5 mM N. In vitro nitrate reductase activity (NRA) ranged from very low (0.03 mumol NO2- g(-1) fw h(-1)) in Russula sp. (wallum) to high (2.16 mumol NO2- g(-1) fw h(-1)) in Gymnoboletus sp. (rainforest) and mirrored the fungi's ability to use nitrate as a N source. All Russula spp. (wallum, Nothofagus and Eucalyptus forests), Lactarills sp, (rainforest) and.4manita sp. (wallum) utilized ammonium and glutamine but had little ability to use other N sources. In contrast,Amanita species (Nothofagus and Eucalyptus forests) grew on all N sources but produced most biomass with ammonium and glutamine. Only Gymnoboletus sp. (rainforest) showed similar growth with nitrate and ammonium as N sources. Fungal N source use was not associated with taxonomic groups, but is discussed in the context of soil N sources in the different habitats.

Mapping glider habitat in dry eucalypt forests for Montreal Process indicator 1.1: Fragmentation of forest types

Accurate habitat mapping is critical to landscape ecological studies such as required for developing and testing Montreal Process indicator 1.1e, fragmentation of forest types. This task poses a major challenge to remote sensing, especially in mixedspecies, variable-age forests such as dry eucalypt forests of subtropical eastern Australia. In this paper, we apply an innovative approach that uses a small section of one-metre resolution airborne data to calibrate a moderate spatial resolution model (30 m resolution; scale 1:50 000) based on Landsat Thematic Mapper data to estimate canopy structural properties in St Marys State Forest, near Maryborough, south-eastern Queensland. The approach applies an image-processing model that assumes each image pixel is significantly larger than individual tree crowns and gaps to estimate crown-cover percentage, stem density and mean crown diameter. These parameters were classified into three discrete habitat classes to match the ecology of four exudivorous arboreal species (yellowbellied glider Petaurus australis, sugar glider P. breviceps, squirrel glider P. norfolcensis , and feathertail glider Acrobates pygmaeus), and one folivorous arboreal marsupial, the greater glider Petauroides volans. These species were targeted due to the known ecological preference for old trees with hollows, and differences in their home range requirements. The overall mapping accuracy, visually assessed against transects (n = 93) interpreted from a digital orthophoto and validated in the field, was 79% (KHAT statistic = 0.72). The KHAT statistic serves as an indicator of the extent that the percentage correct values of the error matrix are due to ‘true’ agreement verses ‘chance’ agreement. This means that we are able to reliably report on the effect of habitat loss on target species, especially those with a large home range size (e.g. yellow-bellied glider). However, the classified habitat map failed to accurately capture the spatial patterning (e.g. patch size and shape) of stands with a trace or sub-dominance of senescent trees. This outcome makes the reporting of the effects of habitat fragmentation more problematic, especially for species with a small home range size (e.g. feathertail glider). With further model refinement and validation, however, this moderateresolution approach offers an important, cost eff e c t i v e advancement in mapping the age of dry eucalypt forests in the region.

Intraindividual allometric development of aerobic power in 8- to 16-year-old boys

Purpose: The aims of this study are two-fold: first, to analyze intraindividual allometric development of aerobic power of 73 boys followed at annual intervals from 8 to 16 yr, and second, to relate scaled aerobic power with level of habitual physical activity and biological maturity status. Methods: Peak (V) over dot O-2 (treadmill), height, and body mass were measured. Biological maturity was based on age at peak height velocity (PHV) and level of physical activity was based on five assessments between 11 and 15 yr and at 17 yr. Interindividual and intraindividual allometric coefficients were calculated. Multilevel modeling was applied to verify if maturity status and activity explain a significant proportion of peak (V) over dot O-2 after controlling for other explanatory characteristics. Results: At most age levels, interindividual allometry coefficients for body mass exceed k = 0.750. Intraindividual coefficients of peak (V) over dot O-2 by body mass vary widely and range from k' = 0,555 to k' = 1,178. Late maturing boys have smaller k' coefficients than early maturing boys. Conclusion: Peak (V) over dot O-2 is largely explained by body mass, but activity level and its interaction with maturity status contribute independently to peak (V) over dot O-2 even after adjusting for body mass.

Control of ovulation with a GnRH agonist after superstimulation of follicular growth in buffalo: fertilization and embryo recovery

The potential to use a GnRH agonist bioimplant and injection of exogenous LH to control the time of ovulation in a multiple ovulation and embryo transfer (MOET) protocol was examined in buffalo. Mixed-parity buffalo (Bubalus bubalis; 4-15-year-old; 529 13 kg LW) were randomly assigned to one of five groups (n = 6): Group 1, conventional MOET protocol; Group 2, conventional MOET with 12 It delay in injection of PGF(2alpha); Group 3, implanted with GnRH agonist to block the pre-ovulatory surge release of LH; Group 4, implanted with GnRH agonist and injected with exogenous LH (Lutropin(R), 25 mg) 24 h after 4 days of superstimulation with FSH; Group 5, implanted with GnRH agonist and injected with LH 36 h after superstimulation with FSH. Ovarian follicular growth in all buffaloes was stimulated by treatment with FSH (Folltropin-V(R), 200 mg) administered over 4 days, and was monitored by ovarian ultrasonography. At the time of estrus, the number of follicles greater than or equal to8 mm. was greater (P < 0.05) for buffaloes in Group 2 (12.8) than for buffaloes in Groups 1 (8.5), 3 (7.3), 4 (6.1) and 5 (6.8), which did not differ. All buffaloes were mated by AI after spontaneous (Groups 1-3) or induced (Groups 4 and 5) ovulation. The respective number of buffalo that ovulated, number of corpora lutea, ovulation rate (%), and embryos + oocytes recovered were: Group 1 (2, 1.8 +/- 1.6, 18.0 +/- 13.6, 0.2 +/- 0.2); Group 2 (4, 6.1 +/- 2.9, 40.5 +/- 17.5, 3.7 +/- 2.1); Group 3 (0, 0, 0, 0); Group 4 (6, 4.3 +/- 1.2, 69.3 +/- 14.2, 2.0 +/- 0.9); and Group 5 (1, 2.5 +/- 2.5, 15.5 +/- 15.5, 2.1 +/- 2.1). All buffaloes in Group 4 ovulated after injection of LH and had a relatively high ovulation rate (69%) and embryo recovery (46%). It has been shown that the GnRH agonist-LH protocol can be used to improve the efficiency of MOET in buffalo. (C) 2002 Elsevier Science Inc. All rights reserved.

Testing landscape metrics as indicators of habitat loss and fragmentation in continuous eucalypt forests (Queensland, Australia)

Landscape metrics are widely applied in landscape ecology to quantify landscape structure. However, many are poorly tested and require rigorous validation if they are to serve as reliable indicators of habitat loss and fragmentation, such as Montreal Process Indicator 1.1e. We apply a landscape ecology theory, supported by exploratory and confirmatory statistical techniques, to empirically test landscape metrics for reporting Montreal Process Indicator 1.1e in continuous dry eucalypt forests of sub-tropical Queensland, Australia. Target biota examined included: the Yellow-bellied Glider (Petaurus australis); the diversity of nectar and sap feeding glider species including P. australis, the Sugar Glider P. breviceps, the Squirrel Glider P. norfolcensis, and the Feathertail Glider Acrobates pygmaeus; six diurnal forest birds species; total diurnal bird species diversity; and the density of nectar-feeding diurnal bird species. Two scales of influence were considered: the stand-scale (2 ha), and a series of radial landscape extents (500 m - 2 km; 78 - 1250 ha) surrounding each fauna transect. For all biota, stand-scale structural and compositional attributes were found to be more influential than landscape metrics. For the Yellow-bellied Glider, the proportion of trace habitats with a residual element of old spotted-gum/ironbark eucalypt trees was a significant landscape metric at the 2 km landscape extent. This is a measure of habitat loss rather than habitat fragmentation. For the diversity of nectar and sap feeding glider species, the proportion of trace habitats with a high coefficient of variation in patch size at the 750 m extent was a significant landscape metric. None of the landscape metrics tested was important for diurnal forest birds. We conclude that no single landscape metric adequately captures the response of the region's forest biota per se. This poses a major challenge to regional reporting of Montreal Process Indicator 1.1e, fragmentation of forest types.

Growth hormone receptor and insulin-like growth factor-I immunoreactivity in osteoclast-like cells during tooth eruption in the toothless (Osteopetrotic) rat following treatment with colony-stimulating factor-1

In the toothless (tl/tl) osteopetrotic rat, teeth form but fail to erupt. Treatment of tl/tl rats with colony-stimulating factor-1 (CSF-1) activates bone resorption by osteoclasts, permits tooth eruption, and up-regulates the immunoreactivity of bone marrow mononuclear cells to growth hormone receptor (GHr) and insulin-like growth factor (IGF)-I. This study examined the distribution of tartrate-resistant acid phosphatase (TRAP) and immunoreactivity for GHr and IGF-I in osteoclast-like cells located on the alveolar bone margin, adjacent to the lower first molar crown, in 14-day-old normal and tl/tl rats, following treatment with CSF-1. Osteoclast-like cells demonstrated a positive reaction for TRAP, GHr, and IGF-I in all groups. However, in tl/tl tissue, osteoclast-like cells were generally negative for GHr. There was no significant difference in the total number of TRAP, GHr, and IGF-I-positive osteoclast-like cells on the adjacent bone margin in normal, normal treated with CSF-1, and tl/tl rats. CSF-1 treatment of the tl/tl rat significantly increased the total number of osteoclast-like cells, which were positive for TRAP (p < 0.001), GHr (p < 0.05) and IGF-I (P < 0.01).

Genetic variation in foliar carbon isotope composition in relation to tree growth and foliar nitrogen concentration in clones of the F-1 hybrid between slash pine and Caribbean pine

The objectives of this study were: (1) to quantify the genetic variation in foliar carbon isotope composition (delta(13)C) of 122 clones of ca. 4-year-old F-1 hybrids between slash pine (Pinus elliottii Engelm var. elliottii) and Caribbean pine (Pinus caribaea var. hondurensis Barr.,et Golf.) grown at two field experimental sites with different water and nitrogen availability in southeast Queensland, Australia, in relation to tree growth and foliar nitrogen concentration (N-mass); and (2) to assess the potential of using delta(13)C measurements, in the foliage materials collected from the clone hedges at nursery and the 4-year-old tree canopies in the field, as an indirect index of tree water use efficiency for selecting elite F-1 hybrid pine clones with improved tree growth. There were significant differences in foliar delta(13)C between the nursery hedges and the 4-year-old tree canopies in the field, between the summer and winter seasons, between the two experimental sites, and between the upper outer and lower outer canopy positions sampled. This indicates that delta(13)C measurements in the foliage materials are significantly influenced by the sampling techniques and environmental conditions. Significant differences in foliar delta(13)C, at the upper outer canopy in both field experiments in summer and winter, were detected between the clones, and between the female parents of the clones. Clone means of tree height at age ca. 3 years were positively related to those of the upper outer canopy delta(13)C at both experimental sites in winter, but only for the wetter site in summer. There were positive, linear relationships between clone means of canopy delta(13)C and those of canopy N-mass, indicating that canopy photosynthetic capacity might be an important factor regulating the clonal variation in canopy delta(13)C. Significant correlations were found between clone means of canopy delta(13)C at both experimental sites in summer and winter, and between those at the upper outer and lower outer canopy positions. Mean clone delta(13)C for the nursery hedges was only positively related to mean clone stem diameter at 1.3 m height at age 3 years on the wetter site. The clone by site interaction for foliar delta(13)C at the upper outer canopy was significant only in summer. Overall, the relatively high genetic variance components for foliar delta(13)C and significant, positive correlations between clone means of foliar delta(13)C and tree growth have highlighted the potential of using foliar delta(13)C measurements for assisting in selection of the elite F-1 hybrid pine clones with improved tree growth. (C) 2002 Elsevier Science B.V. All rights reserved.

Composition and quality of harvest residues and soil organic matter under windrow residue management in young hoop pine plantations as revealed by solid-state C-13 NMR spectroscopy

Solid-state C-13 nuclear magnetic resonance (NMR) with cross-polarisation (CP) and magic-angle-spinning (MAS) was used to: (a) examine the changes in carbon (C) composition of windrowed harvest residues during the first 3 years of hoop pine plantations in subtropical Australia; (b) assess the impacts of windrowed harvest residues on soil organic matter (SOM) composition and quality in the 0-10 cm soil layer. Harvest residues were collected from 0-, 1-, 2- and 3-year-old windrows of ca. 2.5 m width (15 m apart for 0-, 1- and 2-year-old sites and 10 m apart for 3-year-old site). Soils from the 0 to 10 cm soil layer were collected from the 1-, 2- and 3-year-old sites. The 13C NMR spectra of the harvest residues indicated the presence of lignin in the hoop pine wood, foliage and newly incorporated organic matter (NIOM). Condensed tannin structures were found in the decay-resistant bark, small wood and foliage, but were absent in other residue components and SOM. The NMR spectra of small wood samples contained condensed tannin structures because the outer layer of bark was not removed. NIOM showed a shift from foliage-like structures (celluloses) to lignin-type structures, indicating an incorporation of woody residues from the decomposing harvest residues. Suberins were also present in the small wood, foliage and bark. The 13C CP NMR spectra of SOM indicated that in areas where windrows were present, SOM did not show compositional changes. However, an increase in SOM quality under the windrows in the second year after their formation as characterised by the alkyl C/O-alkyl C (A/O-A) ratio was mainly due to inputs from the decomposition of the labile, readily available components of the windrowed harvest residues. (C) 2002 Published by Elsevier Science B.V.