THE EFFECT OF COVER CROP ON SOIL CARBON AND SOIL WATER RETENTION IN TOPOGRAPHICALLY DIVERSE TERRAIN By Jordan Beehler A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree of Crop and Soil Sciences Master of Science 2016 ABSTRACT THE EFFECT OF COVER CROP ON SOIL CARBON AND SOIL WATER RETENTION IN TOPOGRAPHICALLY DIVERSE TERRAIN By Jordan Beehler Farmers must consider real-world landscape variability to maximize yields and minimize environmental impacts when using cover crops in corn-based cropping systems. In the Midwest corn belt much of the variability farmers encounter is due to the topographical diversity of the undulating landscape. The objectives of this study are to explore the interactive effects of the presence or absence of cover crops and topography (summit, slope, and depression positions) on total soil organic carbon, on its labile form, particulate organic carbon, carbon dioxide emissions, cover crop decomposition, soil water retention, and crop growth. A cereal rye cover crop was established each fall after the main crop (corn and soybean) harvest from 2011-2015 at two experimental sites, Kellogg Biological Station (Kellogg) and Mason, which have loam and fine sandy loam soils, respectively. Main crop and cover crop growth were both higher in the depression, but did not differ in topographical position and main crop growth did not differ in the presence/absence of the cover crop. In the absence of the cover crop, topography affected particulate organic carbon but not in the presence of the cover crop (p<0.1). Decomposition and carbon dioxide emission followed the trend depression>summit>slope, but no one variable accounted for the distribution of particulate organic carbon. Total organic carbon and water retention were not affected by the cover crop in this study. A long-term study may reveal additional significant changes in the presence of a cover crop not detectable in a five-year study. iii ACKNOWLEDGMENTS I would like to thank, first and foremost, my family and friends for all their help and support. Also a big thank you to the undergraduate students, graduate students, and my advisor for their support and guidance in the field and in the office and for sharing many laughs with me along the way. iv TABLE OF CONTENTS LIST OF TABLES.. LIST OF FIGURES. KEY TO ABBREVIATIONS CHAPTER ONE Overview Importance of cover crops in corn-based systems Labile carbon pools as the first stage of cover crop driven improvement in soil health Effect of topography on decomposition of cover crop residue Effect of topography on soil respiration Importance of soil water retention Implications Hypotheses.12 CHAPTER TWO Site description Experimental design Plant yield biomass and carbon content sampling.. Particulate organic carbon Rye decomposition CO2 emissions Total soil organic carbon Water retention. Statistical analysis. CHAPTER THREE - RESULTS0 Overview Aboveground biomass..20 Particulate organic carbon.22 Rye decomposition CO2 emissions Total soil organic carbon24 Water retention CHAPTER FOUR DISCUSSIO Conclusions APPENDIX.. LITERATURE CITED.2 v LIST OF TABLES Table 1. Overview of field operations and timing for research plots at Kellogg and Mason Table 2. Soil organic carbon concentration (%) at the four studied depths of the two experimental sites Table 3. Volumetric water content at all pressures at the two studied depths of Kellogg and Mason in the Table 4. Soil texture at Mason and Kellogg for depths 0-10 cm, 10-20 cm, 20-40 cm, and 40-60 cm at three topographical positions (depression, slope, summit) measured at the beginning of the 37 Table 5. Summary of nitrogen fertilizer type, vi LIST OF FIGURES Figure 1. Flowchart depicting hypothesized and observed effects of cover crop presence on carbon cycling in this study Figure 2a. 4 Figure 2b. Digital elevation map of Mason .. Figure 3. Figure 4. Average total carbon (kg ha-1) of rye biomass across topographical position Figure 5. Figure 6. Particulate organic carbon across both experimental sites at 0-10 cm depth for Figure 7. ...23 Figure 8. CO2 emissions from 2012, 2013, and 2014 across topographical positions in Figure 9. Volumetric water content (water retention) (cm3 cm-3) at field capacity (-0.33 bar) for all ye Figure 10. Flowchart depicting hypothesized and observed effects with the factors in Figure 11. Monthly average air temperature and precipitation at Kellogg and Mason for each year of the 39 vii KEY TO ABBREVIATIONS POC Particulate organic carbon CO2 Carbon dioxide C Carbon 1 CHAPTER ONE - INTRODUCTION Overview Farmers must consider real-world problems and variability to maximize yields and minimize environmental impacts when using cover crops in corn-based cropping systems. In the Midwest corn belt, much of the variability that farmers encounter is due to the topographical diversity of the undulating landscape. This topographical diversity is a controlling factor for many soil properties, including the distribution of soil organic carbon. Enhancing a rotation by the inclusion of a cover crop has the potential to accumulate soil organic carbon (West and Post, 2002; Follet, 2001). Cereal rye (Secale cereale L.) is especially suited for this purpose because it is a high biomass grass and is shown to increase soil organic carbon (Kuo and Jellum, 2000; Kaspar et al., 2005; Reicosky and Forcella, 1998). Even though rye is beneficial for soil organic carbon accumulation, it can have negative effects, for example, decreasing the soil moisture or nitrogen availability for the subsequent main crop, or inhibiting emergence or growth of the main crop through allelopathic effects (Duiker and Curran, 2004; Wagner-Riddle et al., 1994; Munawar et al., 1990). The advantages and disadvantages of the use of rye as a cover crop in regards to soil organic carbon improvement are not homogeneously spread across landscapes with topographical diversity. Cover crops, such as rye, can grow at variable rates across topographies, commonly with better biomass accumulation in the flat, depression positions (Munoz et al., 2014). Thus, one can expect that greater soil organic carbon benefits from a rye cover crop can occur in topographical depressions. Yet, recent observations from topographically diverse agricultural fields in Michigan demonstrated that the magnitude of the cover crop effect on soil organic carbon was higher in the summit and slope positions, as compared to depressions (Ladoni et al., 2016). This indicates that other factors, besides the overall amount of the aboveground cover crop biomass production, influence the magnitude of the soil organic carbon benefits in topographically diverse agricultural landscapes. 2 Spatial distribution patterns of soil moisture and temperature are among the factors that likely contribute to soil organic carbon accrual. Topographical differences control soil moisture and soil temperature, both of which regulate the environment for microbial growth and activity as well as provide the environment for crop growth. These factors greatly affect the decomposition of freshly added plant residue. In addition to the topographical controls on physical properties effecting soil microbes, the presence of cover crops in rotation promotes microbial activity. This effect can further enhance the role of topographical differences in the environment for crop growth and decomposition (Wickings et al. 2016). Decomposition rates of a rye cover crop at different topographical positions have not been extensively studied, but it is possible to hypothesize the anticipated effects. Prime decomposition conditions are difficult to obtain in topographically diverse terrain. The summits and slopes tend to be too hot and dry, while depressions are too cool and wet. Despite this, decomposition often happens at a faster rate in depression positions. This may be caused by an increase in microbial activity by up to 55% in depression positions relative to summits and slopes (Wickings et al., 2016). This could be due to many reasons, one of which is the spatial distribution of nitrogen in topographically diverse terrain. Nitrogen, in the form of nitrate, is a mobile nutrient which leads to an accumulation of available nitrogen in the depression positions. Higher biomass accumulation and faster decomposition rates of labile organic matter in depression areas leads to a larger magnitude and faster turnover of carbon relative to summits and slopes. This has a few implications regarding the sequestration, or protection, of carbon in the soil. Carbon can be sequestered in the soil by three main methods; physical stabilization, attachment to silt and clay particles, and biochemical transformation (Six et al., 2002). A faster turnover of carbon in depression positions may lead to greater physical stabilization of soil organic matter by the promotion of 3 increased soil structure through the formation of microaggregates, which protect the soil organic matter (Six et a., 2002). The other two methods of soil sequestration, attachment to silt and clay particles and biochemical transformation, rely on a slower carbon turnover to most effectively protect carbon in the soil. Physiochemical protection, or the sorption of carbon to existing soil surfaces, relies on many factors which are promoted through a slower turnover of carbon including decomposition resulting in smaller C particles, reorientation of particles, and chemical bonding (Jastrow et al., 2007). Biochemical transformation is when soil organic carbon is transformed from an organic to a mineral form of carbon. This is a biologically driven process that is initiated by decomposers who release molecules which can then react to form stable carbon (Jastrow et al., 2007). Since these methods are promoted by slower C turnover, they could be more prominent in the summit and slope positions and not in the depression position. There are many complications in studying impacts of topographical gradients on soils. For example, topographical controls on soil moisture and soil temperature add variability that impact many stages of the carbon cycle in relation to sequestering soil carbon, most of which have not been extensively studied. Root exudates, underlying geological differences, and the nutrient flow path along the topographical gradient may affect cover crop growth and decomposition. In addition, the potential for movement of soil and plant material through erosion and deposition may have a large impact not only on plant growth and decomposition through the availability of nutrients, but the distribution (or redistribution) of soil carbon along the topographical gradient. In addition, topography plays a role in the distribution of nutrients, especially nitrogen, which can impact the variability of biomass accumulation. Nitrogen, most commonly in the form of nitrates, collect at the footslope (depression) leading to higher above ground biomass in that position than in the shoulder (slope) position (Mbonimpa et al., 2016). Biomass growth variability is not only spatially along 4 the topographical gradient but vertically within the plant. Plants are able to direct the allocation of energy to the growth of above or below ground biomass based on soil resources. For example, when there are limited soil resources, plants promote belowground biomass growth but when nutrient availability is high, growth is favored in the above ground biomass (Tateno et al., 2004). Changes in atmospheric carbon, along with other factors, are leading to a modification in climate which is resulting in greater variability of precipitation patterns (Trenberth, 2011). This variability means that precipitation patterns are less predictable, but we can expect that historically wet areas will become more wet, while historically arid areas will become more dry (Dore, 2005). As a result of the altered precipitation patterns, it is expected that soil moisture patterns along topographical gradients, which are normally drier in the upslope position and wetter in the downslope position, will be exacerbated (Western et al., 1999). The objectives of this study are to explore the effect of rye cover crop on soil carbon accrual and changes in soil water retention characteristics in a topographically diverse agricultural landscape. Figure 1 outlines the hypothesized influences of cover crops on soil carbon processes in a corn-based cropping system that were considered in this study, the solid arrows are an observed effect and the dotted arrows are a hypothesized effect. It also indicates, by the outlined boxes, the measurements that were conducted in order to explore the interacting contribution of topography and cover crop use on soil carbon processes. The effects of a rye cover crop have the potential to impact the amount of long-term soil organic carbon but only by going through other processes. 5 Figure 1. Flowchart depicting hypothesized and observed effects of cover crop presence on carbon cycling in this study. Importance of cover crops in corn-based systems The use of cover crops, like rye, is a beneficial management practice in corn-based cropping systems in the Midwest. Cover crops, which are grown when the soil would otherwise be fallow, have many positive effects on the biological, chemical, and physical properties of the cropping system in general, such as suppression of weeds and pests, improvements of soil and water quality, and stimulation of nutrient cycles (Snapp et al., 2005). A cereal rye (Secale cereale L.) cover crop, used in this study, is appropriate when looking to increase biomass production and decrease soil erosion in the fall and winter in Midwest corn and soybean rotations (Raimbult et al., 1989; Bruin et al., 2005). Biomass production and protection from soil erosion are two important characteristics of rye; not only does it facilitate the addition of organic matter through biomass, it also decreases (but does not eliminate) the movement of organic matter through erosion. Both characteristics can facilitate keeping organic matter in the cropping system, which is commonly used as an indicator of soil health (Reeves, 1997). 6 Continuous crop growth and tillage without returning plant residue to the soil can deplete soil organic matter which is why, when returned and incorporated into the soil of the cropping system, cover crops can increase soil health. The absence of soil organic matter leads to reduced soil health and productivity by way of reduced microbial action, increased compaction, and less stable aggregates (Allison, 1973; Van Doren et al., 1976; Campbell and Souster, 1982; Cambardella and Elliott, 1993). Cover crops, especially high biomass, non-leguminous grasses like cereal rye, can provide an additional input of organic carbon into corn-based cropping rotations which can reverse those effects (Sainju et al., 2002; Dabney et al., 2001; Papadopoulos et al., 2006). Labile carbon pools as the first stage of cover crop driven improvement in soil health I hypothesize that the additional input of aboveground biomass from the rye cover crop leads to greater soil organic carbon levels (Figure 1). However, initially I expect to see increases in soil labile carbon pools, e.g., particulate organic carbon. This hypothesis is based on the relatively slow response of total soil organic carbon to management practices (Plaza-Bonilla et al., 2014; Haynes, 2000). Labile carbon, like particulate organic carbon, can show changes on the scale of a few days to a few years (Nascente et al., 2013). Particulate organic matter is an indicator of the amount of stable carbon that will be protected and stored in the soil (Post and Kwon, 2000; Six et al., 2000). Therefore, by using particulate organic carbon as a metric in this five-year study, the changes in soil carbon due to the cover crop can be more easily detected For this study, particulate organic carbon is defined as the active fraction of the total organic carbon pool that consists of organic carbon particles between 0.053 and 2 mm in size which react quickly to changes in management practices (Cambardella and Elliot, 1992; Sequeira et al., 2011; Ladoni et al., 2015). Even though it is highly variable, particulate organic carbon is reflective of management changes (Ladoni et al., 2015). This fraction of the organic carbon pool is also especially prone to 7 redistribution related to downhill erosion and eventually collects in the depression areas, but does not always lead to higher concentrations in the intermediate slope positions when compared to the summits (Dungait et al, 2013). This pattern is likely because of particle size distribution, since high organic matter content is usually correlated with fine material which because of its size and mobility is more likely found in depressions, rather than summits or slopes (Dungait et al., 2013; Parton et al., 1987). Particulate organic carbon can constitute 42-74% of the total organic carbon pool and is preferentially lost during cultivation (Chan, 2001). This loss can be observed through the decomposition rates and carbon dioxide emissions in and from the soil. Effect of topography on decomposition of cover crop residue I hypothesize that greater decomposition of rye residue after rye incorporation in spring will take place in topographical depressions (Figure 1). Greater decomposition will happen in topographical depressions because of more ideal conditions for decomposition including soil moisture and soil temperature leading to an increase in microbial activity. However, the magnitude of topographical effect will vary from year-to-year depending on temporal patterns in spring temperature and precipitation. Decomposition rates and the factors that control them are the major influences on carbon inputs into the system (Henriksen and Breland, 1999). The type and composition of the decomposing litter exerts less control on the rate of decomposition of the labile and long-term pools of organic carbon than other controlling factors, such as soil moisture and temperature (Adair, 2008). The main controlling factors of the rate of decomposition are soil moisture and climate variables, such as rainfall and temperature (Eijsackers and Zehnder, 1990; Andren, et al. 1993; Adair, 2008). In addition, soil texture influences decomposition because in general, coarser textured soils have a slower decomposition rate (Rovira, 2002). However, regardless of controlling factors, rye decomposition was reported to be the highest the week after incorporation because the input of the rye organic matter produced a short-term 8 boost of microbial growth and activity (Lundquist, 1999). This short-term burst of microbial action, lasting about a week after incorporation, happens during the time between the termination of the rye cover crop and the planting of the subsequent main crop. Effect of topography on soil respiration Greater organic inputs and their faster decomposition typically result in greater soil respiration, that is greater CO2 emissions from soil. Thus, increases in CO2 emissions due to cover crop use can be regarded as precursors of subsequent increases in soil organic carbon. I hypothesize that greater increases in soil CO2 emissions due to cover crop use will be observed in topographical depressions because of their greater biomass inputs and faster decomposition (Figure 1). Carbon dioxide emissions from the soil are mainly controlled by soil moisture and climate variables, such as temperature and rainfall which drive microbial action, as well as soil organic matter content (Raich and Potter, 1995; Trumbore et al., 1996; Srivastava et al., 2012; Harrison-Kirk et al., 2013). A preceding study at the two experimental sites used in my work examined effects of topography, cover crop presence, and tillage on soil CO2 emissions during growing seasons of 2012-2013. Across all topographical positions the presence of cover crops resulted in greater CO2 emissions. The magnitude of that effect differed between summits, slopes, and depressions, with the greatest increase occurring in depressions and the smallest in summits, and was inversely proportional to the biomass input from the cover crop (Negassa et al., 2015). The effects of topography and cover crops on CO2 emissions have not been extensively published, so this study aims to add to this area of impact on the soil carbon cycle in corn-based cropping rotations. 9 Importance of soil water retention Topography also plays a role in the spatial variability of soil physical properties, which are themselves highly correlated to soil organic carbon content (Moore et al., 1992; Tromp-van Meerveld and McDonnell, 2006; Romano and Palladino, 2002). One such property of significant value to plant growth is soil water retention which is strongly related to soil organic carbon. On a per volume basis an organic matter increase significantly affects the physical properties of the soil by changing bulk density and in turn, water retention (Hudson, 1994; Adams, 1973; Rawls, 1983; Gupta and Larson, 1979). An increase in organic matter in sandy soils, like the soils in this study, can lead to an increase of water retention (Rawls et al., 2003). Organic matter content in the soil also affects soil aggregation, a key component of soil structure and water retention capability (Franzluebbers, 2002). The effect of organic matter content on water retention can be observed most clearly at field capacity, which can have an effect on cover crop and main crop growth (Jong et al., 1983). By understanding the effects of topography and cover cropping on water retention we can better understand the long-term implications on soil carbon and soil health. Implications It is important to better understand where and how organic carbon cycles in the soil, which holds roughly three times as much carbon as the atmosphevegetation, since scientists continue to observe increases in atmospheric carbon dioxide leading to climate change that will impact agricultural systems, especially in the Midwest corn belt (Lal, 2004). Some agricultural lands have lost one-half to two-thirds of their original total soil organic carbon pool, so by gaining a better understanding of the dynamic turnover and flux of carbon in agricultural systems enhanced with cover crops we will be able to sequester carbon back in to agricultural lands (Lal, 2004). By increasing the soil organic carbon pool in a corn cropping system by 1 ton/ha, corn yield may 10 increase by up to 20 kg/ha, making carbon sequestration not only environmentally advantageous but also potentially financially profitable (Lal, 2004). Increasing soil organic carbon can also increase infiltration and soil fertility, decrease erosion, minimize compaction, and impede pesticide movement, all of which are signs of a healthy, productive soil (Reicosky, 2003). Although this may seem like a very straightforward chain of events grow a cover crop, add more organic carbon to the soil, and then reap the benefits - spatial variability from topographical diversity that impacts cover crop growth and the subsequent impacts on the other variables observed in this study is not extensively researched (Munoz et al., 2014). Real-world, field-scale topographical diversity is often under-represented in agricultural research and experiments. To minimize variability in scientific research, plots chosen are often on flat land with no topographical variance (Munoz, 2014). Despite this, cover crop growth and biomass production can be influenced by the distribution of soil water, which is highly controlled by topography (Munoz et al., 2014; Corre et al., 2002). Topography also affects the distribution of nutrients, like nitrogen, and soil temperature which impacts the growth and microbial environment along the topographical gradient (Moore et al., 1993; Kang et al., 2003; Zhu and Lin, 2011; Bennie et al., 2008; Yimer et al., 2006). The differing environments between topographical positions that drive the growth of cover crop biomass provide highly variable inputs of organic carbon. In addition, topography significantly affects the redistribution (erosion and deposition) of soil across the landscape which is correlated with the distribution of soil organic carbon after the initial carbon input from the cover crop (Ritchie et al., 2007). Both modeling and field studies have found a strong correlation with slope (i.e. topographical position) and soil organic carbon distribution (e.g. Moore et al., 1992; Terra et al., 2004; Kunkle et al., 2011). An increasing popularity in precision agriculture makes studying topographical differences in corn-soy cropping rotations even more pertinent. By adding topographical variables into precision ag management, farmers will be able to be more effective in their use of the land and more efficient in 11 their productivity. Adding topography as a factor in observing the impacts of a rye cover crop in a corn-based cropping rotation makes this study novel and necessary because topography indirectly or directly impacts every variable involved in production agriculture as well as soil carbon sequestration. 12 Hypotheses The principal hypothesis of this study is that topography will be the main cause of observable differences in soil carbon, particulate organic carbon, rye decomposition, CO2 emissions, and water retention because topography plays an important role in spatial variation of biomass growth, distribution of soil organic matter, and soil physical properties (Dharmakeerthi et al., 2005). It is also hypothesized that cover crop effects will be observed through a greater input of aboveground biomass, greater magnitude of carbon flux through CO2 emissions and rye decomposition, and a larger labile carbon pool (Figure 1). Adding cover crops as an additional source of carbon into topographically diverse terrain adds a new and not extensively studied layer of complexity into the understanding of carbon dynamics in corn-based cropping systems in the Midwest. 13 CHAPTER TWO MATERIALS AND METHODS Site description The two experimental sites are located at the Kellogg Biological Station (Kellogg) and the Mason Research farm in southwest Michigan and mid-Michigan, respectively. Both sites average around 76 cm of rainfall in addition to 76 cm of snowfall annually. Temperatures at both sites average -3C in the winter and 20C in the summer. The Mason site, located in Ingham County, is situated near the Michigan State University (MSU) campus. The Kellogg site, located in Barry County, is on the grounds of the Kellogg Biological Station, a large off-site research campus operated by MSU. Soils at the Mason site are classified as Marlette fine sandy loams (Oxyaquic Glossudalfs) and soils at the Kellogg site are classified as Kalamazoo loams (Typic Hapludalfs). A complete description of soil texture at all depths and topographical positions can be found in the Appendix. Experimental design Both experimental sites have contrasting topographies, which enabled incorporating topographical position as a studied factor. At Kellogg, blocks 1 and 2 are slope, blocks 3 and 4 are depression, and blocks 5 and 6 are summit (Fig. 2a). At Mason, blocks 1 and 2 are summit, block 3 and 4 are slope, and block 5 and 6 are depression (Fig. 2b). At each of the sites the experimental design was a split-split plot with whole plot factor, topographical position, in a randomized complete block design with two replications, sub-plot factor, main crop (corn and soybean), and sub-sub-plot factor cover crop (presence and absence) (Fig. 3). The study was conducted from 2011 to 2015. 14 Figure 2a. Digital elevation map of Kellogg site. Figure 2b. Digital elevation map of Mason site Whole plots (blocks) are outlined and Whole plots (blocks) are outlined and numbered. numbered. Figure 3. Representation of treatment splits within each block. Each block is replicated twice within the whole plot topography. Each block (thick outer box) has both corn (diagonal lines) and soybean (vertical lines). Each crop is then split into two halves: with and without winter rye cover crop (thick inner box). Field management and data collection Field sites were managed by our research group with help from farm assistants associated with each research site. Conventional practices were used in plowing, planting, fertilizing, spraying herbicides, and harvesting consistent with practices in the region. An overview of field operations and the timing of those operations can be found in Table 1. 15 Task Projected Date to Complete Task Rye Cover Crop and Weed Sampling First week of May Rye Cover Crop Termination First week of May (approximately 2 weeks before corn planting) Spring Fertilizer and Lime Application Third week of May (before tillage) Tillage Third week of May Plant Corn Last week of May Plant Soybean Last week of May Post Emergence Herbicide Spray Mid-June Corn Harvest First or second week of October Soybean Harvest First or second week of October Plant Rye Cover Crop Second week of October (after harvest) Table 1. Overview of field operations and timing for research plots at Kellogg and Mason. The cereal rye cover crop was established each fall in early October, beginning in 2011.From 2011-2013, rye was sown at 112 kg/ha with a John Deere 15 foot no-till drill. In 2014-2015 rye population was increased to establish a better stand and was sown at 145 kg/ha. Nitrogen fertilizer was applied based on yearly soil tests and recommendations from the Michigan State Soil and Plant Nutrient Laboratory. A summary of nitrogen fertilizer application can be found in the Appendix. Rye was terminated by herbicide then was chisel plowed, followed by a soil finisher to establish a seed bed for the subsequent main crop. Weed growth was terminated by the use of herbicides at least once per growing season with additional weed control by herbicide if necessary. Three composite soil cores were taken in the spring of 2011, 2013, and 2015 from depths 0-10 cm, 10-20 cm, 20-40 cm, and 40-60 cm. Total soil organic carbon was measured at all depths for those three sampling events. Particulate organic carbon samples were collected separately once a year in the 16 spring from 2012 to 2014 at the 0-10 cm depth. Water retention samples were also collected separately in the spring of 2011, 2013, and 2015 at depths of 0-10 cm and 10-20 cm. Carbon dioxide emissions were taken during the main crop growing season in 2012, 2013, and 2014. The litter bag decomposition study was done during the growing season of 2015 and 2016. Plant biomass, both cover crop and main crop, were collected in relation to growth and harvest or termination. Data collected from these sites was based on standardized protocols outlined by Kladivko et al., 2014. Detailed descriptions of the specific data collection/measurement procedures are provided below. Plant yield biomass and carbon content sampling Biomass of corn and soybean was collected at harvest every year. Whole plants of corn and soybean in previously assigned yield rows were collected, then the grain was separated from the vegetative biomass and both were weighed to obtain measurements for yield and vegetative biomass. The remaining plant material in the field was harvested mechanically. For corn, average yield was calculated using the grain collected from the yield rows. For soybean, yield was recorded by the combine. Rye biomass was collected immediately before the termination of the cover crop every year. Three random samples per plot were sampling square was cut at the soil surface to obtain an aboveground biomass. The rye was weighed, dried at 40 C, and then weighed again to obtain a dry biomass measurement. After being dried, corn, soybean, and rye were finely ground in preparation for combustion analysis. The carbon content of the ground plant material was found using a Costech ECS 4010 CHNSO Analyzer (Costech Analytical Technologies, Inc, Valencia, CA). 17 Particulate organic carbon Approximately 20 g of soil was subsampled from the sieved and air dried soil cores to use for particulate organic carbon analysis. Particulate organic carbon was chemically dispersed using a 5% sodium hexametaphosphate solution. The soil solution was then passed through a 53-micron sieve. The contents of the sieve were oven dried at 60 C and ground in the 8500 Shatterbox (Spex Sample Prep, Metuchen, NJ) (Cambardella and Elliott, 1992). Final carbon analysis of the sample was done using a Costech ECS 4010 CHNSO Analyzer (Costech Analytical Tech Inc., Valencia, CA). Rye decomposition Rye was collected from each field site prior to termination by pulling the plant, gathering both aboveground and belowground biomass, to use for decomposition by litter bag methodology. Litter bags were 10 cm by 10 cm and constructed with 0.028 cm plastic mesh. The rye was oven dried at 40C and a random subsample of 5 to 7 grams of hand-cut dried rye, measuring approximately 5-7 cm in length, was placed in each litter bag. In 2015, approximately 60 litter bags were buried about 10 cm deep in plots with cover crop at the Mason site. In 2016, approximately 340 litter bags were buried about 10 cm deep in plots with and without cover crop at the Mason and Kellogg site. Litter bags were removed at three time points, 1, 3, and 5 weeks, at which time the decomposed rye was then carefully cleaned of soil and other debris by hand dusting, at which time the decomposed rye was oven dried at 40C and weighed (Alef, 1995). Commonly in litterbag decomposition studies, a portion of the contents are ashed after decomposition (Nadelhoffer et al., 1999) but no portion of decomposed material was ashed in this study. 18 CO2 emissions Carbon dioxide gas samples were taken biweekly from 8:00 AM to 12:00 PM during the growing season of 2012, 2013, and 2014 from cover and no cover plots in corn and soybean across all topographical positions at both sites. Aluminum static flux chambers (0.375 m x 0.75 m) were manufactured locally (Parkin and Venterea, 2010). The chambers consisted of the anchor, which was installed in the field after planting of the main crop, and a removable cap with a vent tube and sampling port. Once the chamber and cap were in place, CO2 concentration was sampled from the surface soil every two minutes for at least 14 minutes by using an infrared Photoacoustic Spectroscopy (PAS) (INNOVA Air Tech Instruments, Ballerup, Denmark). Soil temperature using a pocket thermometer (Taylor Precision Products, Oak Brook, Illinois) and soil moisture using time domain reflectance (IMKO HD-2 IMKO GmbH, Ettlingen, Germany) were measured at the time of CO2 sampling at three points around the chamber. Calculations from the CO2 concentration, soil temperature, and soil moisture were completed to obtain the rate of CO2 C (CO2) emissions (Iqbal et al., 2013). Total soil organic carbon Soil samples were collected using a Giddings hydraulic probe (Giddings Machine Company, Windsor, CO) (7.6 cm in diameter). Approximately 0.7 kg of soil from each depth increment was wet sieved to pass a 2 mm sieve and air dried. A small subsample, approximately 5 g, was then ground to a fine powder using an 8500 Shatterbox (Spex Sample Prep, Metuchen, NJ) in preparation for flash combustion analysis by the Carlo Erba EA 1108 (CE Elantec Inc., Lakewood, NJ). Combustion analysis of soil samples was completed by the USDA-ARS in St. Paul, Minnesota under the direction of Dr. Kurt Spokas. Samples were sent to that lab so that values could be corrected for small amounts of inorganic carbon at the Mason site through the acetic acid neutralization method, which was an unnecessary correction for the Kellogg site. 19 Water retention Soil samples for water retention analysis were collected one time per year in the spring before corn or soybean planting in brass rings 5.5 cm in diameter and 3 cm in height. There was no continuous measurement of soil moisture dynamics. Pressure plates were used to extract water at pressures of 0.05 bar, 0.1 bar, 0.33 bar, 1 bar, and 3 bar. The water content at 15 bar was measured by drying subsamples of the soil used from the pressure plate method and placing them in a desiccator above an oversaturated KCl solution (500 g KCl/1 L water) for two months. For both methods, the difference in mass was used to calculate water content. Statistical analysis Statistical analysis was performed using MIXED procedure of SAS 9.4 (SAS Institute Inc., Cary, NC). The statistical model for the analysis included presence/absence of cover crop, topographical position, and the interaction between them as fixed factors. The previous crop, corn or soybean, was not used because it introduced too much error into the model to be a useful interaction in this study. Block was nested in topography at each site was included in the model as the random factor and was used as an error term for testing the main effect of topography. When the interaction between topography and cover crops was statistically significant (p<0.1), we used slicing to assess the effect of cover crop presence at each topography level and the effect of topography at each level of the cover crop factor. When slicing effects were statistically significant (p<0.1), comparisons between the means were conducted using t-tests. Because of high variability of the collected field data we reported the results that were statistically significant at both 0.05 and 0.1 levels of significance. For the analysis of total soil carbon, 2011 total organic soil carbon data from both Mason and Kellogg was added to the model and used as a covariate when assessing 2015 values from both sites together. 20 CHAPTER THREE RESULTS Overview The results are presented in the order consistent with the diagram of the hypothesized influences of cover crops on soil carbon processes that were considered in this study (Fig. 1). It starts with topographical and cover crop effects on main crop biomass, followed by topography effect on rye biomass and C inputs, then topographical and cover crop effects on particulate organic carbon (POC), litter decomposition, and CO2 emissions from soil, and finally on total soil organic C and soil water retention. Aboveground biomass Total carbon from aboveground biomass was split into two categories, rye cover crop and main crop (both corn and soybean). The average total carbon in the rye cover crop showed different trends at the two sites. At Kellogg, the slope exhibited a statistically significant higher average total carbon in rye cover crop biomass than the summit and depression positions. At Mason, the depression exhibited a higher average total carbon of rye biomass, but there were no statistically significant differences between the topographies (Fig. 4). 21 Figure 4. Average total carbon (kg ha-1) of rye biomass across topographical position for each site. Letters represent statistically significant differences among the topographical positions at each site (p<0.1). Average total carbon from the main crop aboveground biomass did not differ between sites. It also did not differ in the presences or absence of the rye cover crop. Main crop total carbon was statistically higher in the depression position than in the summits and slopes (Fig. 5). Figure 5. Average main crop (corn and soybean) total carbon for all years. Letters mark statistically significant differences between topographical positions (p<0.1). 22 Particulate organic carbon Overall, the depression positions tended to have more particulate organic carbon than summits and slopes, however the magnitude of the differences varied in cover and no cover treatments (Fig. 6) In the absence of the rye cover crop, POC was significantly lower in summit and slope positions compared same numeric pattern as those in in the no cover treatment, but the differences were not statistically significant. Overall, presence of cover crop tended to result in greater POC, however, the size of the effect varied depending on topography. The presence of the cover crop resulted in significantly higher particulate organic carbon in the slope and summit tdepression positions was not statistically significant. Figure 6. Particulate organic carbon across both experimental sites at 0-10 cm depth for no cover and cover plots. Letters within each cover treatment represent statistically significant differences among the topographical positions (p<0.1). Asterisks mark the topographical positions where the differences between the presence and absence of the cover crop was statistically significant (p<0.1). 23 Rye decomposition At each topographical position (summit, slope, and depression) in both 2015 and 2016, each time point was statistically different (p<0.1), indicating the occurrence of decomposition. In 2016, there were no differences in mass loss between the presence and absence of a cover crop. In 2015, the study was only completed in the presence of a cover crop, so no comparison can be made. In 2015, topographical differences only existed at week 5, where the depression position had significantly higher mass loss than the summit and slope positions (Fig. 7). In 2016, decomposition followed the trend depression>summit>slope, which was statistically significant at week 1 and week 3 (Fig. 7). Figure 7. Rye decomposition (%) at Kellogg and Mason for 2015 and 2016. Letters indicate topographical differences in each year for each sampling date (p<0.1). CO2 emissions Carbon dioxide emissions differed significantly between the topographical positions (p<0.1). For both sites and all years studied, the amount of CO2 released followed the trend slope summit > slope but there were no significant differences between the presence and absence of the cover crop. Kaspar et al. (2006) also observed no significant differences in total soil carbon between the presence and absence of a rye cover crop. Simulated experiments have also found that the presence of a rye cover crop does not increase total soil organic carbon compared to the absence of the cover crop (Basche et al., 2016). Since changes in total soil organic carbon may not be seen until 7 to 10 years or more after changing management practices, like the addition of a cover crop into rotation, the results from the literature and my study are consistent and expected (Duiker and Lal, 1999; Al-Kaisi et al., 2005). This study is shorter than the time it takes to detect significant changes in the long-term pool of carbon. The topographical differences observed in the study were expected as they are controlled by factors with a stronger effect on carbon distribution than the presence of a cover crop. These factors include soil physical properties, such as texture, temperature, and soil water distribution, as well as the erosion and deposition of soil particles across the topographical gradient (Rovira, 2002; Eijsackers and Zehnder, 1990; Andren et al., 1993). For long-term carbon storage, it is important to consider the topographical effects and focus on locations where the effects of cover crops can be used to their highest potential. 34 Conclusions General trends in carbon content of the main crop, the cover crop, particulate organic matter, decomposition, and CO2 emissions cannot by themselves explain the soil carbon results. As shown in Figure 10, the magnitude of the effect of the cover crop on these factors are only significant in the slope and summit positions in relation to particulate organic carbon and the depression position in CO2 emissions. Overwhelmingly, there is no significant difference between the presence and absence of the cover crop. However, by comparing my study with the literature, a few conclusions can be drawn: 1) Better rye growth and a smaller gap in main crop yield between cover and no cover indicate that the largest effect in terms of carbon input from the presence of a rye cover crop is in the summit position; 2) Higher rates of decomposition and soil respiration (CO2 emission) lead to a faster turnover of labile carbon, which diminished the impact of cover crops on soil carbon, especially in the depression position; 3) Erosional deposition and redistribution of particulate organic carbon may lessen the effect of the cover crop on the spatial distribution of carbon along the topographical gradient and is an important factor to remember when studying soil carbon on topographically diverse terrain. For long-term carbon storage, it is important to consider the topographical effects and focus on locations where the effects of cover crops can be used to their highest potential. This means that by focusing on cover crop establishment on summits and slopes, cover crops could have the most impact on long-term carbon storage in Midwest corn cropping systems. 35 Figure 10. Flowchart depicting hypothesized and observed effects with the factors in this study. The magnitude of the effect of cover crop on measured variable shown by brown arrows. 36 APPENDIX 37 APPENDIX Site Topographical Position Depth (cm) Sand (%) Silt (%) Clay (%) Kellogg Depression 0 - 10 57.5 11.5 31.0 Kellogg Depression 10-20 46.2 18.7 35.0 Kellogg Depression 20 - 40 48.0 23.0 29.0 Kellogg Depression 40 - 60 50.1 22.6 27.3 Kellogg Slope 0 - 10 64.6 14.0 21.5 Kellogg Slope 10-20 68.2 16.1 15.8 Kellogg Slope 20 - 40 72.1 13.0 14.9 Kellogg Slope 40 - 60 76.7 11.7 11.7 Kellogg Summit 0 - 10 57.3 8.4 34.4 Kellogg Summit 10-20 54.1 10.6 35.4 Kellogg Summit 20 - 40 65.4 21.5 13.1 Kellogg Summit 40 - 60 82.9 14.0 3.1 Mason Depression 0 - 10 69.0 7.1 24.3 Mason Depression 10-20 65.0 7.0 28.4 Mason Depression 20 - 40 73.8 11.0 15.1 Mason Depression 40 - 60 70.6 14.3 15.3 Mason Slope 0 - 10 62.6 12.3 25.5 Mason Slope 10-20 61.0 14.4 24.9 Mason Slope 20 - 40 53.9 25.0 21.5 Mason Slope 40 - 60 54.6 24.0 21.4 Mason Summit 0 - 10 66.8 9.3 24.1 Mason Summit 10-20 61.8 13.5 25.1 Mason Summit 20 - 40 53.9 26.9 19.3 Mason Summit 40 - 60 56.4 23.5 20.3 Table 4. Soil texture at Mason and Kellogg for depths 0-10 cm, 10-20 cm, 20-40 cm, and 40-60 cm at three topographical positions (depression, slope, summit) measured at the beginning of the study in 2011. 38 Site Year Date Fertilizer Form Fertilized Crop Fertilizer Application Type Fertilizer Nitrogen Applied (kg ha-1) Kellogg 2011 7/13/2011 liquid corn injected between rows 28% UAN 135 Kellogg 2012 5/15/2012 liquid corn injected in every row 28% UAN 28 Kellogg 2012 6/29/2012 liquid corn injected between rows 28% UAN 135 Kellogg 2013 7/22/2013 solid corn broadcast urea 150 Kellogg 2014 6/12/2014 solid corn broadcast urea 46 Kellogg 2014 7/17/2014 solid corn broadcast urea 46 Kellogg 2015 5/27/2015 solid corn broadcast urea 30 Kellogg 2015 7/6/2015 solid corn broadcast urea 90 Mason 2011 7/11/2011 liquid corn broadcast between row 28-0-0 151.2 Mason 2012 5/17/2012 solid corn broadcast urea 46 Mason 2012 7/3/2012 liquid corn injected between rows 28% UAN 135 Mason 2013 7/22/2013 solid corn broadcast urea 150 Mason 2014 7/16/2014 solid corn broadcast urea 150 Mason 2015 5/28/2015 solid corn broadcast urea 30 Mason 2015 7/14/2015 solid corn broadcast urea 90 Table 5. Summary of nitrogen fertilizer type, timing, and amount for both sites. 39 Figure 11. Monthly average air temperature and precipitation at Kellogg and Mason for each year of the study. A B 40 Figure 11 (contd) C D 41 Figure 11 (contd) E 42 LITERATURE CITED 43 LITERATURE CITED Adair, E. C., Parton, W. J., Del Grosso, S. J., Silver, W. L., Harmon, M. E., Hall, S. A., ... & Hart, S. C. (2008). Simple threepool model accurately describes patterns of longterm litter decomposition in diverse climates.Global Change Biology, 14(11), 2636-2660. Adams, W. A. (1973). The effect of organic matter on the bulk and true densities of some uncultivated podzolic soils. Journal of Soil Science, 24(1), 10-17. Alef, K., & Nannipieri, P. (1995). Methods in applied soil microbiology and biochemistry. Academic press. Al-Kaisi, M. M., & Yin, X. (2005). Tillage and crop residue effects on soil carbon and carbon dioxide emission in cornsoybean rotations. Journal of Environmental Quality, 34(2), 437-445. Allison, F. E. (1973). Soil organic matter and its role in crop production. Elsevier Scientific Publishing Company. Andrén, O., Rajkai, K., & Kätterer, T. (1993). Water and temperature dynamics in a clay soil under winter wheat: influence on straw decomposition and N immobilization. Biology and fertility of soils, 15(1), 1-8. Angers, D. A., Bolinder, M. A., Carter, M. R., Gregorich, E. G., Drury, C. F., Liang, B. C., ... & Martel, J. (1997). Impact of tillage practices on organic carbon and nitrogen storage in cool, humid soils of eastern Canada. Soil and Tillage Research, 41(3), 191-201. Basche, A. D., Archontoulis, S. V., Kaspar, T. C., Jaynes, D. B., Parkin, T. B., & Miguez, F. E. (2016). Simulating long-term impacts of cover crops and climate change on crop production and environmental outcomes in the Midwestern United States. Agriculture, Ecosystems & Environment, 218, 95-106. Bennie, J., Huntley, B., Wiltshire, A., Hill, M. O., & Baxter, R. (2008). Slope, aspect and climate: spatially explicit and implicit models of topographic microclimate in chalk grassland. Ecological Modelling, 216(1), 47-59. Cambardella, C. A., & Elliott, E. T. (1992). Particulate soil organic-matter changes across a grassland cultivation sequence. Soil Science Society of America Journal, 56(3), 777-783. Cambardella, C. A., & Elliott, E. T. (1992). Particulate soil organic-matter changes across a grassland cultivation sequence. Soil Science Society of America Journal, 56(3), 777-783. Cambardella, C. A., & Elliott, E. T. (1993). Carbon and nitrogen distribution in aggregates from cultivated and native grassland soils. Soil Science Society of America Journal, 57(4), 1071-1076. Cambardella, C. A., Moorman, T. B., Parkin, T. B., Karlen, D. L., Novak, J. M., Turco, R. F., & Konopka, A. E. (1994). Field-scale variability of soil properties in central Iowa soils. Soil science society of America journal,58(5), 1501-1511. 44 Campbell, C. A., & Souster, W. (1982). Loss of organic matter and potentially mineralizable nitrogen from Saskatchewan soils due to cropping.Canadian Journal of Soil Science, 62(4), 651-656. Chan, K. Y. (2001). Soil particulate organic carbon under different land use and management. Soil use and management, 17(4), 217-221. Corre, M. D., Schnabel, R. R., & Stout, W. L. (2002). Spatial and seasonal variation of gross nitrogen transformations and microbial biomass in a Northeastern US grassland. Soil Biology and Biochemistry, 34(4), 445-457. Da Silva, J. M., & Alexandre, C. (2005). Spatial variability of irrigated corn yield in relation to field topography and soil chemical characteristics.Precision Agriculture, 6(5), 453-466. Dabney, S. M., Delgado, J. A., & Reeves, D. W. (2001). Using winter cover crops to improve soil and water quality. Communications in Soil Science and Plant Analysis, 32(7-8), 1221-1250. De Bruin, J. L., Porter, P. M., & Jordan, N. R. (2005). Use of a rye cover crop following corn in rotation with soybean in the upper Midwest. Agronomy Journal, 97(2), 587-598. Dharmakeerthi, R. S., Kay, B. D., & Beauchamp, E. G. (2005). Factors contributing to changes in plant available nitrogen across a variable landscape. Soil Science Society of America Journal, 69(2), 453-462. Dore, M. H. (2005). Climate change and changes in global precipitation patterns: what do we know?. Environment international, 31(8), 1167-1181. Duiker, S. W., & Curran, W. S. (2004). Rye Cover Crop Management IN Corn. In Tillage Conference FOR Sustainable Agriculture (p. 208). Duiker, S. W., & Lal, R. (1999). Crop residue and tillage effects on carbon sequestration in a Luvisol in central Ohio. Soil and Tillage Research, 52(1), 73-81. Dungait, J. A., Ghee, C., Rowan, J. S., McKenzie, B. M., Hawes, C., Dixon, E. R., ... & Hopkins, D. W. (2013). Microbial responses to the erosional redistribution of soil organic carbon in arable fields. Soil Biology and Biochemistry, 60, 195-201. Eijsackers, H., & Zehnder, A. J. (1990). Litter decomposition: a Russian matriochka doll. Biogeochemistry, 11(3), 153-174. Follett, R. F. (2001). Soil management concepts and carbon sequestration in cropland soils. Soil and Tillage Research, 61(1), 77-92. Franzluebbers, A. J. (2002). Water infiltration and soil structure related to organic matter and its stratification with depth. Soil and Tillage Research,66(2), 197-205. 45 Gale, W. J., Cambardella, C. A., & Bailey, T. B. (2000). Root-derived carbon and the formation and stabilization of aggregates. Soil Science Society of America Journal, 64(1), 201-207. Gessner, M. O., Swan, C. M., Dang, C. K., McKie, B. G., Bardgett, R. D., Wall, D. H., & Hättenschwiler, S. (2010). Diversity meets decomposition.Trends in ecology & evolution, 25(6), 372-380. Gupta, S., & Larson, W. E. (1979). Estimating soil water retention characteristics from particle size distribution, organic matter percent, and bulk density. Water resources research, 15(6), 1633-1635. Harrison-Kirk, T., Beare, M. H., Meenken, E. D., & Condron, L. M. (2013). Soil organic matter and texture affect responses to dry/wet cycles: effects on carbon dioxide and nitrous oxide emissions. Soil Biology and Biochemistry, 57, 43-55. Haynes, R. J. (2000). Labile organic matter as an indicator of organic matter quality in arable and pastoral soils in New Zealand. Soil Biology and Biochemistry, 32(2), 211-219. Henriksen, T. M., & Breland, T. A. (1999). Nitrogen availability effects on carbon mineralization, fungal and bacterial growth, and enzyme activities during decomposition of wheat straw in soil. Soil Biology and Biochemistry,31(8), 1121-1134. Hudson, B. D. (1994). Soil organic matter and available water capacity.Journal of Soil and Water Conservation, 49(2), 189-194. Iqbal, J., Castellano, M. J., & Parkin, T. B. (2013). Evaluation of photoacoustic infrared spectroscopy for simultaneous measurement of N2O and CO2 gas concentrations and fluxes at the soil surface. Global change biology, 19(1), 327-336. Jiang, P., & Thelen, K. D. (2004). Effect of soil and topographic properties on crop yield in a north-central cornsoybean cropping system. Agronomy Journal, 96(1), 252-258. Jong, R. D., Campbell, C. A., & Nicholaichuk, W. (1983). Water retention equations and their relationship to soil organic matter and particle size distribution for disturbed samples. Canadian Journal of Soil Science, 63(2), 291-302. Kang, S., Doh, S., Lee, D., Lee, D., Jin, V. L., & Kimball, J. S. (2003). Topographic and climatic controls on soil respiration in six temperate mixedhardwood forest slopes, Korea. Global Change Biology, 9(10), 1427-1437. Kaspar, T. C., Parkin, T. B., Jaynes, D. B., Cambardella, C. A., Meek, D. W., & Jung, Y. S. (2006). Examining changes in soil organic carbon with oat and rye cover crops using terrain covariates. Soil Science Society of America Journal, 70(4), 1168-1177. 46 Kaspar, T. C., Parkin, T. B., Jaynes, D. B., Cambardella, C. A., Meek, D. W., & Jung, Y. S. (2006). Examining changes in soil organic carbon with oat and rye cover crops using terrain covariates. Soil Science Society of America Journal, 70(4), 1168-1177. Kladivko, E. J., Helmers, M. J., Abendroth, L. J., Herzmann, D., Lal, R., Castellano, M. J., ... & Basso, B. (2014). Standardized research protocols enable transdisciplinary research of climate variation impacts in corn production systems. Journal of Soil and Water Conservation, 69(6), 532-542. Kravchenko, A. N., & Bullock, D. G. (2000). Correlation of corn and soybean grain yield with topography and soil properties. Agronomy Journal, 92(1), 75-83. Kunkel, M. L., Flores, A. N., Smith, T. J., McNamara, J. P., & Benner, S. G. (2011). A simplified approach for estimating soil carbon and nitrogen stocks in semi-arid complex terrain. Geoderma, 165(1), 1-11. Kuo, S., & Jellum, E. J. (2000). Long-term winter cover cropping effects on corn (Zea mays L.) production and soil nitrogen availability. Biology and fertility of soils, 31(6), 470-477. Ladoni, M., Basir, A., Robertson, P. G., & Kravchenko, A. N. (2016). Scaling-up: cover crops differentially influence soil carbon in agricultural fields with diverse topography. Agriculture, Ecosystems & Environment, 225, 93-103. Ladoni, M., Kravchenko, A. N., & Robertson, G. P. (2015). Topography Mediates the Influence of Cover Crops on Soil Nitrate Levels in Row Crop Agricultural Systems. PloS one, 10(11). Lal, R. (2004). Soil carbon sequestration impacts on global climate change and food security. science, 304(5677), 1623-1627. Lal, R. (2004). Soil carbon sequestration to mitigate climate change.Geoderma, 123(1), 1-22. Loomis, R. S., & Connor, D. J. (1992). Crop ecology: productivity and management in agricultural systems. Cambridge University Press. Lundquist, E. J., Jackson, L. E., Scow, K. M., & Hsu, C. (1999). Changes in microbial biomass and community composition, and soil carbon and nitrogen pools after incorporation of rye into three California agricultural soils. Soil Biology and Biochemistry, 31(2), 221-236. Moore, I. D., Gessler, P. E., Nielsen, G. A., & Peterson, G. A. (1993). Soil attribute prediction using terrain analysis. Soil Science Society of America Journal, 57(2), 443-452. Moore, T. R., De Souza, W., & Koprivnjak, J. F. (1992). Controls on the sorption of dissolved organic carbon by soils. Soil Science, 154(2), 120-129. Munawar, A., Blevins, R. L., Frye, W. W., & Saul, M. R. (1990). Tillage and cover crop management for soil water conservation. Agronomy Journal,82(4), 773-777. 47 Muñoz, J. D., Steibel, J. P., Snapp, S., & Kravchenko, A. N. (2014). Cover crop effect on corn growth and yield as influenced by topography.Agriculture, Ecosystems & Environment, 189, 229-239. Nascente, A. S., Li, Y. C., & Crusciol, C. A. C. (2013). Cover crops and no-till effects on physical fractions of soil organic matter. Soil and Tillage Research, 130, 52-57. Negassa, W., Price, R. F., Basir, A., Snapp, S. S., & Kravchenko, A. (2015). Cover crop and tillage systems effect on soil CO 2 and N 2 O fluxes in contrasting topographic positions. Soil and Tillage Research, 154, 64-74. Papadopoulos, A., Mooney, S. J., & Bird, N. R. A. (2006). Quantification of the effects of contrasting crops in the development of soil structure: an organic conversion. Soil Use and Management, 22(2), 172-179. Parkin, T. B., & Venterea, R. T. (2010). USDA-ARS GRACEnet project protocols, chapter 3. Chamber-based trace gas flux measurements.Sampling protocols. Beltsville, MD p, 1-39. Parton, W. J., Schimel, D. S., Cole, C. V., & Ojima, D. S. (1987). Analysis of factors controlling soil organic matter levels in Great Plains grasslands.Soil Science Society of America Journal, 51(5), 1173-1179. Plaza-Bonilla, D., Álvaro-Fuentes, J., & Cantero-Martínez, C. (2014). Identifying soil organic carbon fractions sensitive to agricultural management practices. Soil and Tillage Research, 139, 19-22. Raich, J. W., & Potter, C. S. (1995). Global patterns of carbon dioxide emissions from soils. Global Biogeochemical Cycles, 9(1), 23-36. Raimbult, B. A., Vyn, T. J., & Tollenaar, M. (1991). Corn response to rye cover crop, tillage methods, and planter options. Agronomy journal, 83(2), 287-290. Rawls, W. J. (1983). Estimation soil bulk density from particle size analysis and organic matter content. Soil Science,135(2), 123-125. Rawls, W. J., Pachepsky, Y. A., Ritchie, J. C., Sobecki, T. M., & Bloodworth, H. (2003). Effect of soil organic carbon on soil water retention.Geoderma, 116(1), 61-76. Reeves, D. W. (1997). The role of soil organic matter in maintaining soil quality in continuous cropping systems. Soil and Tillage Research, 43(1), 131-167. Reicosky, D. C. (2003). Conservation agriculture: Global environmental benefits of soil carbon management. In Conservation Agriculture (pp. 3-12). Springer Netherlands. Reicosky, D. C., & Forcella, F. (1998). Cover crop and soil quality interactions in agroecosystems. Journal of Soil and Water Conservation,53(3), 224-229. 48 Ritchie, J. C., McCarty, G. W., Venteris, E. R., & Kaspar, T. C. (2007). Soil and soil organic carbon redistribution on the landscape. Geomorphology,89(1), 163-171. Romano, N., & Palladino, M. (2002). Prediction of soil water retention using soil physical data and terrain attributes. Journal of Hydrology, 265(1), 56-75. Rovira, P., & Vallejo, V. R. (2002). Mineralization of carbon and nitrogen from plant debris, as affected by debris size and depth of burial. Soil Biology and Biochemistry, 34(3), 327-339. Sainju, U. M., Singh, B. P., & Whitehead, W. F. (2002). Long-term effects of tillage, cover crops, and nitrogen fertilization on organic carbon and nitrogen concentrations in sandy loam soils in Georgia, USA. Soil and Tillage Research, 63(3), 167-179. Sequeira, C. H., & Alley, M. M. (2011). Soil organic matter fractions as indices of soil quality changes. Soil Science Society of America Journal,75(5), 1766-1773. Snapp, S. S., Swinton, S. M., Labarta, R., Mutch, D., Black, J. R., Leep, R., ... & O'Neil, K. (2005). Evaluating cover crops for benefits, costs and performance within cropping system niches. Agronomy Journal, 97(1), 322-332. Srivastava, M., Sharma, S. D., & Kudrat, M. (2012). Effect of Crop Rotation, Soil Temperature and Soil Moisture on CO 2 Emission Rate in Indo-Gangetic Plains of India. International Journal of Agriculture and Forestry, 2(3), 117-120. Terra, J. A., Shaw, J. N., Reeves, D. W., Raper, R. L., Van Santen, E., & Mask, P. L. (2004). Soil carbon relationships with terrain attributes, electrical conductivity, and a soil survey in a coastal plain landscape. Soil science,169(12), 819-831. Timlin, D. J., Pachepsky, Y., Snyder, V. A., & Bryant, R. B. (1998). Spatial and temporal variability of corn grain yield on a hillslope. Soil Science Society of America Journal, 62(3), 764-773. Trenberth, K. E. (2011). Changes in precipitation with climate change.Climate Research, 47(1-2), 123-138. Tromp-van Meerveld, H. J., & McDonnell, J. J. (2006). On the interrelations between topography, soil depth, soil moisture, transpiration rates and species distribution at the hillslope scale. Advances in Water Resources,29(2), 293-310. Trumbore, S. E., Chadwick, O. A., & Amundson, R. (1996). Rapid exchange between soil carbon and atmospheric carbon dioxide driven by temperature change. Science, 272(5260). Unger, P. W., & Kaspar, T. C. (1994). Soil compaction and root growth: a review. Agronomy Journal, 86(5), 759-766. 49 Van Doren, D. M., Triplett, G. B., & Henry, J. E. (1976). Influence of long term tillage, crop rotation, and soil type combinations on corn yield. Soil Science Society of America Journal, 40(1), 100-105. Wagner-Riddle, C., Gillespie, T. J., & Swanton, C. J. (1994). Rye cover crop management impact on soil water content, soil temperature and soybean growth. Canadian journal of plant science, 74(3), 485-495. Wander, M. M., Bidart, M. G., & Aref, S. (1998). Tillage impacts on depth distribution of total and particulate organic matter in three Illinois soils. Soil Science Society of America Journal, 62(6), 1704-1711. West, T. O., & Post, W. M. (2002). Soil organic carbon sequestration rates by tillage and crop rotation. Soil Science Society of America Journal, 66(6), 1930-1946. Western, A. W., Grayson, R. B., Blöschl, G., Willgoose, G. R., & McMahon, T. A. (1999). Observed spatial organization of soil moisture and its relation to terrain indices. Water resources research, 35(3), 797-810. Wickings, K., Grandy, A. S., & Kravchenko, A. N. (2016). Going with the flow: Landscape position drives differences in microbial biomass and activity in conventional, low input, and organic agricultural systems in the Midwestern US. Agriculture, Ecosystems & Environment, 218, 1-10. Yimer, F., Ledin, S., & Abdelkadir, A. (2006). Soil organic carbon and total nitrogen stocks as affected by topographic aspect and vegetation in the Bale Mountains, Ethiopia. Geoderma, 135, 335-344. Zhu, Q., & Lin, H. (2011). Influences of soil, terrain, and crop growth on soil moisture variation from transect to farm scales. Geoderma, 163(1), 45-54.