## 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 H_{2}O (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 m^{2} 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 (n_{min} = 3, n_{mean} = 9.2, n_{median} = 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 (cm^{3}) 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 (cm^{3}) 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)].

*Z*-score was calculated from the difference of raw means and a correlation-adjusted pooled standard deviation:

*s*is the raw sample standard deviation, and

*r*is the Pearson’s correlation coefficient between paired bed types within a seeding rate subplot. Only paired data within subplots were used (n = 68 pairs used due to previously described data removal). The cumulative probability distribution greater than

*Z*in a unit normal distribution was calculated using Excel software [=norm.dist(

*Z*,0,1)] to determine CL. CL was similarly calculated for the effect of high seeding rate, after averaging sub-subplot (bed) data within each seeding rate and removing incomplete pairs of subplots (n = 34 pairs used).

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).

Carrot taproot yield parameters measured in three hybrids: ‘Cupar’, ‘Envy’, and ‘Naval’. The number and total weight of marketable taproots, number and total weight of large-size taproots (≥ 127 mm long, ≥ 25 mm diameter), and the percentage of all taproots that were marketable or misshapen (forked, bent, or cracked) were measured. Violin plots illustrate distribution of data, data points indicate predicted mean, and bars represent ±se of predicted mean. Data with different letters indicate statistical differences based on Tukey’s honestly significant difference (α = 0.05), and lack of letters indicates lack of statistical significance.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Carrot taproot yield parameters measured in three hybrids: ‘Cupar’, ‘Envy’, and ‘Naval’. The number and total weight of marketable taproots, number and total weight of large-size taproots (≥ 127 mm long, ≥ 25 mm diameter), and the percentage of all taproots that were marketable or misshapen (forked, bent, or cracked) were measured. Violin plots illustrate distribution of data, data points indicate predicted mean, and bars represent ±se of predicted mean. Data with different letters indicate statistical differences based on Tukey’s honestly significant difference (α = 0.05), and lack of letters indicates lack of statistical significance.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Carrot taproot yield parameters measured in three hybrids: ‘Cupar’, ‘Envy’, and ‘Naval’. The number and total weight of marketable taproots, number and total weight of large-size taproots (≥ 127 mm long, ≥ 25 mm diameter), and the percentage of all taproots that were marketable or misshapen (forked, bent, or cracked) were measured. Violin plots illustrate distribution of data, data points indicate predicted mean, and bars represent ±se of predicted mean. Data with different letters indicate statistical differences based on Tukey’s honestly significant difference (α = 0.05), and lack of letters indicates lack of statistical significance.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Size of large-size carrot taproots (≥ 127 mm long, ≥ 25 mm diameter) measured in three hybrids: ‘Cupar’, ‘Envy’, and ‘Naval’. Violin plots illustrate distribution of data, data points indicate predicted mean, and bars represent ±se of predicted mean. Bars may be obscured by the mean. Data with different letters indicate statistical differences based on Tukey’s honestly significant difference (α = 0.05), and lack of letters indicates lack of statistical significance.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Size of large-size carrot taproots (≥ 127 mm long, ≥ 25 mm diameter) measured in three hybrids: ‘Cupar’, ‘Envy’, and ‘Naval’. Violin plots illustrate distribution of data, data points indicate predicted mean, and bars represent ±se of predicted mean. Bars may be obscured by the mean. Data with different letters indicate statistical differences based on Tukey’s honestly significant difference (α = 0.05), and lack of letters indicates lack of statistical significance.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Size of large-size carrot taproots (≥ 127 mm long, ≥ 25 mm diameter) measured in three hybrids: ‘Cupar’, ‘Envy’, and ‘Naval’. Violin plots illustrate distribution of data, data points indicate predicted mean, and bars represent ±se of predicted mean. Bars may be obscured by the mean. Data with different letters indicate statistical differences based on Tukey’s honestly significant difference (α = 0.05), and lack of letters indicates lack of statistical significance.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

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).

