Rain Forest Fragmentation and Environmental Dynamics on Nosy Be Island (NW Madagascar) at 1300 cal BP Is Attributable to Intensified Human Impact

Madagascar houses one of the Earth’s biologically richest, but also one of most endangered, terrestrial ecoregions. Although it is obvious that humans substantially altered the natural ecosystems during the past decades, the timing of arrival of early inhabitants on Madagascar as well as their environmental impact is still intensively debated. This research aims to study the beginning of early human impact on Malagasy natural ecosystems, specifically on Nosy Be island (NW Madagascar) by targeting the sedimentary archive of Lake Amparihibe, an ancient volcanic crater. Based on pollen, fungal spore, other non-pollen palynomorph, charcoal particle and diatom analyses combined with high-resolution sediment-physical and (in)organic geochemical data, paleoenvironmental dynamics during the past three millennia were reconstructed. Results indicate a major environmental change at ca. 1300 cal BP characterized by an abrupt development of grass (C4) dominated and fire disturbed landscape showing the alteration of natural rain forest. Further, increased soil erodibility is suggested by distinct increase in sediment accumulation rates, a strong pulse of nutrient input, higher water turbidity and contemporaneous increase in spores of mycorrhizal fungi. These parameters are interpreted to show a strong early anthropogenic transformation of the landscape from rain forest to open grassland. After ca. 1000 cal BP, fires remain frequent and vegetation is dominated by forest/grassland mosaic. While natural vegetation should be dominated by rain forest on Nosy Be, these last results indicate that human continuously impacted the landscapes surrounding the lake. At a local scale, our data support the “subsistence shift hypothesis” which proposed that population expansion with development of herding/farming altered the natural ecosystems. However, a precise regional synthesis is challenging, since high-resolution multi-proxy records from continuous sedimentary archives as well as records located further north and in the hinterland are still scarce in Madagascar. The lack of such regional synthesis also prevents precise comparison between different regions in Madagascar to detect potential (dis)similarities in climate dynamics, ecosystem responses and anthropogenic influences at the island’s scale during the (late) Holocene.


