Introduction

1On a vu comment l'aridité pouvait être une contrainte au développement, au Nord-Mexique comme dans d'autres régions. La connaissance de la répartition des précipitations est rendue difficile dans ces régions peu peuplées et de fait équipées d'un réseau de mesure d'autant plus lâche que de nombreuses stations ont été arrêtées ces dernières années.

2La répartition des précipitations annuelles sur la partie centrale de la Sierra Madre occidentale (correspondant aux régions hydrologiques n° 36 et 10) telle qu'on peut l'apprécier sur la figure 21 montre :

• une nette opposition entre les versants ouest et est de la sierra, le versant ouest étant bien plus arrosé mais présentant une bien plus forte variabilité spatiale, fruit d'un relief extrêmement escarpé (Nouvelot et Descroix, 1996 ; Descroix et al., 1997 ; Descroix et al., 2001) ; n régulière des précipitations de la ligne de crête (où elles sont de l'ordre de 800 à 1 000 mm) vers la Laguna et l'est du bassin Nazas Aguanaval (où elles sont localement inférieures à 200 mm par an) ;
• une diagonale de très fortes précipitations, qui se sépare du littoral pour se rapprocher de la ligne de crête, du sud vers le nord ; on atteint des valeurs de plus de 1 500 mm par an à moins de 80 km du littoral au sud, alors que les maxima, proches de 1 200 mm, sont à plus de 150 km de la côte au nord ;
• l'apparition d'une zone littorale aride au nord, où l'on se rapproche des déserts côtiers du Sonora et de Basse-Californie (moins de 300 mm à Topolobampo).

3La figure 22 donne les valeurs du coefficient de variation des précipitations interannuelles moyennes pour le bassin Nazas Aguanaval et une grande partie du versant ouest de la sierra. Plus ce coefficient est élevé, plus l'irrégularité interannuelle de la précipitation est forte. On constate que l'irrégularité croît très nettement de la côte vers l'intérieur des terres. Ici, les zones où le Cv est le plus faible correspond aux zones les plus pluvieuses du versant ouest de la sierra ; les plus forts Cv se rencontrent dans les zones arides du sud du désert de Chihuahua, c'est-à-dire dans les dépressions endoréiques du Bolsón de Mapimí.

La spatialisation des précipitations

4Dans les secteurs peu instrumentés, comme l'est en général le nord du Mexique, il apparaît un jour ou l'autre un besoin de connaître plus ou moins précisément les hauteurs de précipitation en vue de dimensionner des ouvrages (barrages, canaux, ponts, etc.) ou de faire des prévisions de pluies extrêmes et donc de risques naturels (crues, inondations, etc.). Aussi, comme le réseau est lâche, on cherchera à déterminer deux types d'informations complémentaires à celles fournies par le réseau, et qui permettent d'interpoler spatialement la validité de l'information acquise, afin de pallier l'existence de zones sans données :

• la première est la définition de zones homogènes en fonction de l'information pluviométrique recueillie sur le réseau existant ;
• la seconde est la distance de corrélation des données acquises.

