The study was conducted in the state of Bahia, located in the NEB, between the latitudes 8° 32’ and 18° 21’ S and longitudes 37° 20’ and 46° 37’ W (Fig. 1). The following data were used:
Daily meteorological data of minimum and maximum air temperature, average relative humidity, insolation, and average wind speed at 10 m collected in weather stations during a 10-year period (January 1, 1993 to December 31, 2002) for the municipalities studied were downloaded from the Meteorological Database for Teaching and Research (BDMEP) of the National Institute of Meteorology (INMET).
ERA5-Land reanalysis data on air temperature at 2 m height, dew point temperature at 2 m height, solar and thermal surface radiation, u and v components of wind at an hourly frequency were obtained from the ECMWF website for the period between January 1, 1993, and December 31, 2002, with a horizontal resolution of 0.1° x 0.1°. As ERA5 reanalysis data do not provide good rainfall estimates for the tropics (Gleixner et al. 2020; Lavers et al. 2022), daily rainfall data for the period between January 1, 1993, and December 31, 2002, were retrieved from the CPC/NOAA Precipitation Project, with a spatial resolution of 0.5º x 0.5º (Chen et al. 2008), because this dataset has been a useful reference to improve other data sources (Cui et al. 2017; Sadeghi et al. 2019). Due to the different spatial resolutions, the ERA5-Land reanalysis data were scaled to a 0.5° x 0.5° grid.
The Python 3.8 programming language and the Visual Studio Code editor were used to process all reanalysis data and calculate the following variables: relative air humidity, wind magnitude, radiation balance, and reference evapotranspiration.
Potential evapotranspiration, which can be estimated by several methods, is one of the key variables to the application of AZCR based on water balance. In the case of agricultural crops, the Food and Agriculture Organization (FAO) recommends the use of reference evapotranspiration. Thus, in this study, ET0 was estimated by the Penman-Monteith-FAO method proposed by Allen et al. (1998).
The program for crop water balance used in an EXCEL™ environment was created by Rolim et al. (1998) and the method for estimating the reference evapotranspiration was proposed by Allen et al. (1998), enabling the use of the variables degree-days-based crop coefficient (Kc), as described in Table 1, and variable available water capacity (AWC) of the soil, all organized in ten-day periods.
Table 1
Crop coefficient (Kc) for the early, medium, and late cycles with 100, 120 and 140 days, respectively. Values represent ten-day averages adjusted for the state of Bahia.
Cycle | Crop coefficient (ten-day values) |
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 |
Early | 0.40 | 0,42 | 0.63 | 0.90 | 1.12 | 1.15 | 1.15 | 1.14 | 0.95 | 0.67 | - | - | - | - |
Medium | 0.40 | 0.40 | 0.50 | 0.73 | 0.96 | 1.14 | 1.15 | 1.15 | 1.15 | 1.09 | 0.87 | 0.65 | - | - |
Late | 0.40 | 0.40 | 0.44 | 0.62 | 0.81 | 1.01 | 1.14 | 1.15 | 1.15 | 1.15 | 1.15 | 1.02 | 0.83 | 0.64 |
The mean AWC value in each quadrant of the study was calculated for the early, medium, and late cycles with 100, 120 and 140 days from the AWC values of each municipality of Bahia available at the National Water Agency, ANA (2021), as shown in Fig. 2:
In this study we analyzed the cycles classified as early (Group I), medium (Group II) and late (Group III) for the maize crop according to Ordinance number 170 of July, 2020, of the Ministry of Agriculture, Livestock and Supply of Brazil, which approves the AZCR for the 1st harvest of maize crop in the state of Bahia. To calculate the crop water balance, the maize cycle was divided into four phases, namely, germination/emergence (phase I), vegetative growth/development (phase II), flowering/grain filling (phase III), and physiological maturation/harvest (phase IV), respectively. The duration of each phase and the average size of the cycles are presented in Table 2.
Table 2
Average duration of cycles and respective phenological phases in days
Cycle | Phase I | Phase II | Phase III | Phase IV | Average cycle (days) |
Early | 15 | 35 | 30 | 20 | 100 |
Medium | 15 | 45 | 40 | 20 | 120 |
Late | 15 | 55 | 50 | 20 | 140 |
Source: Adapted. MAPA, 2020.
The water requirement satisfaction index (WRSI) developed by FAO (Frère and Popov, 1986) was used as a selection criterion for the start of the planting in each simulation of the crop water balance. This index consists of the ratio of actual to maximum evapotranspiration of the crop at a frequency of 80% established by Assad et al. (1998) and indicate a high (WRSI > 0.55), intermediate (WRSI between 0.45 and 0.55), or low (WRSI < 0.45) suitability of a site for planting maize, considering then that a site is suitable for planting maize when WRSI ≥ 0.45 (Santos et al., 2012).
Crop water balance simulations were carried out in ten-day periods. Thus a total of 36 periods per year were considered in the selection of dates suitable for sowing.
As described in the objectives, the study was divided into two stages: first, the validation of data from ERA5-Land reanalysis and the CPC/NOAA Precipitation Project by comparison with actual values from INMET weather stations; and second, the use of the validated data to elaborate the AZCR for maize cultivation in the state of Bahia as a grid with horizontal spatial resolution of 0.5° x 0.5° covering the entire state.
2.1 Validation
2.1.1 Weather stations
Five INMET weather stations in the state of Bahia were selected for the validation of ECMWF ERA5-Land reanalysis data and CPC/NOAA Precipitation Project because they were the ones that presented the lowest numbers of flaws among the stations of the state. The climatological aspects of the municipalities where the five stations are located regarding biomes and climate type according to Köppen-Geiger classification are shown in Table 3.
Table 3
Municipalities, average annual air temperature, annual precipitation, biome, climate classification according to Köppen-Geiger classification (As – tropical savanna climate with dry summer; Aw - tropical savanna climate with dry winter; BSh - hot semiarid climate; Cwb - subtropical highland climate), and mesoregion of the state of Bahia
City | Air temperature (ºC) | Precipitation (mm year-1) | Biome/ Climate Classification | Mesoregion |
Alagoinhas | 25.41 | 1147.96 | Atlantic Forest and Caatinga/As | Northeast Bahia |
Bom Jesus da Lapa | 26.55 | 814.55 | Cerrado and Caatinga/Aw | São Francisco Valley |
Cipó | 26.69 | 578.17 | Caatinga/BSh | Northeast Bahia |
Correntina | 23.17 | 935.37 | Cerrado/Aw | Far west of Bahia |
Vitória da Conquista | 21.43 | 735.44 | Atlantic Forest/Cwb | South-Central Bahia |
As scarce and incomplete weather station data pose a problem for time series analysis and applications in several regions of Brazil, the few locations where reliable information is available are used as a basic tool to validate reanalysis data.
To create a calendar indicating the periods suitable for growing maize under rainfed conditions, simulations of the crop water balance were performed using INMET weather station data and reanalysis data separately. In addition, a simulation was performed with precipitation data from the weather stations and air temperature and reference evapotranspiration derived from reanalysis data. Reference evapotranspiration was calculated from the variables air temperature and relative air humidity (maximum, minimum and average), average wind speed, solar radiation, and radiation balance. This is a “hybrid” approach and it was adopted because the state of Bahia has a good spatial distribution of pluviometric gauges, with daily and routine observations of precipitation only, but a scarce spatial distribution of weather stations.
The spreadsheet with the BHCult generates results in graphs of reference evapotranspiration (ET0), actual evapotranspiration (ETr), water storage, AWC, crop coefficient (Kc), and water balance history with the periods of surplus and deficit.
2.1.2 Parameters analyzed for validation
ERA5-Land reanalysis data and CPC/NOAA Precipitation Project a were validated and their reliability evaluated by comparing temperature (maximum and minimum), relative humidity, wind speed, and reference evapotranspiration estimates to weather station data. For this purpose, the statistical parameters root mean square error (RMSE), standard error of the estimate (SEE), mean absolute percentage error (MAPE), and mean absolute error (MAE) obtained in the comparisons of ERA5-Land reanalysis data and CPC/NOAA precipitation data with observational data were used in the analysis.
2.2 Agricultural zoning of climate risk for the state of Bahia
After the statistical validation of the data from ERA5-Land reanalysis and the CPC/NOAA Precipitation Project, a spatialization of agricultural zoning comprising the entire state of Bahia was proposed. For this, the ERA5-Land reanalysis data were organized in means grouped every 0.5º x 0.5º to equate the spatial resolution of the CPC/NOAA Precipitation Project throughout the territory of Bahia, within the extreme parallels of 8.25º S and 18.75º S and 37.25º W and 47.25º W, totaling an area with 232 grid points.