INTRODUCTION
Environmental change is a consequence of various kinds of causes. For one thing, natural factors are driving climate and environmental dynamics, however mostly on longer (millennial/orbital) time-scales. These natural, long-term dynamics are recently challenged by a global overprint of direct human impact on all compartments of nature, most obviously seen in increasing rates of loss in biodiversity. Also during past millennia, human impact has triggered environmental change and biodiversity loss, however only on local or regional scale. In this context the conservation of natural ecosystems and biodiversity is one of the greatest challenges of the 21st century (IPCC, 2021). The key to estimate future developments of ecosystems and to improve conservation strategies is, to study and disentangle past human impact and natural factors of climate and environmental dynamics. This is especially important in "closed" and "isolated" environments such as islands which generally harbor endemic and vulnerable biodiversity (Nogué et al., 2017(Nogué et al., , 2021Whittaker et al., 2017;Wood et al., 2017). In tropical to subtropical regions, islands and their ecosystems were frequently impacted by human colonization (Burney, 1997;Kahn et al., 2015;Gosling et al., 2017;Castilla-Beltrán et al., 2019;de Nascimento et al., 2020). Without written archives, paleorecords (e.g., from lake sediments, speleothems, peatbogs) are one of the best tools to assess human impact in the past (Battarbee and Bennion, 2011;Garcin et al., 2018a;Wang et al., 2019). However, the distinction between natural and human induced ecosystem shifts within these archives can be challenging. Human activities, such as deforestation or artificial fire, which lead to a loss of trees can easily be mixed up with natural effects like very dry climatic conditions (e.g., Clist et al., 2018;Garcin et al., 2018b;Maley et al., 2018;Bayon et al., 2019). Madagascar, characterized by one of the world's biologically richest ecosystems, does not deviate from this rule. Early human colonization and its impact on natural environmental changes of this island has been intensively debated for decades (Burney et al., 2003;Virah-Sawmy et al., 2010;Dewar et al., 2013;Anderson et al., 2018;Anderson, 2019;Godfrey et al., 2019;Crowley et al., 2021). Whereas initial settlement of Madagascar prior to 2000 cal BP still remains poorly resolved (Dewar et al., 2013;Anderson et al., 2018;Hansford et al., 2018;Douglass et al., 2019), recent studies on speleothems from northwestern Madagascar suggested that a major ecosystem shift, recorded at ca. 1250-1200 cal BP, was the result of an abrupt transition in human subsistence strategy (Burns et al., 2016;Voarintsoa et al., 2017;Railsback et al., 2020). The "subsistence shift hypothesis, " that consists to the transition from hunting/foraging to herding/farming, has been proposed by Godfrey et al. (2019) to explain the observed changes. Specifically, this transition would have affected habitats by strongly increasing fire frequency which triggered a rapid expansion of grassland at the expense of forest. In addition, population expansion and hence increased anthropogenic impacts, as well as introduction of livestock (Dewar and Wright, 1993) is assumed to have led to destruction of natural wildlife habitats and finally may have caused the extinction of Madagascar's megafauna (Crowley, 2010). Recent studies on the herbivore's diet support various competition factors between endemic and introduced fauna superimposed by environmental change as the main reasons for endemic extinctions (Hixon et al., 2021b,a). Despite these significant advances, disagreements still remain between studies from different regions, regarding the timing and intensity of changes, and the respective interrelation between humans and climate (e.g., Burney, 1987;Matsumoto and Burney, 1994;Virah-Sawmy et al., 2010;Anderson et al., 2018;Godfrey et al., 2019). In order to improve our understanding of this critical and complex transition, it is a necessity to increase spatial resolution of paleoenvironmental study sites to provide new data. Furthermore, the reanalyzes of previously studied sites with more detailed analyses are important as they can improve chronologies, add new proxies or increase the spatial or temporal data-resolution.
In this context, the present study focuses on Lake Amparihibe from Nosy Be island first studied by Burney (1999) and Burney et al. (2003). Previous results showed an increase in fire activity with a rise of Sporormiella at 1130 ± 50 14 C BP which corresponds to a calibrated age [Calib 8.20 (Stuiver and Reimer, 1993) using the SHCal20 curve (Hogg et al., 2020)] of 1000 +180 / −90 cal BP. This was interpreted as human transformation of the local landscape with the development of agriculture (e.g., rice) and first introduction of livestock (cattle, goats and sheep) on Nosy Be by Islamized Indian Ocean traders (Burney et al., 2003;Wright and Radimilahy, 2005). This site therefore represents a key record to illustrate human induced environmental change for the northwestern region of Madagascar and is frequently used as a benchmark for regional comparisons (e.g., Railsback et al., 2020). However, chronological control of this study is quite poor (two radiocarbon dates), and detailed study of sedimentological changes and past vegetation dynamics have not been published. To refine previous analyses and close these obvious gaps, Lake Amparihibe was once more targeted, especially as this site also exhibits excellent conditions to provide a high-resolution record of paleoenvironmental change. Lake Amparihibe represents a deep crater lake (ca. 50 m water depth) making it impervious against desiccation. Thus, it is anticipated that sediments within the lake have continuously recorded past environmental changes in a region where deep lakes are rare, particularly at low elevation. The location of the lake in the northwestern region, close to the coast, makes it also ideal to study initial human impact, since occupation and the establishment of settlements are anticipated to have occurred quite early in this region (Radimilahy, 1997;Dewar et al., 2013;Douglass et al., 2019).
Thus, pollen, fungal spores, other non-pollen palynomorphs (NPPs), macro-charcoal particles and diatom analyses combined with high-resolution sediment-physical and (in)organic geochemical data from Lake Amparihibe will provide an excellent base for comparisons with high-resolution paleoclimate records obtained from speleothems in northwestern Madagascar (Burns et al., 2016;Scroxton et al., 2017;Voarintsoa et al., 2017;Wang et al., 2019;Li et al., 2020;Railsback et al., 2020). Such multiproxy comparisons of different geoarchives are essential for better distinguishing natural from human induced ecosystem shifts in the past.

ENVIRONMENTAL SETTING
The small island of Nosy Be (320 km 2 ) is located approximately 10 km off the northwestern Madagascar coast ( Figure 1A). Nosy Be island is formed by a sedimentary bedrock pierced by volcanic vents and plutonic intrusions during the late Cenozoic and Quaternary (14-0.5 Ma; Cucciniello et al., 2016). This activity resulted in an important volcanic complex surrounding the peak of Mt. Passot (329 m a.s.l., western Nosy Be) with lava flows and tuff rings filled by nine freshwater lakes. The largest of these maar lakes, Lake Amparihibe (-13.32 • N, 48.21 • E, 71 m a.s.l., 163 ha, Figure 1B), a sub-circular crater of 1.5 km diameter with steep slopes at the lake shore and ca. 50 m water depth ( Figures 1B,C), is the target of this study. Its catchment with a size of 713 ha is drained by three major channels, two of them passing the two maar lakes east of the lake (Antsidihy and Maintimaso) and one originating from a third maar lake (Bemapaza). Located on altitudes of 236 m a.s.l. (Bemapaza), 94 m a.s.l. (Antsidihy) and 89 m a.s.l. (Maintimaso) all the lakes are connected hydraulically to Lake Amparihibe forming a lakecascade and serve as primary sediment traps for a big portion of the catchment (Figure 1B).
Climate is tropical sub-humid in northwestern Madagascar, the region surrounding Nosy Be and Ampasindava bay and it corresponds to the most humid area of the west coast of Madagascar. Mean annual temperature is around 25 • C without major seasonal variability. Due to its location, precipitation at Nosy Be is affected by the proximity to the ocean (Mozambique Channel, Indian Ocean) with mean annual precipitation ranging around 2000 mm (Hijmans et al., 2005). Precipitation rates show a strong seasonality, with a maximum during austral summer (November to April) including strong summer monsoon rainfalls accompanied by frequent tropical storms or cyclones, related to the southern position of the Inter Tropical Convergence Zone. During austral winter, rainfalls are less and generally related to trade winds.
With a dry season almost absent, the natural vegetation is a subhumid rain forest called Sambirano rain forest (Koechlin et al., 1974). This forest is part of the East-West rain forest corridor distributed between northwestern and northeastern Madagascar ( Figure 1A). Despite some differences in species composition compared to the forest developed in eastern Madagascar, the main plant families are similar (Koechlin et al., 1974). Among the main tree families are Euphorbiaceae, Rubiaceae, Araliaceae, Ebenaceae, Sapindaceae, Anacardiaceae, Elaeocarpaceae, Lauraceae, Myrtaceae and among the most represented shrub are Ochnaceae, Erythroxylaceae, Myrsinaceae and Celastraceae. Liana and epithetic plants are also well represented by Apocynaceae, Fabaceae, Acanthaceae and Combretaceae. Some families are highly diversified and include many different genera and species, among them are Urticaceae and Melastomataceae. However, due to human impact, natural vegetation on Nosy Be island is now characterized by a secondary complex vegetation including about 14% of rain forest, 41% of wooded savanna, 12% of grasslands and 33% of crops mosaics (Mayaux et al., 2000). Typical trees of secondary forest, such as Trema orientalis are now common in the remaining Sambirano rain forest (Rasoanaivo et al., 2015). Close to the study site, west of Lake Amparihibe, the vegetation is dominated by a wooded grassland-bushland mosaic. To the east, vegetation assemblage shows a succession following an altitude gradient up to the peak of Mt. Passot, with a wide spectrum of species from shrub bushland to rain forest (Koechlin et al., 1974;Moat, 2007).

Pre-site Survey/Hydro-Acoustics
Prior to sediment sampling a shallow seismic survey was conducted to map basin morphology, sediment thickness and structures. A parametric sediment echosounder (SES 2000light, Innomar Technologie) with a primary frequency of 100 kHz and a secondary frequency range of 4 -15 kHz, depending on the desired resolution and penetration depth was applied. Positioning was provided by GPS (Navilock NEO M8U Multi GNSS UDR Receiver). The echosounder was calibrated to a sound velocity of 1460 ms −1 for depth measurements. With these settings, a dense grid of altogether 25 km of profiles were acquired. Final processing of the data was done with the software ISE 2.95 (Innomar Technologie).

Sampling, Water Data, and Chronology
Fieldwork was carried out in November 2016 and based on the shallow seismic survey, two positions within the lake were identified as the most suitable for coring using a gravity coring system (06/2021 1 ). AMP16-1 (-13.323917 • ; 48.209889 • ) with a total length of 164.5 cm was taken at a water depth of 45 m, ∼300 m apart from AMP16-2 which is 168 cm long and obtained at the deepest part of the lake (-13.322106 • ; 48.211828 • , 50 m water depth) ( Figure 1B). Concurrently, vertical in situ FIGURE 1 | (A) Distribution map of natural biomes in Madagascar (Olson et al., 2001). (B) Close up topography of western Nosy Be including the catchment of Lake Amparihibe (AMP), lake bathymetry, coring positions (red stars) and the location of the seismic profile in the lake (purple line) shown in panel (C), as well as an archeological site (black asterisk; Wright and Radimilahy, 2005) (data source: DEM data derived from SRTM 1 arc-second (06/2021 -https://earthdata.nasa.gov/); bathymetric data based on interpolation of seismic profile lines; map creation with QGIS 3.1). (C) Typical seismic image of Lake Amparihibe. The uppermost ca. 1.5 m of the deposits show clearly layered sediments, further penetration is only possible within acoustic windows, which show well-layered sediments down to at least 10 m sediment depth. The limited penetration results from widely distributed gas in the sediment most likely due to decomposition of organic matter. A clear basement reflector was not detected.
water parameter measurements were conducted using a WTW MPP 3430 multi parameter probe for pH, dissolved oxygen, conductivity and temperature next to AMP 16-1 ( Figure 1B). Water transparency was assessed with a Secchi disk. Cores were stored under cool (4 • C) and dark conditions at the laboratory facilities (Physical Geography department, Friedrich-Schiller-University Jena, Germany) until further processing. Cores were split, photographed and described following standard protocols. Magnetic susceptibility was scanned in 2 mm steps with three replicate measurements using a MS2E surface scanning sensor (Bartington Instruments). Based on the initial description, both cores showed different sediment accumulation rates with AMP 16-1 obviously reaching farther back in time. However, since previous analyses of AMP 16-1 found hints for a reworked upper part of the core (at least uppermost 20 cm; Haberzettl et al., 2021; Supplementary Figure 1), the two cores were aligned and combined to a composite master sequence ("AMP-16") using lithological description, patterns of magnetic susceptibility, and specific marker layers (Figure 2 and Supplementary Figure 1). The most outstanding marker layer, expressed as a 6 mm thick orange-brown fine-grained layer, exhibiting the maximum value in magnetic susceptibility at 99.4 cm core depth in core AMP 16-1, as well as its equivalent layer at 156 cm depth in core AMP 16-2 (Supplementary Figure 1) was used as transition depth from AMP 16-2 (=upper part of the composite) to AMP 16-1 (=lower part of the composite).
Due to the lack of plant macro-remains, the age model relies on AMS radiocarbon dating of three bulk sediment samples from each core, performed at either Beta Analytic Inc. (Miami, FL, United States) or the Poznań Radiocarbon Laboratory (Poznan, Poland) ( Table 1). The age of the sediment surface is anticipated to be modern, and was thus set to the year of coring (AD 2016 = −66 cal BP). The age model was performed as a function the composite depth with the RStudio software with the R-package "Bacon" (V. 2.5.6) (Blaauw and Christen, 2011) using SHCal20 calibration curve (Hogg et al., 2020).

Sedimentological and Geochemical Analyses
A total of 46 samples were taken at intervals of 2 to 18 cm depending on sediment layering for grain size analyses. Sample aliquots of ∼1 g were treated with H 2 O 2 (10%, 30%, heated to 80 • C for 2 h) to remove organic matter. Subsequently samples were repeatedly washed with deionized water and centrifuged until pH-neutrality. Measurements were carried out using a laser diffraction particle size analyzer (Beckman Coulter LS 13320) applying the Fraunhofer optical model for equivalent particle size diameter calculation. Each sample was measured in seven runs until signal reproducibility. Statistical indices (median size, sorting, particle size fractions) were calculated using a modified version of the MS Excel R Macro Gradistat V4.5s (Blott and Pye, 2001).
A total of 28 samples (ground to particle size < 40 µm) were analyzed at a coarser resolution for Total Carbon (TC), which equals total organic carbon (TOC) since all samples were carbonate free (no reaction with HCl), as well as Total Nitrogen (TN), bulk δ 13 C and δ 15 N and biogenic silica (BiSi). For carbon and nitrogen analyses, as well as their stable isotopic signature aliquots of 15 mg were packed into tin-boats and analyzed by an element analyzer (Vario EL Cube, Elementar) coupled to an IRMS (Isoprime Vision, Elementar). Based on triplicate measurements precision is 0.21% for TC and 0.01% for TN. T(O)C and TN were used to calculate the molar C/N. Bulk organic δ 13 C and δ 15 N were measured against certified standards (L-Prolin, EDTA and USG65) and reported in standard δ notation ( ) against Vienna Pee Dee Belemnite (VPDB) and Air, respectively. Relative errors based on triplicate measurements are 0.05 for δ 13 C and 0.31 for δ 15 N.
The same set of ground samples was analyzed for biogenic silica following a modified protocol of Ohlendorf and Sturm (2008). Sample aliquots of 40 mg were dissolved using 10 ml NaOH (1 mol l −1 ) supported by 3 min ultrasonic application and subsequent heating at 100 • C for 2 h. Si was measured at sample aliquots (pH set between 1 and 3) using an inductively coupled plasma-optical emission spectrometer (ICP-OES 725-ES, Varian). To account for non-biogenic Si, Al was measured as well, and minerogenic Si was subtracted assuming a Si/Al ratio of 1/1 for dissolved clay minerals.
A total of 21 samples were taken for leaf-wax n-alkanes analyses at sampling intervals between 2 and 18 cm. 10 ml dichloromethane (DCM): methanol (MeOH) (9:1) was added to sample aliquots of 1 to 3.4 g dry ground sediments and extraction was carried out using an ultrasonic bath in three 15 min-cycles. The resulting total lipid extract was separated over aminopropyl pipette columns. The n-alkanes were eluted with 4 ml hexane.
Measurement was carried out using gas chromatography with flame ionization detection (GC-FID) (7890B-GC System, Agilent Technologies) equipped with a HP5MS column (30 m, 320 µm, 0.25 µm film thickness). For identification and quantification an external n-alkane standard (n-alkane mix nC 21 -nC 40 , Supelco) was used. n-Alkane concentration was calculated as sum of nC 21 to nC 35 . Contributions from aquatic plants can be estimated with the Paq index (C23 + C25)/(C23 + C25 + C29 + C31) following Ficken et al. (2000). The average chain length (ACL), for example based on C 27 , C 29 , C 31 and C 33 , can be used to characterize the input from the terrestrial higher plants, with longer chains indicating more input from grasses compared to deciduous trees (Poynter et al., 1989;Zech et al., 2010Zech et al., , 2013aSchäfer et al., 2016). We here also adopt an ACL' based on C 29 , C 31 and C 33 , as well as an ACL based only on C 31 and C 33 .

Paleoecological Analyses
Paleoecological analyses include palynological analyses [counting of pollen grains, fern spores and non-pollen palynomorphs (NPPs)], counting of sedimentary macro-charcoal particles (>125 µm) and counting of diatoms. All analyses were performed on the composite master sequence AMP-16 except for diatoms performed solely on core AMP 16-2.

Palynological Analyses
For pollen and spore extraction, 22 samples of 0.5 cm 3 samples at intervals between 5 and 10 cm were prepared following a standard chemical protocol (Faegri and Iversen, 1975;Daniau et al., 2019). Pollen grains and fern spores were counted using a light microscope (Zeiss) at ×400 magnification, and a minimum sum of 300 terrestrial pollen grains were counted for each sample. The references consulted for identification of pollen and spores included several atlases (Straka and Freidrich, 1989;Gosling et al., 2013;Schüler and Hemp, 2016;Rasoloarijao et al., 2019), the online African Pollen Database (06/2021 2 ) and the reference collections of University of Göttingen (06/2021 3 ) and ISEM (06/2021 4 ). Percentages of pollen and fern spores were computed on the terrestrial pollen sum and the pollen record was subdivided into significant zones by applying a constrained cluster analysis by sum of squares performed with terrestrial pollen (CONISS, Grimm, 1987). Fungal spores were counted on pollen slides at an interval varying between 10 and 30 cm. A minimum fungal spore sum of 100-300 was counted for each sample, except for the depth 96 cm where only 73 spores were found. Fungal spore and other NPP identification are based on NPP Image Database (Shumilovskikh L. S. et al., 2021). Statistical analyses were performed using the RStudio software (Version 1.2.1335) with the package "Rioja" (Version 0.9-21) (Juggins, 2015).

Sedimentary Macro-Charcoal Particles
For macro-charcoal counting a total of 220 contiguous sediment samples of 0.5 cm 3 were retrieved. The extraction of sedimentary FIGURE 2 | Transition between the sediment cores AMP 16 1 and 16-2, digital image, lithological sketch, age-depth model and sediment accumulation rate (SAR) of composite sequence AMP-16. The composite was established based on parallelization of visual marker layers as well as magnetic susceptibility patterns of the two sediment gravity cores AMP 16-1 and AMP 16-2 (Supplementary Figure 1). Age-depth relation is based on a Bayesian model (Blaauw and Christen, 2011) with medians and the 2σ error ranges of five AMS calibrated 14 C ages using the SHCal20 dataset (Hogg et al., 2020). Centimeter-scale SAR variability are an artifact that does not provide a centimeter-scale prediction of SAR but merely a range of possible SAR in different core sections. * 1 year of sampling; * 2 inverse age/outlier, rejected from age modeling; Radiocarbon dates were measured on total organic matter of bulk sediment. Age calibration is computed by running the R-package "Bacon" (V. 2.5.6) with RStudio software with the (Blaauw and Christen, 2011) and by using SHCal20 calibration curve (Hogg et al., 2020).
macro-charcoal particles followed the standard protocol from Stevenson and Haberle (2005). Macro-charcoal particles were counted using a digital microscope at 5.0 magnification (Zeiss) and wooden and grass-based charcoal were distinguished following morphological description from Enache and Cumming (2006). Macro-charcoal concentration (particles cm −3 ) and macro-charcoal influx (particles cm −2 year −1 ) were computed for each sample.

High Resolution Pollen and Macro-Charcoal Particles
Counting From Core AMP16-1 To reconstruct the ecological shift occurring at the transition between the two sediment cores with high precision, Poaceae pollen and macro-charcoal particle were counted on 15 contiguous additional samples from core AMP16-1 that encompasses this transition between the depth 110 and 96 cm. These detailed analyses performed every centimeter (for pollen and macro-charcoals) allow to trace environmental changes with high precision, without any potential disturbance from the switch between the two sediment cores in the composite master sequence AMP16 (Figure 2 and Supplementary Figure 1). For these contiguous samples, the same previously mentioned chemical treatment for pollen and macro-charcoal particle extraction were used. To obtain Poaceae percentages, a minimum of 100 terrestrial pollen grains were counted among which only Poaceae was identified.

Diatom Analyses
Diatom analyses were performed every 6-10 cm. Sample preparation and concentration assessment were conducted following standard procedures (Kalbe and Werner, 1974;Battarbee and Bennion, 2011). If possible, a minimum of 400 valves was identified in each sample. In samples with very low diatom concentration, analyses were stopped when 1000 microspheres had been counted. A simple preservation index was tested by counting at least 100 specimens of Aulacoseira granulata and differentiate between intact and broken valves (ratio I/B). Permanent slides for light microscopy were prepared using Naphrax R as a mounting medium. Slides analyses was performed using a Leica DM 5000 B light microscope with differential interference contrast, equipped with a ProgRes R CT5 camera, under oil immersion at ×1000 magnification. Identification is mainly based on (Krammer and Lange-Bertalot, 1986, 1988, 1991a, Houk and Klee (2007); Houk et al. (2010), Lange-Bertalot et al. (2017) and other relevant taxonomic publications (e.g., Klee et al., 2000). Stratigraphic diatom zones (DZ) were defined by hierarchical cluster analysis of taxa occurring with > 1% using PAST ver. 4.04 (Hammer et al., 2001). Down-core changes in diatom results were visualized using C2 software version 1.7.7. (Juggins, 2007).

Bathymetry and Seismic Survey
The basin morphology ( Figure 1B) is characterized by steep slopes, both on shore and under water. An asymmetric shape in W-E direction might result from a former mass movement at the western side of Lake Amparihibe. A rather shallow-water area is notable in the NE. It may, alternatively, represent an alluvial fan due to river inflow and deposition of sediments from the catchment or result from a fossil mass movement deposit, probably due to a slope failure from a crater wall of the lake located to north/east of Lake Amparihibe.
Seismic survey of sedimentary structures shows typically well layered and undisturbed deposits within the upper ca. 1.5 m followed by a reflector that does not allow further penetration ( Figure 1C). This behavior can be explained by natural gas within the sediment that prevents the acoustic waves from further penetration. Only rare acoustic windows containing less gas give further insight into the sedimentary structure and show that the deeper infill consists of well layered sediments as well, at least down to ca. 8 m sediment depth. This represents a minimum thickness of the sedimentary infill, as gas occurrence prevents further penetration here as well. Thus, a hard-rock basement has not been detected so far.

Lake Water Characteristics
Lake Amparihibe is alkaline (pH ∼8.7) and dissolved oxygen values of surface waters are about 7.9 mg L −1 (∼103%), while vertical water parameter measurements confirm anoxic bottom waters (Supplementary Figure 2). Water temperature is high and decreases with depth from 29.3 • C to 27.2 • C. Temperature depth profile shows that the lake was thermally stratified in November 2016. Electrical conductivity is around 245 µS cm −1 and Secchi depth about 2.5 m.

Sedimentology and Chronology
Sediment core parallelization using magnetic susceptibility patterns of both cores resulted in an almost perfect fit with a correlation coefficient of r = 0.83 (Supplementary Figure 1). According to this fit, the composite record has a total length of 164.5 cm. The lower part of the composite is represented by the lower section of AMP 16-1 from 164.5 to 99.4 cm, whereas the upper 99.4 cm are represented by the compressed core AMP 16-2 from 156.2 to 0 cm (Figure 2 and Supplementary Figure 1).
The composite record generally exhibits dark to light brownish sediments with various light grayish to reddish intercalated 1 mm-to 15 mm-thick layers. Sediments are mainly composed of clayey silts with only little portions of sand (max. 5.6%, see Supplementary Figure 3). From 164.5 to 99.5 cm the record shows general dark brownish colors and only few intercalated grayish and reddish faint layers. Six of the intercalated layers are outstanding since they are quite distinct, exhibit reddish color and show higher proportions of sand. Between 99.5 and 64 cm lamination is much more distinct, layer thickness varies between 1 and 5 mm, and the overall sediment color is brighter. Between 64 and 21.5 cm sediment color and texture differ from the lower parts. Color changes to light brown and particles are distinctively finer with higher proportions of clay. Aside of an indistinct lamination, an 8.7 cm thick grayish layer with a dark, coarse base a reddish center and a fining upward trend is outstanding from 50.8 to 42.1 cm composite depth. In the uppermost unit, the color changes to homogenous slightly darker brown.
The outstanding layers with a minimum thickness > 4 mm, showing both a particle size fining upward trend and peaks in magnetic susceptibility, were interpreted as Event Related Deposits (ERDs) (Figure 2). Assuming a quick (hours to days) sedimentation as turbidity currents, these ERD layers were omitted for age modeling. One of the dated samples shows an inverse age, since it was taken from the base of an ERD, and was thus rejected from age-modeling (Figure 2 and Table 1). The model is thus based on five radiocarbon ages in stratigraphic order and reveals a basal age of 2940 cal BP. This chronology refines a paleomagneticsupported age model previously published for AMP 16-1 (Haberzettl et al., 2021). Sediment Accumulation Rate (SAR) is variable ranging between 0.1 and 4 mm a −1 . At the lowermost part of the record (164.5 to 97 cm) SAR is low ranging between 0.1 and 0.6 mm a −1 . At 97 cm composite depth, SAR dramatically increases to values ranging around 1.3 mm a −1 and remains on this high level until 61.2 cm, showing some major peaks of up to 4 mm a −1 . Between 61.2 and 34.5 cm, a slight declining trend to 0.8 mm a −1 on average is notable. At 34.5 cm a distinct shift to low values ranging around 0.4 mm a −1 occurs. A very slight increment tendency is visible toward the top of the record with SARs up to 0.5 mm a −1 , while within the uppermost 0.5 cm of the record the value increases to 2 mm a −1 , which is likely an artifact of age modeling.
Based on lithological properties, the age-depth model, as well as the SAR, four sedimentary units (phases) were assigned, i.e., unit I from the base of the record to 99.5 cm (2940 -1300 cal BP), unit II from 99.5 to 64 cm (1300 to 1000 cal BP), unit III from 64 to 35 cm (1000 -820 cal BP), and unit IV from 35 cm to the top of the record (820 cal BP -today).

Magnetic Susceptibility
The composite magnetic susceptibility pattern shows values between 86 and 7035 × 10 −6 SI (Figure 3) with noticeable distinct narrow peaks matching the light grayish and reddish ERDs (Supplementary Figure 1). These peaks also represent the maximum values in the record. In unit I (2940 to 1300 cal BP) values are generally on a low level of around 600 × 10 −6 SI. Values drop down from around 1800 × 10 −6 SI in unit II (1300 to 1000 cal BP) to ca. 1000 × 10 −6 SI at 1000 cal BP in unit III (1000 to 820 cal BP). They distinctly increase at 890 cal BP to 3400 × 10 −6 SI and decline steadily thereafter until the top of unit III (820 cal BP) reflecting the uppermost ERD. Unit IV (820 cal BP -today) shows further declining magnetic susceptibility values with some internal variation to 800 × 10 −6 SI at 500 cal BP. After another broad peak centered at ca. 430 cal BP showing values of ca. 1600 × 10 −6 SI, the remaining record reveals fairly stable values ranging around 1100 × 10 −6 SI. However, the distinct minimum during the youngest part of the record (after 10 cal BP) is likely attributed to measurement edge effects.

Particle Size Characteristics
Grain size distributions in all samples are characterized as bito polymodal and poorly sorted (Figure 3 and Supplementary  Figure 3). Median grain size ranges between 4.9 and 16.4 µm. Unit I (2940 to 1300 cal BP) displays the occurrence of coarser particles between 6.6 and 16.4 µm, with a general fining trend. Particle size median (D50) is distinctively smaller 4.9 and 9.4 µm after 1300 cal BP and shows no grading trend and only minor variability. Particle sorting (σ), however generally described as poor, is quite variable over the entire covered time span. In unit I, deposits are very poorly sorted. From 1300 to 820 cal BP (unit II and III) sorting tends to be better, but highly variable and remains quite stable afterward in unit IV (820 cal BP -today). Skewness (S k ) in the distribution of all samples is positive, but comparable to the sorting showing internal variability. In the older deposits (unit I, 2940 to 1300 cal BP) S k reveals lower values with some variability and an increasing trend. In unit II and III, from 1300 to 820 cal BP, S k is on an intermediate level with high internal variability and remains on a high positive level toward the top of the record (unit IV).

Geochemical Characterization
T(O)C, TN as well as the C/N ratio follow the units of the AMP-16 record, but unit II and III (1300 to 820 cal BP) are merged (Figure 3). In general, TC varies between 2.1 and 13.4%, and TN shows contents from 0.2 to 1%. Stable, high amounts of TC and TN are observed within unit I (2940 to 1300 cal BP). In unit II (1300 to 1100 cal BP), TC and TN drop distinctly and remain low and relatively stable, particularly during unit IV (after 820 cal BP). The C/N ratio drops from relatively high values, partly exceeding 16 in unit II (at 1200 cal BP), to values below 14 in the units above. in unit IV (after 820 cal BP) (Figure 3). Stable values for δ 15 N of around 1.8 occur in unit I (2940 to 1300 cal BP) as well as in untis II and IV (after 1000 cal BP), but reveals a distinct enrichment of up to 4.3 in unit II (Figure 3). Biogenic silica (BiSi) contents range between 2.6 and 23.7% (Figure 3). The amount shows a decreasing trend during unit I and unit II until 1100 cal BP, when it reaches its minimum. Contents are dramatically increasing to the maximum thereafter and are again decreasing at ca. 700 cal BP, however remaining on a constant intermediate level of ca. 12% toward the top.
n-Alkanes occur in all investigated samples, and the chromatograms show the odd over even predominance typical for plant-derived leaf-waxes (e.g., Schäfer et al., 2016). Because of the low concentrations of the even n-alkanes and insufficient peak separation from other aliphatic compounds, we quantified only the more abundant odd n-alkanes. The sum of the odd n-alkanes (nC 21 to nC 35 ) reaches 20 µg g −1 dry sediment in unit I (2940 to 1300 cal BP), drops below 10 µg g −1 in unit II and III (1300 to 820 cal BP), and increases again over the last few centuries (Figure 4). The Paq is low at ca. 0.1 in unit I, increases to values > 0.2 in unit II and III, and exceeds values > 0.4 in unit IV. The ACL mirrors this trend, starting with values around 30.5, dropping in unit II and III, and reaching values < 29.5 in unit IV. Excluding C 27 , or even C 27 and C 29 , for the calculation of the ACL and ACL , respectively, of course increases the values, but it also reduces and even removes the trend toward shorter chain length. The remaining feature of the ACL is a pronounced maximum in unit II.

Palynological Analyses
A total of 93 different pollen and spore taxa and 89 NPP types have been identified. Based on the cluster analysis (CONISS) the pollen diagram is subdivided into three significant pollen zones (PZ, Figure 5).
PZ-2 (100-27 cm/1300-600 cal BP/units II, III and IV): AP diversity decreases as well as AP percentages that drop abruptly to 30%, then increase progressively up to 40% at the end of the PZ. Among AP, only percentages of Trema FIGURE 4 | Left plane: Leaf-wax n-alkane indices of composite record AMP-16. Sum of odd n-alkanes (nC21 to nC35), Paq, and average chain length ACL, ACL (nC39-33), and ACL (nC31,33). Right plane: mean n-alkane distribution patterns in respect to the sedimentological units (note that unit II and III were merged due to very similar distribution patterns).
increase slightly (up to 5-10%). Moraceae/Urticaceae (ca. 15-30%) and Macaranga/Mallotus (ca. 5%) still remain among the dominant AP. Poaceae pollen increase abruptly up to 60% then decrease progressively and reach 50% at the end of the PZ showing an opposite trend to AP. Other NAPs remain with low values. Ericaceae percentages show a slight increase (up to ca. 5%) and slight percentage decreases are recorded for Acalypha (<2%) and Begoniaceae (<2%). Percentages of aquatic pollen are characterized by an important increase in Potamogeton, a slight increase in Cyperaceae and Lemna is not recorded anymore. A steady occurrence of fern spores is recorded but with low values (< 3%) and the algae Botryococcus shows an increase at ca. 20% from ca. 1000 cal BP. Considerable changes in fungal spores are indicated by a sharp increase in Coniochaeta ligniaria associated with Delitschia, Cercophora and Sordaria. Spores of pyrophilous fungi Gelasinospora and Neurospora occur, while diversity of lignicolous assemblage decreases. Plant pathogens are presented by Lasioplodia, Tetraploa aristata, Dictyoarthrinium cf. sacchari and rusts uredospores.
PZ-3 (27-0 cm/from 600 cal BP/unit IV): The progressive increase in AP which started in PZ-2, continues in this zone with AP reaching 70%. AP are mainly dominated by Moraceae-Urticaceae (20-30%), Macaranga/Mallotus (5-10%) and Trema (20%). This zone is also characterized by the occurrence of new AP taxa such as Mimosa (up to 5%) or Myrtaceae (up to 3%) while several taxa characteristic of the first zone such as Podocarpus, Weinmannia, Ilex or Meliaceae/Sapotaceae are almost not recorded anymore. Among the NAP, only Poaceae pollen indicate significant changes characterized by a progressive decreasing trend from 45% to 25%. At ca. 10 cm depth an increase in the aquatic pollen Cyperaceae (up to 8%) is recorded. The other aquatic pollen, fern spores and Botryococcus do not show marked changes. PZ-3 is characterized by a further increase in diversity of coprophilous assemblages by Arnium, Apiosordaria and Podospora anserina-type. Chlamydospores Glomus-type occur constantly. Assemblage of lignicolous fungi change to Curvularia, Podospora curvispora, Canalisporium pulchrum, Megalohypha aqua-dulces, Caryospora, Zopfiella lundquistii. Plant pathogens diversity increase by several Spegazzinia species. Diversity of other spores decreases.

Sedimentary Macro-Charcoal Particles
Prior to 1300 cal BP, macro-charcoal particles occur in every sample, but concentrations and influx remain very low (Figures 6A,B). Following morphological descriptions from Enache and Cumming (2006), most macro-charcoal particles are wood-based. At ca. 1300 cal BP (unit II, 1300 to 1000 cal BP), a marked increase of concentration and influx is recorded with peaks exceeding 300 particles cm −3 and 40 particles cm −2 year −1 , respectively. Grass-based charcoal particles increase at the same time with percentages varying around 30% and peaks exceeding 40%. Subsequently, percentages decrease after 1000 cal BP and remain generally below 30%. From 1000 cal BP (unit III and IV, 1000 cal BP to today), macrocharcoal concentration and influx also tend to decrease. After a last peak above 300 particle cm −3 , between 820 and 750 cal BP, concentrations fluctuate around 100 particles cm −3 . Macrocharcoal influx indicates a progressive decreasing trend from FIGURE 5 | (A) Pollen and fern spore records of composite record AMP16 retrieved from Lake Amparihibe (Nosy Be, Madagascar) plotted on an age scale. Percentages of pollen and fern spores were calculated on the terrestrial pollen sum. Pollen zone (PZ) indicated by dashed lines were identified by a cluster analysis (CONISS, Grimm, 1987). The light-colored pollen curves correspond to exaggerated curves by a factor 5. (B) Fungal spore record plotted on an age scale of composite record AMP-16. Percentages are calculated based on total count of fungal spores. Gray shaded areas demarcate sedimentological units. 1000 cal BP and values stabilize from ca. 750 cal BP varying around 5 particles cm −2 year −1 . This decrease in influx is partly the result of the decrease in SAR (Figure 2).
High Resolution Data From Core AMP16-1 (110-96 Cm Depth; Ca. 1640 to 1200 Cal BP) Between 110 and 100 cm (ca. 1640-1300 cal BP) low percentages of Poaceae (<10%) and macro-charcoal (max. 16 particles cm −3 ) are recorded (Figure 6C). Above 100 cm (<1300 cal BP), Poaceae show a continuous increase and reach 40% at 97 cm (ca. 1200 cal BP). Macro-charcoal particles show two important increases above 100 cm depth. The first increase corresponds to a peak at 99 cm and maximum concentration (368 particles cm −3 , ca.1300 cal BP) is reached in the sample where Poaceae percentage shows the first marked increase. This macro-charcoal peak at 99 cm occurs contemporaneously to a peak in magnetic susceptibility. The second increase starts at 97 cm (ca. 1200 cal BP) and macro-charcoal concentration reaches ca. 580 particles cm −3 at 96 cm.

Initial Environmental Conditions Prior 1300 cal BP
Between 2940 and 1300 cal BP, sedimentological and vegetation data do not show major changes (Figures 3, 5A). The high dominance of tree taxa suggests a large area of forest. Most tree taxa that are described as typical of rain forest in modern and fossil pollen samples from Madagascar are well FIGURE 7 | Diatom record from core AMP 16-2 retrieved from Lake Amparihibe plotted on an age scale with relative abundances taxa (>1%), diatom concentration (DW: dry weight), preservation index (ratio intact to broken valves of Aulacoseira granulata) and abundances of benthic + tychoplanktonic and planktonic taxa. Samples with less than 15 valves counted were excluded from abundance plots. Diatom zone (DZ) indicated by dashed lines were identified by a hierarchical cluster analysis. The lighter colored diatom curves correspond to exaggerated curves by a factor 5. Hatched area indicates zone of lowest diatom concentration and poorer valve preservation. Gray shaded areas on the left demarcate sedimentological units (I to IV).
represented in this period (e.g., Moraceae/Urticaceae, Celtis, Elaeocarpaceae, Macaranga/Mallotus and Podocarpus; Burney, 1988;Gasse and Van Campo, 2001). This suggests that the dense and humid rain forest was the dominant vegetation prior to 1300 cal BP on Nosy Be. The low fire activity during the same time confirms the occurrence of rain forest developed under humid environmental conditions and low disturbance regime ( Figure 6A). Furthermore, the presence of a wide variety of lignicolous fungi indicate the availability of dead wood substrate in the lake catchment as well (Figure 5B).
Sedimentological and geochemical results are consistent with the dominance of rain forest prior 1300 cal BP. Relatively low SAR supports the assumption of dense vegetation stabilizing the topsoils in the catchment (Figure 2). The comparably coarse particle median (Figure 3) might either be indicative for lacking grasses and thus a missing "filter" for larger minerogenic particles flushing into the lake by inflowing water, or coarser particles might point to comparably moist conditions with higher amounts of precipitation and thus higher transport energy of the inflows (Conroy et al., 2008). High organic carbon contents, high C/Nratios, and high amounts of leaf-wax n-alkanes during this phase further suggest high productivity within the catchment and input of terrestrial organic matter (Meyers and Ishiwatari, 1993; Figures 3, 4). The depleted δ 13 C values of the bulk organic matter point to predominantly C 3 vegetation (Figure 3; Farquhar et al., 1989;Meyers, 1994;Diefendorf et al., 2010).
The n-alkane patterns show a dominance of nC 29 , nC 31 , and nC 33 , and the Paq is accordingly very low. This indicates a dominant input from higher terrestrial plants over input from aquatic sources (Eglinton and Hamilton, 1967;Zech et al., 2010). Distinguishing between input from grasses and herbs versus deciduous trees, as has been done in other studies (e.g., Zech et al., 2013b;Schäfer et al., 2016), needs to take into account that common plants of rain forests synthesize high amounts of nC 31 and nC 33 (Gnecco et al., 1989;Abas and Simoneit, 1998;Li et al., 2016;Teunissen van Manen et al., 2019). The high ACL values in comparison to northern mid-latitude locations (e.g., Europe) thus reflect the local rain forest signal, and should not be misinterpreted to indicate dominant input from grasses and herbs.

Ecosystem Shift From 1300 to 1000 Cal BP
From 1300 cal BP, most tree taxa decrease and some of them almost entirely disappear (e.g., Podocarpus, Weinmannia, Ilex, Meliaceae/Sapotaceae, Figure 5A). Inversely, an increase of herbs (mainly Poaceae) is observed. Together with an enrichment in δ 13 C of the bulk organic fraction further suggesting a shift from C 3 to more C 4 vegetation (Farquhar et al., 1989;Meyers, 1994;Diefendorf et al., 2010), this points to the development of a mosaic of grassland and forest (Figure 3). Combined with the increase in macro-charcoal mainly associated with increased grass-based charcoal particles (Figure 6A), our results suggest an abrupt development of grasses dominated and fire disturbed landscape and confirm the alteration of the rain forest. Highresolution paleoecological analyses during that critical transition (see Figure 6C) reveal that grassland developed within few decades and occurs simultaneously with increases of fires. The co-occurrence between the first peak of macro-charcoal particles and the peak of magnetic susceptibility indicate that fires were local to the catchment and led to the observed forest alteration and initiated the increase of soil erosion. Similar observations were made at Laguna Potrok Aike, Patagonia where increased frequency dependent magnetic susceptibility and increased charcoal accumulation rates indicated the destruction of the vegetation cover by fire leading to enhanced soil erosion (Haberzettl et al., 2006).
Vegetation changes are also evidenced by the considerable development of saprophilous/coprophilous fungal spores (Figure 5B), that generally develop on decaying plant remains (e.g., litter) and/or herbivore livestock dung (e.g., Krug et al., 2004;van Asperen et al., 2021). Plant parasitic fungi such as Tetraploa aristata, Spegazzinia sp., Dictyoarthrinium cf. sacchari indicate wide spread of herbaceous monocots and dicots as well while lignicolous fungi such as Xylariaceae decrease after 1100 cal BP . Therefore, this complete rearrangement of lignicolous and plant parasites fungal spores confirms an important forest dieback at the expense of grassland. In addition to vegetation changes, enhanced soil erosion is also indicated by Glomus-type (Figure 5B), representing chlamydospores of the arbuscular mycorrhizal fungi growing symbiotically on plant roots (van Geel et al., 1989).
Sediment physical and organic characteristics are altered at this transition and SAR is increasing by almost one order of magnitude, supporting the pronounced decrease in tree cover associated with increased vulnerability of the soils to erosion. This is also visible in the higher frequency of ERDs occurring within the deposits (Figure 2). Following a shortterm coarsening of the particles between 1300 and 1200 cal BP, the fining of the minerogenic particles accompanied by slightly better sorting and a more positive skewness of the distributions is observed at ca. 1150 cal BP (Figure 3). These changes suggest a selective deposition based on alteration of transport processes, i.e., lower amount of transport energy or shielding (McLaren, 1981;McLaren and Bowles, 1985). In this context, the development of Potamogeton from 1150 cal BP onward as an emergent aquatic plant, which likely covered the shallow (near shore/littoral) water areas, might have acted as sediment trap for larger particles. A lake level increase could have provided shallow water areas for macrophytes growth, since this would inundate the northwestern area of the catchment leading to the formation of a larger shallow water area, with emerging aquatic vegetation acting consequently as source of aquatic biomass and as sink for coarse minerogenic particles. A comparable process has been observed recently using aerial imagery (Supplementary  Figure 4). The lake level increase might also be associated with enhanced nutrient supply, which, in combination, most likely explain the rise of aquatic vegetation and algae which is also supported by a distinct shift to rather low C/N-ratios (Meyers and Ishiwatari, 1993). The exceptional enrichment in δ 15 N during this phase might support the enhanced nutrient supply through the input of isotopically enriched N with either topsoil (Ochiai et al., 2015), or distinctive input from anthropogenic/livestock sources (sewage) (Enters et al., 2010). Alternatively, denitrification of the system due to tenacious anoxic conditions at the lake floor might be an explanation for the high δ 15 N (Enters et al., 2010).
This increase of soil erosion impacted lake ecosystem dynamics by increasing nutrient discharges which probably contributed to the development of the aquatic plant Potamogeton followed by the development of the green algae Botryococcus. A change in lake ecosystem dynamic is further supported by the diatom record characterized by minimum concentration values between 1300 and 1040 cal BP (concentration minimum peaking at 1160 cal BP) and minimum valve preservation at ca. 1100 cal BP (Figure 7). Low concentration values are likely caused by enhanced catchment erosion and input of minerogenic material into the lake with associated dilution of autochthonous material. Additionally, the low valve preservation is suggested to be linked to both, mechanical breakage in a high dynamic (energy) environment, and potentially to an increase in pH (Flower and Ryves, 2009). A distinct shift in the dominating planktonic diatom species from Discostella cf. mascarenica to Aulacoseira granulata agg. further proofs a change in the lake ecosystem following the transitional phase. Discostella mascarenica is a rare species and only little is known on its distribution and ecological preferences. It was described from a small, alkaline, eutrophic freshwater pond on Réunion island, situated 800 km to the east of Madagascar (Klee et al., 2000). Aulacoseira granulata occurs worldwide in a wide range of trophic conditions but is often associated with more eutrophic waters. Moreover, it is frequent under turbulent conditions, capable to survive in higher turbid waters, has high growth demands for silica and has also been connected to phases of increased erosion (Kilham et al., 1986;Costa-Böddeker et al., 2012;Bicudo et al., 2016;Nardelli et al., 2016). Shifts among Discostella and Aulacoseira have been associated with physical restructuring of the water column, whereby small cyclotelloid species like Discostella benefit from a strongly stratified water column and larger-celled Aulacoseira taxa are favored by strong vertical mixing (Rühland et al., 2015). Thus, the change from Discotella cf. mascarenica to Aulacoseria granulata agg. in Lake Amaprihibe supports the assumption of a catchment alteration leading to increased turbulence and turbidity.
The distinctive decrease in leaf-wax n-alkane amount, the shift to shorter chain n-alkane predominance and the marked increase in Paq further supports the emergence of aquatic vegetation (Ficken et al., 2000;Aichner et al., 2010;Castañeda and Schouten, 2011) as the result of a change in nutrient supply to the lake. Decreasing ACL values document the input of shorter homologs. In order to evaluate and reduce a potential bias of aquatic and emergent plants (likely also producing C 27 and C 29 to some degree), the modified ACL and ACL is applied. The increase at 1300 cal BP in these modified indices point to enhanced input of long-chain C 33 likely originating from grasses. This would also support the increase in δ 13 C of the bulk organic fraction, which suggest a shift to higher amounts of C4 vegetation.

Development Toward Modern Vegetation Conditions
After 1000-900 cal BP, a slow and progressive forest recovery is evidenced by slight decrease of Poaceae pollen percentages and slight reduction of fire activity, while pioneer tree taxa (e.g., Macaranga/Mallotus) tend to increase (Figures 5A, 6). The coeval pronounced increase in biogenic silica (BiSi) until 820 cal BP (Figure 3), coinciding with high diatom concentrations, especially of the heavily silicified species Aulacoseira granulata, however indicates an enhanced nutrient input (Ferris and Lehman, 2007). The further lowering in C/N as well as the high Paq afterward also suggest a predominance of aquatic (macrophyte and algae) biomass production (Meyers, 1994(Meyers, , 2003 Figures 3, 4). Forest recovery and the decrease in fire occurrences continues after 600 cal BP, evidenced by still decreasing Poaceae pollen percentages, moderate high charcoal concentration and a slight depletion trend in the δ 13 C record toward C 3 vegetation endmember. However, pollen assemblages reveal a different tree composition than in rain forest prior to 1300 cal BP. Pioneer tree taxa characterized by fast recovery after disturbance (Trema, Macaranga/Mallotus) are more abundant, while several tree taxa do not recover (Podocarpus, Weinmannia, Ilex, Meliaceae/Sapotaceae or Schefflera) and new ones appear (Myrtaceae and Mimosa). This new tree assemblage points to the development toward modern vegetation conditions which are characterized by a secondary complex vegetation including rain forest, wooded savanna and grasslands with crops mosaics (Mayaux et al., 2000). The result of this substantial change of the landscape is also evidenced by a stable intermediate SAR and the deposition of constantly small minerogenic particles with a sorting and a distribution skewness on a stable but very different level, compared to the initial phase (unit I, Figures 2, 3). Magnetic susceptibility on a quite stable intermediate level suggests a constant input of soil material, which is likely originating from farming activities in the catchment of the lake.

Regional Implications in Northwestern Madagascar
Several studies proposed that human arrival in Madagascar occurred during the early or mid-Holocene (Dewar et al., 2013;Hansford et al., 2018). However, timing and modality of these early occupations are still debated. Except for cut marks on bones, no direct human evidence (archeological artifacts) has been described for the early Holocene (Hansford et al., 2018), while for the mid-Holocene, radiocarbon dating are controversial (Anderson et al., 2018;Anderson, 2019). The recent synthesis of Douglass et al. (2019) suggests that hunters/foragers populations inhabited the island at least since 2000 cal BP. Such activities generally consisting of ephemeral occupations are challenging to record with paleoecological archives since environmental impact is assumed to have been relatively low. Indeed, the vegetation surrounding Lake Amparihibe was mainly dominated by rain forest without major vegetation/environmental changes between 2940 and 1300 cal BP ( Figure 8A). Consequently, if humans were inhabiting Nosy Be prior to 1300 cal BP, human activity was not strong enough to alter the vegetation composition, or at least, not strong enough to be detected with the applied paleoenvironmental proxies.

Previous Paleoecological Study From Lake Amparihibe
The previous study conducted on Lake Amparihibe by Burney et al. (2003) sets the date for the first occurrence of livestock and initial human impact to 1000 +180 / −90 cal BP (calibrated age from original publication), by moderating the occurrence of Sporormiella and an increase in sedimentary macro-charcoal particles. Our results suggest a major transformation of rain forest vegetation, to an open, fire disturbed, landscape, but starting somewhat earlier at 1306 ± 70 cal BP ( Figure 8A). The rapid development of grasses (within several decades), coupled with peaks of macro-charcoal particles and saprotrophs/coprophilous fungal spores shows the first detectable human impact in northwestern Madagascar. These results illustrate the development of agriculture likely related to fire for pasture. Human impact has altered the environment and led to a major impact on the sediment accumulation increasing by one order of magnitude (0.13 to 1.3 mm a −1 ) within less than 100 years. The lake system reacted with changes in the trophic status on this sudden release of nutrients transferred with the eroded soils. Surprisingly, we were not able to confirm the predicted increase of Sporormiella in our record. After 1300 cal BP, a spore of Sporormiella has been counted only once on the top core while Burney et al. (2003) show that Sporormiella appear at the same time that macro-charcoal particles increase. Raper and Bush (2009) demonstrated that Sporormiella concentration declines from the shoreline across the first 100 m of a lake surface, highlighting the importance of the lake size and its changes in the interpretation of paleoecological records of Sporormiella abundance. Although both coring locations are in the central part of the lake, slight differences of distances from the shoreline might be an explanation since Sporormiella occur in low numbers (<2%) in Burney et al. (2003). However, evaluation of this effect in Nosy Be requires more research and other factors might also contribute to the different results (e.g., chemical treatments, spore identifications). Despite this difference, the entire assemblage of the saprotrophic/coprophilous fungi in our record strongly suggests presence of grazing around the lake (Shumilovskikh and van Geel, 2020;van Asperen et al., 2021), confirming previous results (Burney et al., 2003). By reaching the same conclusion, but with some difference in the results, we demonstrate that this is crucial to not only focus on one taxon, but to focus on the whole assemblages of fungal spores.

Local and Regional Comparison With Archeological Data
Comparison with archeological data reported by Douglass et al. (2019) reveals that this early human impact documented by our results on Nosy Be fits well with the major increase of radiocarbon dates in Madagascar and the maximum pace of showing macro-charcoal concentration, percentages of Poaceae, saprotrophs/coprophilous fungal spores and diatom from Aulacoseria granulata agg. and the δ 13 C signal. (B) Evidence of initial human occupation on Nosy Be (Burney et al., 2003) and Mahilaka (Crowther et al., 2016). Hiatus in Lake Mitsijo interpreted as a dry period (Matsumoto and Burney, 1994) and δ 13 C and δ 18 O speleothem profiles from Anjohibe cave showing C 3 /C 4 vegetation and dynamic and precipitation changes (Scroxton et al., 2017). (C) Black curve and brown curve represent the number of 14 C dated archeological samples obtained from Douglass et al. (2019) for the whole (Continued) FIGURE 8 | island and north of 16.5 • S. For each 14 C date, 1000 calibrated dates were drawn according to their probability distribution following calibration with the SHcal20 calibration curve (Hogg et al., 2020). Kernel density estimation with a bandwidth of 60 years was applied to each 1000 calibrated dataset and the median curve and the 5% and 95% percentiles computed. Purple curve shows the estimated changes in the effective population size after Pierron et al. (2017). The orange dashed line shows the main environmental change recorded on Nosy Be. The horizontal brown bar with varying saturation intensity shows timelines of pace megafaunal decline in Madagascar . The red shaded areas correspond to the dry periods defined by the δ 18 O measured on the speleothem and varying saturation intensity to red marks high values within these periods.
megafauna extinction (Figure 8C). Demography history based on genomic data also revealed the same pattern with a major population expansion starting between 1250 and 750 cal BP ( Figure 8C; Pierron et al., 2017). On Nosy Be, archeological surveys have shown first villages and hamlets on the western coast and overlooking Lake Amparihibe and Lake Anjavibe (Dewar and Wright, 1993;Wright and Radimilahy, 2005). Settlements from the coast had evidence of fishing activities while the interior settlements must have concerned solely with swidden cultivation of rice and the herding of cattle, goats and sheep (Wright and Radimilahy, 2005). Consequently, the establishment of such village or hamlet at 1300 cal BP, might be the best fit to explain our findings. Unfortunately, from these archeological surveys, no archeological remains (fish, turtle bones or mammal bones, carbonized coconut shells and local ceramics) have been dated on Nosy Be. However, ceramics are indistinguishable from those found during in Mahilaka, which corresponds to the oldest town described on Madagascar (Dewar and Wright, 1993;Radimilahy, 1997; Figure 1A). Mahilaka was an Islamic port located at 40 km south of Nosy Be in Ampasindava Bay and it formed, together with the Comoro islands, an important hub in the Indian Ocean trade network with several thousand inhabitants (Radimilahy, 1997;Crowther et al., 2016). Study of Crowther et al. (2016) dating rice remains in the area revealed that Comoro islands and Mahilaka were occupied at least from ca. 1250 cal BP (8th century) and from ca. 1000 cal BP (10th century), respectively. In this context, significant human impact at 1300 cal BP on Nosy Be island off the mouth off Ampasindava Bay certainly marks initial settlement of these populations in the region. After ca. 600 cal BP, human occupation declines and Mahilaka and surrounding villages in Ampasindava Bay including Nosy Be were abandoned (Dewar and Wright, 1993;Radimilahy, 1997). The reasons for this abandonment remain unclear (Wright and Radimilahy, 2005). Although recovery of forest and decrease of fires evidenced by our study may reflect this population decline, the continuous occurrences of macrocharcoal particles show that Nosy Be has certainly not been entirely deserted.

Climate Variability
To better understand the impact of these populations on natural ecosystems it is crucial to consider the potential contribution of natural factors such as climatic changes and provide robust comparisons with paleoclimate and paleoenvironmental records available in the region. Specifically, several speleothems from northwestern Madagascar help to address past variability of precipitation in this region (Figures 1A, 8B). These records show that the mid-to late Holocene was characterized by a general drying trend punctuated by decadal-to-multidecadal megadroughts or prominent droughts (Wang et al., 2019;Li et al., 2020). A common feature revealed by these speleothems during the past 3000 years is an important megadrought centered around 900 cal BP (Scroxton et al., 2017;Voarintsoa et al., 2017;Wang et al., 2019;Li et al., 2020). Human impact and associated main environmental changes observed on Nosy Be at ca. 1300 cal BP occur earlier than this megadrought (Figure 8).
Consequently, this drought cannot be considered as the main trigger of the ecosystem shift during the past millennium. A similar pattern is also recorded by δ 13 C from the same speleothems showing the replacement of C 3 woodland habitat with C 4 grasslands preceding the megadrought. To explain this transition, Godfrey et al. (2019) proposed the "subsistence shift hypothesis" which suggests a shift from hunting/foraging to herding/farming between 1250 and 1050 cal BP. During this interval, settlement of new immigrant groups, concomitant expansion of human population on Madagascar combined with spread of pastoralism altered natural ecosystems and accelerated the pace of megafaunal decline on Madagascar Li et al., 2020;Hixon et al., 2021b). Our data support this assumption on Nosy Be, that development of herding/farming activities from 1300 to 1000 cal BP fragmented the rain forest. At ca. 900 cal BP, our results indicate the largest ERD recorded during the past 3000 years and the maximum value of δ 13 C of the bulk organic fraction (Figure 2 and Supplementary Figure 1). This 8.7 cm thick layer might be related to a lake level decrease during the megadrought and it could explain the observed reworked material and the outlier radiocarbon date (Table 1). However, one must remain careful with the interpretation of this layer because humans were already settled in the vicinity of the lake and distinguishing the effect of human impact and climate change on ecosystems is challenging in this context. Additional study sites or precipitation-related proxy are therefore necessary to confirm this hypothesis. After ca. 500 cal BP, more humid conditions are evidenced by speleothems (Scroxton et al., 2017) and may have also contributed to the observed forest recovery on Nosy Be ( Figure 5). However, except for the well constrained megadrought at 900 cal BP, speleothem records are partly showing contradicting precipitation trends, despite most of them originate from the same cave (Supplementary Figure 3 in Supplementary Material from Li et al., 2020), which makes them rather debatable.
In northwestern Madagascar, other study provided paleoecological data ( Figure 8B). Lake Mitsinjo sediments ( Figure 1A) recorded significant grassland development and an increase in charcoal particles occurred from ca. 1000-900 cal BP, that was preceded by a hiatus in sediment deposition (Matsumoto and Burney, 1994;Burney, 1999). In Benavony (Figure 1A), a marshy meadow, the first significant peak of macro-charcoal occurred even later at ca. 800 cal BP (Burney, 1999). Age uncertainties may be responsible for these discrepancies but small-scale spatial differences in the onset and intensity of human activities may also have played a role. Furthermore, influence of climatic changes cannot entirely be excluded, since vegetation changes in northwestern Madagascar occurred in the context of a multimillennial drying trend and several megadroughts occurred during the Holocene (Li et al., 2020). Resilience of every vegetation type in this region (dry deciduous forest, semi-deciduous forest, rain forest or montane vegetation) may have been altered differently to these climatic changes and to the increase of human activities. In this context, without additional and high-resolution paleoecological records it will further remain a matter of discussion to which extent megadroughts affected human impact on natural ecosystem alteration during the past millennium.

CONCLUSION
In this study, we investigated vegetation changes, fire regime, and various sedimentological proxies on the lake system of Lake Amparihibe on Nosy Be, an island located in northwestern Madagascar, which serves as a natural geoarchive. The lacustrine sediments provide insight into environmental dynamics of the past three millennia using multiple paleoenvironmental proxies. Paleoecological results indicate an ecosystem shift at 1300 cal BP characterized by a fragmentation of rain forest related to local fire intensification with development of herding and farming on Nosy Be island. This major ecological change also facilitated distinctively enhanced erosion of the local soils and thus strong sediment and nutrient input into the lacustrine system, resulting in significant aquatic community alterations. Our data confirm a first and strong human impact via agriculture and fire clearance starting at 1300 cal BP on Nosy Be. This supports the "subsistence shift hypothesis, " at a local scale, which proposed that population expansion and intensification of human activities from 1250-1050 cal BP significantly altered the natural ecosystems . However, robust comparison to investigate if this significant and early human impact, after 1300 cal BP, altered ecosystems at a larger scale and not only on the coastal sites still requires high resolution multi-proxy records from continuous sedimentary geoarchives as well as records located further north and in the hinterland. Lack of such data also prevents precise comparison in Madagascar between different regions to discern different or common features of climatic changes, ecosystem responses and anthropogenic influences at the island's scale during the late Holocene. Provide paleoecological data beyond the late Holocene is also strongly recommended to document the influence of natural climate change on ecosystems dynamics and to understand the current state of Malagasy biodiversity (Teixeira et al., 2021). According to the shallow seismic survey and the age model Lake Amparihibe could fill this gap and contribute environmental data for at least the past 8000 years. Previous studies already recommended environmental and human history detailed investigations focusing on different regions of Madagascar (Dewar, 1997). This recommendation is still very relevant today and high-resolution paleoenvironmental data such as the ones provided here, or from new sites, are strongly needed to better understand current biodiversity and to help appropriate management and conservation of Madagascar's unique biomes.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
GD, VM, and TK designed the study. GD, J-JR, and KK organized and conducted the fieldwork. AR, ML, TK, MB, TH, LS, KK, and VM were responsible for laboratory analyses and data interpretation. AR and VM jointly wrote the first draft of the manuscript. ML, TK, RZ, CF, LS, LB, GD, and KK added sections to the manuscript. All authors discussed the manuscript and approved the submitted version.

ACKNOWLEDGMENTS
We thank the Regional Directorate of Water, Sanitation and Hygiene in Antsiranana and local maire for supporting the field work as well as the locals to give access to the holy site. We also thank David A. Burney and Henry T. Wright for personal discussion about their previous researches on northwestern Madagascar. We are thankful to Denis Pierron and his co-authors to share their data for the effective population size. Further, we like to thank emer. Roland Mäusbacher (Friedrich Schiller University Jena) and Antje Schwalb (Technische Universität Braunschweig) for co-financing field investigations of this pilot study. We are also thankful to Carlos Wetzel for taking SEM pictures of the Discostella species. Finally, we acknowledge William D. Gosling and Henry Lamb whose comments on the previous version improved the quality of this paper substantially.