Microscopic charcoals in ocean sediments off Africa track past fire intensity from the continent

Fires in Africa account for more than half of global ﬁ re-carbon emissions but the long-term evolution of ﬁ re activity and its link to global climate change remains elusive. Paleo ﬁ re records provide descriptive information about ﬁ re changes through time, going beyond the range of satellite observations, although ﬁ re regime characteristics are challenging to reconstruct. To address this conceptual gap, we report here the abundance and morphometric data for a large set of microscopic charcoal samples ( n = 128) recovered from surface ocean sediments offshore Africa. We show that in subtropical Southern Africa, large and intense ﬁ res prevailing in open savanna-grassland ecosystems produce a high abundance of small and elongated microcharcoal particles. In contrast, in the forest ecosystems of equatorial and tropical regions of western and central Africa, low-intensity ﬁ res dominate, producing low amounts of squared microcharcoal particles. Microcharcoal concentration and morphotype in marine sediment records off Africa are thus indicative of ﬁ re regime characteristics. Applied to down-core marine charcoal records, these ﬁ ndings reveal that at orbital time-scale intense and large, open grassland-savanna ﬁ res occurred during wet periods in the sub-tropical areas. A strong contribution of ﬁ re carbon emissions during periods of precession and summer insolation maxima in the geological record is thus expected.


F
ire has shaped ecological systems for millions of years and represents one of the most important causes of environmental disturbances on Earth, playing key roles in global biogeochemical cycles, ecosystem patterns and processes [1][2][3][4] .Satellite records demonstrate that more vegetation is burning in Africa than anywhere else on Earth.Fires there account for more than 50% of fire-related carbon emissions worldwide and more than 70% of the global burnt area.Africa hence holds the most fire-prone ecosystems, where fire plays a major role in the global interannual variability of greenhouse gases 4 .
Fire change projections are largely uncertain, especially in Africa 5 where statistical models predict either a decrease 6 or an increase 7 in the same area for the forthcoming century.These empirical relationships are developed using data covering the instrumental period and, therefore, neglect potential changes in interactions and feedbacks between climate, vegetation and fire beyond the range of modern observations.Because the projected climate change far exceeds the time-scale of instrumental or air/ borne observations, process-based models of the coupled vegetation-fire system are used.However, they still need to be evaluated and benchmarked using paleofire data 8 .We rely, therefore, on natural archives to reconstruct long-term fire regime changes under different climate conditions for future scenarios [9][10][11][12][13] .Sedimentary charcoal records provide an estimate of the changes in fire activity or biomass burning and have the potential to offer quantitative measurements directly comparable to simulated values of the fire regime parameters such as fire number, size, frequency, and intensity 11,[14][15][16] .However, despite the efforts of the paleofire international community, we are still missing comprehensive fire-proxy calibrations 8,11,[14][15][16] .Numerous charcoal calibration studies on terrestrial sediment records suggested a positive relationship between charcoal accumulation and one of the physical parameters of the fire regime, i.e., burnt area [14][15][16][17][18][19][20] , intensity 14,15 , fire number 14,21 and severity 22 .Charcoal preserved in marine records is used to reconstruct regional fire activity.It also provides indirect information about hydroclimate and vegetation changes through time 23 as fire is controlled by climate and vegetation.Previous works on the long-term fire history in subtropical African regions highlighted a monsoon control on fire activity at orbital timescales through fuel development.For example, an increase in fire activity was observed in northwestern Africa during the humid early Holocene when high savanna biomass was reconstructed 24 , and in southern Africa during the Late Quaternary when periods of precession and austral summer insolation maxima enhanced grass- 23 and fynbos-fueled fires 25 .In the marine environment, a few published calibration studies have linked charcoal accumulation and morphotypes with burnt area and fire intensity in California 26 and the Iberian Peninsula 27 .In addition, previous work has shown that the distribution pattern of charcoal concentrations correlated with charcoal influx in North Pacific marine sediments 28 , suggesting that charcoal distribution patterns in sediments are not controlled by any dilution effect but instead faithfully reflect the sedimentary flux of charcoal particles.Nevertheless, there are still important limitations in relating fire proxies from natural records with fire regime characteristics.This is due to the different temporal and spatial resolution between natural archives and instrumental records, the diverse charcoal amount produced during combustion, the difficulty of tracing charcoal source area and various taphonomic processes 11 .The first inferences of fire regime characteristics from marine paleofire records were derived using indirect theoretical principles developed on charcoal from lake sediments and experimental studies.A peak in microcharcoal concentration or influx was interpreted as an increase in fire frequency and intensity while the particle shape was used to infer intensity and vegetation type.Recently, it has been demonstrated that charcoal abundance and morphotype in marine surface sediments off Iberia can trace specific fire regimes on land, in particular fires of high intensity in a mixed open Mediterranean woodland vegetation 29 .However, there is no empirical data yet to support paleofire regime inferences made from marine microcharcoal records off Africa.
Here, we aim to explore the relationships between microscopic charcoal (microcharcoal) characteristics preserved in recent marine sediments along the Atlantic margin off (sub)-tropical Africa (Fig. 1a) and the changes in fire regimes in adjacent regions.To this purpose, we determined the abundance and morphotypes of microcharcoal in surface sediments and compared the obtained data with satellite-derived fire information (fire number, size, intensity and type of burnt vegetation) over the recent period.Based on these results, we propose a new framework linking microcharcoal accumulation and morphometric values in marine sediments to burned vegetation types and fire regimes for Africa.We apply this framework to marine charcoal records located off subtropical African regions to discuss changes in the paleofire regimes over a multimillennial timescale.

