Published online 2 October 2006
Published in Crop Sci 46:2446-2452 (2006)
© 2006 Crop Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
FORAGE & GRAZINGLANDS
Alfalfa Fiber Estimation in Mixed Stands and Its Relationship to Plant Morphology
D. Parsons,
J. H. Cherney* and
H. G. Gauch, Jr.
Dep. of Crop and Soil Sci., Cornell Univ., Ithaca, NY 14853
* Corresponding author (jhc5{at}cornell.edu)
 |
ABSTRACT
|
|---|
In New York, most alfalfa (Medicago sativa L.) is grown in mixed stands with grass, and models for estimating neutral detergent fiber (NDF) are not used. The objectives of this experiment were to test the suitability of existing equations for estimating NDF of the alfalfa component of mixed stands, and to better understand how alfalfa is affected by the presence of grass. Stands of first-cut alfalfa and grass (1090% grass) were sampled at two experimental sites and producers' fields in 19 New York counties during May and June 2004 and 2005. A range of plant measurements and environmental characteristics were recorded. The predictive equations for alfalfa quality (PEAQ) and other models were examined for applicability to mixed stand alfalfa NDF estimation. The R2 values ranged from 0.80 to 0.87 and RMSE ranged from 20.5 to 25.2 g kg1. The most biased model was PEAQ, possibly due to the lower cutting height used to generate the PEAQ equation than the cutting height used in this study. For producer fields, a model based on alfalfa height was the best one-variable model, with R2 of 0.84 and RMSE of 22.7 g kg1. Presence of grass increases the number of nodes and increases alfalfa height; however, the relationship between alfalfa height and NDF is not changed, suggesting that the predictive ability of models based on alfalfa height is not affected by the percentage of grass in the sward.
Abbreviations: GCANOPY, height of the grass canopy in the sample area GDD, growing degree day GRASS%, actual percentage of grass in the sample LC, lack of correlation MAXHT, height of the tallest alfalfa stem in the sample area MSD, mean squared deviation NDF, neutral detergent fiber NU, nonunity slope PEAQ, predictive equations for alfalfa quality RMSE, root mean square error SB, squared bias SBC, Schwarz Bayesian criterion
 |
INTRODUCTION
|
|---|
A PROPERLY-TIMED SPRING ALFALFA HARVEST sets the stage for good harvest management throughout the rest of the season. Timing of spring forage harvest should be based on NDF (Cherney et al., 2006). Methods of estimating NDF for use as a harvest decision aid must be quick, simple, and reasonably accurate. Models based on weather, chronological age, and plant morphology (Fick et al., 1994) have been developed to estimate alfalfa NDF. The most commonly used of these are the PEAQ equations (Hintz and Albrecht, 1991). Equations using the tallest stem and maturity of the most mature stem in the sample gave acceptable RMSE compared with more complex methods involving mean stage determinations (Kalu and Fick, 1981; Allen and Fick, 1990; Fick and Janson, 1990). The Wisconsin PEAQ equations have been modified for other regions of the USA, including Ohio (Sulc, 1996; Sulc et al., 1997) and New York (Cherney, 1995). The original PEAQ equations have been evaluated in New York, Pennsylvania, Ohio, California, and Wisconsin (Sulc et al., 1997), and prediction errors were small enough to suggest the PEAQ equations are robust across a range of environments.
Although most of the alfalfa in the USA is grown in pure stands, this is not the case in the northeastern USA. In New York it is estimated that >80% of the alfalfa land area is seeded with mixtures of alfalfa and grasses (Cherney et al., 2006). Added difficulties with mixtures include estimating the proportion of grass in the stand, estimating the NDF of the grass portion, and knowing how the grass portion affects the NDF of the alfalfa. Parsons et al. (2006b) has suggested models for estimating total sward NDF of mixed stands based primarily on alfalfa height and the percentage of grass in the stand. Informal observations of mixed stands in New York suggest that alfalfa height is positively correlated with the percentage of grass in the stand. The PEAQ equations have been tested and recommended for use only with pure stands of alfalfa that contain no grass or weeds (Sulc et al., 1997). Because the equations rely on alfalfa height, if other plants in the sward affect alfalfa height, then the validity of the models are uncertain. As a result, equations for estimation of alfalfa NDF have not been used to estimate the NDF of the alfalfa component of mixed stands. The objectives of this study were to (i) test existing alfalfa NDF models for their ability to estimate alfalfa NDF in mixed stands; (ii) develop alternative models for estimating the alfalfa NDF of mixed stands; (iii) understand how the presence of grass in the mixed stand affects the morphology of alfalfa; and (iv) determine the validity of applying equations for pure alfalfa to mixed stands.
 |
