Open-access PRODUCTIVE RESPONSES FROM BROILER CHICKENS RAISED IN DIFFERENT COMMERCIAL PRODUCTION SYSTEM - PART II: IMPACT OF CLIMATE CHANGE

ABSTRACT

Broiler chickens are homoeothermic animals, i.e., animals capable of maintaining their body temperature within quite narrow limits; therefore, climate change poses a great challenge to poultry. With this in mind, this research aims to evaluate the performance of broilers submitted to different commercial production systems and exposed to different future scenarios, taking into account the climate change trends. To achieve this objective, we developed and validated a fuzzy model able to predict the performance of a broiler as a function of enthalpy along its life stages. This model was developed and validated in part I of this article based on experimental data collected for one year in three aviaries: conventional, negative pressure, and dark house systems. A Mann-Kendall nonparametric test and linear regression analysis were applied to the enthalpy values, which were calculated as a function of the ambient air temperature and relative humidity in order to study the climate change trends. Later, simulations were performed using the fuzzy model for 2025, 2050, 2075, and 2100 future scenarios. Specific improvements were observed when the heating trends coincided with the initial stages of breeding; however, in general, the productive responses of broilers in the different evaluated systems worsened with the climate change trends. Faced with the climate change trends, the responses improved in the order dark house aviary > negative pressure aviary> conventional aviary.

KEYWORDS broiler industry; artificial intelligence; climate change; fuzzy system

INTRODUCTION

The high temperature and air relative humidity values in tropical and semi-tropical climates, such as in Brazil, offer thermal discomfort conditions to chickens, decreasing their productive performance (Baracho et al., 2013; Boiago et al., 2013; Lara & Rostagno, 2013; Castro, 2014; Santos et al., 2014). The breeding environment, in the case of broilers, is one of the major causes of losses in commercial scale animal production. Several authors highlight that, for high genetic potential animals, it is necessary to provide adequate nutritious food and an aseptic environment thermally adjusted to the chicken's needs in order to allow them to express their full potential (Abreu et al., 2012, Almeida & Passini, 2013; Campos et al., 2013b; Nascimento et al., 2014; Tinôco et al., 2014).

Climate change poses a major challenge to the Brazilian poultry, as the chickens, which are homoeothermic animals, i.e., animals able to maintain their body temperature within quite narrow limits, suffer considerable productive losses when the thermal environment exceeds their comfort limits (Baracho et al., 2013; Boiago et al., 2013; Lara & Rostagno, 2013; Castro, 2014; Santos et al., 2014).

Broiler production is negatively affected by increases in gas emissions (Carbon dioxide - CO2, Methane -CH4, Nitrous Oxide - N2O, Chlorofluorocarbons - CFC-11, and CFC-12) responsible for the greenhouse effect (Gomes et al., 2011). Penereiro et al. (2012) reported that around 40% of the capitals of Brazilian states experienced an average temperature rise trend; four of them are located in the North region, two in the North-East region, two in the Midwest region, one in the South-East, and one in the South, the latter being Curitiba, Paraná, Brazil. Silva & Streck (2014) evaluated the climate change for the city of Santa Maria, Rio Grande do Sul, Brazil, through its monthly average heat index using the Mann-Kendall non-parametric test. Their results revealed that, in the months of March, at 9 am, and January, March, April, and November, at 3 pm, there is a significant positive trend, while the months of September, at 3 pm, and August, at 9 pm, have a negative trend. The most adequate method to study climate change in climatological series is the Mann-Kendall test (Goossens & Berger, 1986).

Therefore, the projection of broiler performance in places exposed to climate change is essential for improving the maintenance mechanisms responsible for keeping the thermal environment within the comfort intervals for genetic screening in order to explore the genetic potential of the lineage. For this purpose, it is necessary to develop environment control algorithms (mathematical models) that can be embedded in microcontrollers. Among the possible models to be developed, those based on artificial intelligence, specifically on the fuzzy set theory, have been adequate in animal comfort research (Gates et al., 2001; Castro et al., 2012; Ponciano et al., 2012; Campos et al., 2013a; Aborisade & Stephen, 2014; Ferraz et al., 2014; Xiang-Jie, 2014; Julio et al., 2015; Mirzaee-Ghalehv et al., 2015; Schiassi et al., 2015; Zare Mehrjerdi et al., 2015).

