Modeling framework
We used an open source, fully documented, and publicly available medium run applied general equilibrium (AGE) model28 with explicit treatment of subnational land markets divided in Agroecological Zones (AEZ), nicknamed GTAP-AEZ22. The GTAP-AEZ framework is based on decision nests at which agricultural producers decide on land cover conversions (Supplementary Information Fig. S2), for example, from pastures to cropland, and then on the allocation of individual crops within the cropland. As producers in different AEZs are connected through land, labor, and capital markets, competition among land uses, and supply chains, the GTAP-AEZ model is ideally suited to study within-country changes in land use across AEZs. Moreover, through an explicit treatment of international trade flows, the GTAP-AEZ framework allows for tracking the effects of regional policies on land use patterns in other countries. The GTAP modeling framework has played an important role in understanding market-mediated effects of important policies including biofuels, soybean, and oilseed markets, and the effects of productivity gains during the Green Revolution25,27,29,30.
We updated the standard GTAP-AEZ model to include a nesting structure that separates the decision to convert forest to agricultural land from the decision to convert pasture to cropland, which is justified by the observation that deforested lands transition first into pastures, and then onto cropland31. This nesting structure applies to all the regions. We also adopted regional elasticities of transformation, from natural covers to agricultural land, and between pastures and cropland, calibrated based on recent historical changes31. We further updated the income elasticities of demand for agricultural and food products to reflect the latest work in this area32. A critical assumption underlying the GTAP-AEZ framework is the productivity of marginal, hitherto, uncultivated lands, as it determines the extensive margin of land expansion. Another key assumption in the GTAP-AEZ model is the response of yields to changes in commodity and input prices33. For both we use the assumptions in the original GTAP-AEZ model22. Given the uncertainty regarding these parameters, we conduct extensive systematic sensitivity analysis of our results to alternative parametric configurations (Supplementary Information S2).
Underlying the model there is a database that consistently represents production, consumption, and trade patterns of 140 regions and 57 sectors in year 201134. To make solution times and model output manageable, we aggregated the model into 11 regions: Brazil, Bolivia, Argentina, Paraguay, Rest of Latin America, US-Canada (North America), European Union (28 countries), China, Malaysia and Indonesia, sub-Saharan Africa, and the rest of the world. We also collapsed the 57 commodity sectors into 18 sectors (i.e., paddy rice, wheat, coarse grains, oilseeds, raw sugar, grazing livestock, non-grazing livestock, forestry, extractive industries, processed livestock, vegetable oils, processed rice, processed sugar, other processed food, chemicals, manufactures and services). For Brazil, we considered the GTAP aggregate oilseed commodity as soybeans because soybeans account for more than 96% of oilseed production in Brazil35. The database is complemented with data on agricultural land rents by land use and natural land covers at the level of Agroecological Zones (AEZ), also representative of 201136.
Spatial Footprint Scenarios (SFS) and Market Share Thresholds of Zero Deforestation Corporate Policies (ZDCPs)
The SFS (Fig. 1) are designed to assess how much deforestation would be avoided by implementing different configurations of the company- and importing country-led ZDCPs. Except for the ASM, most ZDCPs–either voluntary or imposed—have not been implemented6,8,11,17. Therefore, the SFS are counterfactual, non-observed states of the world. We estimate changes in deforestation and other economic outcomes as the difference between the counterfactual SFS and a baseline (as explained below). The baseline includes patterns of land use, land cover, and other economic outcomes obtained by letting the model simulate the changes in equilibrium as the economy responds to a set of drivers of land use during the period 2011-2016, without any land restriction to land expansion in Brazil. The baseline is therefore a plausible representation of the world in the absence of any type of ZDCPs.
We chose 2011 as our initial year because this is the most recent year for which we have the consistent snapshot of the world economy, including land use and land cover—GTAP Database and AEZ Database, V934,36 —needed to calibrate the GTAP-AEZ model. We stop in 2016 as we seek to understand a “medium run” time horizon whereby global production, consumption and trade have enough time to adjust to an equilibrium that includes all the economic effects of the modeled ZDCPs. The ASM was in place during our period of study, and it should be considered a baseline relative to further hypothetical policy developments analyzed here. The drawback of including the ASM in the baseline is that we would not be able to report the ASM outcomes, which, by virtue of its pioneering status, is a natural benchmark of future ZDCPs in Brazil’s soybean sector. For this reason, we exclude the ASM from the baseline by not imposing land restrictions in the Amazon.
The SFS we evaluate are as follows:
-
Amazon Soy Moratorium. This scenario uses the spatial footprint in the Brazilian Amazon of the companies that implemented the moratorium in 2006. These companies are: Abc Industria, ADM, Amaggi, Bunge, Cargill, Louis Dreyfus, Seara, Fiagril, Nidera, Noble, Cofco, Baldo, Imcopa, Agrex, CHS, Coamo, Engelhart CTP, Gavilon, Glencore, Invivo, Marubeni, Multigrain, Nova Agri, Olam, Perdue, Sodrugestvo, Timbro, and Selecta. Other companies that are part of the ASM do not export soybeans from the Amazon are: Binatural, JBS, Oleos menu, Agribrasil, and Culturale. The duration of the ASM has been extended indefinitely14.
-
Global voluntary ZD Commitments + Soy Moratorium. This scenario includes the ASM companies above and adds all Global ZDCPs as if they were implemented in 2011. The Global ZDCPs have been adopted by a subset of the companies that agreed on the ASM. These are (pledge year and in parentheses): ADM (2015 company pledges); Amaggi (2017 company pledges), Bunge (2015 company pledges) Cargill (2014 New York Declaration on Forests); Louis Dreyfus (2018 company pledges); Cofco (2019 Soft Commodities Forum); Glencore (2019 Soft Commodities Forum); and Denofa do Brazil (2014 New York Deforestation of Forests).
-
Import restrictions imposed by the European Union (EU). Agriculture-driven deforestation has become an increasingly polarizing issue between the EU and Brazil, and is a central issue in a potential trade agreement between the EU and the MERCOSUR, a trade bloc agreement among Argentina, Brazil, Paraguay, and Uruguay37. We therefore also explore the effects of the adoption by the EU of mandatory rules currently considered by the European Parliament that would de facto impose ZDCPs on the companies sourcing soybeans from Brazil8. We consider 155 traders exporting to EU plus Switzerland and the United Kingdom that would only procure their soybeans from areas already converted to agriculture prior to 2011. The EU countries appearing as importers consist of Belgium, Bulgaria, Croatia, Cyprus, Denmark, Finland, France, Germany, Greece, Hungary, Ireland, Italy, Latvia, Lithuania, Netherlands, Poland, Portugal, Romania, Slovakia, Slovenia, Spain, and Sweden.
-
Hypothetical import restrictions imposed by China, assumed to be similar to those being currently evaluated by the EU 17 . We consider 140 traders that export to China including Hong Kong (in addition to the ASM) that would only procure their soybeans from areas already converted to agriculture prior to 2011.
-
Hypothetical import restrictions imposed by both China and the EU. These simulations provide an upper bound estimate of the ZCDPs.
In the SFS 3-5 (scenarios 3-5 in the main text) we assume that if a company supplying either EU or Chinese markets decides to produce deforestation-free soybean to preserve market share in one destination, the company will apply those restrictions to their entire supply chain. In other words, we do not allow for different supply chains from the same trader when exporting to different destinations. This is a realistic assumption as companies consider supply chain differentiation extremely costly due to the unwieldy procedures that would be needed to monitor, trace, and certify production21,38.
Market Share Thresholds
Uncertainties exist regarding the critical market share (i.e., the percentage of total regional market share held by corporations with ZDCPs) needed to discourage farmers from selling soybeans produced in recently cleared land to non-committed traders6. If no soybean buyers within a region have a ZDC, farmers have no incentives to avoid forest for soybean clearing. Alternatively, if only committed traders buy soybeans within a region, producers should be forced to comply with the ZDCP land use restrictions to sell their soybeans. In many Brazilian regions, traders with and without ZDCPs purchase soybeans, thus producers with soybeans that are not deforestation free can typically sell their products. We posit that with increasing regional ZDCP market share, the difficultly of selling non-compliant soybeans increases. At some critical market share threshold, farmers may be completely disincentivized from producing soybeans on non-compliant lands that were recently deforested due to the difficulties in selling their product 39,40.
We bound the uncertainty about the competition structure needed to ensure compliance through three market share thresholds built using the export market shares of all the companies active in the Brazil’s soybean market from the Trase v2.4 database 41:
-
The most restrictive market share thresholds requires that at least 75% of the soybeans exported from a given municipality are bought by companies with voluntary ZD commitments. In this scenario, 10% of the area under soybeans in Brazil is subject to the ASM (Fig. 1f). By adding pledged global voluntary ZD commitments in other biomes to the ASM, 27% of Brazil’s soybean area would be under agreements to halt forest conversion for soy production (Fig. 1g).
-
A less conservative market share thresholds requires an export threshold of 50%. Under this scenario, the area under soybeans that is affected by global voluntary ZD commitments under the current pledges amount to 48% of Brazil’s total soybean area (Fig. 1g).
-
The least restrictive scenario requires at least one committed company to be present in the municipality (>0% of market share covered by voluntary ZD commitments). Under global voluntary ZD commitments, 75% of Brazil’s soybean area would be subject to ZDCPs (Fig. 1g).
Land Cover Definitions
We use two different definitions of forests in Brazil to accommodate different biome characteristics and ZDCP targets 42. Definition A was exclusively based on mapped forest cover. Definition B included natural grasslands outside the Amazon Biome, which may have high conservation value and are included in some traders ZDCP definition [e.g., “Transforms our supply chain to be deforestation free while protecting native vegetation beyond forests.” 43]
-
Forest definition “A”: Forest is defined as forest only, as mapped by PRODES for the Amazon 44 and by Mapbiomas for other biomes [Mapbiomas v4 45, classes 1, 2, 3, 4, and 5]. Forest area was derived excluding forest regrowth, with forest base year of 2006 45. We used PRODES for the Amazon biome, because PRODES deforestation maps define the baseline for ASM monitoring, implementation, and enforcement. We used Mapbiomas outside the Amazon Biome. To our knowledge, Mapbiomas provides the most accurate and consistent large scale land use and land cover classification for Brazil.
-
Forest definition “B”: Forest is defined as forest only in the Amazon biome, as mapped by PRODES, and forest and grasslands in all other Brazilian biomes, as mapped by Mapbiomas.
Land Use and Land Cover Databases: In addition to the data on soybean export market share and forests, we gathered municipality-level data on agricultural land cover from Mapbiomas v4 45: total cropland [classes 18-20] and pasture area [class 15], soybean area (ha), soybean production (tonnes), areas with both maize and soybean (ha), maize second harvest area (ha), calculated as the area of second harvest maize that is greater or equal to the area of soy harvested 46, and cattle headcount (heads)47. These data were used to build the different versions of the GTAP-AEZ database, as explained below.
The supply-side spatial footprint scenarios, market share thresholds, and forest definitions give rise to thirty different databases with land use and land cover in each Brazilian municipality (five SFS * three market share thresholds* two Forest Definitions = thirty databases.) We use these land use/land cover databases to build biome-specific distributions of land cover and soybean production with and without ZDCPs that can be used to calibrate the counterfactual experiments using the GTAP-AEZ model.
Model calibration
We use the databases discussed in Supplementary Information S1 to recalibrate the GTAP-AEZ model so that Brazil is split into biomes instead of the standard AEZs. This requires rebuilding the original GTAP and GTAP-AEZ databases. The algorithm to split Brazil’s agricultural output values into biomes proceeds as follows. For each ZDC spatial configuration, market share threshold, and forest definition, we use the following algorithm
-
Overlay the AEZ map used in the GTAP-AEZ database36 on a municipality-level map of Brazil48. In case that a municipality is split across more than one AEZ, assign the municipality to the AEZ with the largest intersection.
-
Overlay a biome map over the AEZ and municipality maps for Brazil. The biomes generally encompass several AEZs and the same AEZ can occur in different biomes. Biomes other than the Amazon and Cerrado are in an “Other” category. In case that a municipality is split across more than one biome, we assign the municipality to the biome with the higher ZDC implementation stringency, prioritizing the Amazon, second the Cerrado, and all other.
-
Each municipality receives a unique id for each biome-AEZ combination, for example: AEZ5 becomes AEZ5-Amazon, AEZ5-Cerrado, and AEZ5-Other.
-
Compute compliant market share thresholds for each municipality (using pledges as of 2017), and then categorize the biome-AEZ ids into compliant and non-compliant based on the SFS. For example: AEZ5 becomes AEZ5-Amazon-ZDC, AEZ5-Cerrado-ZDC, and AEZ5-Other-ZDC, and AEZ5-Amazon-Non-ZDC, AEZ5-Cerrado- Non-ZDC, and AEZ5-Other-Non-ZDC.
-
Use aggregate municipality-level land cover (cropland, pasture, and forest area) and land use (soybean area, soy production, areas with both maize and soy, cattle headcount) to assign land cover and land use areas to each biome-AEZ-market share thresholds level.
-
For all regions other than Brazil, build a conventional GTAP-AEZ database representative of 2011. This step uses a database of land use and land cover areas at the level of AEZs 36,49 to split the country-level output value of relevant products (crops, grazing livestock, and forestry) in the standard GTAP database 34 into AEZs.
-
For Brazil, we use the AEZ-Biome area and production shares created in steps 1-4 to split the aggregate output values of oilseeds, coarse grains, grazing livestock, and forestry into biomes. Each new database represents a counterfactual year 2011 in which some of the area in each AEZ-BIOME was under a ZDCP commitments pledged before 2020. The simulations answer the question: how different area, production, and consumption would have been in 2016 if the pledged ZDC had been in place since 2011.
-
The area of the crops other than oilseeds and coarse grains (paddy rice, etc.) are shared out in each biome-AEZ in proportion to the cropland.
In each experiment we halt land conversion between forest and agriculture in the areas assumed under ZDCs by way of a subsidy that compensates producers for the economic losses of not transforming forests on to agriculture.
Deforestation Leakage Rates
Following the literature on carbon leakage50, we define the deforestation leakage rate as the increase in deforestation in regions without restrictions (No ZDCP) induced by the measures taken in regions with ZD policies (ZDCP) as a percentage share of the absolute value of deforestation in regions with ZDCP. Formally
$$Deforestation Leakage Rate=\frac{{\Delta }Deforestatio{n}_{No ZDCP}}{\left|{\Delta }Deforestatio{n}_{ZDCP}\right|}\times 100.$$
Where \({\Delta }\) denotes the difference between deforestation outcomes in the baseline and counterfactual scenarios, and
$${\Delta }Deforestatio{n}_{Global}={\Delta }Deforestatio{n}_{ZDCP}+{\Delta }Deforestatio{n}_{No ZDCP}.$$
Greenhouse Gas Emissions
Greenhouse gas emissions (GHGs) from the changes in land cover associated with the different experiments are calculated using the also open-source AEZ Emission Factor (AEZ-EF) Model51. The AEZ-EF model closely follows IPCC GHG inventory methods and relies on its default values. The model includes cover-specific (cropland, pastures, and forests) subnational carbon estimates for biomass (above and below-ground), dead organic matter, and soil carbon52. It also includes data on carbon remaining on harvested wood products, non-CO2 emissions, and foregone sequestration. The carbon stock data is combined with assumptions about carbon sequestration from forest growth (foregone if converted), mode of conversion, and CO2 emissions from land clearing using fire, and the fraction of carbon that remains sequestered in wood products during a 30-year time horizon. The AEZ-EF model is designed to estimate land use emissions from land use transitions predicted by comparative static economic models, whereby one starts with a baseline and estimates the resulting final equilibrium. The AEZ-EF model underlies the emission estimates in several analysis of the indirect land use effects of biofuels emissions and land conservation measures27,53.