Crop Rotation and Climate Change Adaptation in Argentina’s Agriculture Sector
Crop Rotation and Climate Change Adaptation in Argentina’s Agriculture Sector
- Ariel R. Angeli, Ariel R. AngeliAsociación Argentina de Consorcios Regionales de Experimentación Agrícola
- Federico E. Bert, Federico E. BertInter-American Institute for Cooperation on Agriculture
- Sandro Díez-Amigo, Sandro Díez-AmigoInternational Finance Corporation
- Yuri Soares, Yuri SoaresInter-American Development Bank
- Jaquelina M. Chaij, Jaquelina M. ChaijAsociación Argentina de Consorcios Regionales de Experimentación Agrícola
- Gustavo D. Martini, Gustavo D. MartiniAsociación Argentina de Consorcios Regionales de Experimentación Agrícola
- F. Martín Montané, F. Martín MontanéAsociación Argentina de Consorcios Regionales de Experimentación Agrícola
- Alejandro Pardo VegezziAlejandro Pardo VegezziInter-American Development Bank
- and Federico SchmidtFederico SchmidtAsociación Argentina de Consorcios Regionales de Experimentación Agrícola
Summary
During the past two decades, extensive agriculture, particularly soybean production, has progressively replaced other crops in Argentina. This transformation was driven by economic, technological, environmental, and organizational factors, such as the increasing demand for agricultural commodities, technological advances, organizational innovations, and climate fluctuations. The expansion of soybean production has brought a substantial increase in agricultural revenue for Argentina. However, the predominance of soybean cultivation poses significant challenges, such as diminished soil fertility, reduction and increased variability in crop yields, ecological imbalance, increased greenhouse gas (GHG) emissions, and vulnerability to climate change.
Crop rotation, particularly balanced crop rotation, may result in very large positive impacts on soybean yields, especially in unfavorable climatic conditions such as those experienced during the La Niña ENSO phase in Argentina. In addition to this positive impact on agricultural productivity and climate adaptation, in some contexts crop rotation may also contribute to the reduction of GHG emissions, increased input energy efficiency, and improved environmental outcomes.
The 2018 Argentinian Association of Regional Consortia for Agricultural Experimentation and Inter-American Development Bank (AACREA-IADB) Integrated Crop Rotation Database compiled and harmonized the information from agricultural diaries kept by Regional Consortia for Agricultural Experimentation (CREA) members in Argentina from 1998 to 2016. This new consolidated data set has replaced previous regional templates, and it is expected to continue to be expanded with new information periodically, offering opportunities for further research on the impact of crop rotation on climate adaptation and on other topics in agricultural and environmental economics.
Keywords
Subjects
- Quantitative Analysis and Tools
- Agriculture and the Environment
The Essential Role of Agriculture in Climate Change Mitigation and Adaptation in Argentina
This article reviews the evolution of Argentina’s agriculture sector, including drivers of change and main challenges, and provides an overview of how crop rotation in particular can facilitate agricultural adaptation to climate change in Argentina. In order to do so, the impact of crop rotation on agricultural yields in more and less favorable climatic scenarios is considered. In addition, the impact of crop rotation on climate mitigation, the environment, and agricultural profitability is also explored.
It is worth noting that, as befits an encyclopedia article, this is not intended to be a research-driven review. Rather, its objective is to provide a rich perspective on the subject that allows readers to subsequently peruse and undertake deeper research on their own. Therefore, the focus of the discussion is to summarize the main insights from the existing body of literature in the relevant areas of agricultural and environmental economics, as well as agronomy. This includes studies published in Spanish, in addition to publications in English.
However, because so far there are no comprehensive large-scale studies of the impact of crop rotation in Argentina, the exposition also relies on the authors’ own analysis of a new nationwide plot-level data set compiled by the Argentinian Association of Regional Consortia for Agricultural Experimentation (AACREA) and the Inter-American Development Bank (IADB)—the 2018 AACREA–IADB Integrated Crop Rotation Database. In order to facilitate the use of this new resource, the structure of the database is described, and sample analyses are showcased. Similarly, key methodological challenges are outlined, and potential avenues for further research are suggested.
The Evolution of Argentina’s Agriculture Sector
The use of land and the structure of agricultural production in developing economies have evolved substantially since the Green Revolution of the 1960s, which witnessed the widespread adoption of modern technologies such as mechanization, irrigation, agrochemicals, and crop high-yielding varieties (HYVs) (Hazell, 2009).
During the past two decades, agricultural production in Argentina has further evolved, with extensive agriculture, particularly soybean production, progressively replacing other crops (e.g., pasture, maize, and cereals) and activities (e.g., livestock farming). This trend could be observed throughout the country, and during the 2010s more than 70% of cultivated land was devoted to soybean crops. However, the transformation was particularly marked in environmentally fragile areas such as the northern provinces (Salta, Chaco, and Santiago del Estero). In these areas, soybean cultivation displaced traditional crops and the native ground cover, and it disrupted customary rotation practices between maize, cereals, and pasture (Magrin et al., 2005).
In the context of Argentina, each agricultural season has two main planting periods: summer and winter. The main summer crops are soybean and maize (and, in some regions, sunflower), and the main winter crops are wheat and barley. According to Satorre and Andrade (2021), the majority of plots in Argentina (~70% as of 2021) are cultivated only using summer crops, without planting any winter crops. Plots planted with winter crops (~30% as of 2021) are typically also planted with a summer crop. However, because the cultivation of a winter crop requires delaying the planting of the summer crop, this usually results in lower yields (i.e., of late summer crops compared to early summer crops).
Although soybean has generally enjoyed higher export prices, since 2015 it has been negatively affected by changes in export taxes, resulting in a slight decline in soybean cultivation. In particular, maize and wheat benefited from the elimination of their export taxes in December 2015, from a previous level of 25%. However, soybean export taxes were only reduced from 35% to 30%. In addition, the elimination of export taxes was accompanied with a “normalization” of the maize and wheat markets, lifting export limits and ensuring that Argentinian farmers could sell their grains without restrictions.
As a consequence of these changes in export taxes, soybean cultivated area decreased from 20 million hectares during the 2015–2016 season to 17 million hectares during the 2019–2020 season. Conversely, the land area devoted to maize and wheat cultivation has expanded very significantly: Maize and wheat cultivated area increased almost 50% to 15.5 million hectares during the 2019–2020 season compared to an average of 10.5 million hectares during the 2015–2016 season (Ministerio de Agricultura, Ganadería y Pesca [MAGYP], 2021) (Figure 1).

