Agricultural Transformations in the Southern Cone of Latin America: Agricultural Intensification and Decrease of the Aboveground Net Primary Production, Uruguay’s Case

The agri-exporting enclaves of the current corporate food regime intensively exploit natural assets in global strategies that govern the local processes. Their multidimensional impacts contribute to the environmental and civilizational crisis. Linked to the agrarian metabolism in its appropriation phase, land use has impacts on local systems. To understand this process, it is necessary to identify the systemic and spatial expression of environmental transformation. The objective of this work was to analyze the dynamic adjustment of aboveground net primary production (ANPP) to agricultural intensification between the years 2000 and 2017, using a land use intensity index and the soils’ productive potential. Agricultural expansion and consolidation are observed, as well as the significant intensification throughout the period in 65% of the country’s area—with differences between regions—associated with soil types. ANPP is higher in areas of low land use intensity and lower in high intensity areas, varying from high to low in soils with low to high productive potential, respectively, and growing throughout the period—depending on the area, with less growth in areas of greater intensity. The appropriation of edaphic wealth puts the systemic functionality at risk and challenges these transforming dynamics, with a strong impact on southern systems.


Introduction
The characteristics of the agri-food system in the production, distribution, and consumption phases work based on a complex network aimed at promoting the global interests of an increasingly small group of large transnational corporations that drive and govern the different links of the globalized agri-food chain [1]. Thus, the system intends to act as a "Food Empire" or "Corporate Empire" [2,3], aiming at unlimited capital expansion and accumulation instead of satisfying the population's food needs. According to Delgado Cabeza [4], in order to achieve such expansion, large corporations need to direct their efforts to erasing social and ecological limits to drive growth and accumulation, by developing strategies with four components: Financialization of the food system, controlling the rules of the food game, using new technologies, and using time-space.
Financialization refers to the increasing importance of financial capital within the operation of the food system, with the control, appropriation, speculation, and commercialization of nature's assets and components within the food sector. Financialization creates the conditions for appropriation by which is now 2.2 billion people [14], and the 821 million undernourished, attributable to the dominant global agri-food system [15]. The direction of these agri-exporting commodity farmers carries the displacement of the use of resources at local and regional levels. These resources are no longer the support of the local food supply, as they now nourish the needs of the global food business. This reality is faced in South America, and it is seen mainly regionally in the continent's warmer areas. The Pampa biome is the group of environmental grasslands in the southern part of the Río de la Plata basin. This biome is one of the most extensive warm grassland regions on the planet. It covers an area of 750,000 km 2 in the territories of Uruguay, Paraguay, Brazil, and Argentina (Soriano et al., 1992), where approximately 35 million people live. Globally, this biome comprises 2.5% of the world's agricultural land, with 0.45% of the world's population. It is a strategic area of the Río de la Plata basin, since it comprises territories of the lower basin occupying 40% of the basin's agricultural land and 35% of the population.
The territories of the basin of the Río de la Plata and, in particular, the Pampa biome, were geopolitically valued by the Spanish and Lusitanian crowns for the first time in the mid-16th century, as a living border territory and gateway into the continent which had to be closed tightly. After years of exhaustive geographical exploration, they decided to create port cities on the basin's main navigable river courses: Buenos Aires, Concepción de Buena Esperanza, and mainly, Asunción del Paraguay, all of them with a privileged geo-strategic position in South America [16]. Thus, given the biome's geographic importance, it was considered by colonizing nations as the final geographic hub in the Río de la Plata basin and South America's warm region during the last decades of the 16th century. Continuing this process, in recent decades, the biome's agricultural area has been used as a specialized production area, with a reduced number of primary production items destined for export. These are a minority sector of very high productivity, which is articulated with an extensive sector maintaining income levels due to the accumulation of land. The logic of the model driven by these sectors is based on the market's hegemony, the predominance of financial capital, and the promotion of transnational agribusiness companies [17]. So far in the 21st century, the agricultural area grew significantly in the Pampa biome, mainly due to the expansion of soybean crops. Argentina went from 27 to 36 million hectares (soybeans 43%), Brazil from 50 to 75 million hectares (soybeans 42%), Paraguay from 1.5 to 4 million hectares (soybeans 70%), and Uruguay from 0.800 to 1.6 million hectares (soybeans 63%) [17].
In this process, soybean crops expanded on natural ecosystems, exerting multidimensional impacts. In Brazil, cultivated and livestock fields replaced the savanna, pushing for a deforestation process which has affected a wide variety of unique ecosystems, from which 43 million hectares have disappeared in the last 30 years [18]. In Argentina, spaces such as the Yungas jungle, a very fragile ecological system which has been deforested by more than 70%, have been affected [19]. In Uruguay, soybeans expanded on natural grasslands, reducing their area and increasing pressure due to overgrazing. Several works report the degradation and loss of soils and/or fertility due to the "export" of natural nutrients in grains which leave from agri-exporting platforms. Regarding Argentina, Pengue estimates that, in 2002/2003, there were 4.5 million tons of nitrogen and phosphorus [20]. Soybean crops are also highly demanding of water: Approximately 1500 liters of water are used per kilogram of soybeans [21]. One of the most important aspects of placing soybean platforms is to have large water sources in the region [4]. In addition, the effects of the intensive use of agrochemicals must be considered, including glyphosate, a broad-spectrum and high-toxicity herbicide associated with the use of resistant transgenic seeds, the use of which has increased in a much greater proportion than soybeans crops [4].
Also, some social impacts are generally seen in the displacement of the population from the rural environment, the progressive impoverishment of the local population, the intensification of migratory movements towards large metropolitan areas or abroad, and the extension of poverty, malnutrition, and hunger [4]. The synergy between social and biophysical impacts compromises the system's functionality regarding the possibility of satisfying social demands in terms of the ability to export biomass. Paradoxically, this occurs in territories and countries specialized in agri-food production [4]. Such production is presented as an agricultural intensification process (AI) leading to the development of national economies [22]. AI should be considered as the process of (self) transformation of environmental systems through greater pressure on their structural and/or functional attributes in the biophysical dimension. This process generates a simpler configuration of production systems, more and more homogeneous and specialized, in which flow speeds increase, the biogeochemical cycles are modified, and the operation of the system is opened to the contribution of increasing numbers of inputs with greater dependence on external sources and a decrease in the overall capacity for internal regulation. As an operational concept, AI has been partially defined, strictly referring to the expansion of agriculture. However, this includes a set of interconnected actions: Increase in the number of crops and yields per area, inputs used, reduction of planned and unplanned components of biodiversity, dependence on the market economy, use of technology and capital, and a set of related negative environmental external factors. Therefore, AI integrates elements related to the management and use of the land and of related ecosystems, affecting the operation of environmental systems which still need to be known. This path has been difficult due to the lack of an operational methodology for space-time measurement and assessment. In addition, because AI affects a set of components of the environmental system, modeling and interpreting all the variables involved in the transformation is complex.
Within this framework, the intention is to contribute to understanding AI through a calculation methodology comprising an index of land use intensity based on the assessment of the main agricultural uses in the territory, and the linking of this index with the behavior of aboveground net primary production (ANPP), depending on the land's productive potential. Since the ANPP is a variable of the behavior of the environmental system comprising structural and/or functional attributes [23][24][25][26] and AI is a process of (self) transformation of these attributes, the analysis of this variable allows us to approach the operation of environmental systems [25].
We aim at contributing to a greater discussion about the impacts of the hegemonic model as a corporate food regime or food empire. Contributing to the debate where the questions of "what", "how", "who", "where", and "for whom" in terms of production are related to an environmental transformation strategy whose impacts become local (soils-ANPP), but are articulated in their social, economic, and environmental dimensions with modifications that, on a global scale, contribute to an unprecedented crisis. This paper proposes as its hypothesis that an emerging attribute of systems such as the ANPP is sensitive to the impact of transformation and responds with an inherent dynamic of productivity adjustment. The aim is to contribute to this analysis by answering five key questions on the spatio-temporal behavior of Uruguayan agricultural systems: What has been the spatial distribution of agricultural use intensity in Uruguayan soils during the first two decades of the 21st century? What has been the temporal evolution of land use intensity in Uruguay during this period? On what types of soils is the agricultural intensification process carried out in the country? How does this intensification process affect the systemic functionality visualized through the ANPP attribute? Is the effect different according to productive potential?

