# Which spatial discretization for which distributed hydrological model?

J. Dehotin, Isabelle Braud

<figure><img src="image_1.png" /></figure>

## To cite this version:

J. Dehotin, Isabelle Braud. Which spatial discretization for which distributed hydrological model?. Hydrology and Earth System Sciences Discussions, 2007, 4 (2), pp.777-829. hal-00330796

HAL Id: hal-00330796

https://hal.science/hal-00330796v1

Submitted on 18 Jun 2008

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

Papers published in *Hydrology and Earth System Sciences Discussions* are under
open-access review for the journal *Hydrology and Earth System Sciences*

Which spatial discretization for which
distributed hydrological model?

J. Dehotin and I. Braud

Cemagref Lyon, Research Unit Hydrology and Hydraulics, 3bis Quai Chauveau, CP 220,
69336 Lyon Cédex 9, France

Received: 12 January 2007 – Accepted: 12 February 2007 – Published: 10 April 2007

Correspondence to: J. Dehotin (dehotin@lyon.cemagref.fr)

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>
Abstract

Distributed hydrological models are valuable tools to derive distributed estimation of water balance components or to study the impact of land-use or climate change on water resources and water quality. In these models, the choice of an appropriate spatial scale for the modelling units is a crucial issue. It is obviously linked to the available data and their scale, but not only. For a given catchment and a given data set, the “optimal” spatial discretization should be different according to the problem to be solved and the objectives of the modelling. Thus a flexible methodology is needed, especially for large catchments, to derive modelling units by performing suitable trade-off between available data, the dominant hydrological processes, their representation scale and the modelling objectives.

In order to represent catchment heterogeneity efficiently according to the modelling goals, and the availability of the input data, we propose to use nested discretization, starting from a hierarchy of sub-catchments, linked by the river network topology. If consistent with the modelling objectives, the active hydrological processes and data availability, sub-catchment variability can be described using a finer nested discretization. The latter takes into account different geophysical factors such as topography, land-use, pedology, but also suitable hydrological discontinuities such as ditches, hedges, dams, etc. For small catchments, the landscape features such as agricultural fields, buildings, hedges, river reaches can be represented explicitly, as well as the water pathways between them. For larger catchments, such a representation is not feasible and simplification is necessary.

For the sub-catchments discretization in these large catchments, we propose a flexible methodology based on the principles of landscape classification, using reference zones. These principles are independent from the catchment size. They allow to keep suitable features which are required in the catchment description in order to fulfil a specific modelling objective. The method leads to unstructured and homogeneous areas within the sub-catchments, which can be used as modelling units. It avoids map

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>
<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

smoothing by suppressing the smallest units, the role of which can be very important in hydrology, and provides a confidence map (the distance map) for the classification. The confidence map can be used for further uncertainty analysis of modelling results. The final discretization remains consistent with the scale of input data and that of the source maps.

We present an illustration of the method using available data from the upper Saône catchment (11 700 km²) in France. We compare the results with more traditional mapping approach, according to the landscape representation and input data scale.

# 1 Introduction

Due to a growing concern about environmental and climate change issues, and the emergence of the concept of sustainable development, the requirements towards hydrology have evolved from the prediction of the water streamflow at a few locations, to the prediction of the water balance components (rainfall, runoff, water storage, transpiration, evaporation, groundwater levels etc.) at every point within a catchment. The consideration of land-use and human-induced modifications of landscapes is a major concern for water management purposes (quantity and quality), especially for flood forecasting, the study of the impact of land management on streamflow, pollutants or sediments transport. To reach these objectives, the explicit treatment of the land-surface heterogeneities is necessary. Furthermore, for some questions such as pollution, the knowledge of the water balance components is not sufficient and fluxes throughout the landscape are required as well as a proper handling of water pathways. In this context, distributed models can be valuable tools as they have the ability to take the landscape heterogeneity into account, and provide distributed output variables. They can be used to test different functioning hypotheses from which simplified and/or predictive models can be derived for more operational purposes. If reliable output variables are expected, model parameters should be estimated a priori from available information. Verification of the model behaviour should be performed not only on the
<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

streamflow at the outlet, but also at intermediate stations and on other variables leading to a multi-objective verification (e.g. Mroczkowski et al., 1997; Beldring, 2002; England et al., 2006; Varado et al., 2006). There was high expectation about the ability of distributed models to take into account changes in the landscape, especially thanks to the increasing availability of high-resolution information. However, the specification of all needed parameters, with a suitable resolution, remains a difficult and uncertain task. Thus, some parameters are often calibrated in practice and the usefulness of such models has been questioned many times due to problems of overparameterization, parameter estimation and validation limitations (e.g. Beven and Binley, 1992; Grayson et al., 1992a, b; Beven, 2001).

One major point of concern today is the definition of an appropriate scale for the modelling units and the choice of a proper representation of hydrological processes at this scale. Although a subject of active research, it still remains an open question, which has not yet been resolved. These questions form part of the challenge proposed by various hydrologists (e.g., Beven, 2002a; Reggiani and Schellenkens, 2003) or in the framework of the Predictions in Ungauged Basins (PUB) decade, initiated by the IAHS (e.g., Sivapalan, 2003). Underlying all these questions is the problem of scales and scaling in hydrology as reviewed by Blöschl and Sivapalan (1995).

In this paper we focus on a discretization strategy which allows a better assessment of the question: which spatial discretization for which distributed hydrological models conceptualisation? The definition of such appropriate spatial scale cannot be given once and for all. The answer should be a function of the catchment under study and of its specificities, and will result from a trade-off between various, sometimes antagonist considerations:

1. What is the objective of the distributed hydrological modelling?

2. Which output variables are required and at which spatial and temporal scales?

3. What are the measured data and at which resolution are they available?
<figure><img src="image_20.png" /></figure>

<figure><img src="image_21.png" /></figure>

<figure><img src="image_22.png" /></figure>

<figure><img src="image_23.png" /></figure>

4. What are the active/dominant hydrological processes on this catchment and what are their functional scales?

5. Which representation of hydrological processes is relevant and at which scale?

6. Which degree of heterogeneity is acceptable within the modelling units?

5

**2 Definition of an appropriate scale for distributed models: an overview**

As said above, the definition of an appropriate spatial scale for hydrological models will result from a trade-off between various considerations. Let's review the various items listed before.

2.1 What is the objective of the distributed hydrological modelling?

10

This question might seem trivial but is often not always well defined by the modellers. Refsgaard et al. (2005) retained it as one of the first item in the list of tasks they identified for performing a modelling study while respecting some quality insurance criteria. Examples of possible modelling objectives are: determination of the components of the water balance of a catchment, quantification of flood or draught risk, evaluation of mitigation solutions for the limitation of the river pollution at a given location, test of functioning hypotheses, search for dominant hydrological processes... For each objective and a given catchment, the required spatial discretization should be different.

15

2.2 Which output variables are required and at which spatial and temporal scales?

20

The definition of the objective of the study leads to the definition of the output variables of interest, and of the spatial and temporal scales at which they are required. For water balance analysis, the outputs are the various components such as runoff, evapotranspiration, groundwater recharge, water storage etc. For pollutant transfer problems,
<figure><img src="image_24.png" /></figure>

<figure><img src="image_25.png" /></figure>

<figure><img src="image_26.png" /></figure>

<figure><img src="image_27.png" /></figure>

outputs can be, integrated fluxes, maximum concentrations over a given time step, du-
rations over which legal thresholds are exceeded. These outputs can be required at the
annual, monthly, daily, hourly time scale; at distributed locations within sub-catchments
or only at the outlet of the whole catchment. It is generally admitted that coarser time
and space scales require coarser modelling units, but there is no clear rule to define
that appropriate scale.

2.3 What are the measured data and at which resolution are they available?

The measured data includes input forcing variables (rainfall and climate forcing), output variables for verification such as streamflow, soil moisture, groundwater levels, surface fluxes, and landscape descriptors such as land use, geology, elevation, which are used for model parameter estimation and definition of the modelling units.

When speaking about observations, Blöschl and Sivapalan (1995) distinguish be-
tween the extent of the data, i.e. the zone other which the data set is collected; the
support of the observation, i.e. the spatial scale at which data are integrated and the
spacing, i.e. the distance (in time and space) between different observations points.
In this discussion, we will focus on the spatial scale, but the temporal scale is obvi-
ously linked (e.g. Skoïen and Blöschl, 2006). For modelling purposes, all input data,
parameters and the verification data are required over the modelling units. Unfortu-
nately, there is often a mismatch between the observations scales and the modelling
scales. The support of in-situ measurements is often punctual. Therefore spatial in-
terpolation using techniques such as kriging is required to derive the values over the
modelling units. Streamflow data are directly integrated over catchment areas and are
more consistent with the modelling scale, but the number of gauging stations is often
much smaller than the number of modelled sub-catchments. New measurements tech-
niques such as scintillometers provide a certain space averaging for sensible and/or
latent heat flux along transects (Green et al., 2001), but do not yet provide values over
sub-catchments.

With the availability of remote sensing data (satellite, planes, drone etc), there was a

5

10

15

20

25
<figure><img src="image_28.png" /></figure>

<figure><img src="image_29.png" /></figure>

<figure><img src="image_30.png" /></figure>

<figure><img src="image_31.png" /></figure>

hope that their spatial nature will help in filling the gap between in-situ measurements and modelling. Existing sensors are more and more able to get data at impressive resolutions, but according to the resolution, some variables do not have the same meaning (e.g. slope, leaf area index etc). Moving from information at a fine resolution (e.g. the trees) to information at a larger scale (e.g. a forest) is not obvious. Thus the fundamental question “which resolution is needed for landscape description to represent hydrological processes?” remains an open question (Puech, 2002). Furthermore, remote sensing measurements are often not directly related to the hydrologic quantity of interest. The sensors often only sample the first cm of the continental surfaces, whereas information on deeper soil layers would be required. Multi-disciplinary research amongst hydrologists (and more generally with researchers in environment) and remote sensing specialists is still needed to progress on these questions.

There is also a paradox in this progress of remote sensing: whereas the continental surface can be described with more and more accuracy, even at the scale of a building, the knowledge of the sub-soil is not progressing so rapidly. For example, there is a lack of knowledge on soil properties. The data support is local, spacing and extend are limited due to the difficulty and cost of soil sampling. To derive reliable maps at larger scales, hypotheses about soil organisation, forming factors are required. In pedology sophisticated classification techniques using geostatistics or fuzzy rules are developed for mapping soil units, but the result remain uncertain (e.g. Burrough et al., 1997; Lagacherie et al., 1997, 2001). Lots of information is available in soil databases. But they often provide rather descriptive information on the soils. Their content is often found disappointing, if not useless by hydrologists who are looking for transfer coefficients. The initiative of Lin et al. (2006) to promote hydropedology as a synergetic discipline between hydrologists and pedologists is promising but will require some years to be fruitful. Progress are also expected on soil characterization with the use of geophysical techniques but their use is still limited and they cannot be deployed routinely yet. Therefore, for some times on, we will have to cope with the paucity of information related to the soil and sub-soil, when performing hydrological modelling. It should not

