Acessibilidade / Reportar erro

Crop-specific technology extrapolation domains for Brazil

ABSTRACT

Climate and soil databases can be used for delineating homogenous environmental zones over large areas of agricultural land to support agricultural research planning and spatial extrapolation of the results of field experiments. In other countries, such framework has proven useful to improve the efficiency and impact of investments in agricultural research. Furthermore, this approach can also contribute to identify key locations with the greatest potential for spatial extrapolation in terms of crop production area with similar climate and soils. In this paper, soil and weather databases covering the main agricultural land in Brazil were developed, and crop-specific technology extrapolation domains (TEDs) were identified for maize, soybean, and sugarcane crops. These databases were useful and viable option for the generation of homogeneous zones maps for different crops and regions. Given the spatial scale of the database applied, it was possible to represent a large area in few TEDs, without losing information regarding climate and soil variability, facilitating the spatial extrapolation of technologies throughout Brazil. Few TEDs are needed to represent a large part of the area of sugarcane (32 TEDs to cover 51% of the area), maize (21 TEDs for 51% of the area), and soybean (17 TEDs for 50% of the total area). Such information may help several areas of the agricultural sector, such as supporting new research that works on areas with large scales, or areas not explored yet.

Key words
agricultural zoning; homogenous zones; soybean; maize; sugarcane

INTRODUCTION

Average crop yields need to increase during the next decades to meet the expected increase in food demand while avoiding crop expansion (Mourtzinis et al. 2018Mourtzinis, S., Edreira, J. I. R., Grassini, P., Roth, A. C., Casteel, S. N., Ciampitti, I. A., Kandel, H. J., Kyveryga, P. M., Licht, M. A., Lindsey, L. E., Mueller, D. S., Nafziger, E. D., Naeve, S. L., Stanley, J., Staton, M. J and Conley, S. P. (2018). Siftin and winnowing: Analysis of farmer field data for soybean in the US North-Central region. Field Crops Research, 221, 130-141. https://doi.org/10.1016/j.fcr.2018.02.024
https://doi.org/10.1016/j.fcr.2018.02.02...
). Brazil has one of the world’s largest reserves of agricultural land, the fastest agricultural expansion during the 21st century, enormous natural resources, and a robust agricultural industry (FAO 2015[FAO] Foo and Agriculture Organization. (2015). Panorama de la inseguridad alimentaria en América Latina y el Caribe: La región alcanza las metas internacionales del hambre. Rome: FAO. [Accessed Dec. 4, 2020]. Available at: http://www.fao.org/3/a-i4636s.pdf
http://www.fao.org/3/a-i4636s.pdf...
).

Around 30% of Brazilian gross domestic product (GDP) depends on the agricultural sector (Martinelli et al. 2010Martinelli, L. A., Naylor, R., Vitousek, P. M and Moutinho, P. (2010). Agriculture in Brazil: impacts, cost and opportunities for a sustainable future. Current Opinion in Environmental Sustainability, 2, 431-438. https://doi.org/10.1016/j.cosust.2010.09.008
https://doi.org/10.1016/j.cosust.2010.09...
), while agricultural research expends a considerable amount of resources from both the private and public sectors. Research on agriculture as a driver of economic development has been widely recognized due to the cross-sectoral linkages between agriculture and nonagricultural economic sectors (Tiffin and Irz 2006Tiffin, R and Irz, X. (2006). Is agriculture the engine of growth? Agricultural Economics, 35, 79-89. https://doi.org/10.1111/j.1574-0862.2006.00141.x
https://doi.org/10.1111/j.1574-0862.2006...
). In general, studies on crop responses to management strategies such as type of cultivars, and use of fertilizers and herbicides are usually restricted to the local scale due to specific characteristics of soil and climate (van Wart et al. 2013Van Ittersum, M. K., Cassman, K. G., Grassini, P., Wolf, J., Tittonell, P and Hochman, Z. (2013). Yield gap analysis with local to global relevance - A review. Field Crop Research, 143, 4-17. https://doi.org/10.1016/j.fcr.2012.09.009
https://doi.org/10.1016/j.fcr.2012.09.00...
). Large amounts of public and private resources are expended yearly for collecting data from experimental field trials, either for testing scientific hypotheses or for evaluating new commercial products before releasing them onto the market (Pardey et al. 2006Pardey, P. G., Alaston, J. M., Chan-Kang, C., Magalhães, E. C and Vosti, S. A. (2006). Internationa and institutional R&D spillovers: attribution of benefits among sources for Brazil’s new crop varieties. American Journal of Agricultural Economics, 88, 104-123. https://doi.org/10.1111/j.1467-8276.2006.00841.x
https://doi.org/10.1111/j.1467-8276.2006...
).

Most of Brazilian agriculture has undergone high rates of expansion over the last four decades in the central area of the country. Even on traditional agricultural regions, however, there is a lack of knowledge on how soil and climate control crop growth and development at large scales, which would allow public and private decision makers to select where experiments should be established for more efficient resources use. For most of Brazil’s agricultural land, high quality climate and soil data at fine spatial resolution are not available for large important areas, and obtaining these data at many locations will be costly, laborious, and requires a long time for field studies and data acquisition. Such a lack data might limit the understanding of crop behavior and potential factors limiting producer yields (Edreira et al. 2017Edreira, J. I. R., Mourtzinis, S., Conley, S. P., Roth, A. C., Ciampitti, I. A., Licht, M. A., Kandel, H., Kyveryga, P. M., Lindsey, L. E., Mueller, D. S., Naeve, S. L., Nafziger, E., Specht, J. E., Stanley, J., Staton, M. J and Grassini, P. (2017). Assessing causes of yield gaps in agricultural areas with diversity in climat and soils. Agricultura and Forest Meteorology, 247, 170-180. https://doi.org/10.1016/j.agrformet.2017.07.010
https://doi.org/10.1016/j.agrformet.2017...
). Therefore, an upscaling method can be a type of “bottom-up” approach from location-specific observed data to larger spatial areas, seeking that the scale chosen should be small enough to minimize variation in climate and crop management practices and large enough to minimize data collection at regional scale (van Ittersum et al. 2013Van Wart, J., van Bussel, L. G. J., Wolf, J., Licker, R., Grassini, P., Nelson, A., Boogaard, H., Gerber, J., Mueller, N. D., Claessens, L., van Ittersumb, M. K and Cassman, K. G. (2013). Use of agro-climatic zones to upscale simulated crop yield potential. Field Crops Research, 143, 44-55. https://doi.org/10.1016/j.fcr.2012.11.023
https://doi.org/10.1016/j.fcr.2012.11.02...
).