Results
Distribution pattern of sedimentary charcoal in seafloor sediment.Microscopic charcoal particles, between 10 and 92 µm in length, were extracted from surface (core-top) marine sediment samples (n = 128; Material and Methods) collected at the African Atlantic margin, from the Gulf of Guinea to the Cape Basin, and on the Indian side, from the Agulhas Plateau to the Limpopo region (Fig. 1a).We analyzed the concentration and morphometric characteristics of microcharcoal, including area, length, width and elongation (length/width ratio), to examine charcoal distribution patterns in ocean sediments and their relationships with fire regimes on land.
Following vegetation burning, microcharcoal particles are released and transported with the fine-grained sediment fraction to the ocean, via fluvial systems or aerial dispersion 29 (SI Section 1).Given the size and weight of microcharcoal particles, ranging between >10 µm and <100 µm in length 30 , they behave similarly to pollen and fine terrigenous particles during transport and sedimentation 31 .In marine sediments proximal to river plumes, fine terrigenous particles are primarily of fluvial origin 32,33 .Only off arid regions, where no hydrographic system is present, the aeolian transport can represent the dominant carrier of fine particles 34 .In the ocean, fine particles become an integral part of the marine snow, reaching the sea floor within a few days or weeks 35,36 .In addition, the latitudinal pollen distribution in marine surface sediments represents faithfully the vegetation distribution on the adjacent continent 32,35,[37][38][39] .Although microcharcoal and pollen share similar fluvial behaviour during fluvial transport, they however behave slightly differently during aerial dispersal due to differing physicalmechanical characteristics 40 .In this study, we therefore evaluated the potential contribution of windblown microcharcoal particles using a dispersion and deposition model (see Material and Methods and SI Section 3) showing that wind and fluvial source areas overlap so that most of charcoal produced in one hydrographic basin is deposited in the same basin and can then be transported to the ocean by the fluvial system.
Consequently, we defined eight regional source areas for sedimentary microcharcoal, roughly based on the proximity to the river system entering the ocean, i.e., Gulf of Guinea, Congo, Angola, which represent the central western Africa region, Namibia, Orange, western South Africa (WSA), Eastern South Africa (ESA) and Limpopo which represent Southern Africa (Fig. 1a).Charcoal preserved in samples off southern Africa (the southern area of study) reflects therefore fire activity of subtropical Africa and charcoal preserved in samples off central western Africa reflects the fire activity of equatorial/tropical Africa.Samples located off the Gulf of Guinea, Orange and WSA share similar, relatively low mean microcharcoal concentration values, below 100,000 nb.g −1 (Fig. 2).Intermediate mean microcharcoal values, between 100,000 and 200,000 nb.g −1 are typical for Angola, Congo but also Namibia and WSA.The highest microcharcoal concentration value is attained by the sample located off Limpopo (Fig. 2).Samples located off Namibia, Orange, WSA and Limpopo show higher mean elongation ratio (>1.8) whereas samples located off Gulf of Guinea, Angola, Congo and ESA coast display values <1.8 (two-tailed t-test, p < 0.05) (Fig. 2).We applied principal component analysis (PCA) on mean ocean charcoal parameters, including concentration, elongation, width, area and length, to identify the association between parameters and potential clustering of charcoal particle types (Fig. 3a).The first two components of the PCA on the mean charcoal (concentration, elongation, width, area and length) per hydrographic basins explain almost 80% of the overall variance of the mean data (Fig. 3a).PC 1 is correlated with area, length and width of charcoal particles, and PC 2 is correlated with charcoal elongation and concentration strength, cos 2 > 0.5, Fig S .5.The PCA results show high heterogeneity in the distribution of samples.However, strong relationships between charcoal abundance and morphotype were established among the different hydrographic basins and between western central and southern Africa (Fig. 3a, Kruskal-Wallis tests, p < 0.05).The Kruskal-Wallis test applied on PC2 (correlated with concentration and elongation) shows evidence for significant differences between basins (Kruskal-Wallis chi-squared = 21.199,df = 7, p-value = 0.003486) as well as a significant difference between the two main regions (Kruskal-Wallis chi-squared = 8.1685, df = 1, pvalue = 0.004262).At sub-continental scale, microcharcoal patterns in marine sediments off Africa are described by two main variables, elongation and concentration.The shape and size of microcharcoal particles group together samples from Namibia, Orange and SAW (southern Africa -southern region), whereas the shape of particles groups samples from Congo, Angola and Gulf of Guinea basins (western central Africa).Samples from western central Africa correspond to squared particles of variable size present in low concentration in the sediment.In contrast, most samples from Southern Africa are characterized by small, elongated particles in high concentration (Fig. 3a).
Charcoal concentrations and morphotypes informing fire regime changes.To identify the potential link between charcoal accumulated in the ocean and fire on land, we compared the detected patterns for marine sedimentary charcoal with various parameters of fire regimes on adjacent continental regions (Fig. 3b).PCA was applied on land fire parameters including number of fires (Nb fires), mean fire radiative power (FRP), area, rate of spread (RoS) and the burnt vegetation types (graminoid -G, mixed -M and trees -T) to characterize the fire regime specific for each hydrographic basin (Fig. 3b).We used the mean FRP in W.m −2 as a measure of fire intensity.The FRP index measures the energy emitted through radiative processes released during the combustion and can be associated with fire intensity throughout the fire-burning process 41 .Principal component axis 1 explains 32%, and axis 2 explains 23% of the variance in all variables.PC 1 opposes FRP and mixed burnt vegetation to burnt trees, while PC2 opposes RoS to burnt graminoids (Supplementary Fig. 5).Samples from Congo, Angola and the Gulf of Guinea tend to cluster on the positive values of PC 1 and are associated with burnt trees.Samples from Orange, SAW and Namibia are found in the negative part of PC 1 and are associated with FRP, burnt mixed, RoS, area and burnt graminoids.The Kruskal-Wallis test shows that fire regimes in each hydrographic basins  2), ocean currents and location of ocean sediment samples analyzed and the studied hydrographic basins with abbreviated names (GG Gulf of Guinea; C Congo; A Angola; N Namibia; O Orange; L Limpopo; SAW South Africa West; SAS South Africa South).GC Guinea Current; SEC South Equatorial Current; EUC Equatorial Under Current; SECC South Equatorial Counter Current; AC Agulhas Current; ARC Agulhas Return Current; EMC East Madagascar Current.b The location of ocean sediment samples (dots) color-coded based on the preferential microcharcoal source area (areas on land) with fire density (calculated from the FRY 2 version as the mean yearly fire number) and major wind direction.The maps were created in ArcMap version 10.8.1.The background map was downloaded from European Space Agency (ESA) Climate Change Initiative (CCI) Land Cover (LC) database 77 and ArcGIS Online.The fire density was generated using the FRY dataset global fire patch database.(Kruskal-Wallis chi-squared = 419.37,df = 7, p-value = 0.000) and between distinct regions, i.e., western central and Southern Africa, are statistically different (Kruskal-Wallis chi-squared = 13.626,df = 1, p-value = 0.0002231).Our assessment of African fire regimes (PCA) indicates that small (local) fires of low intensity occurring in tree-dominated areas mostly prevail in western central Africa, whereas large (more regional) fires of high intensity and high rate of spread involving mixed, graminoid-shrub vegetation are mostly encountered in southern Africa (Fig. 3).Taken together, our results indicate that low amounts of charcoal particles with squared shape are observed in marine sediments proximal to regions dominated by small fires of low intensity in tree-dominated ecosystems, whereas larger amounts of small charcoal particles with elongated shape are observed in proximity to landscapes experiencing large and high-intensity fires in graminoid-mixed ecosystems 32 .

