Water Erosion Processes on the Geotouristic Trails of Serra da Bocaina National Park Coast, Rio de Janeiro State, Brazil

: Conservation units are strategic territories that have a high demand for public use, as they protect attractions of great scenic beauty, geodiversity sites, and numerous leisure areas. However, when carried out in an intensive and disorderly manner, tourist activity in these areas tends to catalyze environmental degradation, triggering, for example, water erosion processes caused by intensive soil trampling on the trails. In this sense, the aim of this study was to determine the soil’s physicochemical characteristics, and to spatiotemporally monitor the microtopography of those areas degraded by erosion along two trails on Serra da Bocaina National Park coast of the Paraty Municipality. The findings verified that intensive trampling, the values of some soil physicochemical characteristics, and the specific meteorological conditions of the coastal region of this protected area were factors that contributed significantly to the evolution of erosion features monitored on these trails. Finally, strategies for appropriate management and recovery actions for these degraded areas are proposed in order to not only stop the erosive processes and re-establish the local ecosystem balance, but also avoid accidents involving the numerous tourists who visit the coastal region.


Introduction
Water erosive processes cause soil degradation through the breakdown, transport, and deposition of particles by the action of water from raindrops and runoff, which corresponds with ~72% of the land affected by soil erosion in the world [1].Its occurrence is intensifying due to climate change and human actions; therefore, it is one of the main environmental and socioeconomic challenges that society faces today.When induced by anthropic actions, these processes, in addition to impacting biogeochemical cycles, promote biodiversity loss, declines in agricultural productivity, carbon storage, increases in hunger, food insecurity, poverty, and social inequality [2][3][4][5][6].
Soil trampling on trails for public use, for example, is a human activity that causes soil degradation and occurrence of these processes, especially when these actions are carried out in an intensive and disorderly way.These activities corroborate with the alteration of the characteristic physicochemical and biological changes in the soil, which, in turn, modify the hydraulic and edaphic dynamics, since, in addition to reducing organic matter levels, porosity values, and water infiltration rates, they cause an increase in bulk density, runoff volume, and soil loss due to erosion [7][8][9][10].
There are many networks of local, regional, national, and international trails maintained and built to offer leisure and recreation opportunities, especially in protected areas, such as conservation units (UCs).These territories have experienced considerable growth in the number of visitors, as their trails are spaces that allow access to areas of leisure and scenic beauty.Consequently, these places have become highly susceptible to soil degradation through water erosive processes [11][12][13][14][15].
Fonseca Filho [16], for example, monitored soil compaction and erosion on trails in the UCs of Minas Gerais State and identified intensive trampling as one of the causes of trail degradation.Rangel et al. [17] also demonstrated the influence of human trampling on soil erosion in two UCs in southeastern Brazil.Figueiredo and Martins [18] discussed the mechanisms that lead to trail degradation in protected areas, pointing out that a lack of planning, management, and maintenance is associated with the occurrence of erosion processes in these areas.Lima and Guerra [10] studied the Morro Dois Irmãos Trail, which is part of a protected area in Rio de Janeiro Municipality, and identified several erosion features along its path, pointing to the physical-chemical characteristics of the soil together with the intensive trampling by people as the main causes of the development of rills.Lima et al. [9] also identified the evolution of erosion features on trails on the Serra da Bocaina National Park (PNSB) coast, indicating intense and disorderly visitation as one of the factors that corroborates the dynamics of this behavior.
Monitoring and analyzing the microtopography and soil physical-chemical characteristics in time and space on trails with water erosion processes, especially those located inside UCs, are extremely useful for the management and planning of these protected areas, because surveys can not only support decision making by their managers in actions that aim to identify and minimize these processes, but also help in optimizing management practices in degraded areas, in addition to guiding the re-establishment of their ecosystem balance through erosion processes in these territories [14,19].
In this context, this study aims to conduct a spatial survey of the soil's physical and chemical characteristics on trails on the coast of the Serra da Bocaina National Park (PNSB) in Paraty Municipality of Brazil, in addition to evaluating the space-time evolution of water erosion processes through monitoring associations between soil microtopography and rainfall data.This protected area, located in one of the largest stretches of continuous and conserved forests in the Atlantic Forest, has numerous geotouristic attractions in its territory, which, in turn, makes it a scene of high demand for tourists, with an intense flow of people on its trails during the summer and holidays (vacations).Its trails and beaches receive an average of around eight thousand people per day [20], making the high demand for visits an agent that triggers the environmental degradation of the trails [9,17,[20][21][22].