5

10

15

20

25
<figure><img src="image_32.png" /></figure>

<figure><img src="image_33.png" /></figure>

<figure><img src="image_34.png" /></figure>

<figure><img src="image_35.png" /></figure>

be forgotten when combining these data with the detailed data derived from remote sensing of continental surfaces.

From this analysis of data in hydrology, it is obvious that the information is often not available with the desired accuracy. Therefore, hydrologists should perhaps change their way of thinking: instead of looking for data for a given model, they should perhaps try to find the relevant formulation of hydrological processes for given available data in order to avoid overparameterization and uncertainty.

## 2.4 What are the active/dominant hydrological processes on this catchment and what are their functional scales?

A lot of former distributed hydrological models were based on Hortonian scheme for the runoff. In such models, the catchment was subdivided into so-called isochrones surfaces. The hydrograph separation with isotopic techniques showed later that this mechanism was uncommon at a flood event scale (Crouzet, Hubert et al., 1970, cited by Gineste, 1997). Frequently, the hydrograph is mainly composed with the waters present in the soil before the rainy event (Grésillon, 1994). This convinces hydrologists of the complexity across scale of flow generation processes inside a catchment.

The question of process scale in hydrology has been reviewed by Blöschl and Sivapalan (1995) who proposed a diagram of characteristic spatial and temporal scales for the various processes in a log-log plot (Blöschl and Sivapalan, 1995, Fig. 2). Lots of research has been dedicated to the determination of these characteristic scales. One of the recent examples is provided by Skoïen et al. (2003) and Skoïen and Blöschl (2006) who analysed rainfall, streamflow, groundwater level and soil moisture records from Austria and Australia, using geostatistical tools. They were able to determine characteristic time and scales for rainfall and runoff of respectively one day and one month, while showing that groundwater levels were not stationary. In space, they found that rainfall was almost fractal without characteristic scales whereas runoff appeared non-stationary but not fractal. This data analysis provided data evidence of the consistency of the Blöschl and Sivapalan (1995) scale diagram. The latter provides guidance for the
<figure><img src="image_36.png" /></figure>

<figure><img src="image_37.png" /></figure>

<figure><img src="image_38.png" /></figure>

<figure><img src="image_39.png" /></figure>

definition of appropriate spatial and temporal discretizations of the various hydrological processes that are included into a particular model.

However, the range of scales for a given process is still large and the dominant processes can change with scale. This is especially evident for the rainfall-runoff relationship for which various authors have shown a decrease of the runoff coefficient with increasing catchment scale for various hydrologic and climatic contexts (e.g., Bergkamp, 1998; Braud et al., 2001; Cerdan et al., 2004). A downward approach of model complexity (Klemes, 1983), based on data analysis can help in the formulation of a conceptual model of the rainfall-runoff relationship, leading to a parsimonious model using parameters derived from available data. This concept has been recently applied by Jothityangkoon et al. (2001) and Eder et al. (2003) for a semi-arid and an alpine catchment respectively. They propose models of the rainfall-runoff relationships at the annual, monthly and daily time scale, by progressively increasing the model complexity until a good reproduction of the data behaviour was obtained. The two cases studies show that, according to climate and catchment characteristics, very different models can emerge, with different dominant hydrological processes in both cases. Macropores, preferential flow, re-infiltration, variability of land cover, the influence of micro-topography leading to the concentration of runoff into small channels, are some of the factors being able to explain threshold effects and the observed differences in dominant hydrological processes across scales (e.g. Bergkamp, 1998; Lehmann et al., 2006; Sidle, 2006; Zehe et al., 2005).

Several concepts have been proposed to describe and explain the variability of landscape characteristics such as organization, hierarchy or fractal behaviour, leading to the definition of various characteristic scales (e.g. Vogel and Roth, 2003; Lin et al., 2006). However, there is still a long way to go before a clear picture can be drawn. Nevertheless, when defining the hydrological modelling units, this multiplicity of scales must be acknowledged and, as much as possible, taken into account.
<figure><img src="image_40.png" /></figure>

<figure><img src="image_41.png" /></figure>

<figure><img src="image_42.png" /></figure>

<figure><img src="image_43.png" /></figure>

## 2.5 Which representation of hydrological processes is relevant and at which scale?

The representation of a process within a model implies the choice of two complementary elements. The first one is the resolution of the model (the spatial scale we are discussing in this paper) and the process conceptualisation. If we borrow the vocabulary from the atmospherical community, the choice of the resolution of the model will allow the separation between the processes which are represented explicitly (i.e. for which a prognostic variable with an evolution equation is defined) and the sub-grid scale processes which will be parameterised<fnref n="1" /> (i.e. for which no evolution equation is solved and for which simplified representations are adopted and added to the prognostic equations). According to the resolution of the model, some processes can be represented explicitly or parameterised. By adopting the same vocabulary in hydrology we would avoid the quarrel on the nature of conceptual or physically based modelling. A process would be conceptual (or parameterised) at one resolution and physically based (explicitly resolved) at another resolution.

Plot scale studies allow the derivation of physical equations, extensively used in hydrology such as the Richards equation for saturated/non-saturated water flow, the Boussinesq equation for 2-D groundwater flow, the Saint-Venant equation for river or overland flow. At this scale, several parameters in these equations such as retention curve, hydraulic conductivity or surface roughness can be estimated from measurements. When moving to larger scales, it is often assumed that the form of the equation remains valid. Then it is necessary to derive so-called effective parameters at those scales (e.g. Blöschl and Sivapalan, 1995). For this purpose the aggregation-disaggregation modelling (Blöschl and Sivapalan, 1995) approach to identify the functional relationship at a larger scale from results at smaller scales can be used (see an example in Viney and Sivapalan, 2004).

<footnote n="1">¹In this context, the word “parameterization” should not be confounded with the estimation of parameters for which it is often used in hydrology. It would be equivalent to conceptualisation in the hydrology jargon.</footnote>
<figure><img src="image_44.png" /></figure>

<figure><img src="image_45.png" /></figure>

<figure><img src="image_46.png" /></figure>

<figure><img src="image_47.png" /></figure>

Several distributed hydrological models are based on a regular mesh over which point scale laws are extended and where effective values of the parameters must be determined. Examples are MIKE-SHE (Abbott et al., 1986a, b), ECOMAG (Motovilov et al., 1999), TOPKAPI (Ciarapica and Todini, 2002). Some authors contest this approach, referred to as a “reductionist” approach (e.g. Gottschalk et al., 2001), arguing that the equation becomes a parameterisation of the process, since parameters cannot be estimated from field measurements (e.g. Beven, 2002b). The choice of the grid size is not always rationalised taking into account the processes that are represented, but seems rather the result of commodity and data resolution. One exception can be found in Beldring et al. (1999) and Motovilov et al. (1999) using the ECOMAG model in the framework of the NOPEX project. They determined the size of the mesh from analysis of averaging properties of point groundwater and soil moisture measurements obtained using a dedicated sampling strategy with nested spacing. Of course, data required for such a study are seldom available.

One of the main criticisms about square elements is their poor handling of heterogeneity because continental surface is not organised in pixels. The task of parameter estimation is therefore more difficult. To overcome the problem, some authors proposed different approaches such as meshes based on iso-contours of elevation in THALES (Grayson et al., 1992a) or TOPOG (Vertessy et al., 1994); and more recently meshes based on Triangular Irregular Networks (TINs) (Ivanov et al., 2004a; Vivoni et al., 2004). The latter offer a good compromise between efficiency and accuracy as shown by the performances of the tRIBS model, developed on this irregular geometry (Ivanov et al., 2004a, b). Other authors tried to define more “hydrological” modelling units. Contributive zones are based on the concept of hydrologic similarity, and can be defined using for instance the topographic index of Beven and Kirkby (1979). Within these areas, it is assumed that the catchment response is similar. The concept is used in TOPMODEL (Beven and Kirby, 1992). In order to represent land-use heterogeneity, some authors have introduced the concept of Hydrological Response Units (HRUs) (e.g., Flügel, 1995), used in the Soil Water Assessment Tool (SWAT) model (Neitsch et
<figure><img src="image_48.png" /></figure>

<figure><img src="image_49.png" /></figure>

<figure><img src="image_50.png" /></figure>

<figure><img src="image_51.png" /></figure>

al., 2005). HRUs represent a sub-catchment scale discretization composed of a unique combination of land cover, soil and land management. One of the drawbacks is that the HRUs mapping induces merging of smaller units into larger ones by applying smoothing filters. From a hydrological point of view, it may result in a loss of information, as some major hydrological processes can be localised on very small units. Illustrations are re-infiltration of runoff at the bottom of hillslopes in the Sahel (Seguis et al., 2002), runoff decrease due to hedge networks (Viaud et al., 2005).

Another example of hydrological spatial discretization is the concept of Representative Elementary Area (REA) proposed by Wood et al. (1988), looking for characteristic spatial scales, beyond which the geographical locations of features could be neglected and the distribution taken into account using statistical distributions. Fan and Bras (1995) questioned the universality of the concept, especially because flow routines and the hierarchical structure of the river network were not taken into account in the analysis. This drawback is overcome with the concept of Representative Elementary Watershed (REWs) proposed by Reggiani et al. (1998, 1999, 2000) and extended by Tian et al. (2006). REWs form the elementary modelling units divided into several zones corresponding to the various hydrological processes. Global mass, momentum and energy balance laws are formulated at the sub-catchment scale. The corresponding equations remain unchanged whatever the scale (e.g. for REWs defined at various Strahler order). On the other hand fluxes between REWs and their zones (saturated, non-saturated, overland, concentrated and river flow) must be defined for each scale. Sub-catchment scale variability can be parameterised in the derivation of these fluxes. The strength of the approach is therefore to translate the general problem of model formulation into the problem of derivation of closure relationships at various scales (Reggiani and Schellekens, 2003). Lee et al. (2005), Reggiani and Rientjes (2005) or Zhang et al. (2005) have provided various formulations of these closure relationships.
<figure><img src="image_52.png" /></figure>

<figure><img src="image_53.png" /></figure>

<figure><img src="image_54.png" /></figure>

<figure><img src="image_55.png" /></figure>

## 2.6 Consequences for the definition of modelling units and process representation: a practical proposition