Discussion
Some relationships between charcoal accumulation or morphotypes and individual fire parameters have been previously documented at local scales in lake records from various climatic and vegetation settings, i.e, boreal 22 , temperate 14 to grassecosystems 18,42 , as well as in marine sediments 26,27 .More specifically, a similar approach 27 albeit at a smaller spatial scale and restricted to temperate and Mediterranean fire ecosystems of the Iberian Peninsula, reported higher charcoal concentration and elongation values in marine samples associated with fires spreading in open Mediterranean woodlands in Iberia 27 .Here, we demonstrate at a larger sub-continental scale and for environmental settings ranging from tropical to subtropical climate, tropical to desert biomes, that concentration and elongation of charcoal in marine sediments off Africa can be used to identify different fire regimes.At this sub-continental scale, it is difficult to argue that transport and/or sedimentation processes mostly control the observed distribution pattern of sedimentary charcoal in seafloor sediments.
Instead, previous studies have shown that the amount of charcoal is generally controlled by production source characteristics such as fuel type, biomass, fire size and intensity 43,44 .A possible explanation for the observed distribution pattern of charcoal morphotypes (concentration, shape and size) could result from different transport mechanisms between regions influenced by grain-size sorting during eolian transport and others dominated by river discharge, or by surface ocean circulations.However, all regions combine the influence of large river basins (Congo and Orange rivers) and intense ocean currents (BCC and AC, Fig. 1a), and prevailing southeasterly trade winds.Another possibility lies in a higher sediment export concentrating charcoal in preferential areas.However, both the Benguela upwelling system and the Congo fan are high sedimentation accumulation systems.Consequently, it seems unlikely that sediment transport processes control the observed latitudinal distribution pattern of charcoal in seafloor sediment.
The high charcoal concentrations observed in the marine samples off Southern Africa are associated with large fires of high-intensity spreading in graminoid-mixed vegetation i.e., open grassland-savanna ecosystems.Here, the graminoid-mixed vegetation is composed of a mixture of shrubland, grassy semi-desert and arid-fertile savanna, fynbos (winter rainfall shrubland) and vlei/dambo (seasonally-dry wetlands), leading to exceptionally high-intensity fires 45 .Rare, intense and large fires dominate the African savannas at its driest range 46 , with fuel availability limited for several years but sufficient and continuous during the rare wettest years to support large fires.Fires of high intensity consume more vegetation and burned large areas in southern Africa 15 .Such a relationship between high charcoal accumulation and high-intensity fires has been identified locally in several water bodies from the boreal region 22 and from the Kruger National Park 15 and in marine sediments off Iberia 27 .Our findings are consistent with the previous work suggesting that charcoal accumulation in sedimentary records reflects a shift in fire regimes rather than a change in a single fire parameter such as fire size, intensity, or frequency 30 .
Higher biomass availability in forested ecosystems was suggested to be the main parameter explaining higher charcoal accumulation in lacustrine sediments 47 or elemental carbon accumulation in marine proximal sediments 48 .However, the fuel load is low in subtropical African grasslands and savannas 49 , thereby arguing against control of overall biomass availability on the concentration of sedimentary charcoal in nearby sediments.Conversely, it has been shown that savanna fires consume leaves at 95% but wood only at 5 to 50% 10 .Therefore, the fuel material combusted by fires might exert a strong control on charcoal amount.In this vein, we propose that the intensity of wildfires in Africa plays a major role in controlling the distribution patterns of associated sedimentary charcoal particles.This finding is fully consistent with some experimental studies showing a positive correlation between the amount of charcoal particles and fire intensity 15,22,31,50,51 .In addition, based on the combustion continuum, the size of charcoal is reduced from cm to micron scale with increasing temperature 52 .We thus infer that the high proportion of small-sized (microscopic) charcoal particles observed in Southern Africa results from high intensity fires due to the high burning temperature of fine fuel (grasses and leaves) which produces high fire emission (gases) 4,49 .Our results also demonstrate that the burning of herbaceous fuels produces much smaller particles than the woody fuels.However, our findings contrast with the experimental work suggesting that low-temperature surface fires produce more elongated particles than hightemperature fires 53 .This opposite trend is not surprising considering the different biogeographical regions (Africa vs Siberia) corresponding to vegetation types/plant associations (grasslandsavanna vs boreal forests) in which fire behaviour and charcoal production and type are expected to differ significantly.
Experimental studies showed that charcoal particles produced by graminoid and graminoid-shrub burning (median elongation ~3) are more elongated than particles derived from the burning of tree and (sub)shrub (median elongation ~2) 54 .Results on the mean morphology of charcoal in our marine sediment samples suggest that southern African fires spread in mixed-graminoid vegetation (high mean elongation values of charcoal samples), whereas western central Africa is marked by forest fires (low mean elongation of charcoal samples) (Fig. 3a and b).This confirms our ocean-land comparison in terms of burnt vegetation type.Our data present much smaller mean elongation values (~<2.5)than those synthetized by experiments 55 .We argue that this difference arises from the fact that 54 compiled elongation data from the burning of a single type of vegetation (grassland and woodland).However, this is rather unrealistic for in-situ environments where fires spread across different vegetation and species types 55 .In samples from central and southern Africa charcoal assemblage (at sample level) is likely composed of particles coming from different burnt vegetation types and species.Another possible explanation for the observed smaller mean elongation values in our sedimentary charcoal samples could be related to the physical breakdown of charcoal into smaller grainsize with increasing distance of fluvial transport.
Thus, we propose a mean elongation classification at charcoal assemblage level to infer burnt vegetation type from charcoal in marine sediments (Supplementary Fig. 8).A value of <1.8 ± 0.3 representing tree burning is obtained by calculating the mean charcoal elongation for samples belonging to central western Africa.A value of >2.1 ± 0.5 obtained from Southern Africa charcoal samples tracks the burning of graminoids-mixed vegetation.Mean charcoal elongation from central western Africa is significantly different from values from southern Africa (t-test, p < 0.05) (Supplementary Fig. 8).The threshold value of ̴ 2 obtained here using the length/width ratio is in line with results from 41 tropical grassy ecosystems 41 and 53 from grass prairie ecosystems of North America 53 .These studies showed that a width/length ratio of 0.5 or less (an inverse length/width ratio >2) characterize the burning of grassland-dominated vegetation, while values >0.5 represent the combustion of ligneous vegetation (including both shrubs and trees).
At sub-continental scale, our results clearly show that microcharcoal assemblages in marine sediments depend on the type of fire regimes on adjacent continental margins.As fires are expected to change predictably along climate-environmental gradients 56 such a relation provides a valuable source of information about fire-biome shifts over different time scales.
Previous work on paleofire in subtropical Africa suggested that at orbital scale the monsoon controls the fire activity through fuel development.In a marine sediment record off southwest Africa (MD96-2098 23 ), higher microcharcoal concentration and more elongated particles were reported for glacial periods and during austral summer insolation maxima between 170 ka to 35 ka.These changes were interpreted as reflecting increasing fire activity during cool and wet periods associated with increasing grassland biomass.Under the same configurations (glacial times and austral summer insolation maxima), another study 25 showed high microcharcoal concentrations associated with fynbos expansion in a 300 ky-long (Site U1479) marine sedimentary record off south-western South Africa, which they interpreted as reflecting fynbos fires.In northern Africa higher microcharcoal concentration was reported during the early humid Holocene in core GeoB7920-2 associated with grass-savanna expansion 24 and interpreted as higher fire occurrence.
The modern land-ocean relationship proposed in this study helps interpreting the three marine charcoal records in terms of past fire regime changes (Fig. 4).In the marine sedimentary core off Namibia 23 , the high microcharcoal concentration associated with high elongation values at times of monsoon enhancement (austral summer insolation maxima) were interpreted as an increase in grass fires of low intensity.Our mean elongation classification applied on the long sedimentary core MD96-2098 suggests that in addition to grass fires 23 , mixed-fuel fires were also occurring.Based on our modern fire land data, contrary to previous interpretations, we suggest that the fires in southern Africa were of high intensity.We propose therefore an alternative interpretation and suggest that those increases in subtropical regions of Africa indicate a fire regime characterized by large fires of high intensity spreading in graminoid-mixed vegetation.Grass-fueled fires of low intensity were inferred on the basis of data for fire characteristics in sites located outside Africa 23 , under different environmental settings, which could explain the observed bias in interpretation.Looking at southwestern South Africa we propose that over the past 300 ky increased charcoal concentration at times of fynbos expansion revealed intense and large fires but not necessarily frequent fires 25 .In northwestern Africa, our land-ocean relationship suggests that fires were large during the humid early Holocene and not numerous 24 .
In summary, microcharcoal concentration and elongation from ocean sediment samples off the African coast portray the fire regime changes on land.On the one hand, high amounts of small and elongated charcoal track intense and large fires spreading into open grassland-savanna vegetation.On the other hand, low amounts of squared charcoal indicate small fires of low intensity, spreading in tree-dominated ecosystems.Among key processes, fire intensity appears to be the most important driver of charcoal morphotypes and abundance in Africa.Our work has important implications for charcoal calibration efforts developed by the international paleofire community.Our study shows that charcoal preserved in marine sediments off Africa is a valuable proxy to trace past open grassland-savanna fires and inform us about their variability, controls, and response to climatic changes.At the orbital timescale, our results suggest that intense and large open grassland-savanna fires occur during wet periods in subtropical regions of Africa.These intense fires might have been a source of carbon to the atmosphere during periods of precession and summer insolation maxima in the geological record.However, more geological data as well as model simulations, are necessary to quantify the contribution of fires to the atmospheric carbon cycle over orbital timescale.