Study Area
The PNSB, covering ~106,000 ha, is a UC that was created in 1972 in São Paulo State (Areias, Cunha, São José do Barreiro and Ubatuba Municipalities) and Rio de Janeiro State (Paraty and Angra dos Reis Municipalities).It has the objectives of preserving and conserving forest formations and associated refuges; ensuring the maintenance of natural landscapes; protecting biodiversity and water resources; and providing opportunities for scientific research, leisure, and recreation [23,24].This UC is one of the largest protected areas of the Atlantic Forest, the largest continuous extension, and one of the most significant Brazilian remnants of this biome, with several geotouristic attractions (Figure 1).On the PNSB coast, the trails which we will focus on in this study are Waterfall of Stone that Swallows (WSS) (red) and the Caixa D'Aço Natural Pool (CDN) (orange) (Figure 2).On the PNSB coast, the trails which we will focus on in this study are Waterfall of Stone that Swallows (WSS) (red) and the Caixa D'Aço Natural Pool (CDN) (orange) (Figure 2).On the PNSB coast, the trails which we will focus on in this study are Waterfall of Stone that Swallows (WSS) (red) and the Caixa D'Aço Natural Pool (CDN) (orange) (Figure 2).These trails are located in the relief of the mountains and hills, under soils and surface residuals, with rocky outcrops of granite Paraty-Mirim and colluvial-alluvial sediments.In general, they have a base saturation of less than 50%, which, in turn, gives them a dystrophic character, in addition to high acidity and high aluminum content.Also, according to mapping carried out at a scale of 1:250,000 [25], the soils that predominate the trails are dystrophic Haplic Cambisols.These classes, as well as Litholic Neosols, also occur in the mountainous region of the PNSB, on the escarpments and on the coastal strip between Serra do Mar escarpment and the plains.Oxisols, on the other hand, occur at the top and on smoother slopes [9,23,26].
The Swallows trail has a length of ~600 m, with an average slope of approximately 12.4% and an altitude range between 0 and 46 m (Figure 3), It gives access to the WSS (Figure 4).The CDN Trail is also ~600 m long, with an average slope of around 14.3% and an altitude range of 0 to 63 m (Figure 3), and it provides access to the CDN (Figure 4).These trails are located in the relief of the mountains and hills, under soils and surface residuals, with rocky outcrops of granite Paraty-Mirim and colluvial-alluvial sediments.In general, they have a base saturation of less than 50%, which, in turn, gives them a dystrophic character, in addition to high acidity and high aluminum content.Also, according to mapping carried out at a scale of 1:250,000 [25], the soils that predominate the trails are dystrophic Haplic Cambisols.These classes, as well as Litholic Neosols, also occur in the mountainous region of the PNSB, on the escarpments and on the coastal strip between Serra do Mar escarpment and the plains.Oxisols, on the other hand, occur at the top and on smoother slopes [9,23,26].
The Swallows trail has a length of ~600 m, with an average slope of approximately 12.4% and an altitude range between 0 and 46 m (Figure 3), It gives access to the WSS (Figure 4).The CDN Trail is also ~600 m long, with an average slope of around 14.3% and an altitude range of 0 to 63 m (Figure 3), and it provides access to the CDN (Figure 4).Geologically, the rocks are from the Rio Negro Magmatic Arch (orthogneisses) and the Parati-Mirim Granite Suite (granites), as well as later fluvio-marine sediments (Cenozoic).Regarding climate, this is influenced by the orographic effect of the Serra do Mar, as the compartmentalization of the relief and the altimetric unevenness of the mountain range.The mountain range extends from sea level to altitudes exceeding 2000 m and acts as a barrier to the air masses responsible for the behavior of meteorological phenomena.The escarpments of Serra do Mar and its mountainous plateau make passage difficult for frontal systems (cold fronts) coming from the South Atlantic/Antarctica, and, as such, are responsible for regional rainfall, generating temporal space-time discontinuities in these phenomena [24,27].Therefore, the Serra do Mar slopes facing the ocean (south) generate a strong seasonality in the precipitation on PNSB coast, making the frontal systems (cold fronts) and the rains frequent in the summer, while winter precipitations are minimal [23,24,27].

Methodology
The determination of the soil physicochemical properties (e.g., granulometry, porosity, bulk density, pH, and organic matter content) were chosen because they are related to erodibility and soil use and management, given that their values are influenced and/or modified by intensive trampling of the ground, mainly by the flow of tourists where the surveyed trails are located [9,10,17,22].
Soil samples were collected at depths between 0 and 20 cm at strategic points along the trails: in places on the beds of the trails where there are erosive features, that is, areas that suffer trampling, and those areas immediately adjacent (edge) to the trails, where there is no passage of visitors (Figures 5 and 6).Therefore, it was possible to infer the impact of trampling by people and to spatially compare soil quality.
Undisturbed samples were collected to determine bulk density and porosity, and disturbed samples to determine pH, granulometry (texture), and organic matter content.Geologically, the rocks are from the Rio Negro Magmatic Arch (orthogneisses) and the Parati-Mirim Granite Suite (granites), as well as later fluvio-marine sediments (Cenozoic).Regarding climate, this is influenced by the orographic effect of the Serra do Mar, as the compartmentalization of the relief and the altimetric unevenness of the mountain range.The mountain range extends from sea level to altitudes exceeding 2000 m and acts as a barrier to the air masses responsible for the behavior of meteorological phenomena.The escarpments of Serra do Mar and its mountainous plateau make passage difficult for frontal systems (cold fronts) coming from the South Atlantic/Antarctica, and, as such, are responsible for regional rainfall, generating temporal space-time discontinuities in these phenomena [24,27].Therefore, the Serra do Mar slopes facing the ocean (south) generate a strong seasonality in the precipitation on PNSB coast, making the frontal systems (cold fronts) and the rains frequent in the summer, while winter precipitations are minimal [23,24,27].