The picture drawn from the review of the various items above might result quite confusing. Contrarily to atmospheric sciences, the definition of a unique scale, separating processes being represented explicitly from those that must be parameterised, is difficult due to the hierarchical nature of the river network and the landscape complexity across scale. Furthermore characteristic scales are different for the various processes. Therefore we propose to adopt an approach, as quoted by Leavesley et al. (2002), requiring that, for a given catchment, the question should not be “which model is most appropriate for a specific set of criteria?” but “what combination of process conceptualisations is most appropriate?”. This approach is consistent with the downward approach mentioned previously and the recognition of the “uniqueness of place” as stated by Beven (2003). It also allows building a specific model for a specific objective, taking into account the availability of data. This pleads for the use of multi-scale hydrological framework, where the processes are develop as independent components, using the facilities provided by Object Oriented Modelling and, if possible with their characteristic time and space scales. They are then coupled through adequate tools provided by the modelling environment (for recent reviews about environmental computing frameworks, see for instance Argent, 2004; Krause et al., 2005).

In order to represent landscape heterogeneity efficiently according to the modelling goals, we propose a flexible methodology for catchment meshing based on the REW discretization. It uses nested discretizations, starting from a hierarchy of sub-catchments, linked by the river network topology. If consistent with the modelling objectives, the active hydrological processes and data availability, sub-catchment variability can be refined using finer discretizations (Representative Elementary Columns, RECs) in order to redefine the estimation of exchange fluxes within the REW. The discretization can take into account different geophysical factors such as topography, land-use, geology, pedology, but also hydrological discontinuities such as ditches and
<figure><img src="image_56.png" /></figure>

<figure><img src="image_57.png" /></figure>

<figure><img src="image_58.png" /></figure>

<figure><img src="image_59.png" /></figure>

hedges, etc., in order to represent sub-catchment variability. These RECs can also be discretized to take into account the vertical structure of soil profiles (Haverkamp et al., 2004). Conservation laws can be solved on the obtained elementary volumes, with various degrees of complexity. At a fine resolution, we can end up with standard partial differential equations. The conservation laws can be solved using finite volumes methods, which are consistent with the exchange fluxes approach. At larger scales, simplified representations can be derived.

Practically, for small catchments, the landscape objects such as agricultural fields, buildings, hedges, river reaches can be represented explicitly, as well as the water pathways between them (e.g., Moussa et al., 2002; Carluer and De Marsily, 2004, Branger, 2007). For larger catchments, such detailed representation is not feasible and simplifications are necessary. The discretization methodology should remain flexible enough to fit into the modular concepts exposed above and as objective as possible concerning the consistency of scales and the simplifications of the patterns. As a practical solution, we suggest to extend the principle of landscape classification, used is soil mapping for the definition of "hydro-landscapes" proposed by Winter (2001). It allows defining different levels of complexity in landscape representation that can be associated with different levels of complexity in hydrological processes description. The factors used in the classification can be adapted to the dominant hydrological processes and the size of the final units remains consistent with the resolution of the available data. The methodology for hydro-landscape classification is described in Sect. 3 and an illustration on the upper Saône catchment (11 700 km²) is proposed in Sect. 4. Comparison with traditional methods is also performed in Sect. 4, with a discussion of the interest and limitations of the approach. Conclusions and perspectives are given in the last section.
<figure><img src="image_60.png" /></figure>

<figure><img src="image_61.png" /></figure>

<figure><img src="image_62.png" /></figure>

<figure><img src="image_63.png" /></figure>

## 3 Discretization of a catchment using landscape classification technique

We borrowed the principles of our landscape classification from those used in soil mapping (Burrough et al., 1997) and more specifically from the method of Robbez-Masson (1994, 1995); Robbez-Masson et al. (1999); Lagacherie et al. (2001). The latter is based on the definition of reference zones and an analysis of the neighbourhood composition at each location.

### 3.1 Hydro-landscape definition and classification requirements

Hydro-landscapes can be defined as areas where hydrological processes can be considered as homogeneous. Their delineation can take into account what will be referred to as factors below, influencing hydrological processes, e.g. slope, land use, geology, pedology, etc.

Hydro-landscapes can be considered as an extension of the HRU concept, as we added the following requirements for the delineation methodology. First, the method should avoid the quite arbitrarily procedures consisting in removing small areas from the map resulting from the overlay of the factors maps. Second, the method should allow better control on the errors arising from the overlay of maps at different resolutions, in an objective and quantifiable way.

### 3.2 Principles of hydro-landscape delineation

The different steps of the method are the following:

1. Definition of the hydrological modelling objectives,

2. Identification of the available data and of their resolution. Note that the procedure could stop here if the data are not sufficient to fulfil the objectives. Various maps describing several catchment characteristics can be use as factors (slope, soil, lithology, land-use...). The classification is a raster-based method. Thus all map files must be in a raster format and are resampled with the same resolution.
<figure><img src="image_64.png" /></figure>

<figure><img src="image_65.png" /></figure>

<figure><img src="image_66.png" /></figure>

<figure><img src="image_67.png" /></figure>

3. Choice of relevant hydrological processes and of their representation according to the available data;

4. Definition of the *p* factors expected to be influential on these hydrological processes, according to the chosen representation,

5. Simplification of each *p* factor map into a number *n*<sub>*p*</sub> of classes. Their combination using GIS leads to a multivariate map of the Π*n*<sub>*p*</sub> combined factors. One class is therefore a unique combination of the *p* factors (Fig. 1);

6. Definition of the reference zones on the study catchment. They are areas representing a unique combination of the retained factors. It may be a specific area of interest. These reference zones are characterized using the combining factor map and the neighbourhood descriptor defined above.

7. Choice of a neighbouring window (size and shape) and a descriptor of the composition of the combined factor within this neighbourhood;

8. All points in the landscape are characterized by their neighbourhood window composition.

9. Mapping of the whole catchment using a pixel-by-pixel analysis, i.e. an allocation of each pixel to one of the reference zones according to a distance criterion between the descriptor of its neighboring window composition and that of the reference zones;

10. Estimation of the distance map;

11. Iteration from steps 7 to 10 until a resolution ensuring a given accuracy and consistency with the input data is obtained.

We will present an illustration of the whole methodology in section 4 and will now detail steps 4 to 10.

5

10

15

20
<figure><img src="image_68.png" /></figure>

<figure><img src="image_69.png" /></figure>

<figure><img src="image_70.png" /></figure>

<figure><img src="image_71.png" /></figure>

## 3.3 Definition of landscape factors (step 4 to 5)

In natural sciences, there is not a unique definition of a “landscape” amongst disciplines. In agriculture, agro-landscape refers to an ensemble of fields that are classified according to natural vegetation, wooded zones, river network, topography, soil surface characteristics (Girard, 2000). In soil sciences, pedology is considered as the results of several interacting factors (climate, geology, slope land use parameters) which are used to define soil-landscapes. In ecology, eco-regions are defined as land and water extends including distinct natural community.

For hydrology, the factors that will be retained in the analysis will depend on the modelled hydrological processes. To model infiltration, factors such as soil surface characteristics, soil types, management practices will be influential. To simulate runoff, we can consider topography, topographical index, and soil surface characteristics. For evaporation modelling, land-use, orientation, groundwater levels (geology), snow melt (topography, orientation) can be taken into account. Once all needed factors are identify, their corresponding layers/maps superposition using GIS gives a composed picture of the landscape. This composed picture defines various combinations of the landscape factors characterizing the spatial organization of water dynamic in the catchment.

## 3.4 Definition of reference zones (step 6)

The notion of reference zone has been borrowed from soil classification where Favrot (1989) proposed the use of reference sectors for small natural region characterisation. Reference zones are extensively surveyed areas that are supposed to contain all the soil classes of the region to be mapped. For the final cartography, one unique soil type characterized each reference sector. The derivation of accurate classifications and the quantification of uncertainties has led to lots of research using techniques such as geostatistics, fuzzy sets, conditional probabilities (e.g. Burrough et al., 1997), the discussion of which is beyond the scope of the paper. We propose to extend the notion
<figure><img src="image_72.png" /></figure>

<figure><img src="image_73.png" /></figure>

<figure><img src="image_74.png" /></figure>

<figure><img src="image_75.png" /></figure>

of reference area for hydrological purposes. We define a reference zone for hydrology
as an area with a certain degree of homogeneity relevant in the characterisation of
hydrological processes. It may correspond to a unique combination of factors or to a
dominant combination of factors. For a modelling of evapotranspiration and runoff such
combination can be “deciduous forests over steep slopes” or “cultures over moderate
slopes”. For modelling saturated zones and their influence on runoff, they can be “high
topographic index” etc. Man influence can be included through management practices:
“drained area”, “urban area”, etc. References zone can also be a particular natural
region where water dynamics are specific (e.g. karstic areas, vineyard). There is no
limitation in the number of reference zones considered in the analysis. Therefore, the
modeller has a lot of flexibility in the mapping of the landscape heterogeneity (com-
plex/simple description), according to the data available, their scale and its objectives.

If there is a good knowledge of the catchment, the reference zones can be quite eas-
ily defined and delineated using the a priori knowledge of field hydrologists. The latter
can know the location of areas prone to saturation, or with high slopes, etc. For larger
catchments or catchments where it is not possible to perform intensive field surveys,
the task is more difficult. In this case, only the factors maps can provide the available
information. These maps can be used to define possible reference zones, according
to traditional/standard knowledge about hydrological processes. In this case, a simpli-
fied classification of factors can be used and the reference zones can be defined using
a statistical and spatial analysis of the multivariate map (see details in Sect. 4). The
relevance of such delineation can then be confirmed by a specific field survey.

3.5 Mapping of modelling units (steps 7 to 10)

Neighbourhood definition and characterization of all point in the catchment (step 6 and 7)

"Landscape is what is around you" (Robbez-Masson, Foltête et al., 1999). With this sentence these authors argue that the integration of spatial neighbourhood is neces-
<figure><img src="image_76.png" /></figure>

<figure><img src="image_77.png" /></figure>

<figure><img src="image_78.png" /></figure>

<figure><img src="image_79.png" /></figure>

sary to describe the landscape. With this mapping approach, each point in the land-
scape is characterized by the composition of a neighbourhood window around the point
(contextual analysis). The modeller must define the size and the shape of a neighbour-
hood window (e.g. ellipsis neighbourhood window on Fig. 1). The latter determines
the resolution of the final units. All points in the catchment are characterized using the
composition of their neighbourhood window.

The characterization of each point is performed (inside the neighbourhood window) using a descriptor on the multivariate image (in pixels). The descriptor may be a histogram, a mean or a standard deviation. This descriptor is calculated for each of the points of the multivariate image. In the example provided in the left of Fig. 2, each colour of the pixels characterizes a specific combination of factors and we chose the histogram as a descriptor. To derive it we consider all the pixels inside the neighbourhood window (here a square) and count the number of each combination class (i.e. each colour) in this window. This histogram is constructed for each point of the map.