MATERIALS AND METHODS
|
|---|
Sampling of Mixed Stands
Alfalfagrass mixed stands were sampled in the spring at two experimental sites and 150 producers' fields in 19 New York counties during May and June 2004 and 2005, as described in Parsons et al. (2006b). The experimental sites were the Cornell University Caldwell Field Research Farm (42°27'1'' N, 76°27'35'' E, 276 m, 02% slope) near Ithaca, NY, and Mount Pleasant Research Farm (42°27'36'' N, 76°22'12'' E, 520 m, 06% slope) near Dryden, NY. The soil at Caldwell Field was a Niagara silt loam (fine-silty, mixed, active, mesic Aeric Endoaqualfs) and the soil at Mount Pleasant was a Mardin silt loam (coarse-loamy, mixed, active, mesic Typic Fragiudepts). The experimental design at each site was a randomized complete block design with four blocks. Each block included three different alfalfagrass species mixtures at two grass seeding rates, and one plot of pure-sown alfalfa, giving a total of 28 plots. The term pure-sown is used here to denote a stand that was planted as pure alfalfa and which may or may not have been invaded by grass species. Each plot measured 2.7 by 6 m (16.2 m2) with 0.15 m between plots and 0.3-m alleys between blocks. Plots were seeded on 19 May 2003 at Caldwell Field and 23 May 2003 at Mount Pleasant. All plots were seeded at Caldwell Field with Hytest 340PLH alfalfa at 13.4 kg ha1 and at Mount Pleasant with Hytest 104PLH alfalfa at 13.4 kg ha1 using a Brillion seeder (Brillion Farm Equipment, Brillion, WI). Grass plots were seeded using a Carter seeder (Carter Mfg., Brookston, IN). Richmond timothy (Phleum pratense L.) was seeded at 3.4 and 6.7 kg ha1, Okay orchardgrass (Dactylis glomerata L.) was seeded at 4.5 and 9.0 kg ha1, and Rival reed canarygrass (Phalaris arundinaceae L.) was seeded at 6.7 and 13.4 kg ha1. Seeding rates were based on pure live seed. Lime, P, and K were applied based on soil test results.
Producer fields and plots were sampled when alfalfa height reached or exceeded 30 cm. To define a portion of the field or plot as the sample area, an area of
1 m2 was visually identified in 2004, and in 2005 a circular quadrat of comparable area was used. The data collected and variable abbreviations are summarized in Table 1. Height of the tallest alfalfa stem in the sample area was measured to the terminal bud (MAXHT). The alfalfa maturity categories of Kalu and Fick (1981) were used to assign a numerical value to the most mature stem in the sample area (MAXSTAGE). The major grass species was recorded (GSPECIES). The height of the tallest grass tiller in the sample area was measured by fully extending the leaf (GMAXHT). The average grass canopy height of the sample area was measured with no extension of leaves (GCANOPY). The developmental stage of the most mature grass tiller in the sample area was determined using the staging system of Moore and Moser (1995) (GMAXNDX). Determination of the index number using this system requires knowledge of the total number of leaves or nodes that will appear before reaching the next development stage. Because this system requires prior knowledge of development norms for each member species, a simplified grass staging system described by Parsons et al. (2006b) was also used (GMAXSTG).
View this table:
[in this window]
[in a new window]
|
Table 1. Descriptions of variables evaluated as potential predictors of mixed alfalfa and grass neutral detergent fiber.
|
|
Time of sampling was recorded and converted to a decimal number (TIME); for example, 2:30 pm was converted to 14.5. The percentage of grass in the sample area was visually estimated (GEST). A representative sample of 500 to 750 g of alfalfa and grass was hand clipped from the sample area at a height of 10 cm, an approximation of typical cutting height in New York. Date of sampling was transformed to day of the year (DOY), the number of days from the beginning of the year. The altitude of the field was recorded (ALTF), as were the geographic co-ordinates. Co-ordinates of the fields were overlayed with the co-ordinates of all New York weather stations using Manifold (Enterprise Edition 6.50, CDA International Ltd., San Mateo, CA). Voronoi cells were created to determine the nearest weather station for each field, thus enabling the calculation of individualized growing degree days (GDDs). Accumulated growing degree days were calculated using both base 0°C (GDD0) and base 5°C (GDD5). Accumulation of GDDs was initiated when the mean temperature exceeded the base for five consecutive days. The altitude of the nearest weather station (ALTWS) was used as a potential explanatory variable.
Samples were separated into alfalfa and grass fractions, oven dried at 60°C until a constant dry weight was reached, and the actual percentage of grass in the sample (GRASS%) was calculated. Samples with GRASS% < 10 or > 90 were not used for further analysis. Samples were ground to pass through a 1-mm screen. Samples (0.25 g) were analyzed for NDF concentration using the procedure described by Van Soest et al. (1991), using the ANKOM (Macedon, NY) fiber analyzer with filter bags. Sodium sulfite and
-amylase were used in the NDF procedure.
Alfalfa Morphology Experiment
Measurements of alfalfa morphology were obtained from the Caldwell Field experimental plots described above. The low and high seeding rates defined in the original experiment were not distinguished as separate treatment levels. The plots were harvested on 20 May and 8 June 2005. From each plot, single 625-cm2 quadrats were harvested with hand shears to a height of 10 cm. Alfalfa and grass were separated and the grass portions were oven dried to a constant weight at 60°C to determine the concentration of dry matter in each sample. The five tallest alfalfa stems from each sample were selected for further measurements and the remaining alfalfa portions were oven dried to a constant weight at 60°C to determine the concentration of dry matter. The following measurements were made on the five tallest stems. Alfalfa stem height was measured to the terminal bud. Alfalfa stem diameter was measured with digital callipers at 2.5 cm above the cutting height. Easily separable nodes above the cutting height were counted. Lengths of the first two internodes at the 20 May harvest, and the first four internodes at the 8 June harvest, were measured. Leaf-to-stem mass ratios were determined by hand separating leaf from stem for the combined five tallest stems, oven drying at 60°C, and weighing. Replicates of leaf and stem fractions were combined and analyzed for NDF concentration. The procedure described by Van Soest et al. (1991) was followed using the ANKOM fiber analyzer with filter bags. Sodium sulfite and
-amylase were used in the NDF procedure.
Predictive Equations and Statistical Analysis
The original PEAQ equation for NDF estimation developed by Hintz and Albrecht (1991) is NDF = 168.9 + 2.7(MAXHT) + 8.1(MAXSTAGE). Cherney and Sulc (1997) developed an equation for NDF estimation in New York fields, based only on alfalfa height. This equation, labeled NYPQ, is NDF = 122.7 + 3.1(MAXHT). Parsons et al. (2006a) developed a similar equation, also based only on alfalfa height. This equation, labeled NYHT, is NDF = 67.7 + 4.1(MAXHT). In addition, Parsons et al. (2006a) developed an equation based on alfalfa height and GDDs. This equation, labeled NYGD, is NDF = 68.9 + 3.4(MAXHT) + 0.14(GDD5). Model validation was performed by regressing actual NDF values on the predicted NDF values. A number of parameters were used in evaluating the equations because no single statistical test can adequately describe the goodness of fit of the model. The coefficient of determination (r2 or R2) and RMSE were determined. Root mean square error has the same units as the variable predicted, and in model validation is the prediction error of the model. All regression equations were tested for intercepts at the origin (a = 0) and unitary coefficients (b = 1). Kobayashi and Salam (2000) presented reasons why these parameters are not entirely satisfactory for model evaluation and promoted the use of mean squared deviation (MSD) and its components as more informative parameters. The MSD reflects discrepancy between a model and the data and is a direct measure of predictive success. Gauch et al. (2003) proposed the partitioning of MSD into the components of squared bias (SB), nonunity slope (NU), and lack of correlation (LC), to provide further insight into model performance. The three components have distinct meanings and simple geometric interpretation, with SB relating to translation, NU relating to rotation, and LC relating to scatter.
For model development, the PROC RSQUARE variable selection procedure of SAS, v. 9.1 software (SAS Inst., 2003) was used to identify models that maximized the coefficient of determination (r2 or R2) and minimized RMSE and the Schwarz Bayesian criterion (SBC). In model development, RMSE is the calibration error of the model. The SBC is a better indicator of predictive accuracy than R2 or RMSE. PROC GLM in SAS, v. 9.1 software (SAS Inst., 2003) was used to determine whether the resultant F value from the addition of the second explanatory variable in a two-variable model was statistically significant.
Alfalfa morphology data were analyzed using PROC GLM in SAS, v. 9.1 software (SAS Inst., 2003). Plant species, seeding rates, and harvest dates were considered fixed effects and replicates were considered random effects. Dunnett's test was used to compare pure alfalfa (the control) with each alfalfagrass mix. Principal component scores were derived using PROC PRINCOMP in SAS, v. 9.1 software (SAS Inst., 2003), and standardized by variable. A biplot for the first two principal components was constructed using the method of Gabriel (1971).
 |