Van Wart et al. (2013)Van Wart, J., van Bussel, L. G. J., Wolf, J., Licker, R., Grassini, P., Nelson, A., Boogaard, H., Gerber, J., Mueller, N. D., Claessens, L., van Ittersumb, M. K and Cassman, K. G. (2013). Use of agro-climatic zones to upscale simulated crop yield potential. Field Crops Research, 143, 44-55. https://doi.org/10.1016/j.fcr.2012.11.023
https://doi.org/10.1016/j.fcr.2012.11.02...
conducted a study on a global scale to seek an appropriate method using upscaling methodology based on the approach of homogeneous climate zones (CZ) and their climatic components, which can be used as a robust guide for analyzing large areas of agricultural interest. Edreira et al. (2018)Edreira, J. I. R.; Cassman, K. G., Hochman, Z., van Ittersum, M. K., van Bussel, L., Claessens. L and Grassini, P. (2018). Beyond the plot: technology extrapolation domains for scaling out agronomic science. Environmental Research Letters, 13, 054027. https://doi.org/10.1088/1748-9326/aac092
https://doi.org/10.1088/1748-9326/aac092...
developed a framework for the US agriculture, called technology extrapolation domains (TED), by using CZs and soil databases within a spatial analysis framework for identifying homogenous zones for crop growth and development in order to support agricultural research planning. Such a framework is useful to improve the efficiency and impact of investments in agricultural research, the more so because it is possible to customize it by adding other variables, including socioeconomic variables that are relevant for out-scaling and adoption of technology. Furthermore, such an approach can also contribute to identify key locations with the greatest potential for spatial extrapolation in terms of crop production area with similar climate and soils, as well as to select the number of sites needed to cover a particular region to be commercially addressed by a company. Finally, such a spatial framework to delineate regions in which crop production technologies perform similarly may overcome the limitations of extrapolation of the results of agronomic field experiments (Edreira et al. 2018Edreira, J. I. R.; Cassman, K. G., Hochman, Z., van Ittersum, M. K., van Bussel, L., Claessens. L and Grassini, P. (2018). Beyond the plot: technology extrapolation domains for scaling out agronomic science. Environmental Research Letters, 13, 054027. https://doi.org/10.1088/1748-9326/aac092
https://doi.org/10.1088/1748-9326/aac092...
).

Brazil is still lacking a spatial framework to guide technology transfer in cropping systems. In this paper, soil and weather databases covering the main agricultural land in Brazil were developed, following the protocols of Edreira et al. (2018)Edreira, J. I. R.; Cassman, K. G., Hochman, Z., van Ittersum, M. K., van Bussel, L., Claessens. L and Grassini, P. (2018). Beyond the plot: technology extrapolation domains for scaling out agronomic science. Environmental Research Letters, 13, 054027. https://doi.org/10.1088/1748-9326/aac092
https://doi.org/10.1088/1748-9326/aac092...
. The noncrop specific CZs worldwide provided by van Wart et al. (2013)Van Wart, J., van Bussel, L. G. J., Wolf, J., Licker, R., Grassini, P., Nelson, A., Boogaard, H., Gerber, J., Mueller, N. D., Claessens, L., van Ittersumb, M. K and Cassman, K. G. (2013). Use of agro-climatic zones to upscale simulated crop yield potential. Field Crops Research, 143, 44-55. https://doi.org/10.1016/j.fcr.2012.11.023
https://doi.org/10.1016/j.fcr.2012.11.02...
, which were based on large scale weather gridded data, were refined by using crop-specific relatively more detailed CZs based on locally measured weather data. The specific objectives were to identify the TEDs for maize, soybean, and sugarcane crops in Brazil, in order to offer a spatial guide for policy makers and private companies to strategically plan agricultural research and experimental results extrapolation.

MATERIAL AND METHODS

Study region and database

A database for soil and climate was developed covering 21 out of 26 states of Brazil and comprising the main agricultural area of the country. For climate, we used a grid of weather stations with well distributed daily data with an average resolution of 10 min (Fig. 1). The climate database used was constituted by 3,159 weather stations distributed across the main agricultural states of Brazil and organized as average monthly values of rainfall, and maximum and minimum air temperature. This database was based on daily weather data provided by several sources such as the National Institute of Meteorology (INMET), National Water Agency (ANA), University of São Paulo (USP), Brazilian Agricultural Research Enterprise (Embrapa), Agronomic Institute of Paraná (IAPAR), Department of Water and Electrical Energy (DAEE), and Agronomic Institute of Campinas (IAC). The weather data can be accessed in a spatial aggregated form in shapefile format at https://www.esalq.usp.br/gepema/teds/.

Figure 1
Spatial distribution of soil data points used for estimating the total available soil water and the weather stations used for creating the homogenous CZs.

Some weather stations have only rainfall data available, it being necessary to estimate maximum, minimum, and average temperatures. For those, the methodology proposed by Alvares et al. (2013a)Alvares, C. A., Stape, J. L., Sentelhas, P. C and Goncalves, J. L. M. (2013a). Modeling monthly mean air temperature for Brazil. Theoretica and Applied Climatology, 113, 407-427. https://doi.org/10.1007/s00704-012-0796-6
https://doi.org/10.1007/s00704-012-0796-...
was used, which correlates temperatures with geographical coordinates (latitude and longitude) and altitude in Brazil, according to Eq. 1:

T i = a o + a 1 φ + a 2 γ + a 3 h + a 4 φ γ + a 5 φ h + a 6 γ h + a 7 φ 2 + a 8 γ 2 (1)

where: Ti is minimum, maximum, and average temperature, monthly temperature (i = 1, 2, 3..., 12), φ is latitude in decimal degrees, is longitude in decimal degrees, h is altitude in meters, and to are the coefficients of the multiple regression equation.

Following Edreira et al. (2018)Edreira, J. I. R.; Cassman, K. G., Hochman, Z., van Ittersum, M. K., van Bussel, L., Claessens. L and Grassini, P. (2018). Beyond the plot: technology extrapolation domains for scaling out agronomic science. Environmental Research Letters, 13, 054027. https://doi.org/10.1088/1748-9326/aac092
https://doi.org/10.1088/1748-9326/aac092...
, the analysis of climatic features of each region was based on four variables: (i) annual total growing degree-days (GDD), which determines the length of the crop growing season; (ii) annual temperature seasonality, which differentiates between temperate and tropical climates; (iii) aridity index, which defines the water limitation in rainfed crop growing; and (iv) total available water (TAW) in the rootable soil depth, which determines the ability of the soil to supply water to support crop growth during rain-free periods, in which a root depth of 120 cm for sugarcane, 100 cm for maize, and 60 cm for soybean (Allen et al. 1998Allen, R. G., Pereira, L. S., Raes, D and Smith, M. (1998). Crop evapotranspiration - Guidelines for computing crop water requirements - FAO Irrigatio and drainage paper 56. Rome: FAO. [Accessed Dec. 4, 2020]. Available at: https://www.scscourt.org/complexcivil/105CV049053/volume3/172618e_5xAGWAx8.pdf
https://www.scscourt.org/complexcivil/10...
) was used. Therefore, each TED would correspond to a combination of these four variables.

To calculate the growing degree-days (GDD) for sugarcane, a base temperature at monthly time step data suggested by Marin et al. (2014)Marin, F. R and Jones, J. W. (2014). Process-Based simple model for simulating sugarcane growt and production. Scientia Agricola 71, 1-16. https://doi.org/10.1590/S0103-90162014000100001
https://doi.org/10.1590/S0103-9016201400...
was used, adopting a temperature base of 19.0 °C for sugarcane Marin et al. (2014)Marin, F. R., Ribeiro, R. V and Marchiori, P. E. R. (2014). How can crop modelin and plant physiology help to understand the plant responses to climate change? A case study with sugarcane. Theoretica and Experimental Plant Physiology, 26, 49-63. https://doi.org/10.1007/s40626-014-0006-2
https://doi.org/10.1007/s40626-014-0006-...
as follows:

For conditions in which the TM > Tb > Tm, Eq. 2 was used:

G D D   =   ( T M T b ) 2   +   ( T M 25 ) 2 2 ( T M t m )   x   f   x   N D (2)

and for Tb < Tm:

G D D   =   T M T   +   ( T M T b ) 2   +   ( T M 25 ) 2 2 ( T M t m )   x   f   x   N D (3)

in which,

f   =   N N - 24 2 (4)

where TM is average maximum temperature (°C), Tm is average minimum temperature (°C), Tb is base temperature (°C), N is the length of the days in hours, and ND is number of days.

For soybean and maize, GDD were calculated according to Licker et al. (2010)Licker, R., Johnston, M., Foley, J. A., Barford, C., Kucharik, C. J., Monfreda, C and Ramankutty, N. (2010). Mind the gap: how do climat and agricultural management explain the ‘yield gap’ of croplands around the world? Global Ecolog and Biogeography, 19, 769-782. https://doi.org/10.1111/j.1466-8238.2010.00563.x
https://doi.org/10.1111/j.1466-8238.2010...
, adopting a temperature base of 14.0 °C for soybean (Pascale et al. 1973Pascale, A. J., Villegas, J. A and Medina, L. F. (1973). Aptitud agroclimática del Noroeste Argentino para el cultivo de la soja. Revista Agronómica del Noroeste Argentino, 10, 173-203.) and 8.0 °C for maize (Warrington and Kanemasu 1983Warrington, I. J and Kanemasu, E. T. (1983). Corn growth response to temperatur and photoperiod I. Seedling Emergence, Tassel Initiatio and Anthesis. Agronomy Journal, 75, 749-754. https://doi.org/10.2134/agronj1983.00021962007500050008x
https://doi.org/10.2134/agronj1983.00021...
) as follows in Eq. 5:

G D D = Σ j = 1 365   max ( 0 . T i - t b ) (5)

where Ti is average temperature (°C) and Tb is base temperature (°C).

Temperature seasonality was calculated as the standard deviation of the average monthly temperatures × 100 (van Wart et al. 2013). The aridity index was calculated by the division between the mean annual precipitation and reference evapotranspiration (ETo), which was based on the approach developed by Camargo et al. (1999)Camargo, A. P., Marin, F. R., Sentelhas, P. C and Picini, A. G. (1999). Ajuste da equação de Thornthwaite para estimar a evapotranspiração potencial em climas áridos e superúmidos, com base na amplitude térmica. Revista Brasileira de Agrometeorologia, 7, 251-257. (Eq. 6), which in turn was proven to be well calibrated for the whole of Brazil:

E T o   =   0 . 01   Q o 2 . 45 . ( 1 . 08 . T M     0 . 36 . T m ) (6)

where ETo is the reference evapotranspiration (mm·day-1), Qo is extraterrestrial solar irradiation (MJ·m-2·d-1), TM is monthly average maximum temperature (°C), and Tm is average minimum temperature (°C).