Materials and Methods
The methodological strategy comprised the spatio-temporal analysis of the ANPP using the Normalized Difference Vegetation Index (NDVI) based on the intensity of land use and its agricultural suitability. The enumeration area (EA) was used as a spatial scale, a unit of greater detail in the systematization of the official information collected in the agricultural censuses, which divide the Uruguayan territory into 637 areas (average surface area of 27,000 hectares).

Uruguay: Territorial Transformation
Uruguay's land territory is located in the warm area of South America, its latitude is between 30 • 05 and 34 • 58 south, and its longitude is between 53 • 11 and 58 • 26 west (Figure 1). The country's total area is over 17 million hectares (176,000 km 2 ); the climate is mild and humid, with annual accumulated rainfall ranging from 900 mm to 1500 mm between east and west, and average temperatures ranging between 17 • C and 20 • C. Most of the agricultural area is covered by natural grasslands and extensive livestock use (over 12 million hectares) [27]. Next, rain-fed crops-cereal and oilseed crops-located in the western region of the country (mainly soybeans, with approximately 1 million hectares), and irrigated crops in the eastern region of the country (mainly rice, with approximately 200,000 hectares) stand out. In addition, an important area is destined for afforestation of exotic species (approximately 1 million hectares), located mostly in four main regions [28]. and oilseed crops-located in the western region of the country (mainly soybeans, with approximately 1 million hectares), and irrigated crops in the eastern region of the country (mainly rice, with approximately 200,000 hectares) stand out. In addition, an important area is destined for afforestation of exotic species (approximately 1 million hectares), located mostly in four main regions [28]. In recent decades, the intensification process in Uruguay has sped up with the expansion of soybeans and afforestation affecting 2.1 million hectares (12.8% of the national territory) [30]. Regarding soybeans, in the last decade, the cultivated area grew steadily; it is currently over one million hectares, and is becoming the main agricultural sector [30]. The growth of agriculture was driven by the profitability linked to a cycle of high prices of commodities. This process has been promoted, defined, and led by new Uruguayan agricultural stakeholders. The articulation between progress, science, development, and capitalism spreads and self-accelerates [31] in the process of the ongoing modernization of the Uruguayan agricultural space [32]. Together with these transformations, new stakeholders have emerged, thus pressuring the different agrarian development models [33].
Agriculture reached its maximum area in the 1950s with 1.5 million hectares, decreasing its importance from then on [34]. In the year 2000, a minimum of 0.4 million hectares was reached. At that time, the discussion focused on the possibility of continuing agricultural production in Uruguay [35]. The expansion of soybean crops, which went from 12,000 hectares in 2000/2001 to over 1 million hectares in the 2012/2013 harvest, led to the significant increase of the area destined for agricultural activity. The cultivated area of other crops such as sorghum, corn, and wheat also increased during this period, delving into the country's agri-exporting nature. In the case of soybeans, Uruguay exports approximately 95% of the country's production as grain [36].
Expanding the area also carried significant technological changes with the use of direct sowing in virtually all crops, characterized by the use of transgenic seeds (soybean and corn) and the In recent decades, the intensification process in Uruguay has sped up with the expansion of soybeans and afforestation affecting 2.1 million hectares (12.8% of the national territory) [30]. Regarding soybeans, in the last decade, the cultivated area grew steadily; it is currently over one million hectares, and is becoming the main agricultural sector [30]. The growth of agriculture was driven by the profitability linked to a cycle of high prices of commodities. This process has been promoted, defined, and led by new Uruguayan agricultural stakeholders. The articulation between progress, science, development, and capitalism spreads and self-accelerates [31] in the process of the ongoing modernization of the Uruguayan agricultural space [32]. Together with these transformations, new stakeholders have emerged, thus pressuring the different agrarian development models [33].
Agriculture reached its maximum area in the 1950s with 1.5 million hectares, decreasing its importance from then on [34]. In the year 2000, a minimum of 0.4 million hectares was reached. At that time, the discussion focused on the possibility of continuing agricultural production in Uruguay [35]. The expansion of soybean crops, which went from 12,000 hectares in 2000/2001 to over 1 million hectares in the 2012/2013 harvest, led to the significant increase of the area destined for agricultural activity. The cultivated area of other crops such as sorghum, corn, and wheat also increased during this period, delving into the country's agri-exporting nature. In the case of soybeans, Uruguay exports approximately 95% of the country's production as grain [36].
Expanding the area also carried significant technological changes with the use of direct sowing in virtually all crops, characterized by the use of transgenic seeds (soybean and corn) and the increase in the use of inputs. Between 2003 and 2013, herbicides and agrochemicals in general increased fourfold [37]. The traditional rotation systems changed, reducing the importance of combined agricultural cycles with pastures and livestock, to move on to cycles of continuous agriculture. This also determines a process of agricultural "summering." The agricultural intensification process in Uruguay has created a link between agricultural production and soil deterioration processes, which decreases the potential for agricultural production. In particular, physical processes of soil erosion are identified as well as, mainly, the decrease of the organic carbon content in soils [28,[38][39][40][41].

NDVI-MODIS Time Series Data
The NDVI is one of the most used indexes for the assessment, tracking, and monitoring of numerous vegetation coverings and, mainly, for ANPP estimates due to its time and cost efficiency. This index combines the electromagnetic bands of red (R) and infrared (IR) (NDVI = (IR-R)/(IR+R)) [42], and has a positive relationship with the fraction of the photosynthetically active radiation absorbed by vegetation and, therefore, with primary productivity [43,44]. For this reason, it has been an indicator widely used worldwide for monitoring vegetation and, in general, changes in land cover.
The Moderate Resolution Imaging Spectroradiometer (MODIS) on board the Aqua and Terra satellites provides a 16 day composition of the NDVI at a moderate spatial resolution (250 meter pixel), generating 23 compositions per year. This product has been used in numerous investigations (e.g., [45][46][47][48]) due to its high quality, scale, and spatial and temporal cover, and the availability of updated and free data. In this study, 414 images of the MYD13Q1 product, provided free of charge by NASA EOSDIS, were used to cover the 2000-2017 period. Thus, a stack was built integrating pixels from the MODIS database from the entire Uruguayan territory during the 2000-2017 period.
NDVI values were averaged for each EA and 16 day period. Only pixels of the NDVI-MODIS images located fully inside the EA were used.

Land Use Suitability
A soil suitability index for agricultural (SSIA) rain-fed use was created. Numerical values were assigned to the conceptual categories defined by the MGAP [49], where 0 corresponds to unsuitable soils and 1 to soils with very high aptitude for agriculture. The MGAP's proposal involved determining crop suitability from the following variables: Erosion resistance, water storage capacity, pH, sodicity, rockiness and drainage; these were georeferenced in the soil charter with a scale of 1:20,000 [50]. This information was processed, systematized, and standardized at the EA scale. Categories defined by MGAP [49] were weighted using the following criteria: 0 for unsuitable, 0.25 for marginally suitable, 0.5 for moderately suitable, 0.75 for suitable, and 1 for very suitable. The SSIA value was the result of the weighted sum of each category of EA's relative area (Equation (1)), with theoretical values between 0 (100% of EA area with zero agricultural suitability) and 1 (100% of the EA area very suitable for agriculture). Finally, values were reclassified into four suitability categories according to SSIA values-low (<0.2), medium (0.2 ≤ x < 0.3), high (0.3 ≤ x < 0.5) and very high (≥0.5)-using the natural breaks classification method, as this method allows the grouping of previously arranged data and defines cuts where differences in values between groups are maximized regarding intragroup differences [51].
where SSIA-EA = soil suitability for agricultural use in EA i; Sni = area of land use suitability category n in EA i; SEAi = total area in EA i; Wn = land use weighting n.

Land Use Intensity
The land use intensity index (LUII) developed by Diaz et al. (2018) [52]-the most accurate land use information available-was used to evaluate land use intensity changes. For the first two temporary cuts, years 2000 and 2011, we worked with information collected in the general agricultural censuses [27,52]; in 2017, this information was updated using satellite images. Land use categories developed by Diaz et al. (2018) [52] were used, grouping land use categories surveyed in the 2000 and 2011 censuses.
The LUII was created to analyze land use intensity as a pressure indicator on natural assets based on official census data. It was carried out through the consultation of specialists on geography, agronomy, natural resources management, and ecology, through the weighting technique Analytical Hierarchical Process (AHP), developed by Saaty (1980) [53], which allowed the weighting of each land use and cover. After verifying the consistency of evaluations by decision makers using the consistency ratio and according to Saaty (1980) [53]-who states that if the ratio is less than 0.1, the evaluation is acceptable-20 evaluations were considered by Diaz et al. (2018) [54]. With this background as a basis, the weights suggest that, according to their decreasing intensity, land uses were as follows: Orchards (0.29), cereals and oilseeds (0.25), citrus, fruit and vineyards (0.17), afforestation (0.13), planted pastures and fodder (0.07), improved natural fields (0.04), natural fields (0.03), and wetlands and natural forests (0.02). The LUII was constructed on an EA scale and is the result of the weighted sum of each category's relative area (Equation (2)).
where IEA = intensity in EA i; S = area of land use suitability category n in EA i; SEA = total area in EA i; W = land use weighting n. Based on this proposal, the index was estimated for 2017. In this case, due to the lack of official information on agricultural censuses, the index was calculated with a supervised pixel-based classification of satellite images [55] from LANDSAT 8OLI for the year 2017 regarding the entire country, using the minimum distance method. For the supervision phase, 100 field data distributed throughout the territory were used to define the training areas. Accuracy was estimated with another 100 field data, resulting in 92% of certainty. This approach allowed the updating of the estimated areas of land uses of greater territorial extension (cereals and oilseeds, exotic afforestation, natural field, improved natural field, and wetlands and native forest).

Data Processing and Analysis
All spatial information was systematized and analyzed through the development and use of a geographic information system (GIS). Thus, it was possible to systematize MODIS satellite information, SSIA, and LUII at EA scale.
First, the existence of spatial autocorrelation was measured using the Moran's I (MI) [56]. After verifying the spatial autocorrelation, a "hot spot analysis" was performed in order to analyze the spatial distribution of land use intensity. This analysis uses the optimized Gi* statistic from Getis-Ord [57] to identify hot spots and cold spots (EA) in the context of neighboring entities (EA). Hot spots are entities with high values surrounded by entities with high values. When these cases are statistically significant, a z value is generated, which, at a higher value, indicates the grouping of entities with greater intensity. The opposite situation occurs when the value of the entity and its neighboring entities are low (cold points), generating negative z values.
Given the absence of normal distribution (Shapiro-Wilk test) and homogeneity of variance (Levene's test) in the data series, non-parametric statistical analyses were used. To analyze the relationship between SSIA and LUII, Kendall tau rank correlation coefficient (τ) analysis was applied.
The comparison of the NDVI values between the different categories of land use intensity and land use suitability (LUII and SSIA, respectively) was analyzed using the Kruskal-Wallis test (H) and later, as a post hoc test, the Mann Whitney test corrected (U), a non-parametric test for sample comparison.
For the purpose of analyzing the temporal behavior of the NDVI in the entire data sample and for each SSIA and LUII category, the Kendall tau rank correlation coefficient (τ) was used, analyzing the correlations between the dates (consecutively numbered) and the NDVI values. This test was used because the data does not meet the normality or independence assumptions.
Data analysis was first carried out for the entire period (2000-2017), and then for the 2000-2010 and 2011-2017 sub-periods, in order to consider the intercensal period and the current period. A significance level of 0.05 (p value) was established in all statistical tests.

Land Use Intensity
Changes in intensity show differences at the spatial level. In the southwest area, there is evidence of greater intensity of rural land use ( Figure 2). Agricultural production is distributed in four areas: South, southwest coast, northwest coast, and east, as observed in the spatial distribution of LUII for the year 2000 ( Figure 2). In 2011 and 2017, these areas remained unchanged, and the expansion of the agricultural frontier from the west coast towards the center of the country was evidenced, along with a very strong consolidation in the southwest coast (Figures 2 and 3), due to the expansion of rain-fed crops. In addition, the decreasing intensity in the east is highlighted due to a loss of relative importance of rice crops. The main spatial change took place in the first decade and was consolidated in the 2011-2017 period.
For 2000-2017, three large territorial processes of agricultural expansion in this period were clearly identified: (a) Consolidation of the west coast region as a key cereal and oilseed crops area, and the decrease in relative importance (at country level) of the southern horticultural area; (b) expansion towards the center of the country; (c) loss of relative importance (stagnation) of the southeast rice area, with an expansion of rain-fed crops in the eastern middle plains.

Land Use Intensity
Changes in intensity show differences at the spatial level. In the southwest area, there is evidence of greater intensity of rural land use ( Figure 2). Agricultural production is distributed in four areas: South, southwest coast, northwest coast, and east, as observed in the spatial distribution of LUII for the year 2000 ( Figure 2). In 2011 and 2017, these areas remained unchanged, and the expansion of the agricultural frontier from the west coast towards the center of the country was evidenced, along with a very strong consolidation in the southwest coast (Figures 2 and 3), due to the expansion of rain-fed crops. In addition, the decreasing intensity in the east is highlighted due to a loss of relative importance of rice crops. The main spatial change took place in the first decade and was consolidated in the 2011-2017 period.
For 2000-2017, three large territorial processes of agricultural expansion in this period were clearly identified: a) Consolidation of the west coast region as a key cereal and oilseed crops area, and the decrease in relative importance (at country level) of the southern horticultural area; b) expansion towards the center of the country; c) loss of relative importance (stagnation) of the southeast rice area, with an expansion of rain-fed crops in the eastern middle plains.
A presence of spatial autocorrelation for the LIUU values was detected for the years 2000 (MI = 0.85, z = 60.5, p < 0.0001), 2011 (MI = 0.41, z = 28.9, p < 0.0001), and 2017 (MI = 0.40, z = 28.2, p < 0.0001). During this period, the LUII showed an important growth. An increase in land use intensity values in 56% of the enumeration areas stand out, which represents 65% of the country's area (Figure 3).   When considering the country, the average LUII value for the year 2000 (0.2) showed a 12% growth compared to 2011 (0.22) and a 15% growth compared to 2017 (0.23) ( Table 1). These results indicate that, during said period, there were two major stages: a) 2000-2011, with significant agricultural expansion towards the center of the country, with an increase in land use intensity of 64% of the EA, and b) 2011-2017, with a consolidation of intensification, where, although intensification was only about 3.5%-three times less than in the previous period-intensification was of 71% of the total enumeration areas, and the rest remained the same. LUII distribution between such dates was positively and significantly correlated, although with significant variations between the two periods, 2000-2011 (τ = 0.64; p < 0.001) and 2011-2017 (τ = 0.95; p < 0.001), as a result of territorial expansion in the first period and consolidation in the second period.

