DOI: 10.5281/zenodo.6079943
title: Raw data for "Assessing Cover Crop and Intercrop Performance Along a Farm Management Gradient" (2022)

description: This dataset accompanies the publication "Assessing Cover Crop and Intercrop Performance Along a Farm Management Gradient" by Stratton et al. in the journal Agriculture, Ecosystems, and Environment (2022). https://doi.org/10.1016/j.agee.2022.107925


METHODS:


We conducted our experiment between May 2018 and December 2019 on 14 farms in the eastern coastal highlands region of Santa Catarina, Brazil. The mean altitude of sites was 467 m (+/- 161 m).  Eastern Santa Catarina has a subtropical climatic pattern, with mean annual rainfall ranging from 1,500-1,700 mm (Wrege et al., 2012). While 2018 had typical weather patterns for the region, 2019 was a dry year, particularly during the spring months (Appendix B, Table B.1). All farms were located in the Colonial Serrana Catarinense soil microregion, one of 16 designated microregions in the state of Santa Catarina (EMBRAPA, 2004). Primary soil types in our study site are associations of dystric Cambisols and haplic Acrisols (typic Dystrocryepts and typic Paleudults in the USDA Soil Taxonomy), which tend to be moderately to highly acidic, with limited soil nutrient availability and moisture retention (EMBRAPA, 2004; IUSS Working Group WRB, 2015; USDA, 2010). To support crop production, farmers in the region typically apply lime (calcium and magnesium carbonate) to agricultural fields to increase soil pH from <5.5 to 6 (Comissão de Química e Fertilidade do Solo - RS/SC, 2016). Exact farm locations within the region are not given and farmer identities have been anonymized.


Experimental design


The fully factorial experiment had six treatments (Figure 2): (1) cover crop + pea-cucumber intercrop, (2) cover crop + pea monocrop, (3) cover crop + cucumber monocrop, (4) fallow + pea-cucumber intercrop, (5) fallow + pea monocrop, and (6) fallow + cucumber monocrop. Due to the timing of farm recruitment, only conventional and transitioning farms participated in the first year of cover cropping (2018); agroecological farms were added to the study during the vegetable intercropping period of 2018 and had their first round of cover cropping in 2019. The cover crop mixture treatment was designed to emulate traditional practices in the region, as well as to include functionally complementary legume and grass species: common vetch (Vicia sativa L.) and black oat (Avena strigosa Schreb). We also selected vegetables with distinct ecological functional traits, such that intercropping represented an increase in functional diversity relative to mono-cropped vegetables. Snow peas are N-fixing legumes with a vining, upright structure and a deep root system, whereas cucumbers are low-lying, non-legume cucurbits that provide groundcover and have a relatively shallow, extensive root system.


 


Cover crop treatments consisted of two adjacent 50 m2 plots in each field, one of which was planted with the cover crop mixture; the other served as a weedy fallow control. In 2018, the cover crop mixture seeding rate was 72 kg/ha black oat and 36 kg/ha common vetch. Due to poor vetch performance in mixtures at this rate, we increased the vetch seeding rate to 60 kg/ha in 2019, maintaining the black oat rate from 2018. Cover crop seeds were inoculated with the Brazilian strain Rhizobium etli (SEMIA 384; source: FEPAGRO) at 4 g/kg vetch seed prior to planting. Cover crops were grown until peak flowering, and then cover crops (and weeds in the fallow) were incorporated into the soil by rototiller (n = 7 farms) or by hand hoeing (n = 7 farms), based on farms’ available machinery, between September 5-10 in 2018 and September 10-18 in 2019 (approximately one week following cover crop sampling on each farm).


 


Vegetables were planted two weeks following cover crop and weed biomass incorporation within a period of 7-10 days across sites. Harvest dates were spaced such that crops were growing for approximately the same period across farms. The 50 m2 plots were each divided into three intercrop treatments with a ~1 m2 pathway between each treatment, for a total of 6 treatments randomly assigned to plots per 100 m2. We planted a climbing variety of snow peas (Pisum sativum subsp. sativum var. macrocarpum, “Torta de flor roxa”) and pickling cucumber (Cucumis sativa L. var. Pepino HT 05) in intercrops and in their respective monocrops, using a replacement design (i.e., equivalent crop densities in all treatments). Snow pea seeds were inoculated with Rhizobium leguminosarum var. viceae (SEMIA 3007/BR 619, source: UFSC ENR/CCA) at a rate of 4 g/kg directly prior to planting. There were five rows of crops per treatment, with only the three middle rows harvested to limit edge effects. In-row spacing was 60 cm for cucumber and 20 cm for peas, with 60 cm between rows in both intercrops and monocrops. Cucumbers were grown as starts for 2.5 weeks before planting, and peas were planted from seed on the same planting date as cucumber starts.


 


In the summer between January and May 2019 all fields were planted to a sunflower (Helianthus annuus L.) crop, which was incorporated into the soil during flowering approximately two weeks prior to cover crop planting in 2019. Because we sought to understand the effects of crop diversification given existing water and nutrient limitations on working farms, the experiment was entirely rainfed and legume N fixation was the sole external N source.


Soil sampling and analysis