Effect size of raised beds on carrot taproot yield and size parameters; 100% (vertical yellow line) represents the predicted mean of flat beds, circles represent predicted mean difference of carrots grown on raised beds, and error bars are the ±95% confidence intervals of the difference from flat beds. Treatments with error bars that do not overlap the flat bed mean indicate no difference between flat and raised beds (closed circles, contrast *P* > 0.05). Flat bed means for each parameter are indicated parenthetically. Large-size taproots were ≥ 127 mm long, ≥ 25 mm diameter. Cylindricity is the volume filled by the taproot compared with a cylinder of the same diameter and length as the taproot.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Effect size of raised beds on carrot taproot yield and size parameters; 100% (vertical yellow line) represents the predicted mean of flat beds, circles represent predicted mean difference of carrots grown on raised beds, and error bars are the ±95% confidence intervals of the difference from flat beds. Treatments with error bars that do not overlap the flat bed mean indicate no difference between flat and raised beds (closed circles, contrast *P* > 0.05). Flat bed means for each parameter are indicated parenthetically. Large-size taproots were ≥ 127 mm long, ≥ 25 mm diameter. Cylindricity is the volume filled by the taproot compared with a cylinder of the same diameter and length as the taproot.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Effect size of raised beds on carrot taproot yield and size parameters; 100% (vertical yellow line) represents the predicted mean of flat beds, circles represent predicted mean difference of carrots grown on raised beds, and error bars are the ±95% confidence intervals of the difference from flat beds. Treatments with error bars that do not overlap the flat bed mean indicate no difference between flat and raised beds (closed circles, contrast *P* > 0.05). Flat bed means for each parameter are indicated parenthetically. Large-size taproots were ≥ 127 mm long, ≥ 25 mm diameter. Cylindricity is the volume filled by the taproot compared with a cylinder of the same diameter and length as the taproot.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

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.

Percentage of harvested carrot taproots in each marketable size category depending on hybrid, seeding rate, and bed shape. Hybrids with the same letter do not have a different percentage of large-size taproots based on Tukey’s honestly significant difference (α = 0.05). ***Treatment effect on the percentage of large-size taproots is different from the reference level (*P* < 0.001).

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Percentage of harvested carrot taproots in each marketable size category depending on hybrid, seeding rate, and bed shape. Hybrids with the same letter do not have a different percentage of large-size taproots based on Tukey’s honestly significant difference (α = 0.05). ***Treatment effect on the percentage of large-size taproots is different from the reference level (*P* < 0.001).

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Percentage of harvested carrot taproots in each marketable size category depending on hybrid, seeding rate, and bed shape. Hybrids with the same letter do not have a different percentage of large-size taproots based on Tukey’s honestly significant difference (α = 0.05). ***Treatment effect on the percentage of large-size taproots is different from the reference level (*P* < 0.001).

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

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).

Effect size of high seeding rate on carrot taproot yield and size parameters; 100% (vertical yellow line) represents the predicted mean of low relative seeding rate, circles represent predicted mean difference of high seeding rate, and error bars are the ±95% confidence interval of the difference from low seeding rate. Treatments with error bars that do not overlap the low seeding rate mean indicate no difference between low and high seeding rate (closed circles, contrast *P* > 0.05). Low seeding rate means for each parameter are indicated parenthetically. Large-size taproots were ≥ 127 mm long, ≥ 25 mm diameter.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Effect size of high seeding rate on carrot taproot yield and size parameters; 100% (vertical yellow line) represents the predicted mean of low relative seeding rate, circles represent predicted mean difference of high seeding rate, and error bars are the ±95% confidence interval of the difference from low seeding rate. Treatments with error bars that do not overlap the low seeding rate mean indicate no difference between low and high seeding rate (closed circles, contrast *P* > 0.05). Low seeding rate means for each parameter are indicated parenthetically. Large-size taproots were ≥ 127 mm long, ≥ 25 mm diameter.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Effect size of high seeding rate on carrot taproot yield and size parameters; 100% (vertical yellow line) represents the predicted mean of low relative seeding rate, circles represent predicted mean difference of high seeding rate, and error bars are the ±95% confidence interval of the difference from low seeding rate. Treatments with error bars that do not overlap the low seeding rate mean indicate no difference between low and high seeding rate (closed circles, contrast *P* > 0.05). Low seeding rate means for each parameter are indicated parenthetically. Large-size taproots were ≥ 127 mm long, ≥ 25 mm diameter.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

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 GDD_{5}; 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 GDD_{5}) is described as midlate. ‘Naval’ (103–115 d, 1599–1794 GDD_{5}) 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).

Soil penetration resistance at 7.6 or 15.2 cm from the surface, measured by a soil penetrometer. Soil was either flat beds or ≈7.6-cm raised beds. Data are from July 2019 and 2020. Violin plots illustrate distribution of data, circles indicate predicted mean, and error bars represent the 95% confidence intervals of the difference between bed types. Within-depth significance level (*P* value) of contrasts between predicted means of each bed type is shown.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Soil penetration resistance at 7.6 or 15.2 cm from the surface, measured by a soil penetrometer. Soil was either flat beds or ≈7.6-cm raised beds. Data are from July 2019 and 2020. Violin plots illustrate distribution of data, circles indicate predicted mean, and error bars represent the 95% confidence intervals of the difference between bed types. Within-depth significance level (*P* value) of contrasts between predicted means of each bed type is shown.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Soil penetration resistance at 7.6 or 15.2 cm from the surface, measured by a soil penetrometer. Soil was either flat beds or ≈7.6-cm raised beds. Data are from July 2019 and 2020. Violin plots illustrate distribution of data, circles indicate predicted mean, and error bars represent the 95% confidence intervals of the difference between bed types. Within-depth significance level (*P* value) of contrasts between predicted means of each bed type is shown.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

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).

