Socioeconomic and environmental effects of soybean production in metacoupled systems | Scientific Reports –

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 11, Article number: 18662 (2021)
Metrics details
Human–environment interactions within and across borders are now more influential than ever, posing unprecedented sustainability challenges. The framework of metacoupling (interactions within and across adjacent and distant coupled human–environment systems) provides a useful tool to evaluate them at diverse temporal and spatial scales. While most metacoupling studies have so far addressed the impacts of distant interactions (telecouplings), few have addressed the complementary and interdependent effects of the interactions within coupled systems (intracouplings) and between adjacent systems (pericouplings). Using the production and trade of a major commodity (soybean) as a demonstration, this paper empirically evaluates the complex effects on deforestation and economic growth across a globally important soybean producing region (Mato Grosso in Brazil). Although this region is influenced by a strong telecoupling process (i.e., soybean trade with national and international markets), intracouplings pose significant effects on deforestation and economic growth within focal municipalities. Furthermore, it generates pericoupling effects (e.g., deforestation) on adjacent municipalities, which precede economic benefits on adjacent systems, and may occur during and after the soybean production takes place. These results show that while economic benefits of the production of agricultural commodities for global markets tend to be localized, their environmental costs tend to be spatially widespread. As deforestation also occurred in adjacent areas beyond focal areas with economic development, this study has significant implications for sustainability in an increasingly metacoupled world.
Almost all places around the world are currently interconnected through the growing flows of people, goods, capital, and information. While such interconnections have existed for a long time, they have increased drastically in recent decades1,2,3. For example, while over the past three decades the world’s total food production has doubled4, a growing fraction of this expansion is associated with commodities produced for global markets6. Thus, current dynamics of land-use and land-cover change are being shaped by new forces, such as the ever-increasing demand for globalized commodities5, with the result that globalization is now more influential than ever before1. This poses unprecedented challenges for the sustainable management of natural resources around the world. The international soybean trade offers an excellent lens through which to explore such challenges.
Soybean has a long history of cultivation for human food and animal feed in China, but since the 1970s, it has been transformed into a commodity crop with industrial uses, particularly in the livestock sector7. This sector has grown considerably in response to a dietary shift towards more animal protein due to the increase in economic affluence in China over the last three decades8. As a result of this dietary shift, starting in 2002 soybean (and its derivatives) started to be considered an industrial, rather than an agricultural commodity in China, which allowed the liberalization of their importation9. The result is that soybean imports, predominantly from the US and Brazil, increased almost exponentially7.
The drastic increase in soybean production and export has brought not only socioeconomic effects10,11, but also serious environmental consequences, as it has both directly and indirectly induced land-use and land-cover changes in Brazilian biomes12,13,14,15. Thus, the international trade of soybean constitutes a complex human–environment interaction that can be conceptualized as a metacoupled system16. Under this conceptual framework, distant, adjacent, and local interactions may occur simultaneously and complementarily17,18. While the metacoupling framework has been used to evaluate the environmental outcomes induced by distant and adjacent interactions such as trade and has been successfully applied to many systems worldwide18,19,20,21,22,23,24,25, no empirical studies have been conducted to assess the economic and environmental effects of intracouplings (human–environment interactions such as soybean production within a system) and pericouplings (human–environment interactions among adjacent systems16) within the context of a dominant telecoupling process (Fig. 1). To fill this important knowledge gap, this study focused on all the municipalities within the state of Mato Grosso, Brazil (a major producer of soybean for international and national markets26), and evaluated socioeconomic and environmental outcomes of the soybean production at the producing municipalities and on neighboring municipalities.
Conceptualization of the telecoupled system, with soybean trade as the major flow. Telecoupling is conceptualized as the flow of soybean between Mato Grosso, Brazil and major international markets (China and the European Union), as well as domestic markets (other states) within Brazil. In 2018 ca. 56% of the state’s total soybean produced went to China and the European Union. In exchange, Mato Grosso received flows of capital and information. Maps were created by the Authors using the software QGIS 3.16.3 (
Our rationale reflects the behavior of soybean producers in Brazilian agricultural areas27,28, in which land-use decisions are driven by producers’ expectations of economic benefits based on previous years of commodity price and market demand, which in turn stimulates further land conversion28. In this regard, our study goes beyond the municipality boundary of the focal production system27 to assess how the aggregated pattern of individual producers in one municipality drives intracoupling and pericoupling effects. In regions such as Mato Grosso state in Brazil, soybean production in a given year requires preparation during previous years27,28. Furthermore, once soybean production gets started in an area, other dynamics also occur, which may lead to different socioeconomic and environmental outcomes29. Therefore, we hypothesize that while soybean production triggers land change processes (e.g., deforestation) not only within but also in neighboring municipalities, associated economic benefits only accrue within the soybean producing municipality. In this study we refer to “intracoupling effects” as those occurring within municipalities while “pericoupling effects” as those occurring in adjacent municipalities.
Previous studies have addressed indirect land use change by taking pixel-based approaches to demonstrate the patterns in which soybean replace pastureland, while leading to the expansion of pastureland in forested areas30,31; or how soybean production promotes economic benefits at varying geographical levels29,32. However, this study addresses two crucial components of complex metacoupled systems. The first one refers to the interactive relationship between land-use and social process within intracoupled systems, while the second one refers to the pericoupling effects between adjacent systems, both at the municipality level. Here we understand an intracoupled system as a place (e.g., a municipality) that exhibits reciprocal human-nature interactions and feedbacks, while pericoupled systems as adjacent systems connected through flows of people, information, resources, and capital, among others16. Such interactions include the use of natural resources to meet human needs (e.g., soybean production) and its socio-environmental effects. Considering the global flows of various natural resources, intracoupled and pericoupled systems may interact with telecoupled systems1, which together form metacoupled systems—i.e., an assemblage of coupled systems at different distances and spatial configurations16. Understanding the interactions among these coupled systems and their multiple effects is crucial for guiding informed decision-making processes related to socioeconomic development in agricultural frontiers.
Mato Grosso exhibited a significant agricultural expansion, although it was not spatially uniform throughout the state (Fig. 2a). For instance, while municipalities with the largest soybean production in 2000 also exhibited the largest production in 2018, the municipalities that exhibited a significant increase in their soybean production between 2000 and 2018 tended to be clustered around the largest producers (Fig. 2a and Supplementary Fig. 1). This pattern is crucial for understanding the influence of neighboring municipalities when examining the spatiotemporal effects of intracouplings and pericouplings on environmental and socioeconomic outcomes, such as deforestation and the change in Gross Domestic Product (GDP) per capita, respectively (Fig. 2a,b).
Dynamics of deforestation, Gross Domestic Product (GDP) per capita, and soybean production in Mato Grosso, Brazil, between 2000 and 2018. (a) Maps represent the changes observed between 2000 and 2018 in forest cover (i.e., deforestation), GPD per capita, and soybean production in each of the 141 municipalities comprising the state of Mato Grosso, Brazil. (b) Line graphs show the trends (2000–2018) of these variables aggregated to the state level. Maps were created by the Authors using the software QGIS 3.16.3 (
We developed two models to assess the effects of soybean production on environmental (deforestation model) and socioeconomic [Gross Domestic Product (GDP) per capita model] variables within (intracoupling) and in adjacent (pericoupling) municipalities (see “Methods“). Results of the deforestation model show significant intracoupling and pericoupling effects of the changes in soybean production on deforestation (Fig. 3a) at five-year intervals (i.e., moving window, see “Methods“), after controlling for other biophysical and socioeconomic variables [e.g., slope, forest cover, GDP per capita, area under pastureland, and second-crop production (Supplementary information 1)]. These results suggest that throughout the 2000–2018 period, soybean production was significantly and positively correlated with deforestation not only within municipalities (i.e., intracoupling effect) but also across neighboring municipalities (i.e., pericoupling effect), with few exceptions (Fig. 3a and Supplementary information 1). The deforestation model also revealed that the increase in pastureland areas at five-year intervals exhibited a significant positive relationship with deforestation but only within the same municipality and not in adjacent municipalities (Supplementary information 1). Although the increase in pastureland exhibited a significant negative relationship with deforestation in adjacent municipalities during some time periods (e.g., 2010–2014), the intracoupling effect contributed with around 93% of the total significant relationship between the increase in pastureland and deforestation. The observed significant negative pericoupling effect of pasture area increase in few deforestation model subsets (Supplementary information 1) may be explained by the expansion of pastureland in previous years. Hence, the large availability of pastureland during a given period in a set of municipalities discourages the increase in their neighbors, thus negatively affecting the deforestation trend.
Results of spatial Durbin Error Models (SDEM) for fifteen subsets of model simulations between 2000 and 2018 using a moving window of five years. (a) SDEM to account for intracoupling and pericoupling effects of soybean production on deforestation across municipalities in the state of Mato Grosso, Brazil. (b) SDEM to account for intracoupling and pericoupling effects of soybean production on Gross Domestic Product (GDP) per capita across municipalities in the state of Mato Grosso, Brazil.
After controlling for variables such as slope, deforestation, pastureland, agricultural labor, maize second-crop production, cattle herd, and forest cover (Supplementary information 2), our results of the GDP per capita model show that soybean production was significantly and positively related with the growth of GDP per capita throughout most of the 2000–2018 period, but only within municipalities (i.e., intracoupling effect) (Fig. 3b and Supplementary information 2). These results suggest that soybean production within a municipality was associated with the economic growth of the same municipality, while showing no association (or a weak association) with the economic conditions of neighboring municipalities (i.e., non-pericoupling effect—only 15% of the total positive effects of soybean production associated with an increase on GDP per capita were from soybean production in neighboring municipalities—Fig. 3b). Furthermore, a Mann–Whitney U test revealed that compared to non-soybean producing municipalities, the GDP per capita in municipalities producing soybean was higher throughout the 2000–2018 (Supplementary Fig. 2). Our model results also revealed that during the five-year intervals in which soybean production did not exhibit a significant positive relationship with the growth of GDP per capita within municipalities, the variable ‘maize second-crop production’ did (Supplementary information 2). This suggests that an increase in the production of maize as a second-crop significantly complements the economic growth within a municipality.
Considering the results for Global Moran’s I tests for spatial dependence on the GDP per capita models, which indicated non-spatial dependence (except for years 2005, 2006, and 2007—i.e., p-values < 0.05), ordinary least square (OLS; Supplementary information 3) regression were applied and results confirmed the significance (p-values < 0.05) of soybean production growth as a major variable positively correlated with the local GPD per capita growth at the municipality level. Additionally to our spatial Durbin GDP per capita models based on five-year moving windows, two different sets of seven- and ten-year moving windows were applied to verify if pericoupling effects of GDP per capita growth would be different from zero, but no evidences of such effects were found.
Our empirical analyses demonstrate that in addition to a strong telecoupling process (i.e., the trade of soybean for national and international markets, Fig. 1) that has been reported to cause significant socioeconomic11,29,32 and environmental12,13,14,15 effects, both intracouplings and pericouplings also exert significant effects (Fig. 4). These interdependent systems include flows of not only soybean but also of information, resources, capital and people, among others27,28,33, all of which influence land-use decisions not only locally but also in adjacent municipalities.
Conceptualization of the intracoupled and pericoupled systems. Intracouplings and pericouplings are conceptualized as the flows of people, information, resources, and capital both within and across neighboring municipalities in response to the soybean market. The pericoupling flows, in turn, have effects on land change dynamics (e.g., pasture conversion, deforestation) in adjacent municipalities.
In contrast to traditional land change analyses that often treat socioeconomic and environmental outcomes by adopting a specific time interval (e.g., defined within a single period)34,35,36, our analyses evaluate deforestation and economic growth as dynamic processes through the use of a moving window approach, hence yielding fifteen subsets of model simulation. Furthermore, through the use of a Spatial Durbin Error Model (SDEM), these analyses provide a different perspective from previous studies that focused on general indirect land change processes pushing deforestation into far distant areas37,38. Hence, our study empirically addressed local (intracoupled) and adjacent (pericoupled) systems simultaneously by taking the municipality as the unit of analysis, and additionally applying the same modeling approach for both processes—economic growth and deforestation—while putting them into the same theoretical metacoupling framework.
The deforestation model shows that the expansion of soybean production was significantly related with deforestation not only within municipalities but also across neighboring municipalities. This is a key empirical finding as it demonstrates that the expansion of soybean production in a given municipality has pericoupling effects that influence land change dynamics in adjacent municipalities within the State of Mato Grosso, hence increasing their likelihood of experiencing deforestation. This land change process is especially reinforced by indirect land-use change in which the expansion of pastureland precedes the expansion of soybean13,34, with soybean areas replacing pastureland areas that had driven deforestation in neighboring municipalities. Furthermore, we demonstrated that neighboring municipalities are related regarding land-use decisions, evidencing that agglomeration economies in fact play an important role not only at the local level but for entire frontier regions39. As has been previously stated16, our analysis demonstrates that human-nature interactions within a system do affect adjacent systems. This result may be particularly explained by the flow of information influencing land-use decisions among neighboring localities—e.g., through the propagation of information in frontier regions40. Among major information flows, we highlight the positively correlated local economic growth promoted by soybean expansion (observed by the GDP per capita model), lower land prices in neighboring municipalities with large areas of arable land available41, lower integration with international markets and hence not committed with supply chain agreements for sustainability, thus facilitating deforestation42, the scarcity of arable lands in former producing municipalities pushing producers to expand into new adjacent areas, and the proximity to more developed hubs of production, which facilitate logistics and trade43.
In contrast, the GDP per capita model shows that the expansion of soybean production was significantly related with economic growth but only within municipalities and not across neighboring municipalities (i.e., no significant pericoupling effects associated with economic growth were found). During the period between 2000 and 2018, the production value of the soybean contributed, on average, to around 27% of municipalities’ GDP44, which demonstrates a strong contribution to the State’s economy. This result empirically confirms that while large agribusiness development hubs greatly benefit from the production of agricultural commodities, adjacent municipalities do not equally benefit45. We highlight that producers in Mato Grosso may purchase farms and open areas for soybean production in municipalities adjacent to well-established production centers, which reinforce economic development in those well-established centers, but not in the neighboring municipalities. This occurs because in these well-established centers producers can more easily: find storage facilities and trading offices; purchase machinery, seeds, and technological infrastructure; and find professional consultants45. Thus, municipalities adjacent to large soybean production centers, which constitute pericoupled systems since they are interconnected through many flows, including people (e.g., farmers), information, and capital, among others, experience significantly lower economic benefits than adjacent agribusiness development hubs45.
In addition to the broad distribution of environmental costs (e.g., deforestation within and in neighboring municipalities) associated with quite localized economic benefits (e.g., GDP per capita growth within but not in neighboring municipalities) of soybean production found in this study, Mato Grosso has experienced the entrance of large international groups purchasing land for soybean production since 200046,47 (which is negatively perceived by 95% of the producers from APROSOJA—the largest association of soybean producers in Mato Grosso48). Thus, Mato Grosso may be experiencing a similar outcome as other agricultural frontier areas such as those in the Chaco and Chiquitano woodlands of Bolivia and Paraguay, respectively. In those areas, new comers from adjacent countries (e.g., Brazil, Argentina) produce agricultural commodities that degrade the local environment without favoring the local economy and instead export the economic gains to investors in distant locations40. Hage et al.47 also highlight the risks of the economic gains from such production systems being exported to geographically distant systems.
Previous studies have shown that soybean production plays an important role in the economic development of soybean producing regions31,32,49, although with varying effects according to the region considered. For instance, soybean production exhibited a positive effect on socioeconomic metrics such as the Gini coefficient for income in the new agricultural frontier of the Amazon biome and MATOPIBA (northeast region within the Cerrado biome49), while in Mato Grosso no such effect was found32. The authors32 attribute these regionalized positive effects to regions such as MATOPIBA because those areas are economically underdeveloped and with a large gap in logistics and infrastructure. This situation has improved over the last couple of decades, hence improving the Gini and other metrics of human wellbeing. Similar to these previous studies, which did not evaluate the effects of soybean production on neighboring localities, our results show that the economic development brought by soybean production is limited to local areas and do not propagate to neighboring municipalities as was observed for deforestation. As noted by Lopes et al.50, the soybean development in MATOPIBA has promoted ecological exclusion of customary users (e.g., local people, smallholders) while taking over diverse former land uses, but also excluding those communities (e.g., indigenous, smallholders) from the socioeconomic outcomes of this agricultural activity. However, a Mann–Whitney U test (Supplementary Fig. 2) showed that soybean producing municipalities have significantly (p-values < 0.05) higher GDP per capita compared to non-producing municipalities, suggesting that over the long run neighboring municipalities affected by negative pericoupling effects (such as deforestation) may also benefit economically but over time periods longer than five years. This trend is also reinforced by our OLS results (Supplementary information 3) in which soybean growth was significantly correlated with increases in GDP per capita at the local production level. Hence, from a complex coupled human–environment systems perspective, our results show that the environmental costs of soybean production (measured as land changes) in both local and adjacent systems are driven by economic expectations, which constitute a feedback from previous economic outcomes (economic growth associated with soybean production), while the economic development of adjacent systems is conditioned on actual production (i.e., soybean yields) if those previous land changes are followed by soybean production (Fig. 4).
The metacoupling framework16 constitutes a suitable approach for both conceptually and empirically evaluating the influence of proximal and distant interactions on socioeconomic and environmental outcomes. Under the influence of telecoupling processes (Fig. 1), Tobler’s first law of geography, which states that “everything is related to everything else, but near things are more related than distant things”51, may be appropriate to explain pericoupling effects on land change processes in adjacent municipalities (Fig. 4) but not on socioeconomic development. Consequently, while socioeconomic development pushed by the agribusiness may be localized, the environmental costs are more broadly distributed. Therefore, the analysis of intracouplings and pericouplings within the context of a strong telecoupling process provides relevant implications for sustainability.
While the soybean agribusiness is in fact a motor of economic and social development31,32,50, our study highlights new implications of its associated environmental impacts (e.g., deforestation) and economic benefits (e.g., local GDP per capita growth). If the soybean dynamics influence the amount of newcomers in municipalities neighboring large agribusiness centers (i.e., municipalities with higher levels of agribusiness infrastructure45), then it is necessary to ensure the economic development at the newcomers’ local level (i.e., pericoupled systems), instead of favoring economic growth of the large agribusiness centers in detriment to the new emerging production areas. Therefore, if economic development does not spillover across broad regions, Federal and State legislators should design redistributive policies aimed at the development of those fronts of production (i.e., the pericoupled systems), hence returning the potential economic gains to local populations in the emerging production areas. In addition, public policies such as the farm credit should be designed in tandem with other policies (e.g., entrepreneurship policies52) to stimulate farmers to purchase production inputs at local levels (instead of purchasing in more developed agribusiness centers) and incentivize the growth of small-to-medium local businesses to meet producers’ demand, while also providing other services to their families while attracting more entrepreneurs.
Our spatial Durbin local econometric models allowed to simultaneously assess the effects of a telecoupling process on local (intracoupled) and adjacent (pericoupled) systems. Nevertheless, this methodological approach has limitations with respect to the geographical scope (i.e., looking at Queen 1st-order of contiguity adjacent systems) and thus require complementary studies using different geographic distances37, spatial structures53,54,55, and modelling approaches. In addition, as the pericoupling effects on economic growth may only manifest when soybean production follows land changes, they may be influenced by significant and long temporal-lag effects. Therefore, future studies exploring economic development in pericoupled systems, not only in our study area but also in many regions around the world, would benefit by time series analyses over longer temporal periods than those analyzed in this study. Such analyses would better inform socioeconomic policies in agricultural development areas. Furthermore, while our study considered the role of agribusiness centers as drivers of economic development, we mainly focused on how the expansion of soybean production in a given municipality influenced not only itself but also its neighbors. Therefore, future studies should evaluate how metacoupled systems are influenced by large socioeconomic players, such as large agribusiness conglomerates.
To assess the effects of intracouplings and pericouplings in a large soybean metacoupled system16 on socioeconomic and environmental outcomes, we used the Spatial Durbin Error Model (SDEM) based standard proposed by Vega and Elhorst56 that model selection should be done under the context of empirical applications. This model was chosen because it accounts for causal variable attributions at the same cross-section (hereafter intracoupling effects) in 141 municipalities (which are used as units of analysis), while also accounting for causal variable attributions on neighboring municipalities (hereafter “pericoupling effects”)57,58. The model provides quantitative estimates of the magnitude of both intracoupling and pericoupling effects and their statistical significance when they are not equal to zero58.
Our deforestation and GDP per capita models were designed to answer whether soybean production drives deforestation and changes in GDP per capita in a constant way or if their significance changes over time. Furthermore, they allow assessing whether the effects of soybean production are contained within the producer municipality or if they influence its neighbors in a contagious path dependence pattern43.
The soybean production of a given year requires previous years of preparation, which involves soil preparation, land-use and land-cover (LULC) changes, access to credit and financial resources, labor, thereby a set of situations to favor its success27,28,33. However, after the crop starts developing in a given area and year, other dynamics are expected to occur as a consequence, which may lead to different socioeconomic and environmental outcomes31. Given Tobler’s first law of geography51 we hypothesized that the effects of the soybean production in a region will influence neighboring areas through path dependencies. Hence, the intracoupling effects of economic development based on soybean production leads to land changes (e.g., deforestation) in neighbor municipalities (Fig. 4) in the following years, as new producers will expect future gains from their soybean production. Therefore, we should expect pericoupling effects for land changes but not for economic development, which is conditioned by the actual production at the focal system. We chose municipalities as our unit of analysis, given that this constitutes the lowest unit of analysis currently available. While previous studies have shown indirect effects of soybean production on deforestation in frontier areas37,38 they did not evaluate cross-sectional units of analysis, hence did not look at intracoupling and pericoupling effects simultaneously and on both land change and economic development processes. The publicly available dataset applied in this research is presented in Supplementary Table 1.
Based on a temporal moving window of five years, our models assess the dynamic effects of contemporary changes in soybean production, land use and socioeconomic dimensions on deforestation and GDP per capita. For this purpose, the t5 represents the dependent variable value of five-year periods (Δt5 = 5 years deforestation or GDP per capita change), while the explanatory variables xt0, xt1, xt5, where t0 representing time invariant variables, t1 the variable value at the beginning of the period, and t5 accounting for the change in value over the five-year period (Δt5). This approach has the following advantages of conventional methods: (1) allows assessing the temporal variability, (2) incorporate stochastic process, and (3) and deals with cointegration. Hence, our model approach explores temporal variability within a time period of analysis (e.g., 2000–2018) instead of looking at the entire period of analysis59 avoiding issues of convolution and non-stationarity60. The approach also allows one to incorporate stochastic processes into non-stationary time series, permitting the study of contemporary effects—i.e., short-term dynamics. Furthermore, by choosing moving windows of 5-year interval (15 in total for each model), this approach allowed one to smooth short-term fluctuations while also allowing the observation of long-term trends. Although many models consider stationarity in time series, stochastic (non-stationary) processes may occur and sometimes yield spurious regressions61,62. Since over the period of evaluation (2000–2018) the data time series of the present study may be non-stationary and lack cointegration among variables, our moving window approach is justified, since it allows cointegration between non-stationary variables, thus obtaining non-biased linear regression models63.
Since 2017, data on the Brazilian LULC is systematically produced by MapBiomas (, a multi-institutional initiative aiming at developing annual national LULC data at 30-m spatial resolution from 1985 to the present, through the automatic classification of satellite imagery. The present study uses the collection v.4.1, which has a global accuracy of 90%64. The original data encompass twenty-seven LULC classes (Level 2) and for the purpose of this study, the class ‘natural forest cover’, key to calculating deforestation (i.e., dependent variable), is represented by the sum of ‘Forest formation’ and ‘Savanna formation’ from MapBiomas. Deforestation was then calculated by the total deforestation in a moving window (Δt5) at the municipality level. When deforestation is not observed in a given municipality and period (sometimes showing growth of natural forest cover), we set the value to zero, hence modeling just deforestation dynamics.
The Gross Domestic Product (in Brazilian Real) is a measure of the total monetary value added by the sectors of agriculture, industry, services and public administration at the municipality level and calculated on an annual basis65. The dependent variable ‘GPD per capita change’ (i.e., Δt5) is calculated by the municipality’s GPD divided by the respective number of inhabitants in the same year (GDPpercapita = GPD/inhabitants). Population statistics were obtained from the IBGE Population Estimates service66. In addition to the five-year moving windows, for the GDP per capita we also built models with windows of 7 and 10 years to assess for temporal-lag effects.
Soybean production was considered a key explanatory variable of interest in our models, as it represents a major agribusiness activity in the state of Mato Grosso—i.e., key to understanding agricultural expansion and land changes. Hence, this variable is represented by ‘soybean production’ measured in tons. The initial ‘natural forest cover’—t1 (hectares) of each deforestation period is also included. The class ‘pasture’ (in hectares) used in the model corresponds to the ‘pasture’ class from the MapBiomas dataset. While agricultural intensification, a process expected to spare land from agricultural expansion and hence decrease deforestation trends67,68,69 has being observed by previous studies69,70, rebound effects may imply that higher yields from agricultural intensification instead of sparing land for conservation would foster more crop expansion (through deforestation and pasture conversion) given the great potential for economic returns71. Hence, our model accounts for ‘maize second-crop production’ (tons). Although indirectly, this approach allows evaluating the impacts of agricultural intensification70. Our model also considered GPD per capita since economic development and growth may foster land conversion and environmental degradation72. Topography was also considered an important factor as it drives large-scale agricultural production73. In particular, the variable ‘slope’ represents the mean value (percentage of slope) for each municipality. Slope information was derived from the Shuttle Radar Topography Mission (SRTM) elevation maps developed by Weber et al.74. All variables used in each model subset represent the change at time Δt5 and the value at t1, except ‘slope’, which is time invariant (t0) and maize second-crop, which is only represented by Δt5 because as long as this crop is dependent on soybean production, their initial time will always be spatially correlated.
In addition to the variables described in the deforestation model, the GDP per capita model also incorporated two additional variables: Formal ‘agricultural labor’ (number of laborers employed in the agribusiness) and ‘cattle herd’ (number of animals). The agricultural labor was obtained from the Annual List of Social Information (RAIS in Portuguese), a database of the Ministry of Labor and Employment (MTE in Portuguese)75. The cattle herd was obtained from the IBGE Municipal Cattle Survey, which is a systematic annual census of properties in the Brazilian countryside76.
To evaluate pericoupling effects of municipalities on its neighbors, we used the Spatial Durbin Error Model (SDEM), which accounts for the spatial lag of independent variables77. The SDEM subsumes the Spatial Error model (SEM, Eq. 1) and the Spatial Lag of X model (SLX, Eq. 2), and described by (Eq. 3)57:
where β2 is the cross-partial derivative or indirect effects (i.e., pericoupling) of the neighboring regions in the W matrix, and β1 the own-region direct effect (i.e., intracouplings). When results from cross-partial derivatives are not equal to zero, it implies impact on neighboring municipalities without endogenous feedback effects, which indicates true local spatial lag58. Hence, we deal with cross-sectional dependence with a balanced data of 141 cross-sectional units (i = municipalities). ɛ is the disturbance vector of random errors from the regression model, which considers an autoregressive spatial process in the error term for the spatial model, where λ the autoregressive coefficient of the error terms; while u is a vector of independent and identically distributed error terms78,79,80. As ordinary least-squares regression (OLS) is unsuitable for spatial regression models because it assumes independence among observations60,81, the proposed SDEM relies on the maximum likelihood estimation77, suitable to estimate the significance and magnitude of spatial lags59,81. Therefore, spatial lags in the proposed model approach constitute the pericoupling effects of independent variables X of a given municipality j over the dependent variable y in the neighboring municipality i. This rationale resembles the flows and mutual feedbacks between municipalities as contextualized in Fig. 4. The spatial weights matrix (W) at the municipality level was created by defining a neighbor based on the Queen 1st-order of contiguity approach, adjacency matrix82,83. As the SDEM does not yield traditional correlation estimates such OLS regression, the independent explanatory variables’ effects are measured by the average impacts (direct, indirect, and total) and their statistical significances estimated by the Bayesian Markov Chain Monte Carlo method, using the random walk Metropolis algorithm set with 500 simulations83. For further model interpretation, see LeSage77.
To all subsets of the deforestation and GDP per capita models (fifteen subsets each), all the explanatory variables were tested by the variation inflation factor (VIF). The study considered a variable critical for multicollinearity if VIF > 5.084. To each moving window (model subset) we applied an OLS model and tested its residuals with the Global Moran’s I85 to check overall spatial autocorrelation where the null-hypothesis indicates no spatial autocorrelation. A cluster and outlier analysis based on the Anselin Local Moran’s I86 was then applied to evaluate clustering patterns over the spatial distribution among municipalities for the overall changes of the three major variables of interest (i.e., GDP per capita, deforestation, and soybean production). To test the assumption of non-stationarity in the time series, we applied the augmented Dickey-Fuller (ADF) test over our dependent and independent variables (performed with aggregated data at State level)87. The null-hypothesis indicates non-stationarity, which means the variable data has unit root and the time series data follows a certain trend. To complement the analysis, the Engle-Granger cointegration test was applied by using the ADF assessed cointegration between dependent and independent variables88. To check model restrictions the Likelihood Ratio test89 was also applied to each of the models’ subsets, to verify if the SDEM has a better fit than an OLS. This test assumes that restricted coefficients are like zero (null-hypothesis), indicating that a less complex model is more suitable, thus suggesting that restrictions are true and a simpler model (e.g., OLS) may have a better fit. Finally, the non-parametric Mann–Whitney U test was applied (normal distribution of data was rejected by a Shapiro–Wilk test using statistical significance threshold of p-value < 0.0190) to the GDP per capita in 2000, 2005, 2010, 2015, and 2018, to verify significant statistical differences in GDP per capita values between soybean producing and non-producing municipalities. Statistical tests’ results for model specifications are presented in Supplementary Information 4.
Liu, J. et al. Framing sustainability in a telecoupled world. Ecol. Soc. 18(2), 26 (2013).
CAS  Article  Google Scholar 
Meyfroidt, P., Rudel, T. K. & Lambin, E. F. Forest transitions, trade, and the global displacement of land use. Proc. Natl. Acad. Sci. 107, 20917–20922 (2010).
ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 
Lambin, E. F. & Meyfroidt, P. Global land use change, economic globalization, and the looming land scarcity. Proc. Natl. Acad. Sci. 108, 3465–3472 (2011).
ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 
Food and Agriculture Organization. FAOSTAT data. Rome: FAO 2020;
Liu, J. Consumption Patterns and Biodiversity. The Royal Society (2020).
UNComtrade. United Nations Statistics Division, 2019;
Oliveira, G. D. L. & Schneider, M. The politics of flexing soybeans: China, Brazil and global agroindustrial restructuring. J. Peasant Stud. 43, 167–194 (2016).
Article  Google Scholar 
York, R. & Gossard, M. H. Cross-national meat and fish consumption: Exploring the effects of modernization and ecological context. Ecol. Econ. 48, 293–302 (2004).
Article  Google Scholar 
Solot, I. B. The Chinese agricultural policy trilemma. Perspectives 7, 36–46 (2006).
Google Scholar 
Garrett, R. D. & Rausch, L. L. Green for gold: Social and ecological tradeoffs influencing the sustainability of the Brazilian soy industry. J. Peasant Stud. 43, 461–493 (2016).
Article  Google Scholar 
Silva, R. F. B. et al. The Sino-Brazilian telecoupled soybean system and cascading effects for the exporting country. Land 6(3), 53 (2017).
Article  Google Scholar 
DeFries, R. S., Uriarte, M., Rudel, T. & Hansen, M. Deforestation driven by urban population growth and agricultural trade in the twenty-first century. Nat. Geosci. 3, 178–181 (2010).
ADS  CAS  Article  Google Scholar 
Macedo, M. N. et al. Decoupling of deforestation and soy production in the southern Amazon during the late 2000s. Proc. Natl. Acad. Sci. 109, 1341–1346 (2012).
ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 
Nepstad, D. et al. Slowing Amazon deforestation through public policy and interventions in beef and soy supply chains. Science 344, 1118–1123 (2014).
ADS  CAS  PubMed  Article  PubMed Central  Google Scholar 
Gasparri, N. I. & de Waroux, Y. L. P. The coupling of South American soybean and cattle production frontiers: New challenges for conservation policy and land change science. Conserv. Lett. 8, 290–298 (2015).
Article  Google Scholar 
Liu, J. Integration across a metacoupled world. Ecol. Soc. 22, 29 (2017).
Article  Google Scholar 
Liu, J. et al. China’s environment on a metacoupled planet. Annu. Rev. Environ. Resour. 43, 1–34 (2018).
CAS  Article  Google Scholar 
Herzberger, A., Chung, M. G., Kapsar, K., Frank, K. A. & Liu, J. Telecoupled food trade affects pericoupled trade and intracoupled production. Sustainability 11, 2908 (2019).
Article  Google Scholar 
Dou, Y. et al. Understanding how smallholders integrated into pericoupled and telecoupled systems. Sustainability 12(4), 1596 (2020).
Article  Google Scholar 
Tromboni, F. et al. Macrosystems as metacoupled human and natural systems. Front. Ecol. Environ. 19(1), 20–29 (2021).
Article  Google Scholar 
Zhao, Z. et al. Synergies and tradeoffs among sustainable development goals across boundaries in a metacoupled world. Sci. Total Environ. 751, 141749 (2021).
ADS  CAS  PubMed  Article  Google Scholar 
Hovis, C. L., Dou, Y., Herzberger, A. & Liu, J. Through the lens of telecoupling and metacoupling: New perspectives for global sustainability. Sustainability 13(5), 2953 (2021).
Article  Google Scholar 
Zhang, J. et al. Sustainability evaluation on the grain to green program in the Hexi Corridor of China: A metacoupled system perspective. Sustainability 13(3), 1498 (2021).
Article  Google Scholar 
Xu, Z. et al. Impacts of irrigated agriculture on food-energy-water-CO2 nexus across metacoupled systems. Nat. Commun. 11, 5837 (2020).
ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 
Carlson, A. K., Taylor, W. W., Rubenstein, D. I., Levin, S. A. & Liu, J. Global marine fishing across space and time. Sustainability 12, 4714 (2020).
Article  Google Scholar 
Estimativas de Comércio Exterior do Agronegócio Brasileiro (AGROSTAT—Ministry of Agriculture, accessed November 2020);
Dou, Y. et al. Land-use changes in distant places: Implementation of a Telecoupling Agent-Based Model. J. Artif. Soc. Soc. Simul. 23, 11 (2020).
Article  Google Scholar 
Dou, Y. et al. Land-use change across distant places: Design of telecoupling agent-based model. J. Land Use Sci. 14, 191–209 (2019).
Article  Google Scholar 
Richards, P., Pellegrina, H., VanWey, L. & Spera, S. Soybean development: The impact of a decade of agricultural change on urban and economic growth in Mato Grosso, Brazil. PLoS ONE 10(4), e0122510 (2015).
PubMed  PubMed Central  Article  CAS  Google Scholar 
Picoli, M. C. A. et al. Impacts of public and private sector policies on soybean and pasture expansion in Mato Grosso—Brazil from 2001 to 2017. Land 9(1), 20 (2020).
Article  Google Scholar 
Song, X. P. et al. Massive soybean expansion in South America since 2000 and implications for conservation. Nat. Sustain. (2021).
Article  PubMed  PubMed Central  Google Scholar 
Martinelli, L. A., Batistella, M., Silva, R. F. B. & Moran, E. Soy expansion and socioeconomic development in municipalities of Brazil. Land 6(3), 62 (2017).
Article  Google Scholar 
Silva, R. F. B., Batistella, M., Moran, E., Celidonio, O. L. M. & Millington, J. D. A. The soybean trap: Challenges and risks for Brazilian producers. Front. Sustain. Food Syst. 4, 12 (2020).
Article  Google Scholar 
Sá, S. A., Palmer, C. & Falco, S. Dynamics of indirect land-use change: Empirical evidence from Brazil. J. Environ. Econ. Manag. 65, 377–393 (2013).
Article  Google Scholar 
Schmidt, J. H., Weidema, B. P. & Brandão, M. A framework for modelling indirect land use change in life cycle assessment. J. Clean. Prod. 99, 230–238 (2015).
Article  Google Scholar 
Magliocca, N. R., Khuc, Q. V., Bremond, A. & Ellicott, E. Direct and indirect land-use change caused by large-scale land acquisitions in Cambodia. Environ. Res. Lett. 15, 024010 (2020).
ADS  Article  Google Scholar 
Arima, E. Y., Richards, P., Walker, R. & Caldas, M. M. Statistical confirmation of indirect land use change in the Brazilian Amazon. Environ. Res. Lett. 6, 024010 (2011).
ADS  Article  Google Scholar 
Richards, P. D., Walker, R. T. & Arima, E. Y. Spatially complex land change: The indirect effect of Brazil’s agricultural sector on land use in Amazonia. Glob. Environ. Change 29, 1–9 (2014).
PubMed  PubMed Central  Article  Google Scholar 
Henderson, J., Godar, J., Frey, G. P., Borner, J. & Gardner, T. The Paraguayan Chaco at a crossroads: Drivers emerging soybean frontier. Reg. Environ. Change 21, 72 (2021).
Article  Google Scholar 
Waroux, Y. L. Capital has no homeland: The formation of transnational producer cohorts in South America’s commodity frontiers. Geoforum 105, 131–144 (2019).
Article  Google Scholar 
IEG/FNP. Agrianual 2017: Anuária da agricultura brasileira (IEG/FNP, São Paulo, 2017).
Google Scholar 
Ermgassen, E. K. H. J. et al. The origin, supply chain, and deforestation risk of Brazil’s beef exports. Proc. Natl. Acad. Sci. 117, 317770–331779 (2020).
Article  CAS  Google Scholar 
Millington, J. D. A., Katerinchuck, V., Silva, R. F. B., Victoria, D. C. & Batistella, M. Modelling drivers of Brazilian agricultural change in a telecoupled world. Environ. Model. Softw. 139, 105024 (2021).
Article  Google Scholar 
Brazilian Institute of Geography and Statistics (IBGE). Pesquisa agrícola municipal (IBGE Sidra—Sistema IBGE de Recuperação Automática, 2019);
Araújo, F. A. O. M. & Boaventura, D. M. R. Municipalities of agribusiness in the State of Mato Grosso: Effects and defects of the activity of modern agriculture. Revista Geográfica Acadêmica 14(2), 100–122 (2020).
Google Scholar 
Contag (Confederação Nacional dos Trabalhadores na Agricultura). Grupos Argentinos Compram Terra para Soja no Mato Grosso. Contag Agência de Notícias (2008).
Hage, F. A. S., Peixoto, M. & Vieira Filho, J. E. R. Aquisição de terras por estrangeiros no Brasil: Uma avaliação jurídica e econômica. Textos para Discussão/Ipea (2012).
Cadore, F. Aquisição de terras por estrangeiros na ótica do produtor rural. Canal Rural (2012).
Silva, R. F. B. et al. Three decades of changes in Brazilian municipalities and their food production systems. Land 9(11), 422 (2020).
Article  Google Scholar 
Lopes, G. R., Lima, M. G. B. & Reis, T. N. P. Maldevelopment revisited: Inclusiveness and social impacts of soy expansion over Brazil’s Cerrado in Matopiba. World Dev. 139, 105316 (2021).
Article  Google Scholar 
Tobler, W. R. A computer movie simulating urban growth in the Detroit Region. Econ. Geogr. 46, 234–240 (1970).
Article  Google Scholar 
Okinyemi, F. O. & Adejumo, O. O. Government policies and entrepreneurship phases in emerging economies: Nigeria and South Africa. J. Glob. Entrep. Res. 8, 35 (2018).
Article  Google Scholar 
Elhorst, J. P., Gross, M. & Tereanu, E. Cross-sectional dependence and spillovers in space and time: Where spatial econometrics and global var models meet. J. Econ. Surv. 35(1), 192–226 (2021).
Article  Google Scholar 
Guan, H. & Li, Q. Spatial spillover effects of economic growth based on high-speed railways in Northeast China. Complexity 2021, 8831325 (2021).
Google Scholar 
Bai, J. & Li, K. Dynamic spatial panel data models with common shocks. J. Econom. 224, 134–160 (2021).
MathSciNet  MATH  Article  Google Scholar 
Vega, S. H. & Elhorst, J. P. The SLX model. J. Reg. Sci. 55, 339–363 (2015).
Article  Google Scholar 
Takagi, D. & Shimada, T. A spatial regression analysis on the effect of neighborhood-level trust on cooperative behaviors: Comparison with a multilevel regression analysis. Front. Psychol. 10, 2799 (2019).
PubMed  PubMed Central  Article  Google Scholar 
LeSage, J. P. What regional scientists need to know about spatial econometrics. SSRN (2014).
Article  Google Scholar 
Tientao, A., Legros, D. & Pichery, M. C. Technology sipillover and TFP growth: A spatial Durbin model. Int. Econ. 145, 21–31 (2016).
Article  Google Scholar 
Zhu, Z. & Wu, Y. Estimation and prediction of a class of convolution-based spatial nonstationary models for large spatial data. J. Comput. Graph. Stat. 19, 74–95 (2010).
MathSciNet  CAS  Article  Google Scholar 
Granger, C. W. J. & Newbold, P. Spurious regression in econometrics. J. Econom. 2, 111–120 (1974).
MATH  Article  Google Scholar 
García, V. R., Gaspart, F., Kastner, T. & Meyfroidt, P. Agricultural intensification and land use change: Assessing country-level induced intensification, land sparing and rebound effect. Environ. Res. Lett. 15, 085007 (2020).
ADS  Article  Google Scholar 
Engle, R. F. & Granger, C. W. J. Co-integration and error correction: Representation, estimation, and testing. Econometrica 55, 251–276 (1987).
MathSciNet  MATH  Article  Google Scholar 
Souza, C. M. et al. Reconstructing three decades of land use and land cover in Brazilian biomes with Landsat archive and earth engine. Remote Sens. 12, 2735 (2020).
ADS  Article  Google Scholar 
Brazilian Institute of Geography and Statistics (IBGE). Produto Interno Bruto dos Municípios (IBGE Sidra—Sistema IBGE de Recuperação Automática, 2018);
Brazilian Institute of Geography and Statistics (IBGE). Estimativas da População (IBGE Sidra—Sistema IBGE de Recuperação Automática, 2020);
Bennett, E. M. Changing the agriculture and environment conservation. Nat. Ecol. Evol. 1, 0018 (2017).
Article  Google Scholar 
Cassman, K. G. & Grassini, P. A global perspective on sustainable intensification research. Nat. Sustain. 3, 262–268 (2020).
Article  Google Scholar 
Pelletier, J., Ngoma, H., Mason, N. M. & Barrett, C. B. Does smallholder maize intensification reduce deforestation? Evidence from Zambia. Glob. Environ. Change 63, 102127 (2020).
Article  Google Scholar 
Garrett, R. D. et al. Intensification in agriculture-forest frontiers: Land use responses to development and conservation policies in Brazil. Glob. Environ. Change 53, 233–243 (2018).
Article  Google Scholar 
Paul, C., Techen, A., Robinson, J. S. & Helming, K. Rebound effects in agricultural land and soil management: Review and analytical framework. J. Clean. Prod. 227, 1054–1067 (2019).
CAS  Article  Google Scholar 
Abel, C. et al. The human-environment nexus and vegetation-rainfall sensitivity in tropical drylands. Nat. Sustain. 4, 25–32 (2021).
Article  Google Scholar 
Jasinski, E., Morton, D. & DeFries, R. Physical landscape correlates of the expansion of the mechanized agriculture in Mato Grosso, Brazil. Earth Interact. 9, 16 (2005).
Article  Google Scholar 
Weber, E., Hasenack, H. & Ferreira, C. J. S. Adaptação do modelo digital de elevação do SRTM para o sistema de referência oficial brasileiro e recorte por unidade da federação (Porto Alegre, UFRGS Centro de Ecologia, 2004);
Fernandes, M. D. E., Gouveia, A. B. & Benini, E. G. Teachers’pay in Brazil: An outlook from the Annual List of Social Information (RAIS). Educ. Pesqui. 38, 339–356 (2012).
Article  Google Scholar 
Brazilian Institute of Geography and Statistics (IBGE). Pesquisa Pecuária Municipal (IBGE Sidra—Sistema IBGE de Recuperação Automática, 2019);
LeSage, J. P. An introduction to spatial econometrics. Revue D’économie industrielle 123, 19–44 (2008).
Article  Google Scholar 
Dantas, R. A., Magalhães, A. M. & Vergolino, J. R. O. Um modelo espacial de demanda habitacional para a Cidade do Recife. Estudos Econômicos 40, 891–916 (2010).
Google Scholar 
Hession, S. L. & Moore, N. A spatial regression analysis of the influence of topography on monthly rainfall in East Africa. Int. J. Climatol. 31, 1440–1456 (2011).
Article  Google Scholar 
Golgher, A. B. & Voss, P. R. How to interpret the coefficients of Spatial Models: Spillovers, direct and indirect effects. Spat. Demogr. 4, 175–205 (2016).
Article  Google Scholar 
Feng, Y., Wang, X., Du, W. & Liu, J. Effects of Ais pollution control on urban development quality in Chinese cities based on spatial Durbin model. Int. J. Environ. Res. Public Health 15, 2822 (2018).
PubMed Central  Article  Google Scholar 
Viña, A., McConnell, W. J., Yang, H., Xu, Z. & Liu, J. Effects of conservation policy on China’s forest recovery. Sci. Adv. 2, e1500965 (2016).
ADS  PubMed  PubMed Central  Article  Google Scholar 
Bivar, R. et al. Spatial Regression Analysis (Package ‘spatialreg’, 2019);
Salmerón, R., García, C. B. & García, J. Variance inflation factor and condition number in multiple linear regression. J. Stat. Comput. Simul. 88, 2365–2384 (2018).
MathSciNet  MATH  Article  Google Scholar 
Cliff, A. D. & Ord, K. K. Spatial Processes: Model and Applications (Pion, 1981).
MATH  Google Scholar 
Anselin, L. Local indicators of spatial association-LISA. Geogr. Anal. 27, 93–115 (1995).
Article  Google Scholar 
Lee, C. & Chang, C. New evidence on the convergence of per capita carbon dioxide emissions from panel seemingly unrelated regressions augmented Dickey-Fuller tests. Energy 33, 1468–1475 (2008).
CAS  Article  Google Scholar 
Agunloye, O. K., Shangodoyin, D. K. & Arnab, R. Lag length specification in Engle-Granger cointegration test: A modified Kouck Mean Lag approach based on partial correlation. Stat. Transit. 15, 559–572 (2014).
Google Scholar 
MacKenzie, D. I. et al. Occupancy Estimation and Modeling: Inferring Patterns and Dynamics of Species Occurrence 2nd edn. (Academic Press, 2018).
MATH  Google Scholar 
Silva, R. F. B., Millington, J. D. A., Moran, E. F., Batistella, M. & Liu, J. Three decades of land-use and land-cover change in mountain regions of the Brazilian Atlantic Forest. Landsc. Urban Plan. 204, 103948 (2020).
Article  Google Scholar 
Download references
We thank the Center for Environmental Studies and Research of the State University of Campinas, Brazil, the Center for Global Change and Earth Observations and the Center for Systems Integration and Sustainability as well as AgBioResearch at Michigan State University, USA, for providing support for this study. We acknowledge the funding support from the National Science Foundation, Grant Numbers 1531086 and 1924111. We also thank FAPESP, which also provided support through the processes 14/50628-9, 15/25892-7, and 18/08200-2. We are very grateful to MapBiomas who generated the land-use and and-cover data used in this paper (Project MapBiomas—Collection v.4.1 of Brazilian Land Cover & Use Map Series, accessed on June 10, 2020 through the link: None of these funding sources are to be held responsible for the opinions and views expressed herein. They are the sole responsibility of the authors. We are very thankful to two anonymous reviewers for their constructive suggestions that improved the manuscript’s clarity.
Center for Systems Integration and Sustainability, Department of Fisheries and Wildlife, Michigan State University, East Lansing, MI, 48823, USA
Ramon Felipe Bicudo da Silva, Andrés Viña & Jianguo Liu
Center for Environmental Studies and Research, State University of Campinas, Campinas, 13083-867, Brazil
Ramon Felipe Bicudo da Silva & Mateus Batistella
Department of Geography, University of North Carolina, Chapel Hill, USA
Andrés Viña
Center for Global Change and Earth Observations, Michigan State University, East Lansing, MI, 48823, USA
Emilio F. Moran
Environmental Geography Group, Institute for Environmental Studies, Vrije Universiteit Amsterdam, de Boelelaan, 1081HV, Amsterdam, The Netherlands
Yue Dou
Embrapa Agricultural Informatics, Brazilian Agricultural Research Corporation, Campinas, 13083-886, Brazil
Mateus Batistella
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
R.F.B.S. and A.V. led the manuscript conceptualization, methodology, data-gathering, software, investigation, formal analysis, and the writing—original draft, writing-review and editing. R.F.B.S created all the maps and figures presented in the article. E.F.M., Y.D., M.B., and J.L. contributed to the research, conceptualization, and the manuscript writing—reviewing and editing.
Correspondence to Ramon Felipe Bicudo da Silva.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit
Reprints and Permissions
da Silva, R.F.B., Viña, A., Moran, E.F. et al. Socioeconomic and environmental effects of soybean production in metacoupled systems. Sci Rep 11, 18662 (2021).
Download citation
Received: 04 June 2021
Accepted: 06 September 2021
Published: 20 September 2021
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advanced search
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
© 2021 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.


Leave a Comment

Your email address will not be published. Required fields are marked *