Prior to the first cover cropping period, we collected a composite sample of 15-20 soil cores (2.5 cm diameter, 20 cm depth) on both the cover crop and fallow sides of each experimental field (n = 28) for analysis of baseline conditions (see Appendix B for full details). Briefly, soil was analyzed for pH, macro- and micronutrients, and soil organic matter (SOM) by the Santa Catarina State Agricultural Agency (EPAGRI) in Ituporanga, Santa Catarina, Brazil, using standard protocols (Comissão de Química e Fertilidade do Solo - RS/SC, 2016). pH was measured with a glass electrode both with and without Sikora’s buffer, and buffered pH is used throughout this paper (Tecnal TEC-11 MP). Soil organic C and total soil N to 20 cm were determined by dry combustion on a Leco TruMac CN Analyzer (Leco Corporation, St. Joseph, Michigan, USA). We measured soil texture (% clay, sand, and silt) using a total dispersion method with sodium hexametaphosphate (Empresa Brasileira de Pesquisa Agropecuaria (EMBRAPA), 1997). Bulk density was estimated from the mass of 10 fresh soil cores per treatment, with subsequent accounting for soil moisture.


 


We measured C mineralization as a baseline indicator of soil microbial activity and biological soil fertility at the start of the experiment, and N mineralization as a response variable following the second year of cover crop treatments. Specifically, using the baseline soil sample, we conducted a short-term (24-hour) C mineralization assay to determine potentially mineralizable C (PMC), which measures the flux of CO2 following re-wetting of previously air-dried, sieved soil using a Li-Cor (Franzluebbers et al., 2000; Hurisso et al., 2016). To measure potentially mineralizable N (PMN), we conducted a two-week aerobic incubation using fresh soil collected at vegetable crop planting in the second year of the experiment (spring 2019), two weeks after cover crop and weed biomass incorporation (Drinkwater et al., 1996; Appendix B.2). PMN was calculated as the difference between extractable soil inorganic N (NH4+ and NO3-) at the start and end of the incubation. We used pre-incubation extractable inorganic N concentration (mg/kg) as a measure of soil inorganic N availability at vegetable crop planting.


 Cover crop sampling and analysis


Cover crop biomass sampling took place from August 28-September 2 in 2018 and September 4-10 in 2019. During peak flowering of both common vetch and black oat, we destructively harvested the aboveground biomass of cover crop mixtures and weedy fallows from two 0.5 x 0.5 m quadrats of each treatment per field. We took care to avoid treatment edges, cut plant material to the soil surface, and separated harvested plant material by species, grouping all weeds together. Aboveground biomass was dried in a forced-air oven at 60 °C for 48 hours. Following grinding in a Wiley mill to 2 mm, % N and C content was determined by dry combustion on an elemental analyzer (Leco, as above). Community-weighted means were calculated for total aboveground biomass C and N in cover crop species and weeds, to determine the overall C and N inputs to soil following incorporation of biomass on each farm. We measured biological N2 fixation in inoculated common vetch from the cover crop phase of the experiment in 2018 and 2019. Vetch N fixation was estimated using the 15N natural abundance method (Shearer and Kohl, 1986), which compares stable N isotope ratios in the legume and reference species (oat monocultures) (Appendix C).


 Vegetable crop sampling and analysis


To capture the full production period of both cucumber and pea crops, yield was measured in two harvests, which were approximately 14 days apart on each farm. Harvest dates ran from November 16-December 6 in 2018 and November 20-December 4 in 2019. We measured yield by weighing all harvestable fruit from three designated, representative row sections (6 plants on average per row) per crop type per treatment. Rows were sampled from the center of each treatment to reduce edge effects. We calculated yield as total crop production (g) per plant harvested in each row. Mean yield for each crop type was calculated as the average of the three harvested rows per treatment on a per-plant basis and was then aggregated to the plot and hectare level based on experimental planting densities. Total N harvested, or “N yield”, was calculated for all treatments by multiplying the % N in each vegetable crop by its yield (kg/ha) after accounting for crop water content. Using plot-level yield data, we subsequently calculated the relative yield total (Land Equivalent Ratio, LER) for intercrop treatments by farm using the standard equation (Vandermeer, 1989) (Table 1). As a relative measure of total crop production per area, when mean LER > 1, intercrops were considered to have “overyielded” compared to their component monocrops. We calculated the LER for N yield (LERN in kg N/ha) using the same formula.


 


At the second vegetable harvest, we destructively sampled whole aboveground crop biomass, including residues and remaining fruits, from the designated experimental rows.  Following the harvest, a minimum of six representative cucumbers per treatment (from different plants) per farm were washed in deionized water, air-dried, sliced, and the middle sections were combined into a homogenized, composite sample of ~100 g and then dried for one week at 60 ºC. All peas from each treatment’s subplot were washed in deionized water, air-dried, de-stemmed, chopped, and each homogenized sample (35-60 g fresh material) was subsequently dried at 60 ºC in a forced-air oven for 48 h to one week, until fully desiccated. Dried vegetable biomass residues were ground using a Wiley mill; vegetable crop samples were ground in a coffee grinder; and all vegetable samples were analyzed for % C and N on a LECO elemental analyzer.


See supplemental material from Stratton et al. 2022 for further detailed information on methods.