Geomorphological studies have helped researchers understand topographical relief. According to [1] , the main objective of geomorphology is to analyze various forms of relief to understand the relationship between past and present processes. Understanding of local relief and the processes that shape, it is important when seeking to responsibly utilize natural resources with minimal damage to the environment. According to [2] , a river system results from the sculpting of active processes that result in a change in relief. Such systems should be studied and analyzed because their formation is “the most active morphogenetic process in sculpting the Earth’s landscape”.

It is important to analyze and understand a watershed’s morphology not only because of its role in the area’s geomorphology but also because watersheds are treated as major geographical units for planning and geo-envir- onmental analysis in Brazil. Morphometric analysis can help determine the capabilities and constraints of a watershed and its suitability for various anthropogenic uses.

In this study, fitness is defined as the capacity of relief to support human activities while minimizing the impact on the watershed. Thus, studying the geomorphological characteristics of a watershed has implications for its occupation and the use and management of its natural resources. The geomorphological characteristics of watersheds are instrumental for academic research and land use planning.

The hydrological behavior of a watershed [3] is determined by its morphological characteristics. To understand the relationships between the hydrological processes in a basin, their features must be expressed quantitatively. Quantitative analysis of morphometric parameters can be useful in assessing watersheds for proper management, conserving natural resources, and supporting the most appropriate uses of their land and water. According to [4] , morphometric analysis provides a quantitative description of a drainage system, which is one of the most important aspects when characterizing a watershed.

[5] finds that the organization of a drainage system largely reflects the complexity of its development and evolution and results from interactions between the local climate, lithology and tectonics. Understanding this organization is important for environmental management and especially the “planned land use, for which the influence of anthropogenic changes in land use and climate on the watershed must be predicted” [6] . Following drainage morphometry studies by [4] [7] - [11] , many studies have used this method to characterize river basins and sub-basins in different geographical regions.

In this sense, morphometric analysis can be considered an essential tool for the qualification of a watershed. The proper interpretation of a drainage network, including the local geology and morphology, can reveal numerous issues related to erosion, which is the principal element of physiographic and geomorphological processes [2] .

Morphometric analysis assumes that quantitative data can serve as parameters for the physical qualification of the basin and its hydrosedimentological processes, characterizing them precisely and revealing their homogeneity, according to [12] . Thus, variables related to the length, height, width, volume, gradient, density, and frequency are used for vulnerability assessments from an environmental geomorphology perspective.

Given the current demand for water resources and growing concern about environmental problems at the local, regional and global scales, morphometric studies can play an important role in consistently characterizing watersheds. Thus, they can support integrated management by quantitatively describing the complex landscapes contained in a drainage basin.

Relief features have cartographic elements that are present in various sizes and shapes. Characterizing the sizes and shapes of these features is important in understanding their genesis and current dynamics, and their cartographic representation allows synthesizing their roles in a landscape. Therefore, geomorphological mapping and classification has been the subject of research and experimentation since the 1960s in various countries, particularly in Europe [13] .

Different land uses (rural, agricultural, industrial and urban) in plains or on hillslopes are elements in a basin’s landscape and scenery. They can have various negative environmental impacts, including vegetation removal, pollution, and silting. Human activities occur at a higher intensity in urban areas because of the concentration and variety of economic activities in such locations that result in the degradation of natural resources [14] .

Geomorphological mapping should consider both taxonomy and scale and aim to characterize various types of relief, grouping them into different classes. In general, maps can be presented as morphostructural domains, geomorphologic regions, geomorphologic units and model types [15] .

When studying relief based on geomorphological units, model types and forms, it should be noted that other elements of the geomorphological system might require special attention. In particular, hydrographic information on land slope and morphometric drainage is important, including relief dissection parameters, such as drainage density, split ratios and the average length of channels. This set of morphometric parameters can be grouped into classes based on the different patterns of the existing relief in a certain geomorphological region.

A morphometric analysis of relief can also help subdivide it into different geomorphological units [16] . Conducted a survey in which the morphometric analysis of the Passa Cinco (SP) river basin allowed the identification of five major geomorphological units. Their methodology was based on the previous studies of slope [17] ; [18] , horizontal dissection [19] ; [20] , vertical dissection, relief [21] , and density and depth of dissection [22] . Their results showed that morphometric analysis can also facilitate understanding the morphological structure of relief and allow for more extensive descriptions of the shapes, ages and genesis of relief.

Therefore, geomorphological characterization is a fundamental basis for environmental studies. This study aimed to analyze morphometric variables in the Maranhão River Basin (DF-GO), Brazil and indicate its geomorphological suitability for various anthropogenic uses.

