EFFECT OF DIGITAL ELEVATION MODEL MESH SIZE ON GEOMORPHIC INDICES: A CASE STUDY OF THE IVAÍ RIVER WATERSHED - STATE OF PARANÁ, BRAZIL

O efeito da dimensão da malha de Modelos Digitais de Elevação sobre os índices geomórficos: O estudo de caso da bacia hidrográfica do rio Ivaí - Estado do Paraná, Brasil

Abstract:

Geomorphometry is the science of quantitative description of land surface morphology by the mean of geomorphic indices extracted from Digital Elevation Models (DEMs). The analysis of these indices is the first and most common procedure performed in several geoscience-related subjects. This study aims to assess the impact of mesh size degradation on different local and regional geomorphic indices extracted for GDEM and TOPODATA DEMs. Thus, these DEMs, having a mesh size of 30 m, were subsampled to 60, 120 and 240 m and then geomorphic indices were calculated using the full resolution DEM and the subsampled ones. Depending on their behavior, these indices are then classified into stable and unstable. The results show that the most affected indices are slope and hydrographic indices such as Strahler order, stream sinuosity and fractal dimension and watershed perimeter, whereas elevation remains stable. It also shows that the effect depends on the presence of the canopy and geological structures in the studied area.

Keywords:
Geomorphic indices; Digital Elevation Model; Mesh Size; Scale

Resumo:

A geomorfometria é a ciência da descrição quantitativa da morfologia da superfície terrestre por meio de índices geomórficos extraídos de Modelos Digitais de Elevação (MDEs). A análise destes índices é um dos primeiros e mais comuns procedimentos executados em estudos relacionados a geociência. Este estudo tem como objetivo avaliar o impacto da degradação da dimensão da malha (mesh size) sobre diferentes índices geomórficos locais e regionais extraídos dos MDEs GDEM e TOPODATA. Estes MDEs, ambos com dimensão original da malha de 30 m, foram reamostrados para 60, 120 e 240 m, sendo em seguida calculados os índices geomórficos, tanto para a dimensão original quanto para as reamostras. Dependendo do seu comportamento, estes índices foram classificados em estáveis e instáveis. Os resultados mostram que os índices mais afetados são a declividade e os índices hidrográficos como a ordem de Strahler, a sinuosidade, a dimensão fractal e o perímetro da bacia, enquanto que a elevação permanece estável. Os resultados também mostram que os efeitos dependem da presença do dossel florestal e das estruturas geológicas na área estudada.

Palavras-chave:
Índices geomórficos; Modelo Digital de Elevação; dimensão da malha, escala

1. Introduction

The analysis of the morphologic aspects of an area is one of the first and most common procedures carried out in hydrological, geomorphological, environmental and territorial planning studies. This analysis can be done through the calculation of several geomorphic indices that can be then used to identify morphologically homogeneous areas. These indices were used in many studies as an auxiliary tool such as in Sodré et al. (2007), Cherem et al. (2008), Bispo et al. (2009), Thommeret (2012) and Valeriano and Rossetti (2017).

The dissemination of global-coverage, free-access Digital Elevation Models (DEMs) supported several research areas related to geomorphological mapping, specifically the quantitative morphometric description and characterization, which has benefited greatly from the high resolution of these DEMs. The use of these DEMs for the derivation of morphologic descriptors have been used in different geomorphological analysis studies to infer interpretations and numeric analyses (Valeriano and Rossetti 2011, 2012).

According to Pike, Evans and Hengl (2008), despite the fact that both classic morphometry - before DEM - and modern morphometry focus on spatial analysis, indices extraction and terrain surface analysis, a clear difference exists between both areas. Indeed, while classic morphometry especially addresses landform analysis, based on hypsometry, slope, frequency of elevation, elevation range, drainage density and topology, and terrain classification, modern morphometry focuses on improving elevation data processing, topography description and visualization, as well as numerical analysis based on digital data expressed by DEMs.