Figure 1. Planted area by crop type—Argentina (2003–2020).
Drivers of Change
The transformation of Argentina’s agriculture sector was driven by economic, technological, environmental, and organizational factors:
The increasing demand for agricultural commodities improved the relative profitability of agriculture in general, and of soybean cultivation in particular.
Technological advances, such as direct seeding and genetically modified crops, drove down the cost and simplified the logistics of agricultural activities.
Organizational innovations in the management of agricultural business led to increased professionalization and economies of scale (e.g., the late 1990s witnessed the rapid development of “sowing pools,” investment funds that provide financial and agronomic management for large-scale agricultural production, contracting land to third parties, using diversification to mitigate risks, and coordinating input purchases in order to minimize costs).
Favorable climate fluctuations resulting in an increase in rainfall, particularly in semi-arid regions, allowed for the expansion of agricultural activities to areas previously considered nonviable.
Increasing Demand for Agricultural Commodities
The increasing demands of a growing world population and a change in the quality of diets have made it necessary to increase agricultural production (Andrade, 2020). The Argentinian Pampas are among the most naturally fertile regions in the world (Calviño & Monzón, 2009), and Argentina has a long tradition in extensive agriculture. Therefore, Argentina is a key player in the quest for global food security.
At the beginning of the 20th century, Argentinian producers already planted 12 million hectares of maize, wheat, and flax (Hora, 2012). The increasing demand for animal protein and biofuels during the 20th century and early 21st century further expanded the market for Argentine grains (Lamers et al., 2008). In combination with other factors, this triggered a process of expansion of grain and oilseed production.
As of 2021, Argentinian producers planted 35 million hectares, predominantly with crops of maize, soybean, and wheat (Satorre & Andrade, 2021). Commodity prices act as the main driver of Argentinian farmers’ decision-making, because the agriculture sector in Argentina is not subsidized, and therefore farmers fully take on the risks inherent to commodity markets.
Technological Advances
Agricultural systems in the Argentinian Pampas underwent significant changes during the last decades of the 20th century and the first decades of the 21st century. This was partly driven by the increased demand of most agricultural commodities produced in Argentina, but also by the evolving local economic and political context, fluctuations in climate and weather, and technological innovations (Bert et al., 2011).
The intertwined effects of these drivers induced significant changes in land use and land tenure patterns in Argentinian agricultural systems (Baldi & Paruelo, 2008). The most remarkable land use change was the expansion—westward and northward—of the area devoted to annual crops, replacing and displacing pastures and native grasslands to marginal areas (Magrin et al., 2005; Viglizzo et al., 2011).
Most of this expansion was tied to the advance of soybean cultivation, which quickly became the dominant crop in most agricultural regions. After its introduction to the country in the early 1970s, soybean production in Argentina reached 11 million tons in 1990, with 5.1 million hectares planted.
The introduction of no-tillage planting and genetically modified herbicide-tolerant soybean varieties in the early 1990s resulted in large cost reductions due to improved weed control, lower energy costs, and simplified agronomic management (Qaim & Traxler, 2005; Trigo & Cap, 2003).
These technological advances resulted in a sustained rapid expansion of soybean cultivation, and as of 2015, Argentinian producers planted 20 million hectares of the crop, yielding 55 million tons.
Organizational Innovations
An additional driver of changes in land use in Argentina has been the rapid expansion of land leasing. As of 2021, approximately half of the area cultivated in the Argentinian Pampas is managed by tenant farmers (Instituto Nacional de Estadística y Censos, 2021). Favorable economic conditions for agricultural production, together with the advent of simplified agronomic systems based on new technologies, have provided incentives for farmers to increase their production scale through land rent and the outsourcing of labor and services. However, advances in technology have also left some farmers out.
The existing evidence indicates that leased agricultural land is managed differently from owned land (Carolan, 2005). Studies of agricultural production in the Argentinian Pampas (e.g., Arora et al., 2016; Senesi et al., 2017) suggest that landowners are more likely to implement agricultural practices that can have a positive medium-term impact. Tenant farmers, on the other hand, are more likely to focus on agronomic practices that allow them to maximize immediate profits, particularly when their land leases are short term.
Given the expansion of land tenancy in Argentina’s agriculture sector, the divergence in incentives between landowners and tenants has resulted in changes in overall agronomic practices.
Favorable Climate Fluctuations
During the 20th century, central-eastern Argentina experienced some of the most significant changes in precipitation trends (Giorgi, 2002). Due to the high cost of irrigation, more than 80% of the agriculture in the Argentinian Pampas is rainfed, and the performance of crops is strongly tied to climate variations (Penalba et al., 2007; Verón et al., 2015). Accordingly, the long-frequency climate variation has contributed to large changes in land use patterns.
A steady increase in both annual and extreme precipitation has been observed in the Pampas since the 1960s (Boulanger et al., 2005; Minetti et al., 2003; Rusticucci & Penalba, 2000). Rainfall variation has been distributed unevenly through the seasonal cycle: Increased precipitation has been concentrated in late spring to summer, whereas the winter season has seen little or no change in precipitation. Furthermore, the increase in rainfall has been particularly marked near the western margin of the Argentinian Pampas, displacing westward the transition to semi-arid regions that represented the traditional boundary of rainfed agriculture and allowing for areas previously considered marginal to become arable (Berbery et al., 2006).
Challenges
The expansion of soybean production has brought substantial benefits to Argentina, as it has been associated with a large increase in agricultural revenue, particularly from exported crops: As of 2021, the soybean value chain accounts for almost 30% of total Argentinian exports (Instituto Nacional de Estadística y Censos, 2021).
However, the predominance of soybean cultivation also poses significant challenges because it can lead to a decrease in the frequency of crop rotation or to outright monoculture practices, with several potential negative effects on the economy, society, and the environment (Altieri, 2009; Kessler et al., 2007; Manuel-Navarrete et al., 2009; Viglizzo et al., 2011).
Among others, the challenges raised by the reduction of crop rotations in the short and medium term include (a) diminished soil fertility, (b) reduction and increased variability in crop yields, (c) ecological imbalance facilitating the expansion of weeds and infestations, and (d) increased greenhouse gas (GHG) emissions and vulnerability to climate change.
Diminished Soil Fertility
Agricultural crop rotation systems, such as those predominant in Argentina, show a lower capacity to sustain soil organic matter balances compared to rotation sequences that mix agricultural crops with pastures (Caride et al., 2012). In turn, within agricultural crop rotation systems, the maintenance of soil fertility (measured in terms of soil organic matter balance and nitrogen stock) is more likely to be achieved with no-tillage systems with rotations that combine different agricultural crops compared to monoculture of a single agricultural crop (Alvarez et al., 2014; Novelli et al., 2011).
Therefore, agricultural systems focused on soybean cultivation, like those used widely in Argentina, are more likely to maintain soil fertility when periodically incorporating alternative crops, such as grasses (maize and wheat) or, ideally, pasture.
However, it is worth noting that the focus of local agricultural rotation systems used in Argentina and other countries is to provide enough nutrients to maximize yields, regardless of whether those are sufficient to compensate for crop uptake (García et al., 2018). In contrast, in worldwide terms, the focus is on limiting the agricultural (and industrial) fixation of nitrogen and the inflow of phosphorus to oceans, so as to avoid transgressing planetary boundaries and causing environmental damage (Rockström et al., 2009).
Reduction and Increased Variability in Crop Yields
The evidence from long-term experimentation indicates that the monoculture of grasses (e.g., maize or wheat) or legumes (e.g., soybean) results in a reduction in crop yields compared with grass–legume crop rotations (Farmaha et al., 2016; Kelley et al., 2003). For example, in Argentina’s agriculture sector, it is considered a well-established fact that the rotation of grasses and legumes is more likely to achieve and sustain high yield levels compared to monoculture (Andrade et al., 2017).
However, the impact of crop rotation systems on yields varies across regions, and their relative profitability depends on exogenous variables such as input costs and taxation regimes. In addition, any medium-term positive impacts may not be considered relevant by producers with a short-term horizon, such as tenant farmers leasing land for a single season.
Ecological Imbalance
The introduction in the 1990s of genetically modified crops with resistance to herbicides resulted in very significant reductions in agricultural costs and greatly simplified agronomic management systems throughout the world.
However, the emergence of herbicide-resistant weeds, together with the mounting evidence regarding the negative impact of herbicide use on ecosystems, makes it necessary to develop alternative weed control strategies such as integrated weed management (IWM) (GROW, n.d.; Stokstad, 2019). Although these systems generally still feature chemical eradication solutions, they emphasize proper herbicide management practices and limit their use, favoring prevention strategies and non-chemical weed control tactics. These include mechanical weed management (e.g., mowing, harvest weed seed control, or robotic weeding), biological solutions (e.g., the introduction of bacteria, fungi, or insects that have a preference for certain weed species), and improved cultural practices (e.g., nutrient management, reduced row spacing, cover crops, or crop rotation) (Scursoni et al, 2019).
In addition, the stress on water resources induced by the monoculture or quasi-monoculture of agricultural crops can result both in increased water scarcity and in excessive soil saturation and more frequent flooding. In particular, monoculture affects, among other things, the physical fertility of soils. For this reason, water infiltration and soil storage capacity tend to be reduced. On the other hand, annual crops, and especially soybean, tend to consume less water. Studies in Argentina show an increase of groundwater levels as a consequence of this and, accordingly, a higher risk of flooding (Jobbagy et al., 2021; Nosetto et al., 2012).
Increased GHG Emissions and Vulnerability to Climate Change
Land use, including agricultural land use, is a significant net source of GHG emissions, contributing more than 20% of anthropogenic emissions of carbon dioxide (CO2), methane, and nitrous oxide combined as CO2 equivalents in 2007–2016. Therefore, climate change mitigation response options related to land use are a key element of most modeled viable climate scenarios (Intergovernmental Panel on Climate Change [IPCC], 2019).
In the specific context of Argentina’s agriculture soybean and maize are the main summer crops. According to Arrieta et al. (2018), the average yields per hectare for the 2012–2013 season were 2.62 tons of soybean and 6.02 tons of maize (Pampas region). In terms of energy efficiency and GHG emissions, on average 0.887 tons of soybean and 0.740 tons of maize were produced per gigajoule of energy, and 6.06 tons of soybean and 5.01 tons of maize were produced per ton of CO2 equivalent (CO2e) emitted (all regions). However, these figures only consider the emissions from cultivation and fertilization, without taking into account carbon changes in soil organic matter.
In overall terms, there are large differences in crop yields, GHG emission, and energy efficiency across regions in Argentina. These are likely driven by distinct climatic patterns, and in particular by differences in mean annual precipitation across regions (Calviño & Monzón, 2009). This is because agricultural production in Argentina is especially sensitive to climate and weather fluctuations. Therefore, agronomic adaptation strategies will play an essential role in increasing climate resilience. These include changing land and cropping practices, the development of improved crop varieties, and changing food consumption patterns (IPCC, 2019).
Crop Rotation and Climate Change Adaptation in Argentina’s Agriculture Sector
Data Sources
Agronomic diaries kept by producers in Argentina provide the opportunity to analyze plot-level data to study the evolution of agricultural practices, including crop rotation, and how they affect agricultural productivity and agriculture’s resilience against changing climate phenomena.
Only a small proportion of agricultural producers in Argentina belong to a technical association, but those who do typically belong to two main organizations: AACREA (Asociación Argentina de Consorcios Regionales de Experimentación Agrícola [Argentinian Association of Regional Consortia for Agricultural Experimentation]) and AAPRESID (Asociación Argentina de Productores en Siembra Directa [Argentinian Association of Direct Seeding Producers]).
AACREA is a national umbrella organization with approximately 2,000 members, composed of regional CREA associations (Consorcios Regionales de Experimentación Agrícola [Regional Consortia for Agricultural Experimentation]), which in turn are composed of local groups of 8–12 producers. Given the technical focus of the association, CREA members are encouraged to keep agricultural diaries so that plot-level records can be shared and analyzed to improve agricultural practices. Although agricultural diaries are generally kept in computer spreadsheet format, facilitating their analysis at the producer level, their format can vary across producers, regions, or even seasons, making aggregation difficult.
In 2018, AACREA partnered with the Inter-American Development Bank (IADB) and its Multilateral Investment Fund (now the IDB Lab) to collect, harmonize, and aggregate agricultural diaries from producers in 11 CREA regions. The consolidated database consists of anonymized plot-level information covering 19.2 million hectares over 19 agricultural seasons (1998–2016).
Figure 2 provides a map of all CREA regions, and Table 1 presents a list of CREA regions included in the 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016) (Asociación Argentina de Consorcios Regionales de Experimentación Agrícola, 2018).