Methodology
The determination of the soil physicochemical properties (e.g., granulometry, porosity, bulk density, pH, and organic matter content) were chosen because they are related to erodibility and soil use and management, given that their values are influenced and/or modified by intensive trampling of the ground, mainly by the flow of tourists where the surveyed trails are located [9,10,17,22].
Soil samples were collected at depths between 0 and 20 cm at strategic points along the trails: in places on the beds of the trails where there are erosive features, that is, areas that suffer trampling, and those areas immediately adjacent (edge) to the trails, where there is no passage of visitors (Figures 5 and 6).Therefore, it was possible to infer the impact of trampling by people and to spatially compare soil quality.
Undisturbed samples were collected to determine bulk density and porosity, and disturbed samples to determine pH, granulometry (texture), and organic matter content.Concerning the depth of these collections, these are justified because they are the main factors that suffer from the impact of trampling during public use activities.
All soil physical-chemical parameters were assessed at the Maria Regina Mousinho de Meis Geomorphology Laboratory of the Geography Department of the Federal University of Rio de Janeiro (UFRJ).Our techniques were in accordance with the methods of Teixeira et al. [28], that is, the texture was assessed by total dispersion of individual soil particles (pipette method); the organic matter content using the indirect method (oxidation); the bulk density by collecting samples in a cylinder of known volume (100 cm 3 ); the mineral density via volumetric flask and ethyl alcohol; the porosity according to its relationship with density values; and the pH using a digital meter.Concerning the depth of these collections, these are justified because they are the main factors that suffer from the impact of trampling during public use activities.All soil physical-chemical parameters were assessed at the Maria Regina Mousinho de Meis Geomorphology Laboratory of the Geography Department of the Federal University of Rio de Janeiro (UFRJ).Our techniques were in accordance with the methods of Teixeira et al. [28], that is, the texture was assessed by total dispersion of individual soil particles (pipette method); the organic matter content using the indirect method (oxidation); the bulk density by collecting samples in a cylinder of known volume (100 cm 3 ); the mineral density via volumetric flask and ethyl alcohol; the porosity according to its relationship with density values; and the pH using a digital meter.
Determination of the texture (granulometry) took place via crushing, sieving, and mixing the soil samples with sodium hydroxide and distilled water [28], while the Concerning the depth of these collections, these are justified because they are the main factors that suffer from the impact of trampling during public use activities.All soil physical-chemical parameters were assessed at the Maria Regina Mousinho de Meis Geomorphology Laboratory of the Geography Department of the Federal University of Rio de Janeiro (UFRJ).Our techniques were in accordance with the methods of Teixeira et al. [28], that is, the texture was assessed by total dispersion of individual soil particles (pipette method); the organic matter content using the indirect method (oxidation); the bulk density by collecting samples in a cylinder of known volume (100 cm 3 ); the mineral density via volumetric flask and ethyl alcohol; the porosity according to its relationship with density values; and the pH using a digital meter.
Determination of the texture (granulometry) took place via crushing, sieving, and mixing the soil samples with sodium hydroxide and distilled water [28], while the Determination of the texture (granulometry) took place via crushing, sieving, and mixing the soil samples with sodium hydroxide and distilled water [28], while the classification was based on the textural triangle of the United States Department of Agriculture (USDA) [29].Wet organic matter was obtained by mixing soil samples with potassium dichromate, silver sulfate, orthophosphoric acid, diphenylamine, and ferrous ammonium sulfate [28].Bulk density was obtained by collecting samples in a cylinder with a known volume (100 cm 3 ), while particle density was obtained by mixing samples with ethyl alcohol [28].Finally, the porosity was determined by the ratio between the particle density (g cm −3 ) and bulk density (g cm −3 ), while the pH was obtained by dissolving the soil in distilled water and its respective reading after equipment calibration [28].
Each physical-chemical parameter was assessed in triplicate.At each site studied on the trails, three collection repetitions were performed, and the presented results refer to the Soil Syst.2024, 8, 24 6 of 18 averages of these triplicates.In addition to these averages, the values of the standard deviation and the coefficient of variation were also calculated in order to obtain the statistical variations using Microsoft Excel software.
The soil microtopography was determined through the installation of erosion bridges (PDE) with stakes from one edge to the other in the cross sections of sites and on the beds of the trails that presented water erosion processes, following the method developed by Shakesby [30] and adopted by Ferreira [31] in studies on water erosion in soils and subsequently applied to trails [9,22,32].Thus, the evolution of erosion features was monitored between August 2021 and March 2022 and between August 2022 and March 2023.Therefore, we were able to obtain significant and contrasting variations in the temporal data on rainfall in the region.It is important to note that the monitored sites were the same sites from which the soil samples were collected to obtain the physical-chemical characteristics (Figures 5 and 6).
The model developed by Shakesby [30] and adapted for trails by Silva and Botelho [32] was used, whereby 50 cm wooden stakes/rods (useful for leveling), 2 m slats (erosion bridge), and a 1 m iron rod (measuring rod) were used with a measuring ruler offering 100 holes (analysis points) at 2 cm intervals.Two stakes were installed at the edges of the cross section (sites) so that they could be leveled between the edges of the trails and the values could be measured (Figure 7).
with a known volume (100 cm 3 ), while particle density was obtained by mixing samples with ethyl alcohol [28].Finally, the porosity was determined by the ratio between the particle density (g cm −3 ) and bulk density (g cm −3 ), while the pH was obtained by dissolving the soil in distilled water and its respective reading after equipment calibration [28].
Each physical-chemical parameter was assessed in triplicate.At each site studied on the trails, three collection repetitions were performed, and the presented results refer to the averages of these triplicates.In addition to these averages, the values of the standard deviation and the coefficient of variation were also calculated in order to obtain the statistical variations using Microsoft Excel software.
The soil microtopography was determined through the installation of erosion bridges (PDE) with stakes from one edge to the other in the cross sections of sites and on the beds of the trails that presented water erosion processes, following the method developed by Shakesby [30] and adopted by Ferreira [31] in studies on water erosion in soils and subsequently applied to trails [9,22,32].Thus, the evolution of erosion features was monitored between August 2021 and March 2022 and between August 2022 and March 2023.Therefore, we were able to obtain significant and contrasting variations in the temporal data on rainfall in the region.It is important to note that the monitored sites were the same sites from which the soil samples were collected to obtain the physicalchemical characteristics (Figures 5 and 6).
The model developed by Shakesby [30] and adapted for trails by Silva and Botelho [32] was used, whereby 50 cm wooden stakes/rods (useful for leveling), 2 m slats (erosion bridge), and a 1 m iron rod (measuring rod) were used with a measuring ruler offering 100 holes (analysis points) at 2 cm intervals.Two stakes were installed at the edges of the cross section (sites) so that they could be leveled between the edges of the trails and the values could be measured (Figure 7).From this monitoring program, graphs were created using Excel software to demonstrate the space-time evolution of erosion features during the four monitoring periods.We aimed to identify sites of removal and accumulation of sediments and organic matter and places where there was intense trampling and removal of organic and inorganic particles, in addition to estimating areas of lost soil along the beds of the trails.After preparing and analyzing these graphs, it was possible to estimate the area of the monitored cross-sections and quantify the evolution of the lost soil areas at the monitoring sites.
Rainfall data were obtained using a manual rain gauge installed at a campsite near the PNSB coast (latitude 23 • 21 ′ 5.28 ′′ S and longitude 44 • 43 ′ 34.47 ′′ W) at a height of ~1.5 m on a wooden platform, away from obstacles, so that it was possible to store and record daily rainfall throughout the monitoring period.Finally, it should be noted that the records of these data correspond to the dates of the beginning and end of the monitoring program of the five erosion bridges (August 2021, March 2022, August 2022 and March 2023), and the readings were recorded daily at 11:00 am.

Rainfall Data
Daily rainfall data allowed total monthly volumes to be calculated (Figure 8).Moreover, rainfall occurred throughout less than half of the total monitoring period (219 days), with daily values higher than the total average, which was 11 mm, being recorded on more than half of these days (Figure 8).It was also possible to identify that this precipitation was concentrated in December (2021), January (2022), April (2022), November (2022), and January (2023), which in turn were the months with the highest averages of rainfall (21 mm, 19 mm, 21 mm, 17 mm, and 17.1 mm, respectively) (Figure 8).In addition, the highest monthly total volumes were also recorded in these months, that is, 647.5 mm in December 2021, 593 mm in January 2022, 630 mm in April 2022, 510 mm in November 2022, and 530 mm in January 2023 (Figure 8).
After preparing and analyzing these graphs, it was possible to estimate the area of the monitored cross-sections and quantify the evolution of the lost soil areas at the monitoring sites.
Rainfall data were obtained using a manual rain gauge installed at a campsite near the PNSB coast (latitude 23°21′5.28″Sand longitude 44°43′34.47″W)at a height of ~1.5 m on a wooden platform, away from obstacles, so that it was possible to store and record daily rainfall throughout the monitoring period.Finally, it should be noted that the records of these data correspond to the dates of the beginning and end of the monitoring program of the five erosion bridges (August 2021, March 2022, August 2022 and March 2023), and the readings were recorded daily at 11:00 am.

