Introduction

Essential for agricultural production, but often limiting, mineable reserves of phosphorus (P) are non-renewable and concentrated in regions with territorial conflicts, adding a geopolitical dimension to P scarcity (Cordell and White 2014). Furthermore, P losses via erosion and leaching are responsible for eutrophication of water bodies and ecosystem degradation (Schoumans et al. 2014). Therefore, reliance on costly P-fertilizer inputs poses a threat to food security. Soil-improving cropping systems such as cover crops and conservation tillage are gaining attention for their potential to enhance overall sustainability of agriculture and P management (Tonitto et al. 2006; Simpson et al. 2011; Scopel et al. 2013; Damon et al. 2014). They play an increasingly important role with respect to the concept of ecological intensification (Bommarco et al. 2013) and agroecology (Altieri 2002; Faucon et al. 2017).

Phosphorus-containing fertilizers are used extensively, but the processes underlying the biogeochemical P cycle are thus far not fully understood (Bünemann et al. 2011), especially under conditions of low P availability (Clarkson 1985; George et al. 2018). Phosphorus is present in soils in both mineral and organic forms with vastly different degrees of availability; only very small amounts of inorganic P (Pi) are present in the soil solution (Pierre and Parker 1927), and this is the form taken up by plants. In agricultural soils, the soil microbial community is increasingly acknowledged as the principal driver of soil P dynamics (Bünemann et al. 2011; Richardson and Simpson 2011), and efforts have been made to include soil microbes into P cycling models (Rengel 2008; Hinsinger et al. 2011; Damon et al. 2014). The relatively large pools of soil organic P (Porg), as a result of the combined action of plant- and microbial-exuded carboxylates for mobilization and enzymes for mineralization, constitute a valuable yet poorly understood resource (George et al. 2018; Menezes-Blackburn et al. 2018). The options for their management and their effects on crop P nutrition status and yield are therefore of paramount interest to both agronomists and farmers.

Cover cropping is the practice of growing plants, usually in the off-season, leaving their biomass on the field to provide various benefits for the agroecosystem, including erosion reduction, soil organic matter (SOM) build-up, weed and pathogen control, and nutrient management. Cover crops are also used to improve the P efficiency of added organic or mineral fertilizers by increasing soil biological activity or uptake and protection of soluble mineral P in strongly P-fixing soils (Kamh et al. 1998; Kuo et al. 2005). In principle, most plant species could be used for these purposes, and the list of plant species that can be used as cover crops is rapidly expanding. However, agronomic requirements of the cover crops (e.g., rapid growth and high biomass production, weed suppression, easy termination) impose some restrictions. Plant species in general, and therefore also when used as cover crops, vary greatly in their biomass production, soil exploration, exudation of P-mobilizing and organic P-mineralizing compounds, as well as their interaction with the rhizosphere microbial community. The variety of strategies for P-acquisition employed by different plant types must be considered, as, e.g., non-mycorrhizal species, while highly efficient at P mobilization, do not necessarily interact strongly with soil microbes (Lambers and Teste 2013; Lambers et al. 2015a). The inclusion of cover crops with special properties such as, e.g., the increase of the inoculum potential by arbuscular mycorrhizal fungi, can be beneficial in agricultural rotations that use domesticated cash crops developed in high-input breeding systems (Plenchette et al. 2005). Conservation agriculture, defined by the combination of cover crops, conservation tillage, and an adequate crop rotation (Hobbs et al. 2008), decreases labor intensity and frequency, enhances soil rest, and benefits soil biota. Here, crop residues are mixed less deeply into the soil than under conventional tillage, modifying soil biological parameters and mineralization dynamics. Other effects of management could result from the termination method used (e.g., spraying or roll-chopping) for winter-hard cover crops (Creamer and Dabney 2002).

The effects of cover crops and conservation tillage on crop yield and soil properties, especially nitrogen (N) dynamics, have been the subject of many studies, several reviews (Dabney et al. 2001; Dreymann et al. 2005; Tonitto et al. 2006; Dahlin and Stenberg 2009), and projects (Crossland et al. 2015). However, P dynamics have, until recently, rarely been addressed. Some reviews on the underlying mechanisms of plant and microbial P acquisition and resulting implications for agricultural management have provided a theoretical foundation for predicting the influence of cover crops on P dynamics (Horst et al. 2001; Guppy et al. 2005; Richardson and Simpson 2011; Richardson et al. 2011; Damon et al. 2014). Interestingly, those publications have largely reflected their regional conditions; climate, soil P content and sorption capacity. Studies of P-fertilization efficiency, often associated with pastures, come mainly from Australia (Rose et al. 2010a; Simpson et al. 2011; McLaughlin et al. 2011; Faucon et al. 2015). Brazilian studies have emphasized conservation agriculture and acid soils with poor P availability (LeMare et al. 1987; Calegari et al. 2013; Balota et al. 2014; Fageria et al. 2016; Varela et al. 2017), while Scandinavian researchers have concentrated on P leaching (Liu et al. 2015; Aronsson et al. 2016). In the USA, cover crops of different species have been investigated (Lal et al. 1978), including their effects on mycorrhizal fungi (Galvez et al. 1995; Zibilske and Makus 2009; Rick et al. 2011; Maltais-Landry 2015), whereas in China and India, studies have often focused on microbial inoculants (Devi et al. 2013; Cui et al. 2015). In some African countries, India, and Mexico, with their prevalent traditional smallholder cropping systems, agroforestry and intercropping have been important topics for research (LeMare et al. 1987; Tarawali et al. 1999; Dinesh et al. 2004; Sileshi et al. 2008; Castillo-Caamal and Caamal-Maldonado 2011; Devi et al. 2013; Tanwar et al. 2014; Parihar et al. 2016). Previous meta-analyses have been concerned mainly with the effects of cover crops on Zea mays (maize) yield in North America (Miguez and Bollero 2005), effects on soil properties and yield in the South American Pampas region (Alvarez et al. 2017), response to woody and herbaceous legumes in sub-Saharan Africa (Sileshi et al. 2008), soil organic carbon (Poeplau and Don 2015) as well as Nitrogen dynamics (Tonitto et al. 2006). P nutrition and dynamics have not been a focus of the meta-analyses in these agro-ecosystems. Consequently, understanding whether and how cover crops can benefit the P nutrition of following main crops is sorely needed.

The aim of our review is to bridge our present knowledge of soil-plant-microbe interactions with the potential of cover crops to stimulate P dynamics in agricultural ecosystems. We begin with a description of P pools and P-acquisition mechanisms of cover crops, and present a conceptual framework of how P dynamics of cover crops and main crops may be linked. Both plants and soil microorganisms are involved in P dynamics of agro-ecosystems; we therefore present the most important mechanisms and pathways for both with respect to cover crops. We reviewed the conceptual framework through an extended meta-analysis that included 240 datasets derived from 25 studies. The focus of the meta-analysis was to elucidate whether

  1. (1)

    P acquisition by cover crops stimulates growth and P uptake of different main crops;

  2. (2)

    cover crops enhance mycorrhizal colonization, short-term storage of P in soil microorganisms, and P mineralization, improving plant- and microbial-driven P uptake of the main crop;

  3. (3)

    site conditions (e.g., fractions of available P) modify P benefit to the main crop; and

  4. (4)

    cover crop management (i.e. species mixtures, tillage intensity, and fertilization) can be used to increase P benefits of the main crop.

Availability of phosphorus

Plant P nutrition is constrained by limited availability, due to physicochemical processes in the soil, of orthophosphate, the form that is taken up by roots. Phosphorus compounds interact strongly with the soil through sorption to particle surfaces (including SOM), slow diffusion into aggregates, and formation of precipitates with cations of calcium (Ca), as well as sorption onto oxides and hydroxides of iron (Fe) and aluminium (Al) (Kelly and Midgley 1943), under alkaline and acid conditions, respectively. These interactions result in low P availability and P-fertilizer efficiency (McLaughlin et al. 2011). Soils with low P-sorption capacity also exhibit the associated hazard of P losses via leaching, whereas P-sorbing soils have problems with P-fertilizer efficiency due to immobilization of added P. As P-fertilizer efficiency is only 10-20% in the short term (Chien et al. 2011), agricultural soils in industrialized countries have commonly received excessive loads of P over decades, often without reaching the soil saturation limit. This accumulated P constitutes a valuable resource that could be accessed by employing appropriate cropping systems. For a recent review about pools of recalcitrant P in agricultural soils and opportunities for mobilization, see Menezes-Blackburn et al. (2018). On the other hand, highly-weathered soils in the tropics are P impoverished, and rich in Fe/Al oxides and hydroxides (Simpson et al. 2011). In these systems, efficient recycling and use of the available inputs is essential. Phosphorus inputs into soil are from weathering of P-containing minerals or mobilization/mineralization of other P pools of low availability, as well as atmospheric deposition by dust, and by fertilizer application. Phosphorus is removed from the system through biomass of harvested crops, by erosion, and by leaching, and is accompanied by accumulation of P forms of low availability (Condron et al. 2005). Erosion, globally the biggest threat to sustainable soil stewardship (Bernoux et al. 2006; Durán Zuazo and Rodríguez Pleguezuelo 2008), is one of the major losses of agricultural P, as it is often associated with particle fractions prone to transport.

Soil P consists of both large, but stable and small, but highly-dynamic pools (Sharpley 1995). Although in heavily fertilized agricultural soils labile Pi may temporarily dominate the plant-available pools (Negassa and Leinweber 2009), usually 30-65% of total P (Pt) is present in organic forms (Harrison 1987; Condron et al. 1990), and more in soils with high SOM content (Borie and Zunino 1983). During pedogenesis, the primary P-bearing minerals are slowly depleted, with highly weathered soils containing almost exclusively occluded P and Porg, which is tightly recycled in the biomass (Smeck 1985). Organic P can contribute to crop nutrition, as plants and microbes can access the more labile Porg by a combination of mobilization with carboxylates and subsequent enzymatic mineralization (Condron et al. 2005; Richardson and Simpson 2011). In addition to phospholipids and nucleic acid-P (both forms account for less than 10% of total soil Porg, but comprise most of the microbial P (Pmic)), inositol phosphates (phytate) accumulate in soil due to their stable nature, and represent the major fraction of Porg (Jones and Oburger 2011). The inositol molecules consist of 1-6 phosphates attached to a C6-ring with ester-bonds, requiring specialized enzymes for breakdown; they also interact strongly with the soil due to their high charge density (Turner et al. 2002; Turner 2007). Due to low substrate availability, low phytase production by roots, and low enzyme-substrate efficiency in soils, only some plant species can access phytate to a limited extent (Menezes-Blackburn et al. 2012). However, several plant species are able to grow with sodium (Na)-hexaphytate as their sole P source in the laboratory (Steffens et al. 2010). In natural soils with endogenous phytates, mobilization via carboxylates and subsequent interactions with microbes appear necessary for plants to use this resource. For a full review of inositol phosphates, see Turner et al. (2007).