Figure 2. Map of CREA agricultural regions in Argentina (2018).
Table 1. CREA Agricultural Regions Included in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016)
Code | Region | Code | Region | |
---|---|---|---|---|
NSF | Norte de Santa Fé | OAR | Oeste Arenoso | |
SFC | Santa Fé Centro | SUD | Sudeste | |
SSF | Sur de Santa Fé | MYS | Mar y Sierras | |
LIS | Litoral Sur | CEN | Centro | |
NBA | Norte de Buenos Aires | CBN | Córdoba Norte | |
OES | Oeste |
Source: Asociación Argentina de Consorcios Regionales de Experimentación Agrícola (2018).
The database was constructed using the Python programming language and PostgreSQL, an open-source relational database management system (RDBMS) that features embedded geographical information system (GIS) support through its PostGIS extension, allowing for easy storage and querying of geographical data. A custom extract–transform–load (ETL) tool was developed in order to automate the data entry and initial validation process from agricultural diaries in spreadsheet format, and the OpenRefine tool was used for basic database cleanup (e.g., imputing some types of missing values, eliminating duplicate values, and correcting spelling).
Table 2 provides a description of the seven main variable categories included in the 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016).
Table 2. Main Variable Categories in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016)
Category | Sample Variables in Category |
---|---|
Producer Data | Producer ID; local CREA group affiliation; regional CREA affiliation |
Plot Data | Plot ID; plot surface area; plot location coordinates |
Agronomic Data | Crop; fertilizer use; previous crop planted |
Seeding Data | Seeding date; soil saturation at seeding; phreatic water level |
Weather Data | Monthly rainfall; El Niño–Southern Oscillation phase |
Crop Data | Crop variety; plant success rate; crop yield |
Loss Data | Weed infestation; frost; hail; drought |
Source: Asociación Argentina de Consorcios Regionales de Experimentación Agrícola (2018).
According to the database, during the 1998–2016 period, approximately half of the land area cultivated by CREA producers was devoted to soybean, whereas the other half produced other crops, varying by region: (a) Southeastern regions (e.g., Sudeste or Mar y Sierras) typically cultivated soybean and feature sunflowers as the main non-grain summer rotation crop, with wheat/barley as the main grain crops used for winter rotation; (b) central-eastern regions (e.g., Oeste, Norte de Buenos Aires, Sur de Santa Fé, Santa Fé Centro, Norte de Santa Fé, or Litoral Sur) generally combined soybean cultivation with grain crops for both summer and winter rotations (maize and wheat/barley, respectively); and (c) central-western regions (e.g., Centro, Córdoba Norte, or Oeste Arenoso) tended to favor summer grain crops (maize) for rotation.
Figure 3 provides a summary of land area by crop in CREA regions included in the 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016). Figure 4 shows the distribution of land area for plots included in the database. Figure 5 graphs the distribution of use for the fertilizers most commonly employed by producers included in the database: nitrogen, sulfur, phosphorus pentoxide, and phosphorus.

Figure 3. Land area by crop and CREA region in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016).
Note: Average cultivated land area from 1998 to 2016 (complete data series starting in 2009). Other summer crops include mainly sunflower. Other winter crops include mainly barley.

Figure 4. Distribution of plot land area in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016) (all crops).

Figure 5. Distribution of fertilizer use by plot in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016) (all crops).
Impact of Climate on Agricultural Productivity
Rainfall in Argentina is highly influenced by the El Niño–Southern Oscillation (ENSO) phenomenon, and its two extreme phases, El Niño and La Niña, are strongly correlated with floods and droughts in areas of the country. The Oceanic Niño Index (ONI), computed by the National Oceanic and Atmospheric Administration of the United States, is used to identify the warm (El Niño) and cold (La Niña) phases of the ENSO phenomenon in the tropical Pacific. The ONI is computed using the running 3-month mean sea surface temperature anomaly for the Niño 3.4 region, and the ENSO extreme phases are defined as three consecutive overlapping 3-month periods above the +0.5ºC anomaly for the warm El Niño phase and below the –0.5ºC anomaly for the cold La Niña phase (Garbarini et al., 2016).
The impact of the ENSO phenomenon in the Argentinian Pampas is particularly strong in the spring and summer, coinciding with the growing season of the most important agricultural crops (Grimm et al., 2000). In neutral years, average precipitation levels during this time of the year tend to be very close to long-frequency climatological values (Garbarini et al., 2016). However, warm ENSO events (El Niño) tend to result in more spring and summer rainfall than in neutral years. Conversely, cold ENSO events (La Niña) tend to result in reduced spring and summer rainfall. Because the majority of agriculture in the Pampas is rainfed, ENSO phases strongly influence summer crop yields (Podestá et al., 1999; Travasso et al., 2003). In particular, soybean yields decrease in response to dry conditions that occur more frequently during the cold La Niña ENSO phase. In contrast, rainier conditions associated with the warm El Niño ENSO phase do not appear to significantly affect soybean yields (Podestá et al., 1999). ENSO phases are not considered to significantly impact yields of winter crops such as wheat and barley.
Impact of Crop Rotations on Agricultural Productivity
Impact of Short-Term Crop Rotation on Agricultural Yields
The impact of short-term crop rotation on agricultural productivity has been studied in many contexts using producer- and plot-level data, generally finding positive impacts. For example, Farmaha et al. (2016) find a positive impact of crop rotation on yields and resource efficiency in the U.S. corn belt when implementing maize–soybean rotations in high-yield irrigation applications.
An analogous analysis can be carried out based on the plot-level data compiled in the 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016) to estimate the impact of crop rotations in Argentina, where soybean has been the most prevalent monoculture or quasi-monoculture during the past few decades. By looking at the summer planting sequence by plot, it is possible to compare soybean yields when soybean was also cultivated during the previous planting season versus when a summer crop other than soybean was rotated in the plot or it was left fallow during the summer.
For simplicity purposes, only summer planting periods in the season are considered in the analysis (i.e., no distinction is made between a plot planted with a winter crop or left fallow during the winter). Also, for comparison robustness purposes, only the yield of early summer soybean crops (i.e., soybean crops planted early in the summer, after the plot was left fallow during the winter) are included in the analysis. This is the most common type of soybean crop planting in Argentina because late summer soybean crops (i.e., soybean crops planted later in the summer, typically following a winter crop such as wheat or barley) usually result in lower yields.
Figure 6 shows the distribution of early summer soybean yields for individual plots included in the 2018 AACREA–IADB Integrated Crop Rotation Database, relative to the average early summer soybean yield during the reference planting season in the local CREA group where the plot is located. A plot with an early summer soybean yield corresponding to the average early summer soybean yield in its local CREA group during the planting season of reference is considered to have a relative yield of 100. Comparing plot-level soybean yields with the average yield during the reference planting season in the relevant CREA local group avoids potential biases caused by season- and group-specific factors (e.g., average local soil characteristics and weather patterns, or differences in technology use and agronomic practices across seasons and CREA local groups). However, this type of comparison is not robust with respect to potential bias arising from plot- or producer-specific factors, such as those that may have affected producers’ choice to plant soybean this or the previous season (e.g., if more skilled producers are more likely to have rotated crops the previous season, or if crop rotation is more likely to be used in plots with poorer soil). Although this may be of no particular concern in the case of the 2018 AACREA–IADB Integrated Crop Rotation Database, given that the technological level and managerial skills of CREA members are likely to be very similar due to cooperation and information sharing. it may be an important consideration in other contexts.

Figure 6. Distribution of early summer soybean relative yields by short-term rotation intensity in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016).
Note: OS: other summer crop/summer fallow in Year 0 followed by early summer soybean in Year 1 (short-term rotation); SS: early or late summer soybean crop in Year 0 followed by early summer soybean in Year 1 (no short-term rotation). Distributional summary statistics: .
Two distributions are graphed in Figure 6: other–soybean (OS), the relative yield of early summer soybean crops in Year 1 for plots left fallow during the summer or where a summer crop different than soybean was planted in Year 0; and soybean–soybean (SS), the relative yield of early summer soybean crops in Year 1 for plots where a soybean (either early or late summer) crop was planted during the previous season. That is, OS and SS are the distributions of early summer soybean relative yields for plots with and without short-term soybean crop rotation, respectively. Both distributions feature comparable dispersions, but the distribution of early summer soybean relative yields with short-term crop rotation (OS) features a higher mean than the distribution without short-term crop rotation (SS). Although simple, this comparison already suggests a positive impact of short-term crop rotation on soybean yields.
The insights provided by the graphical comparison of relative yield distributions can be confirmed by fitting a linear regression model to the plot-level data in the 2018 AACREA–IADB Integrated Crop Rotation Database. Table 3 shows the estimated coefficients for the ordinary least squares (OLS) regression of early summer soybean relative yields y on a constant term c and an indicator variable that equals 1 if soybean was also planted the previous season (no short-term crop rotation), and 0 otherwise (short-term crop rotation).
Table 3. Linear Regression of Early Summer Soybean Relative Yields on Short-Term Soybean Rotation Intensity in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016)a
y = Early Summer Soybean Relative Yields | OLS Coefficient | p Value |
---|---|---|
Constant | 101.85 | 0.000 |
1 = no short-term rotation | –5.56 | 0.000 |
a The estimated regression model is , where denotes the early summer soybean relative yield of plot and is an indicator variable equal to 1 if soybean was also cultivated in plot the previous season, and 0 otherwise.
Source: Authors’ analysis based on data from Asociación Argentina de Consorcios Regionales de Experimentación Agrícola (2018).
The estimated constant term corresponds to the mean early summer soybean relative yield with crop rotation, 101.85. This indicates that, ceteris paribus, plots with short-term crop rotation indeed have mean early summer soybean yields that are 1.85 percentage points higher than the average early summer soybean yield for their respective local group during the reference planting season.
The estimated regression coefficient for , –5.56, is significant with 99% confidence. This indicates that, ceteris paribus, plots without short-term crop rotation have mean early summer soybean relative yields 5.56 percentage points lower than plots with short-term crop rotation.
Although this analysis continues to be very simple (e.g., ignoring any potential omitted variable bias due to plot-level characteristics), it suggests that the positive impact of short-term crop rotation on soybean yields can be quite large in agronomic terms.
Impact of Medium-Term Crop Rotation on Agricultural Yields
If sufficient historical data are available, the analysis of the impact of crop rotation on agricultural productivity can be expanded to estimate the impact of medium-term crop rotation over several planting seasons.
Figure 7 shows the distribution of early summer soybean relative yields in the reference (last in the sequence) planting season for plots in the 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016) for which data are available over periods of four or more consecutive planting seasons, ending in a season during which early summer soybean was planted. A plot with an early summer soybean yield corresponding to the average early summer soybean yield in its local CREA group during the last planting season (in the last consecutive period of four or more seasons) is considered to have an early summer soybean relative yield of 100.