Agricultural Potential Suitability
Of the EAs, 30% showed low agricultural potential suitability, 32% average suitability, 28% high suitability, and 10% very high suitability. Potential soil suitability was significantly higher in the west coast (SSIA was 50% higher than the country's average) and in the southern area (SSIA was 41% higher than the country's average). High potential areas were found in the center of the country and, to a lesser extent, spatially dispersed in the northeast area. Lower potential areas were located to the north and east of the country, and mean potential areas covered the largest area and were dispersed throughout the territory (Figure 4).
EAs of greater land use intensity were located in the areas with the highest SSIA. This link decreased between the year 2000 (τ = 0.56; p < 0.001) and 2011 (τ = 0.50; p < 0.001), and increased When considering the country, the average LUII value for the year 2000 (0.2) showed a 12% growth compared to 2011 (0.22) and a 15% growth compared to 2017 (0.23) ( Table 1). These results indicate that, during said period, there were two major stages: (a) 2000-2011, with significant agricultural expansion towards the center of the country, with an increase in land use intensity of 64% of the EA, and (b) 2011-2017, with a consolidation of intensification, where, although intensification was only about 3.5%-three times less than in the previous period-intensification was of 71% of the total enumeration areas, and the rest remained the same. LUII distribution between such dates was positively and significantly correlated, although with significant variations between the two periods, 2000-2011 (τ = 0.64; p < 0.001) and 2011-2017 (τ = 0.95; p < 0.001), as a result of territorial expansion in the first period and consolidation in the second period.