The direction of the effect of long-term SOM accumulation by cover crops on P availability is not clear. Improved soil physical-chemical parameters (e.g., water-holding capacity, aggregate stability) (Dorado et al. 2003) may increase P availability to crops directly or indirectly (Eichler-Löbermann et al. 2008). Some fractions of SOM may compete with P for binding sites on particle surfaces, decreasing P-sorption capacity (Janegitz et al. 2013). However, large amounts of Porg are in the form of uncharacterized high-molecular-weight organic material (McLaren et al. 2015), linking the accumulation of poorly-available Porg closely to SOM dynamics. This could constitute a problem for the exploitation of these pools (Romanyà et al. 2017), as SOM is needed for soil structure, fertility, and climate change mitigation, and one target of cover cropping is to increase SOM content in soils. However, the fact that C:N:Sulfur (S) ratios are relatively constant across soils, whereas C:P and C:Porg are more variable, may allow increasing available P (Pa) through Porg mineralization without affecting SOM accumulation. A P-priming effect has been described (Randhawa et al. 2005), but remains unresolved due to methodological constraints (Damon et al. 2014).

The most common way to characterize soil P is by using variants of the Hedley fractionation, which determines the amount of Pi and Porg in various soil extracts (i.e. NaHCO3, NaOH, HCl), which are used to represent pools of differing degrees of availability (Hedley et al. 1982; Cross and Schlesinger 1995). However, the net contribution of these operational fractions to uptake by microbes and plants, and therefore the validity of the method for the prediction of plant P uptake, is not straightforward (Negassa and Leinweber 2009; Rose et al. 2010b). Olsen-P (Olsen et al. 1954), using NaHCO3 as extractant, is one of the most widely used methods and often correlates well with yield and P uptake, but was originally developed for calcareous soils. There are many other extractants for soil P testing methods, including Mehlich-III (Mehlich 1984), Bray 1 (Bray and Kurtz 1945), water (van der Paauw 1971), calcium-ammonium-lactate (CAL) (Schüller 1969) and Colwell (1963). Their application differs, even among regions in the same country, due to prevailing soil characteristics, but also due to historical reasons. More recently developed test methods that show promising results such as resin-P and “diffusive gradients in thin-films” (DGT) (Mason et al. 2013) are not yet widely used. One issue is that the color methods used in routine soil testing, mostly molybdate-blue (Murphy and Riley 1962), do not account for Porg (Steffens et al. 2010). This has been justified based on the assumption that Porg seems to play a minor role in plant nutrition under high availability of Pi (Guo et al. 2000), but has consequences for systems dependent on Porg (Dao et al. 2015). The combination of imperfect P-test methods and substantial knowledge gaps in understanding the complex P dynamics in soils constricts a scientific elaboration of general agricultural recommendations (Turner et al. 2005). Another problem is the definition of “plant-available P”, as there are major differences among plant species (Lambers et al. 2006) and even crop varieties (Pang et al. 2018a) regarding their ability to access different soil P pools.

Phosphorus acquisition by cover crops

Plant species differ greatly in their P-acquisition strategies. The ability of cover crops to access poorly-available soil P is often superior to that of cash crops (Fig. 1). Their strategies can be summarized as: 1) exploration of a greater soil volume by an adaptive root architecture and root morphology; 2) mobilization of sparingly-soluble inorganic and organic P forms; and 3) mineralization of Porg. These mechanisms may all be enhanced through interaction with soil microbes.

Fig. 1
figure 1

Strategies and mechanisms for phosphorus (P) acquisition by plants: 1) soil exploration via roots and mycorrhizal hyphae; 2) mobilization of sparingly-soluble inorganic P (Pi) and organic P (Porg) by exudation of H+/OH- and carboxylates; 3) mineralization of Porg by phosphatases. Plant-driven processes have solid outlines, microbial activity is shown by dotted outlines

In addition to the transfer of Pa via cover crop residues to the main crop and to chemical rhizosphere modifications, some studies attribute benefits of cover crops to the subsequent crop to the soil microbial community, which influences P dynamics, both during the cover crop and the main crop phases (Nuruzzaman et al. 2005a; Pypers et al. 2007; Rose et al. 2010a; Mat Hassan et al. 2013).

Plants interact with the soil microbial community by releasing organic compounds into the rhizosphere that serve as substrates and signaling molecules to the microbes, increasing their abundance and activity several-fold (Bünemann et al. 2004; Balota et al. 2014). Our current knowledge indicates that both plants and soil determine microbial community composition (Marschner et al. 2001; Nannipieri et al. 2008). Together with plant roots, microbes are the principal drivers of Pi and Porg solubilization and of the mineralization of more or less recalcitrant Porg in soil. They possess a diverse array of mechanisms to increase P acquisition by plants, including modifications and extension of root systems, allowing roots to access P-rich substrates that are otherwise unavailable to plants due to their location in the soil (i.e. in narrow pore spaces).

Root architecture

The differences in P uptake by cover crops are determined partly by their root architecture, with topsoil exploration and root hair density the most important traits for improved P uptake (Richardson et al. 2011). The kinetic properties of the Pi-uptake system, unlike those of more mobile nutrients such as nitrate, are not a major rate-limiting step in plant P acquisition (Clarkson 1985; Barber 1995). Mycorrhizal fungi play a fundamental role as extensions of the roots, whereas other microorganisms promote root growth and modify root architecture (branching, root hairs) via signaling molecules in the rhizosphere (Hayat et al. 2010). Cover crops with more extensive root systems scavenge P from a larger and deeper soil volume, and make it potentially available for main crops with shallow roots (Dube et al. 2014). Some plants (e.g., Lupinus species), possess specialized root structures, termed cluster roots, that exploit soils with low P availability and potentially enhance P availability for the main crop (Nuruzzaman et al. 2005b; Lambers et al. 2006). They may additionally facilitate P acquisition for neighboring plants (Gardner and Boundy 1983; Horst and Waschkies 1987; Cu et al. 2005).

To date, the most extensive description of cover crop traits related to P uptake is that by Wendling et al. (2016), who classified cover crops into five groups based on shoot biomass and nutrient concentration, comprising species from different families. The main findings of this study were shoot and root traits, rather than taxonomy; species with high nutrient concentrations and high root length density were recommended under high-fertility conditions and from a short-term perspective. However, although biochemical and microbial root P-acquisition strategies were not assessed in this study, they may well have been relevant, especially in systems with low P availability.

Phosphorus mobilization

Soil P-mining strategies enhance desorption and solubilization of sparingly-available Pi and Porg pools, which often limit P availability. Plants and microbes are capable of exuding low-molecular-weight organic anions (carboxylates) to dissolve precipitates and chelate metal cations, both of which make phosphate unavailable; because carboxylates facilitate the release of sorbed P via ligand-exchange reactions (Hinsinger 2001), and block binding sites on soil particles, they increase the concentration of P in solution (Ohno and Crannell 1996). The pH of the soil solution is modified by exudation of H+ or OH-/HCO3-. This, in turn, determines the variable surface charge of minerals and SOM, and may also increase P in solution. Proton release enhances P availability only in calcareous soils, due to the dissolution of Ca-phosphate. However, as the (bio)chemistry of P in soils is very complex, with many processes occurring simultaneously, sometimes in opposite directions, it is difficult to predict the effect of pH changes on P dynamics. A review of Pi bioavailability in the rhizosphere was written by Hinsinger (2001). Despite their name, ‘organic acids’ do not substantially decrease rhizosphere pH, as they are mostly released as organic anions, generally with cations other than protons as balancing ions (Zhu et al. 2005; Roelofs et al. 2008). Exudation of P-solubilizing organic anions differs strongly among plant genotypes and soils (Kamh et al. 1998; Nuruzzaman et al. 2006) and is related to P deficiency and Al toxicity (Richardson et al. 2011). Citrate is a commonly released organic anion and one of the most effective for P mobilization (Jones 1998). It is produced in large quantities by Lupinus albus (Gardner et al. 1983; Dinkelaker et al. 1989; Cu et al. 2005) and other legumes (Kamh et al. 2002), including Cicer arietinum (chickpea) (Veneklaas et al. 2003), Vicia faba (faba bean) (Li et al. 2007), and Trifolium pratense (red clover) (Gerke and Meyer 1995), but also by Brassica napus (canola) (Hoffland et al. 1992). The strategy of dicots (e.g., Fabaceae, Brassicaceae) is to utilize biochemical rhizosphere modification for P mobilization, whereas Poaceae predominantly take up P using their extensive root systems (Maltais-Landry 2015; Schnug and De Kok 2016). The exudation of different organic anions and acidification may be complementary (Gerke and Meyer 1995), but the mechanisms are complex and the result depends strongly on soil chemistry and P level (Oburger et al. 2011). Complicating the system further, microorganisms function as both potential sinks and alternative sources of carboxylates (Deubel et al. 2000), and soil fauna remobilize P from the microbial biomass (Hinsinger et al. 2015). The identification of pH as the principal driver of microbial diversity in soils (Philippot et al. 2009) paired with substantial changes in pH in the rhizosphere led Hinsinger et al. (2009) to the hypothesis that root-induced pH changes shape the structure of the rhizosphere microbial community equally or more importantly than root C deposition.

Microbes may use plant exudates to produce P-solubilizing compounds in the rhizosphere, complementing P mobilization by roots (Schilling et al. 1998). In fact, some root exudates do not function directly in plant nutrient acquisition, but are composed of mobile sugars, which can be used by the rhizosphere microbial community. The critical role of desorption for the mineralization of Porg is receiving increasing attention (Giaveno et al. 2010). Pseudomonas species are among the most frequently studied P-solubilizing bacteria, but also species of Burkholderia, Enterobacter, Pantoea, Bacillus solubilize P (Jorquera et al. 2008). Together with arbuscular mycorrhizal fungi (AMF), P-solubilizing microorganisms are a target for the development of microbial inoculants, although often with limited success in the field (Parray et al. 2016).

The potential of some plants to mine sparingly-available P pools led to optimism about their potential to increase P availability when used as cover crops (Teboh and Franzen 2011; Boglaienko et al. 2014). However, under conditions of P-deficiency, plants do not always respond by releasing organic anions and their effect on P uptake is not consistent (Wang et al. 2016).

The minor direct benefits of cover crop rhizosphere modification (i.e. carboxylate exudation) for subsequent crops (Possinger et al. 2013) may be explained by the short duration of carboxylate exudation (a few weeks) associated with legume roots, limiting their effects on the main crop (Nuruzzaman et al. 2005a). Notwithstanding, the binding of carboxylates to Fe/Al hydroxides could delay microbial mineralization (Jones and Edwards 1998) and reduce soil P-sorption capacity. The ability of microorganisms to access sparingly-available P with benefits to plants may depend on microbial turnover (Richardson and Simpson 2011), as the microbial biomass conserves solubilized P.

The inclusion of carboxylates and other rhizosphere processes could improve plant nutrition models, which have to date failed to predict the actual uptake of P and other low-mobility nutrients, especially under low-input conditions (Hinsinger et al. 2011). Leaf manganese (Mn) concentration, for example, can be used as a proxy for carboxylate concentration in the rhizosphere, providing a tool that may be more reliable than measurement of rhizosphere carboxylates, given their transient nature in the field (Lambers et al. 2015b; Pang et al. 2018a).

Mineralization

