Search

Madagascar is known for its high erosion rates in the central highlands, yet the role of human disturbance versus natural processes is not well understood and is a topic of ongoing debate. At present the necessary quantitative data to couple vegetation dynamics and sediment fluxes over time in Madagascar is scarce. This study aims to provide more insight in vegetation changes and sediment transfers over the last millennia in the Lake Alaotra region, and specifically on the role of human disturbances and natural processes. Our vegetation reconstruction is based on pollen records from two lake sediment cores, covering the last 2600 years. Sediment accumulation rates were calculated from cores obtained from the floodplains, from wetlands surrounding the lake, and from Lake Alaotra itself. Our data show an early opening in the landscape, between 2050 and 1700 cal a BP, with a transition from a wooded grassland or woodland/grassland mosaic towards open grassland and an increase in charcoal accumulation rates. (Indirect) human impact is suggested as the main driver for these vegetation changes. Floodplain and wetland sediment accumulation rates only increase in the last 1000 years and peak in the last 400 years. This increased accumulation can mainly be linked to the increased anthropogenic pressure (grazing and farming activities) that triggered increased lavaka (gullies) activity. No changes in accumulation rate were observed in Lake Alaotra, indicating that most sediments are buffered in the floodplains and wetlands. Overall, our pollen and charcoal data suggest an indirect effect of human disturbance on vegetation shifts whilst strong evidence was found for a direct effect of human disturbance on sediment accumulation through intensified use of the grasslands.

Cosmogenic nuclide analysis of river sediment provides insight into erosion and catchment dynamics. Studies on factors controlling spatial variations in long-term erosion rates have often focussed on tectonically active mountainous areas, where strong linkages with topographic variables like catchment gradient and normalized river steepness have been found. Less is known about rates and controls in tropical areas with deeper soils in tectonically less active regions which are often intensively used by people. Information on long-term erosion rates is crucial if human impact on landscape dynamics is to be understood in these areas.Here, we investigate spatial patterns and controls on 10Be-inferred erosion rates in Madagascar, a moderately seismically active island surrounded by passive margins, with considerable relief and a climate that varies from humid tropical to semi-arid. We use a dataset of 99 detrital in-situ 10Be measurements from a wide range of catchments (combining new measurements with data from the literature), covering more than 30% of the country and including a wide range of topographic, bioclimatic, and geological characteristics. Overall, 10Be erosion rates are very low (2.4 - 51.1 mm kyr-1) but clear differences were found between different geomorphic regions with some of the highest rates on the eastern escarpment while most catchments in the central highlands had extremely low erosion rates. The latter shows that, under (sub-) tropical climax vegetation catchments can be very stable, despite a pronounced topography with convex slopes exceeding 30° and the presence of a thick and erodible regolith mantle covered by a protective laterite. Statistical tests indicate that 27% of the observed variation in 10Be erosion rates is associated with elevation (lower rates for higher catchments) and an additional 18% of the variation is associated with river concavity, seismicity, and lithological erodibility (higher rates for more convex, more seismically active and more erodible catchments). An additional test using random within-between (REWB) analysis, in which different geomorphic regions are also considered as independent variables, shows that the main variations in 10Be-inferred rates between the different regions are linked to river concavity, seismicity and gully abundance, where additional variation within geomorphic regions is linked with seismicity only. The random within-between model explained 73% of the observed variation, suggesting that differences between regions are indeed important, yet are only partly explained by the environmental controls we considered in our analysis. The fact that the association between topographical controls and 10Be-inferred erosion rates is weaker in Madagascar in comparison to tectonically (very) active areas is likely to be related to the long time scale considered and the fact that sediment buffering as well as individual random events may have a more important impact on 10Be-inferred erosion in relatively stable environments such as the ones we studied in Madagascar.