The quality of the geomorphic indices depends on the production technique and its parameters (El Hage et al. 2012), on the terrain morphology and its landcover (El Hage et al. 2017) as well as on the mesh size (Thompson, Bell and Butler 2001, Tang 2003, Kienzle 2004, Erskine et al. 2007, Vaze and Teng 2007, Wu, Li and Huang 2007, El Hage et al. 2010, Vaze, Teng and Spencer 2010, Polidori et al. 2012). The mesh size constitutes one of the main DEM characteristics and it is part of the triple scale concept manifested by dimension, extension and spacing, directly controlling the size of objects that could be extracted from DEMs (Blöschl and Sivapalan 1995). Given the importance of the scale on objects morphology representation, it is essential to study the impact of the mesh size on the information extracted for DEMs (Pike 2000, Pike, Evans and Hengl 2008).

The purpose of this study is to assess the sensitivity of the geomorphic indices extracted from DEMs to the change of the mesh size. Two DEMs are used, GDEM and TOPODATA, which are produced using two different techniques. First, the study area is presented. Then, the adopted datasets and methodologies are explained. Finally, the results are exposed and discussed.

2. Study area

The Ivaí River watershed (Figure 1), located in the State of Paraná (South Brazil), is the second largest watershed in the state, consisting of an area of about 36,533 km², which corresponds to 18% of Paraná territory. The Ivaí River is an important left bank tributary of the Paraná River formed by the confluence of Patos and São João rivers (located on Serra da Boa Esperança, municipality of Ivaí), at an elevation of about 480 m, and it flows at a distance of about 800 km towards its estuary in the Paraná River (Pontal do Tigre, municipality of Icaraíma), at an elevation of about 240 m (Destefani 2005, Meurer, Bravard, and Stevaux 2011).

Due to its heterogeneous geologic, topographic and hydrologic characteristics, Destefani (2005) divided the Ivaí River into three sectors: upper, middle and lower (Figure 1). This subdivision was adopted for this study, and it was extended to the corresponding catchment areas. Thus, 58 sub-watersheds, distributed along the three sectors and which have drainage areas greater than 100 km², were assessed.

Figure 1:
Location of the Ivaí River watershed, its topography and the three sectors.

The upper sector is generally characterized by its convex and concave slopes, flat and encased peaks, narrow and deep valleys, and a higher degree of terrain dissection. The middle sector relies on predominant convex slopes, elongated peaks, encased valleys and a lower degree of dissection compared to the upper sector. The lower sector presents convex slopes, flat peaks, open valleys and low-degree dissection, and it is where the Ivaí river floodplain is located (Maack 2002, Santos et al. 2006).

The studied watershed presents a dendritic drainage pattern. In its upper and middle sectors, however, this pattern is conditioned by geologic structures, forming rectangular channels parallel to those structures (Figure 8) (Maack 2002). The Ivaí River watershed streams, as a whole, have a winding pattern with irregular contours and sudden breaks. According to Destefani (2005), the stream sinuosity shows both higher or lower amplitude and frequency, depending on the local geology and geomorphology.

3. Data and Methods

3.1 Datasets

The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) is one of five remote sensors onboard Terra Satellite, launched by NASA in 1999. The sensor carries three telescopic subsystems that enable the acquisition of images with different spatial resolutions depending on the wavelength range. For instance, the visible and near infrared (VNIR) subsystem has a 15 m pixel size, the short wave infrared (SWIR) has a 30 m pixel whereas the thermal infrared (TIR) has a 90 m pixel (Abrams 2000, Abrams and Hook 2001). Each ASTER scene covers an area of 60 x 60 km (Abrams 2000).

The Global Digital Elevation Model (GDEM) is derived from the VNIR subsystem using nadir and backward images. It has a 30 m mesh size and it was created by processing 1.3 million ASTER scenes, covering the areas situated between 83°N and 83°S latitudes (Abrams and Hook 2001, ERSDAC 2005). These scenes are used to generate relative DEMs (without ground control points) having a vertical accuracy of 10 m, and absolute DEMs (with ground control points) having a vertical accuracy of 7 m (Abrams and Hook 2001).