In addition to mobilization mechanisms described above, the considerable amounts of Porg in soil require, once in solution, enzymatic hydrolysis to become plant-available Pi. In pot experiments, the activity of phosphatase enzymes was three-fold and nine-fold greater in the rhizosphere of Triticum aestivum (wheat) and Lupinus albus (white lupin), respectively, than in bulk soil, and this increased activity was, in both cases, associated with the depletion of soil Porg (Nuruzzaman et al. 2006). However, phosphatases also have a role in recycling P inside cells and recapturing Porg lost from roots or microbial cells (Tarafdar and Jungk 1987; Barrett-Lennard et al. 1993). Experiments with transgenic Trifolium repens have shown that, under laboratory conditions, transgenic expression of phytase and purple acid phosphatase genes from Medicago truncatula increased the plants’ ability to utilize organic P in response to P deficiency (Ma et al. 2009). The use of these techniques is convenient for experimentation, but in the field, the efficacy of single exudation traits appears to be limited in P-deficient soil conditions where the soil does not exactly match the functional requirements of the enzymes of interest (Giles et al. 2017).

Extracellular enzymes interact strongly with soil particles, leading to adsorption and inactivation, but also to protection against degradation (Rao et al. 2000). Adsorption depends on the mineral composition of the soil (Ditterich et al. 2016) and characteristics of the SOM. There are indications that carboxylates may serve a dual role of desorbing P and providing a favorable pH for the phosphatase enzymes, increasing enzymatic activity (Furutani et al. 2017)

Due to sorption processes, the effect of plant-derived phosphatases will be restricted to a few millimeters of distance from the roots. However, mobile rhizodeposited sugars penetrate further into the bulk soil, and can be used by microbes to produce phosphatases, extending the range of Porg mineralization around the roots. Due to complex interactions with soil, microbial degradation, and interception of the products, increased phosphatase activities do not necessarily translate into a more rapid P-uptake rate for plants.

Under natural conditions, the microbial contribution to the mineralization of Porg in the rhizosphere is undisputed; however, it is often difficult to separate the origin of enzymatic activity, as some enzymes, such as acid phosphatases and some phytases, are produced both by plants and microbes (Nannipieri et al. 2011). Diesterases can be produced by plants also, but diesterase activity is mainly related to microbial biomass (Turner and Haygarth 2005; Lang et al. 2017). Alkaline phosphatases, however, some phytases (Azeem et al. 2015) and phosphonate hydrolases (Hunter et al. 2014) are produced only by microbes. Zymography is a promising in situ method for analysis of the two-dimensional distribution of enzymatic activity in soil. In an experiment with Lupinus albus it was combined with 14C imaging, revealing that alkaline phosphatase-producing microorganisms were not dependent on recent rhizodeposition, whereas acid phosphatase activity was concentrated in the direct vicinity of the roots (Spohn and Kuzyakov 2013a). The relative contributions of microbial groups to the activities of the different phosphatases requires further investigation (Turner and Haygarth 2005). In contrast with plants, which take up P exclusively as Pi, microorganisms may be able to take up low-molecular-weight Porg, and protozoa can make use of high-molecular-weight Porg (Jones and Oburger 2011).

The capacities of a plant species to solubilize and mineralize Porg forms may be related to its rhizosphere-associated microbes. For example, the pasture plants Lolium perenne (perennial ryegrass) and Trifolium repens (white clover) have predominantly phytate-mineralizing bacteria in their rhizospheres, whereas in the cereal crops Avena sativa (oat) and Triticum aestivum (wheat), P-solubilizing bacteria dominate in the rhizosphere. Conversely, Lupinus luteus (yellow lupin) shows the lowest proportion of both bacterial types in a Chilean volcanic soil (Jorquera et al. 2008). Larger quantities of phytate-mineralizing and P-solubilizing fungi can be isolated from the rhizosphere of leguminous crops as compared with those of cereals (Gaind and Nain 2015). The ecological interactions between r-strategists in the rhizosphere and K-strategists in the bulk soil may also influence P mineralization (Hunter et al. 2014).

Strategies used by plants vary. There is a rather microbe-independent strategy, as in Lupinus albus, which releases inhibitors that prevent microbial degradation of root-derived carboxylates and phosphatases, intensively changing the chemistry of a small volume of soil around the cluster roots (Weisskopf et al. 2006). However, there are also non-mycorrhizal Brassicaceae with high levels of rhizodeposition, and mycorrhizal plants that scavenge P from a greater soil volume at a lower intensity. Therefore, the P-acquisition strategy of the cover crop influences the mechanisms of P-benefit to the main crop. A selection of common cover crop species and their properties is described in Table 1.

Table 1 Plant species widely used for cover cropping and their properties as described in the literature

Conceptual framework: how do cover crops affect P dynamics?

Among the P dynamics affected by the soil-plant-microbe processes of cover cropping (Fig. 2), the most studied mechanism is the direct uptake of Pa by plants and the transfer of P within the cover crop biomass. The storage of substantial quantities of P, both by high biomass and high P concentrations, is methodologically relatively simple to assess. Phosphorus in the plant biomass is protected from sorption onto the soil (Groffman et al. 1987) or losses by erosion and leaching, but P mineralization needs to be in synchrony with the needs of the main crop. Some cover crops act through positive effects on the soil microbial community (i.e. earlier mycorrhizal colonization, production of enzymes and increased microbial P as a pool for plants), increasing the capacity of the crop-rhizobiome-system to access P (Oberson et al. 2006, Njeru et al. 2014). There is also chemical modification of the rhizosphere, via changes in pH, carboxylate exudation, or phosphatase release, as described above. Biochemical P mobilization would be potentially greatest in soils with a high content of poorly-available P, under the condition that chemical modification of the rhizosphere persists well into the main cropping phase. All these processes occur simultaneously with differing degrees of relative importance depending on the combination of agroecosystem and management.

Fig. 2
figure 2

Pathways of phosphorus (P) transfer and plant-microbial processes affecting P availability by cover cropping. 1) Soil P pools of varying degrees of availability are solubilized and/or mineralized and are immobilized in the microbial biomass. 2) The microbial biomass releases P into the soil solution which 3) ends up in the plant via root or mycorrhizal uptake. Cover crops may additionally possess the capacity to mine P from poorly-available P pools or to produce biochemical rhizosphere modifications to increase P availability. 4) The roots release rhizodeposits that shape the microbial community, eventually leading to increased P mining. 5) The P stored in the cover crop biomass is transferred to the main crop via cover crop residues, which are decomposed by the soil microbial community (6). The soil microbial community (i.e. mycorrhizal fungi) in the main crop phase, enhanced by the cover crops, may possess an increased capacity to mine P for the main crop (7).

The meta-analysis

In order to analyze the general effects of cover cropping on main crop performance in terms of P nutrition, we conducted a meta-analysis. We also assessed more specific effects, such as the multiple ways in which cover crops, interacting with microbes, influence P dynamics and P uptake of the main crop, as well as different cover crop-main crop combinations.

An initial search in 2017 for online available publications using Scopus with the key-words (”phosphorus“ AND "cover crop" OR "green manure" OR "catch crop") yielded 638 matches that were screened by title and abstract. The literature cited in the studies meeting our criteria was also screened, and we expanded the search further using Google Scholar. We selected those studies that reported the effects on main crop yield and P uptake/P concentration, soil P and/or soil biological parameters related to P cycling (phosphatase activity, microbial biomass P, or abundance of AMF) and included a control treatment without cover crops. Phosphorus-mobilizing carboxylates are rarely measured in field studies and could not be included in the meta-analysis. We used only studies with cover crops and main crops grown in rotation, excluding intercropping or living mulch. Greenhouse experiments were excluded, as were agroforestry and grassland studies. Soil biological properties and available P were determined after termination of the cover crop or during growth of the main crop. Experimental factors such as main crop species and/or other factors (e.g., soils, tillage) and data from different years were treated as separate experiments within a study.

The soils included in this meta-analysis were classified according to their P availability, using the descriptions of the field experiments and the results of standard P tests. Datasets from a single field experiment that had been published in several articles (e.g., yield and soil microbiology in different papers) were merged into a single dataset when possible. Details of the studies used (Weerakoon et al. 1992; Medhi and Datta 1996; Boswell et al. 1998; Vanlauwe et al. 2000; Kabir and Koide 2002; Somado et al. 2003; Jensen et al. 2005; Rutunga et al. 2008; Wang et al. 2008; Eichler-Löbermann et al. 2008; Oikeh et al. 2008; Takeda et al. 2009b; White and Weil 2010; Buyer et al. 2010; Rick et al. 2011; Tiecher et al. 2012a, 2012b; Karasawa and Takebe 2012; Njeru et al. 2014; Balota et al. 2014; Maltais-Landry et al. 2015; Karasawa and Takahashi 2015; Mbuthia et al. 2015; Ro et al. 2016; Pavinato et al. 2017) and the extracted data can be found in supplementary material (S1+S2). The data were extracted from the publications using the shareware-tool DataThief III (Tummers 2006) and the open source software Tabula (Aristarán et al. 2017).

We used main crop yield and main crop P uptake as response variables to evaluate the effect of cover crops, because yield is ultimately of interest to farmers. We decided against P concentration, because there is a trade-off between yield and nutrient concentration; that is, high-yielding cropping methods may decrease the concentration of some minerals in the crop (Garibay et al. 1997), whereas a crop with poor field emergence may have a high P concentration. Cover crop biomass and cover crop P uptake were also evaluated to characterize the different cover crop families. To assess the interactions between cover crops, soil microbial community, and the main crops, AMF abundance/colonization, Pmic, and phosphatase activity (alkaline and acid phosphomonoesterase, phosphodiesterase) were first treated as response variables to determine whether or not they were affected by cover cropping, then included in a separate analysis as moderating variables to determine their influence on main crop performance.

The categorical variables soil P availability (high vs low) and climate (tropical vs temperate), and the agronomic factors tillage (inversion tillage vs non-inversion tillage/no-tillage), fertilization (P-fertilized vs unfertilized), and cropping system (conventional vs organic) were also used as moderating variables for the response variables.

The models had the following basic structure with their respective response and moderating variables:

$$ {\log}_e\left(\widehat{y}\right)=m\ \left( fixed\ \right)+ study\ (fixed)+ study: experiment\left( nested, random\right) $$
y :

response variable

m :

moderating variables

The response variables were all loge-transformed to account for different units and scale effects, but back-transformed and reported as percentage change relative to the respective control treatments for graphical visualization. To calculate the relative percentage change by cover crops, the following formula was used:

$$ Y\left(\%\right)=\frac{{\widehat{y}}_{c{over\ crops}_j}-{\widehat{y}}_{control_j}}{{\widehat{y}}_{control_j}}\ast 100 $$
\( {\widehat{y}}_{{cover\ crop}_j} \) :

modeled median or 95% CI, respectively, of the jth cover crop type

\( {\widehat{y}}_{control} \) :

modeled median of the control treatments corresponding to j

Bearing in mind possible interactions, cover crops and main crops were aggregated into phylogenetic families of similar properties (Table S2). Fabaceae and Poaceae cover crops were by far the most studied groups. Lupinus sp. were not included in the Fabaceae group due to the special P-mobilizing properties of this non-mycorrhizal genus (Lambers et al. 2013). Phacelia (Hydrophyllaceae, only tested in one study) was included in the Asteraceae family, due to similarities in their respective mycorrhizal competence and biomass production. Despite promising results, cover crop mixtures, usually consisting of a Poaceae and either a Brassicaceae or a Fabaceae, were seldom assessed.

