In vitro physical mutagenesis of giant reed (Arundo donax L.)

Giant reed (Arundo donax L.) is a C₃ perennial, warm‐season, rhizomatous grass of emerging interest for bioenergy and biomass derivatives production, and for phytoremediation. It only propagates vegetatively and very little genetic variation is found among ecotypes, basically precluding breeding efforts. With the objective to increase the genetic variation in this species, we developed and applied a mutagenesis protocol based on γ‐irradiation of in vitro cell cultures from which regenerants were obtained. Based on a radiosensitivity test, the irradiation dose reducing to 50% the number of regenerants per callus (RD₅₀) was estimated at 35 Gy. A large mutagenic experiment was carried out by irradiating a total of 3120 calli with approx. 1×, 1.5× and 2× RD₅₀. A total of 1004 regenerants from irradiated calli were hardened in pots and transplanted to the field. Initial phenotypic characterization of the collection showed correlated responses of biomass‐related quantitative traits to irradiation doses. Approx. 10% of field‐grown clones showed remarkable morphological aberrations including dwarfism, altered tillering, abnormal inflorescence, leaf variegation and others, which were tested for stability over generations. Clone lethality reached 0.4%. Our results show for the first time that physical mutagenesis can efficiently induce new genetic and phenotypic variation of agronomic and prospective industrial value in giant reed. The methodology and the plant materials described here may contribute to the domestication and the genetic improvement of this important biomass species.     

Screening of giant reed (Arundo donax L.) ecotypes for biomass production under salt stress

Given its high biomass and plasticity, Arundo donax L. is a promising ligno-cellulosic crop for cultivation in marginal lands in temperate climates. In order to test for adaptation to salinity, growth parameters of several A. donax clones were evaluated under two salt regimes in hydroponics. Mild NaCl stress (50 mM NaCl, 5.6 mS cm^?1 EC, for 10 days) failed to discriminate between ecotypes, while a more severe NaCl treatment (150 mM, 18.8 mS cm^?1 EC, for 21 days) enabled the identification of ecotypes maintaining plant growth under high salinity. Among several biometric parameters, 4th leaf width, and shoot and root DW consistently highlighted differences between ecotypes. Gas-exchange parameters also responded to severe NaCl treatment, while the photosystem efficiency was good, regardless of treatment. The results confirm that A. donax can be considered moderately tolerant to NaCl stress, with variation between ecotypes. Our screening protocol identified ecotypes with higher biomass production under severe NaCl treatment and can be useful for preliminary evaluation of NaCl tolerant clones for increasing productivity under salinity. The detected inter-ecotype variability could also be investigated to identify suitable clones for different environments.     

Molecular mapping and location of QTLs for drought-resistance traits in <Emphasis Type="Italic">indica</Emphasis> rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.) lines adapted to target environments

Drought is a major limitation for rice production in rainfed ecosystems. Identifying quantitative trait loci (QTLs) linked to drought resistance provides opportunity to breed high yielding rice varieties suitable for drought-prone areas. Although considerable efforts were made in mapping QTLs associated with drought-resistance traits in rice, most of the studies involved indica × japonica crosses and hence, the drought-resistance alleles were contributed mostly by japonica ecotypes. It is desirable to look for genetic variation within indica ecotypes adapted to target environment (TE) as the alleles from japonica ecotype may not be expressed under lowland conditions. A subset of 250 recombinant inbred lines (RILs) of F₈ generation derived from two indica rice lines (IR20 and Nootripathu) with contrasting drought-resistance traits were used to map the QTLs for morpho-physiological and plant production traits under drought stress in the field in TE. A genetic linkage map was constructed using 101 polymorphic PCR-based markers distributed over the 12 chromosomes covering a total length of 1,529 cM in 17 linkage groups with an average distance of 15.1 cM. Composite interval mapping analysis identified 22 QTLs, which individually explained 4.8–32.2% of the phenotypic variation. Consistent QTLs for drought-resistance traits were detected using locally adapted indica ecotypes, which may be useful for rainfed rice improvement.     

