Abstract
Carrots (Daucus carota L. var sativus) in the midwestern United States are generally raised commercially on sandy-loam or muck soils. Some smaller growers may not have the luxury of an ideal soil type, so raised beds are often recommended to improve taproot growth. To determine specific yield and shape responses of carrots to raised-bed culture in heavy soil, three hybrids (‘Cupar’, ‘Envy’, and ‘Naval’) were grown in 2019 and 2020 on raised or flat beds in clay-loam soil. In addition, two seeding rates (high or low) were incorporated into the split-split plot design. ‘Cupar’ was larger and less misshapen than the other hybrids, but total marketable yields were similar among all three. ‘Naval’ was the most cylindrical. High seeding rate generally resulted in more numerous but smaller marketable roots and no change in total marketable weight. Raised beds caused a 10.6% increase in total marketable weight by allowing more marketable roots classified as “short” (between 76 and 127 mm long, diameter >25 mm) to grow to lengths more than 127 mm long. The average taproot cylindricity was increased and diameter decreased by growing on raised beds. Soil penetration resistance at 7.6 cm was correlated with many of the measured responses.
Carrots are an economically important crop in the upper midwestern United States, with 3493 ha harvested in Minnesota, Michigan, and Wisconsin in 2017; 84% of the acreage was for processing [U.S. Department of Agriculture (USDA) National Agricultural Statistics Service (NASS), 2020]. Soil compaction generally negatively impacts taproot length and yield (Johansen et al., 2015; Pietola, 1995). Typically, carrots are grown on peat/muck or sandy-loam soils that are deeply cultivated to reduce compaction for encouraging straight, uniform taproots. University extension publications and popular press often advise growers and gardeners to plant carrots on raised beds for similar reasons (Delahaut and Newenhouse, 2011; Jauron and Wallace, 2017; Kelley and Phatak, 2009; Smetanka, 2017; Thompson, 2017). Especially on fine-textured mineral soils (not sandy or peat/muck), it is reasonable to assume that raised beds reduce compaction and soil density where the taproot grows. In addition, raised beds can improve water relations in the root zone, especially in wet soils (Bakker et al., 2005). Water relations early in a carrot crop cycle may influence length potential and straightness of the harvested taproot, especially in heavy soils (Brown and Gracie, 2000). A carrot taproot’s ability to avoid areas of low water availability suggests that after germination, drier surface soil in raised beds may encourage longer taproots. However, data to illustrate the specific effects of raised beds on carrots are lacking.
Weed control in carrots can be especially challenging, because carrot seeds are slow to germinate and a canopy is slow to develop. One strategy to reduce weed pressure is to seed at a higher rate, reducing the time to develop a canopy for shading weed seedlings (Colquhoun et al., 2020). However, increased population causes shorter taproots (Thompson, 1969).
Carrot varieties can be classified by their general shape and use characteristics. For example, Chantenay carrots are generally wide and conical-shaped, with medium-length, and are often recommended over other types for heavy or compacted soils. Imperator carrots are long, and generally perform best on sandy-loam or peat/muck soils, although breeders can recommend specific varieties that may perform well on mineral soils. Nantes-type are highly cylindrical and medium-length. Each type has utility in fresh-market or processing industries.
In this study, three carrot hybrids were grown (Chantenay-type, Imperator-type, and Nantes-type) to determine the effects of raised beds and elevated seeding rates, and the interaction between them, on carrot taproot size, shape, and yield in a clay-loam soil.
Materials and Methods
The experiment was planted as a split-split-plot design in 2019 and 2020 at the University of Minnesota Southern Research and Outreach Center in Waseca, MN (lat. 44.08°N, long. 93.52°W). The main plot was carrot hybrid (n = 3; Chantenay-type ‘Cupar’ or Nantes-type ‘Naval’ from Bejo Seeds, Oceano, CA; Imperator-type ‘Envy’ from Seminis Seeds, St. Louis, MO), the first split plot was seeding rate (n = 2, described in the following) and the split-split plot was bed shape (n = 2; flat or raised bed). There were six replicated complete blocks per year.
The soil was a Nicollet clay loam, a fine-textured, mesic Aquic Hapludoll. Soil test properties are described in Table 1. The previous crop was spring wheat with fall fallow. On 29 Apr. 2019 or 22 Apr. 2020, the field was fertilized with phosphorus (triple superphosphate, 0N–20.1P–0K), potassium (KCl, 0N–0P–49.8K), sulfur (gypsum, 18% S), and boron (Granubor; US Borax, Boron, CA; 15% B) as described in Table 1, and cultivated with a rototiller. Raised beds 7.6 to 10.2 cm tall and 36.6 m long with 0.76-m-wide bed tops were formed on 1.5-m centers on 4 May 2019 or 22 Apr. 2020 using a custom bed shaper. To implement the flat bed treatments, a walk-behind tiller (BCS 739; BCS America, Oregon City, OR) was fitted with a hiller/furrower and axle extensions so that the tires would not compact the raised beds. The tiller was set to its most shallow depth, and the furrower was pulled behind the tiller to flatten the raised beds and create the flat bed treatment plots, moving the soil from the raised beds back into the furrows. This was done on 4 May 2019 or 23 Apr. 2020. Flat bed plots were immediately smoothed by hand with a garden rake after leveling. Whole plots (hybrid) were 17.1 m long, evenly divided into two subplots (seeding rate) and four total sub-subplots (bed shape; 4.3 m each). For ≈5 weeks after creation of the bed treatments, the surface soil was allowed to weather into a fine seedbed because cultivation had created large aggregates. Nitrogen (N) fertilizer was added on 11 June 2019 or 28 May 2020 and incorporated with a walk-behind wheel hoe (Table 1). Nitrogen was applied as 25% urea-N and 75% slow-release ESN-N (Nutrien Ag Solutions, Saskatoon, Canada). Fertilizer rates were based on fertility recommendations (Rosen and Eliason, 2005).
Soil and fertility conditions for carrot growth at the Southern Research and Outreach Center (Waseca, MN). Nutrients were all added preplant. The previous crop each year was spring wheat followed by fall fallow. Soil nitrogen and sulfur were not measured.
‘Morex’ barley (Hordeum vulgare L.) was broadcast as a nurse crop at 55 kg⋅ha−1 on 13 June 2019 or 28 May 2020, and also lightly incorporated with a wheel hoe. Carrots were sown the following day. A Nibex 500 Precision Seeder (Andershonstein AB, Vetlana, Sweden) was calibrated for each hybrid and seeding rate (high or low). Seed suppliers recommended 555,750 to 741,000 per ha for ‘Cupar’ and 1,235,000 to 2,099,500 per ha for ‘Envy’ and ‘Naval’. The hybrids used were recommended by the respective seed companies for productivity on mineral soils. A 35% greater population was targeted in the high seeding rate compared with low, and ‘Cupar’ was targeted at 44% of the seeding rate of ‘Envy’ or ‘Naval’. Specifically, ‘Cupar’ was targeted at 592,800 or 800,280 viable seeds per ha, and ‘Envy’ and ‘Naval’ were targeted at 1,333,800 or 1,803,100 viable seeds per ha (low or high within cultivar, respectively). Seeding calculations were based on 1.5-m bed spacing. Supplier-provided values for viability and seed counts were used for planter calibration. Each plot was composed of three planted rows per bed (raised or flat), with outer rows ≈61 cm apart. Each outer row was targeted to contain 37.5% of the total seeding rate per meter of bed, and the inner row was targeted to receive 25% (outer rows were each 150% seeding rate of the inner rows). A 0.6-m area was left between adjacent “hybrid” treatments to accommodate emptying and cleaning the planting boxes between hybrids. Each row was sown in a band using a 65-mm scatter shoe.
Prowl H2O (38.7% pendimethalin; BASF, Florham Park, NJ), AIM EC (22.3% carfentrazone-ethyl; FMC, Philadelphia, PA), Roundup PowerMax (48.7% glyphosate; Bayer, St. Louis, MO), Dakota (26.4% clethodim; Rotam North America, Greensboro, NC), and Lorox (50% linuron; NovaSource, Phoenix, AZ) herbicides were used for weed and barley control, in addition to hand-weeding. Carrots were sprayed twice per year with a total of 2.5 pounds of Kocide 3000 (46.1% copper hydroxide; Certis USA, Columbia, MD) to supply copper and suppress diseases. Additional fungicides and insecticides were used according to label instructions to prevent and control diseases and pests based on weather and scouting. Carrots were drip irrigated as needed with a double row of T-tape on each bed (30-cm emitter spacing).
Carrots were harvested by loosening the soil with a broadfork and digging a uniform 0.8-m-long section of all three planted rows from each sub-subplot (with beds centered at 1.5 m, 1.2 m2 were harvested). All taproot parts, including those that broke during harvest, were included for data collection. Taproots were washed immediately, then stored in a walk-in cooler at 0 to 2 °C for up to 1 week before data collection. Carrots were stored intact in 2019 and with foliage removed at the crown in 2020. Humidity was maintained by manually watering the floor of the cooler. ‘Cupar’ was harvested 3 Oct. 2019 (111 d) and 10 Sept. 2020 (104 d). ‘Envy’ was harvested on 16 to 20 Sept. 2019 (94–98 d) and 22 to 23 Sept. 2020 (116–117 d). ‘Naval’ was harvested on 25 to 26 Sept. 2019 (103–104 d) and 21 to 22 Sept. 2020 (115–116 d).
After removing all foliage, taproots were sized into the following categories, in order: diameter <1.9 cm (unmarketable); diameter ≥ 1.9 cm, but forked, bent, or cracked (unmarketable); diameter ≥2.5 cm and length ≥ 12.7 cm (“large” marketable); 1.9 cm ≤ diameter < 2.5 cm and length ≥ 12.7 cm (“narrow” marketable); diameter ≥ 2.5 cm and 7.6 cm ≤ length < 12.7 cm (“short” marketable); and finally, diameter ≥1.9 cm and length <7.6 cm (unmarketable). The number and weight of carrots in each category were recorded, and the total number and weight were calculated. “Large,” “narrow,” and “short” categories are similar to USDA topped carrot grades U.S. No. 1 Jumbo, U.S. Extra No. 1, and U.S. No. 2, respectively (USDA, 2020). Per-root fresh weight of large-size taproots was calculated from aggregate data.
Taproot length and maximum diameter were measured on a subsample of up to 10 randomly selected large-size taproots from each plot (nmin = 3, nmean = 9.2, nmedian = 10) using a custom-made jig. The number measured was fewer than 10 for plots where fewer than 10 large-size taproots were harvested. Cylindricity was determined as described by Bleasdale and Thompson (1963), by individually dividing taproot fresh weight (g) by the volume (cm3) of a theoretical cylinder with diameter and length of each taproot {g/[π × (cm radius)2 × (cm length)]}. Cylindricity is the percent of volume of the theoretical cylinder filled by the carrot (measured carrot density of 13 taproots was 1.001 ± 0.017 g⋅cm−3 by water displacement at room temperature). Individual length, diameter, and cylindricity values were averaged per plot before analysis.
Soil penetration resistance was measured on 24 July 2019 or 20 July 2020 using a soil penetrometer (AgraTronix model 08180; AgraTronix, Streetsboro, OH). Two sets of penetration resistance measurements per plot were recorded on the same day at 7.6, 15.2, 22.9, 30.5, and 38.1 cm. Weather data were collected on-site by National Weather Service station 218692 (Table 2, Fig. 1).
Temperature and growing degree days accumulated during the 2019 and 2020 carrot-growing seasons.
Statistics.
Two ‘Cupar’ plots were removed from the analysis because the results recorded were clearly inaccurate (>600,000 marketable taproots per ha and <100 g per large-size taproot, for example). Yield, size, and shape responses were subjected to mixed-model analysis of variance (ANOVA) using nlme::lme in R (Pinheiro et al., 2020; R Core Team, 2020). To accommodate the split-split plot design, the random intercept for all models was (1 | year/block/hybrid/seeding rate/bed). Fixed effects of hybrid (n = 3), seeding rate (n = 2), and bed (n = 2) were included as full-factorial effects. Year was included as a random blocking effect and not a fixed effect because differences due to growing season were not of interest. The three categories of marketable carrots (“large,” “narrow,” and “short”) were combined for marketable yield analyses, but “large” carrot yield (number and weight) was also analyzed independently. “Large” taproots as a percentage of total marketable taproots, and forked, bent, or cracked (unmarketable) taproots as a percent of total harvested taproots were analyzed also. Responses were transformed as needed to meet ANOVA assumptions (Tables 3 and 4). Interactions with Type II P values > 0.1 were removed to create reduced models, and the reduced models were compared with the full models using log-likelihood ratio tests (stats::anova; R Core Team, 2020).
Comparisons of a mixed-effect analysis of variance reduced model to a full-factorial model for each of the carrot taproot yield responses measured. Transformation used for each response is indicated, as well as Wald χ2 (Type II) P values for each of the fixed effects remaining in the reduced model: carrot hybrid (n = 3), seeding rate (sr, n = 2), bed shape (n = 2), and interactions. Akaike’s information criterion (AIC) is reported for full and reduced models. Factors with P < 0.05 are highlighted in bold. Common Language effects size (CL) of raised beds or high seeding rate is reported. A log-likelihood ratio test statistic (Λdf), comparing the reduced to the full model, is reported with P value and df. Misshapen taproots include those that were forked, bent, or cracked.
Comparisons of a mixed-effect analysis of variance reduced model to a full-factorial model for each of the carrot taproot shape and size responses measured. A log-likelihood ratio test statistic (Λdf), comparing the reduced to the full model, is reported with P value and df. Akaike’s information criterion (AIC) is reported for full and reduced models. Measurements were made on up to 10 large-size marketable carrots per plot, and arithmetic per-plot mean was analyzed. Cylindricity is g per taproot as a fraction of maximum possible taproot volume (cm3) per carrot. Transformation used for each response is indicated, as well as Wald χ2 (Type II) P values for each of the fixed effects remaining in the reduced model: carrot hybrid (n = 3), seeding rate (sr; n = 2), bed shape (n = 2), and interactions. Factors with P < 0.05 are highlighted in bold. Common Language effects size (CL) of raised beds or high seeding rate is reported.
To compare hybrids, Tukey’s honestly significant difference comparisons were made between predicted means [emmeans::emmeans, Lenth (2020)]. Contrasts between predicted treatment means were also made using emmeans for seeding rate, bed, and significant (P < 0.05) interaction effects. Confidence intervals (CI, ±95%) of the difference between the treatment level (raised bed or high seeding rate) and the reference level (flat bed or low seeding rate) were back-calculated from transformed predicted differences and standard errors of the difference [CI = predicted mean difference ± 1.96 × (standard error of the difference)].
Soil penetration resistance at 7.6 and 15.2 cm was analyzed with mixed-model ANOVA, using bed type (flat or raised) and measurement depth as fixed effects. With two sets of measurements per plot, a within-plot random effect was included in addition to previously described random effects. Residual assumptions could not be satisfied for penetration resistance data beyond 15.2 cm, so three separate Wilcoxon signed-rank tests were performed to test for soil penetration resistance on flat vs. raised beds at 22.9, 30.5, and 38.1 cm (stats::wilcox.test); the two penetration resistance values at each depth from each sub-subplot were averaged, and bed types (sub-subplots) within seeding rate were considered paired. P values from the three tests were Holm-adjusted (stats::p.adjust).
To study direct effects of soil penetration resistance on parameters measured, simple linear models were made regressing normalized responses against penetration resistance at 7.6 cm (stats::lm). Responses were normalized before regression by dividing individual data by the maximum response within hybrid. Pearson correlation coefficients (r) were calculated for estimating response to penetration resistance as well (stats::cor.test).
Results and Discussion
Precipitation was within expected ranges for much of both years (Fig. 1). However, in 2019, two separate 6-day periods with more than 100 mm were recorded. And 4 (nonconsecutive) days in 2020 experienced 39 to 112 mm of rainfall. Temperature and growing degree days from each year are reported in Table 2.
In reporting results from two-level factors (seeding rate and bed shape), we report the mean of the reference level (either low seeding rate or flat bed) and the CI of the difference from the reference level. This is done to illustrate the magnitude and direction of the effect, not just the distribution of the measured responses (raw standard error or CI) at each level of the factor. In addition, a Common Language effects size (CL) is reported for seeding rate and bed effects in Tables 3 and 4. The CL for the effect of high seeding rate on the number of marketable taproots per ha, for example, is 66% (Table 3). This means that a randomly selected sample from the high seeding rate treatment has a 66% chance of having more marketable taproots than a randomly selected sample from the low seeding rate treatment. Similarly, a randomly selected sample from the raised bed treatment has a 34% chance of being larger-diameter than a randomly selected sample from the flat bed treatment (Table 4). The CI of the difference and CL are measures that improve understanding of how meaningful the responses are in the current study, and are tools for understanding this research in the context of similar work (Kramer et al., 2019; Nakagawa and Cuthill, 2007).
A mixed-effect ANOVA model with only hybrid and seeding rate as fixed main effects explained variation in the total number of taproots harvested (marketable and unmarketable) equally as well as a full-factorial model that included bed shape (Λ8 = 7.34, P = 0.50). The reduced model was then compared with a model in which ‘Envy’ and ‘Naval’ results were grouped together, because they were planted at the same seeding rates. This model (reduced and grouped) performed equally as well as the reduced model alone (Λ1 = 0.02, P = 0.88), and AIC (Akaike information criterion) was improved at each step. Compared with ‘Envy’ and ‘Naval’ (928,776 per ha), the relative proportion of total carrots harvested in ‘Cupar’ was 38% to 45% (95% CI); the target ratio was 44% at seeding. At high seeding rate, 23% to 43% more total carrots were harvested than at low seeding rate (517,792 per ha); the target ratio was 35%. Though the relative amounts of harvested plants were within the targeted ranges, the absolute number of harvested taproots was far lower than what was planted. Stand counts and germination were not recorded, so the mechanism of this loss is unknown.
The overall average marketable yield was 34.5 Mg⋅ha−1, which was less than the 54.5 Mg⋅ha−1 reported from Wisconsin in 2019 (USDA NASS, 2020) and less than the yield goal of 44.9 Mg⋅ha−1 described in University of Minnesota fertility recommendations (Rosen and Eliason, 2005). However, carrots in this study were sorted according to fresh-market standards, not processing standards. Regional yield estimates for small-scale fresh-market production are as low as 10 to 15 Mg⋅ha−1 (Delahaut and Newenhouse, 2011), and the national average for fresh-market yield from 2010 to 2015 was 34.3 to 39.5 Mg⋅ha−1 (USDA NASS, 2020). We therefore realized a reasonable fresh-market yield. The number of marketable taproots and large-size taproots was greater in ‘Envy’ and ‘Naval’ compared with ‘Cupar’, but there were no differences in marketable Mg yield among the hybrids (Table 3, Fig. 2) because of the greater large-size taproot diameter and weight of ‘Cupar’ (Table 4, Fig. 3). A greater percentage of taproots were unmarketable in ‘Envy’ and ‘Naval’ (Fig. 2).
The largest impacts of raised beds on yield responses were related to the number and total weight of large-size taproots. There were 10,984 to 54,656 (95% CI) more large-size taproots per ha from raised beds (CL = 65%; Table 3, Fig. 4), and 2.0 to 7.3 (95% CI) more Mg⋅ha−1 of large-size taproots from raised beds (CL = 67%; Table 3, Fig. 4). The increase in number of large-size taproots due to raised beds was largely driven by ‘Naval’ (interaction P = 0.053; Table 3). ‘Cupar’ realized an average gain of 14% more large-size taproots due to raised beds, ‘Envy’ had 4% more, but ‘Naval’ had 41% more large-size taproots on raised beds compared with flat beds (not shown). An increase due to raised beds of 10.6% in total fresh weight of marketable taproots of all sizes was measured (CL = 62%; Table 3, Fig. 4), but this was less substantial than the increase in large-size taproots. There was no impact of raised beds on the total number of marketable taproots (CL = 50%; Table 3, Fig. 4).
Raised beds did not impact average length of large-size taproots (CL = 55%; Table 4, Fig. 4); however, raised beds increased the proportion of large-size taproots (95% CI = 4.7%–13.4% greater percentage of large-size taproots from raised beds, CL = 68%; Table 3, Fig. 5). The fact that the total number of marketable taproots did not increase due to raised beds, but weight of marketable taproots and total weight and number of large-size taproots did increase, led to a reasonable conclusion that “short” marketable taproots (76–127 mm) were able to achieve greater length (≥127 mm) on raised beds compared with flat beds, thus explaining the increase in proportion of “large” taproots. The fact that an increase in length of large-size taproots was not detected from raised beds obscures the increase in overall average length of marketable taproots.
The average fresh weight of large-size taproots grown on flat beds at high seeding rate was lower than fresh weight on flat beds at low seeding rate, but the seeding rate impact was not as great on raised beds (interaction P = 0.004; Table 4, Fig. 6). At low seeding rate, large-size taproot weight was not changed by growing carrots on raised vs. flat beds (Fig. 4). If one is interested in greater fresh weight per large-size taproot, high seeding rate, especially without raised beds, should be avoided. Although 23% to 43% more total taproots were harvested under high seeding rate, only 5% to 30% more marketable taproots were harvested with high seeding rate (CL = 66%; Table 3, Fig. 6). There was no less total marketable weight when seeding rate was high compared with low (CL = 46%; Table 3, Fig. 6), we found more total marketable taproots due to high seeding rate, and there were no differences in the number of large-size taproots due to high seeding rate compared with low seeding rate (CL = 42%; Table 3, Fig. 6). However, a lower percentage of marketable carrots were “large” size when grown at high seeding rate compared with low seeding rate (CL = 21%; Table 3, Fig. 5). We therefore conclude that high seeding rate caused smaller, more numerous marketable taproots. This is also supported by the smaller diameter (CL = 28%) and shorter length (CL = 25%) of large-size taproots at high seeding rate (Table 4, Fig. 6). Previous research has shown similar results. Taproot length was reduced under elevated population density in ‘Amsterdam Forcing’ and ‘Autumn King’ (Thompson, 1969).
It is possible that a limit to length of carrots in our soil was reached, as the maximum plot average of large-size carrots in our study was 195 mm. This was reported in a ‘Cupar’ plot that had the fewest total roots harvested (Fig. 3). There was no difference in length between large-size ‘Cupar’ and ‘Envy’ (Fig. 3), although ‘Cupar’ is reported by the breeder to grow to ≈200 mm and ‘Envy’ is reported to grow to ≈280 mm. ‘Cupar’ did have a lower percentage of large-size taproots compared with ‘Envy’ (Fig. 5), but perhaps growing seasons were not long enough in this study (94–116 d, 1454–1808 GDD5; growing degree days, 5 °C base temperature), especially for ‘Envy’ to reach its full potential. ‘Envy’ is described by the breeder as an early- to midseason hybrid, and ‘Cupar’ (104–111 d, 1671–1691 GDD5) is described as midlate. ‘Naval’ (103–115 d, 1599–1794 GDD5) is also described as midlate and ≈229 mm, but large-size ‘Naval’ were shorter in this study than either ‘Envy’ or ‘Cupar’ (Table 4, Fig. 3).
Overall soil penetration resistance in this study was relatively low at all depths and in both bed types. Eighty-nine percent of soil penetration resistance measurements at 15.2 cm depth were <1 MPa, and average penetration resistance beyond 15.2 cm was 0.77 MPa (not shown). In one study, penetration resistance in noncompacted, 30% clay soil was 1.4 MPa at 7.0 to 17.5-cm depth (Pietola, 1995). In that study, 158-mm taproots were grown with no imposed compaction treatment (cv. Nantes Duke Notabene 370 Sv), which is similar to lengths of large-size carrots recorded in our study (Fig. 3).
Penetration resistance at 7.6 and 15.2 cm soil depth was lower on raised beds compared with flat beds (interaction P < 0.001; Fig. 7). At 7.6 cm, penetration resistance was 33% lower on raised beds compared with flat beds (paired CL for flat-raised = 20%, n = 72). At 15.2 cm, penetration resistance was 9.8% lower (CL = 38%). Wilcoxon signed-rank tests from deeper penetrometer readings revealed no penetration resistance difference between bed shapes at these depths (not shown).
To explore soil compaction as a mechanism of the bed shape effects, yield and shape responses were regressed against penetration resistance at 7.6 cm soil depth. First, responses were normalized within hybrid by dividing each data point by the maximum value within each hybrid. This was done because of significant hybrid effects for some responses, and a correlation coefficient and slope independent of hybrid was desirable. Then, due to correlation between numerous yield responses (Fig. 8), a “combined yield response” was created. This was the average of the normalized responses for Mg marketable taproots per ha, n and Mg large-size taproots per ha, percentage of taproots that were marketable, and (100 − percent of taproots that were misshapen). These responses were all highly positively correlated (0.63 ≤ r ≤ 0.93; Fig. 8).
The combined yield response showed a substantial decrease as soil penetration resistance at 7.6 cm increased (Fig. 9). The number of taproots that were marketable showed a similar, but less severe response (shallower slope) to soil penetration resistance (Fig. 9). Maximizing carrot yield clearly requires limiting compaction, at least in the top 7.6 cm of soil, which is supported by previous research (Johansen et al., 2015). Reducing compaction up to 15.2 cm by creating raised beds is a reasonable mechanism for increasing carrot yield. What is unclear in this study is how raised beds were unable to decrease the percentage of taproots that were forked, bent, or cracked (CL = 49%, P = 0.83; Table 3, Fig. 4) even though soil penetration resistance was positively correlated with the normalized percentage of misshapen taproots (P < 0.001, r = 0.39; data not shown), and compaction substantially increased percentage of forked and split roots in a previous study (Pietola, 1995).
Variation in cylindricity of large-size taproots was largely determined by hybrid and bed shape Table 4). Cylindricity was 61.6% on flat beds and 66.3% on raised beds (CL = 87%; Table 4, Fig. 4). The effect of raised beds on cylindricity was slightly greater in ‘Envy’ and ‘Naval’ compared with ‘Cupar’ (interaction P = 0.081; Table 4), and ‘Naval’ was far more cylindrical on average than ‘Cupar’ or ‘Envy’ (Fig. 10). Nantes-type carrots, like ‘Naval’, are identified by their cylindrical shape and blunt taproot tip. Seeding rate had no effect on cylindricity within ‘Cupar’, and there was a trend for contrasting seeding rate effects in ‘Envy’ and ‘Naval’ (interaction P = 0.039; Table 4, Fig. 10). Cylindricity of large-size taproots decreased slightly as penetration resistance increased, but the correlation was robust (Fig. 9). This response has been noted in previous research (Johansen et al., 2015). The decrease in large-size taproot cylindricity associated with soil penetration resistance and raised beds could be explained by a soil compaction-induced increase in diameter with no change in length or fresh weight of large-size taproots (Table 4, Figs. 4 and 9). We did not study organoleptic properties or chemical constituents of carrots, but one study found compaction-induced reduction in cylindricity was related to reduced quality of color and texture (Taksdal, 1984).
The Nicollet soil used in this study is classified as fine clay loam. Clay loam soils are less compressible than sandy or silty loam soils (Warren and Taylor, 2014) that are commonly used for carrot production in the Midwest. It is possible that more extreme effects of raised beds might be seen on more compressible soils, or on taller raised beds than were made here. However, one previous study found compaction treatments had a more detrimental effect on carrot yield and size in clay soil than in a sandy soil (Pietola, 1995). Soil moisture characteristics were not studied here, but higher water-holding capacity of the clay-loam soil also may have reduced yield compared with coarse-texture or peat/muck soil. A greenhouse study using muck soil found more available water caused shorter carrot length and reduced yield compared with a moderate amount of available water (White, 1992).
The larger diameter and greater taper (lower cylindricity) of ‘Cupar’ made this hybrid easier to harvest by hand compared with ‘Envy’ or ‘Naval’ (personal observation). Soil conditions for harvest in the clay-loam soil were not always ideal (especially in very wet or dry periods), and ‘Cupar’ could sometimes be harvested simply by pulling taproots from loosened soil by their tops. This was not typically the case with ‘Envy’ or ‘Naval’. The latter two tended to break below ground if the soil was not loosened deeply enough, which was less of an issue with ‘Cupar’. This benefit likely extends to other Chantenay-type carrots, with less-cylindrical and larger-diameter taproots. And if this anecdotal evidence holds true, it may be that carrots grown at lower seeding rate are easier to hand-harvest due to reduced cylindricity and increased diameter, with no effects on total yield except for number of marketable roots. However, lower seeding rates may create challenges with weed control (Colquhoun et al., 2020). Increasing the number of rows per bed, rather than plants per row, may be an optimum strategy for weed control and maintenance of taproot size and yield (Colquhoun et al., 2020).
In conclusion, raised beds are an effective measure to increase total yield of carrots in clay-loam soil, and cultivar selection can be a useful tool for maximizing efficiency of producing large-size fresh-market carrots, especially with cultivar-appropriate seeding rates.
Literature Cited
Bakker, D.M., Hamilton, G.J., Houlbrooke, D.J. & Spann, C. 2005 The effect of raised beds on soil structure, waterlogging, and productivity on duplex soils in Western Australia Austral. J. Soil Res. 43 575 585 doi: 10.1071/SR03118
Bleasdale, J.K.A. & Thompson, R. 1963 An objective method of recording and comparing the shapes of carrot root J. Hort. Sci. 38 232 241 doi: 10.1080/00221589.1963.11514074
Brown, P. & Gracie, A. 2000 Factors influencing carrot size and shape. Publication VG97019 Horticulture Austral. Limited Sydney, NSW
Colquhoun, J.G., Rittmeyer, R.A. & Heider, D.J. 2020 A holistic carrot production system for season-long weed management Weed Technol. 34 876 881 doi: 10.1017/wet.2020.84
Delahaut, K.A. & Newenhouse, A.C. 2011 Growing carrots, beets, radishes, and other root crops in Wisconsin: A guide for fresh-market growers UW-Extension publication A3686
Jauron, R. & Wallace, G. 2017 Yard and garden: Planting carrots in the home garden
Johansen, T.J., Thomsen, M.G., Løes, A.-K. & Riley, H. 2015 Root development in potato and carrot crops — influences of soil compaction Acta Ag. Scandinavica B. 65 2 182 192 doi: 10.1080/09064710.2014.977942
Kelley, W.T. & Phatak, S.C. 2009 Commercial production and management of carrots Univ. Ga. Ext. Bul. 1175
Kramer, M.H., Paparozzi, E.T. & Stroup, W.W. 2019 Best practices for presenting statistical information in a research article HortScience 54 9 1605 1609 doi: 10.21273/HORTSCI13952-19
Lakens, D. 2013 Calculating and reporting effect sizes to facilitate cumulative science: A practical primer for t-tests and ANOVAs Front. Psychol. 4 1 12 doi: 10.3389/fpsyg.2013.00863
Lenth, R.V. 2020 Emmeans: Estimated marginal means, aka least-squares means
McGraw, K.O. & Wong, S.P. 1992 A common language effect size statistic Psychol. Bull. 111 2 361 365 doi: 10.1037/0033-2909.111.2.361
Nakagawa, S. & Cuthill, I.C. 2007 Effect size, confidence interval and statistical significance: A practical guide for biologists Biol. Rev. Camb. Philos. Soc. 82 591 605 doi: 10.1111/j.1469-185X.2007.00027.x
Pietola, L. 1995 Effect of soil compactness on the growth and quality of carrots Agr. Sci. Finland 4 144 237 doi: 10.23986/afsci.72611
Pinheiro, J., Bates, D., DebRoy, S. & Sarkar, D. R Core Team 2020 nlme: Linear and nonlinear mixed effects models
Rosen, C.J. & Eliason, R. 2005 Nutrient management for commercial fruit and vegetable crops in Minnesota
R Core Team 2020 R: A language and environment for statistical computing
Smetanka, M.J 2017 Eager to get growing? Many root vegetables can be sown in April
Taksdal, G. 1984 Effect of tractor wheelings on carrot quality Acta Hort. 163 255 260
Thompson, C 2017 Carrots provide benefits for small scale growers
Thompson, R. 1969 Some factors affecting carrot root shape and size Euphytica 18 277 285
USDA 2020 United States standards for grades of topped carrots
USDA NASS 2020 NASS-Quick Stats
Warren, J. & Taylor, R. 2014 Managing soil compaction
White, J.M. 1992 Carrot yield when grown under three soil water concentrations HortScience 27 2 105 106 doi: 10.21273/HORTSCI.27.2.105