Linear mixed models with study as fixed effect and the interaction of study and experiment as random effect were fitted using the package lme4 v1.1-15 (Bates et al. 2015) in R v3.4.3 (R-Core Team 2013) and R-Studio v1.1.423 (RStudio 2013). Graphs were produced with the packages ggplot2 v2.2.1 (Wickham 2009) and cowplot v0.9.2 (Wilke 2017) with estimates from emmeans v1.1 (Lenth 2018) and percentages calculated with plyr v1.8.4 (Wickham 2011). As variance or related parameters were not reported in several studies, the observations were weighted by the number of replicates in each experiment with the weights-statement in the lmer function (all studies had a balanced design). Different models were compared using ML estimation, whereas the final models were fitted with REML. The structure of the fitted models and the F-tests obtained with the package lmerTest v2.0-36 (Kuznetsova et al. 2017) were provided in Table S3 using sjPlot (Lüdecke 2018), sample R code in Table S4.

There was a large variance in and among the published studies due to differences in climate, site conditions, experimental set-ups and management, but also substantial intra-study heterogeneity. Because of missing factor combinations (e.g., not all cover crops were grown with fertilization or tillage), we encountered some difficulties in accounting for interactions of factors. In some cases, we opted to use models with fewer interactions and a higher punctuation by Aikaike’s Information Criterion (AIC) in order to be able to use a greater part of the dataset, under the condition that this did not substantially distort the model output. For the same reason, the moderating variables were tested in separate models. Main crop yield and P content do not represent exactly the same dataset, because not all studies reported both variables. In the studies with wetland rice, only Fabaceae were used as cover crops, so the yield and P uptake of this main crop were calculated separately.

Results of the meta-analysis

Cover crop biomass and P content

From the analysis of the aggregated data from the studies included in the meta-analysis it can be seen that the selection of cover crop was a relevant factor. Cover crop type determined the biomass produced and the tissue P concentration (Fig. 3, Table S3.1: Models 1.1-1.3). Biomass and P concentration were not correlated, resulting in differing C:P ratios: Poaceae cover crops produced the most biomass, but had the lowest P concentration; Polygonaceae had the lowest biomass and Fabaceae and Brassicaceae had rather high P concentrations.

Fig. 3
figure 3

Cover crop dry matter biomass [t ha-1], shoot phosphorus (P) content [kg ha-1] and concentration of P in biomass [g kg-1]. The points represent the modeled median (+/- 95% CI) of the different cover crop treatments. On the left are displayed the number of observations. The letters indicate significant differences among cover crop types with a Tukey post-hoc test (p < 0.05). The corresponding models can be found in Table S3.1

Crop rotation

The integration of cover crops into crop rotations generally increased main crop yields (Fig. 4; Table S3.2: Models 2.1 and 2.2). Main crop yield benefit was determined by main crop species, cover crop type, and their interaction. Maize was most responsive to cover cropping. Other main crops (i.e. Glycine max (soybean) and cereals) tended to respond positively, but the increases were not significant. Wetland rice yields were significantly enhanced by Fabaceae cover crops (Table S3.2: Model 2.3). Brassicaceae, vegetables, and cotton (aggregated as other main crops) were tested in few studies with little response to cover cropping. The interaction cover crop type x main crop type improved the model significantly, although the F-test was not significant. Main crop P uptake was closely related to yields (Fig. S5).

Fig. 4
figure 4

Change in main crop yield and shoot biomass following cover crops from different families. The points represent the modeled median (+/- 95% CI), relative to the respective controls. On the left are displayed the number of observations. The lower-case letters indicate, for a single main crop type with a Tukey post-hoc test (p < 0.05), significant differences among cover crop types (including the control) and the upper-case letters between cover cropping in general and the controls. The corresponding models can be found in Table S3.2

The effect of cover cropping varied also at the species level, as shown for Fabaceae (Fig. 5; Table S3.3: Model 3), with the tropical legumes Lablab purpureus and Mucuna pruriens resulting in the greatest yield increases. Lupinus sp. performed intermediately among the Fabaceae, but the yield increases were not significant.

Fig. 5
figure 5

Change in main crop yield and shoot biomass after different Fabaceae cover crop genera and species. The points represent the percentage change of the modeled median (+/- 95% CI), relative to the respective controls. On the left are displayed the number of observations. The letters indicate significant differences among cover crop species with a Tukey post-hoc test (p < 0.05). Species with only one or two observations were aggregated as “other Fabaceae”: Anthyllis vulneraria, Tephrosia purpurea, Stylosanthes guianensis, Pueraria phaseoloides, Psophocarpus tetragonolobus, Mucuna cochinchinensis, Dolichos lablab, Cassia tora, Canavalia ensiformis, Cajanus cajan, Trifolium sp., Arachis hypogaea. The corresponding model can be found in Table S3.3

Soil biological variables

To understand the mechanisms by which cover crops might stimulate P cycling and yield of the main crop, we explored the soil microbial community after application of different cover crop types. The effects on abundances of AMF, Pmic, and on extracellular P-cycling enzymes (phosphatases) were tested (Fig. 6, Table S3.5: Models 5.1-5.8). Datasets for the soil biological variables included data from seven studies with 60 observations for mycorrhizal abundance, four studies with 53 observations for phosphatase activity, and two studies with 30 observations for Pmic. Abundance of AMF spores and root colonization increased after mycorrhizal cover crops (cover crop mixtures, Fabaceae and Poaceae), but did not change or increased only slightly after non-mycorrhizal cover crops (Brassicaceae and Lupinus sp.). Tillage did not significantly decrease mycorrhizal abundance in the present dataset. Cover cropping generally increased Pmic significantly; with Poaceae, Fabaceae, and Lupinus sp. resulting in the greatest increases, around 25%, but only the effect of Poaceae was significant. Microbial biomass P showed no relationship with main crop yield or P uptake (data not shown). Extracellular phosphatase activity increased around 20% after cover cropping, with Brassicaceae treatments tending to result in the smallest increases over the control, and with Fabaceae, lupins, and Poaceae having the largest effect. Phosphatase activity did not affect main crop growth performance (data not shown).

Fig. 6
figure 6

Soil biological parameters: change in percent abundance of arbuscular mycorrhizal fungi (AMF), microbial biomass phosphorus (P) content, and phosphatase activity as well as available P after different cover crops, relative to the respective controls. On the left are displayed the number of observations. The lower-case letters indicate, for a single main crop type with a Tukey post-hoc test (p < 0.05), significant differences among cover crop types (including the control), and the upper-case letters between mycorrhizal cover crops, nonmycorrhizal cover crops and the controls. The corresponding models can be found in Table S3.5

Standard soil P testing (i.e. Olsen P or similar) was conducted in many studies after cover cropping in order to predict P availability to the main crop. Overall, cover crops had minor effects on the pools measured with these methods, with the exception of Lupinus sp., which increased Pa markedly. There was no evident relationship between Pa after cover cropping and main crop yield or P uptake.

Cover crop effects under different soil P conditions

We intended to explore the effect of cover cropping on soils differing in soil P status by classifying soil Pa into low and high. The majority of studies were conducted in soils with a low Pa. Cover crops had more pronounced effects on main crop performance in these soils compared with systems with abundant labile P (Fig. 7; Table S3.6: Models 6.1 and 6.2). Additionally, under conditions of low Pa, the cover crop benefit was greatest. Cover crop effects on soil microbial parameters were also influenced by soil P status, reflected in a much stronger increase in AMF abundance in soils low in Pa compared with high-Pa soils (Fig. 8; Table S3.7: Model 7).

Fig. 7
figure 7

Main crop yield and phosphorus (P) uptake as affected by cover cropping in soils with low and high available P (Pa). The points represent the percentage change of the modeled median (+/- 95% CI) of the cover crop treatments relative to the controls without cover crops. On the left are displayed the number of observations. The letters indicate significant differences among groups with a Tukey post-hoc test (p < 0.05). The corresponding models are presented in Table S3.6

Fig. 8
figure 8

Effect of cover crops on arbuscular mycorrhizal fungi (AMF) in soils with different P availability. The points represent the percentage change of the modeled median (+/- 95% CI) of the cover crop treatments relative to the controls without cover crops. On the left are displayed the number of observations. The letters indicate significant differences among groups with a Tukey post-hoc test (p < 0.05). The corresponding models are presented in Table S3.7

Discussion

We performed a meta-analysis to explore the importance of different plant P-acquisition strategies and to explain the benefit of cover crops on yield of main crops based on modified plant-microbe interactions during P cycling.

Main crops differ in their response to cover cropping

Our meta-analysis showed that cover crops have the potential to enhance both yield and P uptake of main crops in a variety of agroecosystems and under different management regimes, although the variance is very high. Main crops differ in their ability to profit from the P-benefit of cover crops, and this is related to their P-acquisition strategies. Maize and wetland rice yields increased more than soybean, cereals or vegetables/cotton/Brassicaceae. The response of maize to other cover crop benefits, especially Fabaceae, has also been seen in other meta-analyses (Alvarez et al. 2017), but the modelled high response in the present dataset may be explained by the inclusion of several studies conducted in low-input agroecosystems. The yield response of wetland rice is not directly comparable with that of other main crops due to the practice of flooding and the fact that all studies were conducted with Fabaceae cover crops. The limited benefit of cover crops to some main crops, i.e. canola, vegetables, and cotton, can be attributed in part to the reduced number of trials and possibly to greater management challenges (i.e. cover crop residues interfering with seedbed preparation) compared with other, more robust, arable crops.

Mechanisms underpinning the P benefit of cover-crop families

The main crop response was related to both the cover crop species used and the varying mechanisms of plant-microbial interactions. Fabaceae was overall the most effective cover crop family across all conditions and systems (Figs. 4 and 6). This family combines several of the mechanisms of P benefit: P uptake and carry-over in an abundant biomass with a high P concentration that facilitates release in synchrony with the main crop (Fig. 3), and a lasting effect on the soil microbial community (especially for mycorrhizal abundance and phosphatase activity) (Fig. 6). The N provided by symbiotic fixation provides an additional advantage through the acceleration of residue mineralization.

The separation of Lupinus sp. into a group distinct from other Fabaceae is justified: while Pmic and phosphatase activity were similar, the other (mycorrhizal) Fabaceae were clearly more favorable to mycorrhizal fungi than Lupinus sp. The most striking difference, however, was the effect on the pool of Pa (Fig. 6). The biochemical modification of the rhizosphere of Lupinus sp. increased the abundance of labile P under the main crop also, but its low biomass and low to intermediate P concentration probably limited its benefit to the main crops.

The absence of increases in main crop yield or at least Pa with pure stands of Polygonaceae, mainly Fagopyrum esculentum (buckwheat), was not expected in our meta-analysis, as buckwheat is used as a P-mobilizing species (Boglaienko et al. 2014). Whether or not buckwheat’s potential could be improved with other main crops, in a mixture with other cover crop species, or if its beneficial effects were limited by the low biomass observed in the studies included in the meta-analysis, warrants further investigation.

Asteraceae had intermediate positive effects on the main crops, but were tested exclusively as a cover crop for cereals. The species used in the experiments are all potentially mycorrhizal (Wang and Qiu 2006), but we lack detailed data on the soil biological parameters. The high biomass produced by these species is favorable as long as P concentrations are not too low.

Pure stands of Brassicaceae did not improve P nutrition of the main crop as much as most other families. This was possibly connected to a rather low interaction with the soil microbial community: neither mycorrhizal abundance, nor Pmic or phosphatase activity were increased significantly under the main crop. For species with high biomass production, such as Brassicaceae and Poaceae, P-cycling via the residue pathway is more important than for other cover crops.