Growth and nutrients accumulation potentials of giant reed (Arundo donax L.) in different habitats in Egypt

Arundo donax L. has a high biomass production and a tendency toward community dominance in many habitats and thereby a tolerance to a wide range of environmental conditions. Therefore, the present study investigated the potentiality of A. donax to accumulate nutrients and trace metals in its biomass. Six main habitats (Nile Bank, Drain Bank, Canal Bank, Field Edges, Railways and Roadsides) were recognized. At each habitat, six quadrats (each 1 m²), distributed equally in two sites, were selected for growth measurements (e.g., density, shoot height, diameter, leaf area and biomass), plant and soil analyses. Plants from Nile, Canal and Drain Banks had the highest values of most growth measurements, while those from Railways and Roadsides had the lowest. Canal Bank plants accumulated the highest concentrations of P, Cu and Pb in their leaves; Zn in the stem; and Mg, Cd and Fe in the rhizome. The bioaccumulation factor (BF) of A. donax, for Cd, Fe, Mn and Zn, was greater than 1, while the translocation factor (TF) of most trace metals was less than unity in most habitats. In conclusion, A. donax showed morphological plasticity in response to habitat heterogeneity, and its growth was most vigorous in the riparian habitats. The high BF, as well as the significant positive correlations between trace metals, especially Cd, in soil and plant, renders A. donax a powerful phytoremediator.     

Construction and comparison of three reference‐quality genome assemblies for soybean

We report reference‐quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single‐nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan‐gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40–42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.     

Effects of ecotypes and morphotypes in feedstock composition of switchgrass (Panicum virgatum L.)

Switchgrass (Panicum virgatum L.) is a C4 grass with high biomass yield potential and is now a model species for the Bioenergy Feedstock Development Program. Two distinct ecotypes (e.g., upland and lowland) and a range of plant morphotypes (e.g., leafy and stemmy) have been observed in switchgrass. The objective of this study was to determine the influence of ecotype and morphotype on biomass feedstock quality. Leaf and stem tissues of leafy and stemmy morphotypes from both lowland and upland ecotypes were analyzed for key feedstock traits. The leaf : stem ratio of leafy morphotype was more than 40% higher than the stemmy morphotype in both upland and lowland ecotypes. Therefore, the stemmy morphotype has significant advantages over leafy morphotype during harvesting, storage, transportation and finally the feedstock quality. Remarkable differences in feedstock quality and mineral composition were observed in switchgrass genotypes with distinct ecotypic origins and variable plant morphotypes. Lignin, hemicelluloses and cellulose concentrations were higher in stems than in the leaves, while ash content was notably high in leaves. A higher concentration of potassium was found in the stems compared to the leaves. In contrast, calcium was higher and magnesium was generally higher in the leaves compared to stems. The upland genotypes demonstrated considerably higher lignin (14.4%) compared with lowland genotypes (12.4%), while hemicellulose was higher in lowland compared with upland. The stemmy type demonstrated slightly higher lignin compared with leafy types (P < 0.1). Differences between the ecotypes and morphotypes for key quality traits demonstrated the potential for improving feedstock composition of switchgrass through selection in breeding programs.     

Genotypic variation of cell wall composition and its conversion efficiency in Miscanthus sinensis,a potential biomass feedstock crop in China