RESULTS AND DISCUSSION
|
|---|
Evaluation of Existing Equations
Results of fitting model estimates from PEAQ, NYPQ, NYHT, and NYGD to observed NDF values are contained in Table 2. The range of values for R2 (0.800.87) and RMSE (20.525.2 g kg1) are of comparable magnitude to the PEAQ equations of Hintz and Albrecht (1991). The NYGD equation has the best goodness of fit, with R2 of 0.87 and RMSE of 20.5 g kg1. The other equations have lower goodness of fit. For PEAQ and NYPQ, the b values were significantly >1 and the a values were significantly <0. For NYHT and NYGD, the b values were significantly <1 and the a values were significantly >0. The order of b values ranging from closest to 1 to furthest from 1 is NYPQ, NYGD, NYHT, and PEAQ. Similarly, the order of a values ranging from closest to 0 to furthest from 0 was NYPQ, NYGD, NYHT, and PEAQ.
View this table:
[in this window]
[in a new window]
|
Table 2. Coefficient of determination (R2), root mean square error (RMSE), slope (b), and y intercept (a) derived from fitting equation estimates to observed forage neutral detergent fiber values.
|
|
These results make it difficult to assess the comparative accuracy of the models. For example, the NYGD equation which had the greatest R2 and lowest RMSE had greater deviations in slope and intercept than NYPQ. Thus, it is difficult to compare the strengths and limitations of the models solely using these statistics. Partitioning MSD provides a way to better understand predictive success as a more definitive criterion. In Fig. 1, the order of lowest to greatest MSD is NYGD, NYPQ, NYHT, and PEAQ. The LC for NYGD was lowest (420), whereas LC values for PEAQ (605), NYPQ (634), and NYHT (634) were of similar magnitude. Values for NU were lowest for NYPQ (24) and greatest for NYHT (85), and all were comparatively less than corresponding LC values. The most striking difference between the models was the high SB value for PEAQ (550) compared with NYHT (85), NYGD (43), and NYPQ (24). Simply looking at the results in Table 2 it is difficult to appreciate the bias in the PEAQ model in comparison with the other models. Partitioning of MSD helps to clarify the relative strengths and weaknesses of the models. The relationship between regression lines of model estimates and a 1:1 line is shown in Fig. 2. The PEAQ model appears the most biased in terms of slope and intercept and it is evident that PEAQ deviates from a 1:1 prediction of NDF, particularly at lower values of NDF. Similar results were reported by Parsons et al. (2006a) using the PEAQ model for estimating alfalfa in pure-sown stands in New York. These authors suggested that the PEAQ equation was developed based on a cutting height of 3.8 cm, whereas their study was based on a cutting height of 10 cm, which likely affected bias. Shorter cutting heights would include material of greater NDF near the base of the plant; however, the effect of this extra 6.2 cm would be diluted with increasing plant height. This could account for the observed overestimation of NDF at shorter cutting heights. At a target NDF for alfalfa harvest of 400 g kg1 the divergence between PEAQ and the 1:1 line is minimal, and thus the bias may not be of practical relevance in pure stands. However, in mixed stands where harvest should occur when alfalfa NDF is comparatively less, the bias of PEAQ may be more significant. It may be possible to adjust PEAQ estimations for different cutting heights to give more accurate estimations.