Despite large amounts of P cycled through the biomass of Poaceae cover crops (Fig. 3), yield benefits for the main crops were limited (Fig. 4). Poaceae were most successful in increasing mycorrhizal abundance and microbial P, and they also enhanced phosphatase activity; negative effects on P-related soil biological parameters could therefore be ruled out. Poaceae produced the greatest quantities of biomass, but had the lowest mean P concentration of the cover crops, averaging 2 g P kg-1.

Decomposition of cover crop residues

In spite of their high biomass production and positive effects on soil microbial properties, pure stands of Poaceae were among the least successful cover crop families regarding P benefits for the main crops. This was possibly connected to P immobilization (Eichler-Löbermann et al. 2008), but may have involved other mechanisms as well, as, e.g., incomplete termination or problems with seedbed preparation.

Cover crop biomass and P concentration determine the amount of P cycled through the biomass, which can range between 1 and 30 kg P ha-1, although 3 to 10 kg P ha-1 is more typical, depending on cover crop species and P-availability (Fig. 3). The root:shoot partitioning of P in cover crops is variable, with 16 to 65% of the total plant P in the roots (Franchini et al. 2004). The threshold concentration of P in residues that determines immobilization/mineralization is 2-3 g P kg-1. The P contained in plant residues can be divided into available Pi and a recalcitrant Porg pool. The chemical composition of the plant parts changes with the developmental stage of the plant and with P availability. Cover crop residues are different from harvest residues, as cover crop plants do not reach maturity and Pi is the major pool in the cells (~70%) (Damon et al. 2014).

With tillage, 70-80% of the P in legume cover crop residues is released after 6 months, with roots being slightly more recalcitrant (Talgre et al. 2012). The processes responsible for the decomposition of cover crop residues are related to the mineralization of other pools of Porg in the soil. Cover crops increase phosphatase activity in the soil under a main crop. Although Poaceae, Fabaceae, and lupins tended to increase enzymatic activity more than Brassicaceae did, the plant type seems less important than the practice of cover cropping itself (Fig. 6).

In soils with sufficient P availability, microbial P mineralization is not driven by microbial P requirements, but rather by release of plant-available P as a by-product of carbon mineralization (Spohn and Kuzyakov 2013b). Current modeling approaches assume that microbial biomass releases P upon death, and is connected to the decay of the residue biomass through the availability of C substrates. A single pool for Pmic and residue Porg is used, assuming the same decay coefficients (Damon et al. 2014; Varela et al. 2017). Some authors argue that many studies on the decomposition of cover crop residues have used unrealistically large quantities of finely-ground residues and have not taken into account modification of the rhizosphere by the cover crops (Cavigelli and Thien 2003). Additionally, high concentrations of decomposing legume residues transiently increase the pH in soils with low pH and SOM, potentially leading to increased P availability (Vanzolini et al. 2017). Further field studies on residue decomposition dynamics with tight sampling frequencies are necessary. Different experimental set-ups can shed light on the relative importance of the transfer of P via the (shoot) cover crop biomass, e.g., removing the cover crop shoots, or by applying cover crop residues to previously unplanted soil (Rutunga et al. 2008; White and Weil 2010; Buyer et al. 2010).

The mechanisms by which the soil microbial community determines P dynamics during crop residue decomposition are not fully understood (Maltais-Landry and Frossard 2015), nor are the interactions between the microbial community and particular cover crops with their subsequent decomposition dynamics. The decomposer community adjusts to the cover crop species, as decomposable plant residues are produced over the entire growing period by senesced leaves or dead root hairs. Together with root exudates, this constant input of substrates constitutes a driver for shifts in microbial community structure, increasing the numbers of fast-growing copiotrophic microbes. Due to the strong impact of nutrient availability, phylogenetic diversity decreases in the vicinity of plant roots (Marilley and Aragno 1999); nevertheless, overall species richness in the field is expected to increase due to increased spatial heterogeneity. In a litterbag study in Brazil, mixtures containing Raphanus sativus with Secale cereale or Avena sativa showed peculiar dynamics, with a delayed increase of Pmic associated with RNA that could not be explained by the chemical characteristics of the residues alone (de Oliveira et al. 2017). The soil microbial community also influences suppression of weed germination by cover crop residues through selective decomposition of phytotoxic compounds (Moonen and Bàrberi 2006). The positive effects of cover crops on soil fauna (Blanchart et al. 2006) increase decomposition rates and nutrient cycling, probably persisting into the main crop phase. Although we focused on the soil microbial community, the importance of soil fauna (i.e. earthworms) for the shifts in P dynamics after cover crops should not be underestimated (Roarty et al. 2017).

The relationship between nutrient stoichiometry of the soil, microbial biomass, and plant residues determines microbial colonization and mineralization patterns. Fungi and other microorganisms capable of filamentous growth, unlike unicellular life forms, are capable of translocating nutrients between different compartments (i.e. soil-litter) to compensate for nutrient limitations. This has been demonstrated for C and N (Frey et al. 2003), and seems also plausible for P. Therefore, and bearing in mind the reduced damage to hyphae by reduced soil disturbance, fungal-driven decomposition probably dominates in no-till systems, where residues have less direct contact with the soil. Although no-till significantly increases fungal abundance, the interaction between the factors tillage and cover crop had no significant effect on the ratio of saprotrophic fungi:total bacteria in a long-term experiment under continuous cotton in Tennessee, USA (Mbuthia et al. 2015).

Nutrients other than P also need to be taken into account (Weerakoon et al. 1992), as residue mineralization dynamics and P release can be driven by N availability. An increased supply of N from a leguminous cover crop may permit the main crop to exploit its P-acquisition potential, resulting in increased P uptake. In field experiments, this effect is difficult to control, especially when a significant proportion of nutrients is contained in the cover crop root biomass. Cover crops can cycle substantial amounts of nutrients (potassium (K), magnesium (Mg), and calcium (Ca)) in their biomass (Wendling et al. 2016), and increase the availability of K (Cardoso et al. 2013). Descriptions of negative effects on plant nutrition are scarce and remain hypothetical. In some situations, the biomass pathway described in this paper could increase availability of potentially toxic elements, e.g., Mn, and lead to growth depression (Horst et al. 2001). Other researchers have considered the possibility that high levels of NaHCO3-Pi near the surface could induce copper (Cu) and zinc (Zn) deficiencies in conservation agriculture systems (Dube et al. 2014).

Arbuscular mycorrhizal fungi: early colonization assists crop P uptake by soil exploration

The strong increase, around 50%, in AMF abundance after mycorrhizal cover crops (Fig. 6), is important for the mechanisms of P benefit to the main crop. Most crops, with some notable exceptions, i.e. Brassicaceae, Polygonaceae, and Lupinus sp., can form symbioses with AMF with multiple benefits (Koide and Mosse 2004). AMF hyphae provide some of the functions of root hairs, especially in plant species with thick roots and very few or short root hairs, leading to exploration by plants of an increased volume of soil (Smith et al. 2011). Some functions and mechanisms of AMF symbioses are well known, although several fundamental issues remain unanswered. In AMF-colonized plants, the fungi are usually involved in P uptake with inhibition of a direct pathway via roots, and there may not always be positive growth responses (Smith et al. 2015; Ryan and Graham 2018).

There is no conclusive evidence that AMF-colonized plants are able to take P from soil sources that cannot be accessed by the roots themselves; rather they increase the soil volume from which the same P pools can be acquired (Smith et al. 2015). Direct release of phosphatases by AMF with a significant contribution to plant P uptake is under discussion (Joner et al. 2000). However, the substantial C input from plants through mycorrhizal hyphae extends our concept of a modified rhizosphere to a much greater soil volume, and the microflora of the mycorrhizosphere may play a critical role in P acquisition (Bending et al. 2006).

Mycorrhizal cover crops tended to increase main crop yield and P uptake more than non-mycorrhizal plant species did (Fig. 4), and AMF abundance is positively related to main crop yield and P uptake (S3.9 Models 9.1 and 9.2). With phosphatase activity, this direct relationship was not found. However, methodological deficiencies and the small number of studies may have resulted in a high variance. Due to sorption and stabilization onto soil particles, the enzymatic activity of a soil also reflects the recent history of a soil. This must be taken into account when interpreting the results of the meta-analysis. Samples were taken after termination of the cover crop or under the main crop. Therefore, the enhanced activity relative to that of the fallow control treatments corresponded either to residual phosphatases released by the cover crop roots, to changes in abundance or structure of the soil microbial community, or to an overall substrate-driven increase in phosphatase activity due to P-rich cover crop residues in the soil. Increases in phosphatase activity mirrored the effect of cover crops on microbial biomass P, indicating a potential microbial origin of the enzymes.

A mycorrhizal cover crop can transfer its ability to access P in the soil to the main crop in the form of mycorrhizal inoculum (hyphae or spores in the soil). Although the AMF-plant symbiosis is unspecific and larger plants may be simultaneously colonized by different mycorrhizal fungal species, there are some plant-AMF genotype combinations that are more efficient than others (Jansa et al. 2011). Molecular techniques make it possible to describe AMF diversity associated with specific cover crops (Sharrock et al. 2004), and increases in richness in the main crop have been reported (Ramos-Zapata et al. 2012). However, current knowledge gaps regarding the connection between AMF community assemblage and species function constrains the effective translation of this information into specific cover crop species recommendations.

The build-up of AMF inoculation potential benefits only AMF-competent main crops, and the ability of the main crop to take advantage of earlier increased mycorrhization by previous cover crops determines the P benefit (Bittman et al. 2006). For maize, a positive relationship between mycorrhizal colonization and plant biomass or P content was found in both our aggregated dataset and several single studies (White and Weil 2010; Njeru et al. 2014). Cavigelli and Thien (2003) reported Lupinus albus unexpectedly decreased sorghum P uptake in a pot experiment, although P uptake and biomass was the highest of the tested winter cover crop species; a possible explanation was that lupin was the only non-mycorrhizal crop in the study. However, the lower potential for AMF inoculation cannot have been the only reason, as non-mycorrhizal Brassicaceae performed better, whereas mycorrhizal Poaceae combined poorly.

It is important to bear in mind that non-mycorrhizal crops have evolved special strategies for P acquisition. Non-mycorrhizal families can be broadly classified into Brassicaceae and Proteaceae groups, which evolved in P-rich and severely P-impoverished environments, respectively (Lambers and Teste 2013). Raphanus sativus var. oleiferus (oilseed radish) exudes large amounts of acid phosphatase and other rhizodeposits into the rhizosphere (Kunze et al. 2011); the exceptional P-mining strategies of Lupinus sp. have been described above. The impact of these P-acquisition strategies must be considered when designing site-specific crop rotations that include cover crops.

The meta-analysis also showed an unexpected slight tendency toward increased mycorrhizal abundances, after non-mycorrhizal cover crops. However, species of plant families labeled as “non-mycorrhizal” can be infected at low levels by AMF (Lambers and Teste 2013). Another possibility is that higher herbicide usage in the control treatments decreased AMF abundance because of direct toxic effects (Trappe et al. 1984; Giovannetti et al. 2006) or by fewer weeds acting as mycorrhizal hosts (Oehl et al. 2003).