Figure 7. Distribution of early summer soybean relative yields by medium-term rotation intensity in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016).
Note: Distribution of early summer soybean relative yields for [0, 20] = plots producing a summer crop other than soybean or left fallow during the summer in up to 20% of the last consecutive planting seasons; (20, 40] = plots producing a summer crop other than soybean or left fallow during the summer in more than 20% and up to 40% of the last consecutive planting seasons; (40, 60] = plots producing a summer crop other than soybean or left fallow during the summer in more than 40% and up to 60% of the last consecutive planting seasons; (60, 100] = plots producing a summer crop other than soybean or left fallow during the summer in more than 60% of the last consecutive planting seasons. Distributional summary statistics: .
Given that soybean was the most prevalent crop in Argentina during the time period covered in the 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016), a medium-term crop rotation intensity variable can be constructed based on the percentage of the last consecutive seasons during which a summer crop other than soybean was planted in a plot or it was left fallow during the summer. This intensity variable can then be converted into a multinomial categorical variable or into a set of binomial indicator variables corresponding to each medium-term crop rotation intensity interval.
For example, the indicator variable [0, 20] takes a value of 1 if a summer crop other than soybean was planted in the corresponding plot or it was left fallow during the summer in up to 20% of the last consecutive planting seasons, and a value of 0 otherwise (i.e., if early or late summer soybean was planted during more than 20% of the last consecutive planting seasons, such as in sequence of five planting seasons during which maize was planted instead of soybean during one season, or during which early or late summer soybean was planted in all five seasons). Analogously, the indicator variable (20, 40] takes a value of 1 if a summer crop other than soybean was planted in the corresponding plot or it was left fallow during the summer in more than 20% of the last consecutive planting seasons up to 40%, and a value of 0 otherwise (i.e., if early or late summer soybean was planted during more than 40% of the last consecutive planting seasons, such as in sequence of five planting seasons during which maize was planted instead of soybean during two seasons). The length of the intervals could be chosen arbitrarily, but narrow intervals may cause data sparsity problems, whereas wide intervals may be uninformative. In the case of Argentina, the interval [0, 20] captures plots devoted to soybean cultivation with no or limited long-term crop rotation. The intervals (20, 40] and (40, 60] capture the two crop rotation intensities that are generally considered feasible by soybean producers in Argentina. Finally, the interval (60, 100] captures plots generally devoted to the cultivation of other crops, with limited rotation of soybean cultivation.
As in the analysis of short-term rotations, for simplicity purposes, only summer planting periods in the season are considered in the analysis (i.e., no distinction is made between a plot planted with a winter crop or left fallow during the winter). Similarly, for comparison robustness purposes, only the yield of early summer soybean crops (i.e., soybean crops planted early in the summer, after the plot was left fallow during the winter) are included in the analysis—that is, all the sequences of four or more consecutive seasons included in the analysis end in a season during which early summer soybean was planted.
Four distributions are graphed in Figure 7, corresponding to plots in each of the four defined medium-term crop rotation intensity intervals: [0, 20], (20, 40], (40, 60], and (60, 100]. As in the case of short-term crop rotation intensity, all four early summer soybean relative yield distributions feature comparable dispersions, but mean early summer soybean productivity during the season of reference increases with the intensity of medium-term crop rotation during the previous seasons.
Analogous to the case of short-term crop rotations, the insights provided by the graphical comparison of relative soybean yield distributions can be confirmed by fitting a linear regression model to the plot-level data in the 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016).
Table 4 shows the estimated coefficients for the OLS regression of early summer soybean relative yields y on a constant term c and the crop rotation intensity variable r (i.e., the percentage of the last four or more consecutive seasons for which data are available during which a summer crop other than soybean was planted in the plot or it was left to fallow during the summer). The estimated regression coefficient for r, 0.22, is significantly different from 0 with 99% confidence. This indicates that, ceteris paribus, an additional season of crop rotation in four seasons is estimated to result in early summer soybean yields that are 5.5 percentage points higher than the average early summer soybean yield for the respective local group during the reference planting season.
Table 4. Linear Regression of Early Summer Soybean Relative Yields on Medium-Term Rotation Intensity in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016)a
y = Early Summer Soybean Relative Yields | OLS Coefficient | p Value |
---|---|---|
Constant | 91.86 | 0.000 |
Medium-term rotation intensity (%) | 0.22 | 0.000 |
Constant | 93.97 | 0.000 |
1 = (20,40] medium-term rotation | 3.72 | 0.010 |
1 = (40,60] medium-term rotation | 10.62 | 0.000 |
1 = (60,100] medium-term rotation | 12.18 | 0.000 |
a The first estimated regression model is , where denotes the early summer soybean relative yield of plot and denotes its medium-term crop rotation intensity (i.e., the percentage of the last four or more consecutive seasons for which data are available in which a summer crop other than soybean was planted in plot , or it was left fallow during the summer). The second estimated regression model is , where denotes the early summer soybean relative yield of plot and , and are binomial indicator variables that equal 1 if plot belongs to the corresponding medium-term intensity crop rotation interval, and 0 otherwise.
Source: Authors’ analysis based on data from Asociación Argentina de Consorcios Regionales de Experimentación Agrícola (2018).
This result is very similar to the estimated impact of short-term crop rotation previously obtained by comparing only two seasons. However, the estimated impact of crop rotation on early summer soybean relative yields is even larger when the possibility to rotate during several seasons in the medium term is taken into account: Table 4 also shows the estimated coefficients for the OLS regression of early summer soybean relative yields y on a constant term c and binomial indicator variables , and , which equal 1 if a plot belongs to the corresponding medium-term intensity crop rotation interval, and 0 otherwise. The results indicate that, ceteris paribus, the mean early summer soybean relative yields for plots in which a summer crop other than soybean was planted or which were left fallow during the summer in more than 20% and up to 40% of the last consecutive seasons (i.e., 60% or more, but less than 80% soybean cultivation, moderate crop rotation) are 3.72 percentage points higher than the mean early summer soybean relative yields for plots in which a summer crop other than soybean was planted or which were left fallow during the summer in less than 20% of the last consecutive seasons (i.e., 80% or more soybean cultivation, low or no crop rotation). But, ceteris paribus, the mean early summer soybean relative yields for plots in which a summer crop other than soybean was planted or which were left fallow during the summer in more than 40% and up to 60% of the last consecutive seasons (i.e., 40% or more, but less than 60% soybean cultivation, balanced crop rotation) are estimated to be 10.62 percentage points higher, whereas the mean early summer soybean relative yields for plots in which a summer crop other than soybean was planted or which were left fallow during the summer in more than 60% of the last consecutive seasons (i.e., less than 40% soybean cultivation, predominance of other crops) are estimated to be 12.18 percentage points higher.
This analysis continues to rely on a very simple identification strategy, but it suggests that the positive productivity impact of balanced crop rotation over many seasons in the medium term can be even larger than that of short-term crop rotation over just two seasons.
Panel Data Analysis
So far, the analysis of the impact of crop rotation on soybean productivity has only attempted to avoid potential biases caused by season- and group-specific factors, by examining (early summer) soybean yields relative to the average productivity of plots in the relevant local group during the reference planting season. However, the panel structure of the 2018 AACREA–IADB Integrated Crop Rotation database allows the use of more sophisticated analysis strategies to address potential omitted variable bias arising from non-observable time-invariant plot characteristics (e.g., non-observable soil quality or non-observable producer ability).
Figure 8 summarizes the results of estimating the impact of medium-term crop rotations on relative soybean yields using three panel data models, which leverage the fact that the data set includes observations corresponding to medium-term planting sequences of the same plots at different points in time: (a) A random effects model assumes that time-invariant plot-specific characteristics (i.e., effects) are random variables (orthogonal to any time-varying explanatory variables included in the analysis); (b) a within model assumes that time-invariant plot-specific characteristics are fixed (i.e., non-random) quantities, and it controls for this by examining the deviations from the overall plot-level mean values over time; and (c) a first-difference model also assumes that time-invariant specific characteristics are fixed, but it controls for this by instead examining the evolution of plot-level values from one time period (i.e., planting sequence) to the next. In addition, a pooling model is also considered: This simply pools all observations together in the estimation—even observations that correspond to the same plot at different points in time.

Figure 8. Panel data estimates of impact of medium-term rotation intensity on early summer soybean relative yields in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016).
Note: The y-axis represents the estimated increase in mean early summer soybean relative yields during the last season in a medium-term planting sequence (four or more consecutive planting seasons). The pooling model ignores the panel structure of the data. The random effects model assumes that time-invariant plot-specific effects are random variables. The within model is the standard fixed effects model, where time-invariant plot-specific characteristics are assumed to be fixed quantities, and unbiasedness is achieved by controlling for deviations from the overall plot-level mean values over time. Finally, the first-difference (fd) model also assumes fixed effects, but unbiasedness is achieved by instead controlling for the evolution of plot-level values from one planting sequence to the next.
Given that the pooling model ignores the panel structure of the data, its results are comparable to those obtained so far. The results of estimating a random effects model are also comparable, but this estimation will only be consistent if time-invariant plot characteristics are indeed orthogonal to any time-varying explanatory variables. Otherwise, the results of the within or first-difference fixed effects models suggest that the positive impact of medium-term crop rotation on productivity may be even larger than those estimated so far, with balanced (40–60%) medium-term rotation resulting in an estimated 15–20 percentage point increase in mean (early summer) soybean relative yields compared to no or low medium-term crop rotation intensity (0–20%).
For a detailed discussion of econometric panel data models, see Wooldridge (2010).
Impact of Crop Rotations on Climate Adaptation
Impact of Short-Term Crop Rotation on Overall Agricultural Resilience
Figure 9 graphs early summer soybean absolute yields with and without short-term crop rotation for plots included in the 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016), grouped by CREA local group. Each observation represents the average early summer soybean absolute yield (kg/ha) of plots in a CREA local group for which early summer soybean cultivation in the reference planting season was preceded by soybean (either early or late summer) or by a summer crop other than soybean or summer fallow in the planting season sequence. That is, each observation corresponds to the average early summer soybean yield during a single planting season of two sets of plots managed by producers in a CREA local group: those for which early summer soybean cultivation during the corresponding planting season was preceded by soybean cultivation (either early or late summer) during the previous planting season (x-axis, average CREA local group early summer soybean absolute yield without short-term crop rotation) and those for which soybean cultivation during the same planting season was preceded by cultivation of a summer crop other than soybean or summer fallow during the previous planting season (y-axis, average CREA local group early summer soybean absolute yield with short-term crop rotation).