Material and methods
Different types of corers were used to retrieve the sediments including gravity, piston and box corer during different marine cruises including CONGOLOBE 2011 57 , ZAIANGO set of cruises run between 1998 and 2001 [58][59][60] , RV Tyro 1989 61 and MD105 Images II NAUSICAA 62 .Samples were collected from the top 0 to 1 cm of these sediment cores stored in different repositories (EPOC and IFREMER, France) or obtained from PIs personal repositories.The samples are spread across different depositional settings, from estuary (1 sample), the continental shelf (1 sample), the continental slope (51 samples) to the abyssal plain (75 samples).twenty-one less than 150 years, twenty-four less than 200 years, fifteen less than 500 years and twelve of several hundred years up to late Holocene, and six samples of Holocene age.We were unable to date some core-top samples using the 210 Pb method since some of the marine sediment cores were collected during oceanographic cruises from the 1990s.The ages of 46 core-top samples are based on direct chronologies of the sediment cores.The ages of 13 samples were inferred from closeby cores with direct chronologies.The year of the collection was assigned to one sample from the estuary and to 4 samples from multi-cores that preserved the sediment-water interface (presence of the oxic layer).For the other 4 samples that had no closeby cores with direct chronologies, we assigned a Holocene age based on foraminifera assemblages or benthic oxygen isotopic measurements.For the remaining 60 samples located off the western and central Africa coast, we inferred their ages from published sedimentation rates in multiple cores surrounding our samples 63,64 (please see Data availability).The study of 64 presents means terrigenous fluxes (g/cm 3 /ka) during the Late Quaternary.We calculated the sedimentation rates by dividing the terrigenous flux by the dry density (g/cm 3 ), using a dry density of 0.776 g/cm 3 (estimated using a humid density of 1.5 g/cm 3 , and a 2.65 grain density and a 0.707 porosity).For each sample, we estimated a maximum average sedimentation rate from which we inferred the largest estimated ages (calendar years).We assumed that the range of estimated ages includes age uncertainties.
Studies 65 show that core-top of marine sediments represents the late Holocene and covers few hundred to few thousand years which might be problematic in terms of comparison with more recent observations.However, prior studies have successfully used spatial calibration of marine surface sediment-based proxies to reconstruct different physical and biological ocean and sea-ice properties (for example [66][67][68][69][70] ) or fire regime characteristics 27 suggesting that core-top samples are suitable for our calibration approach.In particular, a study using core-top sediments off Iberia 27 showed that the temporal variability of microcharcoal concentrations in interface cores off Iberia is within the range of modern spatial variability.This supports our comparison of the different datasets as mean environmental conditions on land and in the ocean over the past 100 years, corresponding to these past century biomass burning records, did not radically depart from conditions recorded by the instrumental period.For the development of future calibration studies, it is recommended to focus on sediment samples collected with a multitube corer device that preserves the mud-water interface.
Sample preparation.Charcoal particles are produced at temperatures ranging between 200°to 600 °C52 .They are inorganic carbon, thus inert and resistant to biological and chemical decomposition 29 .The microcharcoal preparation techniques were adapted from previous work 71 , classical palynology methods 72 and petrographic mounted slides 73 .Microcharcoal samples were processed following the sediment preparation protocol previously published 74 slightly modified depending on the composition (organic-, silica-, carbonate-rich) of marine surface sediment samples.The microcharcoal extraction technique was applied to approximately 0.2 g of dried bulk sediment and consisted of the following steps: (a) a chemical treatment of hydrochloric acid (37% HCl) to remove calcium carbonates (CaCO 3 ) until the reaction, i.e., degassing bubbles, stopped; (b) followed by nitric acid (68% HNO 3 ) and hydrogen peroxide (33% H 2 O 2 ) attack on warm water-bath to remove or bleach the organic matter; (c) centrifugation to remove the supernatant before; (d) a chemical attack of hydrofluoric acid (HF) to remove siliceous material followed by centrifugation to remove HF; (e) followed by an HCl treatment to remove colloidal SiO 2 and silicofluorides formed during the HF digestion followed by a centrifugation to remove HCl.For samples with more than 1% organic matter, we added additional steps with alternating nitric acid and hydrogen peroxide chemical attack.This chemical treatment is used to remove carbonates, pyrites, humic material, labile or less refractory organic matter (OM), to bleach non-oxidized OM and to remove silicates.A dilution of 0.1 is applied to the residue.The suspension was then filtered onto a membrane of 0.45 µm porosity.A portion of this membrane was mounted onto a Polymethyl Methacrylate (PMMA) slide with ethyl acetate before gentle polishing with alumin powder.The microscopic charcoal particles were identified under oil immersion using petrographic criteria in reflected light.The microcharcoal quantification was performed in transmitted light using image analysis (Leica DM6000 microscope and LAS software).Only charcoal particles >10 µm were analyzed.Previous work 75 reported some fragmentation of particles when using HNO 3 or H 2 O 2 on modern charcoal and on peat sediment containing charcoal.The microcharcoal concentration (CCnb representing the number of microcharcoal particles per gram of dry weight sediment) and microcharcoal surface area (CCsurf, representing the sum of all surfaces of microcharcoal particles in one sample per gram) are significantly and positively correlated (Supplementary Fig. 7, r = 0.86, p < 0.01), suggesting that potential fragmentation due the chemical procedure 75 is minimal.Both parameters record the same pattern of microcharcoal variability highlighting that CCnb can be used here to inform about these changes.A bias in our charcoal record might come from the use of H 2 O 2 that may bleach or remove partially charred particles resulting from low temperature fires 76 .To minimize the error, we followed their recommendations by using solutions of consistent strength to treat all samples.
Geographic information system (GIS): continent, ocean and fire data.Information about the hydrographic basins was extracted from Natural Earth Data and GRIN database 77 .The characteristics of the hydrographic basins (size, hypsometry, river network) were extracted from the ASTER Global Digital Elevation Model free dataset.
The information regarding vegetation type was extracted from European Space Agency (ESA) Climate Change Initiative (CCI) Land Cover (LC) database 77 .The twenty-eight land cover classes from this database 77 were reclassified into 3 major growth habitat categories: graminoid including non-forested land cover types, mixed including mosaic land cover type and closed including forested types (Fig. 1a, Supplementary Table 2).The other land cover types were included in the not-burning (NB) category (Fig. 1a, Supplementary Table 2).
Fire regimes characteristics, including fire number, fire area, fire intensity (mean fire radiative power), rate of spread (derived from fire size and duration) and the three most important types of burnt vegetation (according to global Land Cover Climate Change Initiative, CCI-LC) were extracted from the global fire patch database FRY (2001-2019) 41,78 , updated with the newly delivered FireCCI51 250 m resolution pixel-level burn dates 79 .We used the mean FRP (Watt/m 2 ) as a measure of fire intensity; FRP measures the energy emitted through radiative processes released during the combustion and can be associated with fire intensity all throughout the fire-burning process 78 .FRP is used for assessing carbon emissions by converting FRP into fuel biomass consumption rate 80 .FRP is also a key component of global fire regime characterization 81 .For ease of interpretation, the most important burnt vegetation classes (L1, L2, L3) provided in FRY v2 were translated into 3 land cover classes according to the vegetation growth habitat they were spreading into (Supplementary Table 2).We used this reclassification to have the same level of information from land fire data and charcoal from the sediment records and examine the relationship between them.The information on fire patches from FRY v2 were aggregated on a 1°× 1°grid for each hydrographic basin by calculating the median of the number of fire patches (>100 ha), burned area, FRP and RoS per vegetation type.
Modelling the dispersion of fire plumes and charcoal particles deposition.The Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT) 82 plume dispersion model and computation of particle concentration was used to simulate the dispersion of fire plumes and deposition of charcoal particles in our attempt to estimate the most likely airborne charcoal source area.Based on the physics of particles, this simulation uses a fixed number of particles which are transported by a wind field, spread by a turbulent component and, in our case, removed by dry and wet mechanisms 82 .For this simulation we used meteorological fields from National Centers for Environmental Prediction's (NCEP) Global Data Assimilation Scheme (GDAS) with 1°× 1°r esolution and 3-hourly data covering the interval 2007 to present.
We tested the Hysplit air dispersal and deposition model in ten fire locations, for three atmospheric levels and three microcharcoal particles sizes that describe the range of our dataset (SI Section 3).Details regarding parametrization and test simulation are provided in SI Section 3. The results were similar in most fire locations showing that the distance traveled by the particles, regardless of their size and atmospheric level, is directly proportional with the injection height and most of the microscopic charcoal stays within a few tens to a few hundreds of kilometers (Supplementary Table 1) from the fire area.From spatial charcoal concentration maps (Supplementary Fig. 3 maps and Supplementary Table 1) wind and fluvial source areas overlap so that most of charcoal produced in one hydrographic basin is deposited in the same basin and can then be transported to the ocean by the fluvial system.Considering the closed land-sea link between pollen assemblages and latitudinal vegetation distribution in Africa 37 , we anticipate microcharcoal deposited in marine sediments come from the closest hydrographic basins and reflect spatial fire regimes distribution.
We therefore used the limits of the hydrographic basins and distance from the river mouth and coast to identify the closest source-areas.Samples located within the Congo system (using the oceanographic limits) were assigned to Congo hydrographic basin.We obtained eight zones as follow, from N to S: Gulf of Guinea (31 samples), Congo (72 samples), Angola (14 samples), Namibia (4 samples), Orange (2 samples), western South Africa (4 samples), Eastern South Africa (4 samples) and Limpopo (1 sample) (Fig. 1a).