The use of cover crops to build up the inoculum potential of beneficial microorganisms, including AMF, has the capacity to considerably improve soil fertility (Galvez et al. 1995; Boswell et al. 1998; Bagayoko et al. 2000; Kabir and Koide 2002; Lehman et al. 2012), although apparently not in all agroecosystems (Sorensen et al. 2005; Higo et al. 2014). In a study in USA, 31 years of Vicia villosa cover cropping decreased mycorrhizal abundance relative to wheat or no cover crop, with some responses associated with high N rates (Mbuthia et al. 2015). Also, in fields with a history of mycorrhiza-enhancing cropping techniques (i.e. rotation dominated by mycorrhizal crops, no-till), cover crops may not increase the mycorrhization of the main crop further (Turmel et al. 2011). In agricultural soils very low in mycorrhizal abundance, cover crops may fail to increase the inoculum potential above a minimum threshold necessary to benefit plant growth (Douds et al. 2011). On the other hand, cover crop mixtures and a strategic AMF-build-up may be especially important in this context (Lehman et al. 2012).

Weeds growing during the off-season may also result in benefits for the main crop; e.g., Taraxacum officinale (dandelion) is a good host for overwintering mycorrhizal fungi (Kabir and Koide 2000). However, some weeds induce negative changes in the microbial community (i.e. a decrease in AMF) and enhance their own competitive advantage over the crops (Wortman et al. 2013). Costs of seeds and labor for cover crop establishment have to be included in evaluating their potential to outperform weeds, and, above all, their easy termination (Wang et al. 2008). Also, management and application form of cover crops determines their effect on AMF, as fresh red clover residues, directly incorporated or used as a mulch layer, result in greater abundance of AMF compared with processed residues (biogas slurry, compost) (Elfstrand et al. 2007). The discovery of significant advantages to seedlings conferred by early establishment of symbioses with mycorrhizal fungi or P-solubilizing bacteria has prompted the development of commercial and on-farm produced inocula for plant growth promotion (Douds et al. 2010). Also, combining cover crops with a simultaneous inoculation of microorganisms has been investigated (Cui et al. 2015).

It is experimentally challenging to separate the direct benefits of improved P availability from other plant-microbial interactions and indirect plant growth-promoting effects. Microbes have multiple effects on plant health via plant pathogen suppression as well as on nutrient cycling and plant nutrition (Bagayoko et al. 2000; Bashan et al. 2013). Cover crop species differ substantially in their root-associated fungal communities (Benitez et al. 2016), suggesting opportunities for management of beneficial and pathogenic fungi, although contradictory evidence has been reported (Turrini et al. 2016). Nutrient uptake and crop yield are often limited by biotic stress. Increased aboveground biodiversity (i.e. by cover cropping) can suppress soil pathogens, thus promoting plant growth. Some cover crop species can be used for specific pest management, whereas others may serve as hosts for pests (Ratnadass et al. 2012). In view of the multiple ecosystem functions of soil microbes, Fester and Sawers (2011) advocated a holistic approach of increasing biodiversity through agronomic management as opposed to a reductionist approach focusing on single species.

Microbial P as a significant pool

Increases in microbial P of around 25% with Poaceae and Fabaceae cover crops (Fig. 6) are worth further discussion, because Pmic constitutes an important pool in soil due to its relatively fast turnover and subsequent availability for plants. However, due to the small number of studies, we did not detect a significant effect of Pmic on main crop yield and P uptake. Interestingly, the mycorrhizal cover crop types tended to enhance Pmic more than the non-mycorrhizal species. The pool of Pmic in agricultural soils typically constitutes 5–70 kg P ha-1, with turnover times of a few months, depending on management and C inputs (Oehl et al. 2001).

Phosphorus in microbial cells is present mostly in the form of nucleic acids, but also as small P-containing esters, free Pi, and phospholipids of cell membranes. The nucleotide content can vary greatly, depending on the growth rate of the cell; surplus P can be stored as polyphosphates (Harold 1966; Stewart and Tiessen 1987). When comparing the nutrient stoichiometry of soils, microbial biomass, Pmic appears to be closely linked to overall microbial biomass (Cleveland and Liptzin 2007). However, in agricultural soils, the microbial C:P stoichiometry may, in response to soil fertility and management, exhibit some plasticity and be affected by the availability of P for the main crop.

Cover-cropping frequency was a greater driver of increases in microbial biomass than compost application, increasing the abundances of Pseudomonas and Agromyces species, including species that are important biological control agents and plant growth-promoting rhizobacteria (Brennan and Acosta-Martinez 2017). In a study in semi-arid Kenya, microbial C, N, and P were strongly increased by different Brachiaria species (Gichangi et al. 2016). In another study, cover cropping was far more important than the measured environmental variables (moisture, temperature, pH) in controlling soil microbial community structure (Buyer et al. 2010). Relative abundance of Gram-positive bacteria was decreased by cover cropping, probably connected to their lower ability to use labile C-inputs compared with other microbial groups. In a field experiment in Sweden, direct incorporation of a red clover crop enhanced and sustained microbial biomass and soil enzyme activities more than did processed forms of green manure applied as biogas slurry or compost (Elfstrand et al. 2007).

The application of isotopic dilution methods has revealed that microbial immobilization and remineralization, rather than mineralization of non-living organic P, represents most of the gross organic P mineralization flux (Bünemann 2015). Sorption/desorption processes dominate in agricultural soils low in microbial biomass. Immobilization/mineralization dynamics are closely linked to overall microbial growth, but trophic interactions should not be overlooked, as amoebae and other bacteria-grazing microfauna are responsible for the remineralization of Pmic (Cole et al. 1977). Cover crops increase the abundance and diversity of microfauna (Blanchart et al. 2006) which restricts long-term microbial immobilization during the decomposition of plant litter. Given the close relationship between Phaseolus vulgaris yield and Pmic, microbial biomass P was proposed as an indicator of soil P availability on P-sorbing Andosols in Japan, instead of the widely used Truog-P (Sugito et al. 2010).

Mobilized and mineralized P is often intercepted by microorganisms before plant roots can take it up (Joner et al. 2000). Although in managed agroecosystems crop yield is the primary objective, some interactions between crops and other organisms, seen as competitive and undesirable in the short term, may have, in principle, longer-term favorable properties (i.e. SOM build-up, pathogen resistance). Competition for available nutrients determines the outcome of the plant-microbe relationship at all levels, from the more opportunistic microbes in the rhizosphere to the mutualists, as even AMF-symbioses can range from being highly beneficial for both partners to ‘parasitic’ (Johnson et al. 1997). In addition to mycorrhizal fungi, there are many other plant-associated organisms that affect P uptake. One example of the very complex interactions between microorganisms and plants are mycorrhiza-helper bacteria, which facilitate the establishment of the symbiosis (Frey-Klett et al. 2007).

Cover crops under different soil-management strategies

The effects of cover crops vary greatly, and some studies report no or even negative effects (Kuo et al. 2005; Takeda et al. 2009b; Rick et al. 2011). This has been due, in some cases, to intrinsic agronomic conditions related to low cover crop biomass or absence of a significant P limitation on experimental plots. Rick et al. (2011) did not find effects of cover crops on labile soil P fractions, wheat biomass, or P concentration, despite differing biomass and P concentrations of the cover crops. However, as neither cover crops nor rock phosphate fertilization resulted in substantial yield responses, it is possible that the lack of positive results was due to a combination of moderately high soil P levels, N limitation, and low precipitation. Also, studies in Brazil on strongly P-sorbing acidic soils, where cover crops often increase labile and moderately labile P pools while decreasing residual P, did not always show the expected increases in terms of P uptake (Almeida and Rosolem 2016). After 3 years in a Brazilian Hapludox under no-till, Avena strigosa (black oat), Vicia sativa (common vetch), and Raphanus sativus (fodder radish) as cover crops cycled P from the non-labile and moderately labile P pools through their biomass without reducing labile P fractions; maize yields were also not affected (Pavinato et al. 2017). Increased P availability was also not consistently translated into improved main crop performance (Pavinato et al. 2017). In other studies, however, although labile P fractions were not affected, yields were increased (Murungu et al. 2011b; Dube et al. 2014), or traditional soil P tests failed to detect the shift in P dynamics by cover crops (Takeda et al. 2009a, 2009b). Reasons for this variability in response are manifold, and reflect both methodological limitations when assessing P dynamics and the diversity of the studies’ designs, as well as abiotic and biotic factors. This points to practical problems in agricultural management that often hamper the successful exploitation of cover crop benefits. Recurrent calls for genetically-engineered crops for improved P efficiency, e.g., in Hunter et al. (2014) overestimate our understanding of the complex rhizosphere processes involved, while existing agricultural management options are underrated.

The importance of site conditions and agricultural management in controlling both growth of cover crops and main crops and the complex mechanisms of their interactions with P explains some of the variation in results of the studies included in the meta-analysis. However, knowledge of the impact of soil P status and different management approaches will not only aid in the interpretation of field experiments, it will also provide many tools for adapting cover crop effects to the specific needs of local agroecosystems. Cover-cropping management can be used according to site conditions in many ways, e.g., by appropriate combinations of cover crop species/mixtures with the crop rotation (e.g., de Oliveira et al. 2017), and by fertilization and tillage (e.g., Mbuthia et al. 2015; Teles et al. 2017). Other management decisions, such as seeding rate (Brennan et al. 2009), seeding and termination date (Nascente et al. 2013), and termination technique (Dorn et al. 2013), extend the opportunities for fine-tuning cover crop performance and P effect. However, as cover crops provide multi-functional tools when enhancing performance for a given function, there can be trade-offs with other functions (Mirsky et al. 2012). For example, on soils with low P availability, a highly-efficient scavenging cover crop with an extensive rooting system that takes up P from the same pools as the main crop could lead to soil P immobilization when the easily-available P is depleted and immobilized in recalcitrant plant residues. Management options could include selection of a different cover crop species/mixture or a later sowing/earlier termination date to reduce the recalcitrance of cover crop residues.

Soil P and fertilization

Phosphorus dynamics, and the effect of cover crops on them, are strongly influenced by the soil, specifically by the size of easily- and sparingly-available P pools. The meta-analysis supports the hypothesis that cover crop P benefits are more evident in soils poor in Pa than on sites with higher Pa (Fig. 7). Also, interactions between cover crops and the soil microbial community are influenced by soil P status: the increase in AMF abundance after mycorrhizal cover crop use is much greater in soils low in Pa (Fig. 8); phosphatase activity shows a similar pattern (data not shown). These results are promising, and the P-mobilizing properties of different cover crop species should be further investigated, specifically in soils high in Pt, but low in Pa. However, soils very low in Pa often also have other fertility problems. Soil quality affects the beneficial effects of cover crops in several ways. Although very fertile or well-managed soils are difficult to improve further with cover crops (Turmel et al. 2011), cover crops require adequate initial soil fertility to effectively improve P cycling and inputs may be necessary to obtain a functional cover crop (Jensen et al. 2005). Examples of useful soil quality parameters are available P, microbial biomass C (Cmic), and the Cmic:C ratio (Koné et al. 2008). Depending on local economic circumstances, a good use of expensive inorganic P fertilizer might be to apply P sources to cover crop legumes; this would improve their ability to benefit from existing soil N as well as residually applied P (Carsky et al. 2001). There are substantial synergies between P application and N fixation by leguminous cover crops (Weerakoon et al. 1992).