As described, after creating the CZ following van Wart et al. (2013)Van Wart, J., van Bussel, L. G. J., Wolf, J., Licker, R., Grassini, P., Nelson, A., Boogaard, H., Gerber, J., Mueller, N. D., Claessens, L., van Ittersumb, M. K and Cassman, K. G. (2013). Use of agro-climatic zones to upscale simulated crop yield potential. Field Crops Research, 143, 44-55. https://doi.org/10.1016/j.fcr.2012.11.023
https://doi.org/10.1016/j.fcr.2012.11.02...
, the CZ were generated for the analyzed area in Brazil by combining with TAW to better represent the environmental conditions for the development and growth of the major crops grown in Brazil (van Bussel et al. 2015Van Bussel, L. G. J., Grassini, P., van Wart, J., Wolf, J., Claessens, L., Yang, H., Boogaard, H., Groot, H., Saito, K., Cassman, K. G and van Ittersum, M. K. (2015). From field to atlas: Upscaling of location-specific yield gap estimates. Field Crops Research, 177, 98-108. https://doi.org/10.1016/J.FCR.2015.03.005
https://doi.org/10.1016/J.FCR.2015.03.00...
). Soil and climate databases were interpolated using the ordinary kriging method with pixel size of 10 min. For the estimates of the TAW in the soil, pedotransfer functions provided by Tomassella et al. (2000)Tomassella, J., Hodnett, M. G and Rossato, L. (2000). Pedotransfer Functions for the estimation of soil water retention in Brazilian soils. Soil Science Society of America Journal, 64, 327-338. https://doi.org/10.2136/sssaj2000.641327x
https://doi.org/10.2136/sssaj2000.641327...
were used to estimate soil hydraulic parameters required, following the approach of van Genuchten (1980),Van Genuchten, M. T. (1980). A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Science Society of America Journal, 44, 892-898. https://doi.org/10.2136/sssaj1980.03615995004400050002x
https://doi.org/10.2136/sssaj1980.036159...
from soil texture, organic carbon content, and soil density, as provided by the soil databases of the Radambrasil Project (Brazil 1973–1986Brazil. Ministério das Minas e Energia. (1973-1986). Projeto Radambrasil. Levantamento de recursos naturais. Rio de Janeiro, Brazil: IBGE.) and Cooper et al. (2005)Cooper, M., Mendes, L. M. S., Silva, W. L. C and Sparovek, G. (2005). National soil profile database for Brazil available to international scientists. Soil Science Society of America Journal, 69, 649-652. https://doi.org/10.2136/sssaj2004.0140
https://doi.org/10.2136/sssaj2004.0140...
(Fig. 1). Compared to the soil database used by Edreira et al. (2018)Edreira, J. I. R.; Cassman, K. G., Hochman, Z., van Ittersum, M. K., van Bussel, L., Claessens. L and Grassini, P. (2018). Beyond the plot: technology extrapolation domains for scaling out agronomic science. Environmental Research Letters, 13, 054027. https://doi.org/10.1088/1748-9326/aac092
https://doi.org/10.1088/1748-9326/aac092...
, soil database can be considered too coarse, which represents a limitation for the study and it should be the focus of future studies to make the soil databases available for Brazil denser and more detailed.

Defining technological extrapolation domains for Brazil intervals and polygons

As mentioned, each TED can be assumed as a combination of CZ with the TAW in the rootable soil depth. Each climate variable had its own classes considering the spatial variability observed throughout the analyzed region of this study and for each crop. After Edreira et al. (2017)Edreira, J. I. R., Mourtzinis, S., Conley, S. P., Roth, A. C., Ciampitti, I. A., Licht, M. A., Kandel, H., Kyveryga, P. M., Lindsey, L. E., Mueller, D. S., Naeve, S. L., Nafziger, E., Specht, J. E., Stanley, J., Staton, M. J and Grassini, P. (2017). Assessing causes of yield gaps in agricultural areas with diversity in climat and soils. Agricultura and Forest Meteorology, 247, 170-180. https://doi.org/10.1016/j.agrformet.2017.07.010
https://doi.org/10.1016/j.agrformet.2017...
, some small areas within CZs were assumed as irrelevant for the study scheme, and two criteria were used to decide when to remove such polygons: those smaller than 350,000 ha and those having areas five times lower than the surrounding CZs.

Validating the technological extrapolation domains for Brazil

