Morphology and morphometry of upland lakes over lateritic crust , Serra dos Carajás , southeastern Amazon region

High-resolution satellite images, digital elevation models, bathymetric and sedimentological surveys coupled with statistical analysis were used to understand the physical environment and discuss their influence on water quality of the five upland lakes of Serra Sul dos Carajás, southeast Amazonia. The lakes have midaltitude ranges (elevation), very small (catchment) and shallow to very shallow (central basins). Based on the length, area and volume, Violão and TI (Três Irmãs)-3 lakes may present large vertical movements of the water due to wind action and weakly stratified waters. Trophic conditions based on depth and shore development (Ld) parameters must be used with caution, since Amendoim Lake is relatively deep, but it is oligotrophic to ultra-oligotrophic. Ld values suggest that the lakes are circular to subcircular and are likely formed by solution process, as also suggested by volume development. TI-2 Lake is only presenting convex central basin and has highest dynamic ratio (DR), thus it may have high sedimentation and erosion rates. Based on the relationship between studied parameters, morphometric index and DR likely influence temperature and dissolved oxygen of waters of TI-2 Lake due to its depth profile and wind-induced surface mixing. Nevertheless, water quality parameters are controlled by catchment characteristics of the lakes.


INTRODUCTION
Understanding physical characteristics of lake system related to morphometric and morphological features is crucial to comprehend its hydrological pattern, structure and ecological functionality (Håkanson 1981, Sperling 1999, Nõges 2009).The lake morphometry is highly influenced by various environmental variables such as nutrient balance, biological productivity and dispersion, thermal stability and stratification of the water column, circulation and sedimentary processes (e.g.Panosso et al. 1995, Beisner et al. 2006, Håkanson 2004, Smith et al. 2005, Genkai-Kato and Carpenter 2005, Fantin-Cruz et al. 2008, Souza Filho et al. 2016).Lake morphometric data are also required for a precise evaluation of sediment accumulation rates, dispersion patterns, sensitivity and resilience due to pollutant input in lake systems (Håkanson 1999).The interactions between these environmental variables may significantly vary depending on the lake characteristics.Therefore, large lakes may be separated into sub ecosystems where different key organisms dominate, while smaller lakes (area < 10 km 2 ) can often be regarded as single ecosystems (Håkanson 1999, Buraschi et al. 2005).The same must be expected for lowland and upland lakes of the Amazon region (Junk et al. 1989, Lopes et al. 2011, Brito et al. 2014, Sahoo et al. 2016a).
Lowland lakes are widespread on floodplains of the Amazon basin, where they are influenced by annual flood pulse of the rivers according to climate seasonality (Junk 1997).Thus, they receive huge amounts of suspended material from different sources, which are likely located hundreds of kilometers distant during high-stand conditions (Lesack and Melack 1995).On the other hand, upland lakes can only be found in scattered sectors of the Brazilian Amazon (e.g.Seis Lagos, Cordeiro et al. 2011and Serra Sul de Carajás, Sahoo et al. 2015, 2016a, b).The upland lakes of the Serra Sul dos Carajás were developed over lateritic crusts located above 600 m in altitude.These lakes were formed by structural and degradation processes of lateritic profiles (Maurity and Kotschoubey 1995), producing sinkhole (solution) lakes and catchment basins hydrologically restricted.However, studies about morphological and morphometric characteristics of upland lakes of Serra Sul de Carajás are surprisingly unknown (Sahoo et al. 2017).Therefore, the present study aims to analyze morphological pattern and morphometric characteristics and their influence on limnology of five upland lakes (Violão, Amendoim, Três Irmãs (TI-1, TI-2 and TI-3) in the Serra Sul de Carajás, southeastern Amazon region.

STUDY AREA
Violão, Amendoim, and Três Irmãs (TI-1, TI-2, and TI-3) lakes are situated in the southern ridge of the Serra dos Carajás into the Carajás National Forest (namely Flona-Carajás), a law-protected area located in the southeastern Amazon region, Brazil (Figure 1a, b).This region represents a metallogenic province on a global scale (Tolbert et al. 1971) and a major Archean tectonic province of the Amazonian Craton (Macambira andLafon 1995, Rämö et al. 2002).The large iron deposits occur associated with the Carajás Formation, a dominant volcano-sedimentary sequence that host significant band-iron formation (Olszewski et al. 1989).An extensively weathering profile was formed due to the laterization of the Carajás Formation under humid climate conditions.This developed fragile zones associated with preexisting fractures, which evolved into sinkhole and lakes (Maurity and Kotschoubey 1995).
The lakes are hydrologically restricted and developed over lateritic crust, which is vegetated by ferruginous montane savanna and bordered by Amazon rainforest (Figure 1c-g).The predominance of montane savanna from Serra Sul de Carajás is mainly related to edaphic conditions (Nunes et al. 2015, Skirycz et al. 2014, Guimarães et al. 2014).The occurrence of extensive lateritic crusts and low water retention allowed the widespread development of plant species adapted to nutrient and hydric stress and hindered the colonization of tree species.Consequently, small patches of forest (capão de mata) occur along structural features or over different lithologies that are much more sensitive to weathering and where a thicker soil horizon was developed, with higher nutrient and water availability.Open forest and forest copse associated to "capão florestal" can also be observed in central portions of the plateau over degraded crusts (Guimarães et al. 2014(Guimarães et al. , 2017)).The Três Irmãs system presents surface connection between the three lakes.
The regional climate is tropical monsoon, with a mean annual temperature of ~ 26 ºC (Alvares et al. 2014).Total annual precipitation ranges from 1420 to 2160 mm, with a total mean of ~ 1450 mm during the rainy season (November to May) and ~ 235 mm during the dry season (June to October).In addition, around 80% of the total annual occur between December to May (Silva Júnior et al. 2017).The pattern of wind velocity and direction in the Serra Sul de Carajás shows slightly higher ranges of wind velocity (~3-4 m/s) with prevailing direction from NE and ENE (Guimarães et al. 2014).

MORPHOLOGICAL SURVEY OF THE CATCHMENT BASINS
Worldview-2 multispectral satellite imagery with 2 m resolution were acquired in May 2013 for morphological mapping of the studied basins.A digital elevation model (DEM) was also developed using LiDAR (Light Detection and Ranging) data acquired in November 2009.The planimetric mapping was referenced to WGS 84 in UTM zone 22 south.The images were orthorectified and the resulting digital values were converted to terrain reflectance by the ATCOR algorithm using PCI Geomatica 13 (Geomatica Algorithm Reference 2003).This dataset was used for interpretation of the morphological features and real-time navigation during bathymetric surveys.

MORPHOMETRIC SURVEY OF THE LAKES
The bathymetric data were acquired using Raytheon DE719E and SyQuest PC 2010 echo sounders operating with a frequency of 200 kHz.The integrated acquisition of bathymetric data, with Differential Global Positioning System (DGPS, model Trimble AG-132 and SPS551, with differential correction from the OmniSTAR satellite ASAT L1 band) was carried out based on Hypack 6.2b software.The digital files were stored in *raw and *xyz extension as follow: GGA (coordinate/s), VTG (velocity and direction/s) and DBS (depth/0.1s).The bathymetric maps were produced from 567,777 sampling points in a total of 89.6 km of bathymetric lines (Table I).The points were interpolated using Kriging geostatistical method in Surfer 11 software (Golden Software 2013).The bathymetric map and digital terrain model were designed for planimetric map with 1 m isobaths of interval.The integration of bathymetric data with LiDAR topography allowed a 3D visualization of the bottom lake morphology and adjacent catchment areas of the studied lakes.Morphometric parameters including size, form and special factors were calculated to the maximum lake levels according to Håkanson (2004) (Table II).The normalized parameters were statically evaluated by multivariate exploratory techniques including factor and cluster analysis using Statistica 8 (STATSOFT Inc. 2009).The main limnological data of the studied lakes are directly related to primary productivity, trophic state and sedimentary process (temperature, dissolved oxygen-DO, total organic carbon-TOC; total iron-Fe (t) , dissolved aluminium-Al (d) , chlorophyll-a-Chl-a, total nitrogen-TN, total phosphorous-TP, Secchi depth-SD and lithological types) available in Sahoo et al. (2017), and morphometric parameters were also statistically evaluated by factor analysis, after raw data normalization.

SEDIMENTOLOGICAL SURVEY OF THE LAKES
Around 120 surficial (0-10 cm) sediment samples were collected using a Van Veen Grab sampler (Figure 1c, d), with geographical coordinates obtained from the DGPS.The sampling points were chosen based on the bathymetric maps of the lakes.These samples were described following Schnurrenberger et al. (2003).

MORPHOLOGY AND MORPHOMETRY OF THE CATCHMENTS AND CENTRAL BASINS
The catchment basin of Violão Lake presents the largest area (1.83 km 2 ) and altitude variability (164 m), while Três Irmãs system that may present surface connection between its three lakes during the rainy seasons (Figure 2a, b) shows the largest perimeter (8.13 km), but the lowest altitude variability (77.5 m).Amendoim Lake has the smallest perimeter and surface area among the studied lakes (Table III).The catchments of these lakes present relatively high declivities (> 20º and maximum of ~ 60º), with steep slopes covered by montane savanna.Slopes with ravines with capão florestal only occur in the catchment of the Violão and Amendoim lakes (Figure 2a).However, patches of forest in the Violão catchment occur over mafic sills.Several caves with guano can also be observed in this area, as well as guano over iron crust in the northwestern border near the patches of forest.The topographic profiles are illustrated in Figure 2c.Violão and Amendoim lakes are separated by an intermediate basin that prevents any surface connection of water among them.The catchment of Amendoim Lake was composed by two active lakes (Figure 2a, c, see section 7-8), but the smaller one located in the western portion was progressively filled by detritic and organic sediments, and nowadays represents a swamp colonized by macrophytes (Guimarães et al. 2017).Once the three lakes of Três Irmãs are connected, the water flow follows the elevation gradient from TI-1 to TI-3 (Figure 2b, c, see section 1-2), forming a small waterfall between them.Patches of forest or open ombrophylous forest, over aluminous crust of likely volcanic origin, are also present in the catchment of the TI-1 and TI-3.
Considering the morphology of the central basins (Figure 3), Violão Lake has an NE-SW elongated guitar-shaped form and washing basin morphology marked by a steep margin, scoured in lateritic outcrops, and flat bottom.Its western portion presents a flat, shallow (< 2.5 m depth) and ellipsoidal area (200 x 100 m) associated to embedded drainage (Souza Filho et al. 2016).This portion is linked by a narrow and elongated channel in NW-SE direction, and connects shallower and deeper zones.Fan morphology can be observed between these zones and in the southeastern portion of the lake (Figure 3).Amendoim Lake has a N-S elongated peanut form, steep margin and relatively flat bottom.Its bottom is colonized by macrophytes (Guimarães et al. 2017).The lake depocenter is located in the southern portion, where a structural high also occur.In addition, fan morphology is present in the southwestern and eastern portion.Differential position of the depocenter can also be observed in the TI-1 and TI-2 lakes of Três Irmãs system, which are located in the northeastern and western portion, respectively (Figure 3).Três Irmãs system has elongated form, and steep and smooth margins near the connection between lakes.The western margins of the TI-3 Lake are also smooth.
Based on the morphometric parameters (Tables II, III and IV), the Violão and TI-3 lakes have the highest values of L max , B max , D max , D mv ,    Line connecting the two most remote points on the shoreline.

Maximum width B max km
The straight line at the right angle to maximum length (L max ) that connects the two most remote extremities on the shoreline without crossing land.

Maximum depth D max m
Greatest known depth in the lakes.

Mean depth D mv m
Ratio between the lake volume (V) and the lake area (A).

Median depth D 50 m
Consists in the depth value covering 50% of the cumulated area of the lake.
Quartile depths D 25 -D 75 m Consists in the depth value covering 25% and 75% of the cumulated area of the lake.

Shoreline length L 0 km
It may be estimated in geoprocessing software.

Measures of area and volume
Area A km 2 It may be calculated by the satellite image interpretation, thus each point in the shoreline is a geographic coordinate.
It is based on the linear and parabolic approximation of the volume.

Special measures and measures of form
Relative depth D rel m Ratio between the maximum depth (D max ) and mean diameter of the lake.

Shore development L d *
Ratio between the shoreline and a perfect circle with area equal to the lake, used to measure the shoreline irregularity.2 (positive loadings on factor 1 and 2; Figure 4b).Amendoim and TI-1 lakes are clustered due to moderate values of the parameters of group 2. TI-2 Lake is isolated from other lakes (Figure 4a), which is related to the lowest value of the parameter of groups 1 (lowest loadings on factor 1) and highest value of group 3 (negative loadings on factor 1and 2) (Figure 4b).The hypsographic and volume curves show that half of the total areas and volumes of Amendoim, TI-1 and TI-2 lakes are located around D mv and D 50 .Half of the total area and volume of Violão Lake are observed at 9 m and 4 m, respectively, and at 5 and 8 m, respectively for TI-3 Lake (Table SI -Supplementary Material; Figure 4c, d).
The relationship between limnological and morphometric data (Tables IV and V) allowed to identify two important principal components with Eigen values > 1, which explain ~ 80% of the total variance of the dataset.The first principal component (PC1) represents 55% and the second one accounts 25% of the total variance.Thus, TOC, TN, Chl-a, pH, TP, Fe (t) , Al (d) are well-correlated with lithological types (group 1; highly positive loading on factor 2). DO and temperature are grouped with MI, DR, ET (group 2; highly negative loadings on factor 1), while SD is positively correlated with V d , L d , L 0 , A, V, D mv , D max , D rel , D 50 , D 25 , D 75 , D wb and S mv (group 3; highly positive loadings on factor 1). Ro did not present significant relationship with studied morphometric parameters.Its highest values are found in Amendoim Lake (Figure 4e).

SURFACE SEDIMENTARY COVER
The margins and shallower portions of Violão Lake are represented by lateritic crust and fragmental peat, respectively (Figure 5).The latter is also occurred at depths of ~ 9 m and near the fan morphologies observed in the morphological analysis.Organic mud occurs in the transition between fragmental peat and slightly oxidized mud that cover most of the central portion of the lake.Mud locally occurs near the outflows of the main drainages.The bottom of Amendoim Lake is mainly represented by organic mud, while mud occurs preferentially near the outflows of the main drainages.Fragmental peat is restricted to the central and deeper portions of this lake.The Três Irmãs system presents lateritic crusts in the margins of the TI-1 and TI-3 lakes, while fragmental peats are predominant in TI-2 Lake, and secondarily in the margins of the TI-1 and TI-3 lakes bordering the lateritic crust.Organic mud occurs in the central portion of the lakes, and in the western portion of the TI-3 Lake (Figure 5).

LAKE MORPHOTYPES
Lakes may present variable form related to their formation and evolution, which are dependent of interactions with environments and climate conditions (Timms 1992).A precise understanding of lake form and function is required to the well interpretation of quantitative morphometric attributes such as size, form and special factors (Håkanson 2004).
The studied lakes may be classified as midaltitude, very small and shallow according to their elevation (above mean sea-level), catchment and central basin areas, and mean depth, respectively (Buraschi et al. 2005, Sahoo et al. 2017).
The Violão and TI-3 lakes with the highest L max , B max and D max may present large vertical movements of the water due to wind action (Fantim-Cruz et al. 2008).As depth is closely related to V, as well as V with A, the high values of these parameters in the Violão and TI-3 lakes as well as their flattened volume curves may indicate restricted solar penetration with stratified waters.However, maximum temperature difference between the surface and bottom water of the studied lakes are very weak (not > 1.5 °C), which indicates weakly stratified waters (Souza Filho et al. 2016, Sahoo et al. 2017).This may be caused by a number of mechanisms, including wind-induced turbulence, high precipitation, and convention mixing might mix the water, resulting in a smaller range of water temperature.On the other hand, TI-2 showed the lowest values of depth and V, likely allowing the activity of autotroph organisms (Morais et al. 2005).
D rel ≤ 5% may suggest complete circulation of lake waters under the presence of relatively strong winds, which allow interchanges between water and bottom surface (Sperling 1999).This is particularly true for lakes with the main axis oriented to the main wind direction (Panosso et al. 1995).In the study site, Violão, TI-1 and TI-3 may be more influenced by winds, since theirs prevailing directions are from NE and ENE (Guimarães et al. 2014).In addition, margins of the Três Irmãs system are smoother and may be more sensible to wind action.However, wind velocities in the studied area are around 3 to 4 m/s, which only produce gentle breeze under such velocities  according to the Beuafort scale (e.g.Saucier 1955), and it may not influence bottom sediment.
Very shallow D wb of Violão, Amendoim, TI-1 and TI-3 lakes, considering D max , D mv and D 50 , are indicative of extensive accumulation areas with restricted areas of physical erosion (e.g.Håkanson 2004).In opposition, TI-2 Lake with similar values of D wb , D mv and D 50 may be more sensible to erosive process.
The highest MI was found in TI-2 Lake, since it has lower D mv and moderate value of L 0 , which indicates that the lake occupies a relatively large area, but it is very shallow and may be highly influenced by hydrodynamic processes.
DR and ET values indicate wind-wave influence on lake bottom and bottom resuspension with DR higher than 3.8.Furthermore, DR of 0.25 suggests that most of the lakes have at least 15% of ET areas (Håkanson and Jansson 1983).However, only TI-2 Lake presented significantly higher DR, which indicates an increased percentage of ET areas.
Considering the surface sedimentary cover, relation of morphometric parameters described above with sedimentological process, it is evident that bottom sediment of all lakes are deposited under very low energy conditions of waters (i.e.suspension process) and they are less influenced by post-depositional processes.This is supported by the geochemical study carried out in Carajás lakes by Sahoo et al. (2015Sahoo et al. ( , 2016b)), where the heavy minerals mainly Si (detrital), Ti and Zr were distributed at the margin of the lake and decreases their concentrations towards the center of the lake.
Low energy flow leads to sorting effect of heavy minerals according to the hydraulic properties and low depositional process (Cingolani et al. 2003).Similarly, despite having high enrichment of Fe and organic matter, the association of trace and REEs with Al-clays in those lake sediments indicates less influence of depositional and post-depositional process (including the erosion and resuspension process) on elemental chemistry (Sahoo et al. 2015(Sahoo et al. , 2016b)).
According to the model of trophic classification based on D mv (Thienemann 1927, Cole 1975), the studied lakes may be grouped as eutrophic (depths < 18 m), with relatively high nutrient availability and euphotic zone.However, trophic conditions may primarily consider geoenvironmental characteristics of the catchment basins such as lithology, soils and vegetation types (Lane et al. 2007).In such case, nutrient input in the Amendoim Lake is significantly low due to catchment lithology, mainly represented by iron-rich crusts, resulting in oligotrophic to ultra-oligotrophic conditions (Sahoo et al. 2016a).
L d > 3 may be indicative of well-developed shore communities less vulnerable to eutrophication (Håkanson 1981), but it is not the case of the studied lakes classified as highly vulnerable (L d from 1.3 to 1.5; Table IV).However, this classification must be viewed with care, since Amendoim Lake is oligotrophic to ultra-oligotrophic (Sahoo et al. 2016a).Thus, it is more suitable the use of L d as a direct proportion of L 0 and A, as well as higher values may be applied to describe more elongated lake forms (Wetzel 2001).L d may also be used to indicate the origin of lakes (Hutchinson 1957, Timms 1992).Therefore, the studied lakes were classified as circular to subcircular.In addition, these lakes have very low degree of shoreline irregularity (Timms 1992), very similar with lakes formed by solution process (Hutchinson 1957 (Thienemann 1927), since the former usually have a greater productive, if they are not limited to light (Magyar et al. 2013).Shallow lakes present large dynamic of ecological conditions as all depths may be affected by sunlight.In addition, they may be holomictic (monomictic and polymictic) with short stable stages (Crisman et al. 2005).
Trophic state (TS) is an indication of a productivity and nutrition levels, which is calculated based on three key water quality parameters such as Chl-a, TP and SD following the modified Carlson index (TSi, Carlson 1977, Lamparelli 2004, Molisani et al. 2010).According to Sahoo et al. (2017), Violão, TI-1 and TI-2 lakes are mainly eutrophic, while Amendoim and TI-3 lakes are oligotrophic and mesotrophic, respectively.Based on the relationship between geomorphology and hydrochemistry of the studied lakes and trophic state variables, TN, pH, Fe (t) and Al (d) are mainly controlled by the catchment characteristics of the lakes, which included lithotypes, drainage system and vegetation types.TOC and TN may also have derived from organic matter of water-living organisms (e.g.Sahoo et al. 2015).Autotrophic organisms may develop in lake waters once organic and inorganic materials relatively enriched in nutrients are transported by surface runoff into the central basins.However, Amendoim Lake has poor-nutrient water, since its catchment basin is only represented by iron crust that is poor in nutrients and very resistant to chemical and physical weathering (e.g.Sahoo et al. 2015Sahoo et al. , 2016a, b), b).Also, its main drainage is related to a filled lake that acts as a nutrient filter.Therefore, sunlight may easily reaches the bottom sediments allowing the growth of macrophytes.
Lakes with great water volume may present high content of DO and consequently high capacity to assimilate organic components (Sperling 1999).
However, MI and DR likely influence temperature and DO of waters of the TI-2 Lake (shallower), since it is the shallowest and hydrodynamically active lake of the studied site.Under this condition, lake water may be mixed by the wind allowing warming and oxygenation of the bottom (Sahoo et al. 2016a).
SD is a rapid way to evaluate the underwater light penetration by visual measure (Luhtala and Tolvanen 2013).Suspended particulate, dissolved matter (inorganic or organic) and chlorophyll along with other photosynthetic pigments of living phytoplankton attenuate light penetration in the water column of the lakes (Kirk 2011).Thus, SD may also indicate thickness of the photosynthetically active water layer (Preisendorfer 1986).In this work, SD represents a mixture of organic and inorganic components.V d , L d , L 0 , A, V, D mv , D max and D rel are associated with horizontal and vertical measurements of the lakes, and well correlated to SD.Thus, larger and deeper lakes such as Violão and TI-3 have high transparency, while TI-2 has the lowest.Nevertheless, Amendoim Lake is relatively deep and has the highest water transparency among the studied lakes, but it is ultra-oligotrophic to oligotrophic (Sahoo et al. 2016a).

CONCLUSION
The upland lakes Violão, Amendoim and Três Irmãs (TI-1, TI-2 and TI-3) are classified as midaltitude, very small and shallow based on elevation, catchment and central basin areas.Morphometric parameters related to length, area and volume suggest that Violão and TI-3 lakes may present large vertical movements of the water due to wind action and weakly stratified waters.Although, depth and shore development (L d ) parameters are directly used to indicate trophic conditions, this may not be applicable for the present study, because the Amendoim Lake is relatively deep but oligotrophic to ultra-oligotrophic.The lakes are circular to subcircular and concave, likely formed by solution process.Convex central basin and highest dynamic ratio are only observed for TI-2 Lake, which suggest higher influence of hydrodynamic.Therefore, the bottom sediment of all lakes were deposited under very low energy conditions of waters (i.e.suspension process) and they are less influenced by post-depositional processes, except LTI-2 Lake.The relationship between geomorphological and limnological characteristics of the lakes, and temperature and dissolved oxygen of waters of the TI-2 Lake are may be associated with morphometric index and dynamic ration, since it is the shallowest and hydrodynamically active lake of the studied site.Nevertheless, lake water quality may be primarily influenced by lithotypes, soil availability, vegetation and drainage system of the catchment basins.

Figure 1
Figure 1 -a) Location map of the study site in the southeastern Amazon region (Imazon 2003) with a detail to b) Carajás National Forest (Flona-Carajás, Landsat image) and a detail to Serra Sul de Carajás and studied lakes; c) WorldView images of the Violão and Amendoim and d) Três Irmãs lakes in the eastern and western portion of the Serra Sul; sediment samples are represented by white dots; e) Aerial photography of the Violão, f) Amendoim and g) Três Irmãs lakes.