Several studies that did not find a consistent effect of cover crops on soil P pools led some authors to argue that the P-mobilizing properties of legumes are more efficient for low-input agriculture on soils with low P availability than for systems where P availability is higher and P limitation weaker (Maltais-Landry et al. 2015). We agree with this line of reasoning and suggest consideration of P reserves to complement the labile P assessed in standard soil P testing when selecting and evaluating cover crops. The classification into high- and low-P soils is a simplification, which becomes clear when treating soils high in Pt but very low in Pa, as in southern China or Andosols in Chile and Japan. On these soils, P-mobilizing cover crops would be of particular interest. Because cover crop P benefit reflects mainly biological processes, these may not be assessed properly with methods designed to predict the effects of inorganic P fertilization. In high-input systems with high P availability, tillage, and additional P fertilization, cover crops often do not increase yields significantly over the control (data not shown).

A positive interaction of cover crops with P fertilization is in part related to the fertilizer type used. Cover crops may increase the availability of the P contained in organic fertilizers or phosphate rock (Ca-P/apatite). In particular, proton- and organic anion-exuding plants such as Fagopyrum esculentum or legumes, but also green manure from Tithonia diversifolia, dissolve calcium (Ca)-P, which is important for organic farmers on calcareous soils (Arcand et al. 2010) and farmers relying on phosphate rock on acid soils (Somado et al. 2003; Ikerra et al. 2006; Oikeh et al. 2008; Opala et al. 2010). However, this effect has not always been found (Rick et al. 2011). Plant-microbial interactions likely play a role in fertilizer effects (Bah et al. 2006). Indications of an interaction between P-fertilizer type (soluble P vs rock phosphate) and cover cropping has also been shown in some studies (Almeida and Rosolem 2016), but was not considered in the present meta-analysis. Plant species respond differently to fertilizer types. Brassica oleracea showed a dramatic response to an otherwise unreactive, Fe-rich, igneous phosphate rock, but this source was ineffective for leguminous cover crops and maize; the following maize yield was not increased by cabbage (Weil 2000). Aluminium phosphate is naturally present in some rock phosphates and soils, but may be inadequate for plants that rely on proton exudation for P-acquisition, as pointed out by Pearse et al. (2006), comparing wheat and lupin in a pot experiment. The interactions between fertilizer type, cropping system and soil are complex (Romanyà and Rovira 2009), and require, therefore, site-specific recommendations.

Some studies have shown a more efficient use of P when added with organic amendments (Eichler-Löbermann et al. 2008; Maltais-Landry et al. 2015), but not others (Takeda et al. 2009a). Under the condition of added C sources, effects of direct competition for Pa by microbial immobilization are at variance with a rate increase in cycling of organic P. The properties of the amendment must be taken into account when selecting the cover crop. Fertilizers with high levels of Pa should be combined with species such as Raphanus sativus oleiformis (oil radish) that maximize uptake and biomass production; species with mobilization traits, such as phacelia would be suitable for sites with less-available P sources (Bachmann and Eichler-Löbermann 2010).

Continuous developments in alternative and innovative cover-crop systems have the potential to increase sustainability of intensive systems such as high-input horticulture (Brennan 2017). In systems with substantial organic inputs (i.e. manure), P availability is usually not a concern, and other cover crop effects are desired (prevention of nutrient leaching, weed suppression, C input). An overlooked property of cover crops regarding P dynamics is connected to the inherent imbalance of the N:P ratio in animal manure: legume cover cropping permits lower rates of manure application by supplying N to subsequent crops (i.e. supporting P-based application rates), reducing the P excess of N-based organic fertilization (Kleinman et al. 2001; Cherr et al. 2006). A special, but not uncommon, situation is the management of large quantities of manure from industrial animal farms by field application in combination with cover crops (Rowe et al. 2006).

Cover crops can be used in watershed management for the reduction of P-runoff (Villamil et al. 2006; Geleta et al. 2006). Legumes are generally appreciated for their ability to mobilize poorly-available P, whereas grasses are more often used as “catch crops”, scavenging the available nutrients and reducing losses (Maltais-Landry et al. 2015). The authors concluded that, in soils with low P-sorbing capacity, P transfer via the cover crop biomass of grasses is more effective than that of legumes, which more strongly modify their rhizosphere. Systems at risk of P losses usually also have problems with N losses (Sharpley and Smith 1991; Aronsson et al. 2016; Brennan 2017), making a focus on grasses instead of legumes meaningful.

Cover crop mixtures

Cover crop mixtures are, in terms of main crop performance, superior to monoculture cover crop species (Figs. 4, 5 and 6); however, most studies used maize as a main crop, which is highly responsive to cover cropping, and also the total number of trials was low. It is difficult, therefore, to draw general conclusions about cover crop mixtures due to the inherent differences in systems depending on the specific components. Mixtures frequently outperform single species in terms of biomass production and P uptake (Li et al. 2007, 2014; Messiga et al. 2016), in addition to the positive influence plant biodiversity exerts on soil biology. Other plants growing in close association with P-mobilizing plants confer additional benefits through intercropping or undersowing, and may also increase access to sparingly-soluble P (Li et al. 2007). Cereal-legume mixtures are among the most widely used and studied (Tarui et al. 2013), due to both their ecological importance and to the availability of practical management expertise with these combinations, because of their use in fodder production. An important benefit of this association is increased N-fixation by the legume driven by the N-demand of the cereal (Høgh-Jensen and Schjørring 1997), N transfer (Brophy et al. 1987), and facilitative interactions via root exudates, as detected between Zea maize and Vicia faba (Li et al. 2016). In cereal-legume mixtures, the cereal component benefited more from the intercropping association than the legumes, reflected in the observed shift in composition towards a higher proportion of grasses over time (Maltais-Landry 2015). However, the legume productivity may also be enhanced by intercropping, because Zea maize improved Fe nutrition of Arachis hypogaea (peanut) in a calcareous soil (Zuo et al. 2000; Zuo and Zhang 2008). In mixtures, there are several trade-offs that have to be balanced; e.g., easily decomposable cover crop residues improve nutrient availability for the main crop, but result in lower weed suppression and less SOM production in comparison with more recalcitrant plant residues (Tarui et al. 2013).

Desirable mixtures of species with complementary functions have a huge potential to increase the benefit of cover cropping, as, e.g., Brassica napus (fodder rape) and Lupinus albus (Little et al. 2004). There is also the advantage of the biofumigant properties of Brassica species and its biomass, as well as the P mobilization of Lupinus sp. (although a mycorrhizal component could enrich the combination). The combined use of a mixture of rye and oat as winter-hardy and winter-killed species also yielded positive results in temperate climates (Kabir and Koide 2002). The transfer of nutrients from dying roots to living roots via AMF, as suggested by Newman and Eason (1989), should also be taken into consideration.

In the meta-analysis, AMF abundance following mixtures containing non-mycorrhizal Brassicaceae tended to be only slightly less than that with pure Poaceae or Fabaceae stands. This suggests that one mycorrhizal partner in a mixture is capable of compensating for the lack of mycorrhizal association of the other. An experiment in mid-Atlantic USA showed that a pure stand of Raphanus sativus var. longipinnatus (forage radish) did not have a negative effect on AMF-colonization, but a pure stand of rye showed a higher rate of colonization than the mixture rye-forage radish. Further research is necessary to determine whether or not Brassicaceae inhibit the extent to which Poaceae are colonized by AMF (White and Weil 2010). As for potentially adverse effects resulting from nonmycorrhizal and mycorrhizal plant species, there are indications that members of the Proteaceae type exhibit fewer competitive interactions than Brassicaceae types (Gardner and Boundy 1983; Lambers and Teste 2013; Lambers et al. 2018).

The plethora of possible combinations and proportions in cover crop mixtures is difficult to handle experimentally. Even studies with more than four species are scarce, although a systematic and simultaneous screening would be important for the selection of appropriate mixture candidates (Horst et al. 2001; Oikeh et al. 2008; Wendling et al. 2016). The selection and testing of innovative cover crop species requires more large-scale projects, such as OSCAR (Crossland et al. 2015), or the 2017 started EU-project REMIX. There is substantial variation, even among genotypes, regarding P-uptake efficiency and subsequent main crop growth (Jemo et al. 2006; Rose et al. 2010a; Pang et al. 2018b). Likewise, site conditions are important, because in an acidic soil, root growth determined P uptake, whereas in an alkaline soil malate exudation was probably the more important parameter (Rose et al. 2010a). With respect to soil biological properties, a truly novel and strategic experimental approach is necessary to assess the multiple combinations of cover crop mixtures.

The systematic use of plant traits that interact with P dynamics, as proposed by Wendling et al. (2016), is already used for other cover crop functions (Damour et al. 2014), and constitutes an important step toward increasing the comparability of trials.

Interactions of cover crops and tillage

Cover cropping tends to have greater positive effects on main crop performance in systems under reduced tillage/no-till than under conventional tillage (Fig. S6). Systems using reduced tillage benefit from the inclusion of cover crops due to weed control, nutrient release, and improved soil structure, together with synergies between both practices regarding soil biological activity. Tillage regime can vary in intensity (shallow/ non-inversion) and frequency (up to no-till), affecting, mainly in two ways, the processes by which cover crops influence soil P dynamics. First is litter distribution; tillage mixes plant residues into the soil, whereas no-till leaves the nutrient-rich residues on the top. Second, tillage regime changes the soil biota, because organisms differ in their sensitivity to soil perturbation.

Cover crops and no-till benefit soil fauna, which are important for residue decomposition, because their feeding activity fragments and buries litter, increasing the surface for microbial decay. Earthworms increase P availability and their interaction with cover crops deserves further attention (Vos et al. 2014). In tropical and arid areas, termites fulfill important functions in nutrient cycling (Rückamp 2011). On the other hand, under conventional tillage, the greater area of plant residue contact with soil increases the decomposition rate, but can lead to P sorption (Tiecher et al. 2012b). The choice of cover crop can have additional effects, e.g., with the release of isothiocyanates by tilled residues of Brassica juncea (Indian mustard) used as biofumigants, possibly increasing negative effects on mycorrhizal-inoculum potential (Njeru et al. 2014). Under no-till, a forage radish cover crop did not negatively affect mycorrhization of maize (White and Weil 2010).

Balota et al. (2014) found that several microbial parameters increased under cover crops, both under no-tillage and under conventional tillage in a Brazilian oxisol. There was a greater relative increase with cover cropping under conventional tillage, although total microbial abundance and activity remained higher under no-tillage when the entire profile was taken into account. The higher metabolic efficiency (lower qCO2–values) indicated that winter cover crops and zero tillage resulted in a more stable system. In the USA, no-till plots increased levels of Gram-positive bacteria, actinomycetes, AMF, and enzymatic activity, whereas tillage enhanced the abundance of saprotrophic fungi and provided a greater total microbial biomass. The cotton yield was greatest after Vicia villosa under no-till (Mbuthia et al. 2015).

One of the principal reasons for the inclusion of cover crops in temperate agroecosystems is the alleviation of soil compaction, which in turn expands the effective rooting zone, benefiting main crops possessing weaker root systems (Calonego and Rosolem 2010). In this context, the property of Raphanus sativus, the so called “tillage radish”, deserves further attention, as it produces a strong taproot that decays during the winter and leaves distinct biopores in the surface soil, potentially leading to locally greater P-availability (White and Weil 2011).