Plant cell walls are composed of cellulose microfibrils embedded in a cross‐linked‐net of matrix polysaccharides and co‐polymerized with lignin. The study presented the genotypic variations of cell wall composition, biohydrogen production, and lignocellulose degradation ratio in a collection of 102 Miscanthus sinensis (M. Sinensis, hereafter) accessions collected from a wide geographical range in China. Significant variations were observed for the determined traits, cellulose content, hemicellulose content, cellulose and hemicellulose degradation efficiency, and biohydrogen yield. The cellulose, hemicellulose, and lignin contents ranged from 30.20–44.25, 28.97–42.65, and 6.96–20.75%, respectively. The degradation ratio of cellulose and hemicellulose varied from 2.08% to 37.87% and from 14.71% to 52.50%, respectively. The feedstock was fermented to produce biohydrogen, and the production varied from 14.59 to 40.66 ml per gram of Miscanthus biomass. The expression profile of three cellulose synthase (MsCesA) genes was initially established to indicate the genotypic difference among the M. sinensis accessions. Pearson's correlations were conducted to reveal the perplexing relationship between the tested traits, biohydrogen yield, cell wall composition and their degradation efficiency. In addition, the relationship pattern, between the test traits and the geographic factors corresponding with the original place, was investigated. The result showed that the significant variation among the M. sinensis genotypes is the result of natural selection in different environments of their original habitats. Improvement in cell wall composition and structure and enhancement of lignocellulose degradation ratio could significantly increase sustainable bioenergy production.     

A salt on the bioenergy and biological invasions debate: salinity tolerance of the invasive biomass feedstock Arundo donax

Arundo donax L., commonly known as giant reed, is promising biomass feedstock that is also a notorious invasive plant in freshwater ecosystems around the world. Heretofore, the salt tolerance of A. donax had not been quantified even though anecdotal evidence suggests halophytic qualities. To test whole-plant and leaf level responses, we established a pot experiment on 80 scions propagated from an A. donax population that has naturalized on the shore of the San Francisco Bay Estuary. To quantify growth and physiological responses to salinity (NaCl), A. donax scions were divided into eight treatments and grown for 60 days across a range of salinities (0–42 dS m⁻¹). Classic growth analysis showed >80% reduction in overall growth at the highest salinities. Yet, there was zero mortality indicating that A. donax is able to tolerate high levels of salt. Declining photosynthesis rates were strongly correlated (R² > 0.97) with decreasing stomatal conductance, which was in turn closely related to increasing salinity. Leaf gas exchange revealed that stomata and leaf limitations of carbon dioxide were three times greater at high salinities. Nonetheless, even when salinities were 38–42 dS m⁻¹ A. donax was able to maintain assimilation rates 7–12 μmol m⁻² s⁻¹. Further, by maintaining 50% relative growth at salinities ~12 dS m⁻¹ A. donax can now be classified as ‘moderately salt tolerant’. A. donax leaf gas exchange and whole-plant salt tolerance are greater than many important food crops (i.e. maize, rice), the bioenergy feedstock Miscanthus × giganteus, as well as some uncultivated plant species (i.e. Populus and Salix) that are indigenous in regions A. donax currently invades. The results of this study have implications for both agronomists wishing to expand A. donax to fields dominated by saline soils, and for others who are concerned about the spread of A. donax with altered stream hydrology or sea-level rise.     

Quantifying plant biomass and seed production in camelina (Camelina sativa (L.) Crantz) across a large range of plant densities: Modelling approaches

To date, there has been little agreement on supporting the hypothesis that how some key vegetative traits of camelina (Camelina sativa (L.) Crantz var. ‘Soheil’) are dependent on plant biomass. Therefore, the main aim of this investigation was to quantify the relationship between the size of camelina plants and seed production across a broad-range of plant densities through modelling approaches. To make a wide range of plant densities, a fan design was used in eight replicates in an experimental field at Sari Agricultural Sciences and Natural Resources University, Iran. To quantify the relation between plant density and other plant traits, a regression analysis was carried out and the coefficient of determination (R²) was considered to evaluate the goodness of fit model. A power model (y = ax^b) could describe well the relationship between plant density (ranged 113–2905 plants m⁻²) and plant biomass, seed production, number of seeds per plant, stem diameter, and siliques number, with the coefficient of determination (R²) values of 0.85, 0.87, 0.65, 0.64, and 0.90, respectively. The harvest indexes were 13.8%–26.9%, depending on plant density. Seed production per plant was positively correlated to the siliques number (r = 0.85), the branch number (r = 0.80), and the seed number (r = 0.99) which could be key components of camelina seed production per plant. Furthermore, no significant correlation was found among plant height, thousand-seed weight, and harvest index with seed production per plant. In conclusion, plant biomass could be considered an important trait to predict plant growth models of camelina. Also, a lower plant density of camelina can be compensated by a greater number of siliques, branches and seeds per plant.     

