The E ﬀ ect of Best Crop Practices in the Pig and Poultry Production on Water Productivity in a Southern Brazilian Watershed

: This study analyzes the relation between Brazilian broiler and pig production and water productivity using recently developed reference guidelines on water footprinting for livestock production systems and supply chains. Di ﬀ erent rainfed crop arrangements, in di ﬀ erent scenarios and producer regions in Brazil, were assessed. Water productivity of broiler feed consumption ranged from 0.63 to 1.38 kg per m 3 water input to rainfed summer maize ( safra ) and from 1.20 to 2.21 kg per m 3 water input to winter maize ( safrinha ) while it ranged from 0.28 to 0.95 kg per m 3 water input to rainfed soy. For pig feed consumption, rainfed maize ranged from 0.68 to 1.49 kg per m 3 water input ( safra) and from 1.30 to 2.38 kg per m 3 water input ( safrinha ) while it ranged from 0.30 to 1.03 kg per m 3 water input to rainfed soy. A potential amount of water saving of 0.0336 km 3 year − 1 and 0.0202 km 3 year − 1 could be attained for producing broiler and pig feed, respectively, depending on the crop rotation and producer region. The results showed that the evapotranspiration of animal feed production represents more than 99% of the total water consumption for broiler and pig production in the study area. The implementation of best crop practices resulted in higher water productivity values of chicken and pork meat and also improved the rainfall water-saving in comparison to conventional agriculture. Hence, the water productivity of the animal production chain in tropical regions demands a close relation to agriculture in order to attain a better understanding and improvement of rainfall water productivity for animal feed production.


Introduction
Water is a crucial resource in agriculture [1], which stands out as the major user of water, accounting for about 70% of the world's freshwater withdrawal [2]. Freshwater demands are increasing worldwide and water scarcity is becoming a threat to the sustainable development of human society [3]. There is a need for widely recognized frameworks to assess the performance of e.g., livestock and livestock products in order to mitigate negative impacts on water resources [1]. Progress has been made recently in the field of assessing water demand in agriculture, and its related impacts on water scarcity, through the development of ISO 14046:2014 [4] "Environmental management-water footprint-principles, requirements and guidelines" and ISO/TR 14073:2016 [5] "Environmental management-water footprint-illustrative examples on how to apply ISO 14046".
This study is the first to apply the guidelines published by the FAO [1] for water productivity assessment on Brazilian livestock production. It aims to contribute to the knowledge related to water demand for animal production based on water productivity assessment and the importance of best practice implementation for raising water productivity of animal feed production in order to achieve a more sustainable production chain in Brazil.

General Approach
Water Productivity (WP) was calculated on farm-scale according to the FAO [1] for pig (WPindirect+direct,pig,Farm) and broiler (WPindirect+direct,broiler,Farm) production performed in the Lajeado Tacongava watershed in Southern Brazil for the year 2019. This study was focused on WP for purchased crops for animal feed. Summer and winter season maize crops, called safra and safrinha, respectively, and soybean produced in specific crop rotations were analyzed, considering the current yield. Potential strategies for raising water productivity, as higher yield and lower soil evaporation, were also analyzed. Evapotranspiration as indirect water consumption was considered to calculate feed water productivity (WPindirect,broiler,Feed; WPindirect,pig,Feed). The direct water consumptions were animal drinking (Qdirect,pig,Animal; Qdirect,poultry,Animal) and water for services (Qdirect,pig,Housing; Qdirect,poultry,Housing).

Study Area and Animal Production System Data
The study area was Lajeado Tacongava watershed ( Figure 1) which is located in the Northeast region of the Rio Grande do Sul State, Southern Brazil. With an area of 149.78 km 2 , the watershed is partially comprised of four cities with tradition in agriculture and livestock production: União da Serra, Serafina Correa, Montauri, and Guaporé. Farms with more than 50 fattening pigs and more than 1000 broilers were investigated. Criteria for farms selection were defined according to the environmental license classification [27] and economic viability. Figure 1 shows the location of farms in the watershed.   Between June 2019 and March 2020, we surveyed 78 farms (37 pig production and 41 poultry farming). Information related to farm-basic data (name, address, location, family workforce details, agricultural activities), environmental farm conditions (type of soil, water supply, water resources located on the farm, wastewater and solid waste management), and animal production features (current herd, the period of the production cycle, breeding systems, type of diet, initial and finishing weight) was collected.
The herd of pigs located on the farms within the Lajeado Tacongava watershed area ranged from 240 to 1150 head farm −1 cycle −1 , with a total herd of 55,071 head year −1 for pig. For broiler, it ranged from 7500 to 85,000 head farm −1 cycle −1 with a total herd of 6,108,600 head year −1 . Animal production is performed in a vertical integration system, where companies define initial and finishing animal weight according to the market demands. The first stages of the animals' life are neither performed on the farms nor in the watershed. Piglets are born on a farm where they spend 28 days in the weaning stage and 35 days in a farm nursery stage until achieving the weight for starting the finishing stage in the last farm [28]. In this paper, the first stage of a pig's life is called pig pre-chain. Broilers, as well, are born in a hatchery and after that, they are transported to the finishing farm. Figure 2 presents the stages, animal performance indicators, and water flows of the two animal production chains addressed in this study.
Water 2020, 12, x FOR PEER REVIEW 4 of 27 Between June 2019 and March 2020, we surveyed 78 farms (37 pig production and 41 poultry farming). Information related to farm-basic data (name, address, location, family workforce details, agricultural activities), environmental farm conditions (type of soil, water supply, water resources located on the farm, wastewater and solid waste management), and animal production features (current herd, the period of the production cycle, breeding systems, type of diet, initial and finishing weight) was collected.
The herd of pigs located on the farms within the Lajeado Tacongava watershed area ranged from 240 to 1150 head farm −1 cycle −1 , with a total herd of 55,071 head year −1 for pig. For broiler, it ranged from 7500 to 85,000 head farm −1 cycle −1 with a total herd of 6,108,600 head year −1 . Animal production is performed in a vertical integration system, where companies define initial and finishing animal weight according to the market demands. The first stages of the animals' life are neither performed on the farms nor in the watershed. Piglets are born on a farm where they spend 28 days in the weaning stage and 35 days in a farm nursery stage until achieving the weight for starting the finishing stage in the last farm [28]. In this paper, the first stage of a pig's life is called pig pre-chain. Broilers, as well, are born in a hatchery and after that, they are transported to the finishing farm. Figure 2 presents the stages, animal performance indicators, and water flows of the two animal production chains addressed in this study. Broilers and pigs were considered mixed between males and females. The mortality rate for broilers was considered 3.5% per cycle [29] and for pigs, it was null in the finishing stage. Feed compositions for fattening broilers (64% maize and 27% soybean meal) [30,31] and pigs (75% maize and 25% soybean) [32,33] were established for different stages of growth as shown in the supplementary material (Table S1). Pig pre-chain feed composition for weaning (75% maize and 25% soybean) and nursery (54% maize and 24.3% soybean) were also defined [28]. Companies provide animal feed to the farms and diets are specific for each stage of life. For this reason, the daily amount of feed and the type of diet were standardized for all farms (Table S1). Broilers and pigs were considered mixed between males and females. The mortality rate for broilers was considered 3.5% per cycle [29] and for pigs, it was null in the finishing stage. Feed compositions for fattening broilers (64% maize and 27% soybean meal) [30,31] and pigs (75% maize and 25% soybean) [32,33] were established for different stages of growth as shown in the supplementary material (Table S1). Pig pre-chain feed composition for weaning (75% maize and 25% soybean) and nursery (54% maize and 24.3% soybean) were also defined [28]. Companies provide animal feed to the farms Water 2020, 12, 3014 5 of 27 and diets are specific for each stage of life. For this reason, the daily amount of feed and the type of diet were standardized for all farms (Table S1). Table 1 shows the input data sources to calculate the water productivity indicators.