Figure 2 -
Figure 2 -Geomorphological maps of the catchment basins of the a) Violão, Amendoim and b) Três Irmãs lakes according Rezende and Barbosa (1972), and Golder (2010), with c) topographic sections of each basin; altitude and distance units are in meters.
relationship between the shoreline lenght (L 0 ) and the mean depth (D mv ).Volume development V d * Ratio between the water volume (V) and the volume of a cone, with the base equal to the water surface area (A) and with a height equal to the maximum depth (D max ).Dynamic ratio DR * DR is a standard morphometric parameter in contexts of lake resuspension based area (A) and mean depth (D mv ).Energy-topography factor ET * It represents the fraction of the lake bed where resuspension is likely to appear, which is dependent of DR.Depth of the wave base D wb m Depth separating transportation area (with resuspension) from accumulation area (without resuspension).Lake bottom roughness Ro * It measures the degree of irregularity of the lake bottom.25 , D 75 , L 0 , A, V, D rel , L d , D wb and S mv , while TI-2 lake shows the lowest values of D max , D mv , D 50 , D 25 , D 75 , V, D rel , V d , L d and S mv , and highest MI.Although TI-1 Lake presents the lowest values of L max , B max , L 0 , A and D wb , it has the highest V d among the studied lakes.Ro has maximum and minimum values for Amendoim and Violão lakes, respectively.This is clearly observed in the multivariate analysis of these data (Figure 4a, b).Two important principal components (PC) could be defined with Eigen values > 1, which explained ~ 88% of the total variance of the dataset where PC1 and PC2 account 73% and 15% of the total variance, respectively.Violão and TI-3 lakes are clustered together (Figure 4a) based on the highest values of the morphometric parameters of groups 1316 MARCIO S. DA SILVA et al.