Figure 9. Average soybean yields with and without short-term crop rotation, by CREA local group and planting season in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016).
y = Average Early Summer Soybean Absolute Yields (kg/ha) by CREA Local Group and Planting Season | OLS Coefficient | p Value |
---|---|---|
Constant | 142.69 (75.96) | 0.061 |
1 = short-term rotation | 0.90 (0.02) | 0.000 |
Note: Analysis excludes observations (unique local group-planting season combinations) comprising less than 20 plots in either the with or the without crop rotation sets. Linear regression is unweighted by observation size (i.e., the linear regression does not take into account the number of plots included in each observation). The estimated regression model is , where denotes the average early summer soybean absolute yield of plot set in CREA local group during planting season , and is an indicator variable equal to 1 if plot set in CREA local group during planting season corresponds to a short-term rotation (other-soybean) planting sequence, and 0 otherwise. A -test rejects the null hypothesis of equality of average absolute yields with and without short-term crop rotation ( value = 0.000).
Therefore, observations in Figure 9 only capture the average absolute yield of plots for which early summer soybean was cultivated during the reference season. Also, each observation represents a unique combination of planting season and CREA local group: Several observations may correspond to the same CREA local group, provided that they correspond to distinct planting seasons. In addition, the diameter of each observation in Figure 9 is proportional to the number of plots included in that observation, relative to the total number of plots included in the analysis. That is, a larger diameter means that the local CREA group corresponding to that observation cultivated soybean in more plots during the corresponding season.
The linear regression analysis of average early summer soybean absolute yields suggests that short-term crop rotation, from one season to the next, had a positive but uneven impact on soybean yields in Argentina. In particular, the estimated y-intercept of the linear regression line, 142.69, is greater than 0, but its estimated slope, 0.90, is less than 1. This indicates that although the overall short-term impact of crop rotation is positive, ceteris paribus, the average impact of short-term crop rotation is smaller for higher yield season–local group combinations (observations up and to the right in Figure 9). In other words, the average positive impact of short-term crop rotation appears to be larger in less favorable season–local group contexts where yields are lower (observations down and to the left in Figure 9): Whereas the t-test null hypothesis of equality in average yields with and without short-term crop rotation is rejected with 99% confidence, at 95% confidence the difference is not significant for season–local group combinations with average yields in excess of 4,000 kg/ha or higher (i.e., starting approximately at that point, the linear regression line’s 95% confidence area overlaps with the 1:1 yield line).
This basic result suggests that the positive impact of short-term crop rotation on soybean yields may be larger in adverse conditions so that crop rotation may contribute to increase the overall resilience of Argentina’s agricultural sector.
Impact of Short-Term Crop Rotation on Agricultural Climate Resilience
What about climate resilience in particular? In order to address this more specific question, the analysis of the overall impact of short-term crop rotation on relative soybean yields can be replicated for different climate scenarios, such as the phases of the ENSO phenomenon. In the context of Argentina, the El Niño phase of the ENSO phenomenon is generally considered to bring more favorable (summer) agricultural conditions than the neutral ENSO phase. Conversely, the La Niña ENSO phase is generally considered to bring less favorable (summer) agricultural conditions than the neutral ENSO phase. Therefore, the analysis of the impact of short-term and medium-term crop rotation during distinct ENSO phases allows the quantification of the effectiveness of crop rotation as a climate adaptation tool.
Analogous to Figure 6, Figure 10 graphs the distributions OS (the early summer relative soybean yields for plots where a summer crop different than soybean was planted or which were left fallow during the summer in the previous season—that is, with short-term crop rotation) and SS (the early summer soybean relative yields for plots where soybean was planted during the previous season—that is, without short-term crop rotation). However, Figure 10 further disaggregates the distributions by ENSO phase, graphing separately the early summer soybean relative yields for planting seasons during the neutral, La Niña, and El Niño ENSO phases.

Figure 10. Distribution of early summer soybean relative yields by short-term rotation intensity and ENSO phase in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016).
Note: OS: other summer crop/summer fallow in Year 0 followed by early summer soybean in Year 1 (short-term rotation); SS: early or late summer soybean crop in Year 0 followed by early summer soybean in Year 1 (no short-term rotation). Distributional summary statistics: .
As it can be observed in Figure 10, the distribution of early summer relative soybean yields features more dispersion during the neutral and La Niña ENSO phases compared to the El Niño ENSO phase. Most important, short-term crop rotation appears to have a large positive distributional impact during the La Niña ENSO phase compared to the neutral and El Niño ENSO phases.
Again, the insights provided by the graphical comparison of early summer soybean relative yield distributions can be confirmed by fitting a linear regression model to the plot-level data in the 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016). Analogous to Table 3, Table 5 shows the estimated coefficients for the OLS regression of early summer soybean relative yields y on a constant term c and an indicator variable that equals 1 if soybean was also planted the previous season (no short-term crop rotation) and 0 otherwise (short-term crop rotation). However, in Table 5, the estimation results are further disaggregated by ENSO phase.
Table 5. Linear Regression of Early Summer Soybean Relative Yields on Short-Term Rotation Intensity by ENSO Phase in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016)a
y = Early Summer Soybean Relative Yields | OLS Coefficient | p Value |
---|---|---|
Neutral ENSO phase | ||
Constant | 101.34 | 0.000 |
1 = no short-term rotation | –6.26 | 0.000 |
La Niña ENSO phase | ||
Constant | 103.54 | 0.000 |
1 = no short-term rotation | –8.18 | 0.000 |
El Niño ENSO phase | ||
Constant | 101.72 | 0.000 |
1 = no short-term rotation | –4.46 | 0.000 |
a The estimated regression model is , where denotes the early summer soybean relative yield of plot and is an indicator variable equal to 1 if soybean was also cultivated in plot the previous season, and 0 otherwise.
Source: Authors’ analysis based on data from Asociación Argentina de Consorcios Regionales de Experimentación Agrícola (2018).
As it can be observed, ceteris paribus, the estimated mean positive impact of short-term crop rotation on early summer soybean relative yields is almost double during the less favorable La Niña ENSO phase compared to the more favorable El Niño ENSO phase (8.18 versus 4.46 percentage points, respectively).
This still basic comparison suggests that the positive impact of short-term crop rotation on soybean yields may be considerably larger in unfavorable climate conditions, such as the drier La Niña ENSO phase in Argentina.
Impact of Medium-Term Crop Rotation on Agricultural Climate Resilience
The analysis of the overall impact of medium-term crop rotation on early summer soybean relative yields can also be replicated for the different phases of the ENSO phenomenon.
Analogous to Figure 7, Figure 11 shows the distribution of early summer soybean relative yields by medium-term crop rotation intensity intervals. However, analogous to Figure 10, Figure 11 further disaggregates the distributions, graphing separately the relative soybean yields for planting seasons during the neutral, La Niña, and El Niño ENSO phases.