However, most fuzzy modeling studies have been individually carried out in one production system or in laboratories, where the environmental conditions are controlled. Thus, the fuzzy modeling conducted under these conditions may have limitations when applied to production systems different from the ones used in the studies. Among the existing broiler production systems, we can highlight the conventional system, the negative pressure system, and the dark house system.

The development of a mathematical fuzzy model based on different commercial production systems and on a given number of broilers raised in these systems shall forecast the broilers’ performance independently of the utilized system in different climate change scenarios.

In this context, the objective of this research was to apply a fuzzy model capable of predicting the performance of broiler chickens raised in different commercial production systems to simulate future scenarios considering possible climatic changes in the State of Santa Catarina.

MATERIAL AND METHODS

The fuzzy model

The fuzzy model applied to the simulations was developed based on data from three aviaries: the conventional, the tunnel with negative pressure, and the dark house breeding systems. The data were obtained from six batches for each system over the course of a year.

The fuzzy model considers the enthalpies as input data in five stages of breeding: Stage 1 (first week of life), stage 2 (second week of life), stage 3 (third week of life), stage 4 (fourth and fifth weeks of life), and stage 5 (from the sixth week of life on). Stages 1, 2, and 3 correspond to the initial breeding phase; stage 4 refers to the growing phase, and stage 5 is the termination phase. The output variables correspond to the productive responses of the broilers, food intake (FI), average weight gain (WG), average feed conversion (FC), and productive efficiency index (PEI).

Mamdani's inference method (Ponciano et al., 2012; Lin et al., 2013; Julio et al., 2015; Schiassi et al., 2015) and the defuzzification method using the center of gravity (Leite et al., 2010) were used. Two hundred and forty-three rules with weighting factors of 1.0 were elaborated (Yanagi Junior et al., 2012; Ponciano et al., 2012; Schiassi et al., 2013; Schiassi et al., 2014). Trapezoidal pertinence curves were utilized to represent the input (Schiassi et al., 2015) and the output variables.

The developed fuzzy model anticipates FI, WG, FC, and PEI with mean standard deviations and percentage errors of 4.16 g and 5.05%, 146.53 g and 8.04%, 0.06 g g−1 and 4.96%, and 24.51 and 12.29%, respectively. Additional details on the development and validation of the fuzzy model are available in Part I of this paper.

Mann-Kendall Test

The Mann-Kendall non-parametric test initially proposed by Sneyers (1975) and the linear regression were applied to the historical enthalpy series (Ávila et al., 2014, Tian et al., 2016). The enthalpy values (Equation 1) were calculated from the mean monthly values of the dry bulb temperature and air relative humidity for the period from 1987 to 2015. The measurements were made at Embrapa Swine and Poultry Agrometeorological Station in Concórdia, Santa Catarina, Brazil.

(1) H = 1 , 006 × t db + W × ( 2501 + 1 , 805 × t db )

Where,

H stands for enthalpy (kJ kg dry air−1);

tdb is the air dry bulb temperature (°C), and

W is the mixing ratio (kg water vapor kg dry air−1).

The mixing ratio was calculated by [eq. (2)] as a function of the current water vapor pressure (ea, kPa) and the local atmospheric pressure (Patm, kPa).

(2) W = 0 , 622 × ( e a P atm )

The Mann-Kendall (MK) test uses the concept of sequential order comparison of the time series value with the remaining values, counting the number of times the remaining terms are above the value analyzed. The S statistic is obtained by expressions 3 and 4.

(3) S = i = 2 n j = 1 i = 1 sinal ( X i X j )

Where,

