Effets des facteurs environnementaux sur la production laitière des vaches Holstein dans la région de Souss-Massa au Maroc

Abdslem Talbi et Said El Madidi

Laboratoire de Biotechnologie et Valorisation des Ressources Naturelles, Facultés des Sciences, Université IBN-ZOHR,
Cité Dakhla, BP 8106, Agadir, Maroc.
Talbi_abdslem@yahoo.fr

Résumé

L’objectif de cette étude est l’évaluation de la variabilité des données du contrôle laitier de 9 194 lactations de 6 307 vaches Holstein de 118 troupeaux contrôlés de 2008 à 2012 dans la région de Souss-Massa, au sud-ouest du Maroc. Les températures moyennes sont de 19°C en hiver et de 27°C en été et la pluviométrie annuelle varie entre 180 et 280 mm par an. Afin d’étudier la variation de la production laitière en fonction de facteurs environnementaux, les données ont été analysées à l’aide d’un modèle linéaire général incluant les effets fixes: rang de lactation, âge au vêlage, année et saison de vêlage.

Les moyennes de la quantité de lait (QL_305j), de la matière grasse (QMG_305j) et du taux butyreux (TB_305j) enregistrées en lactation de référence (en 305 jours) sont égales respectivement à 6 578 ± 1 619 kg, 241 ± 62 kg et 3,7 ± 0,53 %. Le rang de lactation, l’année et la saison de vêlage ont des effets très hautement significatifs (P<0,001) sur (QL_305j) et (QMG_305j). Les saisons de vêlage favorables sont l'hiver et surtout l'automne.

Mots-clés: année de vêlage, rang de lactation, quantité de lait, saison de vêlage

Effects of environmental factors on milk production of Holstein cows in Souss-Massa region of Morocco

Abstract

The objective of this study is to evaluate the variability of a total of 9 194 lactation records of 6 307 Holstein cows from 118 herd kept in the locality of Souss-Massa in southwest of Morocco and controlled between 2008 and 2012. The mean temperatures are 19°C in winter and 27°C in summer and annual rainfall ranges between 180 and 280 mm per year.

The data were analyzed using the least squares approach with a general linear model to evaluate significant effects of lactation number, age at calving, year and season of calving on different milk yield traits. The overall mean of milk yield, fat yield and fat content of Holstein cows in the present study was 6 578 ± 1 619 kg, 241 ± 62 kg and 3,7 ± 0,53 %, respectively. The effects of lactation number, season and year of calving on 305-days milk yield and fat yield were highly significant (P<0,001). Moreover, the suitable calving seasons of milk production are winter and especially autumn.

Keywords: calving season, calving year, lactation number, milk yield

Introduction

Parmi les races pures à hautes performances laitières présentes au Maroc, la race Holstein occupe une place importante et prédominante. Cette race a dominé la production mondiale de lait pour les 25 années passées (Heins et al 2012) et elle est considérée comme la plus productive dans le monde. Cependant, les performances de la vache Holstein varient considérablement dans différentes conditions environnementales. Dans les zones tempérées, le rendement laitier de la race Holstein dépassait largement les 7 tonnes en 305 jours (Albarrán-Portilloa et Pollottb 2011, García-Ispierto et al 2007, González-Recio et al 2004, Knaus 2009, La Bonnardière et al 2007, Pirlo et al 2000, Riecka et Candrák 2011, Wicks et Leaver 2006); il a atteint 9276 kg en 305 jours pour les Prim’Holstein en France au contrôle laitier en 2013 (Prim’Holstein France 2015); mais en régions tropicales et subtropicales la production moyenne est estimée entre 3 et 4 tonnes (Usman et al 2013). Alors que les Vaches Holstein élevées dans les conditions de l'Afrique du Nord produisaient de 5 à 6 tonnes (Ajili et al 2007, Boujenane 2002; 2013, Boujenane and Aïssa 2008, M’hamdi et al 2012).