TOPODATA is a DEM obtained by oversampling the 3 arc-second Shuttle Radar Topography Mission (SRTM) DEM over Brazil. The SRTM mission, which took place in the year 2000, aimed to collect altimetric data for 80% of the Earth’s land surface using SAR interferometry operating at X (λ = 3,1 cm) and C (λ = 5,6 cm) bands. This mission resulted in a global DEM with a mesh size of 1 arc-second (~30 m) for the United States and 3 arc-seconds (~90 m) for the other areas, having a vertical accuracy of 16 m and a horizontal accuracy of 20 m with 90% confidence (Rabus et al. 2003, Ludwig and Schneider 2006, Berry, Garlick and Smith 2007). TOPODATA is generated by oversampling SRTM-DEM from 90 m to 30 m using an algorithm based on kriging, a geostatistical interpolation method. The extraction methodology is further explained in Valeriano and Albuquerque (2010) and Valeriano and Rossetti (2012).

In addition to these DEMs, 1:50,000 and 1:100,000 topographic maps are used. These maps are generated basing of photogrammetric surveys conducted by IBGE (Brazilian Institute of Geography and Statistics) and DSG (the Ministry of the Army - Directorate of Geographic Services) between 1960 and 1970. The maps were then provided free of charge by the ITCG (Institute of Land, Cartography and Geosciences) for the Brazilian southern state of Paraná.

3.2 Methodology

To assess the impact of the mesh size on the geomorphic indices, the DEMs should first be subsampled. Thus, GDEM and TOPODATA are subsampled by keeping one point every two points. Therefore, the first subsampled DEM will have a mesh size of 60 m. Then, the subsampled DEM in his turn is resubsampled by the same method (one point every two), so that its new mesh size is 120 m and 240 m for the last one. This procedure is conducted with the resampling tool of ArcGIS 10.1, and the nearest neighbor interpolation is used, which assign the elevation of the closest pixel to the output pixel. Then, hydrographic and geomorphic indices are extracted from the full resolution DEMs and the subsampled ones and from the topographic maps in order to compare them.

The drainage network extraction and watershed delineation are performed using ArcGIS 10.1, in which the flow direction is determined using D8 algorithm and the extraction chain involves the following steps:

First, the DEM is filled where sinks and peaks that block the flow are removed. Such procedure avoids possible mistakes in the generated drainage network.

Then, the flow direction and flow accumulation are extracted.

Finally, the stream order is calculated and the sub-watersheds are delineated.

Initially, the geomorphic indices were classified into local and regional. Local indices characterize the value at a specific point, such as elevation and slope, whereas regional indices are descriptors of an area of the watershed, such as the form factor (Table 1).

Table 1:
Local and regional geomorphic indices used in the study.

Indices	Type	Description
Elevation	local	Vertical distance between a given point and mean sea level
Slope	local	Angle between the horizontal and the terrain surface
Confluence angle	local	Angle between two channels calculated on stream intersections
Sinuosity	regional	Ratio between stream length and the straight distance between stream ends
Fractal dimension	regional	Characterizes the complexity of streams
Strahler order	regional	Classification of the ramifications and bifurcations of the stream network in a watershed
Drainage area	regional	Total surface of a watershed
Perimeter	regional	Length of the crest line
Form factor or conformation index	regional	Ratio of the basin area to the square of the basin length Kf = A/L2

The calculations have been carried out by ArcGIS and SSM (Self-Similarity Map) software (developed by INRA (Institut National de la Recherche Agronomique)). The latter is dedicated to quantify the geometry and topology of stream networks and, therefore, it enables the quantitative analysis of any channel network (Gaucherel, Salomon and Labonne 2011).

4. Results

Elevation histograms (Figure 2) of the studied watershed show a higher number of pixels for the same elevation value for GDEM compared with TOPODATA. These small differences are related to the production techniques and methods of these DEMs (Polidori, El Hage and Valeriano 2014). This figure also shows that the increase of the mesh size does not affect the elevation histogram for the both DEMs.

Figure 2:
Elevation histograms of the full resolution GDEM (A) and TOPODATA (B) DEMs and the subsampled ones.

