![]() |
![]() |
doi:10.1016/j.agee.2008.06.005 | How to Cite or Link Using DOI |
Permissions & Reprints |
Received 29 January 2008; revised 29 May 2008; Accepted 2 June 2008. Available online 25 July 2008.
Abstract
Soybean (Glycine max) is a booming crop in Brazil. In 2004, the export value was equivalent to 10 billion US $, covering over 10% of total Brazilian exports. Three-quarters of total production leaves the country, mainly to China and the European Union (EU). Soybean cultivation in Brazil is expected to expand further in the coming decades, mainly responding to growing demand in Asia. This will, amongst others, entail transport of vast amounts of nutrients, triggering the need to better study the entire soybean chain. The objective of this study was to estimate and calculate the soybean chain, including five phases: conversion, cultivation, transport and processing, consumption and waste disposal, starting in Brazil, and ending in Brazil, China and EU, using nitrogen (N) as a marker, and looking at three time periods (1993–1995; 1998–2000; 2003–2005). The study revealed that conversion of forest and savanna to pasture and agricultural land entails N losses of 2000–6000 million kg year−1. Removal of N in soybean harvests went up from 1400 million to almost 3000 million kg year−1 between 1993–1995 and 2003–2005. These high values were offset by biological N fixation by soybean and increased adoption of conservation agriculture. N balances in soybean-based agricultural systems became positive after about one decade in the period 2003–2005, thus reducing the soybean-associated global N cascade. Upon crushing, three-quarters of soybeans end up as high-protein soy meal, which is mainly fed to pigs and chickens. Nitrogen in meat, milk and eggs from soy meal-fed animals was estimated at around 20% of N in freshly crushed soy meal. More than half of the lost N can potentially be recycled, although mostly far away from the site of soybean production.
Keywords: Soil fertility; Soybean; Soy meal; Nitrogen stocks and flows; Forest conversion; Conservation agriculture; Animal production; Human consumption; Waste; Brazil; China; European Union
Article Outline
- 1. Introduction
- 2. Materials and methods
- 2.1. Spatio-temporal database
- 2.2. Conversion
- 2.3. Cultivation
- 2.3.1. Nutrient balance
- 2.3.2. Tillage systems
- 2.3.3. IN1 Mineral fertilizer
- 2.3.4. IN2 Organic inputs
- 2.3.5. IN3 Wet and dry deposition
- 2.3.6. IN4 Biological N fixation
- 2.3.7. OUT1 Harvested crop parts
- 2.3.8. OUT2 Removed crop residues
- 2.3.9. OUT3 Leaching
- 2.3.10. OUT4 Gaseous losses
- 2.3.11. OUT5 Erosion
- 2.4. Transport and processing
- 2.5. Consumption
- 2.6. Human waste disposal
- 2.7. Uncertainty analysis
- 3. Results
- 3.1. Conversion
- 3.2. Cultivation
- 3.3. Transport, processing and consumption
- 3.4. Waste disposal
- 4. Discussion
- 5. Conclusions
- Acknowledgements
- References
1. Introduction
With economic growth, urbanization and changing diets, world demand for plant-derived oils and their derivates is soaring. Oils are used inter alia in the food, feed and cosmetics industry, and increasingly as a biofuel. Soybean is one of the major booming oil crops. In addition to oil, it offers a very valuable by-product, i.e. protein-rich soy meal, a raw material for animal feed. Increasing global demand for animal products puts pressure on arable land, given the fact that at least two but often more kilograms of grains are needed to produce one kilogram of meat, i.e. the so-called feed meat ratio (Keyzer et al., 2005), and that an array of environmental problems is looming due to the expanding livestock sector (FAO, 2006). The demand surge largely stems from China. Table 1 shows a ninefold increase in soy imports in 10 years, versus an insignificant home production. This was triggered by China's 2002 WTO membership, which put an end to border tariffs and boosted trade (Van Berkum et al., 2006). Increasing demand has been met with a supply response that is particularly strong in Brazil and Argentina. Country statistics show that in 2004, Brazilian production exceeded 50 million Mg, twice the amount realized in 1997. Area and production increases were particularly strong in the period 2001–2005, following a favorable devaluation of the national currency, the real. Soybean exports in 2004 earned Brazil over 10 billion US $, against 4.2 billion US $ in 2000.
1994 | 1999 | 2004 | |
---|---|---|---|
Production (Mton) | 15 | 15 | 17 |
Area (Mha) | 10 | 9 | 10 |
Export (Mton) | 540 | 214 | 346 |
Import (Mton) | 2654 | 8236 | 24,863 |
Net import (Mton) | 2113 | 8022 | 24,517 |
During the 1960s, soybean cultivation was concentrated in the three states of the Southern Region of Brazil: Rio Grande do Sul, Santa Catarina and Paraná. Later on, expansion to the Central Cerrado Region took place once varieties had been developed that were adapted to low latitudes. Currently, 60% of soybean production is concentrated in the Cerrado Region. The southern area is characterized by smallholder farmers, mostly organized in cooperatives, whereas the central area is characterized by large holdings, with very high levels of mechanization, and mostly organized in large private groups. In addition, climate is more stable in Central Brazil. The latest soybean expansion has taken place in the northern Legal Amazon. Soybean production in the North represents only a small fraction of the currently planted area, but is expanding.
Recognizing that soybean (i) is a highly important commodity for the Brazilian economy, (ii) contributes to the conversion of forest and savanna land, and (iii) travels long distances, ending up in many different foods and feeds worldwide, there is a need to better understand the dynamics of the soybean sector, and how economic and environmental targets can be realized simultaneously. In this study, nitrogen (N) is used as a marker. The important role of N for life, for growth of crops and as an essential element in aminoacids/proteins in our food is obvious. The final judgment is about the trade offs between overloading and loss of N versus the benefits in our food of plant and animal origin. Many studies have looked into the N cycle and the flow of N in the food system, at the global scale ( [Galloway and Cowling, 2002] , [Galloway et al., 2004] , [Galloway et al., 2007] and [Smil, 2002] ), or the national scale (e.g., Antikainen et al., 2005). In this paper we look more specifically into the role of N from soybean in the food system. Soybean is a leguminous crop, capable of fixing atmospheric N through a symbiosis with Rhizobium bacteria in the plant's root nodules. Soybean crops therefore receive no or very small amounts of N fertilizer as a starter.
The objective of this study is to quantify N stocks and flows in and between the compartments of the Brazilian soybean chain, and to find out where in the process N either remains inside or disappears from the food chain. These compartments include forest and savanna conversion, soybean cultivation, transport and processing, animal and human consumption, and waste disposal (Fig. 1).
2. Materials and methods
2.1. Spatio-temporal database
In Brazil, soybean is grown in large and small holdings, and both as a monoculture and in annual rotations with maize, wheat and other crops. For characterizing the soybean farming systems, spatially explicit information was used. To show trends and avoid outliers, average data were taken for the years 1993–1995, 1998–2000, and 2003–2005. Harvested area and production data were obtained from municipality statistics (IBGE, 2006) and related to the country-wide spatial municipality database (n = 5564 for 2003–2005). Besides production data, several spatial data sets were necessary for the N balance calculation. Soil properties were taken from 1:5-million-scale SOTER and WISE databases (FAO et al., 1998; [7] and [8] ). Texture, bulk density, total N and total carbon content of the 0–20 cm soil layer were included for calculation of the N balance. Rainfall data of the CRU TS 2.10 database were used, at a resolution of 0.5 degrees (Mitchell and Jones, 2005). For the calculation of erosion, the slope gradient was required, which was derived from the GTOPO30 digital elevation model (USGS, 1996), with a resolution of 30″ (approximately 1 km). The maps were overlaid in a Geographic Information System, and the mean soil properties, rainfall and slope were calculated for each municipality. The datasets allowed the mapping of soybean land for each time period (Fig. 2a–c), as well as the particular system, i.e. soybean monoculture, soybean/maize and soybean/wheat (Fig. 3a–c).
2.2. Conversion
The expansion of soybean land happens at the expense of pastures, which in turn farmers developed from Cerrado savanna and now increasingly from more densely wooded savanna. Data on decreasing forest and savanna area cannot, however, be fully explained by soybean expansion. To calculate N flows under conversion regimes, a literature review was made. Removals of N are in timber, firewood, in slash and burn of remaining above-ground biomass, and in soil and litter. Land clearing methods range from selective cutting to the destructive savanna clearing system, where a correntao, a giant chain is attached between two bulldozers that move jointly, uprooting not only all vegetation but also topsoil and litter. The totality of this material is then burned (Van Gelder and Dros, 2005).
Mielniczuk et al. (2003) compared data on soil carbon inputs and stocks in long-term experiments in Sweden (Paustian et al., 1992) and in the South of Brazil (Lovato, 2001). In equilibrium, the decomposition rate was 1.05% year−1 in the temperate soil and 3.2% year−1 in the subtropical soil. The carbon to nitrogen ratio of organic matter of Southern Brazilian soils is on average 12.5:1 (Jantalia et al., 2006). A woodland savanna (Cerrado stricto sensu) stores 229 t of C ha−1 in biomass plus soil (1 m depth; Roscoe, 2005). Texture is also an important factor determining carbon stocks in Cerrados soils. Tognon et al. (1998) found a strong linear relationship between texture and carbon stored in soil, varying from 75–87 t of C ha−1 in sandy soil to 182–188 t of C ha−1 in clayey soils. Considering a C:N ratio of 12:1, this would represent an N stock of about 7 t ha−1 in sandy soils and about 15 t ha−1 in clayey soils. Pastures normally develop higher soil N contents than forest and savanna. Neill and Davidson (1999) showed that, from 29 studies reviewed, 19 forest soils had increased their organic matter content in surface when converted into pastures. Other literature data revealed that values for carbon in the upper 1 m of Brazilian forest and savanna soils range from 100 to 200 t ha−1 ( [Matson et al., 1987] , [Cerri et al., 1991] and [Da Silva, 2004] ). Carbon in vegetation ranged from 30–80 t ha−1 in savanna to 140 t ha−1 in forest (Graca et al., 1999). Total N losses in the forest system following clearcutting, based on review data by Matson et al. (1987) are in the order of 1000 kg N ha−1, but data from Cameroon (Kanmegne, 2004) and Madagascar (Brand and Pfund, 1998) show losses as high as 3000 and 5000 kg ha−1 for forest. Bouwman (1995) found in La Selva, Costa Rica, soil C and N stocks of 81,000 and 7600 kg ha−1, respectively. Soil N losses due to decomposition of soil organic matter following forest conversion to pasture were 24% after 3 years, and after 5 years 40%. Schlesinger (1986) gives values of 21% loss of soil C when converting land from tropical rainforest to agriculture (based on 19 sites). Based on this survey, it is assumed that forest and savanna have a ‘starting N stock’ in vegetation and soil, whereas the situation prior to soybean cultivation is the ‘final N stock’, which is the input value for the ‘Cultivation’ phase (Fig. 1). Pasture is considered an intermediate phase. Loss rates are then estimated on the basis of actual deforestation data, and on the percentage loss that can be attributed to soybean expansion. Based on the literature study, and assuming a soil C:N ratio of 12:1, and a relevant topsoil thickness of 20 cm, an initial N stock in forest and savanna systems of 8000–12,000 kg N ha−1 seems plausible, savanna being at the lower end and forest on the upper end. Moreover, it is assumed that 25–50% of forest/savanna system N is lost in the process of conversion to pasture and then to soybean.
2.3. Cultivation
2.3.1. Nutrient balance
The soybean system was analysed, largely using the NUTMON nutrient balance approach ( [Stoorvogel et al., 1993] , [FAO, 2004a] and [Lesschen et al., 2007] ). Table 2 lists the relevant flows. The nutrient balance calculation involved a full cropping cycle of soybean (monoculture) or soybean + maize or wheat (rotation), and the combined balance was thus calculated using the planted areas of each soybean system.
Inputs | Outputs | ||
---|---|---|---|
IN1 | Mineral fertilizer | OUT1 | Harvested crop parts |
IN2 | Organic inputs | OUT2 | Removed crop residues |
IN3 | Wet and dry deposition | OUT3 | Leaching |
IN4 | Biological N fixation | OUT4 | Gaseous losses |
OUT5 | Erosion |
2.3.2. Tillage systems
Agriculture in the 1960s and 1970s in the South was heavily mechanized, leading to strong soil degradation ( [Mielniczuk et al., 2003] and [Jantalia et al., 2006] ). There was a widespread concept among farmers that well-managed soils would have completely clean and smooth surfaces, with residue incorporation and pulverized soils (Jantalia et al., 2006). Comparing the dynamics of soil organic matter in a long-term experiment in the South of Brazil described by Lovato (2001), Mielniczuk et al. (2003) estimated that no-tillage reduced the decomposition rate from 3.2% year−1 to 1.7% year−1. In the Cerrado region of Bahia State, Silva et al. (1994) sampled 220 topsoils (0–15 cm) from soybean fields under conventional tillage and monoculture. The authors observed severe losses from 41% (clayey soils with >30% clay) to 80% (sandy soils with <15% clay) of the initial soil organic matter contents after 5 years of cultivation. However, Freitas et al. (2000) and Roscoe and Buurman (2003) did not observe changes in the SOM stocks (0–40 cm) of a clayey Dark Red Latosol after 25 and 30 years of maize–bean successions with conventional tillage. Llilienfein and Wielcke (2003) reported no significant changes in C content of a clayey Oxisol after 12 years of maize–soybean rotation under conventional tillage. The presence of a high-residue crop such as maize in these experiments may partly explain the persistence of organic matter in these soils. However, the most important factor seems to be the high protection of organic matter in clayey Cerrado soils, which are rich in iron and aluminum hydroxides (Resende et al., 1997).
2.3.3. IN1 Mineral fertilizer
The fertilizer use by crop study for Brazil (FAO, 2004b) estimated fertilizer consumption values per crop and per region (Table 3). These values were estimated for 2002 and were multiplied with a correction factor, based on the relative difference between the fertilizer consumption of Table 3 and the IFA national fertilizer statistics for the calculation years 1993–1995, 1998–2000 and 2003–2005. Total N consumption for these periods was derived from the FAO fertilizer statistics, the values being 1,148,600 Mg N for 1993–1995, 1,738,100 Mg N for 1998–2000 and 2,415,400 Mg N for 2003–2004.
Region | Soybeans | Maize | Wheat |
---|---|---|---|
North | 2 | 10 | 0 |
Northeast | 4 | 22 | 0 |
Central | 7 | 40 | 9 |
Southeast | 7 | 43 | 9 |
South | 9 | 53 | 12 |
2.3.4. IN2 Organic inputs
Although Brazil has the largest livestock herd in the world, as well as large numbers of pigs, sheep and poultry, application of organic manures is limited to horticulture and perennial crops located close to the production sites. In the case of grain crops, the use of organic fertilizers is not very common (FAO, 2004b). Therefore, the organic inputs for soybeans, maize and wheat were set at zero.
2.3.5. IN3 Wet and dry deposition
Trebs et al. (2005) measured dry and wet N deposition for a remote pasture site in the Amazon Basin, and found 7.3–9.8 kg ha−1 year−1. Krusche et al. (2003) also estimated an N deposition rate of about 9 kg ha−1 year−1. According to two chemistry transport models ( [Collins et al., 1997] and [Dentener et al., 2006] ), N deposition for Brazil amounts to 5–10 kg ha−1 year−1. Based on this, input through N deposition was set at 8 kg ha−1 year−1.
2.3.6. IN4 Biological N fixation
Biological N fixation (IN4, kg N ha−1 year−1) is the most important input for a leguminous crop such as soybean. Based on studies of Zotarelli et al. (1998), [Alves et al., 1999] and [Alves et al., 2005] , Hungria et al. (2006) and Araújo et al. (2006), Rhizobium-based N fixation by soybeans was set at 80% of the N in the above-ground plant parts. For non-symbiotic N fixation (based on the input by free living bacteria), the following equation applies (FAO, 2004a):IN4=0.5+0.1×√rainfall
2.3.7. OUT1 Harvested crop parts
Production and harvested areas for soybean, maize and wheat were derived from the PAM (Agricultural Municipalities Production) database (IBGE, 2006). The N content of the harvested crop parts for soybeans was set at 58 g kg−1 harvest product, based on studies by Zotarelli et al. (1998), [Alves et al., 1999] and [Alves et al., 2005] , Hungria et al. (2006) and Araújo et al. (2006). For maize and wheat, the values of FAO (2004a) were used, i.e. 17 and 22 g kg−1, respectively.
2.3.8. OUT2 Removed crop residues
The removal of crop residues from the field is not a common practice in soybean farming systems in Brazil, as almost no alternative uses for them, such as stall-feeding, occur. Therefore, we assumed that all crop residues remain on the field and, therefore, OUT2 was set at zero.
2.3.9. OUT3 Leaching
For N leaching (OUT3, N in kg ha−1 year−1), the regression model of De Willigen (2000) has been used:where P is the annual precipitation (mm); C the clay (%); L the layer thickness, i.e. rooting depth (m); F the mineral and organic fertilizer nitrogen (kg N ha−1 year−1); D the decomposition rate (% year−1); NOM the amount of nitrogen in soil organic matter (kg N ha−1); U the uptake by crop (kg N ha−1 year−1).
The first part of the regression equation determines which fraction of mobile N will be leached, and the second part determines how much mobile N is available. To prevent overestimation of leaching, the first part of the equation was maximized at 1. In the case of soybeans, only the uptake from the soil is taken into account (20%) and not the part which is derived from N fixation. For rooting depth, values of Allen et al. (1998) have been used, which are 0.6 m for soybeans, 0.9 m for maize, and 1.0 m for wheat. Decomposition rate was set at 3.2% year−1 for conventional tillage and 1.7% year−1 for no-tillage (Mielniczuk et al., 2003).
2.3.10. OUT4 Gaseous losses
The regression model developed by FAO (2004a) has been used with an additional factor for crop residues. NH3 volatilisation from field crop residues is set at 4%, which is the volatilisation rate of some ammonium-based fertilizers (Bouwman et al., 2002).OUT4=(0.025+0.000855×P+0.01725×F+0.117×O)+0.113×F+0.04×Rwhere P is the annual precipitation (mm year−1); F the mineral and organic fertilizer nitrogen (kg N ha−1 year−1); O the organic carbon content (%); R the nitrogen in remaining crop residues (kg N ha−1 year−1).
2.3.11. OUT5 Erosion
The RUSLE model was used (Renard et al., 1997), in combination with the spatial rainfall and soil data discussed in Section 2.1. The RUSLE is expressed as:A=R×K×LS×C×Pwhere R is the rainfall erosivity factor (−); K the soil erodibility factor (kg m−2); LS the slope length and slope steepness factor (−); C the crop management factor (−) and P the erosion control practice factor (−).
The R-value was based on the annual erosivity map for Brazil by Da Silva (2004). The K-factor was based on the nomographs of Wischmeier and Smith (1978), but a simplified version was used because of a lack of soil permeability and soil-structure data. The LS factor was calculated using the slope gradient and a default slope length of 22.13 m. According to Morgan (1995), the C-factor for soybeans ranges from 0.20 to 0.50. As the C-factor is significantly influenced by tillage, a distinction was made between conventional tillage and no-tillage. Erosion on no-tillage fields is on average only 25% of erosion under conventional tillage (Roscoe, unpublished data). As a consequence, a C-factor of 0.40 was used for conventional tillage and 0.10 for no-tillage. Table 4 gives the percentages of agricultural area under no-tillage (Roscoe, unpublished data). Besides no-tillage management, no other erosion control practices are used on a large-scale, hence the P-factor is 1. Finally, the loss of N was calculated by multiplying the annual soil loss by the soil N content and an enrichment factor, which was set at 2.3 (FAO, 2004a).
Region | 1993–1995 | 1998–2000 | 2003–2005 |
---|---|---|---|
North | 0 | 10 | 30 |
Northeast | 0 | 10 | 30 |
Central | 5 | 53 | 77 |
Southeast | 21 | 74 | 85 |
South | 21 | 74 | 85 |
2.4. Transport and processing
Once harvested, the soybeans move in different directions (Fig. 4). The part undergoing processing inside Brazil (56%) is crushed to produce soy oil and the protein-rich by-product soy meal. Upon crushing, the chain is only continued for soy meal, soy oil not containing N anymore. Then, 13% of the soy meal is used as feedstuff in Brazil, whereas 30% is exported. Part of Brazilian animal produce is exported (4%). Further assumptions necessary for the calculations of N destinations are given in Table 5.
Process | Remarks/reference |
---|---|
Starting N stock | N in total harvested soybean for the three time intervals (“Cultivation”) constitutes the input N-value for “Transport and Processing” |
Transport and processing losses | Seed use and transport losses between field and factory are estimated at 5%; losses during crushing and further transport are assumed to be 5%. Total losses and seed use between the field and the livestock producers are thus 10%, which is consistent with estimates from EMBRAPA (Maria do Rosario, personal communication) |
Products following crushing | Crushing beans is assumed to produce 20% oil, 75% meal, and 5% processing losses (weight percentages) |
Import/export | Based on the growing importance of China, it is assumed that in 1993–1995, 10% of exported soy went to China and 90% to EU, whereas in 1998–2000 this ratio is assumed to be 30:70%, and in 2003–2005 50:50% |
2.5. Consumption
For the (animal and human) consumption phase, a series of assumptions was made relating to feed conversion efficiency in animals, slaughterhouse efficiency, and human consumption characteristics, including wasting of food. The N losses are in fact not all ‘losses’. Considerable percentages leave the ‘useful product’ chain but can be recycled. This potential re-use is, however, often far away from the site where soybeans are produced, hence causing an accumulation of N at consumption sites. The assumptions are summarized in Table 6.
Process | Remarks/reference |
---|---|
N conversion efficiency in pork and poultry production | Soybean meal is assumed to be completely used in pork and poultry production. Conversion efficiencies (the percentage of the N in feed that is converted to animal tissue) for Brazil and China are assumed to be 21% for pork and 34% for poultry production. For EU, these numbers are 30% for pork and 34% for poultry based on Van der Hoek (1998). The complement of the N conversion is excreted by the animals |
N losses from animal manure | Pork and poultry production is assumed to take place in stables. Ammonia (NH3) losses from animal manure in stables is assumed to be 20% (Bouwman et al., 2002); further losses during the spreading of manure are not taken into account in this analysis |
Recycling of manure | Recycling of pig and poultry manure as organic amendment in agriculture is assumed to be 80% in China and 90% in European countries throughout the period covered by the analysis (Smil, 1999). For Brazil, we assumed that 10% of the manure is not recycled |
N in meat | The N in the consumable meat is calculated from data provided by FAO (2007) on the live weight, dressed carcass weight, retail cuts (52% of live weight for pork and 48% for poultry), and various edible offals (4% of live weight) and non-edible by-products; edible by-products are assumed to be recycled, and non-edible by-products are assumed to be taken out of the food cycle. Protein content of the meat (160 mg g−1 fresh meat for pork and 180 mg g−1 for poultry) and offals (∼180 mg g−1 for both pork and poultry meat) are from FAO (2007); that for bones (∼50 mg g−1 for pigs and 30 for poultry) is from Aerssens et al. (1998). Thus, total N in the live weight is distributed over consumable meat (82% for pork and 85% for poultry) and offals (∼7% for both pork and poultry) |
Consumptive loss | Consumptive loss is defined as the total loss in retailing and wholesaling and consumer and food service losses (from foods forgotten and spoiled in the refrigerator to the uneaten food tossed in the garbage); based on inventories by Kantor et al. (1997) and Bleken and Bakken (1997) these losses are substantial. For the USA, 15% of the meat is lost this way. For Europe, we assumed a 10% loss, for developing countries 5% |
2.6. Human waste disposal
The final part of the chain looks into the human excreta. In our calculations we consider all N in human waste as potential emission to the environment, excluding the N removed in wastewater treatment systems. This is because there is no quantitative information on the fate of the N for other systems such as latrines, septic tanks and other less institutionalized human waste disposal systems. For example in China, many of the human excreta have been recycled, and used to fertilize agricultural systems until the late 1970s (FAO, 1977). However, it is not well known if this practice is still important at present (Zhu and Chen, 2002). Hence, by ignoring recycling of human excreta, we may overestimate the potential emission of human N to the environment.
Depending on the country, part of the population is connected to sewerage systems, and part of the sewage water is treated in wastewater treatment plants. The N removal in treatment installations depends on the type of installation. Assumptions are summarized in Table 7.
Process | Remarks/reference |
---|---|
N excretion by population | Total N excretion is calculated as total intake of N from soybeans through meat consumption. We accounted for N accumulation in the growing human population. Population data are from FAO (2007) |
Human wastewater to sewerage and surface water | Fraction of the population connected to sewerage systems, and N removal in wastewater treatment systems were updated from Bouwman et al. (2005) |
2.7. Uncertainty analysis
An uncertainty analysis was applied to gain insight in the accuracies and possible ranges of the results. For the cultivation phase, we followed the procedure of Lesschen et al. (2007), in which uncertainties were first calculated or estimated for each flow at municipality level (Table 8). Since the results were aggregated to regional and national level the uncertainties have to be adjusted for spatial correlation. Spatial correlation occurs because part of the input data is spatially dependent, e.g., soil and rainfall data. The spatial correlation coefficient was estimated for each nutrient flow and ranged between 0 and 1.
Flow | Variablesa | Calculation | Relative uncertainty (%) |
---|---|---|---|
IN1 | Fertilizer consumption (±10%) | Product | 51 |
Fertilizer use per crop (±50%) | |||
IN3 | Literature values ranging between 6 and 10 kg N/ha | 25 | |
IN4 | Symbiotic N-fixation soybean (±10%) and crop uptake (±22%) ⇒ 24% | Summation | 25 |
Rainfall map (±20%) and regression (±50%) ⇒ 54% | |||
OUT1 | Crop production data (±20%) | Product | 22 |
N content crop products (±10%) | |||
OUT3 | Rainfall map (±20%) and clay content (±20%) ⇒ 28% | Product and summation | 81 |
Fertilizer application (±51%), decomposition rate (±50%), soil N content (±30%) and N uptake (±22%) ⇒ 50% | |||
Regression model R2 = 0.67 (57%) | |||
OUT4 | Rainfall map (±20%), fertilizer application (±51%) and organic carbon content (±20%) ⇒ 35% | Product and summation | 65 |
Regression model R2 = 0.70 (55%) | |||
OUT5 | R: rainfall erosivity factor (±20%) | Product | 81 |
K: soil erodibility factor (±55%) | |||
LS: slope and slope length factor (±50%) | |||
C: crop management factor (±25%) |
For Phase 3 (transporting and processing) and Phase 4 (consumption), a simple sensitivity analysis was done by varying a number of parameters by 25% and calculating the effect on the soybean-derived human N intake.
3. Results
3.1. Conversion
Fig. 5 shows the cumulative area of deforestation, taken from the National Institute of Space Research (INPE), and the total soybean area. The area increased from 11.3 million ha in 1993–1995, to 20.1 million ha in 2003–2005. This increase occurred mainly in the central and northern regions. Until 2000, there is no relation between deforested area and soybean area, but between 2000 and 2004, loss of forest is approximately equal to expansion in soybean area. In reality, the conversion rate of forest/savanna to pasture may, in these years, be approximately equal to the conversion rate between pasture and soybean land. However, these data do not fully explain whether most of the new soybean areas are savanna-derived or forest-derived conversions. Also, soybean expansion has always taken place in land under other arable crops as well.
Taking the rather constant deforestation rate of 2 million ha year−1 (Fig. 5), and attributing 50% of deforestation to soybean cultivation, N losses during Phase 1 are in the order of 2000–6000 million kg year−1 on the basis of the N losses presented in Section 2.2, with 1993–1995 at the lower and 2003–2005 at the higher end of the range. It is difficult to attribute exactly the conversion to vegetation types and time periods, but it is obvious that with the advance of soybean into the central and northern parts of the country, soybean expansion takes place increasingly at the expense of forest.
3.2. Cultivation
Soybean production data have been summarized in Table 9. OUT1 is given for the overall soybean-growing area, also including outputs in maize and wheat. The final row gives the total N production through harvested soybeans, assuming N content in beans of 5.8%. The uncertainty of the total amount of N from soybean production is only 7%, based on the 22% uncertainty for OUT1 and a spatial correlation factor of 0.3.
1993–1995 | 1998–2000 | 2003–2005 | |
---|---|---|---|
Harvested soybean area (million ha) | 11.3 | 13.4 | 20.1 |
Soybean yield (Mg ha−1) | 2.16 | 2.37 | 2.53 |
OUT1—crop products (kg N ha−1) | 101 | 115 | 122 |
Total soybean N (in 1000 Mg) | 1420 ± 94 | 1843 ± 122 | 2949 ± 195 |
N flows were calculated for the different regions and different soybean farming systems portrayed in Fig. 2 and Fig. 3. Fig. 6 shows that surface N balance deficits decrease with time and become surpluses, primarily because of the increased importance of no-tillage, reducing leaching and erosion. Average uncertainty of the N balance is ±9 kg N ha−1, ranging between ±6 kg N ha−1 for the Southeast region and ±16 kg N ha−1 for the North region as calculated for the 1998–2000 period. Soybean is estimated to fix 80% of its needs, and since the crop residues remain in the field, the values of OUT1 and IN4 are about equal (Fig. 7). Based on the total soybean area and N balances calculated for the three time periods in the “Cultivation” phase, there is a total loss of 182 million kg N year−1 in 1993/1995, a loss of 15 million kg N year−1 in 1998/2000, and a gain of 62 million kg N year−1 in 2003/2005.
3.3. Transport, processing and consumption
Based on the assumptions listed in Table 5 and Table 6, flows were calculated during processing, transport and animal and human consumption. Fig. 8 shows the values and flows representing the relevant processes for the three periods. Starting stocks correspond with values given in Table 9, i.e. the output of ‘Cultivation’. About 20% of total N entering Phases 3 and 4 ends up as animal products that are consumed by humans. This percentage decreased slightly between 1993–1995 and 2003–2005 (from 21.1 to 19.8%). In addition, Table 10 shows total N ingested by humans per country/region. Adding up the rows gives totals of 300, 377 and 583 million kg N year−1 for the three periods.
Brazil | China | EU | ||||
---|---|---|---|---|---|---|
Poultry | Pork | Poultry | Pork | Poultry | Pork | |
1993–1995 | 57 | 13 | 6 | 13 | 78 | 134 |
1998–2000 | 74 | 16 | 23 | 50 | 79 | 135 |
2003–2005 | 118 | 26 | 61 | 134 | 90 | 154 |
Table 11 provides a sensitivity ranking for Phases 3 and 4, i.e., on final soybean-derived N intake from meat due to variation of the different input parameters by +25%. It is clear that the feed N conversion in livestock production is the major uncertainty in Phases 3 and 4. This is related to the large fraction of soy meal that goes into livestock production, and the generally low efficiency of feed conversion to meat.
Parameter | Change of N intake (%) | Remarks |
---|---|---|
Transport/processing loss | −2.6 | |
Fraction production to industry within Brazil | +1 | More export of meal, with on the average higher efficiency |
Fraction consumption of processed beans in Brazil | −1.5 | More export of meal, with on the average higher efficiency |
Fraction export of meal to Europe vs. China | +2.1 | |
Use of soybean meal for pigs vs. poultry | −6 | |
N conversion pigs | +13.5 | |
N conversion poultry | +11.5 | |
Fraction offals edible | −2.1 | |
Fraction non-edible offals | −2.7 | |
Fraction retailing, wholesaling and household food waste | −2.7 |
3.4. Waste disposal
Table 12 shows the fate of consumed N. Total consumption (583 million kg N year−1) is the output value of Phase 4 for 2003–2005. European Union (EU) households have by far the largest percentage of sewerage connections (84%). However, even with advanced wastewater treatment in Europe with 51.5% of N in influents being removed, 48.5% of the N entering sewage systems still reaches surface waters. If we consider all households in the EU including those lacking sewage connection, 57% of the soybean N consumed is not removed in wastewater treatment.
Brazil | China | EU | Total | |
---|---|---|---|---|
Consumption | 144 | 195 | 244 | 583 |
Population change (% year−1) | 1 | 1 | 0 | |
N excretion | 142 | 193 | 244 | 580 |
Connection to sewerage (%) | 44 | 19 | 84 | |
Removal in wastewater treatment (%) | 18 | 10 | 51.5 | |
N emission (potential, see text) | 131 | 189 | 138 | 459 |
In Brazil and China less than 10% of the soybean N is removed from sewage water, hence more than 90% is not treated and is a potential emission to the environment, since far fewer households have a sewerage connection than in the EU. The fate of N in non-sewerage outlets, ranging from latrines to septic tanks, has not been considered explicitly due to lack of data. Human excreta that are collected may be used to fertilize agricultural fields, particularly in China. Hence, the potential N emission to the environment may be an overestimation.
4. Discussion
4.1. Limitations of the study
Given the booming market, expansion of the soybean area in Brazil seems economically viable, but it has a social and environmental price. Many negative externalities of the soybean sector are not covered in this study, including plant and animal biodiversity loss, carbon emissions, changing air circulation and increased drought occurrence, unlawful land acquisition, forced migration of forest inhabitants, rural poverty and low employment rates.
On ‘Conversion’, studies exist that look at vegetation carbon, and at soil carbon available in forest and savanna. Fewer look at N, let alone both vegetation N and soil N. Some studies take soil N down to a depth of 1 m, others just consider 0–20 cm. Land can be cleared in many different ways. The fully fledged bulldozer-assisted clearing of cerrado mentioned by Van Gelder and Dros (2005) is much more destructive than the selective cutting that takes place elsewhere. For this study, ‘clearcutting’ is considered to take place, although the actual process may include periods of partial clearing as well as periods under pasture, which often show increases in soil C and N over time, particularly when managed well.
On ‘Cultivation’, N is clearly not a limiting factor as soybean is a leguminous crop, but other nutrients may well be or become limiting. FAO (2004b) shows that soybean in Brazil normally does receive P and other fertilizers, next to tiny amounts of N. These other nutrients are not considered in this study. No attention is paid either to the fact that, currently, most soybean varieties are genetically modified (Brazilian Ministry of Agriculture estimates are 70% of the soybean area under GM soybean in 2005/2006), and based on ‘Round-up ready’ herbicide use, which are imperative for adopting minimum tillage practices. A social externality not taken into account is that the labour/land ratio of large-scale soybean production in central Brazil is very low (1 on 400 ha) compared with small-scale agriculture in the south (80 on 400 ha) (Van Berkum et al., 2006).
Another issue outside the scope of the current study is the increased use of soy oil as a biodiesel. A new market was opened to soybean oil in 2006/2007 as a consequence of the Brazilian National Biodiesel Programme. The programme has allowed the inclusion of 2% of biodiesel in the diesel from petroleum since January 2006. This proportion will become compulsory in 2008 and will be increased to 5% in 2013. The demand for vegetable oil in 2008 is estimated about 1 million t and, in 2013, 2.5 million Mg. According to Petrobrás projections, in 2008, the demand for vegetable oil in this process will be about 500,000 Mg. Therefore, the energy sector will absorb about 1.5–3.0 million Mg of vegetable oil in the coming 2–5 years.
On ‘Transport’, a recently opened port in Santarem and the completion of the road between Cuiaba and Santarem opens up new opportunities for agricultural development, but it may equally pose substantive social, environmental and production risks. Van Berkum et al. (2006) show that export and import figures can differ markedly. For example, ITC/WTO data show that on soy meal trade between Latin America and The Netherlands, Brazil exported 4 million Mg of soy meal, but The Netherlands only imported 1.3 million Mg. Similarly, Argentina export data show 2.3 million Mg shipped to The Netherlands, but the Dutch statistics show imported soy meal from Argentina to have amounted to just over 1.1 million Mg. Although ships do at times change destination (e.g. to France, Belgium or Germany) while moving overseas, the differences are significant.
On ‘Consumption’, there are several causes of uncertainty. The most uncertain parts of the consumption chain are the feed N conversion in livestock production systems and food spoilage and waste. The process that determines a large part of the N loss from the soybean chain is the feed conversion. A small variation causes a large deviation in the end results (Table 11). Estimates for feed N conversion are from large-scale studies, and differences between the estimates for Brazil and China (21% efficiency for pigs) and Europe (30%) indicate that management has an important effect on N lost from the soybean chain. The results shown here show greater efficiency of N capture by animal production than reported previously ( [Smil, 1999] and [Smil, 2001] ). In the previous analysis, the effect of N recycling within the system was not considered, and older data were used to estimate the efficiency of feed N utilization by animals. Since livestock production systems are rapidly changing, particularly in developing countries, where modern units with high levels of management replace existing ones, efficiencies may also increase rapidly in countries like Brazil and China.
In our study roughly 20% of the N in the soybeans is actually consumed as meat by humans. Since we consider only pork and poultry, our estimate exceeds the 14% estimated by Smil (2002) for all livestock including milk and beef cattle.
Estimates of food availability for human consumption from food balance sheets are generally used as an approximate level of food actually consumed. This may work reasonably well in developing countries, but in developed countries this approach can overstate the level of consumption because the amount of food spoilage and waste in catering establishments is rather high. The amount of food actually consumed may be lower than the quantity shown in the food balance sheet depending on the degree of losses of edible food and nutrients in the household, e.g. during storage, in preparation and cooking (which affect vitamins and minerals to a greater extent than they do calories, protein and fat), as plate-waste or quantities fed to domestic animals and pets, or thrown away. Household survey data may prove useful in such instances for adjusting the waste component in the food balance sheet.
On ‘Waste Disposal’, other defecation destinations than sewerage and water treatment are not explicitly dealt with in this study. They may be equally or more constraining to public health, but can be considered as potential emissions of N compounds into the environment.
Quantification of the fate of the N that is lost from the chain and assessment of the environmental impact is outside the scope of this paper. This would require a completely different approach like life cycle analysis (LCA). In actual fact, results of our study may serve as a basis for LCA of the soybean food chain. In this respect, the intention is also not to provide an exhaustive overview of N stocks and flows as was done for Finland by Antikainen et al. (2005).
4.2. Opportunities to make the chain more sustainable
Major reductions in N losses are possible in Phases 1, 4 and 5 (Fig. 1). Although no research was done on solutions, some opportunities stand out that may contribute to striking a more sustainable balance between economic and environmental goals.
On ‘Conversion’, improvements and solutions are more in the field of forest management policies, law enforcement, and human rights protection than in agronomy. Payment for environmental services along the lines of Kyoto Protocol mechanisms may be a way to harmonize economic, environmental and social goals at national and local scale.
As to ‘Cultivation’, Brazil's National Agricultural Research Institute (Embrapa) has indicated that considerable strides can still be made to raise productivity on existing soybean land, to further spread the practice of no-tillage, to restore and cultivate degraded pastures, and to rotate pastures and soybean. Estimates on ‘degraded’ land that could be converted to soybean are between 10 and 30 million ha.
Although trade generates economic benefits, soybeans and soy meal face ‘Transport’ over very wide distances before entering the consumption cycle. Atmospherically fixed N from Brazil ends up in dung piles in Western Europe and China. In the Netherlands, for example, this has for some decades now been regarded as a major environmental problem.
The consumption chain in the EU (soy meal–animal feed sector–farms–slaughterhouse–supermarkets) is increasingly vertically integrated. This makes the sector more efficient, but also a more powerful force vis-à-vis producers and consumers. Reduced consumption of animal products is a way to reduce the tension on the soybean chain. Although not desirable from an economic viewpoint, livestock production with a 20–30% N conversion is quite inefficient when converting plant N to proteins in meat, milk or eggs. Price regulations and campaigns to influence consumer behavior represent a possible way to change human diets.
On ‘human waste disposal’, research and dissemination of good hygienic practices can help in keeping excess N out of surface waters.
On a more general note, the instrument of round tables is used with some degree of success, both for oil palm as well as for soybean. The Round Table on Responsible Soy includes private and public sector stakeholders from Brazil and export countries, social and environmental NGOs, and forest inhabitants. In The Netherlands, a ‘Soy Coalition’ has managed to table the adverse social and environmental aspects of soybean production in parliament. Some major food industries have decided to track and trace soy meal using more specific and measurable indicators to ascertain they are selling products grown in a socially and environmentally benign way.
5. Conclusions
8. Of the N excreted by humans into sewerage systems, 79% is emitted into the environment, whereas the complement is removed from wastewater by treatment.
Acknowledgements
This study has greatly benefited from the guidance of Jan Poulisse (FAO, Land and Water Development Division). Edilberto Sena (Frente em Defesa da Amazonia, Santarem), Siemen van Berkum (LEI, The Hague), Jan Maarten Dros (AIDEnvironment, Amsterdam) and several colleagues of the Brazilian National Agricultural Research Organization (Embrapa) in Dourados, Brazil, are all gratefully acknowledged for providing help and feedback while this study was conducted.
References
Aerssens et al., 1998 J.
Aerssens, S. Boonen, G. Lowet and J. Dequeker, Interspecies differences
in bone composition, density, and quality: potential implications for
in vivo bone research. Endocrinology, 139 (1998), pp. 663–670.
Allen et al., 1998 Allen,
R.G., Pereira, L.S., Raes, D., Smith, M., 1998. Crop
evapotranspiration—guidelines for computing crop water requirements.
FAO Irrigation and Drainage Paper 56. Rome..
Alves et al., 1999 B.J.R.
Alves, W.A.R. Lara-Cabezas, E.A. David and S. Urquiaga, Balanço de N em
soja estabelecida em um Latossolo Vermelho Escuro do triângulo mineiro
em condições de plantio direto e preparo convencional do solo Congresso Brasileiro de Ciência do Solo, 1999Brasília, Brazil (1999).
Alves et al., 2005 B.J.R.
Alves, L. Zotarelli, E.S. Araújo, F.M. Fernandes, J.C. Heckler, A.F.A.
Medeiros, R.M. Boddey and S. Urquiaga, Balanço de N em rotação de
culturas sob plantio direto em Dourados. Boletim de Pesquisa e
Desenvolvimento 07, Embrapa Agrobiologia, Seropédica, Brazil (2005).
Antikainen et al., 2005 R. Antikainen, R. Lemola, J.I. Nousiainen, L. Sokka, M. Esala, P. Huhtanen and S. Rekolainen, Stocks and flows of nitrogen and phosphorus in the Finnish food production and consumption system. Agric. Ecosyst. Environ., 107 (2005), pp. 287–305.
Araújo et al., 2006 E.S.
Araújo, C.P. Jantália, R.M. Boddey, S. Urquiaga and B.J.R. Alves,
Importância do N das raízes da soja para a produtividade da cultura
sucessora e para o balanço de N do sistema Circular Técnica 14.,
Embrapa Agrobiologia, Seropédica, Brazil (2006).
Batjes, 2002 Batjes,
N.H., 2002. Soil parameter estimates for the soil types of the world
for use in global and regional modelling (version 2.0). ISRIC Report
2002/02. ISRIC, Wageningen, The Netherlands..
Batjes et al., 2004 Batjes, N.H., Bernoux, M., Cerri, C.E.P., 2004. Soil data derived from SOTER for studies of carbon stocks and change in Brazil (ver. 1.0; GEFSOC Project). Report 2004/03. ISRIC, Wageningen, The Netherlands..
Bleken and Bakken, 1997 M.A. Bleken and L. Bakken, The nitrogen cost of food production: Norwegian society. Ambio, 26 (1997), pp. 134–142.
Bouwman, 1995 Bouwman,
A.F., 1995. Compilation of a global inventory of emissions of nitrous
oxide. PhD Thesis. Agricultural University Wageningen, Netherlands, 143
pp..
Bouwman et al., 2002 A.F.
Bouwman, L.J.M. Boumans and N.H. Batjes, Estimation of global NH3
volatilization loss from synthetic fertilizers and animal manure
applied to arable lands and grasslands. Global Biogeochem. Cycles, 16 (2002), p. 1024.
Bouwman et al., 2005 A.F. Bouwman, G. Van Drecht, J.M. Knoop, A.H.W. Beusen and C.R. Meinardi, Exploring changes in river nitrogen export the world's oceans. Global Biogeochem. Cycles, (2005), p. 19.
Brand and Pfund, 1998 J.
Brand and J.L. Pfund, Site- and watershed-level assessment of nutrient
dynamics under shifting cultivation in eastern Madagascar. Agric. Ecosyst. Environ., 71 (1998), pp. 169–183.
Cerri et al., 1991 C.C.
Cerri, B. Volkoff and F. Andreux, Nature and behavior of organic matter
in soils under natural forest, and after deforestation, burning and
cultivation near Manaus. Forest Ecol. Manage., 38 (1991), pp. 247–257.
Collins et al., 1997 W.J.
Collins, D.S. Stevenson, C.E. Johnson and R.G. Derwent, Tropospheric
ozone in a global-scale three-dimensional Lagrangian model and its
response to NOx emission controls. J. Atmos. Chem., 26 (1997), pp. 223–274.
Da Silva, 2004 A.M. Da Silva, Rainfall erosivity map for Brazil. Catena, 57 (2004), pp. 251–259.
Dentener et al., 2006 F.
Dentener, D. Stevenson, K. Ellingsen, T. van Noije, M. Schultz, M.
Amann, C. Atherton, N. Bell, I. Bey, L. Bouwman, T. Butler, J. Cofala,
B. Collins, F. Drevet, R. Doherty, B. Eickhout, H. Eskes, A. Fiore, M.
Gauss, D. Hauglustaine, L. Horowitz, I. Isaksen, B. Josse, M. Lawrence,
M. Krol, J.F. Lamarque, V. Montanaro, J.F. Müller, V.H. Peuch, G.
Pitari, J. Pyle, S. Rast, J. Rodriguez, M. Sanderson, N. Savage, D.
Shindell, S. Szopa, K. Sudo, O. Wild and G. Zeng, The global
atmospheric environment for the next generation. Environ. Sci. Technol., 40 (2006), pp. 3586–3594.
De Willigen, 2000 De
Willigen, P., 2000. An analysis of the calculation of leaching and
denitrification losses as practised in the NUTMON approach. Report 18. Plant Research International, Wageningen, The Netherlands..
FAO, 1977 FAO, China: recycling of organic wastes in agriculture, FAO Soils Bulletin 40, Food and Agriculture Organization of the United Nations, Rome (1977), p. 107.
FAO, 1998 FAO,
ISRIC, UNEP, CIP, 1998. Soil and terrain digital database for Latin
America and the Caribbean at 1:5 million scale. Land and Water Digital
Media Series No. 5. FAO, Rome..
FAO, 2004a FAO, Scaling soil nutrient balances—enabling mesoscale approaches for African realities, FAO Fertilizer and Plant Nutrition Bulletin 15, FAO, Rome (2004).
FAO, 2004b FAO, Fertilizer Use by Crop in Brazil, FAO, Rome (2004).
FAO, 2006 FAO, Livestock's long shadow, Environmental Issues and Options, FAO, Rome (2006).
FAO, 2007 FAO 2007. The statistics division. Livestock Statistics. Data available from http://www.fao.org/es/ess/rmlive.asp. FAO, Rome..
Freitas et al., 2000 P.L.
Freitas, P. Blancaneaux, E. Gavinelli, M.C. Larre-Larrouy and C.
Feller, Nível e natureza do estoque orgânico de latossolos sob
diferentes sistemas de uso e manejo. Pesqui. Agropecu. Bras., 35 (2000), pp. 157–170.
Galloway and Cowling, 2002 J.N. Galloway and E.B. Cowling, Reactive nitrogen and the world: 200 years of change. Ambio, 31 (2002), pp. 64–71.
Galloway et al., 2004 J.N.
Galloway, F.J. Dentener, D.G. Capone, E.W. Boyer, R.W. Howarth, S.P.
Seitzinger, G.P. Asner, C.C. Cleveland, P.A. Green, E.A. Holland, D.M.
Karl, A.F. Michaels, J.H. Porter, A.R. Townsend and C.J. Vörösmarty, Nitrogen cycles: past, present, and future. Biogeochemistry, 70 (2004), pp. 153–226.
Galloway et al., 2007 J.N.
Galloway, M. Burke, E. Bradford, R. Naylor, W. Falcon, A.K. Chapagain,
J.C. Gaskell, E. McCullough, H.A. Mooney, K.L.L. Oleson, H. Steinfeld,
T. Wassenaar and V. Smil, International trade in meat: the tip of the
pork chop. Ambio, 36 8 (2007), pp. 622–629.
Graca et al., 1999 P.M.L.
Graca, P.M. Fearnside and C.C. Cerri, Burning of Amazonian forest in
Ariquemes, Rondonia, Brazil: biomass, charcoal formationa and burning
efficiency. Forest Ecol. Manage., 120 (1999), pp. 179–191.
Hungria et al., 2006 M. Hungria, J.C. Franchini, R.J. Campo, C.C. Crispino, J.Z. Moraes, R.N.R. Sibaldelli, I.C. Mends and J. Arihara, Nitrogen nutrition of soybean in Brazil: contributions of biological N2 fixation and N fertilizer to grain yield. Can. J. Plant Sci., 86 (2006), pp. 927–939.
IBGE, 2006 IBGE, Brazilian Institute for Geography and Statistics (IBGE) 2006. Agricultural Municipalities Production Data..
Jantalia et al., 2006 C.P.
Jantalia, B.J.R. Alves, L. Zotarelli, R.N. Boddey and S. Urquiaga,
Mudanças no estoque de C do solo em areas de produção de grãos:
avaliação do impacto do manejo do solo, B.J.R. Alves, Editor, et al.Manejo de sistemas agrícolas: Impacto no seqüestro de C e nas emissões de gases de efeito estufa, Gênesis, Porto Alegre (2006), pp. 35–58.
Kanmegne, 2004 Kanmegne,
J., 2004. Slash and burn agriculture in the humid forest zone of
Southern Cameroon: soil quality dynamics, improved fallow management
and farmers’ perceptions. Ph.D. Thesis, Wageningen University, The
Netherlands..
Kantor et al., 1997 Kantor,
S.L., Lipton, K., Manchester, A., Oliveira, V., 1997. Estimating and
addressing America's Food Losses. USDA Food Review January–April 1997..
Keyzer et al., 2005 M.A.
Keyzer, M.D. Merbis, I.F.P.W. Pavel and C.F.A. Van Wesenbeeck, Diet
shifts towards meat and the effects on cereal use: can we feed the
animals in 2030?. Ecol. Econom., 55 (2005), pp. 187–202.
Krusche et al., 2003 A.V.
Krusche, P.B. De Camargo, C.E. Cerri, M.V. Ballester, L.B.L.S. Lara,
R.L. Victoria and L.A. Martinelli, Acid rain and nitrogen deposition in a sub-tropical watershed (Piracicaba): ecosystem consequences. Environ. Pollut., 121 (2003), pp. 389–399.
Lesschen et al., 2007 J.P.
Lesschen, J.J. Stoorvogel, E.M.A. Smaling, G.B.M. Heuvelink and A.
Veldkamp, A spatially explicit methodology to quantify soil nutrient
balances and their uncertainties at the national level. Nutrient Cycling Agroecosyst., 78 (2007), pp. 111–131.
Llilienfein and Wielcke, 2003 J.
Llilienfein and W. Wielcke, Element storage in native, agri-, and
silvicultural ecosystems of the Brazilian savanna. I. Biomass, carbon, nitrogen, phosphorus and sulphur. Plant Soil, 254 (2003), pp. 425–442.
Lovato, 2001 Lovato,
T., 2001. Dinâmica do carbono e do nitrogênio do solo afetada por
preparos do solo, sistemas de cultura e adubos nitrogenados. Porto
Alegre, Universidade Federal do Rio Grande do Sul. 130 p. (Tese de
Doutorado)..
Matson et al., 1987 P.A. Matson, P.M. Vitousek, J.J. Ewel, M.J. Mazzarino and G.P. Robertson, Nitrogen transformations following tropical forest felling and burning on a volcanic soil. Ecology, 68 (1987), pp. 491–502.
Mielniczuk et al., 2003 J.
Mielniczuk, C. Bayer, F.M. Vezzani, T. Lovato, F.F. Fernandes and L.
Debarba, Manejo de solo e culturas e sua relação com os estoques de
carbono e nitrogênio do solo. Tópicos em Ciência do Solo, 3 (2003), pp. 209–248.
Mitchell and Jones, 2005 T.D.
Mitchell and P.D. Jones, An improved method of constructing a database
of monthly climate observations and associated high-resolution grids. Int. J. Climatol., 25 (2005), pp. 693–712.
Morgan, 1995 R.P.C. Morgan, Soil Erosion & Conservation, (2nd ed.), Longman Group Limited, Essex, UK (1995).
Neill and Davidson, 1999 C. Neill and E.A. Davidson, Soil carbon accumulation for loss following deforestation for pasture in the Brazilian Amazon, R. Lal, Editor, et al.Global Climate Change and Tropical Ecosystems, CRC Press, Boca Raton, FL (1999), pp. 197–211.
Paustian et al., 1992 K. Paustian, W.J. Parton and J. Person, Modeling soil organic matter in organic-amended and nitrogen-fertilized long-term plots. Soil Sci. Soc. Am. J., 56 (1992), pp. 476–488.
Renard et al., 1997 K.G.
Renard, G.R. Foster, G.A. Weesies, D.K. McCool and D.C. Yoder,
Predicting soil erosion by water: a guide to conservation planning with
the revised Universal Soil Loss Equation (RUSLE). USDA Agric. Handbook, 703 (1997), p. 384.
Resende et al., 1997 M.
Resende, N. Curi, S.B. Rezende and G.F. Correa, Pedologia: base para
distinção de ambientes, (2ed edn.), NEPUT, Viçosa (1997), p. 367.
Roscoe, 2005 Roscoe,
R., 2005. Dinâmica da matéria orgânica em solos de cerrado. In: XXX
Congresso Brasileiro de Ciência do Solo, 2005, Recife. Solos,
sustentabilidade e qualidade ambiental. Recife: Embrapa Solos;
Sociedade Brasileira de Ciência do Solo; Universidade Federal de
Pernambuco. 1 CD-ROM..
Roscoe and Buurman, 2003 R. Roscoe and P. Buurman, Tillage effects on soil organic matter in density fractions of a Cerrado Oxisol. Soil Till. Res., 70 (2003), pp. 107–119.
Schlesinger, 1986 W.H. Schlesinger, Changes in soil carbon storage and associated properties with disturbance and recovery, J.R. Trabalka, D.E. Reichle, Editors , The Changing Carbon Cycle. A Global Analysis, Springer Verlag, New York, Heidelberg (1986), pp. 194–220.
Silva et al., 1994 J.E.
Silva, J. Lemainski and D.V.S. Resck, Perdas de matéria orgânica e suas
relações com a capacidade de troca catiônica em solos da região de
cerrados do oeste baiano. Revista Brasileira de Ciência do Solo, 18 (1994), pp. 541–547.
Smil, 1999 V. Smil, Nitrogen in crop production. An account of global flows. Global Biogeochem. Cycles, 13 (1999), pp. 647–662.
Smil, 2001 V. Smil, Enriching the Earth, MIT Press, Cambridge, MA (2001).
Smil, 2002 V. Smil, Nitrogen and food production: proteins for human diets. Ambio, 31 (2002), pp. 126–131.
Stoorvogel et al., 1993 J.J.
Stoorvogel, E.M.A. Smaling and B.H. Janssen, Calculating soil nutrient
balances in Africa at different scales. 1. Supra-national scale. Fertilizer Res., 35 (1993), pp. 227–235.
Tognon et al., 1998 A.A.
Tognon, J.L.I. Dematte and J.A.M. Dematte, Teor e distribuição da
matéria orgânica em Latossolos das regiões da Floresta Amazônica e dos
Cerrados do Brasil Central. Sci. Agric. (Piracicaba), 55 (1998), pp. 343–354.
Trebs et al., 2005 I.
Trebs, L.L. Lara, L.M.M. Zeri, L.V. Gatti, P. Artaxo, R. Dlugi, J.
Slanina, M.O. Andreae and F.X. Meixner, Dry and wet deposition of
inorganic nitrogen compounds to a tropical pasture site (Rondonia Brazil). Atmos. Chem. Phys. Discuss., 5 (2005), pp. 3131–3189.
USGS, 1996 US Geological Survey (USGS), 1996. GTOPO30—global digital elevation model. Sioux Falls, US..
Van Berkum et al., 2006 Van
Berkum, S., Roza, P., Pronk, B., 2006. Sojahandel- en ketenrelaties.
Sojaketens in Brazilië, Argentinië en Nederland. Rapport 5.06.08
Agricultural Economics Research Institute (LEI), The Hague, The
Netherlands (in Dutch)..
Van der Hoek, 1998 K.W. Van der Hoek, Nitrogen
efficiency in global animal production, K.W. Van der Hoek, J.W.
Erisman, S. Smeulders, J.R. Wisniewski, J. Wisniewski, Editors , Nitrogen, the Confer-N-s, Elsevier, Amsterdam (1998), pp. 127–132.
Van Gelder and Dros, 2005 Van Gelder, J.W., Dros,
J.M., 2005. Van oerwoud tot kippenbout. Effecten van sojateelt voor
veevoer op mens en natuur in het Amazonegebied—een ketenstudie (in
Dutch). Profundo, AIDEnvironment, Milieudefensie, Cordaid, The
Netherlands..
Wischmeier and Smith, 1978 W.H.
Wischmeier and D.D. Smith, Predicting Rainfall Erosion Losses, USDA
Agricultural Research Service Handbook 537, USDA, Washington, DC
(1978).
Zhu and Chen, 2002 Z.L. Zhu and D.L. Chen, Nitrogen fertilizer use in China—contributions to food production, impacts on the environment and best management strategies. Nutrient Cycling Agroecosyst., 63 (2002), pp. 117–127.
Zotarelli et al., 1998 Zotarelli,
L., Resende, A., Torres, E., Hungria, M., Urquiaga, S., Boddey, R.M.,
Alves, B.J.R., 1998. O papel da soja para as reservas de N de um
Latossolo Roxo sob plantio direto na região de Londrina, PR. In:
Reunião Brasileira de Fertilidade do Solo e Nutrição de Plantas, 23,
Lavras..
![]() |
![]() |