Pattern and control of biomass allocation across global forest ecosystems

The underground part of a tree is an important carbon sink in forest ecosystems. Understanding biomass allocation between the below‐ and aboveground parts (root:shoot ratios) is necessary for estimation of the underground biomass and carbon pool. Nevertheless, large‐scale biomass allocation patterns and their control mechanisms are not well identified. In this study, a large database of global forests at the community level was compiled to investigate the root:shoot ratios and their responses to environmental factors. The results indicated that both the aboveground biomass (AGB) and belowground biomass (BGB) of the forests in China (medians 73.0 Mg/ha and 17.0 Mg/ha, respectively) were lower than those worldwide (medians 120.3 Mg/ha and 27.7 Mg/ha, respectively). The root:shoot ratios of the forests in China (median = 0.23), however, were not significantly different from other forests worldwide (median = 0.24). In general, the allocation of biomass between the belowground and aboveground parts was determined mainly by the inherent allometry of the plant but also by environmental factors. In this study, most correlations between root:shoot ratios and environmental factors (development parameter, climate, altitude, and soil) were weak but significant (p < .01). The allometric model agreed with the trends observed in this study and effectively estimated BGB based on AGB across the entire database.     

How does plant population density affect the biomass of Ravenna grass?

Ravenna grass, Tripidium ravennae (L.) H. Scholz, is known to produce an abundance of biomass, but how plant density affects its biomass potential remains unknown. The objectives were to determine the effects of plant density on biomass yield; plant growth traits; biomass?carbon, nitrogen, and ash concentrations; heating value; nitrogen removal; and sucrose concentration in leaves and culms. The treatments consisted of five plant densities (1,250; 2,500; 5,000; 10,000; and 20,000 plants per hectare) in a randomized complete block design with four blocks. Plots were nonirrigated, unfertilized, and harvested once during the dormant season each year. Data were collected from 2015?2019. Dependent variables that varied with plant population density (p < .05) were biomass yield, number of reproductive culms per plant, reproductive culm diameter, reproductive culm sucrose concentration, and nitrogen removal with biomass. Biomass yield ranged from 5.6 to 16.3 Mg/ha for plant densities of 1,250–20,000 plants per hectare, respectively. Combined over years, nonlinear regression of the data showed the equation for biomass yield to plateau at 16.2 Mg/ha at a plant density of 10,640 plants per hectare. As plant density increased, the number of reproductive culms per plant, culm diameter, and culm sucrose concentration significantly decreased. At 1,250 plants per hectare, the number of reproductive culms per plant, culm diameter, and culm sucrose averaged 70, 10.2 mm, and 63.2 g/kg, respectively. Nitrogen removed with biomass significantly increased as biomass yield increased with plant density. At a density of 10,000 and 20,000 plants per hectare, the amount of nitrogen removed annually in the harvested biomass averaged 88 kg/ha. The data suggest that 10,000 plants per hectare would produce the greatest annual biomass yields; however, research is needed to determine the nutrient requirement for Ravenna grass to sustain biomass production at that density.     

Genome‐wide association and genomic prediction for biomass yield in a genetically diverse Miscanthus sinensis germplasm panel phenotyped at five locations in Asia and North America