The extent to which the central highlands of Madagascar were once covered by forests is still a matter of debate: while reconstructing past environments is inherently difficult, the debate is further hampered by the fact that the evidence documenting land cover changes and their effects on carbon and sediment dynamics in Madagascar has hitherto mainly been derived from lake coring studies. Such studies provide an integrated view over relatively large areas but do not provide information on how land-use change affects hillslopes in terms of carbon and sediment dynamics. Such information would not only be complementary to lake inventories but may also help to correctly interpret lake sediment data. Carbon stable isotope ratios (δ13C) are particularly useful tracers to study the past dynamics of soil carbon over time spans ranging from years to millennia and thus to understand the consequences of land-use change over such time spans. We analysed soil profiles down to a depth of 2 m from pristine forests and grasslands in the Lake Alaotra region in central Madagascar. Along grassland hillslopes, soil organic carbon (SOC) content was low, from 0.4 % to 1.7 % in the top layer, and decreased rapidly to ca. 0.2 % below 100 cm depth. The current vegetation predominantly consists of C4 grasses (δ13C ∼ −13 ‰), yet topsoil δ13C-OC ranges between −23.0 ‰ and −15.8 ‰, and most profiles show a decrease in δ13C-OC with depth. This contrasts with our observations in the C3-dominated forest profiles, which show a typical profile whereby δ13C values increase slightly with depth. Moreover, the SOC stock of grasslands was ∼ 55.6 % lower than along the forested hillslopes for the upper 0–30 cm layer. δ13C values in grassland and forest profiles converge to similar values (within 2.0 ± 1.8 ‰) at depths below ∼ 80 cm, suggesting that the grasslands in the Lake Alaotra region have indeed developed on soils formerly covered by a tree vegetation dominated by C3 plants. We also observed that the percent of modern carbon (pMC) of the bulk OC in the top, middle and lower middle positions of grasslands was less than 85 % near the surface. This could reflect a combination of (i) the long residence time of forest OC in the soil, (ii) the slow replacement rate of grassland-derived OC (iii) and the substantial erosion of the top positions towards the valley position of grasslands. At the valley positions under grassland, the upper 80 cm contains higher amounts of recent grass-derived OC in comparison to the hillslope positions. This is likely to be related to the higher productivity of the grassland valleys (due to higher moisture and nutrient availability), and the deposition of OC that was eroded further upslope may also have contributed. The method we applied, which is based on the large difference in δ13C values between the two major photosynthetic pathways (C3 and C4) in (sub-)tropical terrestrial environments, provides a relatively straightforward approach to quantitatively determine changing vegetation cover, and we advocate for its broader application across Madagascar to better understand the island's vegetation history.

Over the past decades advanced technology has become available, revolutionizing the assessment of surface topography. At smaller scales (up to a few km²) structure from motion (SfM) algorithms applied to uncrewed aerial vehicle (UAV) imagery now allow sub-meter resolution. On the other hand, spaceborne digital elevation models (DEMs) are becoming increasingly accurate and are available at a global scale. Two recent spaceborne developments are the 12 m TanDEM-X and 30 m Copernicus DEMs. While sub-meter resolution UAV-SfM DEMs generally serve as a reference, their acquisition remains time-consuming and spatially constrained. However, some applications in geomorphology, such as the estimation of regional or national erosion quantities of specific landforms, require data over large areas. TanDEM-X and Copernicus data can be applied at such scales, but this raises the question of how much accuracy is lost because of the lower spatial resolution.Here, we evaluate the performance of the 12 m TanDEM-X DEM and the 30 m Copernicus DEM to i) estimate gully volumes, ii) establish an area-volume relationship, and iii) determine sediment mobilization rates, through comparison with a higher resolution (0.2 m) UAV-SfM DEM. We did this for six study areas in central Madagascar where lavaka (large gullies) are omnipresent and surface area changes over the period 1949-2010s are available. Copernicus derived lavaka volume estimates were systematically too low, indicating that the Copernicus DEM is not suitable to estimate erosion volumes for geomorphic features at the lavaka scale (100 – 105 m²). The relatively coarser resolution of the DEM prevents to accurately capture complex topography and smaller geomorphic features. Lavaka volumes obtained from the TanDEM-X DEM were similar to UAV-SfM volumes for the largest features, while smaller features were generally underestimated. To deal with this bias we introduce a breakpoint analysis to eliminate volume reconstructions that suffered from processing errors as evidenced by significant fractions of negative volumes. This elimination allowed the establishment of an area-volume relationship for the TanDEM-X data with fitted coefficients within the 95% confidence interval of the UAV-SfM relationship. Combined with surface area changes over the period 1949-2010s, our calibrated area-volume relationship enabled us to obtain lavaka mobilization rates ranging between 18 ± 3 and 311 ± 82 t ha-1 yr-1 for the six study areas, with an average of 108 ± 26 t ha-1 yr-1. This does not only show that the Malagasy highlands are currently rapidly eroding by lavaka, but also that lavaka erosion is spatially variable, requiring the assessment of a large area in order to obtain a meaningful estimate of the average erosion rate.With this study we demonstrate that medium-resolution global DEMs can be used to accurately estimate the volumes of gullies exceeding 800 m² in size, where the proposed breakpoint-method can be applied without requiring the availability of a higher resolution DEM. This might aid geomorphologists to quantify sediment mobilisation rates by highly variable processes such as gully erosion or landsliding at the regional scale, as illustrated by our first assessment of regional lavaka mobilization rates in the central highlands of Madagascar.