View larger version (14K):
[in this window]
[in a new window]
|
Fig. 1. Components of mean squared deviation (MSD) of regression models used to estimate the neutral detergent fiber of the alfalfa component of mixed stands. The equations used are PEAQ (Hintz and Albrecht, 1991), NYPQ (Cherney and Sulc, 1997), NYHT (Parsons et al., 2006a), and NYGD (Parsons et al., 2006a). LC, lack of correlation; NU, nonunity slope; SB, squared bias.
|
|
The other models examined (NYPQ, NYHT, and NYGD) showed less bias than the PEAQ equation. For NYPQ, and to a lesser extent NYHT and NYGD, LC constituted the majority of the MSD, indicating that the errors due to SB and NU were comparatively small. These results suggest that these three models successfully estimate the NDF of the alfalfa component of mixed stands.
Development of New Equations
The data were used to develop new equations to estimate the alfalfa NDF of mixed stands. Table 3 shows the variable selection results, grouped according to location. Locations included the Mount Pleasant and Caldwell Field experimental sites, producer fields, and a combination of all locations. For each location, three equations are presented: (i) the best one-variable model, (ii) the best two-variable model including only plant-derived variables in the analysis, and (iii) the best two-variable model including all variables in the analysis. Models were selected on the basis of having the lowest SBC.
For all locations, the ranges of R2 (0.770.95) and RMSE (13.125.2 g kg1) are of comparable magnitude to those of the PEAQ equations of Hintz and Albrecht (1991). In all two-variable models, the addition of the second explanatory variable was statistically significant. The pattern of R2 and RMSE at the Caldwell Field and Mount Pleasant sites were similar. The best one-variable models (Eq. [1] and [4] in Table 3) were based on GDD. The two-variable models, based on the plant-derived variables MAXHT and GCANOPY (Eq. [2] and [5]), had a lower R2 and a greater RMSE. The best two-variable models were based on MAXHT and GDD (Eq. [3] and [6]), with R2 values of 0.85 and 0.95 and an RMSE of 17.8 and 13.1 g kg1.
For the producer field samples, the best one-variable model (Eq. [7]) was based on MAXHT. The best two-variable models were based on MAXHT and MAXSTAGE (Eq. [8]), and MAXHT and GDD5 (Eq. [9]). When all samples were combined, the best one-variable model (Eq. [10]) was also based on MAXHT. The best two-variable models were based on MAXHT and GCANOPY (Eq. [11]), and MAXHT and GDD5 (Eq. [12]).
For both Mount Pleasant and Caldwell Field, a GDD variable was the best single explanatory variable. Fick et al. (1994) detail numerous studies that demonstrate the use of GDD for estimation of alfalfa NDF. Our results confirm this relationship; however, they also show that GDD had less predictive power for the producer fields than for the experimental fields. Both Mount Pleasant and Caldwell Field had a weather station on site, ensuring that GDD data was accurate. In contrast, the nearest weather station for individual producer fields could be a considerable distance away, and the difference in altitude (ALTD) might also be great. Allen and Beck (1996) reported decreased goodness of fit for GDD models as distance from weather stations increased. The excellent results for Caldwell Field suggested that GDD has potential as a tool for NDF estimation if accurate local weather data can be obtained. However, Parsons et al. (2006b) showed that models based on GDD data were less successful when used to predict alfalfa NDF for years on which the model was not built. Therefore, although calibration errors may be low (encouraging use of the model), prediction errors may be inflated when an existing model is used on new data.
The goodness of fit for the Mount Pleasant models is less than for the Caldwell Field models, but comparable with other studies such as the PEAQ equations of Hintz and Albrecht (1991). The poorer results at Mount Pleasant may be reflective of the stands, which were generally uneven in plant distribution and contained more weeds. Although the Caldwell Field plots were better, they were not free of weeds. These observations suggest that equations for estimating NDF of the alfalfa component of mixed stands may be somewhat robust to the presence of weeds. The effect of weeds on alfalfa NDF may depend on both the competitive ability and density of the weeds. Further investigation of this observation would be a practical area of study, particularly in light of the fact that in New York very few stands of alfalfa or alfalfagrass are free of weeds.
In contrast to the experimental plots, MAXHT was the variable that explained the most variation in alfalfa NDF for producer fields. Although goodness of fit for the model with MAXHT and MAXSTAGE is statistically greater (Table 2), we consider this increase to be not biologically significant and thus the extra measurement of alfalfa stage is not warranted. For pure-sown alfalfa stands in New York, Parsons et al. (2006a) found that that addition of alfalfa stage to a model containing alfalfa height was not statistically significant. In this study, producer fields were sampled over a wide geographic area, and included a range of alfalfa and weed densities. These results confirm the robustness of MAXHT and its usefulness in estimating alfalfa NDF not only in pure-sown stands, but also for the alfalfa component of mixed stands.
Alfalfa Morphology in Mixed Stands
Blocks and interactions between harvest date and species were nonsignificant for all variables at a probability level of 0.05. The effects of species and harvest date on alfalfa morphological characteristics are shown in Table 4. Species had a significant effect on alfalfa height. The height of alfalfa growing with each grass species was significantly greater than the height of pure-sown alfalfa. Alfalfa height increased from Harvest 1 to Harvest 2. There was no significant difference in Alfalfa stem diameter either across harvests or across species treatments.
View this table:
[in this window]
[in a new window]
|
Table 4. Treatment means for alfalfa height, alfalfa stem diameter, number of alfalfa nodes, alfalfa internode lengths (IN), alfalfa leaf-to-stem mass ratio (LSR), alfalfa yield, grass yield, and bluegrass yield.
|
|
The number of alfalfa nodes increased from Harvest 1 to Harvest 2. Alfalfa growing with grass had more nodes than pure-sown alfalfa. There was no significant difference in the length of Internode 1 either across harvests or across species treatments. The length of Internode 2 increased from Harvest 1 to Harvest 2. Internode 2 lengths for pure-sown alfalfa were not significantly different from species mixtures. The lengths of Internodes 3 and 4 for pure-sown alfalfa were not significantly different from species mixtures.
Leaf-to-stem mass ratio decreased from Harvest 1 to Harvest 2. Pure-sown alfalfa had a greater leaf-to-stem mass ratio than alfalfa growing with orchardgrass or reed canarygrass, but was not significantly different from alfalfa growing with timothy. Alfalfa yield increased from Harvest 1 to Harvest 2. Yield of alfalfa was greater in pure-sown stands than in orchardgrass or timothy, but was not significantly different from alfalfa growing with reed canarygrass. Grass yield increased from Harvest 1 to Harvest 2. The yield of bluegrass (Poa spp.), a common weed of alfalfa stands, increased from Harvest 1 to Harvest 2. Alfalfa growing in pure-sown stands had a greater yield of bluegrass than alfalfa growing in mixed stands.
Principal component analysis was used to further clarify the relationship between the presence of grass and alfalfa morphology. The first two principal components accounted for 81% of the total variance. Figure 3 is the biplot for the first two principal components and shows the relationship between samples and variables. Combining replicates, there were four species and two harvest dates, giving a combination of eight samples. The eleven variables plotted include leaf-to-stem mass ratio, length of Internode 1, length of Internode 2, grass yield, alfalfa yield, bluegrass yield, alfalfa stem diameter, alfalfa height, number of alfalfa nodes, alfalfa leaf NDF, and alfalfa stem NDF. The cosine of the angle between two points of the same type (i.e., sample or variable) estimates the correlation between the two points; thus points that are clustered are highly correlated. Trends from the principle component analysis are not directly comparable with the analysis of variance results; however some commonalities are expected. There is a cluster of variables that are highly correlated, composed of leaf NDF, stem NDF, nodes, and height (Fig. 3). For example, the correlation between height and nodes is 0.994. The eigenvector for leaf-to-stem mass ratio points in the opposite direction; thus, the four clustered variables are highly negatively correlated to leaf-to-stem mass ratio. As an example, the correlation value between nodes and leaf-to-stem mass ratio is 0.991. This set of variables reflects alfalfa maturation. Moving from left to right in Fig. 3, as the alfalfa plant matures it increases in height, number of nodes, and stem and leaf NDF, whereas the leaf-to-stem ratio decreases.