The choice of the shape of the window allows taking into account some anisotropy in
the catchment. For instance ellipsis neighbourhood window with the major axe oriented
north south can be chosen if specific factor (e.g. topography) has such an orientation.

References zones characterization and mapping procedure (steps 8 to 10)

Similarly to each point of the map, the references zones are characterized by their
composition, using the same descriptor, according to the combination factor map (right
of the Fig. 2).

The mapping consists in the assignment of all points in the landscape to the most
similar reference zone in a statistical sense. Similarity is defined as the minimization
of the distance between the descriptor of each characterized point and those of the
references zones. The distance may be the modal distance, Kolmogorov distance or
Manhattan distance (Robbez-Masson, 1994). See the illustration of the principles in
Fig. 2.
<figure><img src="image_80.png" /></figure>

<figure><img src="image_81.png" /></figure>

<figure><img src="image_82.png" /></figure>

<figure><img src="image_83.png" /></figure>

The result of the mapping is a first map representing a segmentation of the initial multivariate image into elementary landscape units or polygons. The resulting map is thus composed of homogeneous and unstructured areas. A second image is obtained and can be considered as a confidence map for the classification. It quantifies the reliability of the modelling units by providing a criterion of the statistical quality of the classification. If the confidence map is not satisfactory, the classification can be improved by adding more reference zones to get a better representation of the landscape.

Scale and accuracy assessment (step 11)

It is possible to get different units size for the final landscape map, by using different size for the neighbourhood windows. This size must be chosen in consistency with the scale of the input data. The size of the smallest units on the classified map cannot be lower than the finest units on the input maps. An iterative procedure is therefore needed to define a neighbourhood window size, consistent with this first constraint. This iterative procedure will consist in testing several sizes for the neighbourhood window until the constraint is fulfilled (iteration of steps 6 to 10). This ensures consistency of the modelling units with the input data scales.

Once the “optimal” size of the neighbourhood windows is chosen, the classification can be improved using the distance map. The latter provides an idea of the accuracy of the classification. If the distance is large, it means that the similarity of the neighbourhood with the available reference zones is poor. Therefore, new reference zones can be added in the areas with larger distances and used to improve the mapping (iteration of steps 6, 8 to 11). This will reduce uncertainties on the landscape heterogeneity representation and handling.

5

10

15

20
# 4 Illustration of the methodology using data from the upper-Saône catchment in France

In this section, we illustrate the methodology outlined above using data from the upper-Saône catchment upstream the gauging station Lechatelet (11 700 km²). It is located in north-eastern France (Fig. 3). The river rises in the Vosges mountains in Lorraine, and flows south through Burgundy. The elevation ranges from 177 a.m.s.l. at the outlet to more than 1215 a.m.s.l. in the Vosges. The mean annual precipitation varies from 790 mm close to the catchment outlet to 2440 mm in the Vosges. Snowfall constitutes a minor part of the annual precipitation. The coldest month is January (–0.4–2.6°C) whereas the warmest is July (16.1–21.0°C). The average annual potential Evaporation calculated by the Pennman equation varies from 859 mm in the lowest part to 717 mm in the highest part. The highest values appear in July (125–157 mm) and the lowest in January (19–26 mm). This climate results in a runoff regime called Pluvial C regime in the Pardé (1955) classification (Sauquet et al., 2000). This is a regime where the runoff is rain-fed with winter high water and low water in summer due to high evaporation transpiration losses. The mean annual runoff varies from 323 mm in the lowest parts, to more 1451 mm in the Vosges. The main land-use classes are broad-leaved forests, arable land, and pastures. The arable land is mainly located in the south-western parts, whereas the pastures are located in the north-eastern parts. The forests are spread all over the catchment, and in the Vosges the coniferous forest is important. Urban surfaces constitute about 2% of the catchment area with the largest concentration around the town of Dijon. The geology is characterised by limestones in the southern parts and sandstones and granites in the northern parts.

We first present the available data and their resolution. Then, we illustrate the nested discretization procedure and all the steps outlined in Sect. 3. We assume that the modelling objective is to derive over the whole catchment, the water balance components at different temporal scale (e.g. annual, monthly and daily scales). The results of the classification are compared with traditional procedures in a fourth part.

<figure><img src="image_84.png" /></figure>

<figure><img src="image_85.png" /></figure>

<figure><img src="image_86.png" /></figure>

<figure><img src="image_87.png" /></figure>
## 4.1 Available data

The available data for this catchment include:

- 3 h time resolution values of precipitation (rain and snow), potential evapotranspiration, and the maximum and minimum temperatures, distributed on 8×8 km grid for the period 1 August 1981–31 July 1998 and provided by Météo-France. The data were obtained from the SAFRAN (Système d'Analyse Fournissant des Renseignements Atmosphériques à la Neige) analysis (Durand et al., 1993) applied for the whole Rhone basin (including the Saône basin).

- Daily streamflow data from the French hydrographic database (Banque Hydro) for the period 1 August 1981–31 July 1998 at 22 gauging stations. Their catchment areas cover a wide spectrum ranging from 52 to 11 700 km².

- A digital elevation model (DEM) with resolution 200 m, 100 m and 1000 m from the IGN (Institut Géographique National) in France and provided by Water Agency (Agence de l'eau).

- The Corine Land Cover database provided by Institut Français de l'Environnement (IFEN) with a 500 m resolution. The database contains 44 land cover classes organized in three levels (Bossard et al., 2000).

- A soil map from the Soil information system of France from National Institute of Agronomic Research (INRA) with a 1/100 000 resolution (Jamagne et al., 1995). For about one-third of the catchment, another soil database at 1/250 000 resolution was available from the IGCS (Inventaire Gestion et Conservation des Sols) lead by INRA (http://www.gissol.fr/ or http://www.igcs-stb.org/)

- The geology map of France with a 1/100 000 resolution from the BRGM (Bureau de Recherche Géologique et Minières).

5

10

15

20

<figure><img src="image_88.png" /></figure>

<figure><img src="image_89.png" /></figure>

<figure><img src="image_90.png" /></figure>

<figure><img src="image_91.png" /></figure>
- A referential of groundwater systems from SANDRE specification (Secrétariat d'Administration Nationale des Données Relatives à l'Eau) and available at the web site http://sandre.eaufrance.fr/.

In reference with the discussion about data resolution in Sect. 1, we can underline the large heterogeneity in the input maps. Furthermore, meteorological data provide the coarser information, both in space and in time. These data are not used in the definition of the hydro-landscapes units, but the information should be taken into account in the choice of the processes representation, then in the final discretization. The situation for soil data illustrates one of the difficulties when modelling a catchment in France: the sources of data, especially for soils, are not homogeneous.

## 4.2 Determination of the modelling units for a distributed water balance components derivation

In the steps presented in Sect. 3, once the modelling objectives have been defined (step 1) and the data source identified (step 2), the next step is the choice of the processes which will be considered and their representation within the model. Given the time steps and resolutions of the input data, a first step in the modelling is obviously to consider a first discretization at the sub-catchment scale (Sect. 4.2.1) using simple representations such as those presented by Eder et al. (2003). Within a sub-catchment, a large degree of heterogeneity, especially in land-use might remain. Let's suppose that it must be taken into account for an accurate simulation at the chosen scale. It has of course to be proven through simulation, but it is beyond the scope of the paper, where we only wish to illustrate our methodology for landscape discretization, using the available data of the upper-Saône catchment.

### 4.2.1 A first discretization into REWs

The first discretization level was that of the sub-basin (REWs). REWs were determined using DEM analysis and river network structure. Several algorithms were proposed

<figure><img src="image_92.png" /></figure>

<figure><img src="image_93.png" /></figure>

<figure><img src="image_94.png" /></figure>

<figure><img src="image_95.png" /></figure>
<figure><img src="image_96.png" /></figure>

<figure><img src="image_97.png" /></figure>

<figure><img src="image_98.png" /></figure>

<figure><img src="image_99.png" /></figure>

to extract hydrological information from DEM (Peuker and Douglas, 1975; Martz and Garbrecht, 1992). We used the Tardem algorithm (Tarboton, 1997) to derive sub-basin from DEM and a predefined digital river network. This algorithm first performs a detection and treatment of depression zones in the DEM. Then the local direction of out-flows on every cell is calculated using the D8 algorithm. The contributive surfaces are determined for each cell in terms of surface drained through each cell. The sub-basins drained by every link of the network are delimited using a threshold for Strahler order of the river links. Figures 4a and b show a discretization of the Saône catchment in sub-catchment using the first and second Strahler orders. The number of REWs is 341 and 81 for the first and second order respectively, with average surface of 35 and 147 km² respectively.

## 4.2.2 Sub-catchment scale discretization using the landscape classification

### Factors definition

We assumed that the modelling objective was the distributed determination of the components of the water balance at various time scales. One of these components is evapotranspiration and we considered land-use as one of its controlling factors. We assumed the partition of incoming rainfall between runoff and infiltration was controlled by lithology and slope. The use of the topographic index of Beven and Kirkby (1979) was also a possible choice, as well as the direct use of soil-landscape units. In the example described below, we therefore considered three factors: lithology, slope and land-use. The next step was to simplify the information through the definition of classes for continuous data such as slope. In the case of the Saône, we reclassified the original 26 classes of the Corine Land cover map (resolution: 500 m), present in the catchment, into nine classes. This classification was lead by the differentiation in processes representation induced by those land use classes. In the same way we simplified the lithology map according to the parent material and the age of the layers, leading to seven classes for lithology (scale: 1/1 000 000). The slope map was derived from the

5

10

15

20

25
<figure><img src="image_100.png" /></figure>

<figure><img src="image_101.png" /></figure>

<figure><img src="image_102.png" /></figure>

<figure><img src="image_103.png" /></figure>

DEM analysis discussed in Sect. 4.1 and classified into five classes (resolution: 200 m). The various classes for the different factors are given in Table 1. All the maps were re-sampled with a 200 m resolution and were used to define the map of combined factors. The maps of factors are shown in Figs. 5a, b and c. Figure 5d shows the multivariable image of the combined factors obtained after superposition of the three data layers. It is composed of 221 classes of combined factors (amongst the 9×5×7=315 possible classes) and formed the basis for the hydro-landscape mapping procedure. Note in Fig. 5d that within a sub-catchment, the variability of these factors is large, which is one argument for considering this second level of discretization.

10

## Definition of the reference zones

The definition of reference zones may be the most difficult step of this approach. In the ideal case, reference zones should result from a good knowledge on the catchment. Usually, as for the upper-Saône catchment, this is not the case. Therefore, we chose to define them using a statistical and spatial analysis of a simplified combination of factors map. This map was derived using the 9 classes for land-use but only 2 classes for slopes (low and high) and 4 classes for lithology (see Table 1). They were used to define 46 types of reference zones (Table 2). Their location is shown in Fig. 6 and was defined by insuring the representativeness of each reference zone in all areas other the catchment. As much as possible, we chose areas with a certain degree of homogeneity. However, they still keep a histogram of composition for the original 221 combinations of factors.

15

20

## Mapping procedure

For the mapping procedure, we used an available software named CLAPAS developed by Robbez-Masson (1994) (http://sol.ensam.inra.fr/Produits/Asp/ListeProduits.asp?Produit=Clapas). It allows mapping using landscape classification techniques. Several descriptors of the neighbourhood are available: mean, standard deviation,

25
<figure><img src="image_104.png" /></figure>

<figure><img src="image_105.png" /></figure>

<figure><img src="image_106.png" /></figure>

<figure><img src="image_107.png" /></figure>

histogram, matrix of co-occurrence. The software proposes several choices of the
neighbourhood shape: square, ellipse, circle, ring etc. The squared neighbourhood
window is usually used. The choice of the other shapes depends on specific thematic
purposes. We used a square, as particular anisotropy handling was not needed. We
chose to use the histogram of composition as a descriptor of the neighbourhood and
the Manhattan distance to compute the similarity with the composition of the reference
zones. The choice of the Manhattan distance is justified since it is more robust than the
other distances (Kolmogorov, Cramer, etc.) used to perform similarity between vectors
(Robbez-Masson, 1994).

Each point (X) within the catchment was characterized by a specific histogram of the
possible combined factors (p=221) within the neighbourhood window, denoted by Mₓ
= (d1, d2, ..., dp). The histogram calculation was also performed for each k=46 refer-
ences zones. For a reference zone j, the histogram is noted Mⱼ = (m1j, m2j, ..., mpj).
Affectation of a point X to a reference zone j consists in minimizing the Manhattan
distance d(Mₓ, Mⱼ) calculated by Eq. (1); see Fig. 2:

$$
d(M_X, M_j) = \sum_{i=1}^{p} |d_i - m_{ij}| \quad (1)
$$

For the modelling units scale assessment, we varied the neighbourhood windows size,
iteratively from 3 km down to one kilometer.

Resulting hydro-landscapes and distance maps

Figure 7 shows the maps of the hydro-landscapes and their corresponding distance maps for three sizes of the neighbourhood window: 3, 2.2 and 1.4 km. Figure 8 shows the distribution of the surface of the mapped units and Table 3 provides the corresponding statistics.

In the first iteration, we used a neighbourhood windows resolution of 3 km. The resulting map includes homogenous areas with an average surface of 16 km² (standard
<figure><img src="image_108.png" /></figure>

<figure><img src="image_109.png" /></figure>

<figure><img src="image_110.png" /></figure>

<figure><img src="image_111.png" /></figure>

deviation of 99 km²). Fifty percents of the total area was covered with units having a surface lower than 218 km² (top of the Fig. 8). Less than one percent of the total area was covered with modelling unit having a surface lower than 1 km². The corresponding confidence map for this iteration showed larger areas with an unsatisfactory classification (areas in blue on Fig. 7 top right, situated mainly in the Vosges mountainous area). In a second iteration, we tried to refine the mapping by using a smaller neighbourhood window of 2.2 km. The average surface of map units dropped to about 9 km² (standard deviation of 70 km²). Fifty percents of the total area was covered with units having a surface lower than 122 km² (middle of the Fig. 8). Less than one percent of the total area was covered with modelling unit having an area lower than 1 km². The confidence map showed some improvement in the classification (Less areas in blue on the Fig. 7 middle right). In a third step, we tested a window's size of 1.4 km. The map units had an average surface of about 5 km² (standard deviation of 34 km²). Fifty percents of the total area was covered with units having a surface lower than 67 km² (bottom of the Fig. 8). In the latter case, the map units had smaller size than above, three percents of the total area was covered with modelling unit having a surface lower than 1 km². The confidence map showed improvements as compared to the previous iterations.

## Choice of the mapping that fits well with the input data

The use of different sizes for the neighbourhood windows provided modelling units at different scales. The statistical analysis of the surface distribution provided a way to assess the consistency of the results with the scale of the input data. The finest scale of the input data was 1/250 000 and the coarser scale was 1/1 000 000. The final modelling unit could not be more accurate than the finest scale. For a scale of the classified map of 1/250 000, the rule of "quart" states that the modelling units must have a surface larger than 1/4th of 2.5×2.5 km², i.e. 1.6 km². For a scale of the classified map of 1/1 000 000, the units should be larger than 1/4th of 10×10 km², i.e. 6.4 km². The percentage of the cumulated surfaces occupied by the units with surface lower

5

10

15

20

25
<figure><img src="image_112.png" /></figure>

<figure><img src="image_113.png" /></figure>

<figure><img src="image_114.png" /></figure>

<figure><img src="image_115.png" /></figure>

than these two thresholds are provided in Table 3.

A trade off is needed to choose the most appropriate map. The first one (top of Fig. 7) with an average surface of 16 km² seemed too coarse as compared to the scale of the input data. In the third map (bottom of Fig. 7) the average surface of 5 km² is consistent with a scale of 1/250 000 for the classified map. However this feeling of accuracy may be misleading because 18% of the total surface was composed of units with a surface smaller than 6.4 km² and about 6% of the total surface was composed of units with surface smaller than 1.6 km². The better compromise may be to choose the second map (middle, Fig. 7) with an average surface around 9 km². More than 90% of the total surface was covered with landscape units having a surface larger than 6.4 km² and more than 97% of the total surface was covered with landscape units having a surface larger than 1.6 km². Thus, the iterative procedure of this classification approach provided an efficient tool to assess the mapping scale according to the input data scale. Furthermore, the distance map provided an idea of the uncertainties on the modelling units representation, and could be used for uncertainty analysis. However, the input rainfall was available on square grids of 64 km². Therefore, if we consider this data, not taken into account in the classification procedure, the first map resolution could be more consistent with that of the rainfall input data for hydrological modelling.

Other modelling objectives and/or other factors map resolution would produce differ-
ent maps. The presented technique allows a great flexibility in landscape heterogeneity
representation. It may insure a relevant representation of heterogeneity according to
the available data and modelling objectives.

4.2.3 Comparison with the usual mapping techniques for catchment heterogeneity representation

In this section, we present the results of two traditional mapping techniques to repre-
sent the catchment heterogeneity and compare them with those of the classification
methodology presented before.
<figure><img src="image_116.png" /></figure>

<figure><img src="image_117.png" /></figure>

<figure><img src="image_118.png" /></figure>

<figure><img src="image_119.png" /></figure>

Mapping results with the basic smoothing techniques

Smoothing techniques are used to simplify the multivariate factors maps, based on the area of the map units. Usually, the basic technique (we used the ArcView smoothing function) consists in removing areas smaller than those representing the scale of the input data. We used this approach by using areas threshold criteria of 6.4 km² and 1.6 km² corresponding respectively to the two scales chosen for the final landscape map, as defined before. We also performed the mapping by removing the areas lower than 1 km². The corresponding maps are shown in Fig. 8 and the statistics of the surfaces of the mapped units can be found in Table 4.

When areas with surfaces smaller than 6.4 km² were removed, the average surface of the map units was of about 255 km². Less than 0.04% of the total surface was mapped with units having a surface lower than 6.4 km². The mapped units were very coarse and the picture was not very satisfactory. When removing areas lower than 1.6 km², the average surface of the units was of about 32 km² and about 1% of the total surface was covered with units having an area lower than 6.4 km². When the areas lower than 1 km² were removed, the result seemed visually better than the two preceding maps. The average surface was about 14.6 km² and less than 6% of the total surface was covered with map units having an area lower than 6.4 km². The statistics of this map were very close to that of the first map derived from the landscape classification (first line in Table 3).

Mapping by re-sampling and smoothing technique

The final map of hydro-landscapes can be obtained by re-sampling the original multivariate map of the combination of factor in order to get a smoother image. For this re-sampling (we used the resampling function of ArcView), only the desired final factors were considered (in our case the 46 factors retained for the reference zones). Then the same smoothing technique as before was applied to this map. The corresponding maps are shown on Fig. 9 and the statistics of the surface of the units are provided in
Table 4.

When removing areas smaller than 6.4 km², the average surface of the map units was about 64 km² and less than 8% of the total surface was mapped with units having a surface lower than 6.4 km². When removing areas smaller than 1.6 km², the average surface of the map units was about 15 km² and less than 1% of the total surface was mapped with units having a surface lower than 1.6 km².

Comments

In Table 4, we represent the percentages of the 46 initial reference zone encountered in the final map. We observe a large difference between the three mapping technique. The first basic smoothing technique does not allow keeping in the final map a reliable representation of the landscape details. More precisely, the urban zone, wine, agricultural and alluvial areas do not correspond to those present in the input data. When we decreased the minimum area, it was difficult to control the features in the landscape relevant for hydrological modelling, because small units were removed automatically. In any case, with usual mapping techniques, only five to twenty percents of the 46 reference landscapes of Table 2 were represented in the final map (Table 4).

When using re-sampling techniques, the result seemed visually better. But only 28 to 50% of the 46 reference landscapes were represented in the final map.

In the case of landscape classification, between 63 and 90% of the 46 reference landscapes types were represented in the final map (and more when new references were added). This technique provides then better result than the other ones in term of keeping suitable features in the representation of the input data.

4.2.4 Discussion

The usual mapping techniques, because they use a threshold on area, are not efficient in keeping in the landscape suitable features consistent with the available input data and represented hydrological processes. They do not offer any mechanism for map

<figure><img src="image_120.png" /></figure>

<figure><img src="image_121.png" /></figure>

<figure><img src="image_122.png" /></figure>

<figure><img src="image_123.png" /></figure>
<figure><img src="image_124.png" /></figure>

<figure><img src="image_125.png" /></figure>

<figure><img src="image_126.png" /></figure>

<figure><img src="image_127.png" /></figure>

unit scale handling. In any case, there is no quantification of uncertainties of the map units when these methods are used.

The methodology presented in this paper for the derivation of “homogeneous” hydrological modelling units is a first attempt to derive something applicable in practice and allows keeping suitable details in the landscape. It also allows assessing uncertainties on each map units by providing a confidence map. Finally, it allows a large degree of flexibility in the landscape heterogeneity representation, while constraining the hydrologist to formalise the processes he wants to retain and represent within the modelling units.

Nevertheless, we recognise that the method can be improved in many ways for the definition of the classes used in the definition of the combination of factors maps, the reference zones, the choice of the descriptor, the distance used for the classification. A complete sensitivity analysis of the consequences of these choices on the final map should be performed. But this is beyond the scope of the paper where our purpose was to illustrate the methodology outlined in Sect. 3 and initiate discussion about it. This paper does not demonstrate either that the proposed approach is “better” than the traditional use of grid squares or smoothing approach when computing the hydrology. An answer by yes or not is probably not possible because verification data will often be insufficient to discriminate and assess the relevance of the various maps. Nevertheless we think that it can be useful in many ways. By performing sensitivity analysis, it should be possible to get an idea of the refinement below which results – for a given objective and data availability- are insensitive. It can also allow testing hypotheses about catchment functioning or preparing an experimental campaign by highlighting more sensitive zones or defining sampling strategies able to lower model uncertainties. An example of the use of the REC discretization is provided by De Condappa (2005). A catchment of about 50 km² was divided into columns according to the superposition of a land-use map and a soil map. Then 1-D-simulation of water transfer within the various RECS was performed for the derivation of daily water flux below the root zone.

Using this discretization approach we are experimenting new modelling framework
<figure><img src="image_128.png" /></figure>

<figure><img src="image_129.png" /></figure>

<figure><img src="image_130.png" /></figure>

<figure><img src="image_131.png" /></figure>

(Viallet et al., 2006), based on modularity of hydrological processes and appropriate spatio-temporal integration of processes. Once the modelling framework will be built, different tests will provide more precise advice on the suitability of a given catchment discretization for a given objective.

## 5 Concluding remarks

We hope that the ideas expressed in this paper will stimulate research in various domains connected to hydrological modelling. Paraphrasing Beven (2006) about uncertainties, we could consider surface heterogeneity as an opportunity to develop new fields of research instead of considering it a “problem”. As shown in this paper, if we take surface heterogeneity into account, we end up with unstructured meshes, which do not even hold the convexity properties, necessary to use existing classical numerical schemes. There is therefore an opportunity for applied mathematicians to work with hydrologists in order to develop methods applicable on these meshes. One promising example is that of the tRIBS model (Ivanov et al., 2004a) developed using triangular irregular networks, customized to handle hydrological specificities such as catchment boundaries, river network or the requirement to a better discretization along river flat areas where saturated zones can be expected. The inclusion of linear discontinuities (river reaches, hedges, ditches, dikes) is also a promising area of research to derive methods able to i) describe properly these networks, ii) take them into account into hydrological models explicitly at small scales iii) derive simplified parameterisations for use in larger catchment scale models iv) describe the organisation of these networks within the landscape, to be able to take them into account as sub-grid scale parameterisations. For larger catchments, as the Saône catchment presented in this paper, the methodology we proposed is only a first step. We hope that specialists in landscape typology or laws of organization would be interested in taking possession of the subject in order to improve the methodology. The question of spatialization of soil hydraulic properties over the modelling units remains also an open question, deserving research
<figure><img src="image_132.png" /></figure>