Over the past few decades, developments in remote sensing have resulted in an ever-growing availability of topographic information on a global scale. A recent development is TanDEM-X (TerraSAR-X add-on for digital elevation measurements), an interferometric synthetic aperture radar (SAR) mission of the Deutsche Zentrum für Luft- und Raumfahrt, providing near-global coverage and 12 m resolution digital elevation models (DEMs). Moreover, ongoing developments in uncrewed aerial vehicle (UAV) technology have enabled acquisitions of topographic information at a sub-meter resolution. Although UAV products are generally preferred for volume assessments of geomorphic features, their acquisition remains a time-consuming task and is spatially constrained. However, some applications in geomorphology, such as the estimation of regional or national erosion quantities of specific landforms, require data over large areas. TanDEM-X data can be applied at such scales, but this raises the question of how much accuracy is lost because of the lower spatial resolution. Here, we evaluated the performance of the 12 m TanDEM-X DEM to (i) estimate gully volume, (ii) establish an area–volume relationship, and (iii) determine mobilization rates through comparison with a higher-resolution (0.2 m) UAV structure-from-motion (SfM) DEM and a lower resolution (30 m) Copernicus DEM. We did this for six study areas in the Lake Alaotra region (central Madagascar), where lavaka (gullies) are omnipresent and surface area changes over the period 1949–2010s are available for 699 lavaka. Copernicus-derived lavaka volume estimates were systematically too low, indicating that the Copernicus DEM is too coarse to accurately estimate volumes of geomorphic features at the lavaka scale (100–105 m2). Lavaka volumes obtained from TanDEM-X were similar to UAV-SfM volumes for the largest features, whereas the volumes of smaller features were generally underestimated. To deal with this bias we introduce a breakpoint analysis to eliminate volume reconstructions that suffer from processing errors as evidenced by significant fractions of negative volumes. This elimination allowed the establishment of an area–volume relationship for the TanDEM-X data with fitted coefficients within the 95 % confidence interval of the UAV-SfM relationship. Our calibrated area–volume relationship enabled us to obtain large-scale lavaka mobilization rates ranging between 18 ± 3 and 311 ± 82 tha-1yr-1 for the six different study areas, with an average of 108 ± 26 tha-1yr-1 for the full dataset. These results indicate that current lavaka mobilization rates are 2 orders of magnitude higher than long-term erosion rates. With this study we demonstrate that the global TanDEM-X 12 m DEM can be used to accurately estimate volumes of gully-shaped features at the lavaka scale (100–105 m2), where the proposed breakpoint method can be applied without requiring the availability of a higher-resolution DEM. Furthermore, we use this information to make a first assessment of regional lavaka erosion rates in the central highlands of Madagascar.

