Constitution d`un système d`information à référence spatiale
Contents
1. CAPRA H 2002 Simple predictions of instream habitat model outputs for target fish populations Freshwater Biology 47 8 p 1543 1556 PEREZ CORREA MACARENA 2004 D veloppement d une m thode de cartographie de l occupation du sol le long des cours d eau partir de donn es de t l d tection rapport de mast re SILAT ROGERS C F WASSON J G 1997 Cadrage g ographique de la structure et de la dynamique du physique des hydrosyst mes du bassin Rhodanien Fran ais L approche hydro cor gions rapport GIP hydrosyst mes Zone Atelier bassin du Rh ne 03 97 Universit de Saint Etienne CRENAM ENRS UMR 5600 Cemagref Lyon BEA LHQ SOUCHON Y ANDRIAMAHEFA H COHEN P BREIL P PELLA H LAMOUROUX N MALAVOI J R WASSON J G 2000 R gionalisation de l habitat aquatique dans le bassin de la Loire Agence de l eau Loire Bretagne Cemagref Bely LHQ in Ecological method 1978 Southwood TRE London Chapman and Hall p 291 19 N 2 2 w D l o n w S w bn 20 WASSON J G CHANDESRIS A PELLA H 2002 D finition des hydro cor gions de France m tropolitaine Approche r gionale de typologie des eaux courantes et l ments pour la d finition des peuplements de r f rence d invert br s rapport Cemagref Lyon BEA LHQ et MATE DE 190 p WASSON J G CHANDESRIS A PELLA H BLANC L 2004 Typologie de2. partir des informations de relief utilis es ci dessus CO Hvdro cen giens nireau 1 infiuence des roches sur A chimie de leu Eki d ikide oranga ei Eisi pii cotormisss i BON Cas cartoni M dise Nga pris G ologie les cartes existantes l chelle de la France m tropolitaine constituent avant tout une description stratigraphique et structurale du sous sol Elles ne renseignent pas directement la nature des roches et ne sont donc pas imm dia tement exploitables pour une spatialisation des caract ristiques des bassins versants Nature des roches lithologie plusieurs sources de donn es compl mentaires ont t utilis es pour rendre compte de l influence de la litho logie sur le fonctionnement des cosyst mes aquatiques La carte g ologique Bonnefoy et al 1996 du BRGM au 1 1 000 000 sixi me dition a fait l objet d une analyse et d un classement des polygones en fonction de la nature des roches identifi es partir de docu ments cartographiques plus grande chelle et des synth ses nationales du BRGM La classification simple labor e sp cifiquement pour cette tude a pour objectif de d crire l in fluence de la g ologie sur le fonctionnement Constitution d un syst me d information r f rence spatiale dans le contexte de la Directive cadre Herv Pella et al des cours d eau r sistance l rosion perm abi lit propri t s chimiques des
3. transfert de cet outil est envisageable sous r serve d une gestion des copyright appropri e D termination des hydro cor gions Principes m thodologiques La d marche est bas e sur une approche des cendante dont le principe de base consiste utiliser pour la classification les facteurs de con tr le globaux plut t que leurs cons quences au niveau local l chelle du bassin versant les d terminants primaires universellement reconnus du fonctionnement cologique des cours d eau sont la g ologie le relief et le climat Wasson et al 2004 Wasson et al 2002 Des cou ches d information spatialement homog nes portant sur les structures g ophysiques et cli matiques sont utilis es au moyen d un SIG Les HER sont d limit es visuellement en recher chant les discontinuit s des diff rents facteurs Les donn es g ographiques utilis es Les donn es utilis es visent d crire spatiale ment deux grands ensembles de facteurs d une part la structure physique des bassins et d autre part le climat qui s exerce sur ces bassins en tant que d terminant aussi bien de la quantit et de la variabilit des coulements que de la temp rature e La structure physique des bassins est appr hend e partir de la g ologie qui intervient au travers des caract ristiques lithologiques r sistance l ro sion perm abilit propri t s chimiques des eaux qui en d coulent du r
4. laboration de cette plate forme SIG nous a demand un travail important de collecte de validation de mise en coh rence et d archivage des donn es Cette base de donn es g ographiques est com pos e d informations brutes pluviom trie altitude et d informations d riv es den sit de drainage pente calcul es par analyse spatiale Ainsi par exemple le mod le num rique de terrain MNT au pas de 250 m de l Institut g ographique national IGN au format raster est l information de base pour le calcul des cou ches suivantes les pentes le mod le num rique 1 Institut fran ais de l environnement Cemagref UR Biologie des cosyst mes aquatiques 3 bis quai Chauveau CP 220 69366 Lyon Cedex 09 11 N 2 2 w So O l o n w pm S O bn 12 2 Bureau des recher ches g ologiques et mini res 3 Centre d tudes sur les r seaux les trans ports l urbanisme et les constructions publiques Climat M t o France Banque de donn es de la Carte g ologique au 1 1 000 000 6 dition BRGM 1996 Mod le num rique de terrain MNT au pas de 250 m d riv de la BDALTI de l IGN Carte des moyennes mensuelles des pr cipita tions pour la p riode 1961 1990 calcul es par Convention IFEN Convention IFEN Convention IFEN la m thode AURHELY maille 5 x 5 km M t o France Carte des moyennes mensuelles de
5. Ces derni res se d finissent par un ensemble structur d informations d crivant une ressource L objectif principal de ce catalogage est de pouvoir consulter facilement toutes les informations relatives une couche comme par exemple l origine l chelle la pr cision du trac etc Apr s quelques recherches sur les outils existants permettant de formaliser stocker et interroger ces m tadonn es nous avons choisi l outil d velopp par le p le g omatique du CERTU Conforme la norme PR ENV 12657 relative aux m tadonn es des informations localis es le logiciel Reports V2 permet la cons titution et la maintenance de catalogues de don n es g ographiques Celui ci est dot de moyens d importation et d exportation de m tadonn es Les formats de celles ci sont publics afin de permettre d autres logiciels d changer selon ce format Le catalogue pr cise pour chacune des couches les informations relatives leur production leur mode d acquisition leur qualit et enfin leur r f rencement La partie production renseigne sur l organisme qui est l origine de la donn e son r le la personne contact au sein de l orga nisme ainsi que sur l utilisation potentielle qui peut en tre faite et les documents de r f rence La propri t intellectuelle et les restrictions d utilisation sont aussi pr cis es La partie mode d acquisition d taille le cadre dans lequel le laboratoire d tie
6. Ing nieries N 40 p 11 20 D cembre 2004 Constitution d un syst me d information r f rence spatiale dans le contexte de la Directive cadre europ enne sur l eau Herv Pella Andr Chandesris et Jean Gabriel Wasson partir de l exploitation de donn es biologiques repr sentatives de la qualit cologique des cours d eau et de donn es sur l occupation des sols les auteurs expliquent dans cet article comment il est possible de d velopper un syst me d information g ographique destin la d termination typologique des cours d eau L objectif de cette d termination s inscrit dans le projet de r gionalisation des cosyst mes aquatiques en hydro cor gions pour r pondre l objectif principal de bonne qualit cologique des cours d eau de la Directive cadre europ enne sur l eau a Directive cadre europ enne sur l eau DCE pose clairement le principe d une valuation de l tat cologique des milieux aquatiques par la mesure d un cart des conditions de r f rence principale ment sur la base de bioindicateurs La premi re tape de sa mise en uvre en ce qui concerne les eaux de surface continentales passe donc par l tablissement d une typologie permet tant la d finition des conditions de r f rence biologiques mais aussi hydromorphologiques et physicochimiques L article intitul Typologie des eaux courantes pour la Directive
7. atiques travers le filtre de la v g tation terrestre Cette carte a t utilis e comme base de travail et a permis la validation de certains param tres avec toutefois une r interpr tation approfondie pour une r gio nalisation des milieux aquatiques Analyse des pressions Cette analyse consiste valuer le degr d an thropisation des cours d eau de fa on coh rente et homog ne l chelle du territoire national dans un double objectif d une part mettre en vidence les zones les moins anthropis es pour s lectionner des sites de r f rence et d autre part laborer des mod les pressions impacts reliant l alt ration des peuplements aquatiques aux activit s anthropiques rejets polluants agriculture occupation du sol Le couplage entre les donn es biologiques observ es sur des stations ponctuelles et les pressions s exer ant sur les bassins et sur l environnement imm diat des stations corridor rivulaire a n cessit des d veloppements informatiques sp cifiques lan gage objet sous SIG Donn es utilis es Donn es biologiques l indice biologique global normalis IBGN bas sur les invert br s ben thiques est la m thode la plus largement utilis e en France pour valuer la qualit cologique des rivi res Une base de donn es sur les stations IBGN suivies par les DIREN est en cours de d veloppement au sein du laboratoire Actuelle ment plus de 3 500 statio
8. cadre europ enne sur l eau l approche par hydro cor gions Wasson etal 2004 revient sur les objectifs de la typologie expose le cadre conceptuel des hydro cor gions HER ainsi que les principes m thodologiques qui ont conduit leur d limitation et traite la question de leur validation biologique Le pr sent article est compl mentaire de celui cit pr c demment puisqu il aborde l aspect g omatique de cette tude Nous pr sentons tout d abord la mise en place du syst me d in formation g ographique SIG en r ponse aux objectifs de la DCE Ensuite nous voquons le traitement et l analyse des informations g o graphiques pour la d limitation des HER Enfin un exemple d utilisation des donn es pour carac t riser les pressions anthropiques est expos Constitution du SIG Donn es g ographiques l origine du SIG Les donn es utilis es dans le cadre de cette tude sont essentiellement des couches d informations g ographiques acquises par le Cemagref ou mises disposition pour les besoins de l tude par l IFEN dans le cadre de la convention entre le Cemagref et le minist re de l cologie et du D veloppement durable Celles ci sont d taill es dans le tableau 1 page 12 Mise en place du SIG Un SIG l chelle nationale a t cr au sein du laboratoire pour r pondre aux objectifs de r gio nalisation et pour caract riser le milieu naturel et les pressions qui s y exercent L
9. de traitement BASSIN VERSANT DE LA STATION Compte tenu du nombre de stations ponctuelles sur lesquelles nous avons des donn es biologi ques observ es nous avons d velopp un outil sp cifique de calcul de bassin versant Cet outil est bas sur l analyse du mod le num rique de terrain d riv de la BDALTI de l IGN au pas de 250 m sur le d coupage en zones hydrographi ques et enfin sur le trac du r seau hydrographi que de la BD Carthage V 3 0 Une fois la station raccroch e g ographique ment sur le trac du cours d eau auquel elle appartient le traitement consiste d couper le mod le num rique de terrain en fonction de la zone amont probable partir des zones hydrographiques puis d terminer son bassin versant figure 3 Constitution d un syst me d information r f rence spatiale dans le contexte de la Directive cadre Herv Pella et al tant donn que ce calcul de bassin versant est bas sur le MNT de l IGN au pas de 250 m la surface minimale de d tection est de quelques kilom tres carr s En zone de relief peu marqu cette d marche est impossible CORRIDOR RIVULAIRE DE LA STATION Les hypoth ses de travail sont les suivantes le tron on tudi mesure environ 3 km de long soit 1 5 km l amont et l aval de chaque station consid r e la largeur de la zone prendre en compte de part et d autre du cours d eau d pend tout natu rellement de la la
10. eaux figure 1 Cartes g omorphologiques ces cartes au 1 1 000 000 en quatre parties dit es par le GIP Reclus sous la direction de Fernand JOLY individualisent des unit s diff renci es par leur style topographique et secondairement par leur altitude ou leur mat riel Elles constituent un compl ment tr s utile et tr s efficace pour corriger les imperfections inh rentes la carte lithologi que labor e dans le cadre de ce travail Densit du r seau hydrographique la BD Carthage acronyme de Base de Donn es de Cartographie Th matique des Agences de l Eau constitue le r f rentiel spatial commun l ensemble des organismes publics en charge de missions li es la gestion de l eau Elle a t constitu e partir des fichiers num riques de la BD CARTO de l IGN saisis l chelle du 1 50 000 L examen de la carte du r seau issue directement de la BD Carthage apporte un clairage et un compl ment d information quant la nature des formations lithologiques en particulier lorsque des formations superficielles mal renseign es sur les cartes g ologiques viennent modifier la perm abilit des roches sous jacentes Un calcul de la densit de drainage en km km par maille carr e de 2 5 km a permis d diter une carte rendant l information densit du r seau hydro graphique plus directement lisible figure 2 DONN ES CLIMATIQUES Pr cipitations annuelles moyenn
11. elief d crit par l altitude mais aussi de la g omorphologie pentes des terrains pentes des vall es formes du relief de l hydrographie la structure et l organisation du r seau de drainage donnent des informations compl mentaires aux donn es g ologiques perm abilit drainage phr atique e Le climat est appr hend par le r gime des pr cipitations annuelles saison ni res cart saisonnier et maximum d cennal sur 24 heures qui permet de caract riser les r gimes hydrologiques crues violentes les temp ratures maximales les types de v g tation naturelle potentielle ou les tages bioclimatiques qui int grent des infor mations sur les bilans hydriques et thermiques des informations compl mentaires notam ment sur les surplus ou d ficits hydriques les glaciers Autres documents utilis s L utilisation de couches d informations g ogra phiques en format num rique a t compl t e par le recours aux cartes notamment l Atlas des eaux souterraines de la France Margat 1970 les cartes g omorphologiques de la Fran ce coord F Joly GIP Reclus 1988 1987 1992 1993 collection Reclus Mode d Emploi de 1988 1993 la carte des surplus et d ficits hydrologiques de la France Trzpit 1984 scann e et g or f renc e des cartes de la v g tation de la France au 1 200 000 CNRS notamment pour le bassin du Rh ne la car
12. es cette carte r sulte du cumul sur l ann e des moyennes men suelles de pr cipitations p riode 1961 1990 Pr cipitations saisonni res les valeurs men suelles ont t cumul es par saison en groupant d cembre janvier f vrier pour l hiver mars avril mai pour le printemps etc C gences de l eau Densit en krmikrnid 70 03 3 07 E 07 1 1 1 4 Eu 10 F lt d Figure 2 Densit de drainage calcul e partir de la BD Carthage sur une grille maille carr e de 2 5 km France continentale 15 N 2 Q w 1 w l o N w pm Zg Nr bn 16 7 Direction r gionale de l Environnement 8 CORINE Co ORdi nation of INformation on the Environment Cartes de saisonnalit trois cartes exprimant les saisons de pr cipitations dominantes et minimales ainsi que l cart de pr cipitation entre saison ont t d riv es des pr c dentes Ces cartes se sont av r es pertinentes pour identifier les tendances de continentalit et surtout d limiter la r gion m diterran enne au sens climatique Carte des r gions phyto cologiques de Dupias et Rey CNRS 1980 cette carte bas e sur la distribution des associations phyto cologiques naturelles d coupe le territoire national en zones de quatre niveaux hi rarchiques diff rents Les deux premiers niveaux expriment les grandes structures g ologiques et clim
13. idors rivulaires pour caract riser la typologie de chaque station et d velopper des mod les pressions impacts Des d veloppements en cours n cessitant le recours des donn es plus pr cises MNT permettront de d terminer les caract ristiques physiques des tron ons de mani re permettre une extrapolation spatiale des mod les d habi tat poisson d velopp s au sein du laboratoire Lamouroux et Capra 2002 o Remerciements Ce travail a b n fici du soutien financier du minist re de l cologie et du D veloppement durable MEDD Direction de l Eau Il s inscrit dans le programme de la Zone atelier bassin du Rh ne ZABR Nous remercions B atrice Dutang du GIE VALOREZ pour son aide la pr paration du manuscrit Nous remercions galement Nicolas Mengin du GIE VALOREZ pour son travail sur les stations IBGN Constitution d un syst me d information r f rence spatiale dans le contexte de la Directive cadre Herv Pella et al R sum Pour r pondre l objectif principal de bonne qualit cologique des cours d eau de la Directive cadre europ enne sur l eau DCE un projet de r gionalisation des cosyst mes aquatiques a t initi La d termination des hydro cor gions HER qui constituent le cadre adopt pour la typologie des cours d eau a n cessit la mise en place d un syst me d information g ographique SIG d di Cet article pr sente la cr a
14. igure 3 Bassin versant calcul partir de la station et pressions anthro piques d riv es de l analyse de CORINE Land Cover 17 N 2 Q E w So w l o N w gm a O bn 18 D Figure 4 Exemple de calcul des pressions anthro piques d riv es de l analyse de Corine Land Cover l chelle du corridor rivulaire d une station Pressions Corine Land Cover EEE Territoire artificialis dense Agricole intensif Agricole faible impact EN For t _ Espace ouverts Rang de Strahler Ordre 1 MN Ordre 2 N SNS Ordre 5 N 0 1 2 km Conclusion Les travaux du laboratoire sur les bassins versants de la Loire et du Rh ne initi s dans les ann es 90 avaient d j n cessit le d veloppement des syst mes d informations g ographiques Pour r pondre aux objectifs de la Directive cadre le recours aux m mes syst mes s est donc fait dans la continuit des travaux ant rieurs La constitution de la plate forme SIG l chelle nationale telle qu elle est partiellement pr sent e dans ce document a demand un long et impor tant travail de collecte de mise en coh rence et de traitement des donn es g ographiques Cette base de donn e g ographique nous permet de faire le lien entre les stations biologiques les caract ristiques naturelles et les pressions anth ropiques s exer ant sur les bassins versants et les corr
15. ns sont r pertori es et g or f renc es Un travail similaire a t r alis sur les diatom es 832 stations Occupation du sol la couche d information g o graphique utilis e est CORINE Land Cover CLC Elle couvre l ensemble du territoire et repose sur une nomenclature standard hi rarchis e 3 niveaux et 44 postes r partis selon 5 grands types d occupation du territoire territoires arti ficialis s territoires agricoles for ts et milieux semi naturels zones humides et surfaces en eau La p riode d acquisition des images satellitales LANDSAT et SPOT qui ont servi l tablissement de cette couche va de 1987 1994 Pour caract riser les pressions sur le milieu aqua tique des typologies simplifi es pour les postes d occupation du sol sont r alis es pr alablement au traitement Ainsi les 44 postes initiaux ont t reclass s en huit cat gories 1 Territoire artificialis dense regroupant les types CLC 111 12 et 13 2 Territoire artificialis discontinu regroupant les types CLC 112 et 14 3 Agricole intensif regroupant les types CLC 21 22 241 et 242 4 Agricole faible impact regroupant les types CLC 231 243 et 244 5 For ts regroupant les types CLC 31 6 Pelouses landes et friches regroupant les types CLC 32 7 Espaces ouverts regroupant les types CLC 33 8 Zones humides regroupant les types CLC 4et 5 Exemple
16. nt les donn es qu il s agisse d un achat d une convention dans le cadre d une tude en pr cisant la date finale de la convention au del de laquelle les donn es doivent tre supprim es ou d une simple mise Constitution d un syst me d information r f rence spatiale dans le contexte de la Directive cadre Herv Pella et al disposition La partie qualit liste les gammes d chelles d utilisation la pr cision du trac la date de validit des informations ainsi que leur fr quence de mise jour Enfin le r f rencement pr cise le syst me de positionnement et l emprise g ographique Co ts et acc s aux utilisateurs Outre le co t des bases qui ont t acquises dans le cadre de cette tude principalement aupr s de l IGN et de M t o France le temps de travail important pour valider et mettre en coh rence ces donn es entra ne un co t suppl mentaire Ce SIG est mis jour r guli rement l acquisition du MNT 50 m issu de la BD Alti constituerait une volution importante pour des travaux futurs Cette plate forme SIG est actuelle ment un outil de recherche qui peut tre ouvert des projets du m me type Sa mise disposition d autres utilisateurs en dehors du domaine de la recherche serait du ressort d un organisme public d pendant du minist re de l cologie et du D veloppement durable ayant une vocation de diffusion de l information Dans un tel cadre le
17. rgeur du lit mineur Des travaux dans le bassin de la Loire Souchon et al 2000 ont montr la possibilit d extrapoler les largeurs moyennes des cours d eau partir des rangs de Strahler La largeur de la zone prise en compte est donc fonction du rang 100 m pour les cours d eau de rang 1 2 et 3 140 m pour ceux de rang 4 250 m pour ceux de rang 5 600 m pour ceux de rang 6 1 200 m pour ceux de rang 7 et enfin 2 400 m pour ceux de rang 8 Seuls les arcs de m me rang sont retenus pour le calcul de la zone tampon Il est donc possible lorsque la station est proche d une confluence qui change le rang en aval de ne consid rer que la portion du cours d eau de m me rang que celle sur laquelle est situ e la station Le calcul du type d occupation du sol est exprim en pourcentage Ce traitement a t automatis par l criture d un script en langage Avenue sous ArcView 3 2 Un exemple de r sultat est pr sent figure 4 p 18 tant donn que l chelle de travail de Corine Land Cover est le 1 100 000 l analyse pr sent e est en limite d utilisation Ce traitement permet d approcher l occupation dominante du fond de vall e L analyse de la ripisylve ou de l artifi cialisation des berges n cessiterait d utiliser une source de donn e plus pr cise Perez Correa 2004 Ea Toriai arilicialish derse Agricole inlereil Ei Agricole fuble impact EE For M o r i En es lt d F
18. s eaux courantes pour la Directive cadre europ enne sur l eau l approche par hydro cor gion Ing nieries EAT n 40 p 3 10
19. s temp ra Climat tures maximales et minimales pour la p riode 1961 1990 maille 1 x 1 km Aquisition Cemagref Carte des pr cipitations sur 24 h sur une p riode de retour de dix ans maille 1 x 1 km RS Carte des r gions phyto cologiques de Dupias V g tation et Rey 1985 z BD CARTHAGE V 3 0 d riv e de la BDCARTO ee R seau hydrographique de l IGN Aquisition Cemagref Tableau 1 Donn es g ographiques Ces donn es ont t exploit es avec les logiciels Arcview et Arcinfo modifi pour l coulement le plan de la direction d coulement le plan des surfaces amont drai n es la pente des vall es le relief ombr et le r seau hydrographique th orique Cette base compte une cinquantaine de cou ches ce qui repr sente pr s de 10 giga octets Les donn es sont organis es par th mes g olo gie morphologie hydrologie climatologie limi tes administratives occupation du sol pressions polluantes et r gionalisation Elle sert de support effectif tous les travaux de recherche multisites du laboratoire mise en relation de la biologie de la physique et de la chimie des cours d eau diff rentes chelles Elle est utilis e plus lar gement dans le cadre de la zone atelier bassin du Rh ne ZABR Catalogage Les couches d informations g ographiques qui constituent la plate forme SIG tant de prove nances diverses il est important d archiver les m tadonn es
20. te des glaciers Rogers et Wasson 1997 les cartes topographiques touristiques routi res dit es par l IGN 4 Mod le num rique de terrain 5 Produit dit par l Institut g ographi que national 6 Groupement d int r t collectif R seau d tudes des changements dans les localisations et les unit s spatiales Montpellier 13 A 2 S D v D l o n rc KD D bn 14 D Figure 1 Influence des roches sur la chimie de l eau Donn es d riv es et interpr t es partir de la carte g o logique BRGM au 1 1 000 000 6 dition Le traitement des donn es g ographiques STRUCTURE PHYSIQUE Altitude l information du MNT 250 m a t utili s e apr s regroupement en classes d altitudes Pente une carte des pentes l chelle nationale r alis e partir du MNT 250 m a galement t utilis e Relief ombr hillshade ce mode de transformation des donn es du MNT permet de visualiser les formes du relief et d identifier les discontinuit s existantes entre zones de morpho logie diff rentes Pente des vall es partir du MNT un r seau de drainage th orique a t reconstitu Les pentes des tron ons ont t calcul es sur ce r seau Le r sultat permet de visualiser grande chelle les valeurs dominantes de pentes de vall es et facilite ainsi le contr le de la discrimination op r e
21. tion du SIG l analyse des informations g ographiques pour d limiter les HER et un exemple d utilisation des donn es pour caract riser les pressions anthropiques Abstract To answer the main objective of the European Water Framework Directive WFD the good ecological status of water bodies a project of aquatic ecosystem regionalization was initiated The determination ofhydro ecoregions HER which constitute the adopted frame for water bodies typology required the implementation of a dedicated geographical information system GIS This article presents the GIS creation the analysis of geographical informations to delimit the HER and an example of data use to characterize the anthropic pressures Bibliographie BONNEFOY D CHANTRAINE J LEBRET P RABU D 1996 Carte g ologique de la France 1 1 000 000 num rique Guide d utilisation rapport BRGM R39145 33 p BRGM Carte g ologique au 1 1 000 000 6 dition 1996 BRGM Notice de la carte g ologique au 1 1 500 000 dition de 1980 BRGM Cartes g ologiques au 1 250 000 Valence Nice Annecy Marseille Gap Dijon Paris Chalons sur Sa ne Thonon les Bains Corse BRGM Cartes g ologiques diverses DUPIAS G REY P 1985 Document pour un zonage des r gions phyto cologiques CNRS 39 p cartes MARGAT J ed 1970 Atlas des eaux souterraines de la France BRGM DATAR ditions du BRGM Orl ans non pagin LAMOUROUX N
Download Pdf Manuals
Related Search