Ground-based climate data show evidence of warming and intensification of the seasonal rainfall cycle during the 1960–2020 period in Yangambi, central Congo Basin

Meteorological stations are rare in central Africa, which leads to uncertainty in regional climatic trends. This is particularly problematic for the Congo Basin, where station coverage decreased significantly during the last few decades. Here, we present a digitized dataset of daily temperature and precipitation from the Yangambi biosphere reserve, covering the period 1960–2020 (61 years) and located in the heart of the Congo Basin. Our results confirm a long-term increase in temperature and temperature extremes since the 1960s, with strong upward trends since the early 1990s. Our results also indicate a drying trend for the dry season and intensification of the wet season since the early 2000s. Ongoing warming and increasing precipitation seasonality and intensity already have a significant impact on crop yields in Yangambi. This calls for urgent development of climate-smart and dynamic agriculture and agroforestry systems. We conclude that systematic digitization and climate recording in the Congo Basin will be critical to improve much-needed gridded benchmark datasets of climatic variables.


Introduction
Several authors report evidence of a drying trend in central Africa during the last few decades, based on available interpolated gridded climate datasets (Tsalefac et al. 2015;Dezfuli and Dezfuli 2017;Nicholson et al. 2018bNicholson et al. , 2019;;Mabrouk et al. 2022) and satellite products (Diem et al. 2014;Zhou et al. 2014).This drying trend started between 1990 and 2000 and is thought to be one of the most significant worldwide (Malhi and Wright 2004; Extended author information available on the last page of the article 142 Page 2 of 28 Asefi-Najafabady and Saatchi 2013; Zhou et al. 2014;Nicholson et al. 2018b).The longterm drying trend in central Africa was recently supported by in-depth analysis of multiple independent precipitation and satellite-derived vegetation datasets, showing that the dry season length increased by 6.4-10.4days per decade in the period 1988-2013, due to an earlier dry season onset and a delayed dry season end (Jiang et al. 2019).Furthermore, some climate model projections show decreasing precipitation towards the end of the twenty-first century under several RCP scenarios (Fotso-Nguemo et al. 2017).At the same time, wet and dry extreme events are expected to become more frequent and more severe (Fotso-Nguemo et al. 2018;Kendon et al. 2019;Karam et al. 2022).
Nevertheless, large uncertainties on past and future central African climate variability persist because meteorological stations are rare in central Africa (Aguilar et al. 2009;Chaney et al. 2014;Kidd et al. 2017;Nicholson et al. 2018a;Bush et al. 2020;Nicholson 2022).While central Africa was covered by up to 1000 climate stations in the 1960s, 1970s, and 1980s, the network dramatically collapsed afterwards (Nicholson et al. 2018b;Nicholson 2022) (Fig. 1).Especially, the Democratic Republic of the Congo has been poorly covered in recent decades, with only 17 climate stations providing data for the 2000-2016 period.All these stations are located in the South-West of the country, leaving a large gap of more than a million km 2 in the central rainforest area (Fig. 1).This data gap exists due to long-lasting political instabilities and restricted data-sharing policies, which impeded data acquisition, data digitization, and use of the data in a scientific context (Dezfuli and Dezfuli 2017; Nicholson et al. 2018b).Furthermore, satellite products often do not perform well in the tropics (McCollum et al. 2000;Awange et al. 2016;Sun et al. 2018;Nicholson et al. 2019;Igri et al. 2022).This leads to substantial uncertainty in regional long-term climate trends (Washington et al. 2013;Awange et al. 2016).Uncertainties are especially high for rainfall due to high spatial variability in precipitation regimes (Kidd et al. 2017;Bush et al. 2020).Additional to these uncertainties, regional and global climate models show opposing future trends for most of central Africa, suggesting that both past and projected climate change signals may not be robust (Dosio et al. 2019).
Yet, a good understanding of climate variability is of utmost importance to predict potential impact of future climate change in central Africa.Continued long-term drying or intensification of extreme events may likely have large impact on tropical forest composition, functioning, and carbon sequestration (Fauset et al. 2012;Denbow 2013;Zhou et al. 2014;Aguirre-Gutiérrez et al. 2020;Hubau et al. 2020;Bennett et al. 2021).Furthermore, climate change has large impacts on forestry, agriculture, agroforestry and health (Verchot et al. 2007;Diem et al. 2014;Salerno et al. 2019;Tschora and Cherubini 2020).Hence, it is of utmost importance to further investigate central African climate trends by collecting much-needed ground-based data.To do so, the Trans-African Hydro-Meteorological Observatory (TAHMO) is developing a network of modern weather stations across central Africa, which is monitoring climate since 2015 and will continue to do so in the future.Yet, long-term records including historical data are equally important to analyze long-term climate variability and to compare with climate model simulations of past and projected future changes.
Therefore, we visited the climate station of the Yangambi Man-and-the-Biosphere reserve, which is a leading research site in the Democratic Republic of the Congo and is located in the heart of the data-poor area in central Africa (Fig. 1).Although the Yangambi data were not digitized since the 1990s, the station remained operational.We digitized the daily temperature and precipitation observations for the period covering 1960 to 2020 (61 years), we calculated monthly and yearly meteorological indices (describing temperature extremes, rainfall intensity, water demand and availability, drought and seasonality), and we performed long-term and piecewise linear regression to detect long-or short-term trends.We used this analysis to address the following questions: (1) how did temperature and temperature extremes vary in Yangambi (Fig. 2); (2) how did precipitation and rainfall intensity vary in the dry and wet season (Fig. 3); (3) how did water demand and precipitation seasonality vary (Fig. 4)?Finally, we investigated how well observed records of precipitation and temperature from the Yangambi meteorological station compare with those from existing gridded datasets (Fig. 5 and Fig. 6)?