Lavaka (gullies) are often considered as the prime indication of a currently ongoing human-induced environmental crisis in Madagascar's highlands. Yet, lavaka are known to have existed long before human arrival and account for the majority of the long-term sediment input into the highland rivers and floodplains. The role of anthropogenic disturbances in their formation therefore remains highly debated and it is unclear whether lavaka erosion has recently increased. Here,we address these questions by evaluating the dynamics of lavaka in the Lake Alaotra region (central Madagascar). An overall birth to stabilization ratio of 6.1 indicates a rapid lavaka population growth over the period 1949–2010s. Using data on lavaka development we calculated a mean lavaka population age of 410 ± 40 years and estimate that the disequilibrium started at 870 ± 430 cal. BP. Floodplain sedimentation starts to increase around 1000 cal. BP and peaks over the last 400 years, thereby independently confirming this time frame of increased lavaka activity. Lavaka population dynamics modelling shows that a strong increase in environmental pressure over the last centuries is needed to attain current disequilibrium levels. A general drying of the climate since 950 cal. BP in combination with the introduction of cattle and growing human presence around 1000 cal. BP will likely have triggered the increase in lavaka erosion. However, the recent acceleration cannot be explained by climatic changes alone and seems to be linked to increased anthropogenic pressure on the environment. As such, we offer a fresh and quantitatively supported perspective on lavaka dynamics and human impact in central Madagascar, where our methodology can be used in other locations where similar questions on geomorphic equilibrium need to be answered. ispartof: Science Of The Total Environment vol:806 issue:Pt 2 ispartof: location:Netherlands status: Published online

Over the past decades, developments in remote sensing have resulted in an ever growing availability of topographic information on a global scale. A recent development is TanDEM-X, an interferometric SAR mission of the Deutsche Zentrum für Luft- und Raumfahrt providing near-global coverage and high resolution DEMs. Moreover, ongoing developments in unmanned aerial vehicle (UAV) technology has enabled acquisitions of topographic information at a sub-meter resolution. Although UAV products are generally preferred for volume assessments of geomorphic features, their acquisition remains time-consuming and is spatially constrained. However, some applications in geomorphology such as the estimation of regional or national erosion quantities of specific landforms, require data over large areas. TanDEM-X data can be applied at such scales, but this raises the question of how much accuracy is lost because of the lower spatial resolution. Here, we evaluated the performance of the 12 m TanDEM-X DEM to i) estimate gully volumes, ii) establish an area-volume (A-V) relationship, and iii) determine mobilization rates, through comparison with a high resolution (0.2 m) UAV-SfM DEM and lower resolution (30 m) SRTM DEM. We did this for six study areas in the Lake Alaotra region (central Madagascar) where lavaka (gullies) are omnipresent and lavaka surface area changes over the period 1949–2010s are available for 699 lavaka. SRTM derived lavaka volume estimates were systematically too low, indicating that the SRTM DEM is too coarse to accurately estimate volumes of geomorphic features at the lavaka-scale (100–100 000 m2). Lavaka volumes obtained from TanDEM-X were similar to UAV-SfM volumes for the largest features, whereas the volumes of smaller features were generally underestimated. To deal with this bias we introduce a breakpoint analysis to eliminate volume reconstructions that suffer from processing errors as evidenced by significant fractions of negative volumes. This elimination allowed the establishment of an area-volume relationship for the TanDEM-X data that is within the 95 % confidence interval of the UAV-SfM A-V relationship. Our calibrated area-volume relationship enabled us to obtain large-scale lavaka mobilization rates ranging between 18 ± 6 and 289 ± 125 ton ha−1 yr−1 with an average of 102 ± 41 ton ha−1 yr−1. These results indicate that current lavaka mobilization rates are two orders of magnitude higher than long-term erosion rates. With this study we demonstrate that the global TanDEM-X 12 m DEM can be used to accurately estimate volumes of gully-shaped features at the lavaka-scale (100–100 000 m2), where the proposed breakpoint-method can be applied without requiring the availability of a higher resolution DEM. Furthermore, we use this information to make a first assessment of regional lavaka erosion rates in the central highlands of Madagascar.