Statistical analysis.
A set of statistical analyses were run and visualized in Rstudio Version 1.3.1093 83using packages Tidyverse, FactoMineR and Factoextra.Principal Component Analysis (PCA), centered and scaled, was used to assess the potential clustering of charcoal particle types, including their area, length, width, elongation (length-to-width ratio) by hydrographic basins.We used squared cosine (cos2) to determine the quality of representation and strength of the relationship for each variable on the dimension.The same analysis was applied to land fire parameters, including the number of fires, fire radiative power (FRP), area, rate of spread (RoS) and the burnt vegetation types (graminoid -G, mixed -M, trees -T) to characterize the fire parameters specific for hydrographic basins.Given the high number of satellite observations, before running the analysis, the median of fire regime parameters was calculated per grid cell.Kruskal-Wallis rank sum test (KW test, kruskall test() function, R cran) was applied on PC1 values of ocean charcoal and land fire data to determine whether there is a significant difference among hydrographic basin distribution on the PCA driven by charcoal parameters for the ocean data set, and land fire parameters for the land data set.The KW test was also run to determine the statistical differences among regions, western central versus Southern Africa.Bilateral student t-test were performed on the elongation dataset to determine the mean differences between hydrographic basins, t.test() function, R cran.

Fig. 1
Fig. 1 Study area characteristics and ocean sediment sample location.a Study area with main vegetation types following ESA-CCI (Supplementary Table2), ocean currents and location of ocean sediment samples analyzed and the studied hydrographic basins with abbreviated names (GG Gulf of Guinea; C Congo; A Angola; N Namibia; O Orange; L Limpopo; SAW South Africa West; SAS South Africa South).GC Guinea Current; SEC South Equatorial Current; EUC Equatorial Under Current; SECC South Equatorial Counter Current; AC Agulhas Current; ARC Agulhas Return Current; EMC East Madagascar Current.b The location of ocean sediment samples (dots) color-coded based on the preferential microcharcoal source area (areas on land) with fire density (calculated from the FRY 2 version as the mean yearly fire number) and major wind direction.The maps were created in ArcMap version 10.8.1.The background map was downloaded from European Space Agency (ESA) Climate Change Initiative (CCI) Land Cover (LC) database77 and ArcGIS Online.The fire density was generated using the FRY dataset global fire patch database.