The mapping studies have used mapping studies of individualized form without the adequate incorporation of morphometric elements with the proposals of public policies for zoning of the territory. So in this text there is the proposal to analyze the potential of morphometric information for methodological advances in land planning policies.

The Maranhão River Basin is an important unit of study because it covers two states (Goiás and Tocantins) and the Distrito Federal in mid western Brazil. In addition, it is an important catchment area of the Araguaia- Tocantins water system (Figure 1). The river’s source is in Brazil’s Central Plateau and Brasília City. A significant portion of the basin’s headwaters is located in the Environmental Protection Area (APA) of the Cafuringa and Emendadas Water Ecological Station [23] . In January 2002, another portion of the Maranhão River Basin was added to the Central Plateau APA (Federal Decree No. 9468 on 10 January 2002).

This study used drainage and topography data from 1:250,000 maps to conduct a basic morphometric analysis of the basin. All calculations to identify the morphometric indicators followed the recommendations and equations proposed by [2] [4] [8] [24] , which are widely accepted in the literature. The quantities used in this work were the bifurcation ratio (R_b), drainage density (D_d), river density (D_r), extension of the path surface (E_ps), basin area (A), basin length (L), ratio of the average length between two subsequent orders (MLR), and channel tortuosity index (I_s).

The split ratio (LR) is the ratio of the total number of channels in a certain order and the total number of channels in the order immediately above. This ratio must be constant and can never be less than 2 [25] . It is

expressed by the following equation:

where: N is the total number of channels in a particular order; and N_u+1 is the total number of immediately higher order channels.

The drainage density (D_d) describes the total length of channels within a catchment area (in square kilometers). According [2] , this parameter represents the hydrological behavior determined by lithology and geological structures, which control the infiltration capacity and behavior of surface channels. The drainage density has an inverse relationship to the river density; thus, when more channels exist, they are less extensive. This quantity is expressed by the following equation:

where: L_t is the total length of the channels (km), and A is the total basin area (sq/km).

The river density (D_r) reflects the relationship between the number of streams and the area of a basin. This parameter is related to a basin’s ability to generate new fluids. This quantity is expressed by the following equation:

where:

N represents the total number of rivers or channels, and A is the area of the basin (km²).

The extension of the surface route (E_ps) parameter is the average distance traveled by a flood between the interfluve and the permanent river. This variable is important because it represents the link between the hydrological and physiographic development of a drainage network. This quantity is expressed by the following equation:

where:

E_ps is the extent of the surface route, and D_d is the drainage density.

The basin area (A) represents the entire area drained by a river system and is generally expressed in square kilometers. This ratio can be measured using conventional methods or more sophisticated methods using specific software.

The length of the entire channel (L) is distinguished from the length of each order of channels by the letter u. The sum of the lengths of each order of channels is L_u, and L_t represents the total length of all water bodies in a river basin.

The sinuosity index of channels (I_s) relates the true length of the channels with the straight length between their extreme points. This parameter includes the influence of sediment load and the lithological and structural subdivisions [26] . This quantity is expressed by the following equation:

where:

L is the channel length (km);

D_v is the vector distance (km) between the extreme points of the same channel.

Values close to 1 indicate high structural control (high energy), and values above 2 indicate low power; intermediate values generally characterize transitional forms between straight and meandering channels.

In this study, we selected slope classes based on the work of [27] , as shown in Table 1. This table summarizes the relationships between the slope angle, expressed in degrees (˚), and the morphology, erosion and possible anthropogenic activities, providing an objective basis for classifying relief.

We added a classification of six levels of use to the class table [27] : high, high with limitations, moderate, moderate with restrictions, partly restricted and restricted. This classification indicates the amount and type of permissible activities in each class.

Our methodology was based on a spatial analysis of the information layers (PIs) in a geographic information

Text adapted from [27] ; ^*indicates added text or emphasis.

system: a) topography, b) geomorphology, c) drainage channels, d) slope, and e) ecological economic zoning. The data contained in these layers provided in information to support the final analysis.

Table 2 shows morphometric data for the drainage basin based on the Strahler classification of drainage channels, the channel bifurcation ratio (R_b), the average length of the channels in a row (L) and the ratio of the average length of the channels in different hierarchies (MLR).

The basin includes significant numbers of first- and second-order channels, indicating that the drainage headwaters have caused intense dissection of relief. The channel bifurcation ratio is slightly greater than 2, suggesting a relief configuration with a general hilly trend.

Graphing the relationship between the number of channels and the hierarchical order of the streams that constitute the basin (Figure 2) reveals that Horton’s law regarding the number of channels is applicable. According to this law, the number of successively lower-order segments of a given drainage basin tends to form a geometric progression. This progression begins with the single segment of a higher order (the 9th order in this basin) and then increases at a constant ratio. For the scale used in the study (1: 250,000), the r² indicates a 65% correlation between the number of channels and the hierarchical order.