Rainfall Data
Daily rainfall data allowed total monthly volumes to be calculated (Figure 8).Moreover, rainfall occurred throughout less than half of the total monitoring period (219 days), with daily values higher than the total average, which was 11 mm, being recorded on more than half of these days (Figure 8).It was also possible to identify that this precipitation was concentrated in December (2021), January (2022), April (2022), November (2022), and January (2023), which in turn were the months with the highest averages of rainfall (21 mm, 19 mm, 21 mm, 17 mm, and 17.1 mm, respectively) (Figure 8).In addition, the highest monthly total volumes were also recorded in these months, that is, 647.5 mm in December 2021, 593 mm in January 2022, 630 mm in April 2022, 510 mm in November 2022, and 530 mm in January 2023 (Figure 8).When considering the intervals of days between each monitoring period, it was also possible to identify that, from the first monitoring (2021) to the second monitoring period (2022/1), the daily rainfall average was 11 mm, and these precipitations were concentrated within just 85 days.From the second monitoring period (2022/1) to the third (2022/2), despite the daily average rainfall being lower than the total average (11 mm), that is, equal When considering the intervals of days between each monitoring period, it was also possible to identify that, from the first monitoring (2021) to the second monitoring period (2022/1), the daily rainfall average was 11 mm, and these precipitations were concentrated within just 85 days.From the second monitoring period (2022/1) to the third (2022/2), despite the daily average rainfall being lower than the total average (11 mm), that is, equal to 8.3 mm, rainfall was concentrated within only 33 days.From the third monitoring period (2022/2) to the fourth (2023), the daily average rainfall was greater than 13 mm, that is, greater than the total average, with rainfall concentrated within 95 days (Figure 8).In this way, in addition to the rainfall concentration, in specific months, we noted that, between the intervals of the monitoring periods, with the exception of the second to the third, the daily average rainfall was equal to or greater than the total average precipitation during the entire monitoring period.This indicates an ideal scenario for the evolution of erosion features due to splash erosion and the increase in runoff amount on the beds of the trails.

Spatial Variation of Soil Physical-Chemical Characteristics
The survey of the spatial variation in the physical-chemical soil characteristics in the degraded sites due to water erosion processes along the WSS Trail and the CDN Trail on the PNSB coast allowed us to identify disparities between the values of the attributes that were obtained, as well as the relationship between these parameters and the evolution of erosion features.In general, porosity and bulk density values, as well as granulometry, pH, and organic matter content, showed significant relationships with the presence or absence of rills along the two trails.
Analyzing specifically the data from the WSS Trail (Table 1), it was possible to identify that the values of total porosity (%) and organic matter (%) were, without exception, lower on the bed of the trail compared to those obtained on the edges of the trail, given that the average values of these attributes were, respectively, 39% and 0.6% on the bed and 59% and 2.4% on the edges (Table 1).On the other hand, the values of bulk density (g cm −3 ) and pH were, also without exception, higher on the bed of the trail and lower at the edges, so the average values of these attributes were, respectively, 1.1 g cm −3 and 5.3 on the edges and 1.6 g cm −3 and 5.5 on the beds of areas degraded by erosion processes along the trail (Table 1).Regarding the values of soil granulometry, it was possible to identify the domain of the sand fractions, especially fine sand and silt, on the bed of the trail, since the average values were 57, 10, and 24%, respectively.Those on the edges of 57, 6, and 19%, respectively (Table 1).The only exceptions occurred at site A, where the fine sand content was the same in both the soil sample from the edge and the soil from the trail bed (8%), as well as in the silt content, which was lower in the trail bed in comparison with that obtained at the edge (16 and 19%, respectively) (Table 1).It was also possible to identify that the clay contents in the border soil samples were higher than those obtained on the beds, with average values of 24 and 19%, respectively (Table 1).The only exception was at site A, where the clay content was higher in the trail bed samples compared to those obtained at the edges (24 and 14%, respectively) (Table 1).This granulometric composition, rich in sand, indicates a textural classification of soils that vary between sandy loam and sandy clay loam (Table 1).
Regarding the data obtained on the CDN Trail (Table 2) regarding the soil's physical and chemical characteristics, it was possible to identify that the values of total porosity (%) and organic matter (%)were lower on the trail bed compared to those that were obtained at the edges, as was the case on the WSS Trail.This is because on the bed, the average values of these attributes were 33% and 1.0%, respectively, and on the edges they were 59% and 1.4%, respectively (Table 2).The values of bulk density (g cm −3 ) and pH were higher on the trail bed and lower at the edges, so the average values of these attributes, respectively, Soil Syst.2024, 8, 24 9 of 18 were 1.7 g cm −3 and 5.6 on the beds of the trails and 1.1 g cm −3 and 5.4 in the soil along the edges of the trails (Table 2).Regarding the granulometric analysis of the soils, it was noted that, on the bed of the trail, there was a predominance of fractions of silt and sand, especially fine sand.The average values of these fractions were 22%, 58%, and 33%, respectively (Table 2), whereas, at the edges, the values of these granulometric fractions were 21%, 46%, and 24%, respectively (Table 2).It is still possible to observe that clay contents were, without exception, higher on the border than on the bed of the trail; the average values of these attributes were 33% and 20%, respectively (Table 2).These compositions of the granulometric fractions of the soils in this trail indicate the predominance of sand fractions; thus, with the exception of the soil on the edge of site A, the other soils were classified between sandy loam and sandy clay (Table 2).