View larger version (22K):
[in this window]
[in a new window]
|
Fig. 3. Biplot for sample means and variables. Sample points, signified by a letter and number, are composed of mixes of alfalfa and orchardgrass (O), alfalfa and timothy (T), alfalfa and reed canarygrass (R), pure-sown alfalfa (A), Harvest 1 (1), and Harvest 2 (2). The variables include leaf-to-stem mass ratio (LSR), length of Internode 1 (IN1), length of Internode 2 (IN2), grass yield (Grass), alfalfa yield (Alfalfa), bluegrass yield (Bluegrass), alfalfa stem diameter (Diameter), alfalfa height (Height), number of alfalfa nodes (Nodes), alfalfa leaf neutral detergent fiber (LNDF), and alfalfa stem neutral detergent fiber (SNDF).
|
|
The other important trend in Fig. 3 is that of increasing grass yield, which is roughly at a 45° angle to the trend in maturity. There is a moderate correlation between grass yield and the lengths of Internode 1 (0.482) and Internode 2 (0.679). This suggests that with increasing yield of grass, alfalfa may respond by increasing internode length. The correlation is greater with Internode 2. Because the period of extension for Internode 2 is later than for Internode 1, it is possible that the response may be greater as grass yield and the subsequent competition between alfalfa and grass increases.
The negative correlation between grass yield and alfalfa stem diameter (0.393) suggests that alfalfa may also respond to grass competition by decreasing stem diameter, an accompanying response to increasing internode length. In Fig. 3 it appears that this correlation is much stronger. However, Fig. 3 only takes into account Principle Components 1 and 2, whereas correlation matrix values take into account all dimensions. Volenec et al. (1987) reported a decrease in stem diameter as plant population increased in pure-sown stands of alfalfa. Although in this study there was no significant species effect on stem diameter in the ANOVA (Table 4), decreased stem diameter should not be disregarded as a possible response to competition from grass.
Bluegrass is positively correlated with alfalfa yield (0.722) and negatively correlated with grass yield (0.546). This is consistent with observations in New York that bluegrass is widespread in fields of alfalfa, but less prolific where there is a significant sown grass component in the sward.
The location of sample data points (Fig. 3) also contributes to understanding the direction of growth processes in the sward. Markers for the first, less mature harvests are located in the left of the biplot, whereas markers for the second, more mature harvests are in the right. Markers for the first harvest are clustered, signifying that at this stage there is less competition, whereas markers for the second harvest are more dispersed. Particularly for the second harvest, markers for pure-sown alfalfa and alfalfa with reed canarygrass are located toward the top, and are associated with greater yields of alfalfa and bluegrass. This is consistent with observations that reed canarygrass is less competitive with alfalfa than other grasses grown in New York. Markers for alfalfa with orchardgrass and timothy are toward the bottom and signify high grass yield and strong grass competition. These observations are consistent with the typical pattern of growth of these species combinations in New York fields.
 |