Figure 11. Distribution of early summer soybean relative yields by medium-term rotation intensity and ENSO phase in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016).
Note: Distribution of early summer soybean relative yields for [0, 20] = plots producing a summer crop other than soybean or left fallow during the summer in up to 20% of the last consecutive planting seasons; (20, 40] = plots producing a summer crop other than soybean or left fallow during the summer in more than 20% and up to 40% of the last consecutive planting seasons; (40, 60] = plots producing a summer crop other than soybean or left fallow during the summer in more than 40% and up to 60% of the last consecutive planting seasons; (60, 100] = plots producing a summer crop other than soybean or left fallow during the summer in more than 60% of the last consecutive planting seasons. Distributional summary statistics: ; ; .
As it can be observed in Figure 11, the distribution of early summer soybean relative yields for all medium-term crop rotation intervals again features more dispersion during the neutral and La Niña ENSO phases compared to the El Niño ENSO phase. Also, as in the case of short-term crop rotation, medium-term crop rotation appears to have a large positive distributional impact during the La Niña ENSO phase compared to the neutral and El Niño ENSO phases. This is particularly noticeable for moderate (20–40%) and balanced (40–60%) medium-term rotation intensities.
Analogous to Table 4, Table 6 shows the estimated coefficients for the OLS regression of early summer soybean relative yields y on a constant term c and (a) the crop rotation intensity variable r (i.e., the percentage of the last four or more consecutive seasons for which data are available in which a summer crop other than soybean was planted in the plot or it was left fallow during the summer) and (b) binomial indicator variables , and , which equal 1 if a plot belongs to the corresponding medium-term intensity crop rotation interval and 0 otherwise. However, in Table 6 the estimation results are further disaggregated by ENSO phase.
Table 6. Linear Regression of Early Summer Soybean Relative Yields on Medium-Term Rotation Intensity by ENSO Phase in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016)
y = Relative Soybean Yields | OLS Coefficient | p Value |
---|---|---|
Neutral ENSO phase | ||
Constant | 88.56 | 0.000 |
Medium-term rotation intensity (%) | 0.27 | 0.000 |
Constant | 90.43 | 0.000 |
1 = (20, 40] medium-term rotation | 5.81 | 0.010 |
1 = (40, 60] medium-term rotation | 14.61 | 0.000 |
1 = (60, 100] medium-term rotation | 15.11 | 0.000 |
La Niña ENSO phase | ||
Constant | 90.70 | 0.000 |
Medium-term rotation intensity (%) | 0.35 | 0.000 |
Constant | 91.75 | 0.000 |
1 = (20, 40] medium-term rotation | 10.05 | 0.013 |
1 = (40, 60] medium-term rotation | 21.18 | 0.000 |
1 = (60, 100] medium-term rotation | 20.80 | 0.000 |
El Niño ENSO phase | ||
Constant | 97.80 | 0.000 |
Medium-term rotation intensity (%) | 0.06 | 0.044 |
Constant | 99.1 | 0.000 |
1 = (20, 40] medium-term rotation | –0.42 | 0.835 |
1 = (40, 60] medium-term rotation | 2.23 | 0.274 |
1 = (60, 100] medium-term rotation | 3.99 | 0.073 |
Note: The first estimated regression model is , where denotes the early summer soybean relative yield of plot and denotes its medium-term crop rotation intensity (i.e., the percentage of the last four or more consecutive seasons for which data are available in which a summer crop other than soybean was planted in plot , or it was left fallow during the summer). The second estimated regression model is , where denotes the early summer soybean relative yield of plot and , and are binomial indicator variables that equal 1 if plot belongs to the corresponding medium-term intensity crop rotation interval, and 0 otherwise.
Source: Authors’ analysis based on data from Asociación Argentina de Consorcios Regionales de Experimentación Agrícola (2018).
The estimated regression coefficients indicate that although moderate medium-term crop rotation (20–40%) did not have a significant impact during El Niño ENSO phases, it did have a large significant impact during neutral and La Niña ENSO phases: Compared to no or low medium-term crop rotation intensity, ceteris paribus, the estimated mean early summer soybean relative yields with moderate medium-term crop rotation during neutral and La Niña ENSO phases were 5.81 and 10.05 percentage points higher, respectively. The impact is even more pronounced for balanced medium-term crop rotation (40–60%): Compared to no or low medium-term crop rotation, ceteris paribus, the estimated mean early summer soybean relative yields with balanced medium-term crop rotation during neutral and La Niña ENSO phases were 14.61 and 21.18 percentage points higher, respectively. However, the impact of balanced medium-term crop rotation during El Niño ENSO phases was estimated at just 2.23 percentage points compared to no or low medium-term crop rotation. Also, the impact of increased medium-term crop rotation intensity on mean early summer soybean relative yields appears to taper off as the cultivation of other crops becomes dominant (60–100% medium-term crop rotation).
Although with the available data precision in the estimation is low, resulting in wide confidence intervals for the coefficients, these results suggest that moderate (20–40%) and, especially, balanced (40–60%) medium-term crop rotation may be a very effective tool to increase climate resilience in Argentina’s agriculture sector, softening the impact of unfavorable climate conditions such as the drier Neutral and La Niña ENSO phases.
Impact of Crop Rotations on Climate Change Mitigation and the Environment
The discussion has so far focused on the impact of short-term and medium-term crop rotations on agricultural yields and also on their potential to facilitate climate adaptation through increased agricultural resilience. However, the distinct agronomic practices associated with crop rotation can also have an environmental impact.
Table 7 showcases five agronomic models that can be used to estimate the impact of agricultural production on several environmental aspects: (a) GHG emissions, (b) input energy efficiency, (c) soil phosphorus balance, (d) soil organic carbon balance, and (e) pesticide effective lethal dose.
Table 7. Selected Climate Change Mitigation and Environmental Agronomic Indicators
Indicator | Model | Units | Goal |
---|---|---|---|
GHG emissions | CO2e kg/ha | Minimize | |
Input energy efficiency | Minimize energy ratio of products and inputs | Mcal (products)/Mcal (inputs) | Maximize |
Soil phosphorus balance | P kg/ha | Maximize | |
Soil organic carbon balance | AMG (Andriulo et al., 1999) | SOC kg/ha | Maximize |
Pesticide effective lethal dose | Pesticide median lethal dose (LD50) | Pesticide LD50/ha | Minimize |
CO2e, CO2 equivalent; SOC, soil organic carbon.
Figure 12 summarizes the results of applying these five agronomic models to the 2018 AACREA–IADB Integrated Crop Rotation Database. Spider charts are employed in order to provide a quick visual comparison of environmental outcomes by rotation intensity for each CREA region included in the Asociación Argentina de Consorcios Regionales de Experimentación Agrícola (2018) database (100 = positive maximization/minimization outcome achieved). That is, data points closer to the center of the spider graph correspond to worse environmental outcomes, whereas data points closer to the edge of the spider graph correspond to better environmental outcomes.

Figure 12. Estimated achievement percentage of mitigation and environmental goals, by CREA agricultural region and crop rotation intensity in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016).
Note: (a) GHG emissions, (b) input energy efficiency, (c) soil phosphorus balance, (d) soil organic carbon balance, and (e) pesticide effective lethal dose. 100 = 100% goal achievement. See Table 1 for region codes. See Table 7 for model details.
An increase in crop rotation intensity is consistently associated with lower GHG emissions in all regions (Figure 12a). This is primarily due to the improvement in organic carbon balances in the soil associated with an increase in rotation intensity (Figure 12d), which more than compensates for the higher GHG emissions associated with fertilizer and pesticide use in the cultivation of grasses such as maize, compared to the cultivation of soybean.
However, the estimated success of crop rotation in achieving better environmental outcomes in the remaining three dimensions is more mixed because the increase in grass cultivation associated with higher crop rotation intensities is associated with worse phosphorus balance in the soil (Figure 12c) and increased pesticide use in some regions (Figure 12e). In addition, an increase in crop rotation intensity is generally associated with lower input energy efficiency (Figure 12b) due to the higher input requirements of grass cultivation compared to soybean crops.
These results suggest that the rotation of other crops in soybean cultivation in Argentina has the potential to contribute to a reduction in GHG emissions and improved soil organic carbon balances but that in some regions this would occur at the expense of increased soil phosphorus imbalances, increased pesticide use, and/or lower input energy efficiency.
Impact of Crop Rotations on Profitability
Finally, it is important to consider how the implementation of higher crop rotation intensities may impact the profitability of agricultural producers. Profitability will determine whether higher crop rotation intensities are economically sustainable and/or whether there may be a case for public support of crop rotation in view of the potential positive externalities.
However, the profitability of crop rotation is likely to vary across regions and economic contexts. For example, the taxation of agricultural exports in Argentina has a major impact on profitability, and tax rates and policies have varied significantly over time. Similarly, commodity prices, including those of agricultural products, can experience high volatility and cyclicality.
Table 8 describes three sample scenarios: (E1) export taxes for all crops but favorable treatment of maize and wheat, compared to soybean (Argentina’s context before December 2015 tax export reform); (E2) lower export taxes for soybean, but no export taxes for maize or wheat (Argentina’s context after December 2015 tax export reform); and (E3) no export taxes for any crops (hypothetical free market context). In addition, for each scenario, two subscenarios are considered: (a) ignores benefits from medium-term crop rotation (i.e., crop rotation intensity only affects the regional crop mix but not crop yields or fertilizer use, which are always considered to be the average for the CREA region) and (b) includes benefits from medium-term crop rotation (i.e., crop rotation intensity affects the crop mix, crop yields, and fertilizer use, based on actual data from the Asociación Argentina de Consorcios Regionales de Experimentación Agrícola [2018] database for each medium-term crop rotation intensity interval and CREA region).
Table 8. Description of Economic Scenarios for Profitability Estimation
Scenario | Prices (Net of Taxes) | Export Taxes | Crop Yields | Fertilizer Use | Input Costs | |
---|---|---|---|---|---|---|
E1 | a | FAS 2005–2014 (All) | 35% (S); 25% (M, W) | Region average | Region average | Region average |
b | “ | “ | Actual | Actual | Region average | |
E2 | a | FAS 2005–2014 (S) FOB 2005–2014 (M, W) | 30% (S); 0% (M, W) | Region average | Region average | Region average |
b | “ | “ | Actual | Actual | Region average | |
E3 | a | FOB 2005–2014 (All) | 0% (all) | Region average | Region average | Region average |
b | “ | “ | Actual | Actual | Region average |
M, maize; S, soybean; W, wheat. Actual, from Asociación Argentina de Consorcios Regionales de Experimentación Agrícola (2018) database (for each crop rotation intensity interval); region average, from complementary data.
The profitability of agricultural production in each of these three scenarios is then estimated for the four medium-term crop rotation intensity intervals defined (i.e., 0–20%, 20–40%, 40–60%, and 60–100%), by CREA region. Profitability is defined as gross margin/direct costs. The gross margin is defined as revenues net of commercial expenses and direct costs. In turn, direct costs include labor, inputs, insurance, and land leases.
As shown in Figure 13, the preferential tax treatment of non-soybean crops (maize and wheat) in scenario E2 results in higher relative profitability of increased crop rotation intensities compared to scenarios E1 or E3, where the more symmetric tax treatment of soybean and non-soybean crops lowers the profitability of increased crop rotation intensities. Also, the estimated profitability of increased crop rotation intensities is generally higher when their positive impact on yields and fertilizer use is considered, again particularly in scenario E2, in which non-soybean crops enjoy a preferential tax treatment. Finally, estimates for agricultural profitability vary substantially from region to region. For example, the estimated profitability of increased crop rotation intensities appears to be generally higher in regions in which non-soybean crops are typically rotated in during winter (e.g., Mar y Sierras or Sudeste) compared to regions in which non-soybean crops are typically rotated in during summer (e.g., Centro or Oeste Arenoso).