Agricultural Potential Suitability
Of the EAs, 30% showed low agricultural potential suitability, 32% average suitability, 28% high suitability, and 10% very high suitability. Potential soil suitability was significantly higher in the west coast (SSIA was 50% higher than the country's average) and in the southern area (SSIA was 41% higher than the country's average). High potential areas were found in the center of the country and, to a lesser extent, spatially dispersed in the northeast area. Lower potential areas were located to the north and east of the country, and mean potential areas covered the largest area and were dispersed throughout the territory (Figure 4).
EAs of greater land use intensity were located in the areas with the highest SSIA. This link decreased between the year 2000 (τ = 0.56; p < 0.001) and 2011 (τ = 0.50; p < 0.001), and increased slightly in 2017 (τ = 0.51; p < 0.001). High correlations indicate that the most intensive uses are mostly located on soils with greater agricultural suitability. In particular, at the beginning of the analyzed period, agriculture was developed in soils of greater suitability. Subsequently, the expansion towards the center of the country was carried out on soils with less potential, and in 2017, the trend was towards the use of soils with greater agricultural potential. slightly in 2017 (τ = 0.51; p < 0.001). High correlations indicate that the most intensive uses are mostly located on soils with greater agricultural suitability. In particular, at the beginning of the analyzed period, agriculture was developed in soils of greater suitability. Subsequently, the expansion towards the center of the country was carried out on soils with less potential, and in 2017, the trend was towards the use of soils with greater agricultural potential.