The slope histograms (Figure 3), show that an increase in the mesh size induces a disappearance of steep slopes and low slope angles become more abundant, which is an expected result since the terrain surface becomes smoother. This can have an important effect in studies that search to map unstable areas. Indeed, the extent of unstable regions that are characterized by the presence of steep slopes, tend to be underestimated with the increase of the mesh size as these slopes disappear (El Hage et al. 2010, Polidori el al. 2012). Moreover, the impact of subsampling is more pronounced on GDEM than on TOPODATA. This could be due the strong effect induced by the presence of canopy in GDEM that increases the number of steep slopes in the DEM, specifically with small mesh sizes (Polidori and Simonetto 2014). This effect diminishes with the increase of the subsampling ratio, and with a mesh size of 240 m, the histograms of the two DEMs become very similar.

Figure 3:
Slope histograms of the full resolution GDEM (A) and TOPODATA (B) DEMs and the subsampled ones.

To assess the impact of mesh size on the hydrography, stream networks were extracted from both DEMs. Some studies evaluated the drainage network extracted from GDEM and they concluded some imperfections in the stream morphometry (Hosseinzadeh 2011) and a good positional accuracy of the streams (Polidori et al. 2014). Figure 4 shows that there is a loss in the topologic and morphologic details of streams obtained from GDEM and TOPODATA compared with the drainage extracted from the topographic map at a scale of 1: 50,000. The difference is more significant in the upper and lower sectors for both GDEM and TOPODATA.

Likewise, the subsampling has also an effect on the drainage network for both DEMs, and this is visually clear as shown in Figure 5. This impact is manifested by a reduced number of streams having small orders and a decrease in the morphological details with slight differences between GDEM and TOPODATA.

Figure 4:
Drainage network extracted from topographic maps and from GDEM and TOPODATA in (A) the upper sector; (B) the middle sector (C) and the lower sector.

Figure 5:
Drainage networks extracted from GDEM (A) and TOPODATA (B) in the lower sector.

Figure 6 shows the relation between Strahler order (Strahler 1957) and the number of streams in the lower sector for the different mesh sizes. The streams extracted from full resolution GDEM are more numerous the those extracted from full resolution TOPODATA for all orders, which also could be an effect induced by the canopy. However, as the mesh size increase, the number of streams decreases for all the orders and for both DEMs, and the difference between the two DEMs becomes smaller.

Figure 6:
Number of streams for each Strahler order in the lower sector.

In general, the sinuosity of the streams in Ivaí River watershed is moderate, and their extent and frequency depend on local geology and geomorphology. Figure 7 represents the mean value of sinuosity and fractal dimension for the different Strahler orders and different mesh sizes. This figure shows a big difference between the different sectors with respect to the mesh size. Basing on Figure 8, a relation can be inferred between this difference and the presence of a complex system of fractures in the basin. Indeed, the upper sectors is characterized by a very dense system of structural lineaments with three main directions: NE, NW and EW (the lineaments are extracted using a method based on satellite imagery and geophysical measurements), the middle sector has a lower density and lower sector has the least density of lineaments that do not have a main direction (Wildner et al. 2006).

The sinuosity and fractal dimension of the first two Strahler orders increases with the increases of mesh size and decrease for the higher orders in the upper and middle sectors. However, in the lower sector, the increase of mesh size does not significantly affect these indices (Figure 7). The low value of sinuosity and fractal dimension for the first two orders in the case of full resolution DEMs is due to the presence of an intensive network of faults in the upstream of the upper and middle sectors (Figure 8). Indeed, these faults are linear objects and their sinuosity and fractal dimension values tend to be 1. An increase in the mesh size causes deformation of these linear structures and then an increase in their sinuosity and fractal dimension.

Figure 7:
Mean sinuosity and mean fractal dimension for all Strahler orders in the upper sector (A) the middle sector (B) and the lower sector (C).

The confluence angles in the upper and middle sectors are conditioned by local geology as shown in the Figure 9, where it is possible to visualize the influence of lineaments in the streams. Figure 10 represents the confluence angles generated with different mesh sizes. In general, the mean value of confluence angles decreases when the mesh size increases. This is partly due to the intensive presence of linear structures that with NE-NW direction that create confluence with near right angles (this is why the mean value is high). The subsampling produces a smoothing in the relief and the angles between these structures tend to decrease as a result.