Madagascar is an island characterized by a sharp vegetation gradients in landscape and the extent to which the central highlands were once covered by forest is still under debate. Stable carbon isotope ratios have been an important tracer to understand vegetation shifts in a landscape over time because forest plants (following the C3 photosynthetic pathway) and grasses (following the C4 pathway for the majority of tropical and subtropical species) show a different degree of isotope fractionation. The main aim of this study is to provide insight into the past vegetation history of the grasslands in the Lake Alaotra region (central Madagascar). This region is under high anthropogenic pressure as it is the most important rice-producing region in Madagascar. While much of the Lake Aloatra catchment consists of grassland-covered hills, with a high density of lavaka (large erosional features), pristine forests are located just east of the lake and these form the most western part of a larger rainforest-covered region extending up to the east coast of the island. Soil profiles were sampled along the hillslope gradient at both forested and grassland sites, whereby carbon stocks were quantified and δ13C values of soil organic carbon (OC) were measured to find evidence for past forest vegetation in the grassland sites. The soil organic carbon content of grassland soil profile was extremely low, from 0.4 to 1.8% in the top layer and rapidly decreasing to 0.2 % below 100 cm depth. The current vegetation predominantly consists of C4 grasses (δ13C ~-13 ‰), yet soil δ13C-OC ranges between -25.9 and -16.6‰, and most profiles show a decrease in δ13C-OC with depth, in contrast to observations in the (C3-dominated) forest profiles, which show a typical enrichment in 13C with depth. δ13C values in grassland and forest profiles converge to similar values (within 2.1 ±1.8 ‰) at depths below ~80cm, suggesting that the grasslands in the Lake Alaotra region have developed on formerly forested soils. Moreover, the soil OC stock of grasslands was ~55.6% lower than along the forested hillslopes for the upper 30 cm layer. Our results are consistent with the hypothesis that a vegetation change has occurred in the lake Alaotra region, and offer a promising avenue to expand this approach on a wider scale to help understand the vegetation cover change in the central highland of Madagascar.

To understand the driving processes for changes in African ecosystems and related problems such as soil degradation, it is crucial to gain insight in the relative importance of human disturbance and climate change. Madagascar is known for its particularly high erosion rates in the central highlands, yet the role of human disturbance versus natural processes are not well understood and is a topic of ongoing debate. Recent studies have challenged the traditional view that the currently observed intense erosion processes and sediment fluxes in Madagascar are mainly driven by recent large-scale deforestation. However, at present almost no quantitative data is available to couple vegetation dynamics and sediment fluxes over time in Madagascar. This study aims to provide more insight in landscape changes (vegetation changes, sediment mobilization and deposition) in central Madagascar, and in the specific role of man and climate. The study focuses on the 1800 km² catchment of Lake Alaotra, located ca 200 km northeast of Antananarivo. Lake Alaotra is formed in a graben system in the highlands of Madagascar, and is the largest freshwater lake of the country (400 km²). A pollen record from the lake was used to reconstruct regional vegetation changes. Radiocarbon dates of extracted pollen provide a detailed chronostratigraphic framework. Augerings and radiocarbon dates from floodplains and marshes in the catchment were used to reconstruct the sediment deposition history. The pollen record and charcoal data shows the vegetation changes over the last 3000 years. The main observed shift in vegetation is a transition from a woodland/grassland mosaic towards an open grassland, starting ca 1850 years ago, which coincides with the onset of human activities. Data on floodplain sedimentation show an increase in accumulation rates in the last 600 years, from ca 1 mm yr-1 to ca 30 mm yr-1, which can be linked to increased hillslope erosion rates during that time period. The sedimentation wave, however, does not reach Lake Alaotra nor the surrounding marshes as floodplains act as a buffer. Overall, this study provides a spatial and temporal integrated reconstruction of vegetation changes in the Lake Alaotra catchment and the link with sediment mobilization and deposition, thereby providing a better understanding of environmental changes in central Madagascar and its driving forces.