(4) sinal ( Xi Xj ) = { 1 , for (Xi Xj) < 0 0 , for (Xi Xj) = 0 1 , for (Xi Xj) > 0

The S statistic tends to normality for high n, with mean and variance given by Expressions 5 and 6, respectively.

(5) E ( S ) = 0
(6) Var(S)= 1 18 [ n ( n 1 ) ( 2 n + 5 ) ]

Where,

n is the size of the time series.

Therefore, the ZMK statistic test is represented by Expression 7.

(7) Z MK = { S 1 Var(S) , if S > 0 0 , if S = 0 S + 1 Var(S) , if S < 0

By means of the Zmk value, we can verify the statistically significant trend in the time series, and this statistic is used to test the null hypothesis that no trend exists. From the positive or negative values of Zmk, it can be inferred if the data trend is increasing or decreasing. The significance level adopted was 5%, and significance was determined using the p-value test.

Simulations

The simulations were performed with the help of the MATLAB's Fuzzy Toolbox® software, 7.13.0.564 (R2011b) version, in which the entire modeling was designed (Part I of this article). The damping coefficient of the barns was calculated by means of the experimental enthalpy data measured inside the aviaries and in the external environment through the ratio between the external and internal enthalpies.

The years 2025, 2050, 2075, and 2100 were chosen as the future scenarios. The enthalpy increases for each future scenario were obtained by multiplying the angular coefficients of the linear regression applied to the enthalpy data series, for the months when climate change trends appear, by the number of years from 2016 on. Subsequently, the previously calculated enthalpy increase was added to the external enthalpy experimental data of each breeding stage, and the result was multiplied by the damping coefficient of each aviary. Then, the values of FI, WG, FC, and PEI were simulated using the fuzzy system developed and validated in part I of this article. The simulations were performed only for the batches raised in the months in which a trend of climate change was detected (February and June); the lots are identified in Table 1.

TABLE 1
Batches from each evaluated system and their respective evaluation periods.

RESULTS AND DISCUSSION

The analysis of the Mann-Kendall trend test for the monthly average enthalpy shows that, for the city of Concórdia, there was an increasing tendency for the months of February and June, as shown in Figure 1. There was no climate change trend for the other months. Several authors have pointed out that climate change results in more frequent, more intense, and longer heat waves (Renaudeau et al., 2012, Skuce et al., 2013, Sossidou et al., 2014).

FIGURE 1
Time trend of enthalpy in February (A) and June (B) for the city of Concórdia, SC, Brazil, in the 1987-2015 period.

In a study conducted by Vale et al. (2016), in which the effect of heat waves on the thermal and aerial environment for broilers was simulated, it was noticed that temperature rises linearly increased the ammonia concentration and bed temperature. These authors also noticed that temperature rises linearly reduced the oxygen concentration in the air.

According to the heating trend scenarios for February and June, we can see, in Table 1, that the batches affected by these months were batches 3, 4, and 6, and the productive performance results obtained by the fuzzy model for the batches affected by these months are shown in Table 2. It can be seen that, for the 2025 scenario, all batches from the evaluated systems maintained the same productive performances.

TABLE 2
Stages and batches affected by the climate change future scenarios and their respective food intake (FI, g), mean weight gain (WP, g), feed conversion (FC, g g−1), and productive efficiency index (PEI) values predicted by the fuzzy model as functions of enthalpy and broilers’ life stage for the future scenarios.

The Conventional aviary displayed the worst productive responses in the 2050, 2075, and 2100 scenarios for batches 3 and 6. Batch 4 exhibited the worst FC and PEI productive responses in the 2050 and 2075 scenarios. In turn, FI and WG in the 2050, 2075, and 2100 scenarios, and FC and PEI, in the 2100 scenario, showed improved productive responses for the respective batches.

The Negative Pressure aviary had the worst productive responses in the 2050, 2075, and 2100 scenarios for batches 3 and 6. In batch 4, the FC and PEI results remained unchanged up to the 2075 scenario, while the FI and WG responses remained unchanged for all scenarios.

For the dark house aviary, the productive responses of almost all batches remained unchanged or improved until the 2050 scenario, with the exception of the PEI for 2050, which fell from 393 to 392. We can see an improvement in batch 4, in all productive responses, for the 2050 and 2075 scenarios.

Studies carried out by Vale et al. (2010) and Vale et al. (2016) show that the increase in temperature associated with the increase in ammonia concentration and reduction in oxygen availability during periods of heat waves may justify the increase in mortality in heat wave events or even offer reasonable explanations for the lower performance of broilers, as observed in this study.

In general, for all the evaluated confinement systems, batch 4, raised from February to March (Table 1), achieved the best results when exposed to the heating trends of the future scenarios. This was the only batch able to improve the productive responses to the heating trends for the dark house and conventional aviaries or to keep them unchanged for a longer time for the negative pressure aviary. These results can be explained as follows: the heating tendency coincided with the batch's initial production phase, when the broilers require a higher environment temperature condition (Menegali et al., 2013), and because they experience a cold stress condition, the warming trend of future scenarios brings an improvement in the production environment for these stages.

According to Cassuce et al. (2013), chickens submitted to cold stress during their growth stage change their behavior, which affects their physiological and metabolic functions, hence decreasing their productivity.

In batches 3 and 6, in which the warming trend affected the final production stages, we can see a worsening of productive results from the middle of the twenty-first century, when chickens suffered heat stress and, as a consequence, a drop in final batch performance, impairing the activity's profitability (Carvalho, 2012).

Observing the annual means of the productive responses of all 6 batches (Table 3) in the three barns, we noticed that the broilers raised in the dark house system achieved the best results when faced with climate change trends, followed by those in the negative pressure system and then those in the conventional system, findings that were expected due to the different control levels of the systems.

TABLE 3
Feed intake (FI, g), mean weight gain (WG, g), feed conversion (FC, g g−1), and productive efficiency index (PEI) values as functions of enthalpy and broilers’ life stages predicted by the fuzzy model for future scenarios.

Even with occasional improvements due to the effects of the heating trends of these months on the breeding stage, it is observed that there was a worsening in productive responses in the different evaluated systems. This reinforces the need for artificial thermal conditioning systems and the use of materials with greater thermal inertia, as well as their correct management, which would allow for a reduction in the thermal environment variations inside the aviaries and, consequently, the problems related to thermal discomfort (Santos et al., 2014).

According to the Broiler Performance and Nutrition Supplement (Cobb-Vantress, 2015), the cumulative feed conversion for male broilers at 42 days of age is around 1.667 g g−1. This study showed that, for all future scenarios, the estimated mean feed conversion values are below the value stipulated by the lineage index (Cobb-Vantress, 2015).

CONCLUSIONS

In light of the expected climate change trends for the months of February and June, with monthly average enthalpy increase for the state of Santa Catarina, the response results improved in the order dark house aviary > negative pressure aviary > the conventional aviary.

There were occasional improvements when the heating trends coincided with the initial stages of breeding; however, in general, there was a worsening of the productive responses in the different evaluated systems. Thus, the use of materials with higher thermal resistance and inertia values, and artificial thermal conditioning systems, as well as their correct management, can lead to a reduction in thermal environment variations inside the aviaries and, consequently, problems related to thermal discomfort.

ACKNOWLEDGEMENTS

The authors thank FAPEMIG, CAPES, CNPq and EMBRAPA Swine and Poultry for their support of this research.

REFERENCES

  • Aborisade DO, Stephen O (2014) Poultry house temperature control using Fuzzy-PID controller. International Journal of Engineering Trends and Technology 11(6):310-314.
  • Abreu PG de, Abreu VMN, Coldebella A, Hassemer MJ, Tomazelli IL (2012) Medidas morfológicas em função do peso e da idade da ave, por meio de imagens. Revista Brasileira de Engenharia Agrícola e Ambiental 16(7):795-801.
  • Almeida EA, Passini R (2013) Conforto térmico em modelos reduzidos de casas de frangos de corte, sob diferentes tipos de materiais de cobertura. Engenharia Agrícola 33(1): 19-27.
  • Ávila LF, Mello CR, Yanagi SNM, Sacramento Neto OBS (2014) Tendências de temperaturas mínimas e máximas do ar no Estado de Minas Gerais. Pesquisa Agropecuária Brasileira 49(4):247-256.
  • Baracho MS, Cassiano JA, Nääs IA, Tonon GS, Garcia RG, Royer AFB, Santana MR (2013) Ambiente interno em galpões de frango de corte com cama nova e reutilizada. Agrarian 6(22):473-478.
  • Boiago MM, Barba H, Souza PA, Scatolini AM, Ferrari FB, Giampietro-Ganeco A (2013) Desempenho de frangos de corte, alimentados com dietas contendo diferentes fontes de selênio, zinco e manganês, criados sob condições de estresse térmico. Arquivo Brasileiro de Medicina Veterinária e Zootecnia 65(1):241-247.
  • Campos AT, Castro JO, Schiassi L, Yanagi Junior T, Pires MFÁ, Mattioli CC (2013a) Prediction of free-stall occupancy rate in dairycattle barns through fuzzy sets. Engenharia Agrícola 33(1):1079-1089.
  • Campos AT, Klosowski ES, Sousa FA, Ponciano PF, Navarini FC, Yanagi Junior T (2013b) Eficiência de sistema de aquecimento auxiliar para aviários, com base nos índices de conforto térmico. Bioscience Journal 29(3):703-711.
  • Carvalho LSS (2012) Nutrição de poedeiras em clima quente. Revista Científica Eletrônica de Medicina Veterinária 18:1-15.
  • Castro JO (2014) Avaliação e modelagem do desempenho de codornas japonesas em postura submetidas a diferentes ambientes térmicos. Tese Doutorado, Universidade Federal de Lavras.
  • Castro JO, Veloso AV, Yanagi Junior T, Fassani EJ, Schiassi L, Campos AT (2012) Estimate of the weight of Japanese quail eggs through fuzzy sets theory. Ciência e Agrotecnologia 36(1): 108-116.
  • Cassuce DC, Tinoco IDF, Baeta FC, Zolnier S, Cecon PR, Vieira MDA (2013) Atualização da temperatura de conforto térmico para frangos de corte de até 21 dias de idade. Engenharia Agrícola 33(1):28-36.
  • COBB (2015) Suplemento de nutrição e desempenho do frango de corte: Cobb-Vantress, 14 p. Available in: http://www.cobb-vantress.com/languages/guidefiles/fa217990-20c9-4ab1-a54e-3bd02d974594_pt.pdf Accessed: Jan 20, 2016.
    » http://www.cobb-vantress.com/languages/guidefiles/fa217990-20c9-4ab1-a54e-3bd02d974594_pt.pdf
  • Ferraz PFP, Yanagi Junior T, Julio YFH, Castro JO, Gates RS, Reis GM, Campos AT (2014) Predicting chick body mass with artificial intelligence-based models. Pesquisa Agropecuária Brasileira 49(7):559-568.
  • Gates RS, Chao K, Sigrimis N (2001) Identifying design parameters for fuzzy control of staged ventilation control systems. Computers and Electronics in Agriculture 31(1):61-74.
  • Gomes RCC, Yanagi Junior T, Lima RR, Yanagi SNM, Carvalho VF, Damasceno FA (2011) Predição do índice de temperatura do globo negro e umidade e do impacto das variações climáticas em galpões avícolas climatizados. Ciência Rural 41(9):1645-1651.
  • Goossens C, Berger A (1986) Annual and seasonal climatic variations over the northern hemisphere and Europe during the last century. Annales Geophysicae 4:385-400.
  • Julio YFH, Yanagi Junior T, Pires MFA, Lopes MA, Lima RR (2015) Fuzzy system to predict physiological responses of Holstein cows in southeastern Brazil. Revista Colombiana de Ciências Pecuárias 28(1):42-53.
  • Lara LJ, Rostagno MH (2013) Impacto do estresse térmico sobre a produção de aves. Animals 3(2):356-369.
  • Leite MS, Fileti AMF, Silva FV (2010) Desenvolvimento e aplicação experimental de controladores fuzzy e convencional em um bioprocesso. Revista Controle & Automação 21(2):147-158.
  • Lin CS, Yeh PT, Chen DC, Chiou YC, Lee CH (2013) The identification and filtering of fertilized eggs with a thermal imaging system. Computers and Electronics in Agriculture 91:94-105.
  • Menegali I, Tinoco IFF, Carvalho CCS, Souza CF, Martins JH (2013) Comportamento de variáveis climáticas em sistemas de ventilação mínima para produção de pintos de corte. Revista Brasileira de Engenharia Agrícola e Ambiental 17(1):106-113.
  • Mirzaee-Ghalehv E, Omid M, Keyhani A, Dalvand MJ (2015) Comparison of fuzzy and on/off controllers for winter season indoor climate management in a model poultry house. Computers and Electronics in Agriculture 110:187-195.
  • Nascimento GR, Nääs IA, Baracho MS, Pereira DF, Neves DP (2014) Termografia infravermelho na estimativa de conforto térmico de frangos de corte. Revista Brasileira de Engenharia Agrícola e Ambiental 18(6):658-663.
  • Penereiro JC, Ferreira DHL, Garcia VFV (2012) Caracterização estatística de tendências climáticas em séries anuais de temperatura média e precipitação para as capitais Estaduais do Brasil. Revista GEINTEC 3(4): 1-15.
  • Ponciano PF, Yanagi Junior T, Schiassi L, Campos AT, Nascimento JWB (2012) Sistema fuzzy para predição do desempenho produtivo de frangos de corte de 1 a 21 dias de idade. Engenharia Agrícola 32(3):446-458.
  • Renaudeau D, Collin A, Yahav S, Basilio V de, Gourdine JL, Collier RJ (2012) Adaptation to hot climate and strategies to alleviate heat stress in livestock production. Animal 6(5):707-728.
  • Santos GB, Sousa IF, Brito CO, Santos VS, Barbosa RJ, Soares C (2014) Estudo biológico das regiões litorâneas, agreste e semiárida do estado de Sergipe para a avicultura de corte e postura. Ciência Rural 44(1):123-128.
  • Schiassi L, Melo NSM, Tavares GF, Souza ÍP, Araújo HB, Della Giustina C (2013) Modelagem fuzzy em parâmetros de bem-estar humano. Nativa 1(1):8-12.
  • Schiassi L, Yanagi Junior T, Damasceno FA, Saraz JAO, Amaral AG (2014) Thermal-Acoustic Comfort Index for Workers of Poultry Houses Using Fuzzy Modeling. International Journal of Engineering Research and Applications 4(9):60-64.
  • Schiassi L, Yanagi Junior T, Reis GM, Abreu LHP, Campos AT, Castro JO (2015) Modelagem fuzzy aplicada na avaliação do desempenho de frangos de corte. Revista Brasileira de Engenharia Agrícola e Ambiental 19(2):140-146.
  • Silva SD, Streck NA (2014) Tendências das séries históricas do índice de calor no município de Santa Maria - RS. Ciência Rural 44(8):1360-1366.
  • Skuce PJ, Morgan ER, van Dijk J, Mitchell M (2013) Animal health aspects of adaptation to climate change: beating the heat and parasites in a warming Europe. Animal 7(2):333-345.
  • Sneyers R (1975) Sur l'analyse statistique des séries d'observations. Genève, Secrétariat de l’ Organisation Météorologique Mondiale.
  • Sossidou EN, Tsiplakou E, Zervas G (2014) Options for managing livestock production systems to adapt to climate change. Journal of Earth Science and Engineering 4(1):15-427.
  • Tian Q, Prange M, Merkel U (2016) Precipitation and temperature changes in the major Chinese river basins during 1957-2013 and links to sea surface temperature. Journal of Hydrology 536(1):208-221.
  • Tinôco IFF, Souza CF, Baêta FC, Coelho DJR, Mendes MASA (2014) Ambiência e Instalações na Avicultura de Postura Brasileira - Avanços e Perspectivas. Animal Busineess Brasil 4(14):6-9.
  • Vale MM, Moura DJ, Nääs IA, Curi TMRC, Lima KAO (2016) Effect of a simulated heat wave in thermal and aerial environment broiler-rearing environment. Engenharia Agrícola 36(2):271-280.
  • Vale MM, Moura DJ, Nääs IA, Pereira DF (2010) Characterization of heat waves affecting mortality rates of broilers between 29 days and market age. Brazilian Journal of Poultry Science 123(4):279-285.
  • Xiang-Jie N (2014) Research on the temperature control algorithm of the poultry farm. Applied Mechanics and Materials 602:1206-1209.
  • Yanagi Junior T, Schiassi L, Abreu LHP, Barbosa JA, Campos AT (2012) Procedimento fuzzy aplicado à avaliação da insalubridade em atividades agrícolas. Engenharia Agrícola 32(3):423-434.
  • Zare Mehrjerdi MR, Ziaabadi M, Ayatollahi Mehrgardi A, Dayani O (2015) Comparison of fuzzy and on/off controllers for winter season indoor climate management in a model poultry house. Journal of Livestock Science and Technologies 3(1):34-40.

Publication Dates

  • Publication in this collection
    Jan-Feb 2019

History

  • Received
    11 July 2018
  • Accepted
    29 Oct 2018
location_on
Associação Brasileira de Engenharia Agrícola Associação Brasileira de Engenharia Agrícola - SBEA, Departamento de Engenharia - FCAV/UNESP, Via de Ac. Prof. Paulo Donato Castellane, KM 05, CEP: 14884-900 , Phone: +55 (16) 3209-7619, WhatsApp: +55 (16) 98118-8978 - Jaboticabal - SP - Brazil
E-mail: revistasbea@sbea.org.br
rss_feed Acompanhe os números deste periódico no seu leitor de RSS
Acessibilidade / Reportar erro