Figure 8:
Structural lineaments for Ivaí River watershed (Minepar, 2006).

Figure 9:
River confluences in Ivaí River watershed in the upper sector (A) the middle sector (B) the lower sector (C) and the floodplain area (D).

Figure 10:
Mean values of confluence angles extracted from GDEM and TOPODATA.

The drainage area shows slight difference between full resolution and subsampled DEMs and also when compared to topographic maps (Figure 11 - A). A small difference is shown in the lower sector, where the terrain is very rugged and the subsampling can cause the disappearance of some terrain details and consequently the sub-basin boundaries differs (Figure 12).

By comparing the values of the perimeters extracted from topographic maps and DEMs for all the watershed sectors studied (Figure 13), an inverse relation that exist between the mesh size and the perimeter. Indeed, the more the mesh size, the less the perimeter and this is due to a decrease in the surface roughness, where noise and small landforms disappear. It is important to highlight that the perimeter values for the 30 m mesh size in both DEMs are slightly higher than those of the topographic map. This may be due to the presence of canopy that has a clearer effect on 30 m DEM than on the subsampled ones (Polidori and Simonetto 2014), specifically in the case of GDEM. Moreover, the perimeters calculated by GDEM are always higher than those calculated by TOPODATA except the last mesh size, which also implies the effect of the canopy that tends to disappear with a mesh size higher than 120 m.

Figure 11:
Drainage area (A) and form factor (B) calculated for topographic maps, GDEM and TOPODATA DEMs for the upper, middle and lower sector sub-watersheds.

Figure 12:
Boundaries of a sub-basin extracted from GDEM (A) and TOPODATA (B) in the lower sector.

Figure 13:
Perimeter of the sub-basins extracted from the topographic map, GDEM and TOPODATA in (A) upper sector; (B) middle sector (C) lower sector of the watershed.

The form factor (Figure 11-B) undergoes a small change in value with the increase of the mesh size. For example, the drainage area in the upper sector is almost stable for all DEMs, but the form factor presents a difference in value, which comes from a difference in the perimeter of the sub-watershed (Figure 13).

5. Discussion and conclusion

This study aimed to assess the impact of subsampling on different geomorphic indices extracted from DEMs. Local and regional indices were evaluated and the results show that the mesh size produces a significant effect on some indices like slope and watershed perimeter and low effect on other indices like drainage area and elevation. Table 2 represents the obtained stable and unstable indices.

The effect of mesh size depends on presence of canopy in the studied area. Indeed, the canopy modify the morphology on the DEM and the geomorphic indices that arise from are affected specifically with a small mesh. According to Polidori, El Hage and Valeriano (2014), TOPODATA carries the advantage of extending SRTM geomorphologic properties to a rather local scale, since the oversampling method of SRTM from 90 m to 30 m is based on geomorphologic assumptions and that is why the effect of the canopy is low in TOPODATA. The increase of the mesh size decreases the effect of the canopy in GDEM and with a mesh size higher than 120 m, the difference between GDEM and TOPODATA becomes very small. The impact also depends on the presence of structures in the studied area. Indeed, the lineaments are better represented in full resolution DEMs, where a subsampling tends to modify these linear features, which affects the geomorphic indices, specifically the morphometry of the stream network.

Table 2:
Stable and unstable geomorphic indices.

Stable	Unstable
	Moderately	Highly
Elevation (local)	Confluence angle (local)	Slope (local)
	Drainage area (regional)	Strahler order (regional)
	Form Factor (regional)	Sinuosity (regional)
		Fractal dimension (regional)
		Perimeter (regional)

The drainage network extracted from the 30 m GDEM and TOPODATA showed a difference in the stream network compared with that extracted from the topographic map. In the floodplain areas, the movement of water is hardly detected using a DEM as it is controlled by very small topographic features, specifically if the mesh size is not fine enough to represent these small features such as narrow channels connecting main rivers and floodplains (Yamazaki et al 2012, Saksena and Merwade 2015). Moreover, the presence of noise and canopy in the DEM highly affect the drainage network in these regions, specifically with small mesh sizes (Lane, James and Crowell 2000, Polidori and Simonetto 2014).