Spatiotemporal Dynamics of Soil Microtopography
Soil microtopography monitoring of degraded sites with erosive features along the WSS Trail and the CDN Trail (Figures 9-13) indicated a space-time evolution of erosion processes, as well as deposition of organic and inorganic materials on the beds of these trails.This may be associated with intensive trampling and soil compaction, removal and deposition of vegetation cover and soil organic matter, and the rainfall indices recorded in the region.The evolution of the cross-sectional area at the monitored site was evident when observing that the soil surface from the last monitoring (black line) was in a lower position than the soil surface during all other periods (orange, blue, and yellow lines) at all monitored sites (Figures 9-13).
period was ~4 cm greater than that in the first period (Figure 9).This space-time evolution corroborates the increase in the total cross-sectional area between the monitoring periods, since this area was 0.265 m 2 in the first monitoring period (2021), 0.284 m 2 in the second period (2022/1), 0.295 m 2 in the third period (2022/2), and 0.304 m 2 in the last period (2023).Therefore, this indicates an increase of approximately 0.039 m 2 in the cross-sectional area from the beginning of monitoring.Figure 10, which shows the soil microtopography at site B of the WSS Trail, also indicates a space-time evolution of the monitoring period erosion feature.In the central portion of this rill, for example, it was possible to identify that in the first monitoring period (2021), the maximum soil depth surface was ~65 cm, while in the last monitoring period (2023), this depth was greater than 70 cm (Figure 10).This space-time evolution was reinforced by the increase in the total cross-sectional area between the monitoring periods, since the total area was 0.327 m 2 in the first monitoring period (2021), 0.336 m 2 in the second period (2022/1), 0.258 m 2 in the third period (2022/2), and 0.372 m 2 in the last one (2023), i.e., an increase of about 0.044 m 2 in the total cross-sectional area from the beginning of the monitoring period (Figure 10). Figure 11, which shows the microtopography of the soil at site C of the WSS Trail, also shows a space-time evolution of the monitored erosion feature, especially in the lateral portions of the cross section (Figure 11).The exception was found in the central portion of this erosive feature.This occurred because the soil loss was severe at this site, and, in turn, reached the rocky layer underlying the soil layer, as can be seen in Figure 7B.Although the central portion of this rill did not indicate a space-time evolution at this point due to the rock outcrop, which is more resistant to the erosion process, it was possible to identify lateral soil loss in the other portions of this cross section.We considered the total cross-sectional area at this site, which was 0.177 m 2 in the first monitoring period, 0.192 m 2 in the second, 0.205 m 2 in the third, and 0.210 m 2 in the last, that is, an increase of 0.032 m 2 from the beginning of the monitoring period (Figure 11). Figure 11, which shows the microtopography of the soil at site C of the WSS Trail, also shows a space-time evolution of the monitored erosion feature, especially in the lateral portions of the cross section (Figure 11).The exception was found in the central portion of this erosive feature.This occurred because the soil loss was severe at this site, and, in turn, reached the rocky layer underlying the soil layer, as can be seen in Figure 7B.Although the central portion of this rill did not indicate a space-time evolution at this point due to the rock outcrop, which is more resistant to the erosion process, it was possible to identify lateral soil loss in the other portions of this cross section.We considered the total cross-sectional area at this site, which was 0.177 m 2 in the first monitoring period, 0.192 m 2 in the second, 0.205 m 2 in the third, and 0.210 m 2 in the last, that is, an increase of 0.032 m 2 from the beginning of the monitoring period (Figure 11). Figure 12 indicates the data associated with the monitoring of soil microtopography at site A of the CDN Trail.It is possible to identify that, in the central portion of this cross section, that is, inside the erosion feature, the maximum depth during the first monitoring period (2021) was approximately 14 cm.In the last period, this depth was approximately 4 cm greater in relation to the first.This space-time evolution is corroborated by data on the total area of the cross-section, which was 0.336 m 2 in the first monitoring period, 0.355 m 2 in the second, 0.37 m 2 in the third, and 0.381 m 2 in the fourth.This indicates an increase of approximately 0.045 m 2 in the total area from the beginning of the monitoring period.Figure 13, which presents the soil microtopography data from site B of the CDN Trail, also shows a space-time evolution of the monitored erosion feature that is similar to the behavior of the other monitoring sites,.At this site, it was possible to identify a significant space-time evolution of soil loss by exposing the roots, the greater distance between the lines that indicate the soil surface, or the increase in the total value of the cross-sectional area (Figure 13).In the first monitoring period, for example, the total cross-sectional area was approximately 0.259 m 2 , while in the second monitoring period, it was 0.284 m 2 .In the third monitoring period, the total value of this area was 0.307 m 2 , and in the last monitored period, 0.318 m 2 .This indicates an increase of about 0.059 m 2 from the beginning of the first monitoring period.

Discussion
When jointly considering the obtained data, it was possible to identify that these factors are related to each other, since certain physical-chemical parameters, for example, when subjected to certain dynamics of use and precipitation, tend to condition the evolution of erosive features such as those that were evidenced during the soil microtopography monitoring period (Figures 9-13) [9,14,17,18].
Knowing that water erosion is a process of mobilization, transport, and deposition of soil sediment through the action of rainwater, the concentration of rainfall in the PNSB region, especially in specific periods, tends to intensify aggregate breakdown and the detachment of soil particles from the bed of the trails, corroborating the evolution of the  Figure 13, which presents the soil microtopography data from site B of the CDN Trail, also shows a space-time evolution of the monitored erosion feature that is similar to the behavior of the other monitoring sites,.At this site, it was possible to identify a significant space-time evolution of soil loss by exposing the roots, the greater distance between the lines that indicate the soil surface, or the increase in the total value of the cross-sectional area (Figure 13).In the first monitoring period, for example, the total cross-sectional area was approximately 0.259 m 2 , while in the second monitoring period, it was 0.284 m 2 .In the third monitoring period, the total value of this area was 0.307 m 2 , and in the last monitored period, 0.318 m 2 .This indicates an increase of about 0.059 m 2 from the beginning of the first monitoring period.