For validating the Brazilian TEDs of each crop, a survey of yield at municipality scale provided by the Brazilian Institute of Geography and Statistics (IBGE, http://www.sidra.ibge.gov.br) was used. Data from the last 10 years were used and the yield data for each TED component (climate and soil zones) were spatially aggregated, assuming that within each TED the yield variability is lower than the variability among TEDs. Moreover, the growing area data for each crop were used to define the number of TEDs necessary to represent the fraction of total area for maize, soybean, and sugarcane. Finally, comparisons among yield means associated with selected TED across Brazil were performed using the Tukey’s test (p < 0.05).

RESULTS AND DISCUSSION

Homogenous climate and soil zones

The GDD maps were closely related to temperature spatial variation in Brazil, with areas close to the Equator having higher GDD and no limitation for crop development through the year (Fig. S1). In high latitude areas in the south of Brazil, however, the maps showed some limited availability of GDD for growing crops in specific periods of the year, because of temperature seasonality. Given the wide range of latitudes where maize, soybean, and sugarcane have been cultivated in Brazil, the number of classes were considered adequate to represent the environmental spatial variability of the country. Such variability implies the need to select the right genotypes for each production zone, which is demonstrated by the diversity of genotypes offered to growers by private companies and public institutions in Brazil (Bustos et al. 2016Bustos, P., Caprettini, B and Ponticelli, J. (2016). Agricultural productivit and structural transformation: evidence from Brazil. American Economic Review, 106, 1320-1365. https://doi.org/10.1257/aer.20131061
https://doi.org/10.1257/aer.20131061...
). The differentiation in terms of GDD is also important for improving crop management, such as weeds, diseases, and pest control (López-Ovejero et al. 2005López-Ovejero, R. F., Garcia, A. G., Carvalho, S. J. P., Christoffoleti, P. J., Dourado Neto, D., Martins, F and Nicolai, M. (2005). Using thermal units for estimating critical period of weed competition in off-season maize crop. Journal of Environmental Scienc and Health, Part B, 40, 1-11.).

Related to the GDD maps (Fig. S1), but not with the same meaning, annual temperature seasonality maps show the variation of temperature throughout the year, in which regions near to the Equator line presented lower temperature variation than subtropical regions in southern Brazil (Fig S2a). The aridity index is an indicator of the level of dryness of the climate at a given location. Thus, areas with low rainfall and high temperature presented high degree of dryness (or low aridity index). On the other hand, areas with high rainfall and low temperature indicate a lower dryness degree (or a higher aridity index). The aridity index map (Fig. S2b) represented the spatial variability of Brazil well in terms of water availability, contrasting the humid regions in southern Brazil with the semiarid zones in northeast Brazil. This spatial distribution conforms with the natural biome distribution in Brazil and the climate classification of the country (Alvares et al. 2013bAlvares, C. A., Stape, J. L.; Sentelhas, P. C., Gonçalves, J. L. M and Sparovek, G. (2013b). Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift, 22, 711-728. https://doi.org/10.1127/0941-2948/2013/0507
https://doi.org/10.1127/0941-2948/2013/0...
).

By intersecting GDD maps for each crop, seasonal temperature, and aridity index, maps of crop-specific climatic zones (CZs) were created. In total, in Brazil, there are 85 CZs for maize, 84 CZs for soybean, and 91 CZs for sugarcane in the main agricultural states of Brazil (Fig 2). Compared with the number of CZs provided by van Wart et al. (2013)Van Wart, J., van Bussel, L. G. J., Wolf, J., Licker, R., Grassini, P., Nelson, A., Boogaard, H., Gerber, J., Mueller, N. D., Claessens, L., van Ittersumb, M. K and Cassman, K. G. (2013). Use of agro-climatic zones to upscale simulated crop yield potential. Field Crops Research, 143, 44-55. https://doi.org/10.1016/j.fcr.2012.11.023
https://doi.org/10.1016/j.fcr.2012.11.02...
for Brazil, a higher number of CZs was found due to the higher resolution of the weather database used for mapping the GDD (Fig. S1), aridity index (Fig. S2a), and temperature seasonality (Fig. S2b). van Wart et al. (2013)Van Wart, J., van Bussel, L. G. J., Wolf, J., Licker, R., Grassini, P., Nelson, A., Boogaard, H., Gerber, J., Mueller, N. D., Claessens, L., van Ittersumb, M. K and Cassman, K. G. (2013). Use of agro-climatic zones to upscale simulated crop yield potential. Field Crops Research, 143, 44-55. https://doi.org/10.1016/j.fcr.2012.11.023
https://doi.org/10.1016/j.fcr.2012.11.02...
used lower resolution gridded data and provided a worldwide framework for defining homogenous CZs, and that is the reason for the relatively lower resolution of their map.

Figure 2
Crop-specific climate zones for maize (a), soybean (b) and sugarcane (c) in the main agricultural States of Brazil.

Total available soil water (TAW) is basically related to the type of soil and physical characteristics, in which soils with higher clay content showed higher storage capacity, while less structured soils and sandy soils present less water storage capacity (Fig S3). The maps show soils in important productive regions of Brazil storing large amounts of water, and this could explain the good performance of agriculture in these regions. Intersecting CZs with TAW (Fig. 3) resulted in a total of 199 TEDs for maize, 201 TEDs for soybean, and 245 TEDs for sugarcane. This is the number of zones required to cover the main agricultural states of Brazil for the different crops and is essential for the efficient evaluation of a new technology to ensure that the number of field experiments can be reduced by covering larger scales based on climate and soil type. This number of TEDs is not considered high for a large country such as Brazil (Fig. 3). Rather, it is reasonable for rationalizing the number and location of experiments, and would be useful for providing spatial guidance to reduce the number of experiments to a minimum while representing Brazilian environmental diversity. For instance, considering that soybean has been cultivated on c.a. 35 million ha in the last season in Brazil and that such a large area could be represented by only 136 TEDs (Fig. 3), this information will help policy makers and companies to better plan research resources, retaining experimental sites that are well located and moving other experimental sites to locations that are not currently represented. Edreira et al. (2018)Edreira, J. I. R.; Cassman, K. G., Hochman, Z., van Ittersum, M. K., van Bussel, L., Claessens. L and Grassini, P. (2018). Beyond the plot: technology extrapolation domains for scaling out agronomic science. Environmental Research Letters, 13, 054027. https://doi.org/10.1088/1748-9326/aac092
https://doi.org/10.1088/1748-9326/aac092...
found 620 and 1,140 TEDs for moderate and high resolution, respectively. Comparatively, the low number of TEDs found for Brazil is mainly related to the low resolution of the soil data used for mapping the crop-specific TAW (Fig. S3).

Figure 3
Crop-specific technology extrapolation domains for maize (a), soybean (b) and sugarcane (c) in the main agricultural states of Brazil.

Planning the location of experimental sites based on the ted approach

Considering the representativeness of each TED in terms of cultivated area, it was calculated how many TEDs would be needed to represent a fraction of national area for each crop. For instance, only 21 TEDs would be necessary to cover 50.5% of the total maize area in Brazil (Fig. 4a). Soybean showed the most favorable scenario for saving resources by using the TED approach, as only 17 TEDs would be necessary to represent 50% of the total soybean area in Brazil (Fig. 4b). For sugarcane, 50 TEDs would be necessary to cover 69% of the sugarcane growing area in Brazil (Fig. 4c). However, a total of 100 TEDs would be needed to cover more than 94% of the area, while only 32 TEDs would cover more than 50% of the total sugarcane area in Brazil (Fig. 4c).

Figure 4
Crop area coverage as a function of number technology extrapolation domains (TEDs). TEDs were sorted fromlargest to smallest according to their 2018 harvested maize (a), soybean (b), and sugarcane (c) area (IBGE, http://www.sidra.ibge.gov.br).

Therefore, it is easy to wonder how much time and resources could be saved by using the TED approach in Brazilian agriculture. These resources are mainly related to agricultural researchers (reduction of work and costs for studies and tests in new areas), companies (implementation of existing technologies to other areas), and farms (benchmarking of other regions with similar climate and soil conditions), using this approach of location specific to larger spatial areas (Edreira et al. 2018Edreira, J. I. R.; Cassman, K. G., Hochman, Z., van Ittersum, M. K., van Bussel, L., Claessens. L and Grassini, P. (2018). Beyond the plot: technology extrapolation domains for scaling out agronomic science. Environmental Research Letters, 13, 054027. https://doi.org/10.1088/1748-9326/aac092
https://doi.org/10.1088/1748-9326/aac092...
). This means a massive reduction of experimental sites and costs, as environmental variability can be captured well by only a few TED and experiments should be planned to follow this approach.

Validating the TED approach for Brazil

As expected, soils with lower TAW presented low average yield when compared to those with higher TAW (Table 1). For example, counties with similar TAW but located in different CZs exhibited significant differences in average yield. Similarly, counties located in the same CZ, but with different TAW exhibited increasing yields in TEDs with greater TAW (Table 1). Moreover, the Tukey’s test showed that means differed from each other for most TAW classes within each selected CZ, this evidenced that the approach used for classifying CZ and TAW was able to adequately discriminate the environmental differences among different TEDs. Furthermore, the TED maps are also consistent with a well-validated simulation model that accounts for the effects of rainfall, temperature regime, and soil water availability on the crop physiology of sugarcane (Marin et al. 2008Marin, F. R., Lopes-Assad, M. L., Assad, E. D., Vian, C. E and Santos, M.C. (2008). Sugarcane crop efficiency in two growing seasons in São Paulo State, Brazil. Pesquisa Agropecuária Brasileira, 43, 1449-1455. https://doi.org/10.1590/S0100-204X2008001100002
https://doi.org/10.1590/S0100-204X200800...
; Marin and Carvalho 2012Marin, F. R and Carvalho, G. L. (2012). Spatio-temporal variability of sugarcane yield efficiency in the state of São Paulo, Brazil. Pesquisa Agropecuária Brasileira, 47, 149-156. https://doi.org/10.1590/S0100-204X2012000200001
https://doi.org/10.1590/S0100-204X201200...
; Marin and Jones 2014Marin, F. R and Jones, J. W. (2014). Process-Based simple model for simulating sugarcane growt and production. Scientia Agricola 71, 1-16. https://doi.org/10.1590/S0103-90162014000100001
https://doi.org/10.1590/S0103-9016201400...
, Vianna et al. 2020Vianna, M. S., Nassif, D. S. P., Carvalho, K. S and Marin, F. R. (2020). Modelling the trash blanket effect on sugarcane growt and water use. Computer and Electronics in Agriculture, 172, 105361. https://doi.org/10.1016/j.compag.2020.105361
https://doi.org/10.1016/j.compag.2020.10...
), maize (Souza et al. 2019Souza, T. T., Antolin, L. A. S., Bianchini, V. J. M., Pereira, R. A. A., Silva, E. H. F. M and Marin, F. R. (2019). Longer crop cycle lengths could offset the negative effects of climate change on Brazilian maize. Bragantia, 78, 622-631. https://doi.org/10.1590/1678-4499.20190085
https://doi.org/10.1590/1678-4499.201900...
), and soybean (Silva et al. 2019Silva, E. H. F. M., Gonçalves, A. O., Pereira, R. A., Fattori Júnior, I. M., Sobenko, L. R and Marin, F. R. (2019). Soybean irrigation requirement and canopy-atmosphere coupling in Southern Brazil. Agricultural Water Management, 218, 1-7. https://doi.org/10.1016/j.agwat.2019.03.003
https://doi.org/10.1016/j.agwat.2019.03....
). Based on the publicly available weather and soil data in Brazil, it was therefore considered that the crop-specific TED framework was able to capture the influence of soil type on yield within the same CZs in Brazil.

Table 1
Average maize, soybean, and sugarcane yield across selected climate zones (CZ) and total available water (TAW) within each CZ. Means followed by same lower-case letter in the row, and upper-case letter in the column do not differ by Tukey’s test (p < 0.05).

CONCLUSION

Weather and soil databases were found to be a useful and viable option for the generation of homogeneous zones maps for different crops and regions. It was possible to optimize the number of sites in Brazil to represent a large area in a few TEDs, without losing the variations in weather and soil conditions of a crop, facilitating technology transfer to different Brazilian regions. A low number of TEDs are needed to represent a large part of the acreage of sugarcane (32 TEDs to coverage 50.8% of area), maize (21 TEDs to coverage 50.5% of area), and soybean (17 TEDs to coverage 50.0% of acreage).

This study can contribute to several areas of the agricultural sector, such as supporting new research that works in areas with large scales, areas not explored yet, or areas that do not have climate and soil data. The shapefile delineating the TED spatial framework for Brazil is available for downloading at: https://www.esalq.usp.br/gepema/teds/.

Future studies should focus on the use of higher resolution soil data for improving the TAW map quality. Furthermore, the development of TED for other crops would also benefit the Brazilian agricultural industry.

ACKNOWLEDGMENTS

Not applicable.

REFERENCES

  • [FAO] Foo and Agriculture Organization. (2015). Panorama de la inseguridad alimentaria en América Latina y el Caribe: La región alcanza las metas internacionales del hambre. Rome: FAO. [Accessed Dec. 4, 2020]. Available at: http://www.fao.org/3/a-i4636s.pdf
    » http://www.fao.org/3/a-i4636s.pdf
  • Allen, R. G., Pereira, L. S., Raes, D and Smith, M. (1998). Crop evapotranspiration - Guidelines for computing crop water requirements - FAO Irrigatio and drainage paper 56. Rome: FAO. [Accessed Dec. 4, 2020]. Available at: https://www.scscourt.org/complexcivil/105CV049053/volume3/172618e_5xAGWAx8.pdf
    » https://www.scscourt.org/complexcivil/105CV049053/volume3/172618e_5xAGWAx8.pdf
  • Alvares, C. A., Stape, J. L., Sentelhas, P. C and Goncalves, J. L. M. (2013a). Modeling monthly mean air temperature for Brazil. Theoretica and Applied Climatology, 113, 407-427. https://doi.org/10.1007/s00704-012-0796-6
    » https://doi.org/10.1007/s00704-012-0796-6
  • Alvares, C. A., Stape, J. L.; Sentelhas, P. C., Gonçalves, J. L. M and Sparovek, G. (2013b). Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift, 22, 711-728. https://doi.org/10.1127/0941-2948/2013/0507
    » https://doi.org/10.1127/0941-2948/2013/0507
  • Bustos, P., Caprettini, B and Ponticelli, J. (2016). Agricultural productivit and structural transformation: evidence from Brazil. American Economic Review, 106, 1320-1365. https://doi.org/10.1257/aer.20131061
    » https://doi.org/10.1257/aer.20131061
  • Brazil. Ministério das Minas e Energia. (1973-1986). Projeto Radambrasil. Levantamento de recursos naturais. Rio de Janeiro, Brazil: IBGE.
  • Camargo, A. P., Marin, F. R., Sentelhas, P. C and Picini, A. G. (1999). Ajuste da equação de Thornthwaite para estimar a evapotranspiração potencial em climas áridos e superúmidos, com base na amplitude térmica. Revista Brasileira de Agrometeorologia, 7, 251-257.
  • Cooper, M., Mendes, L. M. S., Silva, W. L. C and Sparovek, G. (2005). National soil profile database for Brazil available to international scientists. Soil Science Society of America Journal, 69, 649-652. https://doi.org/10.2136/sssaj2004.0140
    » https://doi.org/10.2136/sssaj2004.0140
  • Edreira, J. I. R., Mourtzinis, S., Conley, S. P., Roth, A. C., Ciampitti, I. A., Licht, M. A., Kandel, H., Kyveryga, P. M., Lindsey, L. E., Mueller, D. S., Naeve, S. L., Nafziger, E., Specht, J. E., Stanley, J., Staton, M. J and Grassini, P. (2017). Assessing causes of yield gaps in agricultural areas with diversity in climat and soils. Agricultura and Forest Meteorology, 247, 170-180. https://doi.org/10.1016/j.agrformet.2017.07.010
    » https://doi.org/10.1016/j.agrformet.2017.07.010
  • Edreira, J. I. R.; Cassman, K. G., Hochman, Z., van Ittersum, M. K., van Bussel, L., Claessens. L and Grassini, P. (2018). Beyond the plot: technology extrapolation domains for scaling out agronomic science. Environmental Research Letters, 13, 054027. https://doi.org/10.1088/1748-9326/aac092
    » https://doi.org/10.1088/1748-9326/aac092
  • Licker, R., Johnston, M., Foley, J. A., Barford, C., Kucharik, C. J., Monfreda, C and Ramankutty, N. (2010). Mind the gap: how do climat and agricultural management explain the ‘yield gap’ of croplands around the world? Global Ecolog and Biogeography, 19, 769-782. https://doi.org/10.1111/j.1466-8238.2010.00563.x
    » https://doi.org/10.1111/j.1466-8238.2010.00563.x
  • López-Ovejero, R. F., Garcia, A. G., Carvalho, S. J. P., Christoffoleti, P. J., Dourado Neto, D., Martins, F and Nicolai, M. (2005). Using thermal units for estimating critical period of weed competition in off-season maize crop. Journal of Environmental Scienc and Health, Part B, 40, 1-11.
  • Marin, F. R., Lopes-Assad, M. L., Assad, E. D., Vian, C. E and Santos, M.C. (2008). Sugarcane crop efficiency in two growing seasons in São Paulo State, Brazil. Pesquisa Agropecuária Brasileira, 43, 1449-1455. https://doi.org/10.1590/S0100-204X2008001100002
    » https://doi.org/10.1590/S0100-204X2008001100002
  • Marin, F. R and Carvalho, G. L. (2012). Spatio-temporal variability of sugarcane yield efficiency in the state of São Paulo, Brazil. Pesquisa Agropecuária Brasileira, 47, 149-156. https://doi.org/10.1590/S0100-204X2012000200001
    » https://doi.org/10.1590/S0100-204X2012000200001
  • Marin, F. R and Jones, J. W. (2014). Process-Based simple model for simulating sugarcane growt and production. Scientia Agricola 71, 1-16. https://doi.org/10.1590/S0103-90162014000100001
    » https://doi.org/10.1590/S0103-90162014000100001
  • Marin, F. R., Ribeiro, R. V and Marchiori, P. E. R. (2014). How can crop modelin and plant physiology help to understand the plant responses to climate change? A case study with sugarcane. Theoretica and Experimental Plant Physiology, 26, 49-63. https://doi.org/10.1007/s40626-014-0006-2
    » https://doi.org/10.1007/s40626-014-0006-2
  • Martinelli, L. A., Naylor, R., Vitousek, P. M and Moutinho, P. (2010). Agriculture in Brazil: impacts, cost and opportunities for a sustainable future. Current Opinion in Environmental Sustainability, 2, 431-438. https://doi.org/10.1016/j.cosust.2010.09.008
    » https://doi.org/10.1016/j.cosust.2010.09.008
  • Mourtzinis, S., Edreira, J. I. R., Grassini, P., Roth, A. C., Casteel, S. N., Ciampitti, I. A., Kandel, H. J., Kyveryga, P. M., Licht, M. A., Lindsey, L. E., Mueller, D. S., Nafziger, E. D., Naeve, S. L., Stanley, J., Staton, M. J and Conley, S. P. (2018). Siftin and winnowing: Analysis of farmer field data for soybean in the US North-Central region. Field Crops Research, 221, 130-141. https://doi.org/10.1016/j.fcr.2018.02.024
    » https://doi.org/10.1016/j.fcr.2018.02.024
  • Pardey, P. G., Alaston, J. M., Chan-Kang, C., Magalhães, E. C and Vosti, S. A. (2006). Internationa and institutional R&D spillovers: attribution of benefits among sources for Brazil’s new crop varieties. American Journal of Agricultural Economics, 88, 104-123. https://doi.org/10.1111/j.1467-8276.2006.00841.x
    » https://doi.org/10.1111/j.1467-8276.2006.00841.x
  • Pascale, A. J., Villegas, J. A and Medina, L. F. (1973). Aptitud agroclimática del Noroeste Argentino para el cultivo de la soja. Revista Agronómica del Noroeste Argentino, 10, 173-203.
  • Silva, E. H. F. M., Gonçalves, A. O., Pereira, R. A., Fattori Júnior, I. M., Sobenko, L. R and Marin, F. R. (2019). Soybean irrigation requirement and canopy-atmosphere coupling in Southern Brazil. Agricultural Water Management, 218, 1-7. https://doi.org/10.1016/j.agwat.2019.03.003
    » https://doi.org/10.1016/j.agwat.2019.03.003
  • Souza, T. T., Antolin, L. A. S., Bianchini, V. J. M., Pereira, R. A. A., Silva, E. H. F. M and Marin, F. R. (2019). Longer crop cycle lengths could offset the negative effects of climate change on Brazilian maize. Bragantia, 78, 622-631. https://doi.org/10.1590/1678-4499.20190085
    » https://doi.org/10.1590/1678-4499.20190085
  • Tiffin, R and Irz, X. (2006). Is agriculture the engine of growth? Agricultural Economics, 35, 79-89. https://doi.org/10.1111/j.1574-0862.2006.00141.x
    » https://doi.org/10.1111/j.1574-0862.2006.00141.x
  • Tomassella, J., Hodnett, M. G and Rossato, L. (2000). Pedotransfer Functions for the estimation of soil water retention in Brazilian soils. Soil Science Society of America Journal, 64, 327-338. https://doi.org/10.2136/sssaj2000.641327x
    » https://doi.org/10.2136/sssaj2000.641327x
  • Van Bussel, L. G. J., Grassini, P., van Wart, J., Wolf, J., Claessens, L., Yang, H., Boogaard, H., Groot, H., Saito, K., Cassman, K. G and van Ittersum, M. K. (2015). From field to atlas: Upscaling of location-specific yield gap estimates. Field Crops Research, 177, 98-108. https://doi.org/10.1016/J.FCR.2015.03.005
    » https://doi.org/10.1016/J.FCR.2015.03.005
  • Van Genuchten, M. T. (1980). A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Science Society of America Journal, 44, 892-898. https://doi.org/10.2136/sssaj1980.03615995004400050002x
    » https://doi.org/10.2136/sssaj1980.03615995004400050002x
  • Van Ittersum, M. K., Cassman, K. G., Grassini, P., Wolf, J., Tittonell, P and Hochman, Z. (2013). Yield gap analysis with local to global relevance - A review. Field Crop Research, 143, 4-17. https://doi.org/10.1016/j.fcr.2012.09.009
    » https://doi.org/10.1016/j.fcr.2012.09.009
  • Van Wart, J., van Bussel, L. G. J., Wolf, J., Licker, R., Grassini, P., Nelson, A., Boogaard, H., Gerber, J., Mueller, N. D., Claessens, L., van Ittersumb, M. K and Cassman, K. G. (2013). Use of agro-climatic zones to upscale simulated crop yield potential. Field Crops Research, 143, 44-55. https://doi.org/10.1016/j.fcr.2012.11.023
    » https://doi.org/10.1016/j.fcr.2012.11.023
  • Vianna, M. S., Nassif, D. S. P., Carvalho, K. S and Marin, F. R. (2020). Modelling the trash blanket effect on sugarcane growt and water use. Computer and Electronics in Agriculture, 172, 105361. https://doi.org/10.1016/j.compag.2020.105361
    » https://doi.org/10.1016/j.compag.2020.105361
  • Warrington, I. J and Kanemasu, E. T. (1983). Corn growth response to temperatur and photoperiod I. Seedling Emergence, Tassel Initiatio and Anthesis. Agronomy Journal, 75, 749-754. https://doi.org/10.2134/agronj1983.00021962007500050008x
    » https://doi.org/10.2134/agronj1983.00021962007500050008x

SUPLEMENTARY FILES

Figure S1
Spatial variation of crop-specific growing degree-days for maize (a), soybean (b) and sugarcane (c) of the main agricultural states of Brazil.
Figure S2
Spatial variability of temperature seasonality (a) and aridity index (b) in the main agricultural states of Brazil.
Figure S3
Spatial variation of crop-specific total available water (mm) in soil for maize (a), soybean (b) and sugarcane (c) in the main agricultural states of Brazil.

Edited by

Section Editor:

Gabriel Constantino Blain

Publication Dates

  • Publication in this collection
    15 Jan 2021
  • Date of issue
    2021

History

  • Received
    03 July 2020
  • Accepted
    06 Oct 2020
Instituto Agronômico de Campinas Avenida Barão de Itapura, 1481, 13020-902, Tel.: +55 19 2137-0653, Fax: +55 19 2137-0666 - Campinas - SP - Brazil
E-mail: bragantia@iac.sp.gov.br