1Landscape evolution is the result of various forces: tectonics, erosion, deposition, which in turn, are governed by complex processes related to lithospheric movements, climatic conditions and rock properties. A better understanding of landscape evolution can be assisted by numerical tools that make it possible to explore the complex non-linear and interacting processes shaping landscapes. Significant progress has been made in quantifying the evolution of drainage-basin landscapes by physical processes such as mechanical weathering, sediment transport on hillslopes and channels, and channel incision into bedrock. For many of these processes, we now have geomorphic transport laws that describe their long-term impact on topography (see review by Dietrich et al., 2003). In contrast, few quantitative models exist to describe chemical geomorphic processes and their role in shaping karst terrain. In this paper, we study the evolution of cockpit karst landscapes by means of a simulation model.
2Karst is terrain that possesses distinctive hydrology and landforms as a result of high rock solubility and well-developed secondary porosity (Ford and Williams, 1989). Often, but not always, found in areas underlain by carbonate rocks such as limestone and dolomite, karst terrain appears very distinct to “normal” fluvial landscapes with a lack of established surface drainage networks. Major surface karst landforms include enclosed depressions such as dolines, which are the major routeways by which water is translated from the surface to the subsurface zone. Several conceptual models of karst landscapes stress the importance of the epikarstic zone (i.e. the uppermost zone of weathered bedrock below the soil) which is usually heavily fissured and is the location for most dissolutional activity (Williams, 1983). In the past, a clear climatic control on karst landscapes was recognised, and although a simplistic link between climate and relief is now less accepted, temperate karst areas are often characterised by polygonal networks of dolines whilst humid tropical karsts often consist of networks of dramatic hills. In Jamaica and other areas, conical hills are separated by deep dolines (called cockpits) making up what is known as “cockpit karst” (fig. 1).
Fig. 1 – An oblique aerial view of cockpit karst terrains in Jamaica (Photo by P. Lyew-Ayee).
Fig. 1 – Vue aérienne oblique d’un paysage de karst en cockpits dans la région de « Cockpit Country » en Jamaïque (photo de P. Lyew-Ayee).
3Over the last ten years, much effort has been focused on developing models that can simulate landscape evolution on large time and space scales (e.g. Willgoose et al., 1991; Braun and Sambridge, 1997; Tucker et al., 2001). Most such models incorporate both short-range processes such as hill slope diffusion, and long-range processes such as fluvial transport. A partial list of current models includes SIBERIA (Willgoose et al., 1991), GOLEM (Tucker and Slingerland, 1994), DELIM (Howard, 1994), CASCADE (Braun et Sambridge, 1997), CAESAR (Coulthard et al., 1998), ZSCAPE (Densmore et al., 1998), and CHILD (Tucker et al., 1999) and they are reviewed in Coulthard (2001).
4Although modelling of subsurface runoff and relief generation in karst has resulted in considerable research (Groves and Howard, 1994; Clemens et al., 1997; Kaufman and Braun, 1999, 2000), very few models explore rock dissolution processes by surface water and surficial landscape evolution (Ahnert and Williams, 1997; Kaufman and Braun, 2001).
5Starting from a squared two-dimensional network and simple erosion laws, Ahnert and Williams (1997) recreated tower karst landscapes. Kaufman and Braun (2001) demonstrated how carbon dioxide-enriched surface runoff shapes a landscape typical of large karstic valleys by the chemical dissolution of calcite contained in rocks. They used the CASCADE model for their study, and included an erosion factor related to the dissolution of carbonate rocks. Fleurant et al. (2007) used the CHILD model to simulate the evolution on a single cockpit.
6A wide variety of karst landscapes exist whose morphologies depend on a number of variables (White, 1984): chemical variables such as temperature, pressure of carbon dioxide and rainfall; physical variables such as relief and rainfall; geological variables such as tectonics, characteristic features of soluble rocks, stratigraphy and lithology. The variations of the different components listed above over a range of timescales result in different landscapes: doline karst, fluvio-karst or tower karst. Nevertheless, a common characteristic of karst landscapes is that they are mainly composed of a series of closed depressions and hills (Trudgill, 1985).
7Cockpit karst is a form of polygonal karst usually thought to occur only under humid tropical climatic conditions. Polygonal karst (Ford and Williams, 1989) is characterised by a network of closed depressions covering an entire landscape; cockpit karst is a variant of this where the closed depressions are deep, large and often star-shaped with residual hills found between them. In order for cockpit karst to form several requirements must be met, i.e. abundant runoff, massive and pure limestones, and a regional system of fractures which allow the development of underground drainage (White, 1984). Sweeting (1972) hypothesizes that annual rainfall of over 1500 mm is required to form cockpit karst. The type site for cockpit karst is found in the Cockpit Country in Jamaica, where the landscape is underlain by White Limestone, a massive and exceptionally pure unit (Lyew-Ayee, 2004). Karstic processes are characterised by the predominance of mass transport of solutes over other transport processes (White, 1984). This is why the integration of carbonate rock dissolution processes is essential for the modelling karstic phenomena.
8The objective of the present paper is to examine the topographic, hydrologic, and geochemical conditions under which cockpit karst can develop. We describe the integration of chemical dissolution processes of carbonate rocks in the computer model CHILD. Integrating chemical erosion into the model will allow us to test some hypotheses of cockpit karst landscape development in Jamaica, and will also permit further testing of karst landscape evolution theories.
9We will consider briefly the functioning principles of the CHILD model as well as the modifications which were necessary to simulate rock dissolution processes. These modifications build upon the work presented in Fleurant et al. (2007) and are now tested at the regional scale. The simulated landscape will be compared to a real cockpit karst landscape through three morphometric criteria typical of this type of landscape, and recently highlighted by Lyew-Ayee (2004): average slope, relative relief and scaling characteristics. Such a comparison will make it possible to provide both ranges of validity for the various parameters of dissolution models and put forward hypotheses of scenarios for cockpit formation.
10The CHILD model (Channel-Hillslope Integrated Landscape Development) is a numerical model of landscape evolution (Tucker et al., 1999). The evolution of topography over time is simulated through the interaction and feedback between surface flows, rock erosion and the transport of sediments. CHILD is a numerical simulation method used for modelling a large number of processes related to the geomorphology of catchment basins. As it has been well-described in previous work, we will consequently only present a small number of these features. In this computer model, the evolution of topography over time results from the combination of several factors:
11These various factors are numerically solved using a particular spatial framework: Voronoï polygons. The first three factors are present in most models of landscape evolution; a detailed description can be found in Tucker et al. (2001). The last factor, concerning rock dissolution, will be discussed in a subsequent section. Here we focus on dissolution processes of limestone, the other processes are neglected. Several studies (e.g. White, 1984 or Lyew-Ayee, 2004) and the presence of thick equatorial forest strengthen the assumption of the dissolution erosion predominance over mechanical erosion. So equation (1) reduces to:
12Implications of uplift and mechanical erosion are discussed in several papers (e.g. Pfeffer, 1997). Cockpit karst landforms in Jamaica are restricted to areas of the White Limestone Group. This white limestone was first exposed to erosion by uplift in the Upper Miocene and no traces can be discerned today of a surface drainage system in the cockpit areas. Thus, we assume that landscape development from the Miocene onwards has been dominated by solutional processes, without any significant contribution from either other denudation processes or tectonics.
13The spatial framework of the CHILD model is composed of an irregular network of points and is used to obtain a reliable discretization of the topographic surface. These points, also called nodes, are connected to form a mesh of triangles known as a Triangulated Irregular Network, or TIN (Braun and Sambridge, 1997; Tucker et al., 2001). This gridded model is built up using Delaunay triangulation and offers several advantages over a raster-based approach (Braun and Sambridge, 1997; Tucker et al., 2001). From the numerical point of view, triangulation makes it possible to solve equations involving partial derivatives using the finite-volume method (Versteeg and Malalasekera, 1995). Each node of the triangulation is associated with a Voronoï cell (fig. 2). Such a method can thus be used to solve physical equations related to geomorphologic processes, with exchanges of mass between nodes calculated at the Voronoi polygon interfaces.
Fig. 2 – Schematic illustration of the spatial framework in CHILD model. The study area (16 km2) is discretized into 40000 Voronoï cells. Each cell is 400 m2 large.
Fig. 2 – Illustration schématique du mode de discrétisation du modèle CHILD. La zone d’étude a une surface de 16 km2 et est découpée en 40000 cellules de Voronoï. Chaque cellule de Voronoï représente une surface de 400 m2.
14The temporal framework of the CHILD model makes it possible to consider the significant time disparities resulting from geomorphologic processes: the evolution of topography spreads over thousands of years, whereas rainfall frequency can range from a few minutes to a few days. Rainfall variability at a geological time scale is simulated using a stochastic method based on a Poisson storm-arrival model and exponentially distributed storm intensity and duration, which gives the distribution of rainfall frequency and intensity (Tucker and Bras, 2000). Equations are solved iteratively, according to an alternation of rainy and dry periods. Storm intensity and duration are assumed to be independent. The probability density functions for storm intensity, duration and interstorm interval are given by:
15Surface runoff is modelled using an algorithm that routes water from one cell to the next following the path of steepest descent. In the case of an enclosed depression, water either evaporates or is forced to find an outlet, using a lake-filling algorithm (Tucker et al., 2001). Local runoff generation, at one node, can be modelled using several functions: Horton runoff or excess runoff on saturated soils.
16Although calcite is slightly soluble in deionised water, most karst denudation occurs through the dissolution of calcite by water acidified by carbon dioxide. In the case of open systems, carbonates are affected by climatic events (rain in the studied case) and other flowing waters (rivers for example). The dissolution of carbonate rocks will then depend on the flux of water in contact with carbonates, and on the hydrodynamic nature of flows (Dreybrodt, 1988).
17The calcite dissolution process in an open system can then be expressed as follows (Trudgill, 1985):
18CaCO3 + CO2 + H2O CaCO3 + H2CO3 + Ca2+ +2HCO3-.
19It should be noted that, for stoichiometric reasons, one mole of Ca2+ consumes one mole of CO2. Consequently, the absolute flow value of CO2 is equal to that of Ca2+.
20The karst denudation process corresponds to the erosion of carbonates through dissolution. Several models have been put forward to quantify this erosion (Corbel, 1959; Ford, 1981). White's (1984) denudation model (DR) is the most comprehensive and widely used; it speculates that infiltration waters and carbonates are in equilibrium:
21where K1 (mol/m3), K2 (mol/m3), KC (mol2/m6) and KH (mol/m3/atm) are equilibrium constants which only depend on temperature, gX is the activity coefficient of X and PCO2 (atm.) is the pressure of carbon dioxide in open system, i.e. under atmospheric conditions. The DR values may vary according to the temperature and CO2 pressure. The values of PCO2 pressure depend on local factors: the atmospheric conditions and soil biological activity. The relation between PCO2 and the actual evapotranspiration rate (E) and therefore the geographical variation of carbon dioxide pressure in soils is (Brook and Hanson, 1991):
22In Jamaica for example, E ≈ 1 m/y, which gives PCO2 ≈ 0.0176 atm. Such values are consistent with the values Smith et al. (1976) put forward for tropical areas: PCO2 ≈ 0.01 atm.
23The relation (5) combines rock, equilibria and climatic factors. Each of the three climatic variables (temperature, carbon dioxide pressure and effective rainfall) has a different impact on denudation rate. In order of importance, the key controls are: effective rainfall (linear relation), carbon dioxide pressure (cube root relation) and temperature, which has a non-linear influence on chemical constants. As mentioned previously, White's model (1984) gives the maximum value of karst denudation for the whole studied area. Therefore, if the model (5) is applied, it is assumed that the karst denudation rate is uniform across the studied karst. For typical values of cockpit karst terrains DR ≈130 mm/ky (Smith et al., 1976).
24As the karstification process is complex it is difficult to establish clear and quantifiable relations between its various components: calcite dissolution, geologic properties of carbonates, rock structure, and climatic variations. There is a considerable literature on the qualitative description of these processes (e.g. Jennings, 1985; Sweeting, 1972; Williams, 1985), but little modelling work has so far been produced quantifying these. Although dissolution plays an essential role in the development of karst, it is not the only process involved with collapse, and subsidence brings material down into closed depressions (Jennings, 1985). Many authors also highlight the close relation between surface denudation and underground karstification processes, but do not provide quantitative data to link the two (Dreybrodt, 1988).
25These linkages are strongly suggested by the epikarst zone (e.g. Ford and Williams, 1989). The epikarst is the uppermost zone of highly weathered karst bedrock below the soil. Runoff waters passing through the soil become enriched in carbon dioxide, producing a maximum dissolution rate within the epikarst zone. Most dissolution is accomplished within 10 m of the surface as a result. This dissolution within the epikarst zone leads directly to surface lowering, as soils subside into the enlarged voids created. Thus, the higher the subsurface dissolution rate, the more rapidly topography lowers over time. Therefore, in the model presented in this paper, the rate of dissolution of topography is taken to be proportional to the underground dissolution.
26Equations of the model express a simple linear relationship between surface denudation rate and underground dissolution, thus subsurface flow:
27where zi (m) is the altitude of the node i, DRmax (m/s) and DRmin (m/s) are respectively the maximum and the minimum denudation rates on the study area, qi (m3/s) is the subsurface flow going through the Voronoï cell i, qmax (m3/s) and qmin (m3/s) are respectively the maximum and the minimum subsurface flow on the study area at time t (s). Equation (7) is illustrated by figure 3: the more subsurface flow increases, the more denudation rate increases.
Fig. 3 – Principle of equation (7). Denudation rate of the topography DR is taken to be proportional to the underground flow q (subcutaneous zone).
Fig. 3– Illustration du fonctionnement de l’équation (7). L’érosion karstique DR est directement en relation avec le flux de subsurface q (zone épikarstique).
28It is assumed that development of cockpit landforms begins on a smooth gently undulating surface (X = Y = 4000 m). The average altitude of these depressions and summits is Z0 = 600 m (fig. 4).
Fig. 4 – Boundary conditions. Initial topography is a smooth 600 m high surface with random sinks (-10 cm). The density of these depressions is 25 km-2.
Fig. 4 – Conditions aux limites. La topographie initiale est une surface lisse de 600 m d’altitude sur laquelle sont réparties aléatoirement des dépressions (-10 cm). La densité de ces dépressions est d’une pour 25 km2.
29In fact, following the limestone sediment deposition that occurred during the Oligocene and Eocene, the cockpit area of Jamaica is thought to have undergone a tectonic elevation of the ground, from subaqueous to approximately Z0 = 600 m above sea level (Lyew-Ayee, 2004). Consequently, the karstification process is presumed to have only started during the Pliocene, when climatic conditions were favourable (Pfeffer, 1997).
30The limestone plateau is discretized into 40000 Voronoï cells. Hydraulic conductivity of each Voronoï cell is 10-7 m/s. This value is typical of slow initial karstification processes (Singhal and Gupta, 1999). This hydraulic conductivity controls the partitioning of surface runoff and subsurface flow. Initially, depressions, cockpits or sinks are created using a random uniform function in the study area (16 km2). The density of sinks is approximatively d = 25 km-2 (fig. 5).
Fig. 5 – Example of a simulated cockpit karst landscape. Here, the climatic factor is 2 and the simulation time is 4 My. The DEM of a simulated landscape can be used to perform analysis. In this study, three analyses are carried out: average slope, relative relief and scaling characteristics.
Fig. 5 – Exemple d’un paysage de karst en cockpit simulé. Pour cette simulation, le facteur climatique est de 2 et le temps de simulation est de quatre millions d’années. Le MNT simulé, peut ainsi être analysé : calculs de la pente moyenne, du relief relatif moyen, des propriétés d’échelle.
31This average value was calculated by Lyew-Ayee (2004) for the whole cockpit country in Jamaica. In the model, the value of rainfall is <P> = 2 m/y. Many field studies (e.g. Sweeting, 1972) clearly show that average rainfall has to be greater than 1 m/yr to produce cockpit karst landscapes. Moreover, actual average effective rainfall in Jamaica is about 2 m/yr (Lyew-Ayee, 2004).
32The local maximum and minimum denudation rates were estimated by Smith et al. (1972) in the Cockpit karst areas:
33where (CaCO3) (ppm) is the minimum or maximum hardness values of CaCO3 which were founded in the cockpit karst areas and rCaCO3 (kg/m3) is volumetric weight of calcite. For values of the parameters given in the table 1, one obtains DRmin = 48 mm/ky and DRmax = 240 mm/ky. These denudation rates are local values of the limestone dissolution. DRmin occurs when subsurface flow is low whilst DRmax holds at higher values of subsurface flow. The average value of these denudation rates is 144 mm/ky which is very close to the typical value uses for cockpit karst terrains (DR = 130 mm/ky from Smith et al., 1976).
Table 1 – Values of the model's parameters.
Tableau 1 – Description et valeur des paramètres utilisés pour les simulations du modèle.
34The objective of simulations is to test the mathematical model (7), in order to assess its capacity to recreate geomorphologic processes typical of cockpit karst landscapes of the cockpit type.
35There has been considerable research effort expended on morphometric analysis of karst landscapes aiming at producing quantitative data with which to test landscape evolution models (e.g. Brook and Hanson, 1991). Such research has been vastly aided by the advent of digital datasets from which DEMs can be produced. Lyew-Ayee (2004) used a suite of morphometric techniques to characterise and analyse and area covering 84 km2 of the cockpit karst in Jamaica. He proposed several parameters, which could be used to distinguish cockpit karst from other relief types, of which we use the following: average slope, relative relief and scaling properties. Slopes are calculated from the DEM using a moving window – a 33 cell kernel – which is displaced throughout the grid. To calculate slopes, a third-order finite difference estimator is used that is based on the eight outer points of the moving window. Determination of relative relief needs to identify sinks and summits within the landscape. Then relative relief may be determined using sink and summit points. A cluster analysis of both the average slope and the average relative relief, that is the distance between depressions and summits of cockpits, was used (Lyew-Ayee, 2004) as a crucial measure in distinguishing cockpit and non-cockpit karst landscapes. The scaling properties give a picture of the scale of variations of the horizontal and vertical characteristics of the landscape (Weissel et al., 1994). These scaling properties can be studied by conducting neighbourhood analysis of the grid. The range of elevation within a specified horizontal radius (from 10 m to 4000 m) is calculated. Results of the scaling characteristics are plotted as a log-log plot. Data show a major scaling break in the plot of vertical relief versus horizontal scale at around 100 m in the X axis. This value reveal the average spacing between sinks and summits which is significant with those observe in initial geometry conditions of the model (fig. 6).
Fig. 6 – Evolution of average slope versus relative relief.
Fig. 6 – Évolution de la pente moyenne en fonction du relief relatif moyen du paysage simulé.
The simulation 1 is computed with
La simulation 1 correspond à un facteur climatique de
36Then, data clearly show two scales to the landscape: before and after the scaling break. Data can thus be fitted by a linear equation before and after this breaking scale: z = c(x,y)m where c is an amplitude factor and m is a scaling exponent (tab. 2).
Table 2 – Morphometric properties of the Cockpit Country study area (Lyew-Ayee, 2004).
Tableau 2 – Propriétés morphométriques des trois zones d’étude dans la région des « Cockpit Country » (Lyew-Ayee, 2004).
Name of the ‘Cockpit Country’ area
|
Slope (in°)
|
Relatif relief (in m)
|
c (before x = 100 m)
|
c (after x = 100 m)
|
m (before x = 100 m)
|
m (after x = 100 m)
|
Barbecue bottom
|
31
|
74
|
0.08
|
11.21
|
1.45
|
0.44
|
Quickstep
|
24
|
54
|
0.099
|
9.01
|
1.33
|
0.45
|
Windsor
|
26
|
57
|
0.097
|
9.77
|
1.32
|
0.41
|
37Starting from these boundary and initial conditions, and applying the model (7), is it possible to simulate the landforms of a cockpit karst at regional scale?
38The model needs the input of many variable values (tab. 1) most of which are drawn from published sources (e.g. Smith et al, 1972 or Lyew-Ayee, 2004). The only parameters which can be modified are the storm duration Tr and the interstorm period Tb. The duration ratio
39Simulations were run over a time period of 10 million years, based on published suggestions that karstification has been active here since the Pliocene (Pfeffer, 1997). The average hillslope and relative relief of the test-cockpit are computed every 10,000 years in order to note its evolution over time. Therefore, it is possible to compare the simulated values with the real average measurements of the test-cockpit in the studied area of Jamaica. The simulated landscapes (figure 5 is an example for have been analysed in the same way as a DEM of a real cockpit karst landscape, with data collected for average slope, relative relief and scaling properties.
40The results of the model simulations are shown in figures 6 to 9. Firstly, we can observe that the relative relief, the average slope and the scaling characteristics of the simulated landscape are quite similar to those observed in the field. These simulations lead to a key conclusion: using the chosen values of the model parameters, it is possible to recreate the morphometry of a cockpit karst landscape.
Fig. 7 – Evolution of average slope versus relative relief.
Fig. 7 – Évolution de la pente moyenne en fonction du relief relatif moyen du paysage simulé.
The simulation 2 is computed with
La simulation 2 correspond à un facteur climatique de
Fig. 8 – Evolution of average slope versus relative relief.
Fig. 8 – Évolution de la pente moyenne en fonction du relief relatif moyen du paysage simulé.
The simulation 3 is computed with
La simulation 3 correspond à un facteur climatique de
ttt
41This result is very consistent with many authors' observations (e.g. Sweeting, 1972) that humid tropical climates, characterised by year-round abundant rainfall, are needed to produce a cockpit karst landscape.
42Results of the scaling analysis are presented on figure 9.
Fig. 9 – Scaling properties of the simulated landscape.
Fig. 9 – Les propriétés d’échelle des paysages simulés.
These virtual landscapes are simulated for several climatic factors:
Ces paysages virtuels sont simulés pour différentes valeurs du facteur climatique :
43For each simulation, only three scaling curves are plotted corresponding to the most relevant times. Field data from the Jamaican cockpit karst clearly show a scaling break at around 100 m (average distance between sinks and summits), whereas these two scales (before and after 100 m) are less distinct within the simulated landscapes. The reason is probably due to the standard deviation of the average distance between sinks and summits. In the simulations, this distance is about 109 m, and the standard deviation is large (40 m), so that most of the distances range between 40 m and 180 m when the range of the distance in the field is from 80 m to 120 m only (Lyew-Ayee, 2004).
44Concerning the scale less than 100 m, one can notice that simulations begin only at 20 m on the horizontal scale which is the grid resolution of the model. The scaling characteristics of the simulated landscapes, before and after the scaling break, show a relatively good general trend. But, if we go into the details one can notice (fig. 10) that relative (%) and absolute (m) errors are lower after the scaling break. This is probably due to the size of the model grid: the size of the Voronoï cells is, in average, 20 m. As the average radius of a cockpit is 109 m, we only have 5 Voronoï cells to model the hillslope landform of a cockpit. Beyond the horizontal scale (X,Y) = 100 m, the landscape is analysed at a regional scale, so the grid size of the model (20 m) becomes less important.
Fig. 10 – Model errors. Example of relative (black curve in %) and absolute (grey curve in m) errors between the "best" simulated curve (Simulation 1: t = 1.2 My) and data (Windsor).
Fig. 10 – Erreurs du modèle. Exemple des erreurs relatives (courbe noire en %) et absolues (courbe grise en m) entre la meilleure des simulations (simulation 1 pour le temps t = 1.2 millions d’années) et les données de terrain (Windsor).
45The behaviour of the scaling characteristics curves versus the climatic factor is the same as that of the slope/relative relief index. Increasing the climatic factor, increases the time needed to obtain a landscape with similar landforms.
46Simulations carried out in the model are compared to real cockpit karst landforms from the Cockpit Karst Country in Jamaica. Morphometric characteristics can be deduced from simulated karst and compared with those measured from DEMs of real karst landscapes. Three morphologic characteristics based on Lyew-Ayee' work (2004) were chosen to compare simulated and real karst: average slope, relative relief and scaling characteristics.
47The similarity between the simulated and real karst morphology constitutes an essential stage in the validation of the hypotheses put forward on their genesis. Though it is a necessary stage, it is probably not a sufficient one, and further tests need to be carried out to test the validity of our model. The present study highlights how essential it is to take spatial anisotropy of denudation into account when studying the morphometric evolution of cockpit karst.
48The spatial anisotropy of denudation corresponds to the heterogeneity of many variables, e.g. hydraulic conductivity or depth of soil, which are frequently observed in karstic systems and that we have encapsulated in a simple linear model. Further work could usefully be done in the field in quantifying different aspects of these spatial variations, and using them to refine our model if necessary.
49Considering that underground dissolution has been shown to be a significant factor in the evolution of surface karst geomorphology, as highlighted in the epikarst, the simulations obtained using our model are encouraging. The spatial-anisotropic dissolution model reveals that the lowest topography locations on the hillslope receive the largest quantity of infiltrating water, and vice versa. Consequently, dissolution of the topography is more active downslope.
50The aim of the present study on the morphology of cockpit karst was to simulate their evolution processes on the basis of various hypotheses on their genesis. The model was based on a spatial anisotropy of dissolution processes which has already been established as being necessary to simulate karst landscape development at the scale of individual cockpits by Fleurant et al. (2007). Denudation rate of the topography was linked with underground dissolution through the concept of the epikarst zone.
51Based on the epikarst processes, and considering that the denudation rate of the surface topography is proportional to underground dissolution, we have been able to develop a cockpit karst landscape evolution model at regional scale. The model is based on a number of simplifying assumptions and could be improved in several areas if more field data become available. To compare real and simulated landscapes we have used three types of morphometric parameters: average slope, relative relief and scaling characteristics (Lyew-Ayee, 2004). However, several other parameters could also be used to describe and compare modelled output with cockpit karst terrains and particularly geostatistic analysis or aspect analysis.