<figure><img src="image_133.png" /></figure>

<figure><img src="image_134.png" /></figure>

<figure><img src="image_135.png" /></figure>

investment. A closer collaboration between pedologists and hydrologists, as proposed by Lin et al. (2006) is a necessary step towards the study of such questions.

*Acknowledgements*. The authors thank L. Bargeot and J. Gourmelon from ENESAD in Dijon for fruitful discussion about soil classification and for their assistance in using the Clapas software.

5

References

Abbott, M. B., Bathurst, J. C., Cunge, J. A., O'Conell, P. E., and Rasmusson, J.: An introduction to the European Hydrological System-Système Hydrologique Européen “SHE” 1: History and philosophy of a physically based distributed modeling system, J. Hydrol., 87, 45–59, 1986a.

10

Abbott, M. B., Bathurst, J. C., Cunge, J. A., O'Connell, P. E., and Rasmussen, J.: An introduction to the European Hydrological System - Système Hydrologique Européen, "SHE", 2, Structure of a physically-based, distributed modelling system, J. Hydrol, 87, 61–77, 1986b.

Argent, R. M.: An overview of model integration for environmental applications—components, frameworks and semantics, Environ. Modell. Software, 19(3), 219–234, 2004.

15

Beldring, S.: Multi-criteria validation of a precipitation-runoff model, J. Hydrol., 257, 189–211, 2002.

Beldring, S., Gottschalk, L., Seibert, J., and Tallaksen, L. M.: Distribution of soil moisture and groundwater levels at patch and catchment scales, Agric. Forest Meteorol., 98–99, 305–324, 1999.

20

Bergkamp, G.: A hierarchical view of the interactions of runoff and infiltration with vegetation and microtopography in semiarids shrublands, Catena, 33, 201–220, 1998.

Beven, K.: How far can we go in distributed hydrological modeling, Hydrol. Earth Syst. Sci., 5(1), 1–12, 2001.

25

Beven, K.: Towards an alternative blueprint for a physically-based digitally simulated hydrological modelling system, Hydrol. Processes, 16, 189–206, 2002a.

Beven, K.: Towards a coherent philosophy for modelling the environment, Proc. R. Soc. London A, 458, 2465–2484, 2002b.

Beven, K.: On environmental models of everywhere on the GRID, Hydrol. Processes, 17, 171–174, 2003.

30

Beven, K.: On undermining the science, Hydrol. Processes, 20, 3141–3146, 2006.
<figure><img src="image_136.png" /></figure>

<figure><img src="image_137.png" /></figure>

<figure><img src="image_138.png" /></figure>

<figure><img src="image_139.png" /></figure>

Beven, K. and Binley, A.: The future of distributed models: model calibration and uncertainty prediction, Hydrol. Processes, 6, 279–298, 1992.

Beven, K. and Kirkby, M. J.: A physically-based variable contributing area model of basin hydrology, Hydrol. Sci. Bull., 24, 43–69, 1979.

Blöschl, G. and Sivapalan, M.: Scales issues in hydrological modeling – a review, Hydrol. Processes, 9(3–4), 251–290, 1995.

Bossard, M. Feranec, J., and Otahel, J.: Corine land cover technical guide – Addendum 2000, Technical report No 40, European Environment Agency, Copenhagen, Denmark, 2000.

Branger, F.: Prise en compte des discontinuités hydrologiques dans le transfert de pesticides : utilisation d'une plate-forme de modélisation hydrologique pour la construction de modules spécifiques et l'intégration d'un modèle de bassin versant. Application au bassin versant de la Fontaine du Theil., Université Joseph Fourier, Grenoble I, Grenoble, France, 2007.

Braud, I., Vich, A., Zuluaga, J., Fornero, L., and Pedrani, A.: Vegetation influence on runoff and sediment yield in the Andes region: observation and modelling, J. Hydrol., 254, 124–144, 2001.

Burrough, P. A., van Gaans, P. F. M., and Hootsmans, R.: Continuous classification in soil survey: spatial correlation, confusion and boundaries, Geoderma, 77(2–4), 115–135, 1997.

Carluer, N. and De Marsily, G.: Assessment and modelling of the influence of man-made networks on the hydrology of a small watershed: implications for fast flow components, water quality and landscape management, J. Hydrol., 285(1–4), 76–95, 2004.

Cerdan, O., Le Bissonnais, Y., Govers, G., Lecomte, V., van Oost, K., Couturier, A., King, C., and Dubreuil, N.: Scale effect on runoff from experimental plots to catchments in agricultural areas in Normandy, J. Hydrol., 299, 4–14, 2004.

Ciarapica, L. and Todini, E.: TOPKAPI: a model for the representation of the rainfall-runoff process at different scales, Hydrol. Processes, 16(2), 207–229, 2002.

De Condappa, D.: Etude de l'écoulement d'eau à travers la zone non-saturée des aquifères de socle à l'échelle spatiale du bassin versant. Application à l'évaluation de la recharge au sein du bassin versant de Maheshwaram, Andhra Pradesh, Inde., Université Joseph Fourier, Grenoble, France, 22 Avril 2005, 355 pp, 2005.

Durand, Y., Brun, E., Mérindol, L., Guyomarch, G., Lesaffre, B., and Martin, E.: A meteorological estimation of relevant parameters for snow models, Ann. Glaciol., 18, 65–71, 1993.

Eder, G., Sivapalan, M., and Nachtnebel, H. P.: Modelling water balances in an Alpine catchment through exploitation of emergent properties over changing time scales, Hydrol. Pro-
<figure><img src="image_140.png" /></figure>

<figure><img src="image_141.png" /></figure>

<figure><img src="image_142.png" /></figure>

<figure><img src="image_143.png" /></figure>

cesses, 17, 2125–2149, 2003.

Engeland, K., Braud, I., Gottschalk, L., and Leblois, E.: Multi-objective regional modelling, J. Hydrol., 327(3–4), 339–351, 2006.

Fan, Y. F. and Bras, R. L.: On the concept of a representative elementary area in catchment runoff, Hydrol. Processes, 9, 821–832, 1995.

Favrot, J. C.: Une stratégie d'inventaire cartographique à grande échelle: la méthode des secteurs de reference, Sciences du sol, 27, 351–368, 1989.

Flügel, W. A.: Delineating Hydrological Response Units by Geographical Information System analyses for regional hydrological modeling using PRMS/MMS in the drainage basin of the river Bröl in Germany, Hydrol. Processes, 9(3–4), 423–436, 1995.

Gineste, P.: Télédétection spatiale et modélisation du fonctionnement hydrologique des petits bassins versants. Proceedings of the international workshop, organized by LCT – CEMA-GREF – ENGREF ORSTOM. Water report 16 - FAO, 1997.

Girard, M. C.: Glossaire de science du sol. Institut National Agronomique de Paris-Grignon, 2000.

Gottschalk, L., Beldring, S., Engeland, K., Tallakse, L., Saelthun, N. R., Kolberg, S., and Mo-tovilov, Y.: Regional/macroscale hydrological modeling: a Scandinavian experience, Hydrol. Sci. J. – Journal des Sciences Hydrologiques, 46, 963–982, 2001.

Grayson, R. B., Moore, I. D., and McMahon, T. A.: Physically based hydrologic modelling 1 : A terrain-based model for investigative purposes, Water Resour. Res., 28(10), 2639–2658, 1992a.

Grayson, R. B., Moore, I. D., and McMahon, T. A.: Physically based hydrologic modelling 2. Is the concept realistic, Water Resour. Res., 28(10), 2659–2666, 1992b.

Green, A. E., Astill, M. S., McAneney, K. J., and Nieveen, J. P.: Path-averaged surface fluxes determined from infrared and microwave scintillometers, Agric. Forest Meteorol., 109(3), 233–247, 2001.