The Malagasy highlands are scattered with large inverse teardrop-shaped gullies called lavaka, which are by many considered as the prime indication of a currently ongoing human-induced environmental crisis. Yet, these gullies are known to have existed long before human arrival on the island, resulting in the highly debated role of anthropogenic disturbances on their formation. Here, we assess the dynamics of 700 lavaka in the lake Alaotra region from 1949 to the 2010s by using historical aerial pictures and present day satellite imagery. An overall birth to stabilization ratio of 6.1 indicates a currently rapid growing lavaka population. Observed growth-, birth- and stabilization rates allowed us to calculate a mean lavaka population age of 410 ± 40 years, and estimate that the current crisis started at 943 ± 430 cal. yr BP. This timeframe corresponds well with the “subsistence-shift”, where people move from hunting and foraging to farming and herding practices upon the introduction of cattle in the region. These findings were integrated into a novel, temporally explicit lavaka population model - building upon the observed lavaka growth-, birth- and stabilization rates and lavaka size distributions - where different environmental pressure scenarios were tested. Modelling outcomes show that the currently observed lavaka crisis largely results from a rapid increase in environmental pressure over the last centuries, likely caused by the combined effects of deforestation and overgrazing related to human population growth and the introduction of cattle. With this study we show the potential of an integrated data-modeling approach, where demographic concepts are applied to geomorphological features, allowing to link their evolution with past anthropogenically driven environmental changes.

Madagascar shows very high erosion rates and, in some regions, the landscape is dotted by major gullies called “lavaka”. Lavaka are also very frequent in the surroundings of Lake Alaotra, a large, shallow lake in the Malagasy highlands. A central question with respect to landscape evolution in Madagascar is to what extent human impact has triggered environmental change in terms of vegetation, erosion rates (lavaka formation) and sediment dynamics. Sedimentary archives in lakes such as lake Alaotra can be of great help to resolve this question provided that we understand how sediment and carbon are mobilised and transported through the landscape. In this study, we traced pathways of sediment and carbon fluxes through this eroding landscape, from the eroded hillslopes over various sediment deposition zones (floodplains, reservoirs, marshes..) to Lake Alaotra. Detailed profiles taken along convex hillslope transects (grasslands and primary forest), in the marsh peat, floodplains and lake were analyzed for carbon and nitrogen content, texture, and stable carbon isotope ratios (δ13C). Along the grassland hillslopes, soil OC content is extremely low, from 0.4 to 1.8% in the top layer and rapidly decreasing to 13C ~-13 ‰), yet soil δ13C ranges between -24 and -18‰, and most profiles show a decrease in δ13C with the depth – in contrast to observations in most C3-dominated systems. Contrary to our expectations, Lake Alaotra was found not to be a major sink of hillslope-derived sediments and/or carbon. Sediment cores from different parts of the lake have high OC contents (5 to 18%) and contain only minor amounts of sand, the dominant grain size class on the hillslopes. The high OC content of the lake sediments, in combination with data on C/N ratios and δ13C indicate that the OC in the lake sediments is mainly derived from the surrounding marshes and in situ primary production rather than from terrestrial C eroded from the catchment. Floodplains are likely a key sink for soil-derived sediments: similar to hillslope soils, sediment profiles in the floodplains show a low %OC and relatively high δ13C values ranging between -21 and -14‰.We conclude that most of the detritic sediments and carbon mobilised on the hillslopes through erosion do not reach lake Alaotra, even though erosion rates in the landscape are extremely high. Studying sedimentary profiles in the lake may provide information on environmental change (e.g. through changes in carbon contents and/or characteristics) but is insufficient to understand the entire sediment and carbon cascades in the Malagasy landscape.