According to Fernández (2011), the studied DEMs have limitations in detailed drainage network extraction (scale > 1:50,000). With TOPODATA, though, it is possible to generate coherent drainage network compatible with a scale of 1:100,000; although it still requires editing work, especially in flat areas as floodplains because of a limitation in the drainage extraction algorithms. Nevertheless, in low slope angles, GDEM and TOPODATA are compatible with scale of 1:50,000, whereas in steeper slopes, like the upper and middle sectors, the compatible scale is 1:100,000 or smaller (Moura, Bias and 2014).

Flow direction extraction algorithm (D8) has some limitations and may, therefore, cause some errors. Algorithms such as D8 (unidirectional) gives better results on converging topographies. In diverging topographies, however, it generates errors since it restricts the direction of the flow towards a single pixel among its eight neighboring pixels (Fernández 2011, Fernández et al. 2012). There is also an influence caused by other factors such as DEM elevation error that propagates in the geomorphic indices (Tang et al. 2001, El Hage et al. 2017).

It is important to emphasize that the watershed in question is strongly controlled by a series of structural lineaments, which is prominently reflected on the hydrographic network and, consequently, on the local and regional indices, especially in the upper and middle sectors (Santos 2015). Another watershed located in a smoother terrain and missing intense structural lineaments could provide different results. Further research is needed to profoundly understand the impact of mesh size on the geomorphic indices over these structures.

In this study, the influence of scale change on geomorphic indices has been assessed. However, to have a better idea on how the mesh size change affect the geomorphology of the DEM, the quality of the aforementioned indices should be assessed using ground control data, and this will be investigated in future research. It is also important to point out that the results obtained depend on the extraction algorithms adopted in the study, where other algorithms would give slight differences

ACKNOWLEDGEMENTS

We thank to CAPES/COFECUB (the Coordination for the Improvement of Higher Education Personnel) for the grant and ESGT-CNAM (Ecole Supérieure des Géomètres et Topographes, Le Mans - France) for the technical support.

REFERENCES

• Abrams, M. 2000 The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER): data products for the high spatial resolution imager on NASA’s Terra platform. International Journal of Remote Sensing 21(5), pp. 847-859.
• Abrams, M. and Hook, S. 2001. ASTER user handbook. Version 2 Jet Propulsion Laboratory/EROS Data Center: California.
• Berry, P.A.M., Garlick, J.D. and Smith, R.G. 2006. Near-global validation of the SRTM DEM using satellite radar altimetry.Remote Sensing of Environment, 106 (1), pp.17-27.
• Bispo, P. C. et al. 2009. Variáveis geomorfométricas locais e sua relação com a vegetação da região do interflúvio Madeira-Purus (AM-RO). Acta Amazonica, 39(1), pp. 81-90.
• Blöschl, G. and Sivapalan, M. 1995. Scale issues in hydrological modelling: A review. Hydrological Processes, 9(4), pp. 251-290.
• Cherem, L. F. S. 2008. Análise morfométrica da Bacia do Alto do Rio das Velhas PhD. Universidade Federal de Minas Gerais
• Destefani, E. V. 2005. Regime Hidrológico do Rio Ivaí - PR Master. Universidade Estadual de Maringá.
• ERSDAC - Earth Remote Sensing Data Analysis Center. 2005. ASTER user’s guide part III DEM Product (pdf) Available at: Available at: https://unit.aist.go.jp/igg/rs-rg/ASTERSciWeb_AIST/en/documnts/users_guide/part1/pdf/Part3D_1.1E.pdf (Accessed 01 April 2017).
  » https://unit.aist.go.jp/igg/rs-rg/ASTERSciWeb_AIST/en/documnts/users_guide/part1/pdf/Part3D_1.1E.pdf