Grésillon, J. M.: Contribution à l'étude de la formation des écoulements de crue sur les petits bassins versants. Approches numériques et expérimentales à différentes échelles. Habilita-tion à Diriger des Recherches, Université Joseph Fourier, Grenoble, 1994.

Haverkamp, R., Angulo-Jaramillo, R., Braud, I., Debionne, S., De Condappa, D., Gandola, F., Roessle, S., Ross, P. J., Sander, G., Vachaud, G., Varado, N., Viallet, P., and Zin, I.: POWER - Planner Oriented Watershed modelling system for Environmental Responses-. Final Integrated report AgriBMPwater, Project UE n° EVK1-CT-1999-00025, 26 pp., 2004.

5

10

15

20

25

30
<figure><img src="image_144.png" /></figure>

<figure><img src="image_145.png" /></figure>

<figure><img src="image_146.png" /></figure>

<figure><img src="image_147.png" /></figure>

Ivanov, V. Y., Vivoni, E. R., Bras, R. L., and Entekhabi, D.: Catchment hydrologic response with a fully distributed triangulated irregular network model, Water Resour. Res., 40(11), W11102, doi:10.1029/2004WR003218, 2004a.

Ivanov, V. Y., Vivoni, E. R., Bras, R. L., and Entekhabi, D.: Preserving high-resolution surface and rainfall data in operational-scale basin hydrology: a fully-distributed physically-based approach, J. Hydrol., 298(1–4), 80–111, 2004b.

Jamagne, M., Hardy, R., King, D., and Bornand, M.: La base de données géographiques des sols de France. Etude et Gestion des Sols, 2, 153–172, 1995.

Jothityangkoon, C., Sivapalan, M., and Farmer, D. L.: Process controls of water balance variability in a large semi-arid catchment: downward approach to hydrological model development, J. Hydrol, 254, 174–198, 2001

Klemes, V.: Conceptualisation and scale in hydrology, J. Hydrol., 65, 1–23, 1983.

Krause, P., Kralisch, S., and Flügel, W. A.: Preface: Model integration and development of modular modelling systems, Adv. Geosci., 4, 1–2, 2005, http://www.adv-geosci.net/4/1/2005/.

Lagacherie, P., Cazemier, D. R., van Gaans, P. F., and Burrough, P. A.: Fuzzy k-means clustering of fields in an elementary catchment and extrapolation to a larger area, Geoderma, 77(2–4), 197–216, 1997.

Lagacherie, P., Robbez-Masson, J. M., Nguyen-The, N., and Barthes, J. P.: Mapping of reference area representativity using a mathematical soilscape distance, Geoderma, 101(3–4), 105–118, 2001.

Leavesley, G. H., Markstrom, S. L., Restrepo, P. J., and Viger, R. J.: A modular approach to addressing model design, scale, and parameter estimation issues in distributed hydrological modelling, Hydrol. Processes, 16(2), 173–187, 2002.

Lee, H., Sivapalan, M., and Zehe, E.: Representative Elementary Watershed (REW) approach: a new blueprint for distributed hydrologic modeling at the catchment scale: development of closure relations, in: Predictions ungauged streamflow in the Mackenzie river basin: today's techniques and tomorrow solutions, edited by: Spence, C., Pomeroy, J., and Pietroniro, A., Canadian Water Resources Association (CWRA), Ottawa, Canada, 165–218, 2005.

Lehmann, P., Hinz, C., McGrath, G., Tromp-van Meerveld, H.J. and McDonnell, J.J., 2006. Rainfall threshold for hillslope outflow: an emergent property of flow pathway connectivity. Hydrology and Earth System Sciences Discussion, 3: 2923-2961.

Lin, H., Bouma, J., Pachepsky, Y. A., Western, A. W., Thompson, J., van Genuchten, M. T.,

5

10

15

20

25

30
<figure><img src="image_148.png" /></figure>

<figure><img src="image_149.png" /></figure>

<figure><img src="image_150.png" /></figure>

<figure><img src="image_151.png" /></figure>

Vogel, H.-J., and Lilly, A.: Hydropedology: synergetic integration of pedology and hydrology, Water Resour. Res., 42, W05301, doi:05310.01029/02005WR004085, 2006.

Martz, L. and Garbrecht, J.: “Numerical definition of drainage network and subcatchment areas from digital elevation models.”, Comp. Geosci., 18(6), 747–761, 1992.

Motovilov, Y. G., Gottschalk, L., Engeland, K., and Rodhe, A.: Validation of a distributed hydrological model against spatial observations, Agric. Forest Meteorol., 98–99, 257–277, 1999.

Moussa, R., Voltz, M., and Andrieux, P.: Effects of the spatial organization of agricultural management on the hydrological behaviour of a farmed catchment during flood events, Hydrol. Processes, 16(2), 393–412, 2002.

Mroczkowski, M., Raper, G. P., and Kuczera, G.: The quest for more powerful validation of conceptual catchment models, Water Resour. Res., 33(10), 2325–2335, 1997.

Neitsch, S. L., Arnold, J. G., Kiniry, J. R., and Williams, J. R.: Soil and Water Assessment Tool, SWAT 2000. Theoretical documentation, 542 pp. http://www.brc.tamus.edu/swat/doc.html, 2005.

Parde, M.: Fleuves et rivières, Collection Armand Colin, Paris, 224 pp., 1955.

Peuker, T. and Douglas, D. H.: “Detection of surface-specific points by local parallel processing of discrete terrain elevation data.”, Comput. Graphics Image Process, 4, 375–387, 1975.

Puech, C.: Utilisation de la télédétection et des modèles numériques de terrain pour la connaissance du fonctionnement des hydrosystèmes, Institut National Polytechnique de Grenoble, Grenoble, France, 105 pp, 2002.

Refsgaard, J.C., Henriksen, H. J., Harrar, W. G., Scholten, H., and Kassahun, A.: Quality assurance in model based water management – review of existing practice and outline of new approaches, Environ. Modell. Software, 20(10), 1201–1215, 2005.

Reggiani, P. and Rientjes, T. H. M.: Flux parameterization in the Representative Elementary Watershed approach: application to a natural basin, Water Resour. Res., 41, W04013, doi:04010.01029/02004WR003693, 2005.

Reggiani, P. and Schellekens, J.: Modelling of hydrological responses: the representative elementary watershed approach as an alternative blueprint for watershed modelling, Hydrol. Processes, 17, 3785–3789, 2003.

Reggiani, P., Sivapalan, M., and Hassanizadeh, S. M.: A unifying framework for watershed thermodynamics : balance equations for mass, momentum, energy and entropy, and the second law of thermodynamics, Adv. Water Resour., 22(4), 367–598, 1998.

Reggiani, P., Sivapalan, M., Hassanizadehb, S. M., and Gray, W. G.: A unifying framework for
<figure><img src="image_152.png" /></figure>

<figure><img src="image_153.png" /></figure>

<figure><img src="image_154.png" /></figure>

<figure><img src="image_155.png" /></figure>

watershed thermodynamics: constitutive relationships, Adv. Water Resour., 23(1), 15–39, 1999.

Reggiani, P., Sivapalan, M., and Hassanizadeh, S. M.: Conservation equation governing hill-slope response : exploring the physical basis of water balance, Water Resour. Res., 36(7), 1845–1863, 2000.

Robbez-Masson, J. M.: Reconnaissance et délimitation de motifs d'organisation spatiale. Application à la cartographie des pédopaysages. PhD Thesis, Ecole Nationale Supérieure Agronomique de Montpellier. 191 pp, 1994.

Robbez-Masson, J. M., Foltête, J. C., Cabello, L., et al.: "Prise en compte du contexte spatial dans l'instrumentation de la notion de paysage. Application à une segmentation assistée." Revue internationale de géomatique 9: 173-195, 1999.

Sauquet, E., Krasovskaia, I., and Leblois, E.: Mapping mean monthly runoff pattern using EOF analysis, Hydrol. Earth Syst. Sci., 4, 79–93, 2000, http://www.hydrol-earth-syst-sci.net/4/79/2000/.

Sauquet, E.: Mapping mean annual river discharges: geostatistical development for incorporating river network dependency, J. Hydrol., 331(1–2), 300–314, 2006.

Seguis, L., Cappelaere, B., Peugeot, C., and Vieux, B.: Impact on Sahelian runoff of stochastic and elevation-induced spatial distributions of soil parameters, Hydrol. Processes, 16(2), 313–332, 2002.

Sidle, R. C.: Field observation and process understanding in hydrology: essential components in scaling, Hydrol. Processes, 20, 1439–1445, 2006.

Sivapalan, M.: Prediction in ungauged basins: a grand challenge for theoretical hydrology, Hydrol. Processes, 17, 3163–3170, 2003.

Skoien, J. O. and Bloschl, G.: Catchments as space-time filters – a joint spatio-temporal geostatistical analysis of runoff and precipitation, Hydrol. Earth Syst. Sci., 10, 645–662, 2006, http://www.hydrol-earth-syst-sci.net/10/645/2006/.

Skoien, J. O., Bloschl, G., and Western, A. W.: Characteristic space scales and timescales in hydrology, Water Resour. Res., 39(10), 1034, doi:10.1029/2002WR001736, 2003.

Tarboton, D. G.: A new method for the detemrination of flow directions and contributing areas in Grid Digital Elevation models, Water Resour. Res., 33(2), 309–319, 1997.

Tian, F., Hu, H., Lei, Z., and Sivapalan, M.: Extension of the Representative Elementary Watershed approach for cold regions via explicit treatment of energy related processes, Hydrol. Earth Syst. Sci., 10, 619–644, 2006,

5

10

15

20

25

30
<figure><img src="image_156.png" /></figure>

<figure><img src="image_157.png" /></figure>

<figure><img src="image_158.png" /></figure>

<figure><img src="image_159.png" /></figure>

5

10

15

20

25

30
<table><caption>**Table 1.** Table defining the factors and the delimitation of classes for the whole classification.</caption><tr><td>Land Use</td></tr><tr><td>Urban areas</td></tr><tr><td>Open spaces</td></tr><tr><td>Agricultural areas</td></tr><tr><td>Vineyards and fruit trees</td></tr><tr><td>Pastures</td></tr><tr><td>Sparsely vegetated areas</td></tr><tr><td>Broad leaved-forest</td></tr><tr><td>Coniferous forest</td></tr><tr><td>Water body</td></tr><tr><td>Slope</td></tr><tr><td>Very low slope</td></tr><tr><td>Low slope</td></tr><tr><td>Moderate slope</td></tr><tr><td>High slope</td></tr><tr><td>Very high slope</td></tr><tr><td>Lithology</td></tr><tr><td>Old sediments</td></tr><tr><td>Sediments</td></tr><tr><td>Chalk</td></tr><tr><td>Alluvium</td></tr><tr><td>Sandstone</td></tr><tr><td>Tertiary chalk</td></tr><tr><td>Marl and clayey chalk</td></tr><tr><td>Sand and clay</td></tr></table>