ANPP Behavior
The average NDVI was significantly higher in areas of lower land use intensity, followed by those of medium intensity and those of higher intensity with the lower values (H = 117, p < 0.0001). The coefficient of variation was higher in areas of greater intensity (9.7%) than of lower intensity (7.2%). The same variability pattern was found for both sub-periods.
NDVI median analysis for the entire period showed significant differences regarding soils' productive potentials (H = 313, p < 0.0001). Values ranged from highest to lowest NDVI in soils according to their suitability; in general, soils of lower agricultural suitability show higher NDVI values, and soils of greater suitability show the lowest NDVI values (U = 31877, p < 0.0001). The same pattern was observed for the two sub-periods considered, although in the last period, the differences between marginal and unsuitable soils were not significant ( Table 2).

ANPP Behavior
The average NDVI was significantly higher in areas of lower land use intensity, followed by those of medium intensity and those of higher intensity with the lower values (H = 117, p < 0.0001). The coefficient of variation was higher in areas of greater intensity (9.7%) than of lower intensity (7.2%). The same variability pattern was found for both sub-periods.
NDVI median analysis for the entire period showed significant differences regarding soils' productive potentials (H = 313, p < 0.0001). Values ranged from highest to lowest NDVI in soils according to their suitability; in general, soils of lower agricultural suitability show higher NDVI values, and soils of greater suitability show the lowest NDVI values (U = 31877, p < 0.0001). The same pattern was observed for the two sub-periods considered, although in the last period, the differences between marginal and unsuitable soils were not significant ( Table 2). In areas with very highly suitable soils, the NDVI did not show any significant trends throughout the period. A trend was identified in the remaining areas. In these cases, the NDVI's growth was greater, as the soil showed less agricultural suitability.
The NDVI within the group of agricultural areas showed a growing trend throughout the period (τ = 0.26, p < 0.001), which was explained by the growing trends of low intensity areas (τ = 0.27, p < 0.001) and of medium intensity areas (τ = 0.17, p < 0.001). On the other hand, in high intensity areas, the trend was not significant (p > 0.05).
In the first sub-period, no significant trends were identified for the group of areas. These differences were only identified for high intensity areas showing a decreasing trend in the NDVI (τ = −0.19, p < 0.001). In the second sub-period, a growing trend was seen (τ = 0.22, p < 0.001), which was explained by the growing trend of low intensity areas (τ = 0.22, p < 0.001) and of medium intensity areas (τ = 0.20, p < 0.001). On the other hand, in high intensity areas, the trend was not significant.
The analysis of the three areas according to their level of agricultural intensification showed, for the 2000-2017 period, an increasing trend of the NDVI in the three areas, although it was higher in areas which remained stable in their land use intensity, and showed lower values in areas where intensification was stronger (Table 3). This variation was different according to each sub-period. In the first sub-period, it was significant in areas of material changes showing a decreasing trend. In the second sub-period, it was increasing in all cases, and growth was greater in the areas where intensification was stable (Table 3).