• El Hage, M. et al. 2010. Impact of DEM reconstruction parameters on topographic indices. International Archives of Photogrammetry, Remote Sensing and Spatial Information Sciences, 38(3B), pp. 40-44.
• El Hage, M. et al. 2012. Evaluation of elevation, slope and stream network quality of SPOT DEMs. ISPRS Annals of Photogrammetry, Remote Sensing and Spatial Information Sciences, 1(2), pp. 63-67.
• El Hage, M. et al. 2017. Effect of Image-Matching Parameters and Local Morphology on the Geomorphological Quality of SPOT DEMs. Photogrammetric Record, In Press, doi:10.1111/phor.12196.
  » https://doi.org/10.1111/phor.12196
• Erskine, R.H. et al. 2007. Digital elevation accuracy and grid cell size: effects on estimated terrain attributes. Soil Science Society of America Journal, 71, pp. 1371-1380.
• Fernández, D. C. J. 2011. Avaliação de algoritmos e modelos digitais de elevação para extração da drenagem Master. Instituto Nacional de Pesquisas Espaciais.
• Fernández, D. C. J. et al. 2012. Extração automática de redes de drenagem a partir de Modelos Digitais de Elevação. Revista Brasileira de Cartografia, 64 (5), pp. 619-634.
• Gaucherel, C., Salomon, L. and Labonne, J. 2011. Variable self-similar sinuosity properties within simulated river networks. Earth Surface Processes and Landforms, 36(10), pp.1313-1320.
• Hosseinzadeh, S.R. 2011. Assessing the quality of ASTER DEMs for Hydrological Applications. International Conference on Environment Science and Engineering, IPCBEE, Bali Island, Indonesia, 1-3 April 2011, vol.8, pp. 37-41.
• Kienzle, S. 2004. The effect of DEM raster resolution on first order, second order and compound terrain derivatives. Transactions in GIS, 8, pp. 83-111.
• Lane, S. N., James, T. D. and Crowell, M. S., 2000. Application of digital photogrammetry to complex topography for geomorphological research. Photogrammetric Record , 16(95), pp. 793-821.
• Ludwig, R. and Schneider, P. 2006. Validation of digital elevation models from SRTM X-SAR for applications in hydrologic modeling. Journal of Photogrammetry & Remote Sensing, 60, pp. 339-358.
• Maack, R. 2002. Geografia física do Estado do Paraná Curitiba, Imprensa Oficial.
• MINEROPAR - Minerais do Paraná. 2006. Atlas geomorfológico do Estado do Paraná Universidade Federal do Paraná. Curitiba.
• Moura, L. Z., Bias, E. S. and Brites, R. 2014. Avaliação da acurácia vertical de Modelos Digitais De Elevação (MDEs) Nas Bacias do Paranoá e São Bartolomeu. Revista Brasileira de Cartografia , 66(1), pp. 1-14.
• Meurer, M., Bravard, J. P. and Stevaux, J. C. 2011. Granulometria dos sedimentos marginais do rio Ivaí com vistas à compreensão da dinâmica hidrosedimentar montante-jusante. Revista Brasileira de Geomorfologia, 12(1), pp. 39-44.
• Pike, R.J. 2000 Geomorphometry - diversity in quantitative surface analysis. Progress in PhysicalGeography, 24(1), pp.1-20.
• Pike, R.J., Evans, I. and Hengl, T. 2008. Geomorphometry: A Brief Guide. In: Hengl, T.; Reuter, H.I. Geomorphometry: Concepts, Software Elsevier: Oxford, pp. 3-33.
• Polidori, L. et al. 2012. Scale effects in the assessment of digital terrain models and slope maps. Geotunis, 26-30 Mars 2012, Tunis.
• Polidori, L. and Simonetto, E. 2014. Effect of scale on the correlation between topography and canopy elevations in an airborne InSAR product over Amazonia. Procedia Technology, 16, pp. 180-185.
• Polidori, L. et al. 2014. Elaboration du référentiel hydrographique d'Haïti à partir d'un MNT ASTER. Revue Française de Photogrammétrie et de Télédétection, 205, pp. 49-57.
• Polidori, L., El Hage, M. and Valeriano, M. M. 2014. Digital Elevation model validation with no ground control: application to the TOPODATA DEM in Brazil. Boletim de Ciências Geodésicas, 20(2), pp. 467-479.
• Rabus, B. et al. 2003. The shuttle radar topography mission-a new class of digital elevation models acquired by spaceborne radar. ISPRS Journal of Photogrammetry and Remote Sensing, 57 (4), 241-262.
• Saksena, S. and Merwade, V. 2015. Incorporating the effect of DEM resolution and accuracy for improved flood inundation mapping. Journal of Hydrology, 530, pp. 180-194.
• Santos, V.C. 2015. Ambientes de confluência no contexto da rede de drenagem: Exemplo da bacia hidrográfica do rio Ivaí - Estado do Paraná PhD. Universidade Estadual Paulista.
• Santos, L.J.C. et al. 2006. Mapeamento geomorfológico do Estado do Paraná. Revista Brasileira de Geomorfologia , 7(2), pp. 03-12.
• Sodré, R. V. R. et al. 2007. Classificação de bacias de drenagem do alto Jequitaí (Minas Gerais) a partir da análise de principais componentes e análise de grupos. Revista Brasileira de Geomorfologia , 8(2), pp. 73-86.
• Strahler, A. N. 1957. Quantitative analysis of watershed geomorphology.Transactions of the American Geophysical Union, 38 (6), pp. 913-920
• Tang, G. et al. 2001. The impact of resolution on the accuracy of hydrologic data derived from DEMs. Journal of Geographical Sciences, 11(4), pp. 393-401
• Tang, G. et al. 2003. Simulation on slope uncertainty derived from DEMs at different resolution levels: a case study in the Loess Plateau. Journal of Geographical Sciences , 13, pp. 387-394.
• Thommeret, N. 2012. Analyse spatiale de reseaux de ravines hierarchisees a partir de MNT a differentes resolutions - Application aux badlands de draix (Alpes de Haute-Provence PhD. Université Paris 1 Pantheon-Sorbonne.
• Thompson, J.A., Bell, J.C. and Butler, C.A. 2001. Digital elevation model resolution: effects on terrain attribute calculation and quantitative soil-landscape modeling. Geoderma, 100, pp. 67-89.
• Valeriano, M. M. and Albuquerque, P. C. G. Topodata: processamento dos dados SRTM. São José dos Campos: INPE, 2010.
• Valeriano, M. M. and Rossetti, D. F. 2017. Regionalization of local geomorphometric derivations for geological mapping in the sedimentary domain of central Amazônia. Computers & Geosciences, 100, pp. 46-56
• Valeriano, M. M. and Rossetti, D. F. 2012. Topodata: Brazilian full coverage refinement of SRTM data. Applied Geography, 32, pp. 300-309.
• Valeriano, M. M. and Rossetti, D. F. 2011. Análise digital do padrão de coerência da orientação de vertentes em terrenos planos. In: XIII Congresso da Associação Brasileira de Estudos do Quaternário ABEQUA Armação dos Búzios - RJ, Brazil, 9-14 October, 2011.
• Vaze, J. and Teng, J. 2007. Impact of DEM Resolution on Topographic Indices and Hydrological Modelling Results. MODSIM 2007 International Congress on Modelling and Simulation Modelling and Simulation Society of Australia and New Zealand, December 2007.
• Vaze, J., Teng, J. and Spencer, G., 2010. Impact of DEM accuracy and resolution on topographic indices. Environmental Modelling & Software, 25 (10), 1086-1098.
• Wildner, W.; et al. 2006. Texto explicativo dos mapas geológico e de recursos minerais do sudoeste do Estado do Paraná - escala 1:200.000 Brasília: CPRM.
• Wu, S., Li, J. and Huang, G.H., 2007. Modeling the effects of elevation data resolution on the performance of topography-based watershed runoff simulation, Environmental Modelling & Software , 22, pp. 1250-1260.
• Yamazaki, D. et al. 2012. Adjustment of a spaceborne DEM for use in floodplain hydrodynamic modeling. Journal of Hydrology , 436-437, pp. 81-91.

Publication Dates

• Publication in this collection
  Oct-Dec 2017

History

• Received
  2 May 2017
• Accepted
  28 Aug 2017