To improve the efficiency of breeding of Miscanthus for biomass yield, there is a need to develop genomics‐assisted selection for this long‐lived perennial crop by relating genotype to phenotype and breeding value across a broad range of environments. We present the first genome‐wide association (GWA) and genomic prediction study of Miscanthus that utilizes multilocation phenotypic data. A panel of 568 Miscanthus sinensis accessions was genotyped with 46,177 single nucleotide polymorphisms (SNPs) and evaluated at one subtropical and five temperate locations over 3 years for biomass yield and 14 yield‐component traits. GWA and genomic prediction were performed separately for different years of data in order to assess reproducibility. The analyses were also performed for individual field trial locations, as well as combined phenotypic data across groups of locations. GWA analyses identified 27 significant SNPs for yield, and a total of 504 associations across 298 unique SNPs across all traits, sites, and years. For yield, the greatest number of significant SNPs was identified by combining phenotypic data across all six locations. For some of the other yield‐component traits, greater numbers of significant SNPs were obtained from single site data, although the number of significant SNPs varied greatly from site to site. Candidate genes were identified. Accounting for population structure, genomic prediction accuracies for biomass yield ranged from 0.31 to 0.35 across five northern sites and from 0.13 to 0.18 for the subtropical location, depending on the estimation method. Genomic prediction accuracies of all traits were similar for single‐location and multilocation data, suggesting that genomic selection will be useful for breeding broadly adapted M. sinensis as well as M. sinensis optimized for specific climates. All of our data, including DNA sequences flanking each SNP, are publicly available. By facilitating genomic selection in M. sinensis and Miscanthus × giganteus, our results will accelerate the breeding of these species for biomass in diverse environments.     

Diversity and population structure of northern switchgrass as revealed through exome capture sequencing

Panicum virgatum L. (switchgrass) is a polyploid, perennial grass species that is native to North America, and is being developed as a future biofuel feedstock crop. Switchgrass is present primarily in two ecotypes: a northern upland ecotype, composed of tetraploid and octoploid accessions, and a southern lowland ecotype, composed of primarily tetraploid accessions. We employed high‐coverage exome capture sequencing (~2.4 Tb) to genotype 537 individuals from 45 upland and 21 lowland populations. From these data, we identified ~27 million single‐nucleotide polymorphisms (SNPs), of which 1 590 653 high‐confidence SNPs were used in downstream analyses of diversity within and between the populations. From the 66 populations, we identified five primary population groups within the upland and lowland ecotypes, a result that was further supported through genetic distance analysis. We identified conserved, ecotype‐restricted, non‐synonymous SNPs that are predicted to affect the protein function of CONSTANS (CO) and EARLY HEADING DATE 1 (EHD1), key genes involved in flowering, which may contribute to the phenotypic differences between the two ecotypes. We also identified, relative to the near‐reference Kanlow population, 17 228 genes present in more copies than in the reference genome (up‐CNVs), 112 630 genes present in fewer copies than in the reference genome (down‐CNVs) and 14 430 presence/absence variants (PAVs), affecting a total of 9979 genes, including two upland‐specific CNV clusters. In total, 45 719 genes were affected by an SNP, CNV, or PAV across the panel, providing a firm foundation to identify functional variation associated with phenotypic traits of interest for biofuel feedstock production.     

Biomass yields,cytogenetics, fertility,and compositional analyses of novel bioenergy grass hybrids (Tripidium spp.)