Discussion
Differences, which show greater intensity of rural land use in the southwest, are consistent with the spatial pattern of historical agricultural uses of Uruguayan soils, which expands from Montevideo's surroundings towards the west [28]. In the first decade of this study, agricultural expansion from the west coast to the center of the country was confirmed. In this process, it is worth noting that the consolidation of cereal crops took place in soils of high agricultural suitability, but the expansion towards the center during the first period was developed on soils with less agricultural potential. The process of intensification and expansion of the agricultural area follows a similar pattern in all countries of Latin America's Southern Cone [17], and Uruguay does not escape this logic of use and appropriation of the agricultural space's natural assets. Expansion implies an advance on grasslands, hence modifying the system, since it removes the natural tapestry. The consequences of this transformation are linked to various processes, such as loss of the ability to retain carbon in soils [38], degradation [58], and loss of biodiversity associated with the Pampa biome [59]. In addition, there are changes in the "productive profile" of these soils linked to traditional livestock use, which is replaced by agriculture and new social stakeholders [60].
This transformation also results in a decrease in the relative importance of the southern area of fruit and vegetable production at the country level, which stood out until the year 2000. Its importance is linked to national food security and sovereignty, given that food production in this area is primarily destined for the domestic market [61]. Another aspect associated with this transformation is the loss of relative importance of the southeast rice area and the increase of rain-fed crops, which may be interpreted as part of the agricultural expansion of soybeans, but with a change in land use-mainly from rice crops, linked to "flooding" soils in middle plains and to rain-fed crops [62]. This shows a stagnation and relative impact on one of the most important crops in Uruguay's exporting scheme [63].
During this 17 year period, the intensification process expanded and consolidated, on the one hand, in existing agricultural areas, and on the other hand, in new agricultural areas, which shows a process of appropriation of natural assets and territorial transformations. Resulting impacts may be associated not only with the decrease of natural grasslands and their systemic use, but also with its linkage to social and economic transformations affecting traditional livestock systems, together with other impacts associated with the technology package, including increased use of biocides and fertilizers [36].
ANPP as a systemic attribute of environmental systems allowed the capture of part of this transformation. Therefore, it is evident that its growth trend is higher in areas with less intensification, and lower in those with greater intensification. The same behavior is observed regarding the suitability of land use: Higher NDVI levels are found in soils with lower agricultural suitability, and lower values are found in soils with higher suitability. These trends coincide with those found when analyzing the country's general tendency [54].
These systems of greater suitability are also those involving the greatest intensification, which seems to result in a lower ANPP. This behavior can be interpreted as the result of two factors acting simultaneously and reinforcing each other. A first factor, which would be explained by intensive land uses, reduces vegetable soil cover during a large part of the year (e.g., soybeans cover the soil with vegetation only in summer, while natural grasslands do so throughout the year). A second factor, explained by the export of nutrients from insensitive land use areas, is generated mainly in the harvest phase (system outputs) and, therefore, determines a "mobilization" of the productive potential of soils towards the productive results, as identified in local studies [64,65]. The first factor contributes greatly to the explanation of productivity at a specific time, while the second factor contributes to the explanation of the productivity trend over time.
By analyzing this behavior, it may be interpreted that the "mobilization" of the productive potential of soils towards productive results is increasingly lower. Considering this aspect together with the slowdown in cereal productivity [66] implies the understanding that productivity increases must be achieved at the cost of total inputs (fertilizers and biocides), which consequently increases environmental impacts [67].
From this interpretation, maintaining or increasing productivity may only be achieved on the basis of projecting greater expansion and intensification, consolidating a growing spiral in appropriation mechanisms determined by global logics, which are based on deepening the use of technology and capital [68].
The growing NDVI trend throughout the period, in all agricultural areas, is explained by a greater productivity in areas of low and medium use intensity. Although, at the country level, the NDVI may be offset between regions; local productive capacities are different, resulting in more serious consequences to more suitable and more intensively used soils, while expanding on less suitable soils with higher relative productivity, because such soils have undergone less intensification.
Increasing productivity on the basis of expansion, intensification, and greater dependence on inputs is the strategy of agricultural agribusiness, which articulates on a global scale both the input business and the production, marketing, and/or consumption business [69]. This results in local systems of wealth extraction and appropriation, in which the intensive exploitation of resources and local specificities takes place. Thus, the dynamics of local productive systems are governed from a global perspective [70], making evident their impacts which, from the point of view of the Southern Cone, may already be considered as regional [32]. Any appropriation made over the limits of the recovery of ecosystems/landscapes ends up showing in the decrease of the return flow in the short, medium, or long terms, or in the direct or indirect effects of the mechanisms used. The accumulation of decisions and/or actions in space and time for the benefit of maintaining extraction, concentration, and appropriation logics for some, and at the cost of others, may lead to the collapse of the material basis of rural metabolism, causing territorial conflicts in several dimensions.
Results obtained show the expansion of industrial agriculture by increasing the intensive use of inputs (energy, machinery, fertilizers, pesticides) and the greater control of production methods, which coincides with analyses carried out for the region [17].