Discussion
When jointly considering the obtained data, it was possible to identify that these factors are related to each other, since certain physical-chemical parameters, for example, when subjected to certain dynamics of use and precipitation, tend to condition the evolution of erosive features such as those that were evidenced during the soil microtopography monitoring period (Figures 9-13) [9,14,17,18].
Knowing that water erosion is a process of mobilization, transport, and deposition of soil sediment through the action of rainwater, the concentration of rainfall in the PNSB region, especially in specific periods, tends to intensify aggregate breakdown and the detachment of soil particles from the bed of the trails, corroborating the evolution of the Figure 9 shows that, during the first monitoring period of this site (2021), for example, the maximum depth of this rill was approximately 25 cm, while in the last monitoring period (2023), the depth increased by ~8 cm.Furthermore, it should be noted that the width of this rill also showed space-time evolution, since the value in the last monitoring period was ~4 cm greater than that in the first period (Figure 9).This space-time evolution corroborates the increase in the total cross-sectional area between the monitoring periods, since this area was 0.265 m 2 in the first monitoring period (2021), 0.284 m 2 in the second period (2022/1), 0.295 m 2 in the third period (2022/2), and 0.304 m 2 in the last period (2023).Therefore, this indicates an increase of approximately 0.039 m 2 in the cross-sectional area from the beginning of monitoring.
Figure 10, which shows the soil microtopography at site B of the WSS Trail, also indicates a space-time evolution of the monitoring period erosion feature.In the central portion of this rill, for example, it was possible to identify that in the first monitoring period (2021), the maximum soil depth surface was ~65 cm, while in the last monitoring period (2023), this depth was greater than 70 cm (Figure 10).This space-time evolution was reinforced by the increase in the total cross-sectional area between the monitoring periods, since the total area was 0.327 m 2 in the first monitoring period (2021), 0.336 m 2 in the second period (2022/1), 0.258 m 2 in the third period (2022/2), and 0.372 m 2 in the last one (2023), i.e., an increase of about 0.044 m 2 in the total cross-sectional area from the beginning of the monitoring period (Figure 10). Figure 11, which shows the microtopography of the soil at site C of the WSS Trail, also shows a space-time evolution of the monitored erosion feature, especially in the lateral portions of the cross section (Figure 11).The exception was found in the central portion of this erosive feature.This occurred because the soil loss was severe at this site, and, in turn, reached the rocky layer underlying the soil layer, as can be seen in Figure 7B.Although the central portion of this rill did not indicate a space-time evolution at this point due to the rock outcrop, which is more resistant to the erosion process, it was possible to identify lateral soil loss in the other portions of this cross section.We considered the total cross-sectional area at this site, which was 0.177 m 2 in the first monitoring period, 0.192 m 2 in the second, 0.205 m 2 in the third, and 0.210 m 2 in the last, that is, an increase of 0.032 m 2 from the beginning of the monitoring period (Figure 11).
Figure 12 indicates the data associated with the monitoring of soil microtopography at site A of the CDN Trail.It is possible to identify that, in the central portion of this cross section, that is, inside the erosion feature, the maximum depth during the first monitoring period (2021) was approximately 14 cm.In the last period, this depth was approximately 4 cm greater in relation to the first.This space-time evolution is corroborated by data on the total area of the cross-section, which was 0.336 m 2 in the first monitoring period, 0.355 m 2 in the second, 0.37 m 2 in the third, and 0.381 m 2 in the fourth.This indicates an increase of approximately 0.045 m 2 in the total area from the beginning of the monitoring period.
Figure 13, which presents the soil microtopography data from site B of the CDN Trail, also shows a space-time evolution of the monitored erosion feature that is similar to the behavior of the other monitoring sites,.At this site, it was possible to identify a significant space-time evolution of soil loss by exposing the roots, the greater distance between the lines that indicate the soil surface, or the increase in the total value of the cross-sectional area (Figure 13).In the first monitoring period, for example, the total cross-sectional area was approximately 0.259 m 2 , while in the second monitoring period, it was 0.284 m 2 .In the third monitoring period, the total value of this area was 0.307 m 2 , and in the last monitored period, 0.318 m 2 .This indicates an increase of about 0.059 m 2 from the beginning of the first monitoring period.