High biomass yields have been documented for Tripidium spp. (Erianthus spp., Saccharum spp.), but targeted breeding for bioenergy applications has been limited. Advanced, interspecific hybrids between Tripidium ravennae and T. arundinaceum were planted in replicated field plots in 2016. Comparative feedstock evaluations examined biomass yields, cytogenetics, plant fertility, and compositional analyses relative to Miscanthus × giganteus. Dry biomass yields varied as a function of year and accession and increased each year ranging from 3.4 to 10.6, 8.6 to 37.3, and 23.7 to 60.6 Mg/ha for Tripidium hybrids compared to 2.3, 16.2 and 27.9 Mg/ha for M. × giganteus in 2016, 2017, and 2018, respectively. Cytology and cytometry confirmed that Tripidium hybrids were tetraploid with 2n = 4x = 40 (2C genome size = 5.06 pg) and intermediate between T. ravennae with 2n = 2x = 20 (2C genome size = 2.55 pg) and T. arundinaceum with 2n = 6x = 60 (2C genome size = 7.61 pg). Plant fertility characteristics varied considerably with some accessions producing no viable seeds or fewer than that observed for M. × giganteus. Accessions varied significantly for flowering culm number and height and dates of peak anthesis ranging from 14 September to 2 October. Variations in yield and compositional analyses contributed to variations in theoretical ethanol yields ranging from 10,181 to 27,546 L/ha for Tripidium accessions compared to 13,095 L/ha for M. × giganteus. Relative feed value (RFV) indices for winter‐harvested Tripidium accessions varied from 52.8 to 60.0 compared to M. × giganteus with 45.4. RFV for summer‐harvested Tripidium accessions varied from 71.6 to 80.5 compared to M. × giganteus with 61.0. These initial findings for Tripidium hybrids, including high biomass yields, cold hardiness, and desirable traits for multiple markets (e.g., forage, bioenergy, bioproducts), are promising and warrant further development of Tripidium as a temperate bioenergy feedstock.     

Ecological correlates of invasion by <Emphasis Type="Italic">Arundo donax</Emphasis> in three southern California riparian habitats

Arundo donax L. (Poaceae) is an aggressive invader in California’s riparian habitats. Field experiments were conducted to examine invader and site attributes important in early invasion. One hundred A. donax rhizomes were planted along five transects into each of three southern California riparian habitats. Pre-planting rhizome weight was recorded, along with site variables including percent bare ground, litter depth, PAR, soil moisture, soil temperature, incidence of herbivory, native canopy cover, and plant community richness and diversity. A. donax shoot emergence, survival time, and shoot height were recorded for approximately 10 months. The experiment was repeated over three years in different locations within each site. When years and sites were pooled to reveal large-scale patterns, A. donax performance was explained by rhizome weight, soil moisture, bare ground, soil temperature, and herbivory. When each site was considered singly, A. donax was positively correlated with different variables in each location. Species richness was correlated with A. donax performance in only one site. Our results indicate that A. donax establishment in riparian habitats is promoted by both vegetative reproduction and favorable abiotic environmental factors and relatively unaffected by the composition of the native community. The positive response of A. donax to disturbance (bare ground) and high resource availability (soil moisture), combined with a competitive perennial habit suggest that this species takes advantage of a competitive-ruderal life history. The ability of A. donax to respond to different conditions in each site combined with low genetic and phenotypic variation seen in other studies also suggests that a high degree of environmental tolerance contributes to invasion success.     

Host‐associated genomic differentiation in congeneric butterflies: now you see it,now you do not

Ecotypic variation among populations may become associated with widespread genomic differentiation, but theory predicts that this should happen only under particular conditions of gene flow, selection and population size. In closely related species, we might expect the strength of host‐associated genomic differentiation (HAD) to be correlated with the degree of phenotypic differentiation in host‐adaptive traits. Using microsatellite and Amplified Fragment Length Polymorphism (AFLP) markers, and controlling for isolation by distance between populations, we sought HAD in two congeneric species of butterflies with different degrees of host plant specialization. Prior work on Euphydryas editha had shown strong interpopulation differentiation in host‐adapted traits, resulting in incipient reproductive isolation among host‐associated ecotypes. We show here that Euphydryas aurinia had much weaker host‐associated phenotypic differentiation. Contrary to our expectations, we detected HAD in Euphydryas aurinia, but not in E. editha. Even within an E. aurinia population that fed on both hosts, we found weak but significant sympatric HAD that persisted in samples taken 9 years apart. The finding of significantly stronger HAD in the system with less phenotypic differentiation may seem paradoxical. Our findings can be explained by multiple factors, ranging from differences in dispersal or effective population size, to spatial variation in genomic or phenotypic traits and to structure induced by past histories of host‐adapted populations. Other infrequently measured factors, such as differences in recombination rates, may also play a role. Our result adds to recent work as a further caution against assumptions of simple relationships between genomic and adaptive phenotypic differentiation.     