Pearson’s correlation between measured carrot responses. Correlation coefficient (*r*) value is displayed, and also represented by color and circle size. Only correlations with *P* > 0.05 are shown.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Pearson’s correlation between measured carrot responses. Correlation coefficient (*r*) value is displayed, and also represented by color and circle size. Only correlations with *P* > 0.05 are shown.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Pearson’s correlation between measured carrot responses. Correlation coefficient (*r*) value is displayed, and also represented by color and circle size. Only correlations with *P* > 0.05 are shown.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

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).

Response of carrot taproot yield, size, and shape parameters to soil penetration resistance (MPa) at 7.6-cm depth. All values were normalized to % of maximum value within each hybrid before analysis. Combined yield response is the per-plot average of 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), all normalized within hybrid. The width of the red regression line is proportional to (1 − *P* value) for the simple linear model. The 95% confidence interval (CI) for the correlation between penetration resistance and response is reported (*r*), and regression slope (m) 95% CI is reported for regressions with *P* < 0.01.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Response of carrot taproot yield, size, and shape parameters to soil penetration resistance (MPa) at 7.6-cm depth. All values were normalized to % of maximum value within each hybrid before analysis. Combined yield response is the per-plot average of 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), all normalized within hybrid. The width of the red regression line is proportional to (1 − *P* value) for the simple linear model. The 95% confidence interval (CI) for the correlation between penetration resistance and response is reported (*r*), and regression slope (m) 95% CI is reported for regressions with *P* < 0.01.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Response of carrot taproot yield, size, and shape parameters to soil penetration resistance (MPa) at 7.6-cm depth. All values were normalized to % of maximum value within each hybrid before analysis. Combined yield response is the per-plot average of 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), all normalized within hybrid. The width of the red regression line is proportional to (1 − *P* value) for the simple linear model. The 95% confidence interval (CI) for the correlation between penetration resistance and response is reported (*r*), and regression slope (m) 95% CI is reported for regressions with *P* < 0.01.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

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).

Carrot seeding rate effects on taproot cylindricity in three carrot hybrids. Cylindricity is the volume filled by the taproot compared with a cylinder of the same diameter and length as the taproot. Violin plots illustrate distribution of data, circles indicate predicted mean, and error bars represent the 95% confidence interval of the difference between seeding rates within hybrid. Within-hybrid significance level (*P* value) of contrasts between predicted means of each bed type is shown.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Carrot seeding rate effects on taproot cylindricity in three carrot hybrids. Cylindricity is the volume filled by the taproot compared with a cylinder of the same diameter and length as the taproot. Violin plots illustrate distribution of data, circles indicate predicted mean, and error bars represent the 95% confidence interval of the difference between seeding rates within hybrid. Within-hybrid significance level (*P* value) of contrasts between predicted means of each bed type is shown.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

Carrot seeding rate effects on taproot cylindricity in three carrot hybrids. Cylindricity is the volume filled by the taproot compared with a cylinder of the same diameter and length as the taproot. Violin plots illustrate distribution of data, circles indicate predicted mean, and error bars represent the 95% confidence interval of the difference between seeding rates within hybrid. Within-hybrid significance level (*P* value) of contrasts between predicted means of each bed type is shown.

Citation: HortScience horts 56, 6; 10.21273/HORTSCI15823-21

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/SR03118Bleasdale, 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.11514074Brown, 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.84Delahaut, 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.977942Kelley, W.T. & Phatak, S.C. 2009 Commercial production and management of carrots

*Univ. Ga. Ext. Bul.*1175Kramer, 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-19Lakens, 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.00863Lenth, 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.361Nakagawa, 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.xPietola, L. 1995 Effect of soil compactness on the growth and quality of carrots

*Agr. Sci. Finland*4 144 237 doi: 10.23986/afsci.72611Pinheiro, 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 260Thompson, C 2017 Carrots provide benefits for small scale growers

Thompson, R. 1969 Some factors affecting carrot root shape and size

*Euphytica*18 277 285USDA 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