<figure><img src="image_160.png" /></figure>

<figure><img src="image_161.png" /></figure>

<figure><img src="image_162.png" /></figure>

<figure><img src="image_163.png" /></figure>
<table><caption>Table 2. Table defining the reference zones type.</caption><thead><tr><th>Numbers</th><th>Nom de l'Unité paysagère</th></tr></thead><tbody><tr><td>1</td><td>Urban areas</td></tr><tr><td>2</td><td>Open spaces on sedimentary land with low slope</td></tr><tr><td>3</td><td>Open spaces on chalky land with low slope</td></tr><tr><td>4</td><td>Open spaces on alluvium with low slope</td></tr><tr><td>5</td><td>Open spaces on marl with low slope</td></tr><tr><td>6</td><td>Open spaces on chalky land with high slope</td></tr><tr><td>7</td><td>Agricultural areas on sedimentary land with low slope</td></tr><tr><td>8</td><td>Agricultural areas on chalky land with low slope</td></tr><tr><td>9</td><td>Agricultural areas on alluvium with low slope</td></tr><tr><td>10</td><td>Agricultural areas on marl with low slope</td></tr><tr><td>11</td><td>Agricultural areas on sedimentary land with high slope</td></tr><tr><td>12</td><td>Agricultural areas on chalky land with high slope</td></tr><tr><td>13</td><td>Agricultural areas on alluvium with high slope</td></tr><tr><td>14</td><td>Vineyard and fruit tree on sedimentary land with low slope</td></tr><tr><td>15</td><td>Vineyard and fruit tree on chalky land with low slope</td></tr><tr><td>16</td><td>Vineyard and fruit tree on marl with low slope</td></tr><tr><td>17</td><td>Vineyard and fruit tree on sedimentary land with high slope</td></tr><tr><td>18</td><td>Vineyard and fruit tree on chalky land with high slope</td></tr><tr><td>19</td><td>Prairies on sedimentary land with low slope</td></tr><tr><td>20</td><td>Prairies on chalky land with low slope</td></tr><tr><td>21</td><td>Prairies on alluvium with low slope</td></tr><tr><td>22</td><td>Prairies on marl with low slope</td></tr><tr><td>23</td><td>Prairies on sedimentary land with high slope</td></tr><tr><td>24</td><td>Prairies on chalky land with high slope</td></tr><tr><td>25</td><td>Prairies on alluvium with high slope</td></tr><tr><td>26</td><td>Sparately vegetated areas on sedimentary land with low slope</td></tr><tr><td>27</td><td>Sparately vegetated areas on chalky land with low slope</td></tr><tr><td>28</td><td>Sparately vegetated areas on alluvium with low slope</td></tr><tr><td>29</td><td>Sparately vegetated areas on marl with low slope</td></tr><tr><td>30</td><td>Sparately vegetated areas on sedimentary land with high slope</td></tr><tr><td>31</td><td>Sparately vegetated areas on chalky land with high slope</td></tr><tr><td>32</td><td>Broad leaved-forest on sedimentary land with high slope</td></tr><tr><td>33</td><td>Broad leaved-forest on chalky land with high slope</td></tr><tr><td>34</td><td>Broad leaved-forest on alluvium with high slope</td></tr><tr><td>35</td><td>Broad leaved-forest on marl with high slope</td></tr><tr><td>36</td><td>Broad leaved-forest on sedimentary land with high slope</td></tr><tr><td>37</td><td>Broad leaved-forest on chalky land with high slope</td></tr><tr><td>38</td><td>Broad leaved-forest on alluvium with high slope</td></tr><tr><td>39</td><td>Coniferous forest on sedimentary land with low slope</td></tr><tr><td>40</td><td>Coniferous forest on chalky land with low slope</td></tr><tr><td>41</td><td>Coniferous forest on alluvium with low slope</td></tr><tr><td>42</td><td>Coniferous forest on marl with low slope</td></tr><tr><td>43</td><td>Coniferous forest on sedimentary land with high slope</td></tr><tr><td>44</td><td>Coniferous forest on chalky land with high slope</td></tr><tr><td>45</td><td>Forêts de conifères on alluvium with pentes fortes</td></tr><tr><td>46</td><td>Water body</td></tr></tbody></table>

<figure><img src="image_164.png" /></figure>

<figure><img src="image_165.png" /></figure>

<figure><img src="image_166.png" /></figure>

<figure><img src="image_167.png" /></figure>
<table><caption>**Table 3.** Statistics of the surfaces of the mapped units for several sizes of the neighbourhood window.</caption><thead><tr><th></th><th>Average Surface (km²)</th><th>Standard deviation of surface (km²) lower than 1.6 km²</th><th>% of the total surface with units surface lower than 6.4 km²</th><th>% of the total surface with units surface</th></tr></thead><tbody><tr><td>Map 1 (3 km)</td><td>16</td><td>99</td><td>1</td><td>5</td></tr><tr><td>Map 2 (2.2 km)</td><td>9</td><td>70</td><td>2.5</td><td>9.6</td></tr><tr><td>Map 3 (1.4 km)</td><td>5</td><td>34</td><td>6</td><td>18.4</td></tr></tbody></table>

Title Page

Abstract

Introduction

Conclusions

References

Tables

Figures

I ◀

▶I

◀

▶

Back

Close

Full Screen / Esc

Printer-friendly Version

Interactive Discussion
<table><caption>**Table 4.** Synthetic Statistics of different mapping techniques and their capacity to keep in the final map suitable references zones.</caption><thead><tr><th></th><th>Average Area (km²)</th><th>% of the total surface with units surface lower than 1.6 km²</th><th>% of the total surface with units surface lower than 6.4 km²</th><th>% Reference areas kept in the final map</th></tr></thead><tbody><tr><td>Landscape Classification</td><td>9</td><td>2.5</td><td>9.6</td><td>76</td></tr><tr><td>Basic smoothing technique</td><td>14.6</td><td>0.4</td><td>6</td><td>28</td></tr><tr><td>Re-sampling technique</td><td>15</td><td>1</td><td>8</td><td>43</td></tr></tbody></table>

Title Page

Abstract

Introduction

Conclusions

References

Tables

Figures

<figure><img src="image_168.png" /></figure>

<figure><img src="image_169.png" /></figure>

<figure><img src="image_170.png" /></figure>

<figure><img src="image_171.png" /></figure>

Back

Close

Full Screen / Esc

Printer-friendly Version

Interactive Discussion
<figure><img src="image_172.png" /><figcaption>Fig. 1. Illustration of the neighbourhood windows definition.</figcaption></figure>

<figure><img src="image_173.png" /></figure>

<figure><img src="image_174.png" /></figure>

<figure><img src="image_175.png" /></figure>

<figure><img src="image_176.png" /></figure>
<figure><img src="image_177.png" /></figure>

<figure><img src="image_178.png" /></figure>

<figure><img src="image_179.png" /></figure>

<figure><img src="image_180.png" /></figure>

<figure><img src="image_181.png" /><figcaption>**Fig. 2.** Basic principles for the mapping by landscape classification. On the left map is figured a point to be mapped (the black one) with a squared neighbourhood window and the composition histogram within this window. On the right map are shown the reference zones and their composition histograms. In the middle the Manhattan distance is used to search for the minimum distance between the left histogram and those of the reference zones.</figcaption></figure>
<figure><img src="image_182.png" /></figure>

<figure><img src="image_183.png" /></figure>

<figure><img src="image_184.png" /></figure>

<figure><img src="image_185.png" /></figure>

<figure><img src="image_186.png" /><figcaption>Fig. 3. Location of the upper-Saône catchment in France.</figcaption></figure>
<figure><img src="image_187.png" /></figure>

<figure><img src="image_188.png" /></figure>

<figure><img src="image_189.png" /></figure>

<figure><img src="image_190.png" /></figure>

<figure><img src="image_191.png" /><figcaption>**Fig. 4.** Discretization of the catchment in sub-catchments (REW) using the first (Fig. 4a) and second (Fig. 4a) Strahler order.</figcaption></figure>
<figure><img src="image_192.png" /></figure>

<figure><img src="image_193.png" /></figure>

<figure><img src="image_194.png" /></figure>

<figure><img src="image_195.png" /></figure>

<figure><img src="image_196.png" /><figcaption>**Fig. 5.** Maps of (a) land-use, (b) slope and (c) lithology on the upper Saône catchment. (d) Map of the combination of the three previous factors. Superimposed is the discretization of the catchments into sub-catchments at the two Strahler order.</figcaption></figure>
<figure><img src="image_197.png" /></figure>

<figure><img src="image_198.png" /></figure>

<figure><img src="image_199.png" /></figure>

<figure><img src="image_200.png" /></figure>

<figure><img src="image_201.png" /><figcaption>Fig. 6. Map of the reference zones.</figcaption></figure>
<figure><img src="image_202.png" /></figure>

<figure><img src="image_203.png" /></figure>

<figure><img src="image_204.png" /></figure>

<figure><img src="image_205.png" /></figure>

<figure><img src="image_206.png" /><figcaption>**Fig. 7.** Left: hydro-landscapes units and Right: distance map for three sizes of the neighbourhood window; top: 3 km, middle: 2.2 km; bottom: 1.4 km.</figcaption></figure>
<figure><img src="image_207.png" /><figcaption>**Fig. 8.** Fraction of the total surface occupied with units landscapes of surface lower than the surface given in abscissa for three values of the neighbourhood window.</figcaption></figure>

<figure><img src="image_208.png" /></figure>

<figure><img src="image_209.png" /></figure>

<figure><img src="image_210.png" /></figure>

<figure><img src="image_211.png" /></figure>
<figure><img src="image_212.png" /></figure>

<figure><img src="image_213.png" /></figure>

<figure><img src="image_214.png" /></figure>

<figure><img src="image_215.png" /></figure>

<figure><img src="image_216.png" /><figcaption>**Fig. 9.** Mapping of the catchment derived from classical smoothing technique obtained by removing areas lower than (a) 6.4 km², (b) 1.6 km², (c) 1 km².</figcaption></figure>
<figure><img src="image_217.png" /></figure>

<figure><img src="image_218.png" /></figure>

<figure><img src="image_219.png" /></figure>

<figure><img src="image_220.png" /></figure>

<figure><img src="image_221.png" /><figcaption>**Fig. 10.** Mapping of the hydro-landscapes derived by re-sampling into the 46 combined factors retained for the reference zones and removing areas of less than (a) 6.4 km², (b) 1.6 km², (c) 1 km².</figcaption></figure>