5C'est pour parvenir à la définition de régions homogènes que l'on a fait appel une nouvelle fois (cf. « Le climat et l'aléa pluviométrique au Nord-Mexique », p. 129) à la Méthode du vecteur régional (MVR). Cette technique, proposée par Hiez (1986), est basée sur la notion de pseudoproportionnalité des totaux pluviométriques entre stations voisines dont l'ensemble des données permet ainsi l'élaboration d'un vecteur représentatif de la région couverte. Cette élaboration s'appuie sur un principe de maximum de vraisemblance qui part du postulat que l'information la plus vraisemblable est celle qui se répète le plus fréquemment.

6Mais alors que précédemment, on a utilisé cette méthode pour vérifier la validité des données et repérer d'éventuelles erreurs dans les séries, on s'en sert ici pour la détermination par défaut de zones « pluviométriquement » homogènes. C'est-à-dire que l'on considère que si un vecteur est constitué par un ensemble de stations au comportement « proportionnellement » semblable, il définit une région homogène en terme de comportement pluviométrique.

7L'analyse de toutes les données disponibles sur le réseau existant a permis de constituer sept « régions pluviométriques », cinq pour le bassin du Nazas Aguanaval (RH36) et deux pour le versant Pacifique de la sierra (RH10) (Gutierrez, 2003), constitué de nombreux bassins versants dont les eaux s'écoulent directement vers le Pacifique (fig. 23).

Le rôle des facteurs locaux et zonaux dans la distribution spatiale des précipitations

8Il est toujours délicat de rechercher les facteurs locaux pouvant expliquer la répartition spatiale des précipitations. Néanmoins, le premier regard sur la carte des isohyètes (fig. 21) permet de comprendre le rôle primordial de l'altitude et de la distance à l'Océan. Ceci est d'ailleurs paradoxalement bien plus net pour le versant intérieur (oriental) de la Sierra Madre occidentale. On y voit nettement la pluviométrie diminuer en même temps que l'altitude (on passe de 900 mm vers 2 800 m à 200 mm vers 1 100 m), et inversement à la distance au Pacifique. Les deux évolutions sont logiques et se retrouvent sous pratiquement tous les climats. Simplement, on a ici, du fait que la sierra est parallèle au trait de côte (ce qui est dû à la tectonique régionale), une addition des deux phénomènes sur le versant interne.

9Sur le versant Pacifique de la sierra, l'évolution est forcément plus complexe pour la simple raison que l'altitude augmente avec l'éloignement de la côte, ce qui fait que les deux facteurs se contrarient. Mais la complexité est accentuée par :

• le relief, très escarpé, du versant ouest de la chaîne ;
• un fort gradient latitudinal lié aux climats zonaux.

10Le relief impose une très forte variabilité, spatiale et temporelle, des précipitations ; celle-ci serait encore plus marquée si la densité des postes était plus forte, car la variabilité est ici à l'échelle de ces vallées profondes et des lambeaux de plateaux qui les dominent parfois de plus de 1 500 m.

11Par ailleurs, on est ici à des latitudes subtropicales ; on passe de fait, de climats tropicaux, au sud à des climats subtropicaux, au nord, caractérisés par un désert côtier sur le littoral (désert du Sonora) et des secteurs arides continentaux sur l'altiplano nord-central mexicain (désert de Chihuahua). De ce fait, il y a une configuration régionale de la pluviométrie, les secteurs littoraux ou leur immédiat arrière-pays étant les plus arrosés en zone tropicale. Mais 500 km plus au nord, les littoraux sont concernés zonalement par l'aridité côtière des façades ouest de continent aux latitudes subtropicales.

12Du fait de cette complexité, on s'est intéressé ici uniquement au versant interne de la sierra ; on a voulu déterminer les facteurs locaux de la pluviométrie. Pour ce faire, on a réalisé deux analyses statistiques : une Analyse en composantes principales (ACP) et une Analyse factorielle des correspondances (AFC). Ces analyses permettent de classer et hiérarchiser les variables explicatives d'un phénomène (qui joue ici le rôle de variable dépendante vis-à-vis des autres, considérées comme indépendantes). L'ACP permet aussi d'en supprimer quelques-unes en mettant en évidence la redondance des variables. Ces analyses ont été faites à partir des données de précipitation annuelle des 52 stations les plus documentées de la Région hydrologique 36 (bassin du Nazas-Aguanaval).

13L'ACP a été menée avec comme variable dépendante la pluviométrie, comme variables indépendantes l'altitude, la distance à l'océan Pacifique, la latitude, la longitude et la densité de la végétation sur trois aires (0,3, 16 et 94 km²) autour du poste pluviométrique. La figure 24 montre l'espace des variables basé sur les deux premières composantes de l'ACP, qui représentent 75 % de la variance expliquée.

14On remarque que la pluviométrie est le plus corrélée avec l'altitude (positivement) et avec la distance à l'Océan (MER) négativement. La végétation et la longitude (ouest, donc valeurs croissantes vers l'ouest) sont également corrélées positivement, mais à un moindre degré. Ceci, de même que le fait que ces dernières variables sont bien corrélées entre elles, s'expliquent aisément parce que les régions les plus hautes et les plus pluvieuses correspondent aux secteurs les plus proches de la ligne de partage des eaux avec les cours d'eau se jetant dans le Pacifique, c'est-à-dire les plus occidentales.

15Pour pouvoir prendre en compte des valeurs non numérisables, on a mené à bien une AFC, analyse statistique qui est le pendant de l'ACP mais permet d'introduire des valeurs non numériques. Il faut toutefois déterminer des classes pour chacune des variables prises en compte, afin de les hiérarchiser si nécessaire, et de toutes manières, de les distinguer les unes des autres. Il faut noter qu'en établissant des classes aussi pour les variables qui ont des valeurs numériques (déjà prises en compte par l'ACP), on peut aussi les analyser dans l'AFC.

16Les variables étudiées ici sont rassemblées dans le tableau XII. Ce sont celles qui ont servi à réaliser l'ACP, exceptées la longitude et la latitude, mais on y a ajouté le type de site et deux types d'exposition, celle du site (échelle kilométrique) et celle, plus globale, de la région du site (échelle de la vingtaine de km).

17Les résultats sont exposés figure 25 ; celle-ci représente le plan des variables telles que classées (d'après le tableau XII), sur les axes 1 et 2 de l'AFC, qui représentent 60 % de la variance expliquée des variables.

18On peut constater que l'AFC permet de distinguer quatre « régions pluviométriques » différentes :

• une région de pluviométrie très faible (inférieure à 300 mm) ; elle est associée à un éloignement maximal de l'Océan, une altitude et des densités de végétation minimales, et à une situation en plaine ou de fond de vallée encaissée ; cette dernière configuration s'explique par la situation d'extrême abri de certains postes situés pourtant dans la Sierra Madre, mais dans des sites très protégés ;
• une région de faible pluviométrie (300 à 400 mm) associée à une distance à l'océan de 300 à 350 km, un indice de végétation normalisée NDVI (Normalized Difference Vegetation Index) faible à moyen, une position de piedmont ou de col, une altitude comprise entre 1 200 et 1 600 m ;
• une région de pluviométrie moyenne, où l'on a une situation de versant ou de fond de vallée majoritaire, une pluviométrie de 400 à 550 mm, une altitude supérieure à 1 600 m, une distance à la mer comprise entre 250 et 300 km, un indice de végétation moyen à élevé ; cela correspond aux bas versants et aux fonds des deux grandes vallées du haut Nazas, dans la Sierra Madre (Sextin et Ramos) ;
• et une région de pluviométrie élevée (+ de 550 mm), qui correspond aux plus forts indices de végétation, aux plus hautes altitudes, aux zones les plus proches de l'océan Pacifique. Ces stations sont celles situées sur les hauts versants ou sur la ligne de partage des eaux entre les deux versants (RH36 et RH10).

19On constate par ailleurs que hormis la première zone où aucune exposition n'est remarquée (on est au fond du Bolsón de Mapimi essentiellement), les types d'exposition ne permettent pas de classer les régions ; on les trouve, tous, dans les trois autres zones.

20La régionalisation telle que définie par l'analyse statistique n'est pas très différente de celle obtenue grâce à la Méthode du vecteur régional. Il lui manque les régions 6 et 7 du découpage MVR, qui sont absentes ici car on n'a pas considéré le versant Pacifique de la sierra. Par ailleurs, les régions 2 et 3 sont ici découpées presque de la même manière car les plateaux de l'Aguanaval (région 3 de la figure 23) sont à plus haute altitude que la vallée du Nazas (région 2), ce qui lui permet d'en être distinguée par l'AFC qui tient compte de l'altitude ; mais ces régions sont néanmoins grosso modo les mêmes.

Conclusion

21Les régions homogènes restent les mêmes quelle que soit la méthode utilisée, montrant une certaine redondance de l'information (Descroix ef al., 2001 ). Les régions telles que définies, ainsi que les variogrammes destinés à connaître la distance de dé-corrélation (ou, à l'inverse, la distance maximale de validité du signal pluviométrique), semblent donc suffisamment robustes pour permettre une extension de l'information obtenue sur les quelques postes pour lesquels on a de longues séries d'observations, et qui fournissent pour l'essentiel une information de bonne qualité.

22L'une ou l'autre et a fortiori l'ensemble de ces méthodes permettent de pallier la faible densité du réseau et d'obtenir une information spatialisée pour mieux prévoir les quantités précipitées sur les zones peu instrumentées, et régionaliser aussi les risques de crues et d'inondation et à l'inverse, de sécheresses et de pénurie.