Cover crop residues at the soil surface can delay soil warming in spring, and, therefore, the onset of microbial activity; this constrains the use of mulches in cold climates (Sarrantonio and Gallandt 2003). The increased soil water content under mulch layers can also constitute a problem in wet years (White and Weil 2010). Some studies suggest that cover crops have only limited potential to decrease P losses under the wet and cold climatic conditions of northern Europe, as dissolved organic P leaching over winter may outweigh reductions in erosion losses when plant residues experience several freeze-thaw cycles in the field (Bergström et al. 2015; Aronsson et al. 2016; Kirchmann and Wessling 2017). However, cover crops that were incorporated into the soil during winter did not increase P leaching in a study in Belgium (Vanden Nest et al. 2014); similarly, at least in milder Nordic climates, winter-hardy cover crops could be an option.

The termination method

When cover crops do not die off via natural causes (e.g., frost or drought), a termination step is required. The most widespread method is a termination herbicide, usually glyphosate. However, this practice has been criticized because of undesirable side effects on the agroecosystem and the environment (Johal and Huber 2009; Yamada et al. 2009; Mamy et al. 2016). Glyphosate may also reduce main crop yields when applied too close to the main crop seeding date (Nascente et al. 2013) or it may interact with P-fertilizer application (Rose et al. 2018). Termination can also be mechanical, by tillage of varying degrees of intensity, mixing the plant residues with the soil, or by flailing, disking or rolling the shoot biomass, resulting in a mulch layer that covers the soil surface. The use of roll-choppers requires exact timing to be effective, but can be an adequate alternative to glyphosate (Creamer and Dabney 2002; Dorn et al. 2013). The termination step interacts with the temporal dynamics of the residue mineralization, important for synchronization of nutrient release with the requirements of the main crop (Zibilske and Makus 2009; Murungu et al. 2011a; Damon et al. 2014). The achievement of desirable C:N:P ratios with timely termination may be constrained by field accessibility related to soil water content (Odhiambo and Bomke 2007). The cover crop species should be adapted to variations in sowing and termination dates in order to achieve high biomass and nutrient content, which does not necessarily follow a linear pattern (Anugroho et al. 2009). Cover crops are usually terminated between one and a few weeks before planting the main crop in cases where they are not winter-killed. With chemical termination, a certain time for pesticide inactivation needs to be taken into account.

In some rotations, a cover crop (e.g., white clover before maize) may simply be clipped and permitted to regrow as living mulch, conferring high AMF inoculation potential and improving P nutrition (Deguchi et al. 2007). An interesting approach for integrated crop-livestock systems consists of cover crop termination by grazing (Clark 2008), which can increase P efficiency under adequate management (Costa et al. 2014). Nutrient transfer by cut and carry use of green manure cover crops on improved fallow fields may lead to nutrient impoverishment in low-input systems. Shrub species with high litter and root production, such as Tithonia diversifolia, may be advantageous in these situations (Rutunga et al. 2008). Burning of bulky cover crop biomass, as practiced by some resource-limited smallholder farmers, reduces fertility benefits and increases nutrient losses (Oikeh et al. 2008).

Cover cropping: long-term effects and adoption

Contrary to our expectations, the published studies found no strong evidence that cover crop effects on P nutrition of the main crop increase with time after introduction into the cropping system. The few long-term studies available (Kuo et al. 2005; Abdollahi and Munkholm 2014; Mbuthia et al. 2015; Mukumbareza et al. 2015) did not report substantially better results regarding P benefit than short-term experiments over 1 or 2 years. This indicates that cover crop management and selection of appropriate species could be more important than time since adoption of cover cropping, although more long-term trials are warranted. However, all the aforementioned studies and more (Balota et al. 2014) reported increases in microbial biomass and enzymatic activity. Regarding P pools, cover crops increased mainly the Porg and Pmic pools over time (Maltais-Landry et al. 2015; Mukumbareza et al. 2015) which may have been related to overall increases in SOM (Blanco-Canqui et al. 2015). Another study found no changes in soil P fractions after 9 years of cover cropping in a temperate system with high P availability (Kuo et al. 2005).

Cover cropping tends to provide better results in tropical than in temperate climates (data not shown). However, the numerous studies that show important main crop P benefits by cover cropping in temperate systems suggest that perhaps not only climate makes a difference, but also the prevailing agroecosystem. Studies in the tropics often involve low-P soils receiving fertilizer treatments, and are therefore more responsive to cover cropping compared to the northern countries whose mouldboard-ploughed fields are high in Pa. Additionally, given the management challenges of cover cropping, the greater experience with cover crops in the tropical regions may be an additional factor. However, except for regions with the limiting factors of water scarcity (as, e.g., in Mediterranean climates or with further climate change) and the short vegetation periods of cold areas, which strongly influence both plant growth and decomposition dynamics, cover cropping has the potential to be a successful management option in most climates.

Cover cropping is one of many agricultural tools and must be integrated and adapted adequately into the management strategy of an agroecosystem. Evaluation of both the multiple benefits and site-specific management needs of cover crops requires multi-faceted and whole-system approaches in research and extension (Cherr et al. 2006). Cover crop species relevant for the local farming system for other reasons than P nutrition should also be tested. Economic benefits may be indirect, e.g., from the possibility of reduced planting densities after cover crops (Wang et al. 2008). An interdisciplinary or multidisciplinary approach involving farmers and other practitioners at early stages of experimental design may increase the efficiency and practical relevance of scientific studies (Weil and Kremen 2007; Reed 2008; Scopel et al. 2013; Smith et al. 2014).

Despite its many benefits, the use of cover cropping by farmers is often less utilized than is desirable. Barriers to adoption include the following: benefits are site- and soil-specific, establishment and management problems exist, and climatic variability can lead to uncertainty in outcome (Ro et al. 2016). Water use by the cover crop may be a problem under drought stress, resulting in lower yields and P uptake (Turmel et al. 2011). Economic considerations are necessary for adoption of a farming practice: in warmer climates there may be no off-season which means that cover crops compete directly with cash crops for the same space. To compensate for the income loss, the benefit over a continuous cropping system needs to be very high. Easily-manageable multi-purpose cover crops (edible seeds, fodder, wood) would be required to increase adoption by smallholder farmers in developing countries. Also, external factors, such as changes in land markets or novel pests, determine the adoption or abandonment of sustainable practices (Neill and Lee 2001). Reasons for low adoption of cover crops in temperate high-input systems are discussed by Brennan (2017).

The yield and P effects of cover crops hold across both organic and conventional cropping systems, as the results from five studies with 45 observations using organic fertilizers and no chemical weed control were not significantly different from the nine studies with 103 observations using conventional management practices (data not shown). In order to achieve widespread adoption of cover cropping, it is essential to overcome constraints to its adoption (e.g., insufficient experience for specific site conditions, unavailability of machinery, missing management alternatives in case of cover crop failure, and general pressure by markets to externalize costs). Improved management, increased knowledge, and development of instruments to relieve farmers of socially and environmentally unsustainable short-term market pressures, as for example community supported agriculture (CSA), are important (Lass et al. 2003). With rising fertilizer prices and pressure to reduce environmental impacts, cover crops constitute a promising, multifunctional tool for sustainable intensification of agriculture, on the conditions that species selection and management match the agricultural goals (Blanco-Canqui et al. 2015).

Conclusions

Cover crops can be successfully used to stimulate main crop yield and P uptake. However, site conditions and agronomic management lead to varying results. Plants have different P-acquisition strategies, and main crops show varying abilities to take advantage of the cover crop legacy. Cover crops benefit the P nutrition of main crops by different, simultaneous processes: soil P (sometimes from sparingly-available pools) accumulates in the cover crop biomass, and the mineralization of P-rich litter provides available P for the main crop (plant-storage pathway). This pathway is most relevant for cover crops with high biomass such as Poaceae, Brassicaceae, and Fabaceae (Fig. 9). The P concentration of the cover crop biomass determines mineralization dynamics, which may partly explain the limited efficacy of Poaceae cover crops. Cover cropping enhances soil microbial community abundance (Pmic) and activity (extracellular phosphatase activity), and maize, in particular, benefits from increases in AMF abundance in soils with low available P (soil microbe pathway). Poaceae and Fabaceae have the greatest impact on soil microorganisms. Other cover crop species (e.g., Lupinus sp.) are capable of mining P pools, improving soil P availability even during the main crop phase (biochemical rhizosphere modification pathway).

Fig. 9
figure 9

Radar chart summarizing the properties of the cover crop families and their effects on soil. The lines correspond to the calculated quantile moment of each data point (R code provided in S4). Grid lines correspond to the 0, 25, 50, 75, and 100-quantiles of each variable. Asteraceae and Polygonaceae had missing data points and could not be displayed here.

Consideration of the abovementioned P-acquisition mechanisms, the interactions with the crop rotation, and the use of plant traits for the characterization of cover crop species would facilitate the generalization of the results of different studies. Further research is needed to elucidate the relative contributions of the different P-acquisition mechanisms, both to P uptake of the cover crop and their effects on the main crop, in order to optimize combinations. Soil P availability affects the P dynamics of the system and the mechanisms of cover crop benefit: generally, cover crops enhance main crop yield and P uptake in systems low in Pa more than in systems with abundant Pa, and have a greater effect on abundance of AMF.

Cover crop benefits are greater with reduced tillage or no-till. Management determines the success of cover cropping in general, and it is possible to fine-tune P dynamics with appropriate techniques. Cover crops are used on a global scale under varying circumstances with successful examples in many regions and agroecosystems. However, cover crops may not be the optimal solution for P management in all situations. Factors that limit the use of cover crops, i.e. adverse climatic conditions (water scarcity, short growing seasons, freeze-thaw cycles), insufficient soil fertility, or problems with pests can be overcome by cover crop selection and management. However, for conventional high-input systems, though there is little likelihood of yield improvement, there is potential for environmental benefits.

Finally, in many areas of the world, the principal reason to include cover crops into the rotation is erosion control, which is probably the most important global issue regarding P management and soil in general. We need to avoid perpetuating systems in which single characteristics are overemphasized, resulting in significant trade-offs in overall performance and sustainability of agroecosystems. The isolated effect of a management practice (cover crop) on a single nutrient (P), may make sense in situations with one dominating limiting nutrient. Yet a more comprehensive evaluation of the ecosystem is required in most situations (Schipanski et al. 2014). This is especially the case for complex systems that replace technological inputs with ecosystem services of biological components. We need to find a balance that takes advantage of the numerous contributions of cover crops to agroecosystem health.

The effects of cover crops on P uptake of the main crop depend on many factors, offering opportunities for site-specific adoption and optimization of the system, but also restraining general agronomic recommendations. However, we can draw some broad conclusions about the potential for P management by cover-crop management and directions for future investigation:

  • The different mechanisms of cover crop P benefit we have discussed, i.e. P transfer via cover crop residues, organic anion exudates, root-exuded enzymes, and microbial interactions, may happen simultaneously and warrant further investigation.

  • Cover crop biomass determines in many cases the magnitude of its effect, because, in addition to the transfer of P in plant residues to the next crop, it affects the potential for rhizosphere modifications and microbial interactions. Appropriate management of the cover crop is required, acknowledging its importance to the overall return of the rotation. To determine suitable cover crop mixtures and management, more interdisciplinary projects are required.

  • In order to advance our understanding of cover-crop-related effects, we suggest for future research that comparability among studies through the inclusion of appropriate controls and additional data be improved, i.e. biomass and P content of cover crops and main crops, as well as soil biological parameters and soil P pools.

Bearing these suggestions in mind, scientists will be able to join farmers in moving towards soil improving cropping systems, through relevant research on the benefits of cover cropping.