Toutefois, le rendement laitier d’un élevage bovin est le résultat de la combinaison du génotype et des facteurs non-génétiques (Çilek 2009), puisque l’interaction environnement- génotype joue un rôle important dans l'expression de la valeur génétique complète des bovins laitiers (Cienfuegos-Rivas et al 2006, Hammami et al 2009, Ravagnolo et Misztal 2002, Scholtz et al 2010), d’où l'expression du mérite génétique varie d'un état environnemental à un autre et elle est fortement influencée par des facteurs non-génétiques (Javed et al 2001). La production de lait est influencée par le potentiel génétique de l’animal et les facteurs non-génétiques (Albarrán-Portilloa et Pollott 2011); la connaissance de ces facteurs est essentielle pour une gestion efficace et une estimation précise des valeurs d'élevage (Nyamushamba et al 2014, Nyamushamba et al 2013).

Cette étude a été menée pour évaluer la variabilité de la quantité de lait en 305 jours, la quantité de matière grasse et le taux butyreux chez des bovins Holstein élevés dans le sud-ouest du Maroc avec un climat semi-aride et pour étudier les effets de certains facteurs environnementaux sur leurs performances laitières afin de se servir de ces données pour l’amélioration de la gestion des troupeaux laitiers.

Situation géographique et climat

La région de Souss-Massa est située au sud-ouest du Maroc et couvre une superficie de l’ordre de 25 000 km² réparties entre 25 % de plaine et 75 % en montagne. Elle comprend deux principales plaines qui sont la plaine du Souss (4 150 km²) et la plaine de Chtouka-Massa (1 260 km²) (MEMEE 2015), d’altitude comprise entre 0 et 700 m et elle occupe une cuvette entre les montagnes du Haut Atlas du Nord, Anti-Atlas du Sud et l'Océan Atlantique de l’ouest. Sur le plan administratif, la zone du Souss Massa couvre la préfecture d'Agadir Ida-Outanane, la préfecture d'Inezgane-Ait-Melloul, la Province de Taroudant et la Province de Tiznit (MEMEE 2015).

Le climat est semi-aride à aride (Bouchaou et al 2011, Tagma et al 2009), caractérisé par un climat littoral doux dans l'Ouest et semi-continental chaud dans l'Est (Tagma et al 2009). La région est soumise à l'influence des vents froids de l'océan Atlantique (courant des Canaries) et l'influence des vents chauds sub-sahariens (Manssouri et al 2014). Les températures moyennes annuelles varient de 14°C sur le Haut-Atlas au nord à 20°C sur l’Anti-Atlas au sud. Ainsi les stations d’Agadir, Barrage-Abdelmoumen, Barrage Aoulouz, Issen, Taroudant et Barrage Youssef Ben Tachfine ont enregistré, respectivement, 18,3; 20,6; 19,9; 19,3; 19,7 et 19,9°C (ABHSM 2015) par contre les températures moyennes dans la plaine de Souss-Massa sont relativement élevées (19°C en hiver et 27°C en été) et les températures journalières minimale et maximale atteignent respectivement -3°C et 48°C (Bouchaou L et al 2008, Bouragba et al 2011).

Les précipitations sont en général caractérisées par leur faiblesse, leur irrégularité et leur mauvaise répartition allant de 180 mm / an sur la plaine à 600 mm / an dans les montagnes (altitude > 700 m) (Bouchaou et al 2011). Sur les plaines les précipitations diminuent du nord au sud et de l’ouest vers l’est (ABHSM 2015) enregistrant des moyennes de l’ordre de 280 mm sur la plaine de Souss, 265 mm sur la plaine de Massa et 180 mm à Tiznit au sud (MEMEE 2015, ABHSM 2015).

L'économie de la Souss-Massa est principalement basée sur l'agriculture, la pêche maritime et le tourisme. La majorité des activités agricoles sont basées sur l’irrigation avec un potentiel en terres irrigables de l’ordre de 250 000 ha. La superficie irriguée s’étend sur (143 310 ha), située principalement sur les plaines de Souss (114 645 ha) et Massa (27 105 ha), dont près de 60 % est irriguée de manière moderne. L’élevage représente 28% de la production agricole dans la région (MEMEE 2015) et les cultures fourragères, essentiellement le maïs fourrager et la luzerne, occupent une superficie irriguée de 7 500 ha avec une augmentation de 583 ha annuellement (Loulitit 2008).

Matériels et méthodes

Conduite de la reproduction

Le troupeau Holstein de la coopérative COPAG contient environs 80 000 têtes dont 40000 vaches laitières distribuées dans les quatre provinces de la région de Souss-Massa. Ce troupeau est le résultat des génisses nées localement et des génisses introduites dans la région soit d’autres régions du Maroc soit importées de la France, d’Allemagne, du Canada et du pays bas.

