This study aimed to modeling and bioclimatic zoning of ‘Saanen’ goat milk production in northeastern Brazil. Data were obtained from a 14-year dairy control program for 246 ‘Saanen’ goat matrices. Dairy production meteorological influence and mathematical modeling were assessed through statistical analyses such as Pearson’s correlation matrix and path analysis. Bioclimatic zoning for current and future scenarios, with climate changes, were delineated using inverse-square distance (ISD) and spherical model (kriging) interpolations. Results showed maximum air temperature to contribute substantially to goat dairy production; therefore, its modeling was based thereon. Current and future dairy production maps highlighted Maranhão and Piauí as regions of limited ‘Saanen’ milk production. Furthermore, we should emphasize that temperature rises in further scenarios may have significant impact on the potential of dairy farming for the region under study.