Discussion
When jointly considering the obtained data, it was possible to identify that these factors are related to each other, since certain physical-chemical parameters, for example, when subjected to certain dynamics of use and precipitation, tend to condition the evolution of erosive features such as those that were evidenced during the soil microtopography monitoring period (Figures 9-13) [9,14,17,18].
Knowing that water erosion is a process of mobilization, transport, and deposition of soil sediment through the action of rainwater, the concentration of rainfall in the PNSB region, especially in specific periods, tends to intensify aggregate breakdown and the detachment of soil particles from the bed of the trails, corroborating the evolution of the erosive features that were identified (Figures 9-13).In this sense, the impact of the kinetic energy of raindrops when they reach the soil surface causes its sealing, compaction, and detachment through splash erosion.Consequently, this material is lost and transported through runoff, therefore leading to the space-time evolution of the erosion features monitored.This culminates in an increase in the area of the soil surface and its crosssections, as observed during the monitoring periods [9,14,17,18].
It is worth mentioning that Rangel and Guerra [22], as well as Lima et al. [9], have already identified that the evolution of cross-sectional areas on trails during monitoring periods of soil microtopography can also be associated with the presence or removal of organic material/litter by surface runoff.Therefore, this evolution is not strictly associated with the loss of mineral material from the soil.This behavior reinforces a possible influence of physical-chemical characteristics on this evolution, given that the organic matter content, for example, can accentuate or delay the evolution of the erosion processes.
Considering that the beds of the trails are the places with the most trampling and least vegetation cover compared to their edges, the values of bulk density, total porosity, organic matter, and pH, for example, tended to be different from those of soil samples that were collected on the edges, as evidenced in the data in Tables 1 and 2. In the soils on the trail beds, the minimal vegetation cover and the intensive trampling of the soil tend to contribute to the decrease in porosity values and organic matter content, in addition to favoring the increase in bulk density values and soil loss.This explains not only the evolution of the monitored erosion features (Figures 9-13), but also the values of porosity, organic matter, and bulk density (Tables 1 and 2).At the trail edges, the presence of vegetation cover and low trampling tend to reduce bulk density and affect pH values, while also directly affecting total porosity values and organic matter content (Tables 1 and 2).
The lack of vegetation cover on the trail beds also contributes to the dynamics of this behavior, in addition to exposing the soil to the impact of the kinetic energy of raindrops (splash), which in turn contributes to its compaction.Its respective loss also influences the levels of organic matter [8,[33][34][35][36].
As the lack of soil trampling on the edges favors a greater development of vegetation compared to the beds of the trails, the mechanical action of roots during the growth and retraction processes increases the amount of empty space inside the soil that can be filled by air and water, that is, the value of total soil porosity [35][36][37].In addition, because it is composed predominantly of plant and animal residues, organic matter provides nutrients for the endopedonic fauna, which, together with the agents secreted by the roots, favor an increase in microbial activity and total soil porosity [33,36,38].
The soil's characteristic chemical conditions influence the formation of aggregates, retention, and water infiltration.They also increase aggregate stability, the number of pores, and biological activity, in addition to reducing the density [8,14,[38][39][40][41].By acting as a cementing agent that unites soil particles, organic matter forms organic substances that, when released and made available to the soil, bind and agglutinate the particles, increasing aggregate stability and reducing their erodibility [39,42].Therefore, the ability of organic matter to affect the structural integrity of the soil, form aggregates, and stimulate the development of its biota favors aeration and increased permeability, which in turn increases water infiltration, reduces runoff, and increases soil resistance to erosion.Its absence, on the other hand, tends to trigger water erosion processes, remove soil material, and both cause and reduce nutrient retention [8,41,43,44].
In addition to vegetation cover influencing the disposition of organic matter, it also influences the total soil porosity, which in turn controls percolation, infiltration of water inside the soil, the penetration of air, and the movement of roots.The modifications to the distribution, volumes, sizes, and shapes of these pores through use and soil cover, for example, affect permeability, aeration, hydraulic conductivity, and, consequently, the dynamics of water erosion processes [8,[44][45][46].
Considering that porosity values are inversely proportional to their density, and that intensive trampling by people contributes to reducing soil porosity, this behavior implies that the soils on the bed of the trail are more susceptible to erosion compared to the soils on the edge, either due to low porosity values, low levels of organic matter, or high bulk density values.This justifies not only the values found in Tables 1 and 2, but also the evolution of soil microtopography (Figures 9-13).On the other hand, the edge soils, which are not trampled, have low erodibility due to the maximum values of porosity and organic matter and the low values of bulk density, thus favoring the dynamics of the values shown in Tables 1 and 2 [8,45].
Erodibility tends to increase with high rainfall in the PNSB region, which, in turn, tends to be concentrated in specific periods, especially during the beginning of summer and the holiday (vacation) period (December and January) (Figure 13).In addition to the absence of vegetation covering the low levels of organic matter in the trails' soils and, in turn, increasing the instability of its aggregates, it facilitates the rupture of its particles due to the impact of raindrops, thus favoring the formation of crusts in the topsoil layer, reducing water infiltration, and increasing the soil loss via runoff.Therefore, when concentrated in specific periods, the rain tends to intensify the compaction of soils due to trampling by people on the beds of the trails, favoring the formation of preferential paths of water that culminate in the development of erosion features.This influences not only the values of some physical-chemical parameters (Tables 1 and 2), but also the loss of organic and mineral particles, as evidenced in the monitoring period of soil microtopography with the evolution of the cross-sectional area (Figures 9-13).
The pH values were all acidic (pH < 7), with the highest values found at the beds of the trails and the lowest associated with the trail edges (Tables 1 and 2).The fact that the edges have greater vegetation cover and, consequently, higher levels of organic matter may indicate the excretion of certain acidic substances, either by the action of the roots or by the decomposition of vegetable and animal matter.Soils obtained on the beds of the trails, on the contrary, lack vegetation cover and materials for decomposition (Tables 1 and 2).The more acidic soils at the trail edges tend to present colloidal complexes deficient in chemical elements that confer greater stability, that is, they lead to reductions in erodibility and increase the shear resistance of soil particles, since some of these elements result not only in greater biological activity, but also in stability of the aggregates [47][48][49].
Soils on the edges, despite having lower pH values, are less susceptible to erosion processes because, in addition to not suffering from trampling or the direct impact of raindrops, these soils are more porous, less dense, less compacted, and have higher levels of organic matter, which, by producing and releasing humus into the soil, favor their aggregate stability.In contrast, the soils on the trails, despite having less acidic pH values than the soils on the edges, in addition to being subjected to intensive public use and the impact of raindrops, have higher bulk density values, low porosity, and low levels of organic matter.This justifies our findings (Tables 1 and 2), and the evolution of erosion was evidenced by monitoring the soil microtopography (Figures 9-13).
The granulometry and textural classification of the soils indicate the predominance of clayey, medium-texture soil (clay-sandy loam), in which clay and sand fractions predominate in the granulometric composition of the soils.The predominance of coarse sand fractions and clay fractions causes the soils to have low erodibility, considering that the coarse sand fractions, due to their diameters, weights, and fast decanting speeds, make their removal and transport difficult due to the action of water.The clay fractions, in turn, have high aggregation capacity, the colloids and specific surfaces with cohesive strength, and a significant presence of fractions of silt and fine sand, especially in soils on the beds of trails.These favor the formation and loss of soil through erosion, since these fractions are easily removed by low cohesion and insufficient weight in the face of detachment and transport caused by the action of water [50][51][52][53].
Further, these soils have granulometric compositions of solid particles with minimal cohesion, so these characteristics together can be fundamental conditioning factors for soil loss, which became evident due to the evolution of the area of the cross sections recorded during the monitoring of the soil microtopography (Figures 9-13).On the other hand, the soils on the edges had the highest values of porosity and organic matter, the presence of vegetation cover, low values of bulk density and pH, a lack of trampling, and greater presence of granulometric fractions of clay.These factors increase the cohesion of soil particles, hinder soil loss, and elucidate the data (Tables 1 and 2) [8,10,14,17,20].
In this sense, the joint action of soil physical and chemical characteristics, mainly the values of bulk density, total porosity, granulometry, and organic matter content associated with rainfall indices in the coastal region of the PNSB, as well as the intense public use of the trails of this protected area, were assessed.These conditions culminate in the development and evolution of erosion features along the trails, as evidenced in cross sections of the soil microtopography (Figures [9][10][11][12][13].This is because public use tends to influence some of these physical-chemical characteristics, while the dynamics of regional rainfall together with other intrinsic characteristics of the soil (granulometry) favor the evolution of erosion features.

Conclusions
The obtained data indicate the intensive trampling of the soil resulting from public use is associated with the intrinsic characteristics of the location, such as the granulometric fractions of high erodibility and concentrated rains, especially in periods that coincide with greater local visitation.These conditions are leading to the environmental degradation of the trails of this conservation unit, as evidenced by the evolution of the erosion features in the cross-sectional areas monitored during the survey.
Monitoring the soil microtopography has demonstrated a loss of solid organic and inorganic particles which culminates in a degradation process, given that in all the cross sections monitored, there were increases in the areas of erosion features from the beginning of the monitoring.Such behavior, therefore, evidences the surface runoff process, which has a high capacity to remove solid materials.In addition, the reduction in soil protection against the impact of raindrops is enhanced by the trampling of the ground resulting from the intense flow of people along the trails.
The physical-chemical characteristics of the soil also showed an association with the degradation process, given that compaction, evidenced by the values of bulk density, especially on the bed of the trails, contributes to the reduction in water infiltration and to the increase in surface runoff.This, in turn, culminates in soil loss along the trails.Organic matter content also contributes to the dynamics of this process; all results were higher at the edges of the trails compared to those obtained in the respective beds, and these contents affect soil erodibility.
Intrinsic pedological and climatic factors, such as soils with coarse texture and intense rainfall in concentrated periods, also contribute to the evolution of the monitored erosion processes.The presence of sand and silt fractions and low levels of clay in the soils from the bed of the trails, for example, contribute to soil loss, especially when accompanied by compacted soils and by concentrated rainfall in specific periods.
This collective evidence shows the need for management and recovery actions for degraded areas along these trails through strategic management actions and improvement of their physical-chemical attributes, whether through the incorporation of organic matter or the construction of steps.
Finally, it is concluded that the obtained data can be fundamental tools for identifying environmental weaknesses or potential, as well as supporting actions aimed at planning and managing trails for public use in conservation units, especially when these protected areas are in high demand for public use and have intrinsic environmental characteristics that favor the triggering of environmental degradation processes.

Figure 2 .
Figure 2. Location map of the study area on the Serra da Bocaina National Park (PNSB) coast.Author: Lima.

Figure 2 .
Figure 2. Location map of the study area on the Serra da Bocaina National Park (PNSB) coast.Author: Lima.

18 Figure 3 .
Figure 3. Topographic profiles of the trail to the Waterfall of Stone that Swallows (A) and the Caixa D'Aço Natural Pool Trail (B) on the coast of the PNSB.Author: Lima.Figure 3. Topographic profiles of the trail to the Waterfall of Stone that Swallows (A) and the Caixa D'Aço Natural Pool Trail (B) on the coast of the PNSB.Author: Lima.

Figure 3 .
Figure 3. Topographic profiles of the trail to the Waterfall of Stone that Swallows (A) and the Caixa D'Aço Natural Pool Trail (B) on the coast of the PNSB.Author: Lima.Figure 3. Topographic profiles of the trail to the Waterfall of Stone that Swallows (A) and the Caixa D'Aço Natural Pool Trail (B) on the coast of the PNSB.Author: Lima.

Figure 3 .
Figure 3. Topographic profiles of the trail to the Waterfall of Stone that Swallows (A) and the Caixa D'Aço Natural Pool Trail (B) on the coast of the PNSB.Author: Lima.

Figure 4 .
Figure 4. Geotouristic attractions on the trails studied on the PNSB coast: (left) Waterfall of Stone that Swallows and (right) Caixa D'Aço Natural Pool.Photos: Lima.

Figure 4 .
Figure 4. Geotouristic attractions on the trails studied on the PNSB coast: (left) Waterfall of Stone that Swallows and (right) Caixa D'Aço Natural Pool.Photos: Lima.

Figure 5 .
Figure 5. Soil collection sites A, B and C on the trail to the Waterfall of Stone that Swallows on the PNSB coast.The red arrow shows the direction of main flow during periods of intense rainfall.Photos: Lima.

Figure 6 .
Figure 6.Soil collection sites A and B on Caixa D'Aço Natural Pool Trail on the PNSB coast.The red arrow shows the direction of main flow during periods of intense rainfall.Photos: Lima.

Figure 5 .
Figure 5. Soil collection sites (A-C) on the trail to the Waterfall of Stone that Swallows on the PNSB coast.The red arrow shows the direction of main flow during periods of intense rainfall.Photos: Lima.

Figure 5 .
Figure 5. Soil collection sites A, B and C on the trail to the Waterfall of Stone that Swallows on the PNSB coast.The red arrow shows the direction of main flow during periods of intense rainfall.Photos: Lima.

Figure 6 .
Figure 6.Soil collection sites A and B on Caixa D'Aço Natural Pool Trail on the PNSB coast.The red arrow shows the direction of main flow during periods of intense rainfall.Photos: Lima.

Figure 6 .
Figure 6.Soil collection sites (A,B) on Caixa D'Aço Natural Pool Trail on the PNSB coast.The red arrow shows the direction of main flow during periods of intense rainfall.Photos: Lima.

Figure 7 .
Figure 7. Example of the PDE leveling on site A and B the trail to the Waterfall of Stone that Swallows on the PNSB coast.Photos: Lima.Figure 7. Example of the PDE leveling on sites (A,B) the trail to the Waterfall of Stone that Swallows on the PNSB coast.Photos: Lima.

Figure 7 .
Figure 7. Example of the PDE leveling on site A and B the trail to the Waterfall of Stone that Swallows on the PNSB coast.Photos: Lima.Figure 7. Example of the PDE leveling on sites (A,B) the trail to the Waterfall of Stone that Swallows on the PNSB coast.Photos: Lima.

Figure 8 .
Figure 8.Total monthly values of precipitation on PNSB coast during the four monitoring periods.

Figure 8 .
Figure 8.Total monthly values of precipitation on PNSB coast during the four monitoring periods.

Figure 9 .
Figure 9. Soil microtopography at site A of the Waterfall of Stone that Swallows Trail on the PNSB coast.

Figure 9 . 18 Figure 10 .
Figure 9. Soil microtopography at site A of the Waterfall of Stone that Swallows Trail on the PNSB coast.Soil Syst.2024, 8, x FOR PEER REVIEW 11 of 18

Figure 10 . 18 Figure 10 .
Figure 10.Soil microtopography at site B of the Waterfall of Stone that Swallows Trail, on the PNSB coast.

Figure 11 .
Figure 11.Soil microtopography at site C of the Waterfall of Stone that Swallows Trail on the PNSB coast.

Figure 11 .
Figure 11.Soil microtopography at site C of the Waterfall of Stone that Swallows Trail on the PNSB coast.

Figure 12 .
Figure 12.Soil microtopography at point A of the Caixa D'Aço Natural Pool Trail on the coast of the PNSB.

Figure 13 .
Figure 13.Soil microtopography at site B of the Caixa D'Aço Natural Pool Trail on the PNSB coast.

Figure 12 .
Figure 12.Soil microtopography at point A of the Caixa D'Aço Natural Pool Trail on the coast of the PNSB.

Figure 12 .
Figure 12.Soil microtopography at point A of the Caixa D'Aço Natural Pool Trail on the coast of the PNSB.

Figure 13 .
Figure 13.Soil microtopography at site B of the Caixa D'Aço Natural Pool Trail on the PNSB coast.

Figure 13 .
Figure 13.Soil microtopography at site B of the Caixa D'Aço Natural Pool Trail on the PNSB coast.

Table 1 .
Values of the physical-chemical characteristics of soils obtained from the Waterfall of Stone that Swallows Trail on the coast of the PNSB.

Table 2 .
Physical-chemical parameters of the soils obtained from the Caixa D'Aço Natural Pool Trail on the coast of the PNSB.