Figure 4 -
Figure 4 -Statistical analysis of morphometric parameters of the studied lakes according to a) similarities between lakes by cluster analysis, and b) relationship between parameters by factor analysis; B max : maximum width; D max : maximum depth; D mv : mean depth; D 50 : median depth; D 25-75 : quartile depths; L 0 : shoreline length; A: area; V: volume; D rel : relative depth; L d : shore development; MI: morphometric index; V d : volume development; DR: dynamic ratio; ET: energy-topography factor; D wb : depth of the wave base; Ro: lake bottom roughness; S mv : mean slope; c) Hypsographic and d) volume curves considering cumulative areas, volumes and depths of the studied lakes; e) Relationship between limnology and morphological and morphometric parameters of the lakes by factor analysis; D max : maximum depth; D mv : mean depth; L 0 : shoreline length; A: area; V: volume; D rel : relative depth; L d : shore development; MI: morphometric index; V d : volume development; DR: dynamic ratio; Temp: temperature; DO: dissolved oxygen; TOC: total organic carbon; Fe (t) : total iron; Al (d) : dissolved aluminium; Chl-a: chlorophyll-a; TN: total nitrogen; TP: total phosphorous; litho: lithological types; SD: Secchi depth.

TABLE I Sampling points and lines generated during the bathymetric surveys of the studied lakes.
MARCIO S. DA SILVA et al.

TABLE V Main limnological data of the studied lakes related to primary productivity. Trophic state and sedimentary process available in Sahoo et al. (2017). Lakes
Filho et al. 2016ery low V d values reinforced this interpretation (Wetzel 2001).MARCIO S. DA SILVA et al.High Ro of Amendoim Lake is directly related to the presence of a structural high in its southern portion.High S mv as observed for Violão and TI-3 lakes is indicative that both lakes present washing basin morphology (e.g.SouzaFilho et al. 2016).