Understanding the implications of this bifurcation ratio (R_b) for relief and, thus, anthropogenic use is of paramount importance when characterizing the environmental fragility of lower-order waterways. For hydrogeomorphological management, awareness of the intense processes of natural and anthropogenic erosion in the basin’s headwaters is required. This is especially important in the Maranhão River Basin, which includes more than 690 first-order channels, representing 49.3% of the basin’s channels.

Board 1 shows the other major indices that were important inputs for creating the suitability classes: the drainage density (D_d), river density (D_r), and extension of the surface route (E_ps). For the Maranhão River Basin, the drainage density and river density are relatively low, indicating lower relief. Thus, other elements must be observed outside the area for a hydrographic basin of this size.

The morphometric parameters of the main drainage channel can helpcharacterize the river basin. Among these, the sinuosity index (I_s) is the relationship between the channel length and the vector distance from the ends of the main channel. When I_s ≅ 1.0, the channel tends to be straight, whereas when I_s > 2.0, the channel tends to be tortuous.

Intermediate values of I_s indicate transitional, regular and irregular shapes. Sediment load, lithology, geological structures and channel slope all influence the sinuosity of the channels.

Board 1. General morphometric parameters of the Maranhão River Basin’s hydrography.

The longitudinal profile of the main channel shows its slope (Figure 3). The profile is balanced, with no sharp drops in slope. It is gently concave, and the slope is steeper in the drainage headwaters. The profile was divided into three sections in which the I_s values increasing with distance from the headwaters: the upper third has I_s = 1.23, the middle third has I_s = 1.48, and the lower third has I_s = 2.58. For the entire channel, I_s = 2.12.

The geomorphological map of Maranhão River Basin was adapted from the work of [28] who examined the effect of morphostructural processes on the distribution of features in the basin. In this study, we considered features with high degrees of dissection and, thus, higher relief to correspond to areas that should be subject to more restrictive use.

After comparing the PIs, we created a suitability relief map showing the distribution of slope classes associated with the geomorphological map and other spatial information (Figure 4). We examined the data to verify where human activities are now occurring with respect to these classes (Table 3). We used maps of Cerrado vegetation remnants, which included spatial data on anthropogenic uses of the basin (agriculture, cultivation and urban pasture).

Anthropogenic use corresponds to 25.67% of the area of the basin and is distributed among the suitability classes, indicating a high degree of control over the relief (Figure 5). In total, 70% of the anthropogenic use is concentrated in classes 2 and 3. Clearly, the difficulties imposed by high-relief areas impact use because less than 10% of the anthropogenic use occurs in areas classified with restrictions.

The classes proposed in this study were compared with the suitability classes from the Economic Ecological Zoning (ZEE) of the State of Goiás and the Distrito Federal (Table 4). We aimed to verify the possibility of including morphometric indicators in feasibility, economic and ecological studies.

Table 4 shows that the proposed suitability classes exhibit a greater degree of detail than the classes listed in the ZEE. More notably, though, the ZEE suggests that 78% (43% + 35%) of the basin is suitable for intensive

anthropogenic use. In contrast, based on morphometric parameters, this percentage is at most 40% (13.1% + 26.9%). Furthermore, the ZEE indicates that 21% of the area should be considered to be restricted, whereas the morphometric approach includes three levels of restricted activity affecting 22.6% of the area.

We conclude that applying morphometric analysis of the basin to determine the suitability of areas for anthropogenic use represents an improvement over previous methods because it results in a higher level of detail. In addition, it highlights some inconsistencies in the ZEE results.

This study sought to use spatial analysis to determine the suitability of land in a watershed for anthropogenic use. Spatial analysis of the morphometric parameters of a basin can support qualitative analysis, and this method could be incorporated into the zoning process. We suggest that this quantitative approach is necessary in the context of determining subsidies for territorial development and environmental management procedures.

It highlights the contribution of morphometric measures extracted in the scale of 1:250,000 in particular where articulated with the geomorphological features for a proposed new regional configuration of class stated for the land uses according to the degree of impact they can generate.

This text certainly is a landmark for future studies in the same watershed in more detailed scales, which may further contribute to the sustainable planning for the anthropic uses.

The authors thank their institutions Education and Research (UNB and UEPG) the opportunity to develop a collaborative work carried out under the technical cooperation agreement in the 2012 period to 2015.

Valdir A. Steinke,Ercilia Torres Steinke,Mario Diniz Araujo Neto,Maria Ligia Cassol Pinto, (2016) Proposed Relief Map of the Suitability of the Maranhão River Basin, Brazil, for Anthropogenic Use. Journal of Geographic Information System,08,351-360. doi: 10.4236/jgis.2016.83030