CONCLUSIONS
|
|---|
Existing models for pure-sown alfalfa NDF can be used to estimate the alfalfa component of mixed stands. This is significant for producers in New York and other locations where mixed stands are common. Although equations exist to estimate the total NDF of mixed stands (Parsons et al., 2006b), this finding is still pertinent. Previously, producers have been advised against using the PEAQ equation in fields where alfalfa may be dominant but a small amount of grass is present.
Acknowledging that PEAQ was originally calibrated and only advocated for pure stands of alfalfa cut at a 3.5-cm stubble height, an issue the results raise is the bias of PEAQ when it is used for alfalfa not cut at this height. In reality, producers cut at varying stubble heights, based on such factors as available machinery and stoniness of the field. Adjusting PEAQ (and other equations) for expected cutting height could help reduce the potential for biased NDF estimates.
The additional equations presented in this paper offer alternatives for predicting the NDF of the alfalfa component of mixed stands. The equations could be tested in other locations outside New York to assess their predictive ability. The results of this study should increase confidence in MAXHT as the premium explanatory variable for alfalfa NDF in both pure and mixed stands. Samples were taken from a wide geographical area and from fields differing in many respects. No attempt was made to exclude fields from the study because of a few weeds or to restrict the fields selected to having a narrow range of grass species.
The alfalfa morphology results suggest that the presence of grass causes an increase in alfalfa height by stimulating an increase in the number of nodes, consequently increasing NDF. It is possible that the presence of grass may also increase internode length and reduce stem diameter; however, these effects were less evident. Although grass causes a morphological response in alfalfa, the relationship between alfalfa height and NDF remained consistent across the height range in the study. This suggests that an alfalfa plant growing in a pure-sown stand is similar in terms of forage chemistry as an alfalfa plant of the same height growing in a mixed stand. It follows that for practical purposes the PEAQ equation and other alfalfa NDF predictive equations are applicable to the alfalfa component of mixed stands.
 |