农林复合系统中滇龙胆形态和生物量变化研究

以幼龄茶树(Camellia sinensis var.sinensis)、十年茶树、木瓜(Chaenomeles sinensis)、尼泊尔桤木(Alnus nepalensis)、大叶桉(Eucalyptus robusta)间作模式和荒坡栽培滇龙胆为研究对象,测定其株高、茎粗、叶长、根粗、生物量等性状,采用单因素方差分析、多重比较、相关性分析和多元逐步回归分析对其形态和生物量数据进行分析。结果表明:荒坡栽培、滇龙胆与十年茶树间作其株高最高,为(37.32±8.36)cm和(37.31±9.62)cm,与大叶桉间作其株高最低,为(19.08±12.40)cm;荒坡栽培植株茎粗数值最高,为(0.36±0.13)cm,大叶桉间作茎粗数值最低,为(0.23±0.04)cm;与茶树间作植株根系最长,为(18.42±7.23)cm和(17.71±7.34)cm,与尼泊尔桤木间作植株须根数最少,为(7.32±2.23)cm;根粗受栽培模式影响不显著(P0.05);荒坡栽培全株生物量最高,为(14.52±13.37)g,大叶桉间作全株生物量最低,为(2.17±1.42)g。相关性分析显示,株高、茎粗和须根数与全株生物量呈极显著的正相关(R=0.514,P0.01;R=0.510,P0.01;R=0.339,P0.01)。根茎比与全株生物量呈极显著的负相关(R=-0.295,P0.01)。多元逐步回归分析显示,各性状对全株生物量的积累贡献程度不同,株高、茎粗、须根数和根粗是影响滇龙胆生物量积累的主要性状。荒坡栽培光照条件较好,有利于植株生物量的积累;与桉树间作,滇龙胆可能受生物与非生物胁迫的共同影响,使其植株矮小,生物量偏低。6种栽培模式中荒坡、滇龙胆与茶树、木瓜间作栽培为高产模式。研究结果可为农林药用复合种植中物种搭配、时间和空间种植技术优化升级以及管理提供理论依据。     

Carbon and nitrogen dynamics in bioenergy ecosystems: 1. Model development,validation and sensitivity analysis

Biofuel made from conventional (e.g., maize (Zea mays L.)) and cellulosic crops (e.g., switchgrass (Panicum virgatum L.) and Miscanthus (Miscanthus × giganteus)) provides alternative energy to fossil fuels and has been considered to mitigate greenhouse gas emissions. To estimate the large‐scale carbon and nitrogen dynamics of these biofuel ecosystems, process‐based models are needed. Here, we developed an agroecosystem model (AgTEM) based on the Terrestrial Ecosystem Model for these ecosystems. The model was incorporated with biogeochemical and ecophysiological processes including crop phenology, biomass allocation, nitrification, and denitrification, as well as agronomic management of irrigation and fertilization. It was used to estimate crop yield, biomass, net carbon exchange, and nitrous oxide emissions at an ecosystem level. The model was first parameterized for maize, switchgrass, and Miscanthus ecosystems and then validated with field observation data. We found that AgTEM well reproduces the annual net primary production and nitrous oxide fluxes of most sites, with over 85% of total variation explained by the model. Local sensitivity analysis indicated that the model sensitivity varies among different ecosystems. Net primary production of maize is sensitive to temperature, precipitation, cloudiness, fertilizer, and irrigation and less sensitive to atmospheric CO₂ concentrations. In contrast, the net primary production of switchgrass and Miscanthus is most sensitive to temperature among all factors. Nitrous oxide fluxes are sensitive to management in maize ecosystems, and sensitive to climate factors in cellulosic ecosystems. The developed model should help advance our understanding of carbon and nitrogen dynamics of these biofuel ecosystems at both site and regional levels.