Yangambi station and daily meteorological measurements
The study was carried out in the UNESCO Biosphere Reserve of Yangambi.The reserve covers an area of 230,000 ha and is situated ~ 100 km west of Kisangani, at the banks of the Congo River and just north of the equator in the Tshopo Province (Democratic Republic of Congo), between 0° 49′-0° 51′N and 24° 29′-24° 35′E (Fig. 1).Vegetation in Yangambi consists mainly of moist semi-deciduous rainforest, with patches of moist evergreen  (Nicholson et al. 2018b).The two maps represent the geographic position of the stations with digitized climate data for two-decadal periods (1960-1980 and 2000-2020).Each blue dot in the maps represents a climate station; the yellow dot in the center of each map represents the Yangambi climate station (N 00° 49′ 12.4″; E 024° 27′ 22.6″).Green background color represents closed tropical rainforest cover following Global Landcover Map 2000 (Mayaux et al. 2003) 142 Page 4 of 28 rainforest (in river valleys), transition forest, and agricultural land (Mangaza et al. 2021).Four forest types are found in Yangambi: (i) young secondary forest dominated by pioneer species such as Musanga cecropioides R. Br. ex Tedlie and Macaranga monandra Müll.Arg.; (ii) semi-deciduous mixed forest dominated by long-lived pioneer species such as Pericopsis elata (Harms) Meeuwen; (iii) semi-deciduous mixed forest dominated by shadetolerant species such as Scorodophloeus zenkeri Harms; and (iv) evergreen monodominant forests dominated by Gilbertiodendron dewevrei (De Wild.)J. Léonard or Brachystegia laurentii (De Wild.) Louis ex Hoyle (Luambua et al. 2021).Yangambi has the Af climate according to the Köppen classification with a slightly marked dry season, an average annual precipitation of 1837 mm, and mean annual temperature of 25.1 °C.
The Yangambi meteorological station is managed by the INERA (Institut National pour l'Étude et la Recherche Agronomiques).It is positioned on a large plateau to avoid topographical influences.The site of the station is levelled, maintained regularly, and located far from the immediate influence of trees and buildings.The station records precipitation through an analog precipitation gauge (Casella London), which collects rainfall through a horizontal opening of known cross-section, and stores it in a test tube which is graduated (gauged) from 0 to 13 mm.Minimum (Tmin) and maximum (Tmax) temperature are recorded using two Hg thermometers placed inside a well-aerated housing in an open shelter.Meteorological data at Yangambi is recorded daily at 6 am, 9 am, 12, 15 pm, and 18 pm on paper sheets.Supplementary Fig. 1 shows an example of a manual record and pictures of the Yangambi climate station.
The precipitation gauges deployed at Yangambi have not been changed since their installation in 1928 and thermometers not since 1960.After digitization, data quality checks consisted of a visual inspection of the time series for the detection of possible typos, leading to unrealistic precipitation and extreme temperatures.These checks were complemented by direct comparisons with pictures of the original manual documents.As such, we digitized daily minimum temperature (in °C), daily maximum temperature (in °C), and daily precipitation (in mm), for the period between 1 January 1960 until 31 December 2020.There is a hiatus in the dataset from 1 December 1964 to 30 September 1965, when all activities at Yangambi were suspended due to political unrest.Hence, we excluded the year 1965 from our analysis.For December 1964, we estimated daily values of each variable by interpolating between the equivalent calendar days in 1963 and 1965.
To validate the analog measurements, a state-of-the-art meteorological station (Atmos 41, METER Group, Pullman, WA, USA) has been mounted next to the traditional Yangambi instrumentation in September 2018.This digital meteorological station is part of TAHMO, and data can be monitored in quasi real-time at https:// tahmo.org/ clima te-data/.Digital recordings of daily minimum and maximum temperature as well as precipitation cover one full year of data, from 30 September 2018 until 30 September 2019 (365 days).Daily measurements of analog and digital devices are compared in Supplementary Fig. 2 and Supplementary Fig. 3, suggesting a consistent agreement between analog and digital data.

Monthly meteorological indices
To study climate variability, we considered monthly and yearly meteorological indices developed by the Expert Team on Climate Change Detection Monitoring Indices (ETCCDMI) from the World Climate Research Program of the World Meteorological Organization (see Table 1 for index definitions).Indices were calculated on a monthly and/or annual basis using the daily meteorological measurements.All calculations were performed using the R statistical platform, version 3.2.1 (R Core Team 2017).

Temperature
We first calculated general temperature indices.These are the monthly average of the daily minimum temperature (ATN, in °C), the monthly average of the daily maximum temperature (ATX, in °C), the monthly average of the daily mean temperature (ATM, in °C), and the monthly average diurnal temperature range (DTR = ATX-ATN).Furthermore, we calculated monthly lowest (TXmin) and highest (TXmax) maximum daily temperature and monthly lowest (TNmin) and highest (TNmax) minimum daily temperature.

Rainfall
To quantify rainfall, we considered simple indices such as the monthly total precipitation (PTOT, in mm) and the monthly number of wet (or rainy) days (Rd, in days), with a wet day defined as a day with rainfall ≥ 1 mm.To quantify rainfall intensity, we calculated the monthly average rainfall from wet days, known as the simple day intensity index (SDII, in mm).

Standardized precipitation and evapotranspiration index
We then calculated the standardized precipitation evapotranspiration index (SPEI) at different time scales (3, 6, 12, 24, and 48 months).The SPEI index was developed by Vicente-Serrano et al. (2010).This index is the derivative of the SPI, a commonly used indicator to detect meteorological drought (precipitation deficits) (Mckee et al. 1993;Edwards and Mckee 1997).The fundamental difference between these two indicators is that SPI is calculated based on precipitation, while SPEI is based on the difference between precipitation and potential evapotranspiration (Karam et al. 2022) in including the role of temperature.The procedure to calculate the SPEI index involves a climatic water balance, the accumulation of deficit/surplus at different time scales, and adjustment to a log-logistic probability distribution (Vicente-Serrano et al. 2010).For each decade within our record, we calculated the drought duration as the number of months in drought conditions according to the different SPEI values.

Yearly meteorological indices
Apart from monthly meteorological indices, we also calculated and analyzed yearly indices to fully explore the data.Monthly and yearly indices allow using different statistical approaches to test long-term trends.

Temperature
General temperature indices are the annual average of the daily minimum temperature (ATN, in °C), the annual average of the daily maximum temperature (ATX, in °C), the annual average of the daily mean temperature (ATM, in °C), and the annual average diurnal temperature range (DTR).Furthermore, we calculated annual lowest (TXmin) and highest (TXmax) maximum daily temperature and annual lowest (TNmin) and highest (TNmax) minimum daily temperature (see Table 1 for all index definitions).

Temperature extremes
To capture exceptional hot temperature, we quantified the annual number of warm and extreme warm days as TX95p and TX99p (in days per year), which is the annual count of days when Tmax ≥ 95th or 99th percentile of the 1960-2020 daily measured Tmax record.We quantified the number of warm and extreme warm nights as TN95p and TN99p (in nights year −1 ), which is the annual count of nights when Tmin ≥ 95th or 99th percentile of the 1960-2020 daily measured Tmin record.To capture exceptional cold temperature, we quantified the number of cool and extreme cool days as TX5p and TX1p (in days year −1 ), which is the annual count of days when Tmax ≤ 95th or 99th percentile of the 1960-2020 daily measured Tmax record.We quantified the number of cool and extreme cool nights as TN5p and TN1p (in nights), which is the annual count of nights when Tmin ≤ 95th or 99th percentile of the 1960-2020 daily measured Tmin record.

Rainfall
To quantify rainfall, we considered annual total precipitation (PTOT, in mm), annual dry season precipitation (PDRY, in mm), and annual wet season precipitation (PWET, in mm).Dry season months are defined as those with median precipitation below the overall 1960-2020 median monthly precipitation (142 mm month −1 ).Yangambi exhibits two dry seasons alternating with two rainy seasons (Supplementary Fig. 4).A long dry season extends from December to February, followed by a short rainy season from March to May.Then, a short dry season covering the months of June and July is followed by a long rainy season from August to November.

Rainy days
We first counted the annual total number of wet (or rainy) days (Rd, in days year −1 ), with a wet day defined as a day with rainfall ≥ 1 mm.We then counted the number of wet days in the dry season (RdDRY) as the sum of long (RdDRY djf ) and short (RdDRY jj ) dry season rainy days.Finally, we counted the number of wet days in the wet season (RdWET) as the sum of short (RdWET mam ) and long (RdWET ason ) wet season rainy days.

Rainfall intensity
To quantify rainfall intensity, we calculated the annual average rainfall from wet days, known as the simple day intensity index (SDII, in mm).We also calculated dry (SDII_ DRY) and wet (SDII_WET) season simple day intensity index.To quantify extreme precipitation, we considered indices based on the 95th and 99th percentiles, where R95p is the annual count of days when daily rainfall ≥ 95th percentile of the 1960-2020 daily rainfall record and R99p is the annual count of days when daily rainfall ≥ 99th percentile of the 1960-2020 daily rainfall record.R95p and R99p (in days) represent the number of very wet days and extreme wet days respectively.Derived indices are R95pSUM which is the annual precipitation (in mm) from days when rainfall ≥ 95th mm percentile of the 1960-2020 daily rainfall record and R99pSUM which is the annual precipitation (in mm) from days when rainfall ≥ 99th mm percentile of the 1960-2020 daily rainfall record.R95pSUM and R99pSUM (in mm) represent very wet day intensity and extreme rainfall intensity respectively.Finally, we considered R95pTOT which is the percentage of annual precipitation from days when rainfall ≥ 95 th percentile of the 1960-2020 daily rainfall record and R99pTOT which is the percentage of annual precipitation from days when rainfall ≥ 99th percentile.R95pTOT and R95pTOT (in %) represent the very wet day proportion and extreme rainfall proportion respectively.

Water demand and availability
Potential evapotranspiration (PET) is a widely used indicator of the environmental demand for evapotranspiration (Tadese et al. 2020).Here, we calculate PET using Hamon's equation (Lu et al. 2005) where ATM is the annual average of the daily mean temperature (in °C), Ld is the day length (almost constantly 12 h in Yangambi), and ESAT is the saturation vapor pressure: Furthermore, we calculated net water available (NWA) as follows (Tadese et al. 2020): NWA is a measure for the water available for runoff, tissue storage, soil moisture, and groundwater recharge.

Drought and seasonality indices
We calculated annual maximum climatic water deficit (MCWD), which is a commonly used metric of dry season intensity that has been widely applied in analysis of tropical forest responses to climate change (Aragão et al. 2014;Aguirre-Gutiérrez et al. 2020;Bennett et al. 2021).First, we calculated monthly climatic water deficit (CWD) values for each subsequent series of 12 months (complete years).Monthly CWD estimation begins with the wettest month of the first year in the record and is calculated as monthly evapotranspiration (ET i ) minus monthly precipitation (P i ) (Hubau et al. 2020;Bennett et al. 2021).ET i is calculated from monthly PET i and P i using a simple model calibrated using the global eddy covariance flux network and remote sensing data (Sun et al. 2011;Sun 2013;Fang et al. 2016): Then, CWD values for the subsequent 11 months were calculated recursively as: where negative CWD i values were set to zero (no drought conditions).This procedure was repeated for each subsequent complete 12 months.We then calculated the annual MCWD as the largest monthly CWD value for every complete year within the record.Larger MCWD indicates more severe water deficits.
Finally, we calculated the seasonality index (SI) for each year in the record as (Walsh and Lawler 1981): where R i is the total annual precipitation for the particular year i and X in is the actual monthly precipitation for month n.The seasonality index allows classifying climate based on rainfall seasonality using the following categories: SI < 0.19: precipitation spread throughout the year SI = 0.20-0.39,precipitation spread throughout the year but with a definite wetter season; SI = 0.40-0.59,rather seasonal climate with a short dry season; SI = 0.60-0.79,seasonal climate; SI = 0.80-0.99,marked seasonal climate with a long dry season; SI = 1.00-1.19,most precipitation in < 3 months

Long-term linear trends
To analyze long-term linear trends over the entire 1960-2020 dataset, we fitted linear mixed models on the daily and monthly meteorological indices from the Yangambi climate station, following Bush et al. (2020).For temperature indices (ATN and ATX) and SDII, we used a linear mixed effect model via REML (lmer command) from the lme4 R package (Bates et al. 2015).For rainfall (PTOT, Rd) indices, we fitted compound Poisson generalized linear mixed models (cpglmm command) from the cplm R package (Zhang 2013).The cpglmm models are exponential dispersion models using true likelihood-based inferential procedures (i.e., the Laplace approximation and the adaptive Gauss-Hermite quadrature).They are recommended for (positive and continuous) daily or monthly rainfall data which typically contain many zeros (Hasan and Dunn 2010).The lmer and cpglmm models regress the response variable against the year as predictor and include a random effect to account for seasonality and the hierarchical structure of the data.In the models using daily data, day-of-the-year was included as a random effect, that is, by assuming that the intercept can vary randomly among the days of the year.In the models on monthly data, we used month-of-the-year as a random effect.To evaluate the meaningfulness of the slope in the models, we compared model AIC with that of an intercept-only model (representing no long-term If AIC of the model with slope is lower than that of the intercept-only model, we consider the slope significant.To evaluate seasonal trends, models were also run on dry season and wet season subsets of the daily and monthly rainfall indices (PDRY, PWET, SDII_DRY, SDII_WET, RdDRY, RdWET).
Additionally, we performed the nonparametric Kendall's tau test on the yearly meteorological indices to determine if the monotonic long-term linear trend in the time series was significant.We performed Kendall's test on smoothed time series using locally weighted regression (lowess filter) following Peterson et al. (2008).Kendall's tau test has been widely used to compute trends in hydrometeorological series because it is robust to the effect of outliers in the series (Peterson et al. 2008;Aguilar et al. 2009;Chaney et al. 2014).We performed smoothing using the lowess function of the gplot package (Cleveland 1979) and we performed Kendall's test using the kendall-TrendTest function of the EnvStats package in R (Kendall 1975).

Piecewise linear trends
To analyze short-term trends within the 1960-2020 dataset, we fitted linear breakpoint models on the yearly meteorological indices from the Yangambi climate station.We first smoothed each annual index using the rollmean function of the zoo package in R with a 5-year moving-average window.Then, we parameterized a simple linear regression model on the full time series using the lm function in R. Finally, for each index, we parameterized single breakpoint and double breakpoint models using the segmented package in R (Muggeo 2008).For the single breakpoint models, we specified 1990 as the initial breakpoint value; for the double breakpoint models, we specified 1980 and 2000 as initial breakpoints.Starting from these, the segmented function then estimates best-fit breakpoints using the bootstrap restarting algorithm (Wood 2001).
To select the best-fit model for each index, we considered the R 2 value and the P-value of the Davies test (using the Davies test function of the segmented package).The Davies test allows to evaluate whether a model would be significantly improved when adding an extra breakpoint (Davies 2002).Specifically, if the P-value of the Davies test (P davies ) for a certain model is < 0.05, then the model would be improved by allowing another breakpoint.
For each index, we started by considering the simplest model.We selected the nonbreakpoint model if its P davies > 0.05 and R 2 > 0.2.Then, if the non-breakpoint model was rejected based on one or both of these criteria, we selected the single breakpoint model if its P davies > 0.05 and R 2 > 0.2, and if the model converged.Then, if the single breakpoint model was rejected based on these criteria, we selected the double breakpoint model if it converged.
Additionally, we considered the meaningfulness of the breakpoints for selected double breakpoint models.If the two breakpoints were less than 5 years apart, the double breakpoint model was rejected in favor of the single breakpoint model.Supplementary Table 1 shows the slope(s), intercept(s), P-value(s), and breakpoint(s) for the selected model for each index.Slopes with a P-value < 0.05 are considered significant.

Comparison with gridded datasets
To compare the long-term rainfall trends observed in the Yangambi record, we used two gridded monthly rainfall datasets.We extracted Yangambi rainfall data from the Global Precipitation Climatology Centre (GPCC) at 1° resolution (~ 3025 km 2 ) (Schneider et al. 2014).We used the GPCC's most accurate in situ land surface precipitation analysis product by combining the Full Data V7 Product , based on quality-controlled data from 67,200 meteorological stations worldwide that feature record durations of 10 years or longer, with the V6 monitoring product (2011)(2012)(2013)(2014), based on quality-controlled data from 7000 stations.This resulted in a GPCC rainfall record covering 1960-2014.Second, we extracted Yangambi rainfall data from the 0.25° resolution (~ 757 km 2 ) data from the Tropical Rainfall Measurement Mission (TRMM product 3B43 V7) from 1998 to 2020 (Huffman et al. 2007) which provides up-to-date satellite data for the whole tropical realm.Because TRMM only starts at 1998, we combined it with GPCC, following Hubau et al. (2020) and Bennett et al. (2021).The fit for the overlapping time period (1998-2014) was used to correct any systematic difference between GPCC and TRMM (Hubau et al. 2020).This resulted in a TRMM-GPCC rainfall record covering 1960-2020.GPCC and TRMM are the most widely used benchmark datasets to validate other gridded rainfall datasets (Nicholson et al. 2018b;Sun et al. 2018;Harris et al. 2020;Igri et al. 2022;Ndehedehe and Agutu 2022) and to analyze tropical forest responses to climate change (Zhou et al. 2014;Hubau et al. 2020;Bennett et al. 2021).
Comparing observed (Yangambi) and extracted (CRU, GPCC, TRMM) data is a symmetric problem, meaning that we have a pair of Y variables and we want to see how they are related to each other.We therefore performed standardized major axis estimation using the sma command of the smart package in R (Warton et al. 2012).We used robust estimation to ensure that inferences are valid in the presence of outliers, which exist in all these databases.

Long-term linear trends
We find a significant long-term increase in minimum (ATN), mean (ATM), and maximum (ATX) daily temperature, at 0.14, 0.18, and 0.21 °C per decade respectively, using the daily Yangambi climate station data (Table 2).This is confirmed using the monthly (Table 3) and yearly (Table 4) meteorological indices.We also find a long-term significant increase in (extreme) warm nights and (extreme) warm days and a significant decrease in (extreme) cool days (Table 4).All temperature indices reached record breaking values during the last few years (Fig. 2).ATX and ATN reached absolute maxima in respectively 2016 and 2020.ATX increased by 1.6 °C (from 29.2 to 30.8 °C) and ATN by 1.4 °C (from 19.8 to 21.2 °C) since the 1960s.Similarly, the number of warm days (TX95p) and warm nights (TN95p) also peaked in 2016 and 2020, at 46 and 118 days/ nights per year respectively.To compare, 1960 only had 2 warm days and 4 warm Page 11 of 28 142 nights.Finally, the number of cool days (TX5p) decreased to a minimum of 11 per year in 2016 and the number of cool nights (TN5p) decreased to zero in 2019.
We do not find a significant long-term change in total annual precipitation (PTOT) using the daily and monthly data, although the difference in AIC between the slope and the no-slope models is very small (ΔAIC < 2; Table 2 and Table 3).However, Kendall's tau test using annual data smoothed by the lowess model shows contrasting longterm slopes between dry season precipitation (PDRY; decreasing at − 5 mm per year per decade; P = 0.009) and wet season precipitation (PWET; increasing at 17 mm per year per decade; P < 0.001) (Table 4).Declining dry season precipitation and increasing wet season precipitation suggest an increasing precipitation seasonality trend in Yangambi (Fig. 2, Table 4).This is confirmed by long-term trends in net water availability (NWA), a metric that accounts for both precipitation and temperature.NWA is declining significantly in the dry season (NWA_DRY, at − 10.56 mm per year per decade) and increasing significantly in the wet season (NWA_WET, at 8.5 mm per year per decade) (Table 4).
Using the monthly and yearly data, we find a significant long-term decrease in the yearly number of rainy days (Rd, at − 2.8 days per year per decade) and a significant increase in rainfall intensity, with a simple day intensity index (SDII) increasing at 0.48 mm per day per decade (Table 3, Table 4).This is confirmed by significant increases in very wet day frequency, intensity, and proportion (R95p, R95pSUM, R95pTOT; see Table 1 4).Thick colored lines represent best-fit linear regression models (non-breakpoint models in gray, single breakpoint models in orange, double breakpoint models in blue).Regression lines are solid if significant or dashed if non-significant.Index definitions are given in Table 1; slopes and P-values of breakpoint regression models are given in Supplementary Table 1 142 Page 12 of 28 definitions) and extreme rainfall frequency, intensity, and proportion (R99p, R99pSUM, R99pTOT).The number of rainy days is decreasing both in the dry and wet seasons (Rd_ DRY and Rd_WET), and rainfall intensity is increasing in both dry and wet seasons (SDII_ DRY and SDII_WET).This suggests that precipitation intensity is increasing throughout the year.

Piecewise trends: temperature
Selected simple linear models vary among temperature indices.For ATX, the nonbreakpoint model is selected, while for ATN and ATM, the double breakpoint model was selected (Supplementary Table 1, Fig. 2).Using the simple regression model, ATX increased significantly (P < 0.001) and continuously over the full 1960-2020 period, at 0.22 °C per decade, which is comparable to the linear mixed models on the daily (Table 2) and monthly (Table 3) data.ATN and ATM both initially show an increase, then a significant decrease between 1982 and 1993, and finally a very steep, nearly three-decade increase from 1993 until 2020, at 0.47 and 0.34 °C per decade for ATN and ATM respectively (all trends significant with P < 0.001) (Fig. 2).
Selected simple linear models also vary among night-time extreme temperature indices (Supplementary Table 1, Fig. 2).Single breakpoint models were selected for the number of cool nights (TN5p) and extreme cool nights (TN1p).These indices increased in the first half of the record, but decreased significantly at the end of the record (all trends significant with P < 0.001).Especially, the number of extreme cool nights (TN1p) decreased spectacularly since 1994, at − 10.8 nights per decade.In contrast, the number of extreme warm nights (TN99p) increased continuously and significantly (P < 0.001) throughout the 1960-2020 record, at 1.8 nights per decade.A double breakpoint model was selected for the number of warm nights (TN95p), with a spectacular increase since 1994, at 2.1 nights per decade (Fig. 2).Selected simple linear models are remarkably comparable among extreme day-time temperature indices.For each of the indices, the non-breakpoint model was selected.The numbers of cool days (TX5p) and extreme cool days (TX1p) are significantly decreasing throughout the record, at − 1.8 and − 0.6 days per decade respectively (P < 0.001 for both indices; Supplementary Table 1).In contrast, we observe a consistent increase in warm days (TX95p) and extreme warm days (TX99p), at 4.2 and 1.3 days per decade respectively (P < 0.001 for both indices; Supplementary Table 1).
Furthermore, trends are remarkably different among the seasons.The wet season is characterized by a significant long-term increase in total precipitation (PWET, P = 0.013) and rainfall intensity (SDII_WET, P < 0.001), along with a decrease in rainy days (Rd_WET, P < 0.001) (Supplementary Table 1, Fig. 3).Also, very wet day frequency (R95p) and proportion (R95pTOT) are increasing throughout the record (Supplementary Table 1).In contrast, the dry season is characterized by two-decadal trends that are very similar to full-year trends.Specifically, dry season precipitation (PDRY), number of rainfall days (Rd_DRY), and simple day intensity index (SDII_DRY) are decreasing from 1960 to 1979 (breakpoint varies between 1975 and 1980 among indices).Then, all indices are increasing from 1980 to 1998 and decreasing again from 1999 to 2020 (Supplementary Table 1; all P < 0.001).

Piecewise trends: drought and seasonality
The short-term trends observed in the full-year and dry season records are also reflected in the standardized precipitation evapotranspiration index (SPEI) (Vicente-Serrano et al. 2010), which is mathematically similar to the standardized precipitation index (SPI), but it includes the role of temperature.We counted the frequency of dry months by decade in the Yangambi precipitation record, using well-defined drought categories (Supplementary Table 2 and Supplementary Fig. 5).We find that the 1960s had a very low number of dry months in all categories (12%), followed by the 1970s which had a slight increase in dry months (15%); the 1980s and 2010s had a moderately high number of dry months Fig. 6 Comparison between observed records of precipitation and temperature from the Yangambi meteorological station with those from existing gridded datasets.Scatterplots show observed (x-axis) monthly precipitation and mean temperature observations from the Yangambi record (1960-2020) versus extracted observations for the same months from GPCC (green, 1960-2014), TRMM (red, 1998(red, -2019(red, ), and CRU (blue, 1960(red, -2019)).Colored lines represent standardized major axis estimation models; correlation coefficients (R) and slopes are shown in the top right corner.The 1:1 line is shown in black Page 15 of 28 142 respectively (17% in each of the two decades); and finally, the 1990s and 2000s had a very high number of dry months (19% in each of the two decades).This indicates a wetting trend from 1960 to 1970, an increasing trend to moderate drought observed in the 1980s and 2010s, and then a trend to very high drought from 1990s to 2000s.
A drought index based on both mean annual temperature (ATM) and total annual precipitation (PTOT) is net water availability (NWA) (Tadese et al. 2020).NWA is the balance between total water availability (precipitation, PTOT) and water demand (expressed as potential evapotranspiration, PET, which is calculated as a function of ATM).We find that both full-year and dry season patterns in NWA are very similar to those of PTOT (Supplementary Table 1, Fig. 4).Dry season NWA was negative (indicating water shortage) during most periods in the record, except the 1960s and 1990s (Supplementary Table 1).In contrast, wet season NWA was always positive and increased throughout the 1960-2020 record (P = 0.044), while dry season NWA decreased, increased, and decreased again during subsequent two-decadal periods (all P < 0.001).
The rainfall (Fig. 3) and drought (Fig. 4) trends show substantial differences between the dry and wet season.This suggests shifts in precipitation seasonality across the record.To test this, we finally assessed trends in two rainfall seasonality indices: the Walsh and Lawler seasonality index (SI) and the maximum climatological water deficit (MCWD) (Fig. 4, Supplementary Table 1).For both indices, the double breakpoint model was retained, again with breakpoints close to those in the PTOT model, dividing the record in three nearly two-decadal periods.We find that seasonality was significantly increasing, then decreasing, and finally increasing again (all P < 0.001).The seasonality index indicates that the precipitation regime in Yangambi switched between two different categories over the 1960-2020 period.In the wettest decades (1960s and 1990s), SI reached minima of 0.28 and 0.29 (in 1962 and 1996), which indicates that "precipitation was spread throughout the year" (Walsh and Lawler 1981).In the driest decade (1980s), SI reached a maximum of 0.55, which indicates a "rather seasonal precipitation regime."The last two decades also tended to be rather seasonal, with an average SI of 0.40 and 0.44 and a maximum SI of 0.53 (2012 and 2018).
The annual average mean daily temperature record extracted from CRU shows a nonsignificant long-term trend (P = 0.298), in contrast with a significant warming trend in the observed Yangambi record (Table 2, Table 3, Table 4).Piecewise regression in the CRU record does show a significant warming trend between 1975 and 2008, but a drop during the last decade (Fig. 5), which contrasts with a consistent steep warming trend in the observed record since 1995 (Fig. 2).
Slopes of standardized major axis estimation models are close to 1, suggesting good agreement between the observed and extracted monthly precipitation values, for GPCC (slope = 0.79), TRMM (slope = 0.91), and CRU (slope = 0.96) (Fig. 6).This indicates that   high observed values generally correspond to high extracted values.However, moderate correlation coefficients (R = 0.62 for GPCC, R = 0.53 for TRMM) suggest that absolute monthly precipitation values are sometimes very different between the extracted and observed data, although differences exist in both ways (sometimes observed values are higher than extracted values and vice versa).For CRU, the correlation coefficient is high (R = 0.79) but extracted values from CRU tend to be systematically higher than observed values in Yangambi, in the order of ~ 0.5 °C.

Consistent warming trend
Long-term continuous warming over the entire 1960-2020 period in Yangambi (Table 2, Table 3, Table 4) is consistent with the overall temperature change in Western and Central Equatorial Africa, where ATM has increased by about 1 °C in the past 70 years (Dezfuli and Dezfuli 2017).Increasing numbers of (extreme) warm days and nights and decreasing numbers of (extreme) cool days and nights are consistent with reported changes in temperature extremes throughout Central Africa since 1955 (Aguilar et al. 2009;Chaney et al. 2014).

Contrasting regional long-term precipitation trends
The non-significant long-term precipitation trend in Yangambi (Table 2, Table 3) corroborates a non-significant trend over the 1979-2005 period in the same region (Chaney et al. 2014).The short-term decline in PTOT observed in Yangambi between 1960 and 1980 corresponds to a reported continuous rainfall decline over most of the Central African continent, as indicated by analysis of all available gauge stations during that time period (Nicholson et al. 2018a).One of the wettest years in the overall central African gauge station network is 1961, whereas 1983 is among the driest years (Nicholson et al. 2018b).
Table 2 Comparisons of linear regression model results to test for long-term trends in temperature and rainfall indices using "daily" observations from the Yangambi meteorological station .We used linear mixed models (lmer) for temperature indices and compound Poisson generalized linear mixed models cpglmm) for precipitation.For each model, day-of-the-year was included as a random intercept to account for seasonality and the hierarchical structure of the data.For each index, we ran a model without slope (representing no long-term change) and a model with a slope.Here, we present AIC of each model and the difference (ΔAIC) between both models.For each index, the model with lowest AIC is retained (indicated in bold) The short-term increase in PTOT in Yangambi between 1981 and 1999 is also prominent in the interpolated central African gauge dataset, where 1997 is among the wettest years (Nicholson et al. 2018b).Finally, the declining trend in PTOT in Yangambi between 2000 and 2020 corroborates previous reports of significant drying in the Congo Basin during the past two decades (Malhi and Wright 2004;Asefi-Najafabady and Saatchi 2013;Diem et al. 2014;Zhou et al. 2014;Tsalefac et al. 2015;Dezfuli and Dezfuli 2017;Nicholson et al. 2018bNicholson et al. , 2019;;Mabrouk et al. 2022).
We note substantial similarities between our observed piecewise trends in PTOT and piecewise trends in published Congo River's annual discharge data (Mahé 1995;Laraque et al. 2013Laraque et al. , 2020;;Tsalefac et al. 2015;Bola et al. 2022; Ndehedehe and Agutu 2022) (Supplementary Fig. 6).Nevertheless, we also note substantial differences in long-term trends, with significant decline in Congo River's discharge since 1960 versus no significant change in PTOT in Yangambi.This is due to the historically high Congo River discharge levels in the 1960s, which must be related to anomalously high precipitation in regions within the Congo Basin other than the Yangambi region.Similarly, the precipitation pattern in Yangambi contrasts with the strong negative continuous drying trend at − 75 mm per decade over the 1984-2018 period in the Lopé National Park (Gabon) (Bush et al. 2020).Furthermore, it contrasts with strong negative rainfall trends over the northern and southern sector of the Congo Basin (Dezfuli and Dezfuli 2017) and a significant negative trend at − 31 mm per year per decade (Aguilar et al. 2009(Aguilar et al. ) over the 1955(Aguilar et al. -2003 period in the entire central Congo Basin.This confirms strong differences in long-term precipitation trends among central African regions (Dezfuli and Dezfuli 2017).Indeed, recent analysis of gridded products (GPCC and CRU) showed remarkable spatial variability in precipitation change across the entire Congo Basin, with two-decadal trends (decreasing precipitation in 1961-1980, increase in 1981-2000, decrease in 2001-2013) for the Yangambi area similar to those observed in our record (Ndehedehe and Agutu 2022).
Table 3 Comparisons of linear regression model results to test for long-term trends in temperature and rainfall indices using "monthly" indices calculated from the Yangambi meteorological station data .We used linear mixed models (lmer) for temperature and SDII indices.We used compound Poisson generalized linear mixed models (cpglmm) for precipitation and Rd.For each model, month-of-the-year was included as a random intercept to account for seasonality and the hierarchical structure of the data.For each index, we ran a model without slope (representing no long-term change) and a model with a slope.Here, we present AIC of each model and the difference (ΔAIC) between both models.For each index, the model with lowest AIC is retained (indicated in bold)

Increasing precipitation "seasonality"
The increasing seasonality index (SI) and increasing maximum climatic water deficit (MCWD) over the 1998-2020 period (Fig. 4) suggest recent intensification of the seasonal rainfall cycle in Yangambi, with a drier dry season and wetter rainy season.This is confirmed by an increasing trend in wet season precipitation and a decreasing trend in dry season precipitation (Fig. 3).This trend towards intensifying precipitation seasonality in central Africa was recently confirmed by a combination of farmer perceptions, satellitebased estimates, and ground-based station measurements (Salerno et al. 2019).
Our observation of increasing long-term wet season precipitation in Yangambi is consistent with an observed spectacular increase in the annual lightning flash density and the number of stormy days in the central Congo Basin between 2005 (114 stormy days per year) and 2013 (189 stormy days per year) (Soula et al. 2016).Furthermore, our observation of a drying dry season over the last two decades (1999-2020) is consistent with the reported significant increase in dry season length in the Congo Basin, by 6.4-10.4days per decade, attributed to an earlier dry season onset and a delayed dry season end (Jiang et al. 2019).Finally, it is consistent with an observed long-term  increase in the areal extent and intensity of thunderstorms over the Congo Basin, causing reduced moisture availability in the Lower troposphere over the Congo Basin (Raghavendra et al. 2018).
Opposing dry and wet season long-term precipitation trends in Yangambi (Table 4) are also comparable to the situation in Lopé (Gabon), where total precipitation decline was driven by a decline in rainfall in the long dry season (June-July-August-September) but not in the wet seasons.This suggests that increases in intra-annual rainfall variability, more specifically increased contrast between the wet and dry seasons, are a widespread phenomenon in Central Africa (Bush et al. 2020).

Increasing precipitation "intensity"
Our results (Table 4) also corroborate previously reported decrease in rainy days (at − 0.67 days per year per decade), a long-term increase in SDII (at 0.06 mm per day per decade), and increasing very wet day intensity in Central Africa (at 13.6 mm per year per decade) (Aguilar et al. 2009;Chaney et al. 2014).The combination of decreasing number of rainy days (Rd), increasing rainfall intensity (SDII), and increasing very wet day intensity (R95pSUM) confirms model outputs showing that wet extreme events are expected to become more frequent and more severe in Equatorial Africa (Fotso-Nguemo et al. 2018;Kendon et al. 2019;Karam et al. 2022).

Impact of climate variability in Yangambi
Continued warming and increasing precipitation seasonality and intensity in Yangambi (2000-2020) have a significant impact on local livelihoods.A recent study showed that the vast majority of farmers in Yangambi not only perceived changes in temperature and precipitation, but also experienced reduced crop yields and the emergence of new weed species and pests (Mangaza et al. 2021).Increasing impact of climate variability will particularly require implementation of climate-smart agriculture (Laux et al. 2010;Mangaza et al. 2021).Ongoing divergence of dry season (getting drier) and wet season (getting wetter and more intense) climate in Yangambi urges for a dynamic adaptation approach, where farmers introduce small management adjustments based on observed changes in Page 23 of 28 142 productivity in previous years, rather than following prescribed adjustment protocols (Meza and Silva 2009).The tangible impact of climate variability in Yangambi might also call for an increased focus on agroforestry, which is an example of a livelihood system combining various benefits and combining mitigation and adaptation effects (Batsi et al. 2020;Tschora and Cherubini 2020).
Finally, a shift towards a hotter and more seasonal climate may also have important consequences for Congo Basin forest composition and functioning.Several authors report ongoing shifts in species composition due to long-term drying in West Africa (Fauset et al. 2012;Aguirre-Gutiérrez et al. 2020) and in the drier parts of the Amazon (Esquivel- Muelbert et al. 2019).Dying wet-affiliated trees in the forest are replaced by more drought-tolerant (deciduous) species.No analysis using ground-based long-term forest inventory data has been conducted for Yangambi or any other region in the central Congo Basin yet, but analysis of satellite images revealed a long-term decline of Congo rainforest greenness, associated with a shift of species composition towards a larger proportion of deciduous species in the canopy (Zhou et al. 2014).Even though no shifts in forest composition have been reported for Yangambi through direct ground-based observations, increasing precipitation seasonality and decreasing net water availability are likely to lead to drought stress in its forests (Karam et al. 2022), which could push these ecosystems beyond their limits (James et al. 2013).Combined long-term and systematic monitoring of both biodiversity and climate in Yangambi and, more broadly, in the Congo Basin is indispensable to develop early warning systems for ecosystem degradation and collapse.
Our comparison analysis shows that remarkable discrepancies exist between observed climatic variables and those extracted from gridded products.Differences are found both in long-term trends and in absolute monthly values.Recent analysis showed that the correlation between observed and extracted (from TRMM) precipitation values is lowest in the wettest areas of central Africa (Igri et al. 2022).These results, combined with the exceptionally poor station network in the wet central Congo Basin, call for more intense efforts in climate monitoring.

Conclusions
Our newly digitized data from Yangambi represent a unique climate record in the central Congo Basin, a remarkably under sampled and hence enigmatic region in climate research.Our results from Yangambi corroborate recently reported long-term increase in temperature and temperature extremes throughout Central Africa since at least 1960.Short-term trends in the Yangambi records suggest that the increase in temperature and temperature extremes accelerated since the early 1990s.Our results show no long-term change in precipitation and drought indices in Yangambi.Nevertheless, we find long-term increasing intra-annual rainfall variability, both in terms of rainfall amount (decreasing the in dry season, increasing in the wet season) and in rainfall intensity (especially increasing in the wet season).Furthermore, we find highly significant short-term trends in precipitation and drought indices, dividing the record in three almost equally long two-decade periods.The first period (1960)(1961)(1962)(1963)(1964)(1965)(1966)(1967)(1968)(1969)(1970)(1971)(1972)(1973)(1974)(1975)(1976)(1977)(1978) is characterized by decreasing precipitation, decreasing net water availability, and no change in rainfall intensity.The second period  is characterized by increasing precipitation, increasing net water availability, and increasing rainfall intensity.The third period  is characterized by decreasing precipitation, decreasing net water availability, and decreasing rainfall intensity.These short-term trends are prominent in the dry season, but not in the wet season, which shows a significant longterm increase in precipitation, net water availability, and rainfall intensity.The seasonality index confirms that the precipitation regime in Yangambi switched between two different categories over the 1960-2020 period, from "rather seasonal" (1960)(1961)(1962)(1963)(1964)(1965)(1966)(1967)(1968)(1969)(1970)(1971)(1972)(1973)(1974)(1975)(1976)(1977)(1978) to "precipitation spread throughout the year"  and back to "rather seasonal" .
Our results are in line with previous reports of a drying trend in the dry season in central Africa during the last two decades, based on available interpolated gridded climate datasets.Ongoing warming and increasing precipitation seasonality and intensity in Yangambi, with increasing extreme temperature events, drier dry seasons, and wetter rainy seasons, already have a significant impact on crop yields in Yangambi.This calls for urgent development of climate-smart and dynamic agriculture and agroforestry systems.To this end, systematic digitization and climate recording in the Congo Basin will be critical, especially because existing satellite and interpolated climate products are not entirely representative.
Page 25 of 28 142 Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Fig. 1
Fig. 1 Evolution and geographic position of climate stations with available precipitation data in and around the Congo Basin (outlined in thick black line).The top panel represents the evolution of the number of stations with digitized climate data for the last six decades (1960-2020)(Nicholson et al. 2018b).The two maps represent the geographic position of the stations with digitized climate data for two-decadal periods(1960-1980 and 2000-2020).Each blue dot in the maps represents a climate station; the yellow dot in the center of each map represents the Yangambi climate station (N 00° 49′ 12.4″; E 024° 27′ 22.6″).Green background color represents closed tropical rainforest cover following Global Landcover Map 2000(Mayaux et al. 2003) PET = 0.1651 × Ld × 216.7 × ESAT∕(ATM + 273.3)ESAT = 6.108 × EXP(17.26939× ATM∕(ATM + 237.3))NWA = PTOT − PET 142 Page 8 of 28

Fig. 2
Fig. 2 Patterns of temperature and temperature extreme indices in Yangambi, measured over the period 1960-2020.Each open circle represents one yearly record; the thin gray line shows the smoothed 5-year moving-average time series; the thick solid black line shows the lowess smoother (slopes and P-values in Table4).Thick colored lines represent best-fit linear regression models (non-breakpoint models in gray, single breakpoint models in orange, double breakpoint models in blue).Regression lines are solid if significant or dashed if non-significant.Index definitions are given in Table1; slopes and P-values of breakpoint regression models are given in Supplementary Table1

Fig. 3
Fig. 3 Patterns of precipitation (dry and wet season), rainy days, and rainfall intensity in Yangambi, measured over the period 1960-2020.Figure details as in Fig. 2

Fig. 4 142
Fig. 4 Patterns of water demand (PET), net water availability (NWA), maximum climatic water deficit (MCWD), and seasonality (SI) in Yangambi, measured over the period 1960-2020.Figure details as in Fig. 2 cool nights-annual nr of nights when Tmin ≤ 1th percentile of 1960-2020 TN5p Number of cool nights-annual number of nights when Tmin ≤ 5th percentile of 1960-2020 TN95p Number of warm nights-annual number of nights when Tmin ≥ 95th percentile of 1960-2020 TN99p Number of extreme warm nights-annual nr of nights when Tmin ≥ 99th percentile of 1960-2020 TX1p Number of extreme cool days-annual nr of days when Tmax ≤ 1th percentile of 1960-2020 TX5p Number of cool days-annual number of days when Tmax ≤ 5th percentile of 1960-2020 TX95p Number of warm days-annual number of days when Tmax ≥ 95th percentile of 1960-2020 TX99p Number of extreme warm days-annual nr of days when Tmax ≥ 99th percentile of 1960i.e., annual total number of wet days (rainfall ≥ index: average rainfall from wet days (rainfall ≥ 1 mm) SDII_DRY Dry season simple day intensity index SDII_WET Wet season simple day intensity index R95p Very wet day frequency, i.e., annual number of days when rainfall ≥ 95th percentile of the 1960-2020 daily rainfall record R99p Extreme rainfall frequency, i.e., annual number of days when rainfall ≥ 99th percentile of the 1960-2020 daily rainfall record R95pSUM Very wet day intensity, i.e., total annual precipitation from days when rainfall ≥ 95th mm percentile for 1960-2020 R99pSUM Extreme rainfall intensity, i.e., total annual precipitation from days when rainfall ≥ 99th mm percentile for 1960-2020 R95pTOT Very wet day proportion, i.e., percentage of annual precipitation from days when rainfall ≥ 95th percentile of 1960-2020 R99pTOT Extreme rainfall proportion, i.e., percentage of annual precipitation from days when rainfall ≥ 99th percen-

Table 4
Comparisons of linear regression model results to test for long-term trends in temperature and rainfall indices using "yearly" indices calculated from the Yangambi meteorological station data.For each index, we performed a nonparametric test for a monotonic trend based on Kendall's tau statistic, using yearly data