Fig. 2
Fig. 2 Microcharcoal results derived from surface marine sediment sample.Boxplots of mean microcharcoal concentration and mean elongation ratio derived from marine sediment samples in each zone.

Fig. 3
Fig.3Statistical results on microcharcoal assemblages from ocean sediment samples and land fire data.a Results of the principal component analysis on microcharcoal assemblages including length, area, width, elongation and charcoal concentration (Ccnb) with the interpretation of the four spaces defined.b Principal component analysis on land fire parameters (derived from FRY2 version at 1°× 1°grid) including fire number, area (ha), mean fire radiative power (FRP, W/m 2 ), rate of spread (RoS) for major burnt vegetation types (G Graminoids, M Mixed, T Trees).Ellipses with different colors and symbols represent the hydro-basins (please note that Namibia and Limpopo do not form an ellipse on the ocean data due to the restricted number of samples).The legend for the eight hydrographic basins (colors and symbols) are presented on the right side of the chart.Graminoids correspond to grassland-savanna vegetation, and trees correspond to woodland vegetation.

Fig. 4
Fig. 4 Proposed fire regimes scenarios obtained from the analysis of microcharcoal types and concentration produced and deposited in the marine environment off Africa.The forested ecosystem of central western Africa is dominated by small and low-intensity fires producing squared microcharcoal particles in low concentration.Intense and large fires are specific for southern Africa with fires burning graminoid-mixed vegetation and producing small, elongated particles in high concentration.The figure was created in Affinity Designer 1.10.