Le mode de la reproduction le plus utilisé au niveau des fermes adhérentes était l’insémination artificielle accomplie par une équipe d’inséminateurs de la COPAG. Les semences étaient importées de la France, Canada, l’Allemagne et les États Unis de l’Amérique. En cas de l’échec de la troisième insémination artificielle, la saillie était assurée par la monte naturelle par des taureaux géniteurs, choisis sur la base de leur conformation et la production laitière de leurs mères.

Conduite alimentaire

La ration alimentaire distribuée aux vaches laitières diffère selon leur stade physiologique et le niveau de production. Elle est essentiellement composée de fourrages verts ou conservés en foin ou en ensilage dans des silos produits dans les exploitations. Ainsi les cultures fourragères contiennent essentiellement le maïs fourrager et la luzerne avec des superficies irriguées de 5 000 ha et 2 500 ha, respectivement (Loulitit 2008). Ces surfaces connaissent une augmentation annuelle de 433 ha pour le maïs et de 150 ha pour la luzerne. Les animaux reçoivent aussi des aliments concentrés dont la composition a varie en fonction de l’état de la vache (tourteau de soja, tourteau de colza, tourteau de tournesol, les grains de maïs et les pulpes déshydratées de betterave). Les éleveurs utilisent encore des aliments composés fabriquées par une unité interne de production «COPAG AALAF».

Origine des données et analyse statistique

La base de données est issue des résultats du contrôle laitier officiel effectué au sein de la coopérative agricole COPAG dans la province de Taroudant appartenant à la région agricole Souss-Massa au sud-ouest du Maroc. Les données analysées ont été enregistrées entre 2008 et 2012, regroupant 9 194 lactations appartenant à 6 307 vaches Holstein qui ont été nées entre 2003 et 2010 appartenant à 118 troupeaux adhérents à cette coopérative agricole. Ce cheptel de vaches Holstein présente une moyenne d’âge au vêlage de 44 ± 17 mois s’étalant sur une gamme d’âge comprise entre 19 mois et 115 mois et exposant un coefficient de variation très élevé 39 %.

Le contrôle des performances laitières des vaches Holstein dans la région Souss-Massa se base sur le contrôle laitier par un contrôle mensuel réalisé par un agent agréé sur une période de 24 heures, nommé: protocole A4 recommandée par l’ICAR (International Commitee for Animal Recording, Comité International pour le Contrôle des Performances en Élevage)-«A» pour Agent et « 4 » pour la fréquence- (ICAR 2012). Pour l’interpolation des quantités produites on a utilisé la méthode de Fleischmann recommandée par ce Comité International. Les quantités cumulées ont été standardisées à 305 jours à l'exception des dossiers des vaches qui ont été taries avant les 305 jours.

Pour estimer les effets de l'année de vêlage, la saison de vêlage, le rang de lactation et l’âge au vêlage sur des différents paramètres de production laitière, y compris la quantité de lait en 305 jours (QL_305j), la quantité de matière grasse (QMG_305j), le taux butyreux (TB_305j), on a adopté la procédure suivante:

Pour la saison de vêlage des vaches, quatre saisons ont été considérées: hiver (de décembre à février), printemps (de mars à mai), été (de juin à août), et l'automne (de septembre à novembre).

Concernant l’année de vêlage, les lactations ont été divisées en cinq groupes: 2008, 2009, 2010, 2011 et 2012.

A propos de l’âge au vêlage, 6 groupes ont été créés: groupe 1 (âge ≤ 30 mois), groupe 2 (30 mois <  âge ≤ 42 mois), groupe 3 (42 mois < âge ≤ 54 mois), groupe 4 (54 mois < âge ≤ 66 mois), groupe 5 (66 mois  âge ≤78 mois) et groupe 6 (âge > 78 mois).

Pour le rang de lactation, on a subdivisé les lactations allant de la première jusqu’à la sixième lactation en 6 groupes.

Le modèle linéaire général (GLM) a été utilisé pour les analyses des effets selon le modèle suivant: Y_ijkm = m + A_i + B_j + C_k + E_m + e_ijkm

