Sorghum (Sorghum bicolor L.) is the fifth most-cultivated cereal in the world. Currently, more than 35% is grown directly for human consumption  . According to  , it is the staple food for more than 500 million people living mainly in Africa and Asia. In these countries, it supplies 70% of the daily calorie intake, thus playing a fundamental role in food safety    .
Global production of sorghum for 2016/2017, according to the United States Department of Agriculture  , was 63.9 million tonnes. The United States leads world production with a little over 10 million tonnes, in second place is Mexico, with more than 7.0 million tonnes, and in third place, Nigeria, with around 6.0 million tonnes.
Brazilian production in the 2015/2016 season was 1,378,700 tonnes, according to data from the National Company for Supply  . The Central-West and Southeast of Brazil are noteworthy for their production of 780,200 tonnes and 439,200 tonnes respectively. The Northeast, equal to 72% of the semi-arid region of Brazil, had a production of 96,500 tonnes.
According to  , the semi-arid region of the Brazilian Northeast has an agricultural production which is intrinsically dependent on the rainfall in the region, but this regime is quite irregular, both on inter- and intra-annual scales. This uncertainty exposes the rural population, and agricultural production in particular, to a high level of vulnerability to the climatic conditions of the region. As an integral part of the semi-arid region, the State of Ceará is also characterised by high spatial and temporal variability in the rainfall regime  . Despite this, agricultural production in the state is based on rainfed agriculture  .
Data from the National Institute for the Semiarid  show that the semi- arid region of the Brazilian Northeast has around 23.8 million inhabitants, and is extremely vulnerable to variations in climate and the frequency of the droughts that mainly affect family farming, which consists for the most part of plantations under a rainfed system.
According to  , approximately 56% of the total cultivated area in India is under rainfed agriculture. The importance of rainfed agriculture can be evaluated by the fact that it contributes to 40% of the country’s food production.
Therefore, the question of the development of rainfed agriculture is critically important, due to growth and its implications for the subsistence security of a significant number of people living in the semi-arid region of the Northeast. Given the above, the aim here was to analyse the productivity and water-use efficiency of sorghum in a rainfed regime in the semi-arid region of Brazil.
2. Material and Methods
2.1. Location of the Study Area
The research was carried out in an area of the semi-arid region of Brazil (Figure 1) located in the Settlement of Lagoa Grande, which is supported by the National Institute for Colonisation and Agrarian Reform―INCRA. This has a total area of 2992.57 ha, where there are 92 farming families, and a total population of 392 inhabitants. The area of the settlement under rainfed agriculture is inserted
Figure 1. Location of the semi-arid region of Brazil and the area of rainfed agriculture.
in the Chapada do Apodi located geographically at 5˚20' latitude and 38˚0' longitude, at a distance of 30 km from the town of Tabuleiro do Norte, Ceará, Brazil.
The study was carried out with 12 producers, growing rainfed grain sorghum in an area equal to 138 ha of the universe of 250 ha used for agriculture in the above settlement, the study area being divided into 12 lots.
Data for the survey were obtained from spreadsheets and through interviews with the representatives of the family production units, where each producer represented one sample unit of the productivity of the rainfed sorghum.
2.2. Characterisation of Soil and Climate
The climate of the region is hot semi-arid, with average monthly temperatures greater than 18˚C (BSw'h'). The average annual rainfall is 790 mm, in a unimo- dal regime, with maximum rainfall concentrated from February to April. The predominant vegetation represents deciduous spiny forest, dense shrub-like caatinga, open brush-like caatinga, and mixed dicot-palm forest  .
Other characteristics of the climate in the region can be seen in Table 1, com- prising the mean values of a historical series from 2000 to 2016.
The soil in the area of rainfed agriculture of the Lagoa Grande Settlement lo- cated in the Chapada do Apodi is classified as a Red-Yellow Eutrophic Cambi- sol, deep, rich and well-drained. The characteristics physicochemical measured the 0.0 a 0.2 m depth is shown in the Table 2.
2.3. Description of the Production System
The first operation carried out by the farmers in the settlement was the lodging and desiccation of existing vegetation in the area where the sorghum was planted. The herbicide used was Glyphosate at a dose of 2.5 L∙ha−1 applied initially in the preparation phase of the area for planting. The prepared by mixing one litre of the herbicide in 100 litres of water, and applied with a 600 litre Jacto mechanical sprayer, driven by a 75 hp Massey Fergusson 296 tractor.
Table 1. Characteristics of the climate in the study area.
Source:  .
Table 2. Parameters physicochemical of soil of the area of rainfed agriculture.
Source:  .
The second operation was planting, using seeds of the sorghum hybrid “BR 304”. This was carried out under a system of direct seeding, using a model SP Light 3000 Baldanseeder, driven by a 75 hp Valtra 785 tractor. Planting was performed eight days after application of the herbicide, at a depth of 5 to 6 cm, and with 6 to 7 seeds per linear metre, at a spacing between rows of 0.45 m, giving a population of around 133,332 to 155,554 plants ha−1. Germination took place between 5 to 8 days after sowing.
The BR 304 is a simple hybrid, a result of the crossroad of lineage BR 001A with BR 012R, cycle of approximately 125 days. BR 304 is tolerant to the dry stem rot, caused by the pathogen Macrophomina phaseolina, under conditions of water deficiency. This characteristic becomes the BR 304 highly recommended for planting in succession to precose summer crops and late plantings. It is recommended for the Northeast, Midwest, Southeast and South regions, especially for times when rainfall is lower  .
Finally, harvesting, like the other operations, was also carried out mechanically using a CASE IH, model 7130 harvester. All the stages were performed under conditions of soil and climate that were favourable to each agricultural operation, since, according to  , tillage, sowing, cropping practices and harvesting should be carried out under adequate conditions of soil moisture, as they can cause a reduction in soil porosity and in the flow of liquid and gases, increase soil density, and consequently reduce crop productivity.
2.4. Analysed Variables
2.4.1. Productivity of the Sorghum
The mean productivity of the crop was analysed based on the ratio of sorghum production to the area of cultivation, as per Equation (1).
in which: Pme―Mean productivity (kg∙ha−1); p―sorghum production (kg); x― planted area (ha).
2.4.2. Water-Use Efficiency
Water-use efficiency (WUE) was obtained from the ratio of crop production to the volume of water applied, as per  , Equation (2).
in which: WUE―water-use efficiency (kg∙m−3); p―crop production (kg); v― volume of water applied (m3).
2.4.3. Financial Analysis
For the financial analysis, the methodology proposed by the Company for Technical Assistance and Rural Extension-EMATER-DF  was adopted, in which the financial concepts are defined by Equations (3)-(5). Given that total revenue refers to unit price multiplied by the quantity of produced goods.
in which: TR―Total revenue (BRL∙ha−1).
The cost of production is the sum of the fixed and variable costs, i.e. all the expenses generated during cultivation of the sorghum, such as tillage, planting, crop treatments, herbicides, seeds and labour, among others (Equation (4)).
in which: Cp―production costs (BRL); Cf―fixed costs (BRL); Cv―variable costs (BRL).
Net revenue refers to gross revenue minus product returns and taxes paid by the company (Equation (5)).
in which: NR―net revenue (BRL); GR―gross revenue (BRL); PE―general expenses of the production process (BRL); T―taxes (BRL).
2.5. Statistical Analysis
Descriptive statistics were used to analyse the following parameters: mean, standard deviation and coefficient of variation. The data were analysed using the SPSS v 16 software, with the Excel software being used for the graphs and regression analysis of the correlations between the different parameters under evaluation.
3. Results and Discussion
With the production data of the 12 lots of producers at the Grande Lagoa settlement, located in the semi-arid region of Brazil, values can be seen from the descriptive analysis, which vary from 438 kg∙ha−1 to 2143 kg∙ha−1, with a mean value of 919.42 kg∙ha−1 and a coefficient of variation of 49% (Table 3). Although the coefficient of variation had a high value, 75% of the values for productivity are within the range of the mean ± standard deviation (919 ± 469). Of all the lots under evaluation, 67% had lower productivity than the mean for the group, while 18% achieved a value for productivity statistically higher than the mean at a level of 5% probability.
The mean value for productivity found was lower than the values for produc- tivity registered in other Brazilian states, including that of the state itself. Ac- cording to a survey by  , the highest production of sorghum, 4500 kg∙ha−1, was found in the Federal District, and the lowest, 1006 kg∙ha−1, was obtained in the State of Bahia. In the State of Ceará, the average value was 1346 kg∙ha−1, for the 2015/2016 harvest.
Although the production lots received the same type of management, it can be seen from the coefficient of variation that productivity displayed high variability,
Table 3. Descriptive analysis of data for sorghum production.
this being attributed to the spatial variability of the soil attributes, and especially the spatial and temporal variation of the rainfall.
According to  and  , when analysing forage productivity in a Red Latosol they also found high variability (29%), which was attributed to the spatial variability of the soil attributes. On the other hand,  , studying the physical attributes of the soil, found that at both depths analysed, the coefficients of variation were classified as low, medium, high and very high, ranging from 10.7% to 32.9%.
As regards the data for sorghum production and productivity for each of the lots in the settlement (Table 4), it was found that lots 02, 03, 04 and 12 displayed higher values for productivity. This fact is attributed to the time of planting, since, especially in the case of rainfed agriculture, planning the agricultural calendar to take better advantage of the rainfall is a prime factor in the viability of a production system.
The low productivity and high variability are related to the distribution of rainfall during the period the sorghum was under cultivation, since the greatest limitation on agricultural production is not expressed by the total rainfall, but by the uncertainty of how these events are distributed over time, as well as the constant occurrence of periods with no rainfall (hot and dry) during the rainy season. According to  , and  , these hot, dry spells are characterised as an area of climate vulnerability for agricultural sustainability.
According to  and  further point out that although water is the limiting factor in agricultural production, the distribution of rainfall events over time and space is what most affects crop development and agricultural production. As found in a study by  , the low productivity seen in crops of forage sorghum on small farms in the north of the State of Minas Gerais is mainly due to the irregular distribution and low volume of rainfall during the development cycle of the sorghum.
Table 4. Sorghum production and productivity for each lot.
Regarding daily and monthly rainfall occurring in the study area during cultivation, these can be seen in Figure 2, where the highest values of 300 and 280 mm were recorded in March and April respectively, equivalent to 62% of the total rainfall, characterising the irregularity of rainfall distribution in the region.
Planting was concentrated from 1 April to 3 June, when rainfall during the growing season was 857 mm. It is noteworthy that from January to July, rainfall exceeded both the regional average of 790 mm∙year−1  and the historical average for the last 34 years, which was equal to 737 mm  .
Similarly,  point out that the total rainfall in the rainy season is sufficient for the development of rainfed agriculture in the intertropical zone, albeit affected by hot, dry spells; a phenomenon characterised by periods when the rainfall is interrupted during the rainy season.  add that the problem of irregular rainfall in the Northeast results not only from variations in total rainfall, but also from the duration and intensity of that rainfall.
Table 5 therefore, shows the number of dry spells and the length of each such event during cultivation of the sorghum; by means of the planting date, it was possible to identify which lot had the most influence on productivity as a function of these events. Thus, lots that were sown in April were less influenced by these periods, but also because the rainfall depth for March provided greater soil moisture at the time of planting, resulting in consequently higher productivity in lots 04, 12, 03 and 02, where planting was carried out on 20, 26, 27 and 01 April respectively.
On the other hand, lots 05, 06, 07 and 10 had lower values for productivity due to the time of planting, when they were influenced by the length of the dry spells, with six such periods being registered on average during cultivation, most
Table 5. Number of dry spells and duration of each event for the different Phenological stages of the crop in each lot.
Figure 2. Monthly rainfall during the sorghum cycle.
of them of six to seven consecutive days with no rainfall, and one period of 30 days, which contributed significantly to the low values for productivity found in these lots.
According to  , in a study in the State of Pernambuco in the semi-arid region of the Brazilian Northeast, found on average six to seven dry events per rainy season, with an average lack of rainfall of nine to eighteen days. However, longer dry spells of the order of 14 to 36 days were noted.
In addition to determining the number of dry spells that occurred, it was important to analyse the length of each event, as it was thereby possible to observe the high variability in productivity as a function of the spatial and temporal distribution of rainfall for each phenological stage of the crop, as well as the duration of the events for each physiological phase of the sorghum (Table 5).
Thus, it can be also be seen in Table 5 that during the first stage of crop growth (EC1), which runs from planting and germination to panicle initiation, approximately 30 days, there were around two to three dry spells, lasting from five to six days, in each lot; there was therefore no significant influence at this phenological stage for the low indices of productivity. However, the speed of germination, emergence and seedling establishment is important, since the plant is small, has slow initial growth, and at this stage, grain yield can be seriously reduced.
During the next stage (EC2), which includes panicle initiation to flowering, several growth processes, if affected, could compromise crop yield, among them, according to  , leaf area and root system development, dry matter accumula- tion and the establishment of a potential number of seeds, the latter being probably the most critical production component associated with increased yield.
It was found that during the EC2 stage, which goes from day 30 to day 70, there were on average two dry spells. However, it was the length of these events that had the most significant influence on the low productivity seen in lots 05, 06, 07 and 10, where the average length was 24 days. Whereas, for lots 02, 03, 04 and 12, where the duration of these events was around 5 to 6 days, productivity was not significantly compromised.
For  , any aggression to the plants at the EC2 stage, such as the undue application of agrochemicals or an unfavourable climatic event, such as a lack of soil moisture, will affect panicle emergence and compromise the final productivity of the crop.
During the third phase (EC3), which goes from flowering to physiological maturation and lasts a little more than 20 days, the factors considered the most important are those related to grain filling. At this stage, the plant continues to depend on a good level of water in the soil for adequate grain filling, with water deficiency during this period usually causing shrivelled grains and a fall in productivity  . It was found that lots with a low value for productivity were those that had an average of 30 days with no rainfall.
According to  point out that the hot, dry spells are a serious problem when they last longer than ten days without continuous rainfall, and especially when they affect the phenological stages of flowering and formation and of grain filling, stages EC2 and EC3 respectively; they are also harmful when they occur at the time of plant emergence, EC1. Therefore, comparing the data for productivity in Table 3 with the number of hot, dry spells in Table 5, a close relationship was seen between productivity and the occurrence of the dry spells, especially for the duration of these events.
Similarly,  point out that in rainfed production systems, the best use of the water resources is obtained by adjusting the time of sowing to the rainfall regime. According to  , climatic variability can significantly affect crop yield, even when the crop is sown at a suitable time.
Predicting the occurrence of hot, dry spells and their duration is therefore quite important, since it provides additional information for the planning of both rainfed and irrigated agriculture, helping to maximize water-use efficiency in areas under cultivation  .
3.1. Water-Use Efficiency
To calculate the volume of water per kilogram of sorghum produced in the rainfed system, the rainfall depth during the crop cycle and the total production of each lot were considered, and the amount of rainfall determined as a function of the area of each lot, as shown in Table 6.
It can be seen in Table 6 that water-use efficiency in the sorghum crop showed high variability as a natural result of the distribution of rainfall during the crop cycle in each lot, in addition to the variation in production values. The high value for lot 02 for the ratio of grain produced to water applied, of 4.83 m3∙kg−1, is explained by the input of a rainfall depth greater than that required by the crop, which is 300 - 325 mm. The best ratio of 1.13 m3∙kg−1 was seen in lot 04, which displayed maximum productivity.
Table 6. Water-use efficiency for each lot in the semi-arid region.
For  , improving water-use efficiency is a key factor in the continuous increase in crop productivity in arid and semi-arid regions.
Thus, it was found that the best ratio of grain production to water demand shows a strong correlation with water distribution throughout the cycle, and not only with the total volume; according to  , soil moisture conditions have a significant influence on water-use efficiency in arid and semi-arid regions. This fact demonstrates the great challenge of successful rainfed agriculture, where the irregularity of rainfall distribution was seen in this research as a limiting factor of greater weight than the actual volume of rainfall.
In their study,  found large variations in water-use efficiency in a crop of maize for the towns being studied, but found no direct relationship between grain yield and water-use efficiency in maize during the rainy season, or between rainfall and water productivity.
It is therefore necessary to adopt production systems that can adapt to this reality, as well as working together with meteorological agencies, taking advantage of rain forecasts, and adopting planting times in accordance with the information provided by such bodies.
3.2. Financial Analysis
For the financial analysis, the production costs per hectare of the sorghum crop were calculated, as shown in Table 7.
According to the data shown in Table 7, it was found that the cost of production for each hectare of planted sorghum was BRL 554.22, making a financial analysis of the system possible, as shown in Table 8.
According to the data shown in Table 8, it can be seen that only lots 01, 02, 03, 04 and 12 presented positive economic results, especially lots 02, 03, 04 and 12, with a good return on the capital invested, of the order of 11%, 20%, 132%
Table 7. Production costs for the sorghum crop.
Table 8. Financial analysis of the production system.
and 48% respectively. It is worth noting that these four lots were planted at the end of April, thereby suffering less influence from the hot, dry spells, allowing greater productivity of the sorghum, and consequently better economic return.
The low productivity of the sorghum under rainfed conditions shows that the lack of success of production systems in semi-arid regions is not due to total rainfall, but rather the spatial and temporal distribution of the rainfall during the rainy season, since the lowest values for productivity were seen with the occurrence of a greater number and the greater duration of hot, dry spells during cultivation.
The maximum yields of 1032 kg∙ha−1 (lot 02), 1102 kg∙ha−1 (lot 03), 2143 kg∙ha−1 (lot 04) and 1367 kg∙ha−1 (lot 12) were recorded when planting was carried out on 1, 27, 20 and 26 April respectively, when there were on average seven hot, dry spells of a shorter duration of around 8 days, while the lowest values for productivity occurred when planting was outside this period, with an average occurrence of six dry spells of a longer duration of around 17 days.
No direct relationship was found between grain yield and water-use efficiency in sorghum under rainfed conditions, or between rainfall and water productivity, since the best ratio of grain production to water demand shows a strong correlation with water distribution throughout the cycle and not only with the total volume.
Productivity under rainfed conditions is still low in some regions, requiring further research to determine the causes and set up strategies that would increase crop yield. This way, knowledge of the occurrence and duration of hot, dry spells, together with the time of planting of the sorghum, has the potential to approximate actual productivity to the productivity that can be achieved under a rainfed regime.
The authors acknowledge for financial support and the Coordination for the Improvement of Higher Education Personnel (CAPES) and Cearense Foundation of Scientific and Technological Support (FUNCAP). Thanks also go to the Secretary of Agrarian Development of the State of Ceará (SDA) and Association of Producers Irrigated Perimeter Curupati.
Conflicts of Interest
There are no conflicts of interest in present study.
 Mutisya, J., Sun, C., Rosenquist, S., Baguma, Y. and Jansson, C. (2009) Diurnal Oscillation of SBE Expression in Sorghum Endosperm. Journal of Plant Physiology, 166, 428-434.
 Dicko, M.H., Gruppen, H., Traoré, A.S., Castano, G., Menendez, R., Mas R., Voragen, A.J. and Berkel, W. (2006) Sorghum Grain as Human Food in Africa: Relevance of Content of Starch and Amylase Activities. African Journal of Biotechnology, 5, 384-395.
 USDA—United States Department of Agriculture. Foreign Agriculture Service, Production, Supply, and Distribution Database.
 CONAB—National Supply Company. Crop 2015/2016.
 Nelson, D.R., West, C.T. and Finan, T.J. (2009) Introduction to “In focus: Global Change and Adaptation in Local Places”. American Anthropologist, 111, 271-274.
 Guerreiro, M.J.S., Andrade, E.A., Abreu, I. and Lajinha, T. (2013) Long-Term Variation of Precipitation Indices in Ceará State, Northeast Brazil. International Journal of Climatology, 33, 2929-2939.
 Alves, J.M.B., Servain, J. and Campos, J.N.B. (2009) Relationship between Ocean Climatic Variability and Rainfed Agriculture in Northeast Brazil. Climate Research, 38, 225-236.
 INSA—National Institute of the Semiarid. Population of the Semi-Arid.
 Institute of Economic Research and Strategy of Ceará (2015) Ceará in Maps.
 National Institute of Meteorology Climatology (2017) Maps.
 Costa, S.C., Soares, A.A., Sediyama, G.C., Viana, T.V.A. and Moreira, F.V.O. (2011) Variation of Physical and Chemical Parameters of Fruit Fruits in ‘Pacovan’ Subjected to Different Irrigation Blades and Doses Potassium in Chapada do Apodi— Limoeiro do Norte-Ce. Irriga, Botucatu, 16, 82-92.
 Basso, F.C., Andreotti, M., Carvalho, M.P. and Lodo, B.N. (2011) Relationship between Forage Sorghum Yield and Physical Attributes and Organic Matter Content Latosol at the Cerrado. Pesquisa Agropecuária Tropical, 41, 135-144.
 Loomis, R.S. (1983) Crop Manipulations for Efficient Use of Water: An Overview. In: Taylor, H.M., Jordan, W.R. and Sinclair, T.R., Eds., Limitations to Efficient Water Use in Crop Production, American Society of Agronomy, Crop Society of America, and Soil Science Society of America, Madison, 345-374.
 Emater, D.F. (2009) Technical Assistance and Rural Extension of Distrito Federal.
 Pariz, C.M., Carvalho, M.P., Chioderoli, C.A., Nakayama, F.T., Andreotti, M. and Montanari, R. (2011) Spatial Variability of Forage Yield and Soil Physical Attributes of a Brachiaria Decumbens Pasture in the Brazilian Cerrado. Revista Brasileira de Zootecnia, 40, 2111-2120.
 Montanari, R., Zambianco, E.C., Corrêa, A.R., Pellin, D.M.P., Carvalho, M.P. and Dalchiavon, F.C. (2012) Physical Attributes of a Red Latosol Correlated Linearly and Spatially with the Consortium of Pigeon Pea with Millet. Ceres, 59, 125-135.
 Freddi, O.S., Carvalho, M.P., Veronesi Júnior, V. and Carvalho, G.J. (2006) Productivity of Maize Related to Mechanical Resistance to Soil Penetration under Conventional Tillage. Revista de Engenharia Agrícola, 26, 113-121.
 Silva, F.A.S. and Rao, T.V.R. (2002) Regimes Rainfall, Rainy Season and Probability of Occurrence of Dry Spells in the State of Ceará. Revista Brasileira de Engenharia Agrícola e Ambiental, 6, 453-459.
 Barron, J., Rockstrom, J., Gichuki, F. and Hatibu, N. (2003) Dry Spell Analysis and Maize Yields for Two Semiarid Locations in East Africa. Agricultural and Forest Meteorology, 117, 23-37.
 Goenster, S., Wiehle, M., Gebauer, J., Mohamed, A.A., Stern, R.D. and Buerkert, A. (2015) Daily Rainfall Data to Identify Trends in Rainfall Amount and Rainfall Induced Agricultural Events in the Nuba Mountains of Sudan. Journal of Arid Environments, 22, 16-26.
 Neto, O.A., Caxito, A.M., Chaves, F.F. and Pessoa, S.T. (2016) Production Silage Brs Ponta Negra Sorghum in Small Farms in the North of Minas Gerais. Proceedings of 31st Congresso Nacional de Milho e Sorgo “Milho e Sorgo: Inovações, Mercados e Segurança Alimentar”, Bento Gonçalves, 25-29 September 2016, 181-185.
 Foundation of Meteorology and Water Resources (2016) Historic Serie.
 Carvalho, D.F., Faria, R.A., Sousa, S.A.V. and Borges, H.Q. (2000) Spatialization of the Dry Spell Period for Different Levels of Loss of Production in Maize Rio Verde Grande, MG. Revista Brasileira de Engenharia Agrícola e Ambiental, 4, 172-176.
 Silva, L.L., Costa, R.F., Campos, J.H.B.C. and Dantas, R.T. (2009) Influence of Precipitation on Agricultural Productivity in the State of Paraíba. Revista Brasileira de Engenharia Agrícola e Ambiental, 13, 454-461.
 Soler, C.M.T., Sentelhas, P.C. and Hoogenboom, G. (2007) Application of the CSM CERES-Maize Mode for Planting Date Evaluation and Yield Forecasting for Maize Grown Offseason in a Subtropical Environment. European Journal of Agronomy, 27, 165-177.
 Andrade, C.L.T., Amaral, T.N., Silva, D.F., Garcia, A., Hoogenboom, G., Brito, R.A.L., Borges Júnior, J.C.F. and Gomide, R.L. (2009) Using the Ceres Maize Model as a Tool in the Definition of Maize Sowing Strategies: 3—Water Requirement and Productivity. Proceedings of XVI Congresso de Agrometeorologia, Belo Horizonte, 22-25 September 2009, 1528-1532.
 Carvalho, D.F., Oliveira, M.A.A., Sousa, S.A.V. and Carvalho, P.O.L. (1999) Estimates of Occurrences of Dry Spell in Seropédica, Vassouras e Piraí (RJ), and Influence on Yield of Bean (Phaseolus vulgaris L.). Ciência e Agrotecnologia, 23, 323-330.
 Xiao, G.J., Zhang, F.J., Huang, J.Y., Luo, C.K., Wang, J., Ma, F., Yao, Y.B., Wang, R.Y. and Qiu, Z.J. (2016) Response of Bean Cultures’ Water Use Efficiency against Climate Warming in Semiarid Regions of China. Agricultural Water Management, 173, 84-90.
 Tenhunen, J., Roupsard, O. and Rambal, S. (2002) Severe Drought Effects on Ecosystem CO2 and H2O Fluxes at Three Mediterranean Evergreen Sites: Revision of Current Hypotheses. Global Change Biology, 8, 999-1017.
 Correia, C.A., Andrade, C.L.T., Silva, P.P.G. and Magalhães, B.G. (2016) Simulation of Water Productivity for Rainfed Corn in Municipalities of Minas Gerais. Proceedings of 31st Congresso Nacional de Milho e Sorgo “Milho e Sorgo: Inovações, Mercados e Segurança Alimentar”, Bento Gonçalves, 25-29 September 2016, 1200-1203.