Calculation of Water Productivity
WP was calculated for broilers in the fattening stage (WP indirect+direct,broiler,Farm ) and pigs (WP indirect+direct,pig,Farm ) for pre-chain and fattening production stages accounting for direct and indirect water inflows, according to the FAO [1] and presented in Section 2.1.
The boundaries of the system were cradle-to-farm gate which involved feed and meat production at farm level. This study applied a bottom-up approach where farm specific data was used to calculate a detailed WP. In this study, WP was calculated in three references units in order to increase comparability with other studies. They were expressed on a mass basis (WP indirect+direct,broiler,Mass,Farm (Equation (1)); WP indirect+direct,pig,Mass,Farm (Equation (2))) (kg Carcass Weight (CW) m −3 ), food energy basis (WP indirect+direct,broiler,Energy,Farm (Equation (3)); WP indirect+direct,pig,Energy,Farm (Equation (4))) (MJ m −3 ) and monetary basis (WP indirect+direct,broiler,Mon,Farm (Equation (5)); WP indirect+direct,pig,Mon,Farm (Equation (6))) (R$ m −3 ).
Legend to above equations: WP indirect+direct,broiler,Mass,Farm is water productivity of chicken meat produced on mass base (kg Carcass Weight m −3 ); WP indirect+direct,pig,Mass,Farm is water productivity of pork meat produced on mass base (kg Carcass Weight m −3 ); WP indirect+direct,broiler,Energy,Farm is water productivity of chicken meat produced on food energy base (MJ m −3 ); WP indirect+direct,pig,Energy,Farm is water productivity of pork meat produced on food energy base (MJ m −3 ); WP indirect+direct,broiler,Mon,Farm is water productivity of chicken meat produced on a monetary base (R$ m −3 ); WP indirect+direct,pig,Mon,Farm is water productivity of pig meat produced on a monetary base (R$ m −3 ); Q direct+indirect,broiler,Farm is water consumption for broiler in fattening stage production (m 3 year −1 ); Q direct+indirect,pig,Farm is water consumption for pig production in fattening and pre-chain stages (m 3 year −1 ).
WP was also reported with fractions of green and blue water consumed (% blue water/% green water) as defined by the FAO [1]. Water consumption (Q indirect+direct,broiler,Farm (Equation (7)); Q indirect+direct,pig,Farm (Equation (8))) was calculated according to the consumed water, assigned to the generation of output for the stages assessed. Q indirect+direct,broiler,Farm = Q indirect,broiler,Feed + Q direct,broiler,Animal + Q direct,broiler,Housing Q indirect+direct,pig,Farm = Q indirect,pig,Feed + Q indirect+direct,pig,Animal + Q indirect+direct,pig,Housing (8) Legend to above equations: Q indirect_direct,broiler,Farm is the total water consumed for broiler purchased feed production + water consumed for broiler production (m 3 year −1 ) considering broiler fattening stage; Q indirect_direct,pig,Farm is the total water consumed for pig purchased feed production + water consumed for pig production (m 3 year −1 ) considering pig fattening and pre-chain stages; Q indirect,broiler,Feed is the total water consumed (evapotranspiration (ET); fresh matter (FM)) for purchased broiler feed production (m 3 year −1 ). It was based on the ratio of the yield of the field (cropland) for producing broiler feed and the ET from the field (from harvest of the previous crop through to harvest of the crop) [1]; Q indirect,pig,Feed is the total water consumed (ET; FM) for purchased pig (fattening and pre-chain stages) feed production (m 3 year −1 ). It was based on the ratio of the yield of the field (cropland) for producing pig feed and the ET from the field (from harvest of the previous crop through to harvest of the crop) [1]; Q direct,broiler,Animal is the total water consumed for broiler drinking (m 3 year −1 ); Q indirect+direct,pig,Animal is the total water consumed for pig drinking in pre-chain and fattening stages (m 3 year −1 ); Q direct,broiler,Housing is the total water consumed for services (cooling, cleaning) (m 3 year −1 ) for broiler production; Q indirect+direct,pig,Housing is the total water consumed for services (cleaning) (m 3 year −1 ) for a pig in pre-chain and fattening stages production.
Feed crop water productivity WP indirect,Feed of broiler WP indirect,broiler,Feed and pig WP indirect,pig,Feed was calculated according to Equations (9) and (10) Legend to above equations: WP indirect,broiler,Feed is water productivity of broiler feed consumption (kgFM m −3 ); WP indirect,pig,Feed is water productivity of pig feed consumption (fattening and pre-chain stages) (kgFM m −3 ).
Water demand for broiler feed produced for supplying pre-chain was neglected because broilers are transported to the farms a few hours after they were born. WP indirect,Feed was calculated for crops addressed to animal feed in two different crop rotations produced in three sites located in two crop producer regions ( Figure 3 Qindirect+direct,pig,Animal is the total water consumed for pig drinking in pre-chain and fattening stages (m 3 year −1 ); Qdirect,broiler,Housing is the total water consumed for services (cooling, cleaning) (m 3 year −1 ) for broiler production; Qindirect+direct,pig,Housing is the total water consumed for services (cleaning) (m 3 year −1 ) for a pig in pre-chain and fattening stages production.
Water demand for broiler feed produced for supplying pre-chain was neglected because broilers are transported to the farms a few hours after they were born.
WPindirect,Feed was calculated for crops addressed to animal feed in two different crop rotations produced in three sites located in two crop producer regions ( Figure 3): Center West, Primavera do Leste city (Mato Grosso State-MT) and South, Vacaria city (Rio Grande do Sul State-RS) and Cascavel city (Paraná State-PR). Throughout the paper, cities are abbreviated with RS for Vacaria, MT for Primavera do Leste, and PR for Cascavel.   The first crop rotation comprised soybean (fs-soy) and safra (fm-safra) intercropping with fallow and the same crops intercropping with ryegrass (safra: rm-safra; soy: rs-soy) where soy is produced one year and in the next year maize is produced in the same field driving to one crop harvest per year. Safrinha (smf-safrinha) harvested after soybean (smf-soy) intercropping with fallow compose the second crop rotation which represents two crops harvested per year (soy and maize). As suggested by Flach et al., [26], modeling water productivity of double cropping was applied in this study. The last crop rotation was not analyzed in Vacaria-Rio Grande do Sul state, because the climate characteristics are not favorable to safrinha production in the region.
Water consumed for animal feed production (Q indirect,Feed ) was calculated as actual evapotranspiration (ET) of the crops through the modeling software AgroHyd Farmmodel [20], which is based on the FAO's 56 dual crop coefficient method [44]. This method requires calculating (a) the reference evapotranspiration (ET0), (b) the potential crop transpiration (Tc), and (c) the actual transpiration (Tact) from the different datasets for climate, plants, and soil containing regional climate data, plant-specific parameters, and regional soil data. In the calculation presented in this study (d), the actual evapotranspiration of the individual crop is calculated as the sum of Tact and the actual Evaporation (Eact).
(a) With regional climate data ET0, a grass reference surface was calculated using the FAO Penman-Monteith equation. (b) To model Tc, the ET0 was adjusted for the individual crop with plant-specific parameters (e.g., the plant-specific basal crop coefficient (K cb )). Plant-specific parameters are provided in Table S2. (c) The calculation of Tact incorporates the effect of daily water stress due to water-limited conditions by linking the datasets on plants, soil, and climate on Tc. A water stress coefficient (Ks) incorporated water stress and reduced Tc to Tact. To determine the water stress coefficient (Ks), a simple tipping bucket approach was combined with regional soil and precipitation data. The equation for Tact (mm) applied here was: For adjustment on specific climatic conditions, the calculated K cb values were improved using the formula of K cb,adj where RH min is the minimum relative humidity, u 2 the wind speed at 2 m height (m s −1 ), and h the mean plant height during the mid or late season stage (m) for 20% ≤ RH min ≤ 80%. If the amount of soil water drops below a critical value, the crop is water-stressed [45]. To calculate the water stress coefficient, values of total available soil water in the root zone, readily available soil water in the root zone, and the root zone depletion are needed. Ks is given by: where Ks is the water stress coefficient (0-1), Dr is the root zone depletion (mm), TAW is the total available soil water in the root zone (mm), and RAW is the readily available soil water in the root zone (mm). The maximum value of Ks of 1 shows the absence of soil water stress. The total available soil water TAW (mm) can be calculated by the difference between water content at field capacity θ FC (m 3 m −3 ) and water content at wilting point θ WP (m 3 m −3 ). This value is multiplied by the effective rooting deep Zr (mm).
Water 2020, 12, 3014 9 of 27 Soil classification of USDA is available in the modeling system and clay soil was used. For this type of soil, the model considered a soil water content at field capacity of 0.396 (m 3 m −3 ) and soil water content at a wilting point of 0.248 (m 3 m −3 ) [44]. The readily available soil water content is described as p is a tabular value (Table S2) describing the average fraction of TAW that can be depleted from the root zone, without causing moisture stress for the crop. It can be adjusted with the formula In order to determine water availability for evapotranspiration, a root zone depletion Dr was calculated using a daily water balance using a simple tipping bucket approach: where Dri (mm)is the root zone depletion at the end of day i, Dri-1 (mm) the root zone depletion at the end of the previous day i − 1, Pi (mm) the precipitation on day i, Tacti (mm) the actual transpiration on day i, Ii the interception on day i (mm), and DPi (mm) is the water loss out of the root zone by deep percolation on day i. After heavy precipitation or irrigation, the soil water content in the root zone might exceed field capacity. The difference between the content which exceeded the field capacity and the soil water at field capacity is called deep percolation. Deep percolation is given by with Pi as precipitation on day i (mm), Ii for the interception on day I (mm), dPi for deep percolation on day i (mm), Dri-1 for water content in the root zone at the end of the previous day, i − 1 (mm), Iri for irrigation on day i, and Tacti for transpiration on day i. For the instant calculation of the values of DP and Dr for day i = 1 were approximated. The rainfall interception calculation used here was based on the work of von Hoyningen-Hüne [46] and Braden [47]. The approach was implemented in several agro-hydrological models of different complexity for the estimation in particular of the interception for agricultural crops, e.g., the physical-based model SWAP [48]. The authors measured the interception of precipitation for various crops. The general formula for canopy interception proposed is where I is the intercepted precipitation (mm), P is the gross precipitation (mm d −1 ), a is an empirical coefficient (mm day −1 ), and cf is the soil cover fraction (1 × 10 −0 . 385 LAI (-)). For increasing precipitation amounts, the amount of intercepted precipitation asymptotically reached the saturation amount a × LAI. We assumed a = 0.25 (mm day −1 ) for the agricultural crops. The calculation of the actual evapotranspiration of the individual crop ET(mm) was based on the previous three calculation steps (a, b, and c) of the actual transpiration (Tact) and the calculation of the actual evaporation (Eact) according to the procedures of the FAO 56 dual crop coefficient method [44] ET = Tact + Eact (20) The actual evapotranspiration (ET) was considered as water consumed for animal feed production (Q indirect,Feed ). A polygon with an area of 1 ha was defined as the modeling area for each crop analyzed and it was assumed that maize and soybean were produced at least one year before animal consumption (the year 2018).
Indirect water consumption (ET) was calculated considering the mean annual crop yield (CY) for maize and soy produced in each region (RS, PR, and MT) in a period of 10 years (2008-2018), obtained from the Brazilian Statistics Bureau [36], as shown in Table 2. The highest state mean annual crop yield attained for each region in 10 years was assumed to be the potential highest mean annual yield (HY) ( Table 2). The annual mean variation of CY and HY are shown in Table 2 with the standard variation (SE) and their variation along the years analyzed is shown in the supplementary material (Tables S3 and S4, respectively).
The sowing period was defined according to the Climate Risk Agricultural Zoning and the soil type was defined according to a normative instruction established by MAPA-Brazilian Ministry of Agriculture, Livestock, and Food Supply [49]. Soil types for all sites analyzed ( Figure 3) had a clay content higher than 35% [39,40], with a classification of clay soil [50]. Data on sowing and harvest and vegetation period are also shown in Table 2. Climate data, such as precipitation (mm), minimum and maximum temperature ( • C), sunshine hours (h), relative humidity (%), and wind speed (km h −1 ), was taken in monthly mean values in five climate stations from the INMET (National Institute of Meteorology, Brazil) located near the sites analyzed [37]. Data recorded on climate stations located in Mato Grosso State (Poxoréo), Paraná State (Campo Mourão; Maringá), and in Rio Grande do Sul State (Santa Maria; Vacaria) were uploaded in the modeling software. The distribution of precipitation per crop cycle had a large variation around the sites analyzed, where availability for summer crops (safra and soybean) ranged from 1218 mm year −1 to 2679 mm year −1 and for winter crop (safrinha), ranged from 432 mm year −1 to 994 mm year −1 [37]. The annual mean precipitation is presented in Figure 3.
Alvares et al. [38] showed that the analyzed feed crop production regions have different climate classifications according to the Köppen climate classification. Vacaria city, in RS, has a climate classification of Cfb (subtropical climate, without a dry season and with a temperate summer); Cascavel city, in PR, has a climate classification of Cfa (subtropical climate, without a dry season and with a hot summer) while Primavera do Leste city, in MT, has a climate classification of Aw (tropical climate, with a dry winter).
Water consumed for animal drinking (Q direct,broiler,Animal ; Q indirect+direct,pig,Animal ) (m 3 year −1 ) was calculated according to Palhares [41]. Water consumed for services (Q direct,broiler,Housing ; Q indirect+direct,pig,Housing ) accounted for water inflows used for cleaning pig [43] and broiler [42] facilities and for cooling, which was calculated only for broilers [42]. It was considered that the cooling system is activated whenever the temperature inside the stable is higher than 24 • C [29]. According to data on the maximum temperature in the watershed [37], the system was activated for 1420 h in 2019. The water consumption for animal drinking, cleaning, and cooling is shown in Table S5.

Scenarios Analyzed
In order to assess the effect of the implementation of potential agricultural best practices techniques (soil coverage and yields improvement) on the WP indirectly and feed for crops, four scenarios were analyzed: • 1st scenario (CY): WPindirect,Feed was calculated and assessed considering the current mean annual crops yield (CY) for maize and soy of each producer region analyzed (RS, PR, and MT) from 2008 to 2018 (Table 2; Table S3); • 2nd scenario (M50): lower soil evaporation of 50% was applied considering the same mean annual crop yield analyzed in the CY scenario; • 3rd scenario (HY): the highest state mean annual crop yield ( Table 2, Table S4) for maize and soy was applied to the respective city located in the same state for analyzing the impact of crop yields potential improvement in the WP indirect,Feed ; • 4th scenario (HY50%): a restriction of 50% in the soil evaporation was applied for the same highest state mean annual crop yield analyzed in the HY scenario.
The restriction in soil evaporation proposed in the 2nd scenario (M50) and 4th scenario (HY50%) was inserted in the modeling software AgroHyd Farmmodel. The potential water-saving for producing animal feed (Q indirect,Feed ) demanded to pig and broiler feed in the whole watershed was estimated by comparing the results of the 1st scenario with the results of the 2nd scenario, 3rd scenario, and 4th scenario (agricultural best practices techniques). Water demand for the production of diet ingredients such as minerals, vitamins, enzymes, etc, was not included in the account for both broiler and pig production, since it was assumed to be negligible.

Sources of Uncertainty
Uncertain information of input data needs to be carefully evaluated, as it is often highly uncertain due to variability and lack of measured data. The Monte Carlo simulation was applied for uncertainty data analyses, as suggested by the FAO [1], considering the annual mean evapotranspiration (m 3 year −1 ) generated between 2008/18. Ten simulations were run for each year (from 2008 to 2018), resulting in 110 simulations for each crop produced in each crop rotation (fm-safra; fs-soy; rm-safra; rs-soy; smf-safrinha; smf-soy) appreciating the analyzed production regions (RS, PR, MT) respective broiler and pig feed. Simulations were facilitated through the Dell ® Statistic 13 software.

Feed Crop Water Productivity Scenarios and Potential Water Saving
Results of feed crop water productivity of maize (safra, safrinha) and soy produced in the analyzed regions and used to feed broilers (WP indirect,broiler,Feed ) and pigs (WP indirect,pig,Feed ) produced in the Lajeado Tacongava watershed for the 1st scenario (CY) are shown in Table 3. Crop water productivity of maize and soy produced in the analyzed regions is shown in Table S6. Drastig et al. [42] assessed water productivity of maize and soy used for feeding broilers on Brazilian farms, considering transpiration to calculate the water demand for feed production in a dry matter (DM) base. found in this study for Primavera do Leste, which is located in the same State (MT) considering water demand for transpiration and DM base. The same plant parameters for both maize types (safra and safrinha) were applied for Drastig et al. [42] which found a WP indirect,Feed of 1.40 kgDM m −3 for safrinha produced in MT. Specific plant parameters were applied in this study (Table S2) The highest WP indirect,Feed of safra arrangements were found in PR, which were 54.7% (fm-safra) and 44.1% (rm-safra) higher in comparison to the lowest values of safra found in RS for both animal feed production. PR also attained the highest WP indirect,Feed of soy (rs-soy), which were 12.7% and 17.3% higher than the lowest value found in RS for broiler and pig feed production, respectively. Among all safra and soy arrangements analyzed, only WP indirect,Feed of soy (fs-soy) driven to broiler feed was the highest in MT (9.3%) than in PR. Despite both production regions (PR and RS) being located in the south of Brazil, the highest WP indirect,Feed found in PR may result from a combination of favorable climate conditions and the implementation of best crop management techniques that provide an efficient rainfall water use.
In the second maize crop season (safrinha), which has a longer vegetation period of 140 days (Table 2), MT achieved a 10.9% higher yield of safrinha (smf-safrinha) than PR, with a 16.7% lower mean precipitation and lower ET for broiler (12.3%) and pig feed (18.7%) production. Under these conditions, MT attained the highest values of WP indirect,Feed of safrinha for both animal feed productions analyzed (Table 3). During the last decade , the State of MT, which currently stands out as the larger safrinha producer in Brazil, improved safrinha yield by 58%, with a mean yield from 4.0 to 6.4 t ha −1 , while the State of PR attained the improvement of 24% with a mean yield from 4.9 to 6.0 t ha −1 [53]. In order to improve WP indirect,Feed of safrinha in PR, best crop management techniques should be implemented as well as the selection of crops with better genetics for greater water use efficiency. Nascimento et al., [54], hybrids of safrinha have been developed allowing early sowing, more efficiency for water consumption, light, temperature, and nutrient absorption.
In a comparison between maize crops (safrinha: smf-safrinha; safra: fm-safra) produced in MT and PR, the safra mean yield was 65% higher than the safrinha mean yield during the analyzed period. Even with a lower yield than safra, currently, safrinha is responsible for 65% of the Brazilian maize national production and took it to center stage as the main maize crop in Brazil [55]. Safrinha attained a mean higher WP indirect,Feed of 41.2% (PR) and 41.8% (MT) than safra for broiler and pig feed, respectively (Table 3), with a lower median ET ratio (Table S7). One agricultural practice that could increase safrinha yield is the use of irrigation, in which water productivity values could be assessed in further studies.
The future intention of the authors is to build a database of water indicators from farms that are using different farming measures to increase livestock water productivity. That will allow us to evaluate the differences in water productivity between different farm systems. In order to implement agriculture best practices for raising water productivity, factors such as climate conditions, soil water availability, crop quality, and cropping system, should be carefully analyzed. In Brazil, due to the tropical and subtropical climate conditions and rainfed agriculture, soil conservation techniques are key for enhancing the storage of rainfall water [18]. Among all scenarios analyzed, the improvement of crop yield and the reduction of soil evaporation analyzed in the 4th scenario (HY50%), resulted in the highest WP indirect,Feed values for crops produced in all crop rotations and among all scenarios and producer regions analyzed (Table 4). In a higher crop yield scenario assessed (3rd scenario (HY)), WPindirect,Feed of all maize (safra and safrinha) and soy crop rotation arrangements were improved (Table 4) in comparison to the values found in the 1st scenario (CY) ( Table 3). Results of WP indirect,Feed which considered a lower soil evaporation of 50% (2nd scenario (M50%)), which could be attained with the mulching of forage crops or residues of the main crops that remain in the soil after harvesting, also demonstrated the possibility of attaining a higher efficient rainfall water use in all crop rotation and producer regions analyzed (Table 4).
WP indirect,Feed of safrinha (smf-safrinha) and soy (smf-soy) produced in PR was improved by 84% and 43%, respectively, for both animal feed productions in the 4th scenario (HY50%) In PR, the WP indirect,Feed of safra (fm-safra, rm-safra) could be improved by 19.4% and 21.4% for broiler and pig feed production. All safrinha (smf-safrinha) and soy (smf-soy) arrangements analyzed attained higher WP indirect,Feed in MT in the 2nd scenario (M50%), with the highest WP indirect,Feed improvement value found for safrinha (smf-safrinha) which was 18% higher than in the 1st scenario (CY) for both animal feed production.
WPindirect,Feed of safra produced in RS intercropping with fallow (fm-safra) and used for broiler (86%) and pig (87%) feed, and safra intercropping with ryegrass (rm-safra) used for both animal feeds (59%), had the highest improvement between the 4th scenario (HY50%) and the 1st scenario (CY). The lowest WPindirect,Feed of safra (fm-safra; rm-safra) and soy (fs-soy; rs-soy) among all crops arrangements, producer regions and the four scenarios analyzed was found for crops produced in RS for both animal feed production (Tables 3 and 4). Considering forage crops like ryegrass, for instance, can present high water consumption throughout their growth cycle [56], another intercropping species with lower water demands (lower single crop coefficient K c ) could be sown in order to improve WP indirect,Feed of crops produced in PR and RS.
An impressive amount of water-saving (m 3 year −1 ) could be attained depending on the crop rotation, production region, and implementation of best practices. Table 5 shows the potential water-saving (m 3 year −1 ) for producing animal feed (Q indirect,Feed ) demanded to pigs and broilers.  Water savings were calculated as the volumetric difference between water consumed in scenario CY and scenarios (M50%, HY, and HY50%). Higher water-saving could be attained with the implementation of soil evaporation restriction and higher crop yields in the 4th Scenario (HY50%), mainly for crop rotation composed by safra (fm-safra) and soy (fs-soy) produced in RS, with 45% and 46% of potential water-saving for broiler and pig feed production, respectively. The amount of water-saving was higher for crop rotation composed of safra and soy (fm-safra; fs-soy) than safrinha and soy (smf-safrinha; smf-soy). The water scarcity footprint in a watershed level is mainly related to blue water availability, but green water is the unique source to recharge the blue waters. The amount of potential rainfall water saved (green water) will be incorporated in the watersheds by infiltration and run-off, contributing to groundwater and superficial reservoirs recharge as well as improving water storage in the soil.
In MT, for instance, the water-saving for crop rotation composed by safra and soy (fm-safra; fs-soy) was more than 88% higher than water saving for safrinha and soy (smf-safrinha; smf-soy) for both animal feed productions. In a study addressing rainfed crop production in Brazil, the Brazilian National Agency [17] stated that the Center-West region had a mean water deficit, which means the amount of water demand for rainfed agriculture which exceeds the rainfall water availability of 34% (26% in the first stages of crop development and 8% in the last crop development stage) between 2013-2017 mainly for maize and soy production. For the same period, rainfed crops produced in the south of Brazil had a mean water deficit of 18%. In a ranking of the cities with the highest rainfall consumption and water deficit in the year 2017, the production region analyzed in MT (Primavera do Leste city) was in the fourteenth position with a total rainfall consumption of 40.8 m 3 s −1 and a water deficit of 27%, mainly for soybean production [17]. Best practices could be implemented for higher water use efficiency and consequently, higher water productivity.

Water Input and Water Productivity of the Chicken Meat and Pig Pork
Water input for broiler (Q indirect+direct,poultry,Farm ) and pig (Q indirect+direct,pig,Farm ) production in the 1st scenario (CY) is shown in Table S8. Water demand for animal feed production accounted for more than 99% of all water inputs for broiler and pig production considering crop rotation arrangements produced in RS, PR, and MT. Ridoutt et al. [57] showed that the production of animal feed can be considered as one of the major environmental impacts of livestock production. In Germany, water demand for poultry feed production accounted for 90% to 93% [52] while in Tunisia it accounted for 95.7% of all water inputs for broiler production [58]. Water required for pig feed in this study was higher than found by Nakamura and Itsubo [59] and Bai et al., [60], where it accounted for 90% of the entire water input. The sustainability of the broiler and pig production chain is highly dependent on the improvement of rainfall use efficiency (green water).
Broiler and pig industrial production systems are characterized by producing animal feed off the farm. Consequently, the increase of feed water productivity of these products is not in the hands of farmers, since the food companies are the ones who choose breeds and determine ration composition and sourcing. Although Brazil is one of the largest world grain producers, climate and soil conditions, crop yield, transportation, and marketing are key variables for supporting companies' decision-making regarding the definition of the maize and soy source for supplying broiler and pig feed production.
For broilers, water required for animal drinking, cooling, and cleaning accounted for less than 1% of the total amount of total water demand. These three consumptions represented 57.3%, 42.7%, and <0.02% of water demand, respectively (Table S8). A similar water input pattern for Brazilian poultry production was found by Palhares [61], from the years 2000/10, and Drastig et al. [42]. In both studies, water inputs (drinking, cooling, and cleaning) accounted for less than 1%. Krauß et al. [52] found that drinking and cleaning water accounted for less than 1% in German broiler farms.
The same pattern was found for pig production, with water required for animal drinking and cleaning accounting for less than 1% of the total water demand (Table S8). Water demand for cleaning represented 55.1% while drinking water represented 44.9%. Palhares [62] found similar results for pigs produced in the State of Rio Grande do Sul, where drinking water accounted for less than 1%.
Even though it accounts for less than 1%, water required for drinking, cooling, and cleaning is withdrawal from surface or groundwater water sources that could be overexploited and limited dependent on their availability in the watershed. Livestock drinking water represents 8% (166.8 m 3 s −1 ) of the total water demand in the country while it represents 16% in Southern Brazil [63]. Considering Brazil is one of the largest pig and broiler producers in the world, and the number of animals raised every year, the amount of water demand for these activities is highly relevant. According to Palhares et al. [64], there is a need for recognizing the importance of improving the water productivity of livestock products to mitigate the negative impacts on water availability. Another important fact is that Brazilian farms have to pay to use water sources on their farms. This charge has not yet been fully implemented in the whole country but it is a process in progress [65] and shall contribute to farmers enhancing water productivity.
Water productivities of broiler and pig production considering crops produced in the 1st scenario (CY) are shown in Table 6. WP indirect+direct,poultry,Mass,Farm of poultry chicken is higher than WP indirect+direct,pig,Mass,Farm of pig pork. It is related to a lot of genetic and non-genetic aspects, but one of the most important aspects is the higher feed conversion efficiency (kg feed to produce kg meat) of broilers. The poultry production system has the longest history among livestock activities of using knowledge and technologies to improve production efficiency. This makes these animals the best feed-to-food converters, which will have a positive impact on increasing water productivity. Mekonnen et al. [66] highlighted that the results of this kind of study should be compared with values from other studies at three levels: feed conversion efficiency, total feed requirement, and the WP of livestock products.
The water productivity is too variable even for the same animal product because several factors contribute to the final value. Two of the most important are crop productivity and the performance indexes of the animal. Therefore, this great variability is expected and this highlights the importance of local assessments, considering data from each farm and watershed. This is the only way to propose practices and policies according to local conditions. It will result in water efficiency improvements and reduce pressure on water availability.
In terms of revenues (WP indirect+direct,broiler,Mon,Farm ; WP indirect+direct,pig,Mon,Farm ), considering the price charged for animal carcass in Rio Grande do Sul State, it could be possible to raise incomes by 0.94 R$ m −3 and 0.64 R$ m −3 if pig and broiler feed were produced in MT (smf-safrinha; smf-soy) rather than in RS (rm-safra; rs-soy), respectively.
As shown in Figure 4 (Table S9) and Figure 5 (Table S10), water productivity can be improved for mass, energy, and revenues if the best crop management and conservation agriculture techniques were implemented and improved for animal feed production. Water input of broilers and pigs in the farms investigated in the 2nd Scenario (M50%), 3rd Scenario (HY), and 4th Scenario (HY50%) are shown in Tables S11 and S12, respectively. Table 6. Water productivity of broiler chicken production (WP indirect+direct,poultry,Farm ) in the farms investigated in the 1st scenario (CY). RS: Vacaria (city) located in Rio Grande do Sul State; PR: Cascavel (city) located in Paraná State; MT: Primavera do Leste (city) located in Mato Grosso State.

RS
WP indirect+direct,Mass,Farm varied according to the crop rotation system applied for animal feed production. The highest values of WP indirect+direct,broiler,Mass,Farm (0.53 kgCW m −3 ) and WP indirect+direct,pig,Mass,Farm (0.35 kgCW m −3 ) were found in MT for crop rotation composed of safrinha (smf-safrinha) and soy (smf-soy), which were 139% and 132% higher, respectively, than the lowest values considering animal feed produced in RS (rm-safra; rs-soy). An amount of water-saving (Q indirect,Feed ) of 0.0336 km 3 year −1 and 0.0202 km 3 year −1 could be attained if broilers and pigs would be fed with safrinha (smf-safrinha) and soy (smf-soy) produced in MT rather than safra (rm-safra) and soy (rs-soy) produced in RS. This potential amount of water-saving represents 0.07%, 0.06%, and 0.09% of the total water demand for rainfed maize and soy produced in the States of Paraná (PR), Mato Grosso (MT), and Rio Grande do Sul (RS), respectively, for broiler feed production and 0.04%, 0.04%, and 0.05% for pig feed production in the same states [17].
The highest water productivity found by Drastig et al. [42] for Brazilian poultry farms was 0.33 kgCW m −3 . Considering the same water demand conditions for feed production, it was close to WP found in this study considering feed (rm-safra; rs-soy) produced in RS, 0.30 kgCW m −3 and in PR, 0.35 kgCW m −3 . A higher WP of 0.87 kgCW m −3 was found in this study considering animal feed produced in MT (smf-safrinha; smf-soy), which was 163.3% higher than the highest WP found by Drastig et al. [42]. It reinforces the importance of crop rotation and the source of feed (producer region) in the WP of animals produced in the studied watershed.
In Germany, Krauß et al., [52] found a similar WP of 0.30 kgCW m −3 while Renault et al., [67] found a lower WP of 0.24 kgCW m −3 in California, USA. The last lower WP value could be associated with different yield index and crop management techniques applied in the 1990s in the US. The WP of 0.29 kgCW m −3 found in the USA in the period of 2014/16 [66] was almost the same as found in this study for animal feed produced in PR (rm-safra; rs-soy) of 0.28 kgCW m −3 and 45.3% lower than the highest WP of 0.53 kgCW m −3 found for broiler feed produced in MT (smf-safrinha; smf-soy).
In California, USA, Renault et al. [67] found a pig WP of 0.21 kg CW −1 m −3 while Mekonnen et al. [66] found 0.25 kg CW −1 m −3 for the USA, which was close to the mean value (0.22 kgCW m −3 ) found in this study considering all crop rotations and grain production regions analyzed.
Among the analyzed production regions, the state of MT is the most threatened by the potential impacts of global warming [68,69]. In a business-as-usual scenario, the state's agricultural production would drop significantly causing a great increase in the water demand for products that depend on the state crops. Another point to be considered is emphasized by Zhuou et al. [70], i.e., the livestock industry expansion will undoubtedly increase feed production demand and associated water consumption, thus increasing local water scarcity levels. MT is located more than 1800 km away from the farms analyzed in this study and grains should be transported via a highway system. The transport sector is one of the major contributors to global CO 2 emission where Brazil stands out among the top ten transport CO 2 emitters in the world [71]. Thus, it is necessary to assess other environmental potential impacts for supporting decision-making in order to define if MT is the best region for animal feed production for supplying broiler and pig feed demands on farms located in the south among the producer regions analyzed. Animal production in Brazil plays an important social and economic role and investments in researches for improving water productivity of rainfed crops used for animal feed under sustainable crop management techniques are essential.
Results from this study showed that independent of the crop rotation system, water demand for animal feed production has the same percentage of total water demand to produce the livestock product even with differences in structural and input parameters applied to each crop rotation arrangement analyzed. Likewise, results also demonstrated the high impact of the feed water productivity on the WP of an animal product. Therefore, in order to attain more animal protein per drop of water, the livestock sector must work closely with the agriculture sector and best practices should be implemented in order to attain a higher rainfall water productivity. Nakamura and Itsubo [59]; Vellenga et al. [72]; Palhares et al. [73] emphasize that the greatest potential to increase water productivity and efficiency in livestock is proliferous feed production and precise nutritional animal management.

Uncertainty Data Analyses
Results for the uncertainty data analyzed considering evapotranspiration of crops used for animal feed are shown in Figure 6  social and economic role and investments in researches for improving water productivity of rainfed crops used for animal feed under sustainable crop management techniques are essential.
Results from this study showed that independent of the crop rotation system, water demand for animal feed production has the same percentage of total water demand to produce the livestock product even with differences in structural and input parameters applied to each crop rotation arrangement analyzed. Likewise, results also demonstrated the high impact of the feed water productivity on the WP of an animal product. Therefore, in order to attain more animal protein per drop of water, the livestock sector must work closely with the agriculture sector and best practices should be implemented in order to attain a higher rainfall water productivity. Nakamura and Itsubo [59]; Vellenga et al. [72]; Palhares et al. [73] emphasize that the greatest potential to increase water productivity and efficiency in livestock is proliferous feed production and precise nutritional animal management.

Uncertainty Data Analyses
Results for the uncertainty data analyzed considering evapotranspiration of crops used for animal feed are shown in Figure 6  The confidence interval demonstrates the higher variation of evapotranspiration ratio for all crops used for animal feed and crop rotations produced in the producer regions throughout the period analyzed (2008-2018), which demonstrated a high value of uncertainty. Climatic data hold inherent uncertainty due to an overlapping accumulation of uncontrollable sources of variation. Evapotranspiration, for instance, is influenced by a considerable number of uncontrolled random factors, as precipitation, temperature, solar radiation, wind, humidity.
The evapotranspiration ratio was calculated based on climate data which are monitored by climate stations managed by the INMET [37]. Potential breaks for station maintenance or due to malfunctions may have occurred over the period analyzed and might have affected the quality of The confidence interval demonstrates the higher variation of evapotranspiration ratio for all crops used for animal feed and crop rotations produced in the producer regions throughout the period analyzed (2008-2018), which demonstrated a high value of uncertainty. Climatic data hold inherent uncertainty due to an overlapping accumulation of uncontrollable sources of variation.
Evapotranspiration, for instance, is influenced by a considerable number of uncontrolled random factors, as precipitation, temperature, solar radiation, wind, humidity.
The evapotranspiration ratio was calculated based on climate data which are monitored by climate stations managed by the INMET [37]. Potential breaks for station maintenance or due to malfunctions may have occurred over the period analyzed and might have affected the quality of data collection. Missing data on climate grids were automatically corrected by the modeling software AgroHyd Farmmodel for ET calculation and may also represent an inaccurate value.
Input data related to animal feed composition and water flows of drinking and housing were standardized for all farms, as explained in Section 2. Water meters, for instance, used for accounting water flows demanded by animal drinking and housing, are unfortunately rare in the area of study. Hence, it is not possible to have an accurate flow of water used for supplying these demands for each farm. A similar situation should be highlighted related to animal diets. Diets are provided by the companies to the whole farms and formulations are not published. Hence, it is essential to highlight the results found in this study have uncertainty.

Applications of the LEAP-Guidelines
The guidelines published by the FAO [1] to assess water use at the farm-scale by the indicator water productivity have been successfully applied in the first Brazilian case study. This will strengthen further development and adoption of harmonized scientifically robust and practical frameworks for the assessment of water use in livestock production systems and supply chains worldwide. Research is needed for further development and application of the guidelines including calculation of the water scarcity footprint (WSF). According to the guidelines published by the FAO [1], the overall water productivity metric of a production system incorporating indirect water use shall be accompanied by the water scarcity footprint of the analyzed system. Due to the fact that the main focus in this study was to improve the insight about the green water demand for maize and soy rainfed crops, the water scarcity footprint (WSF) was not analyzed. As a next step, the analysis of water productivity and the water scarcity footprint in the Lajeado Tacongava Watershed in Brazil is planned, which will justify the watershed as the study area used in this study.
The use of direct water use vs. direct + indirect flows for the calculation of water productivity has its pros and cons. The use of direct water helps to assess yield or economic return from a farm's livestock production to assist farmers to understand the water flows in their farms and to optimize water use by agronomic measures and farm management at one specific farm location. The inclusion of direct + indirect water flows helps to assess yield or economic return from a farm's livestock production to assist farmers to understand the water flows in their farms and the resulting effects of optimizing water use through agronomic measures and farm management at the specific farm location and in potentially different regions [1]. From the results presented here, it is suggested that the calculation of the water productivity of the pig and poultry supply should include direct and indirect water. It is important to understand the water flows in the different regions from which the feed stems because 99% of the water demand for broiler and pig production in the study area stems from the purchased feed in different regions.
The water productivity method can be applied in the assessment of farming measures to improve rainwater productivity in livestock production. The development of a consistent and coherent approach to assessing water productivity in livestock supply chains could lead to a benchmark for water productivity values for specific species and production intensities in specific regions.
Furthermore, research is needed for further development and application of the methodology including: • To analyze the temporal variation in the water productivity indicators and their range, depending on the farming system.

•
To investigate the effectiveness of single and combined measures of farmers to improve water productivity and the water scarcity footprint and • To analyze the uncertainty considering evapotranspiration combined with the uncertainty of farm-basic data, environmental farm conditions, animal production features, and complementary data, such as maximum animal housing and size of the barns.

Conclusions
Water productivity of broiler chicken and pig pork produced in the Lajeado Tacongava watershed is directly dependent on rainwater productivity for animal feed production which varies according to the crop producer region, crop rotation, maize season crop (safra; safrinha), and agronomic techniques applied (e.g., no-tillage system). The results show that the implementation of best agriculture practices for animal feed production allows for higher rainfall water-savings in comparison to a conventional agriculture system, which has a positive impact on the water productivity of chicken and pork meat. Water productivity of the animal production chain in tropical regions demands a close relation to agriculture in order to attain a better understanding and improvement of rainfall water productivity for animal feed production. Besides, water productivity of animal feed production and water consumed in the farm-livestock sector should also take into account the water availability of the basin to support the decision-making in order to reach a more sustainable animal production chain. Water productivity assessment is also important for supporting government decision-making considering the economic and social importance of animal production and its dependence on water, which is a limiting factor. Management and best practices for crop production regions, crop rotation arrangements, and different types of soil can be assessed in further studies to improve the insight about rainfall water efficiency in order to attain higher water productivity of crops which will reflect in the water performance of animal products. Considering the direct relation between agriculture and livestock, climate change predictions should also be taken into account for future studies. The guidelines published by the FAO [1] can be applied in the assessment of farming measures to improve rainwater productivity in livestock production Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4441/12/11/3014/s1, Table S1: Maize and soybean meal rations and feed intake for WP calculations, Table S2: Plant-specific parameters for the crop-related modeling, Table S3: Current annual mean yield (t ha −1 ) of each city (2008-2018).), Table S4: Highest annual mean yield (t ha −1 ) of each State (2008-2018)., Table S5: Water consumption (m 3 head −1 cycle −1 ) considered for animal drinking, cleaning, and cooling., Table S6: Crop water productivity of maize (safra, safrinha) and soy produced in the Brazilian sites analyzed (2008/18) in the 1st scenario (CY), Table S7: Median values of ET of maize and soy produced in the Brazilian regions analyzed used for broiler and pig feed (2008-2018) in the 1st scenario (CY), Table S8: Water input (Q indirect+direct,Farm ) (m 3 year −1 ) of broiler chicken and pig pork production in the farms investigated in the 1st scenario (CY), Table S9: Water productivity of broiler chicken production (WP indirect+direct,poultry,Farm ) in the farms investigated considering the 2nd scenario (M50%), 3rd scenario (HY), and 4th scenario (HY50%), Table S10: Water productivity of pig pork production (WP indirect+direct,pig,Farm ) in the farms investigated considering the 2nd scenario (M50%), 3rd scenario (HY), and 4th scenario (HY50%), Table S11; Water inputs (Q indirect+direct,Farm ) (m 3 year −1 ) of broilers in the farms investigated in the 2nd Scenario (M50%), 3rd Scenario (HY), and 4th Scenario (HY50%), Table S12: Water inputs (Q indirect+direct,Farm ) (m 3 year −1 ) of pigs in the farms investigated in the 2nd Scenario (M50%), 3rd Scenario (HY), and 4th Scenario (HY50%).