Y_ijkm est la valeur mesurée, m est la moyenne générale, A_i est l’effet de l’année vêlage, B_j est l’effet la saison de vêlage, C_k est celui relatif à l’âge au vêlage, E_m désigne l’effet du rang de lactation et e_ijkm est l’erreur résiduelle.

Pour les comparaisons multiples des moyennes, on a utilisé le test de comparaison de Newman-Keuls. Les analyses statistiques ont été réalisées en utilisant le logiciel SAS version 9.2 (SAS Institute 2008).

Résultats et discussion

Les moyennes et les écarts-types de la quantité de lait, la quantité de matière grasse et le taux butyreux produits en 305 jours de lactation sont présentés dans le tableau 1. La quantité moyenne de lait en 305 jours (QL_305j) est de 6578 ± 1 619 kg, comprises entre 1 969 et 12 855 kg avec un coefficient de variation 24,6 %. Les moyennes de la QMG_305j et le taux butyreux ont égales respectivement à 241 ± 62 kg et 3,7 ± 0,53 %.

Les performances moyennes enregistrées dans cette étude sont en accord avec celles rapportés par plusieurs études en race Holstein dans des conditions climatiques semblables à celle de la région concernée par cette étude mais restent inférieures à celles obtenues dans les pays tempérés. En Turquie, la QL_305j de la race Holstein varie entre 6 810 ± 56 kg et 6 222 ± 36 kg (Bakir et al 2009a, Kaya et al 2003, Sahin et al 2012, Tekerli et Gündoǧan 2005). En Iran (6 428 ± 1 357 kg) (Nilforooshan et Edriss 2004), en Tunisie on a des valeurs moyennes comprises entre 5 808 ± 78 kg et 5 905 ± 1 895 kg (Ajili et al 2007, Bouallegue et al 2013, M'hamdi et al. 2012), alors qu’au Maroc la QL_305j a été signalée entre 5 353 et 6 293 kg (Boujenane 2002, 2013, Boujenane et Aïssa 2008).

L’analyse de la variance (le tableau 2) a montré que les effets de l'année au vêlage, la saison du vêlage et le rang de lactation sont très hautement significatifs (P < 0,001) pour la quantité de lait en 305 jours et celle de la matière grasse. Pour le taux butyreux, un effet très hautement significatif est observé pour la saison du vêlage et l'année au vêlage. Concernant l’âge de vêlage, l’effet est non-significatif (P > 0,05) pour les différents paramètres de production étudiés.

Les effets de la saison de vêlage

L’effet de la saison de vêlage sur la production laitière est très hautement significatif (P<0.001) pour tous les paramètres de production étudiés (tableau 3), ce qui suggère que la production de lait est très sensible aux variations saisonnières. Les quantités de lait et de matière grasse chez les vêlages de l’automne et de l’hiver sont supérieures à celles obtenues en vêlages de printemps et d’été. Pour les vaches ayant vêlé en automne, on a enregistré la QL_305j significativement la plus élevée (6 866 ± 1 699 kg) suivie par celle des vaches vêlées en hiver (6 576 ± 1 636 kg). De la même manière, les vaches vêlées en automne et en hiver ont produit plus de matière grasse: 251 ± 64,7 kg et 240 ± 63 kg, respectivement.

Dans la province de Taoudant, région Souss-Massa, la saison froide et pluvieuse s’étale de novembre à mars et la saison chaude et sèche d'avril à octobre, sachant que la gamme convenable de température chez les bovins Holstein en allaitement est entre 4°C et 24°C (Hahn 1981). Les effets du stress thermique sur le bétail peuvent être observés au-dessus de 24°C et la production de lait diminue nettement après 27°C (Johnson et Vanjonack 1976). Durant la saison pluvieuse et froide, la température moyenne oscille au voisinage de 19°C alors qu’en saison chaude et sèche d'avril à octobre, dans la région étudiée, les températures ambiantes moyennes sont généralement supérieures à 24°C et peuvent s’accroître jusqu’à 48°C sous l’influence des vents chauds sub-sahariens (Bouchaou et al 2011, Bouragba et al 2011). C’est pourquoi la production lors des saisons de l’été et du printemps est inferieure à celle de l’automne et de l’hiver. De même (Usman et al (2013) ont signalé que les conditions tropicales et subtropicales avec une température ambiante au-dessus de 40°C provoque une chute de la production laitière des Holstein.

Les effets de l’année de vêlage

L’analyse de l’effet de l’année de vêlage montre des différences de production d’une année à l’autre avec un effet très hautement significatif (P <0,001) pour tous les paramètres de production étudiés. La comparaison des moyennes montre que la quantité de lait la plus élevée a été obtenue en 2008 (6 788 ± 1 816 kg) qui se détache significativement des autres années. Les valeurs les plus faibles pour QL_305j et QMG_305j sont enregistrées en 2012 et sont égales respectivement à 6 186 ± 1 430 kg et 221 ± 57,3 kg. Ces résultats sont en concordance avec plusieurs études qui ont rapporté une grande fluctuation de QL_305j en fonction des années et soulignent une influence significative de l’année de vêlage (Bakir et Kaygisiz 2013, Bakir et al 2009b, Çilek 2009, García-Ispierto et al 2007, Katok etYanar 2012, M'hamdi et al 2012, Sahin et al 2012, Tekerli et Gündoǧan 2005).

Les effets du rang de lactation

Le rang de lactation influence significativement (P < 0,001) la QL_305j et la quantité de matière grasse, mais son effet sur le taux butyreux est non significatif (P < 0,05). Les vaches primipares ont produit les plus faibles QL_305j et QMG_305j avec 6 144 ± 1 463 kg et 223 ± 53,6 kg, respectivement. Les quantités de lait et de matière grasse ont connu une nette augmentation jusqu’à un maximum de production en 3^ème lactation avec 7 001 ± 1 635 kg et 257 ± 65,9 kg respectivement puis une diminution progressive avec les lactations subséquentes.

Plusieurs études, ont révélé un effet significatif du rang de lactation sur la production de la quantité de lait (Bakir et al 2009b, Çilek 2009, M'hamdi et al. 2012). Le niveau le plus élevé de production a été obtenu en 4^ème lactation (Topaloǧlu et Günes 2010), en 2^ème lactation (Bakir et al 2009b, Çilek 2009, M'hamdi et al 2012) ou en 3^ème lactation (Bakir et Kaygisiz 2013, Sahin et al 2012).

Conclusion

• Cette étude a montré que les performances laitières des vaches Holstein dans la province de Taroudant, région de Souss-Massa au sud-ouest du Maroc ont été influencées significativement par les facteurs non génétiques analysés. Parmi les facteurs défavorables, le stress thermique créé par les températures élevées pendant le printemps et l’été est vraisemblablement l’un des facteurs les plus limitants pour maintenir une production laitière optimale dans la région. Malgré ces facteurs de variations qui limitent le potentiel productif de cette race, les quantités de lait et de matière grasse sont supérieures à celles enregistrées dans d’autres régions au Maroc. Cela peut être le résultat de l'amélioration des conditions d’élevage (bon encadrement des éleveurs en matière de gestion et de santé, alimentation variée et suffisante, meilleure gestion de la reproduction,…). En revanche, l’amélioration du rendement laitier doit optimiser les effets environnementaux limitants en prenant en considération ces facteurs lors de la gestion des élevages. Les éleveurs peuvent par exemple ajuster la saison de vêlage de leurs troupeaux avec les saisons pluvieuses et froides (l’automne surtout et l’hiver). En fin, la grande variabilité des paramètres productifs de la race Holstein dans la région nécessite l’évaluation d’autres facteurs de variation pour une meilleure rentabilité de ces élevages ; tels l’effet troupeau, l’âge au premier vêlage, l’origine des vaches.

Remerciements

Les auteurs remercient Mr. Zaid Annasser directeur de la production animale, les ingénieurs agricoles, les techniciens et les contrôleurs laitiers de la coopérative agricole COPAG qui ont mis les données analysées à leur disposition.

Références bibliographiques

ABHSM 2015 Climat du Souss-Massa. Agence du Bassin Hydraulique du Souss-Massa (ABHSM), Agadir, Maroc, from http://www.abhsmd.ma/spip.php?article63

Ajili N, Rekik B, Gara A B and Bouraoui R 2007 Relationships among milk production, reproductive traits, and herd life for Tunisian Holstein-Friesian cows. African Journal of Agricultural Research, volume 2(2), article E9A007528937, from http://www.academicjournals.org/journal/AJAR/article-abstract/E9A007528937

Albarrán-Portillo B and Pollott G E 2011 Environmental factors affecting lactation curve parameters in the United Kingdom’s commercial dairy herds. Archivos de Medicina Veterinaria, volume 43(2), from http://www.veterinaria.uach.cl/archivos_med_vet/en/papers/?v=138

Bakir G and Kaygisiz A 2013 Milk Yield Characteristics of Holstein Cows and the Effect of Calving Month on Milk Yield. Journal of Natural Sciences, volume 16 (1), article #1017000107, from http://dergi.ksu.edu.tr/article/view/1017000062/1017000107

Bakir G, Kaygisiz A and Cilek S 2009a Estimates of genetic trends, for 305-Days milk yield in Holstein Friesian cattle. Journal of Animal and Veterinary Advances, volume 8 (12), article 2553.2556, from http://medwellpublishing.com/fulltext/?doi=javaa.2009.2553.2556

Bakir G, Kaygisiz A and Cilek S 2009b Milk yield traits of Holstein cattle reared at Tahirova State Farm in Balikesir province in Turkey. Journal of Animal and Veterinary Advances, volume 8 (11), article 2369.2374, from http://medwelljournals.com/fulltext/?doi=javaa.2009.2369.2374

Bouallegue M, Haddad B, Aschi M and Hamouda M B 2013 Effect of environmental factors on lactation curves of milk production traits in Holstein-Friesian cows reared under North African conditions. Livestock Research for Rural Development, Volume 25, Article #75, from http://www.lrrd.org/lrrd25/5/boua25075.htm

Bouchaou L, Tagma T, Boutaleb S, Hssaisoune M and El Morjani Z E A 2011 Climate change and its impacts on groundwater resources in Morocco: the case of the Souss-Massa basin, Chapter 8:129-144. In Treidel H, Martin-Bordes J L and Gurdak J J (eds.). Climate Change Effects on Groundwater Resources: A Global Synthesis of Findings and Recommendations, Volume 27, 414p, CRC Press/Balkema, Leiden, Netherlands.

Bouchaou L, Michelot J, Vengosh A, Hsissou Y, Qurtobi M, Gaye C, Bullen T and Zuppi G 2008 Application of multiple isotopic and geochemical tracers for investigation of recharge, salinization, and residence time of water in the Souss-Massa aquifer, southwest of Morocco. Journal of Hydrology, Volume 352 (3-4): 267-287.

Boujenane I 2002 Estimates of Genetic and Phenotypic Parameters for Milk Production in Moroccan Holstein-Friesian Cows. Revue d'Elevage et de Médecine Vétérinaire des Pays Tropicaux, Volume 55 (1), Article 511853 from http://remvt.cirad.fr/revue/notice_fr.php?dk=511853

Boujenane I 2013 Comparison of Different Lactation Curve Models to Describe Lactation Curve in Moroccan Holstein-Friesian Dairy Cows. Iranian Journal of Applied Animal Science, Volume 3(4), Article 902, from http://ijas.ir/main/modules/content/index.php?id=902

Boujenane I et Aïssa H 2008 Performances de reproduction et de production laitière des vaches de race Holstein et Montbéliarde au Maroc. Revue d'Elevage et de Médecine Vétérinaire des Pays Tropicaux, Volume 6 (3-4), Article 556170, from http://remvt.cirad.fr/revue/notice_fr.php?dk=556170

Bouragba L, Mudry J, Bouchaou L, Hsissou Y, Krimissa M, Tagma T and Michelot J L 2011 Isotopes and groundwater management strategies under semi-arid area: Case of the Souss upstream basin (Morocco). Applied Radiation and Isotopes, Volume 69: 1084-1093.

Cienfuegos-Rivas E, Blake R, Oltenacu P and Castillo-Juarez H 2006 Fertility responses of Mexican Holstein cows to US sire selection. Journal of dairy science, Volume 89(7), from http://www.journalofdairyscience.org/article/S0022-0302%2806%2972352-9/fulltext

Çilek S 2009 Milk yield traits of Holstein cows raised at polatli state farm in Turkey. Journal of Animal and Veterinary Advances, Volume 8 (1), article 6.10, from http://medwelljournals.com/fulltext/?doi=javaa.2009.6.10

García-Ispierto I, López-Gatius F, Bech-Sabat G, Santolaria P, Yániz J, Nogareda C, De Rensis F and López-Béjar M 2007 Climate factors affecting conception rate of high producing dairy cows in northeastern Spain. Theriogenology, volume 67(8): 1379-1385.

González-Recio O, Pérez-Cabal M and Alenda R 2004 Economic value of female fertility and its relationship with profit in Spanish dairy cattle. Journal of dairy science, Volume 87 (9), from http://www.journalofdairyscience.org/article/S0022-0302%2804%2973438-4/fulltext

Hahn G L 1981 Housing and management to reduce climatic impacts on livestock. Journal of Animal Science, volume 52(1), from https://www.animalsciencepublications.org/publications/jas/abstracts/52/1/JAN0520010175?search-result=1

Hammami H, Rekik B and Gengler N 2009 Genotype by environment interaction in dairy cattle. Biotechnologie Agronomie Société et Environnement, volume 13(1), from http://www.pressesagro.be/base/index.php/base/article/view/387

Heins B, Hansen L and De Vries A 2012 Survival, lifetime production, and profitability of Normande × Holstein, Montbéliarde × Holstein, and Scandinavian Red × Holstein crossbreds versus pure Holsteins. Journal of dairy science, volume 95(2), from http://www.journalofdairyscience.org/article/S0022-0302%2812%2900052-5/fulltext

ICAR 2012 International agreement of recording practices. Guidelines approved by the General Assembly held in Cork, Ireland, June 2012. International Committee for Animal Recording, Rome, Italy, from http://www.icar.org/Documents/Rules%20and%20regulations/Guidelines/Guidelines_2012.pdf

Javed K, Mohiuddin G and Akhtar P 2001 Heritability estimates of some productive traits in Sahiwal cattle. Pakistan Veterinary Journal, volume 21, from http://www.pvj.com.pk/abstract/21_3/1.htm

Johnson H and Vanjonack W 1976 Effects of environmental and other stressors on blood hormone patterns in lactating animals. Journal of dairy science, volume 59(9), from http://www.journalofdairyscience.org/article/S0022-0302%2876%2984413-X/abstract

Katok N and Yanar M 2012 Milk traits and estimation of genetic, phenotypic and environmental trends for milk and milk fat yields in holstein friesian cows. International Journal of Agriculture & Biology, Volume 14(2): 311-314.

Kaya I, Uzmay C, Kaya A and Akbas Y 2003 Comparative analysis of milk yield and reproductive traits of Holstein-Friesian cows born in Turkey or imported from Italy and kept on farms under the Turkish-ANAFI project. Italian Journal of Animal Science, volume 2, from http://www.aspajournal.it/index.php/ijas/article/view/ijas.2003.141

Knaus W 2009 Dairy cows trapped between performance demands and adaptability. Journal of the Science of Food and Agriculture, volume 89 (7): 1107-111.

La Bonnardière C, Bonaiti B, Abrial D, Gasqui P, Calavas D, Ducrot C and Barnouin J 2007 Milk yield, age at first calving, and the risk of BSE: an analysis at the farm level in France. Preventive veterinary medicine, volume 78(1): 67-78.

Loulitit M M 2008 Action de développement de la COPAG. Revue Homme Terre et Eau, volume 139, article 5, from http://www.anafide.org/Revu%202008.html

M'hamdi N, Bouallegue M, Hamouda MB, Frouja S, Brar SK and Ressaissi Y 2012 Effects of Environmental Factors on Milk Yield, Lactation Length and Dry Period in Tunisian Holstein Cows. Chapter 8: 153-164, In Narongsak C (ed.). Milk Production-An Up-to-Date Overview of Animal Nutrition, Management and Health", 404p, INTECH Open Access Publisher.

Manssouri T, Sahbi H and Manssouri I 2014 Elaboration of stochastic mathematical models for the prediction of parameters indicative of groundwater quality Case of Souss-Massa Morocco. IJCER International Journal of Computational Engineering Research 4(5),article III-3, from http://www.ijceronline.com/v4i5.html

MEMEE 2015 Ressources en eau en Souss-Massa. Ministère de l'Énergie des Mines de l'Eau et de l'Environnement, Rabat, Maroc, from http://www.water.gov.ma/index.cfm?gen=true&id=13&ID_PAGE=32

Nilforooshan M and Edriss M 2004 Effect of age at first calving on some productive and longevity traits in Iranian Holsteins of the Isfahan Province. Journal of dairy science, volume 87(7), from http://www.journalofdairyscience.org/article/S0022-0302%2804%2970032-6/abstract

Nyamushamba G B, Halimani TE, Imbayarwo-Chikosi V E and Tavirimirwa B 2014 Comparative evaluation of non-genetic factors affecting milk yield and composition of Red Dane and Jersey cattle in Zimbabwe. SpringerPlus, volume 3, from http://www.springerplus.com/content/3/1/88

Nyamushamba G B, Dhliwayo M, Chimonyo M, Tavirimirwa B and Banana N Y D 2013 Non-genetic factors affecting milk yield and composition of Holstein-Friesian cows nested within natural ecological regions of Zimbabwe. Scientific Journal of Animal Science, volume 2(5), Article 102-108, from http://www.sjournals.com/index.php/SJAs/issue/view/113

Pirlo G, Miglior F and Speroni M 2000 Effect of age at first calving on production traits and on difference between milk yield returns and rearing costs in Italian Holsteins. Journal of dairy science, volume 83(3), from http://www.journalofdairyscience.org/article/S0022-0302%2800%2974919-8/abstract

Prim'Holstein France 2015 La Prim'Holstein, from http://primholstein.com/la-prim-holstein/presentation-primholstein/

Ravagnolo O and Misztal I 2002 Effect of heat stress on nonreturn rate in Holsteins: fixed-model analyses. Journal of dairy science, volume 85(11), from http://www.journalofdairyscience.org/article/S0022-0302%2802%2974397-X/fulltext

Riecka Z and Candrák J 2011 Analysis of relationship between production and reproduction traits of Holstein cattle population in the Slovak Republic. Scientific Papers, Animal Science and Biotechnologies, volume 44: 332-336.

Sahin A, Ulutas Z, Adkinson A Y and Adkinson R W 2012 Genetic and environmental parameters and trends for milk production of Holstein cattle in Turkey. Italian Journal of Animal Science, volume 11(3), Article e44, from http://www.aspajournal.it/index.php/ijas/article/view/ijas.2012.e44

SAS Institute Inc 2008 Statistical analysis Software version 9.2, Cary, NC: SAS Institute Inc. USA

Scholtz M, Furstenburg D, Maiwashe A, Makgahlela M, Theron H and Van der Westhuizen J 2010 Environmental-genotype responses in livestock to global warming: A Southern African perspective. South African Journal of Animal Science, Volume 40 (5), from http://www.ajol.info/index.php/sajas/article/view/95890

Tagma T, Hsissou Y, Bouchaou L, Bouragba L and Boutaleb S 2009 Groundwater nitrate pollution in Souss-Massa basin (South-West Morocco). African Journal of Environmental Science and Technology, Volume 3(10), from http://www.ajol.info/index.php/ajest/article/view/56258

Tekerli M and Gündoǧan M 2005 Effect of certain factors on productive and reproductive efficiency traits and phenotypic relationships among these traits and repeatabilities in West Anatolian Holsteins. Turkish Journal of Veterinary and Animal Sciences,Volume 29(1), from http://mistug.tubitak.gov.tr/bdyim/abs.php?dergi=vet&rak=0206-19

Topaloǧlu N and Günes H 2010 Effects of Some Factors on Milk Yield and Components of Holstein-Friesian Cattle in England. Journal of the Faculty of Veterinary Medicine Istanbul, volume 36(1), Article 07, from http://vetjournal.istanbul.edu.tr/--2010-1-en.html

Usman T, Qureshi M S, Yu Y and Wang Y 2013 Influence of various environmental factors on dairy production and adaptability of Holstein cattle maintained under tropical and subtropical conditions. Advances in Environmental Biology, Volume 7(2), Article 366-372, from http://www.aensiweb.com/old/aeb_Feb-_2013.html

Wicks H C F and Leaver J D 2006 Influence of genetic merit and environment on somatic cell counts of Holstein-Friesian cows. The Veterinary Journal, Volume 172(1): 52-57.

Received 10 April 2015; Accepted 26 April 2015; Published 3 June 2015

Go to top