Conclusions
LUII, SSIA, and NDVI are indicators which are complementary and adequate for capturing the spatio-temporal transformations of use intensity and understanding the spatial dynamics of the process based on soils' productive potentials.
In the period analyzed, agrarian intensification presents a spatial strategy that is structured from the west coast and moves towards the center of the country, which consolidates industrial agriculture and agribusiness in Uruguay. At the beginning, the growth of intensive agriculture is a slow process advancing on soils with greater productive potential. Subsequently, growth accelerates, and prioritizes the territorial consolidation of non-agricultural areas. Finally, when the process slows down, it retracts mainly on less suitable soils.
Lower PPNA values in areas with more suitable soils and greater intensification require further study of nutrient reserves in these soils and their management in the future. This knowledge would contribute to avoiding cycles of decreased productivity, decreased plant waste, changes in soil structure, and therefore, lost productivity.
The space-time pattern of primary productivity found in this work shows that productivity decreases in soils with greater agricultural suitability while remaining stable in other soils (or remains stable while increasing in other soils, or increasing less than in other soils). This pattern questions the strategy for the use and appropriation of nature's assets. The overexploitation of soils linked to agricultural businesses endangers primary productivity and generates great uncertainty regarding its medium and long term. Therefore, it is worth asking if we are facing a new impulse of colonialism, acting again as a colonized geographical pole in the Río de la Plata basin and the warm region of South America. In this sense, the region is facing the need to define strategies for its agrarian systems: Becoming functional territories of this logic of appropriation, or building options from the counter-hegemonic agrarian territories, the livelihoods of agriculture, and self-determined, resilient, and fair ways of living.