ACKNOWLEDGMENTS
|
|---|
The authors thank Sam Beer, Kai Ming Zhao, Molly Lebowitz, Jen Beckman, Peter Barney, Aaron Gabriel, Jeff Miller, Bruce Tillapaugh, Rick Faucett, Michael Hunter, Michael Davis, Aysin Bilgili, and Leon Hatch for assistance with harvesting and analysis. This research was supported by a Kieckhefer Adirondack Fellowship.
Received for publication March 2, 2006.
 |
REFERENCES
|
|---|
- Allen, M., and J. Beck. 1996. Relationship between spring harvest alfalfa quality and growing degree days. p. 1625. In Proc. National Alfalfa Symposium, 26th, East Lansing, MI. 45 Mar. 1996. Certified Alfalfa Seed Council, Davis, CA.
- Allen, S.J., and G.W. Fick. 1990. On-farm testing of mean stage by count as a predictor of alfalfa forage quality. p. 185. In Agronomy abstracts. ASA, Madison, WI.
- Cherney, J.H. 1995. Spring alfalfa harvest in relation to growing degree days. p. 2936. In Proc. Natl. Alfalfa Symp., 25th, Syracuse, NY. 2728 Feb. 1995. Certified Alfalfa Seed Council, Woodland, CA.
- Cherney, J.H., D.J.R. Cherney, and D. Parsons. 2006. Grass silage management issues. p. 3749. In Proc. from Silage for dairy farms: Growing, Harvesting, Storing, and Feeding, Harrisburg, PA. 2325 Jan. 2006. NRAES-181. Natural Resource, Agriculture, and Engineering Service, Ithaca, NY.
- Cherney, J.H., and R.M. Sulc. 1997. Predicting first cutting alfalfa quality. p. 5365. In Silage: Field to feedbunk. Proc. from the North American Conf., Hershey, PA. 1113 Feb. 1997. NRAES-99. Northeast Regional Agricultural Engineering Service, Ithaca, NY.
- Fick, G.W., and C.G. Janson. 1990. Testing mean stage as a predictor of alfalfa forage quality with growth chamber trials. Crop Sci. 30:678682.[Abstract/Free Full Text]
- Fick, G.W., P.W. Wilkens, and J.H. Cherney. 1994. Modeling forage quality changes in the growing crop. p. 757795. In G.C. Fahey, Jr. et al. (ed.) Forage quality, evaluation, and utilization. ASA, CSSA, and SSSA, Madison, WI.
- Gabriel, K.R. 1971. The biplot graphic display of matrices with application to principal component analysis. Biometrika 58:453467.[Abstract/Free Full Text]
- Gauch, H.G., J.T.G. Hwang, and G.W. Fick. 2003. Model evaluation by comparison of model-based predictions and measured values. Agron. J. 95:14421446.[Abstract/Free Full Text]
- Hintz, R.W., and K.A. Albrecht. 1991. Prediction of alfalfa chemical composition from maturity and plant morphology. Crop Sci. 31:15611565.[Abstract/Free Full Text]
- Kalu, B.A., and G.W. Fick. 1981. Quantifying morphological development of alfalfa for studies of herbage quality. Crop Sci. 21:267271.
- Kobayashi, K., and M.U. Salam. 2000. Comparing simulated and measured values using mean squared deviation and its components. Agron. J. 92:345352.[Abstract/Free Full Text]
- Moore, K.J., and L.E. Moser. 1995. Quantifying developmental morphology of perennial grasses. Crop Sci. 35:3743.[Abstract/Free Full Text]
- Parsons, D., J.H. Cherney, and H.G. Gauch. 2006a. Estimation of spring forage quality for alfalfa in New York State. Forage Grazingl. doi:10.1094/FG-2006-0323-01-RS (available at http://www.plantmanagementnetwork.org/pub/fg/research/2006/alfalfa/).
- Parsons, D., J.H. Cherney, and H.G. Gauch. 2006b. Estimation of preharvest fiber content of mixed alfalfagrass stands in New York. Agron. J. 98:10811089.[Abstract/Free Full Text]
- SAS Institute. 2003. The SAS system. v. 9.1. SAS Inst., Cary, NC.
- Sulc, R.M. 1996. Equations for predicting quality of alfalfa. p. 115124. In 1996 Proc. Tri-State Dairy Nutrition Conf., Fort Wayne, IN. 1415 May 1996. Ohio State Univ., Columbus.
- Sulc, R.M., K.A. Albrecht, J.H. Cherney, M.H. Hall, S.C. Mueller, and S.B. Orloff. 1997. Field testing a rapid method for estimating alfalfa quality. Agron. J. 89:952957.[Abstract/Free Full Text]
- Van Soest, P.J., J.B. Robertson, and B.A. Lewis. 1991. Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition. J. Dairy Sci. 74:35833597.[Abstract]
- Volenec, J.J., J.H. Cherney, and K.D. Johnson. 1987. Yield components, plant morphology, and forage quality of alfalfa as influenced by plant population. Crop Sci. 27:321326.[Abstract/Free Full Text]