Figure 13. Estimated mean agricultural profitability, by CREA agricultural region and crop rotation intensity in 2018 AACREA–IADB Integrated Crop Rotation Database (1998–2016).
Note: Profitability is defined as gross margin/direct costs. See Table 1 for region codes. See Table 8 for scenario definitions.
These results suggest that although in many cases balanced crop rotations will be the most profitable option for agricultural producers, in certain regions and/or under certain tax regimes, balanced crop rotations may decrease profitability and additional public sector support may be needed in order to incentivize their implementation.
Potential Avenues for Future Research
The basic results showcased in this article suggest that crop rotation, and particularly balanced crop rotation, may result in very large positive impacts on soybean yields, particularly in unfavorable climatic conditions such as those experienced during the La Niña ENSO phase in Argentina. In addition to this positive impact on agricultural productivity and climate adaptation, in some contexts crop rotation may also contribute to the reduction of GHG emissions, increased input energy efficiency, and improved environmental outcomes.
These results are based on sample analyses of the 2018 AACREA–IADB Integrated Crop Rotation Database, which compiled and harmonized the information from agricultural diaries kept by CREA members in Argentina from 1998 to 2016. This new consolidated data set has replaced previous regional templates, and it is expected to continue to be periodically expanded with new information. Therefore, it offers potential opportunities for further research on the impact of crop rotation on climate adaptation and on other topics in agricultural and environmental economics.
The sample analyses showcased in this article can be augmented with more sophisticated identification strategies that address some of the robustness questions discussed (e.g., omitted variable bias from unobserved plot- and producer-level characteristics) and/or that shed more light on the specifics of how crop rotation can be most effectively used to facilitate agricultural climate adaptation (e.g., by further exploring how the basic results are influenced by regional differences in soil characteristics and weather patterns, or which planting sequences are optimal in different contexts).
In addition, further research is needed to refine complementary estimation models, such as the climate mitigation, environmental, and profitability models discussed in this article. For example, multidimensional models can be developed to simultaneously estimate the impact of crop rotation on several key aspects, such as yields and environmental benefits. This would allow for a better understanding of the trade-offs involved, the facilitation of communication to agricultural producers, and the identification of opportunities for public policy interventions (e.g., to incentivize longer land leases that attenuate the impact of the shorter time horizons over which land lessees optimize compared to land owners). Similarly, additional scenarios can be developed, for example, to model the impact of a transition to a low- or zero-emission economy.
Some of these avenues of research are already open, and others will become feasible as the AACREA–IADB Integrated Crop Rotation Database gains temporal and cross-sectional depth. The insights obtained could prove pivotal in facilitating and accelerating climate change adaptation in Argentina’s agriculture sector.
However, ultimately there is a limit as to how far agriculture can adapt to climate change, and a political will to reduce the impact of fossil fuels is essential for long-term food security (Anderson et al., 2020).
Acknowledgments
We thank CREA regional consortia and their members for sharing the agricultural experimentation data on which this article relies. The insightful comments of two anonymous reviewers are also gratefully acknowledged. Funding and institutional support for this research were generously provided by the Inter-American Development Bank’s Multilateral Investment Fund (now the IDB Lab) and AACREA. However, the opinions represented in this article are those of the authors only, and they do not purport to reflect the views of any institution.
Further Reading
- Andrade, F. H. (2020). Los desafíos de la agricultura global. Ediciones INTA.
- Andrade, J. F., & Satorre, E. H. (2015). Single and double crop systems in the Argentine Pampas: Environmental determinants of annual grain yield. Field Crops Research, 177, 137–147. https://doi.org/10.1016/j.fcr.2015.03.008
- Angrist, J. D., & Pischke, J.‑S. (2009). Mostly harmless econometrics: An empiricist’s companion. Princeton University Press.
- Anwar, M. R., Liu, D. L., Macadam, I., & Kelly, G. (2013). Adapting agriculture to climate change: A review. Theoretical and Applied Climatology, 113(1), 225–245.
- Asociación Argentina de Consorcios Regionales de Experimentación Agrícola. (2018). Base de datos integrada de rotaciones agrícolas en Argentina (Convenio AACREA-BID-FOMIN) (v. 2018–02) [Data set].
- Asociación Argentina de Consorcios Regionales de Experimentación Agrícola. (2021). Sustainability indicators system CREA [Data set].
- D’Acunto, L., Andrade, J. F., Poggio, S. L., & Semmartin, M. (2018). Diversifying crop rotation increased metabolic soil diversity and activity of the microbial community. Agriculture, Ecosystems & Environment, 257, 159–164.
- de Abelleyra, D., & Verón, S. (2020). Crop rotations in the Rolling Pampas: Characterization, spatial pattern and its potential controls. Remote Sensing Applications, 18, 100320.
- Farmaha, B. S., Eskridge, K. M., Cassman, K. G., Specht, J. E., Yang, H., & Grassini, P. (2016). Rotation impact on on-farm yield and input-use efficiency in high-yield irrigated maize–soybean system. Agronomy Journal, 108(6), 2313–2321.
- Food and Agriculture Organization. (2004). Fertilizer use by crop in Argentina.
- Hora, R. (2012). La evolución del sector agroexportador argentino en el largo plazo, 1880–2010. Historia Agraria, 58, 145–181.
- Inter-American Development Bank. (2021). Agricultural policies monitoring system (BID Agrimonitor) [Data set].
- Intergovernmental Panel on Climate Change. (2019). Climate change and land: An IPCC special report on climate change, desertification, land degradation, sustainable land management, food security, and greenhouse gas fluxes in territorial ecosystems. In P. R. Shukla, E. Calvo Buendia, V. Masson-Delmotte, H.‑O. Pörtner, D. C. Roberts, P. Zhai, R. Slade, S. Connors, R. van Diemen, M. Ferrat, E. Haughey, S. Luz, S. Neogi, M. Pathak, J. Petzold, J. Portugal Pereira, P. Vyas, E. Huntley, K. Kissick, . . . J. Malley (Eds.). In press.
- Jobbagy, E. G., Aguiar, S., Piñeiro, G., & Garibaldi, L. A. (2021). Impronta ambiental de la agricultura de granos en Argentina: Revisando desafíos propios y ajenos. Ciencia Hoy, 29(173), 35–44.
- Kuntashula, E., Chabala, L. M., & Mulenga, B. P. (2014). Impact of minimum tillage and crop rotation as climate change adaptation strategies on farmer welfare in smallholder farming systems of Zambia. Journal of Sustainable Development, 7(4), 95.
- Liu, H., Xiong, W., Li, Y.‑C., & Yang, D. (2017). Advances of impacts and adaptation of climate change on crop rotations in China. Chinese Journal of Agrometeorology, 38(10), 613–631.
- Manuel-Navarrete, D., Gallopín, G. C., Blanco, M., Díaz-Zorita, M., Ferraro, D. O., Herzer, H., Laterra, P., Murmis, M. R., Podestá, G. P., Rabinovich, J., Satorre, E. H., Torres, F., & Viglizzo, E. F. (2009). Multi-causal and integrated assessment of sustainability: The case of agriculturization in the Argentine Pampas. Environment, Development and Sustainability, 11(3), 621–638.
- Ministerio de Ambiente y Desarrollo Sostenible. (2015). Tercera Comunicación Nacional de Cambio Climático. Dirección Nacional de Cambio Climático, Ministerio de Ambiente y Desarrollo Sostenible.
- National Oceanic and Atmospheric Administration. (2021). El Niño–Southern Oscillation (ENSO).
- Nosetto, M. D., Jobbágy, E. G., Brizuela, A. B., & Jackson, R. B. (2012). The hydrologic consequences of land cover change in central Argentina. Agriculture, Ecosystems & Environment, 154, 2–11.
- Pengue, W. (2005). Transgenic crops in Argentina: The ecological and social debt. Bulletin of Science, Technology & Society, 25, 314–322.
- Rockström, J., Steffen, W., Noone, K., Persson, Å., Chapin, F. S., Lambin, E. F., Lenton, T. M., Scheffer, M., Folke, C., Schellnhuber, H. J., Nykvist, B., de Wit, C. A., Hughes, T., van der Leeuw, S., Rodhe, H., Sörlin, S., Snyder, P. K., Costanza, R., Svedin, U., . . . Foley, J. A. (2009). A safe operating space for humanity. Nature, 461, 472–475.
- Rositano, F., Bert, F. E., Piñeiro, G., & Ferraro, D. O. (2018). Identifying the factors that determine ecosystem services provision in Pampean agroecosystems (Argentina) using a data-mining approach. Environmental Development, 25, 3–11.
- Satorre, E., & Andrade, F. (2021). Cambios productivos y tecnológicos de la agricultura extensiva argentina en los últimos quince años. Ciencia Hoy, 173, 19–27.
- Senesi, S., Daziano, M., Chaddad, F., & Palau, H. (2017). Ownership versus management: The role of farming networks in Argentina. International Food and Agribusiness Management Review, 20, 221–238.
- Stevenson, J. R., Villoria, N., Byerlee, D., Kelley, T., & Maredia, M. (2013). Green Revolution research saved an estimated 18 to 27 million hectares from being brough into agricultural production. Proceedings of the National Academy of Sciences of the United States of America, 110(21), 8363–8368.
- Teixeira, E. I., de Ruiter, J., Ausseil, A.‑G., Daigneault, A., Johnstone, P., Holmes, A., Tait, A., & Ewert, F. (2018). Adapting crop rotations to climate change in regional impact modelling assessments. Science of the Total Environment, 616–617, 785–795.
- Volanti, M., Savarino, F. O., Passarini, F., Vassura, I., & Grosso, S. A. (2021). Environmental analysis of crop rotations through the application of the Grain Unit approach. Ecological Indicators, 121, 107199.
- Wooldridge, J. M. (2010). Econometric analysis of cross section and panel data (2nd ed.). MIT Press.
References
- Altieri, M. A. (2009). The ecological impacts of large-scale agrofuel monoculture production systems in the Americas. Bulletin of Science, Technology & Society, 29(3), 236–244.
- Alvarez, C., Alvarez, C. R., Costantini, A., & Basanta, M. (2014). Carbon and nitrogen sequestration in soils under different management in the semi-arid Pampa (Argentina). Soil and Tillage Research, 142, 25–31.
- Anderson, R., Bayer, P. E., & Edwards, D. (2020). Climate change and the need for agricultural adaptation. Current Opinion in Plant Biology, 56, 197–202.
- Andrade, F. H. (2020). Los desafíos de la agricultura global. Ediciones INTA.
- Andrade, J. F., Poggio, S. L., Ermácora, M., & Satorre, E. H. (2017). Land use intensification in the Rolling Pampa, Argentina: Diversifying crop sequences to increase yields and resource use. European Journal of Agronomy, 82, 1–10.
- Andriulo, A., Mary, B., & Guerif, J. (1999). Modelling soil carbon dynamics with various cropping sequences on the rolling pampas. Agronomie, 19(5), 365–377.
- Arora, P., Peterson, N. D., Bert, F., & Podesta, G. (2016). Managing the triple bottom line for sustainability: A case study of Argentine agribusinesses. Sustainability, 12(1), 60–75.
- Arrieta, E. M., Cuchietti, A., Cabrol, D., & González, A. D. (2018). Greenhouse gas emissions and energy efficiencies for soybeans and maize cultivated in different agronomic zones: A case study of Argentina. Science of the Total Environment, 625, 199–208.
- Asociación Argentina de Consorcios Regionales de Experimentación Agrícola. (2018). Base de datos integrada de rotaciones agrícolas en Argentina (Convenio AACREA-BID-FOMIN) (v. 2018–02) [Data set].
- Baldi, G., & Paruelo, J. M. (2008). Land-use and land cover dynamics in South American temperate grasslands. Ecology and Society, 13(2), 6.
- Berbery, E. H., Doyle, M., & Barros, V. (2006). Tendencias regionales en la precipitación. In V. Barros, R. Clarke, & P. Silva Días (Eds.), El cambio climático en la Cuenca del Plata (pp. 67–92). Consejo Nacional de Investigaciones Científicas y Técnicas.
- Bert, F. E., Podestá, G. P., Rovere, S. L., Menéndez, Á. N., North, M., Tatara, E., Laciana, C. E., Weber, E., & Toranzo, F. R. (2011). An agent based model to simulate structural and land use changes in agricultural systems of the Argentine Pampas. Ecological Modelling, 222(19), 3486–3499.
- Boulanger, J.–P., Leloup, J., Penalba, O., Rusticucci, M., Lafon, F., & Vargas, W. (2005). Observed precipitation in the Paraná-Plata hydrological basin: Long-term trends, extreme conditions and ENSO teleconnections. Climate Dynamics, 24, 392–413.
- Calviño, P., & Monzón, J. (2009). Farming systems of Argentina: Yield constraints and risk management. In V. Sadras & D. Calderini (Eds.), Crop physiology: Applications for genetic improvement and agronomy (pp. 55–70). Academic Press.
- Caride, C., Piñeiro, G., & Paruelo, J. M. (2012). How does agricultural management modify ecosystem services in the Argentine Pampas? The effects on soil C dynamics. Agriculture, Ecosystems & Environment, 154, 23–33.
- Carolan, M. S. (2005). Barriers to the adoption of sustainable agriculture on rented land: An examination of contesting social fields. Rural Sociology, 70(3), 387–413.
- Ciampitti, I., Rubio, G., Picone, L., & García, F. (2009). El fósforo en la agricultura: Mejores prácticas de manejo. In F. Garcia & I. Ciampitti (Eds.), Simposio Fertilidad 2009: Mejores prácticas de manejo para una mayor eficiencia en la nutrición de cultivos (pp. 35–51). IPNI Cono Sur.
- Farmaha, B. S., Eskridge, K. M., Cassman, K. G., Specht, J. E., Yang, H., & Grassini, P. (2016). Rotation impact on on-farm yield and input-use efficiency in high-yield irrigated maize–soybean system. Agronomy Journal, 108(6), 2313–2321.
- Garbarini, E., Skansi, M. M., Gonzalez, M. H., & Rolla, A. (2016). ENSO influence over precipitation in Argentina. In J. A. Daniels (Ed.), Advances in environmental research (pp. 223–246).
- García, F., Grasso, A., Gonzalez San Juan, M. F., Correndo, A., & Salvagiotti, F. (2018). Connecting crop nutrient use efficiency to future soil productivity. Better Crops, 102(4), 8–10.
- Giorgi, F. (2002). Variability and trends of sub-continental scale surface climate in the twentieth century. 717 Part I: observations. Climate Dynamics, 18(8), 675–691.
- Grimm, A. M., Barros, V. R., & Doyle, M. E. (2000). Climate variability in southern South America associated with El Niño and La Niña events. Journal of Climate, 13(1), 35–58.
- GROW. (n.d.). What is integrated weed management?
- Hazell, P. B. R. (2009). The Asian Green Revolution [Discussion paper]. International Food Policy Research Institute.
- Hora, R. (2012). La evolución del sector agroexportador argentino en el largo plazo, 1880–2010. Historia Agraria, 58, 145–181.
- Instituto Nacional de Estadística y Censos. (2021). Comercio exterior. Vol. 5, nº 4. Complejos exportadores: Año 2020. Informes Técnicos, 5(38), 11–14.
- Instituto Nacional de Estadística y Censos. (2021). Censo Nacional Agropecuario 2018: resultados definitivos.
- Intergovernmental Panel on Climate Change. (2006). IPCC guidelines for national greenhouse gas inventories, prepared by the National Greenhouse Gas Inventories Programme. In H. S. Eggleston, L. Buendia, K. Miwa, T. Ngara, & K. Tanabe (Eds.). 2006 IPCC guidelines for national greenhouse gas inventories. IGES.
- Intergovernmental Panel on Climate Change. (2019). Climate change and land: An IPCC special report on climate change, desertification, land degradation, sustainable land management, food security, and greenhouse gas fluxes in territorial ecosystems. In P. R. Shukla, E. Calvo Buendia, V. Masson-Delmotte, H.‑O. Pörtner, D. C. Roberts, P. Zhai, R. Slade, S. Connors, R. van Diemen, M. Ferrat, E. Haughey, S. Luz, S. Neogi, M. Pathak, J. Petzold, J. Portugal Pereira, P. Vyas, E. Huntley, K. Kissick, . . . J. Malley (Eds.). In press.
- Jobbagy, E. G., Aguiar, S., Piñeiro, G., & Garibaldi, L. A. (2021). Impronta ambiental de la agricultura de granos en Argentina: Revisando desafíos propios y ajenos. Ciencia Hoy, 29(173), 35–44.
- Kelley, K. W., Long, J., & Todd, T. C. (2003). Long-term crop rotations affect soybean yield, seed weight, and soil chemical properties. Field Crops Research, 83, 41–50.
- Kessler, J. J., Rood, T., Tekelenburg, T., & Bakkenes, M. (2007). Biodiversity and socioeconomic impacts of selected agro-commodity production systems. Journal of Environment & Development, 16(2), 131–160.
- Lamers, P., McCormick, K., & Hilbert, J. A. (2008). The emerging liquid biofuel market in Argentina: Implications for domestic demand and international trade. Energy Policy, 36(4), 1479–1490.
- Magrin, G. O., Travasso, M. I., & Rodríguez, G. R. (2005). Changes in climate and crop production during the 20th century in Argentina. Climatic Change, 72(1), 229–249. https://doi.org/10.1007/s10584-005-5374-9
- Manuel-Navarrete, D., Gallopín, G. C., Blanco, M., Díaz-Zorita, M., Ferraro, D. O., Herzer, H., Laterra, P., Murmis, M. R., Podestá, G. P., Rabinovich, J., Satorre, E. H., Torres, F., & Viglizzo, E. F. (2009). Multi-causal and integrated assessment of sustainability: The case of agriculturization in the Argentine Pampas. Environment, Development and Sustainability, 11(3), 621–638.
- Minetti, J. L., Vargas, W. M., Poblete, A. G., Acuña, L. R., & Casagrande, G. (2003). Non-linear trends and low frequency oscillations in annual precipitation over Argentina and Chile, 1931–1999. Atmósfera, 16, 119–135.
- Ministerio de Agricultura, Ganadería y Pesca. (2021). Estimaciones Agrícolas [Data set].
- Nosetto, M. D., Jobbágy, E. G., Brizuela, A. B., & Jackson, R. B. (2012). The hydrologic consequences of land cover change in central Argentina. Agriculture, Ecosystems & Environment, 154, 2–11.
- Novelli, L. E., Caviglia, O. P., & Melchiori, R. J. M. (2011). Impact of soybean cropping frequency on soil carbon storage in Mollisols and Vertisols. Geoderma, 167–168, 254–260.
- Penalba, O. C., Bettolli, M. L., & Vargas, W. M. (2007). The impact of climate variability on soybean yields in Argentina: Multivariate regression. Meteorological Applications, 14(1), 3–14.
- Podestá, G. P., Messina, C. D., Grondona, M. O., & Magrin, G. O. (1999). Associations between grain crop yields in central-eastern Argentina and El Niño–Southern Oscillation. Journal of Applied Meteorology, 38, 1488–1498.
- Qaim, M., & Traxler, G. (2005). Roundup Ready soybeans in Argentina: Farm level and aggregate welfare effects. Agricultural Economics, 32(1), 73–86.
- Rockström, J., Steffen, W., Noone, K., Persson, Å., Chapin, F. S., Lambin, E. F., Lenton, T. M., Scheffer, M., Folke, C., Schellnhuber, H. J., Nykvist, B., de Wit, C. A., Hughes, T., van der Leeuw, S., Rodhe, H., Sörlin, S., Snyder, P. K., Costanza, R., Svedin, U., . . . Foley, J. A. (2009). A safe operating space for humanity. Nature, 461, 472–475.
- Rusticucci, M., & Penalba, O. (2000). Interdecadal changes in the precipitation seasonal cycle over Southern South America and their relationship with surface temperature. Climate Research, 16, 1–15.
- Satorre, E., & Andrade, F. (2021). Cambios productivos y tecnológicos de la agricultura extensiva argentina en los últimos quince años. Ciencia Hoy, 173, 19–27.
- Scursoni, J. A., Duarte Vera, A. C., Oreja, F. H., Kruk, B. C., & De la Fuente, E. B. (2019). Weed management practices in Argentina crops. Weed Technology, 33(3), 459–463.
- Senesi, S., Daziano, M., Chaddad, F., & Palau, H. (2017). Ownership versus management: The role of farming networks in Argentina. International Food and Agribusiness Management Review, 20, 221–238.
- Stokstad, E. (2019). Jury verdicts cloud future of popular herbicide. Science, 364(6442), 717–718.
- Travasso, M. I., Magrin, G. O., & Rodríguez, G. R. (2003). Relations between sea-surface temperature and crop yields in Argentina. International Journal of Climatology, 23(13), 1655–1662.
- Trigo, E., & Cap, E. J. (2003). The impact of the introduction of transgenic crops in Argentinean agriculture. AgBioForum, 6(3), 87–94.
- Verón, S. R., de Abelleyra, D., & Lobell, D. B. (2015). Impacts of precipitation and temperature on crop yields in the Pampas. Climatic Change, 130(2), 235–245.
- Viglizzo, E. F., Frank, F. C., Carreño, L. V., Jobbágy, E. G., Pereyra, H., Clatt, J., Pincén, D., & Ricard, M. F. (2011). Ecological and environmental footprint of 50 years of agricultural expansion in Argentina. Global Change Biology, 17(2), 959–973.
- Wooldridge, J. M. (2010). Econometric analysis of cross section and panel data (2nd ed.). MIT Press.