<figure><img src="image_1.png" /></figure>

HAL
open science

Shifting in the global flood timing

Gonghuan Fang, Jing Yang, Zhi Li, Yaning Chen, Weili Duan, Charles Amory,

Philippe de Maeyer

<figure><img src="image_2.png" /></figure>

To cite this version:

Gonghuan Fang, Jing Yang, Zhi Li, Yaning Chen, Weili Duan, et al. Shifting in the global flood timing. Scientific Reports, 2022, 12 (1), pp.18853. 10.1038/s41598-022-23748-y. hal-04924158

HAL Id: hal-04924158

https://hal.science/hal-04924158v1

Submitted on 4 Feb 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_4.png" /></figure>

Check for updates

OPEN

# Shifting in the global flood timing

Gonghuan Fang<fnref n="1" />,²,³, Jing Yang⁴, Zhi Li<fnref n="1" />, Yaning Chen<fnref n="1" />✉, Weili Duan<fnref n="1" />, Charles Amory⁵ & Philippe De Maeyer⁶,⁷

Climate change will have an impact on not only flood magnitude but also on flood timing. This paper studies the shifting in flood timing at 6167 gauging stations from 1970 to 2010, globally. The shift in flood timing and its relationship with three influential factors (maximum 7-day precipitation, soil moisture excess, and snowmelt) are investigated. There is a clear global pattern in the mean flooding date: winter (Dec–Feb) across the western Coastal America, western Europe and the Mediterranean region, summer (Jun–Aug) in the north America, the Alps, Indian Peninsula, central Asia, Japan, and austral summer (Dec–Feb) in south Africa and north Australia area. The shift in flood timing has a trend from – 22 days per decade (earlier) to 28 days per decade (delayed). Earlier floods were found extensively in the north America, Europe and northeast Australia while delayed floods were prevailing in the Amazon, Cerrado, south Africa, India and Japan. Earlier flood timing in the north America and Europe was caused by earlier snowmelt while delayed extreme soil moisture excess and precipitation have jointly led to delayed floods around the monsoon zone, including south Africa, India and Japan. This study provides an insight on the shifting mechanism of flood timing, and supports decisions on the global flood mitigation and the impact from future climate change.

Floods are one of the most dangerous climate-related disasters, and climate change has altered the distribution, intensity and timing of floods worldwide<fnref n="1" />⁻⁷.

There are intensive studies on flooding, most of which focus on the historical trends in flood magnitude and intensity and could not reach a consistent signal in flood magnitude change on global scale⁸⁻¹⁴. The influencing factors in flood magnitude changes have also been investigated extensively including climatic and human factors, e.g., precipitation, soil moisture, snow, dam construction, land-use change, river training¹⁴⁻¹⁸. Recently, there have been more studies on flood timing³,¹⁹⁻²⁵. The spatial pattern of mean flooding date has been identified in Europe³,¹⁹,²²,²⁶, north America²³,²⁷⁻²⁹, Australia³⁰, Brazil³¹, Africa³² or globally¹⁰,²¹.

Shifts in flood timing is often used as a monitor to interpret the mechanisms that cause floods²⁰,³³. The underlying mechanism under exploitation includes the atmospheric conditions and others, and the atmospheric factors are much more influential than the catchment and river system factors¹¹,²³,³⁴. For a specific flood event, the influential factors are well documented in terms of hydro-meteorological or circulation aspect³⁵,³⁶. On the regional or continent scales, these influential factors include short-term heavy rainfall, rainfall on saturated soil, snow melt water, which individually or jointly stimulate the flood events³,³⁷⁻³⁹. For example, the mean flood date has altered especially in the high latitude north hemisphere due to warming induced earlier snowmelt³,¹⁹. In Europe, earlier spring snowmelt floods in northeastern Europe, delayed winter flood around the North Sea and parts of the Mediterranean coast and earlier winter flood in Western Europe were observed and most annual floods are caused by subextreme precipitation with high antecedent soil moisture, not by annual peak rainfall³,¹⁹. In the Upper Austria, one- or seven-day extreme precipitation is usually a better covariate for variations of the flood frequency curve than precipitation on longer time scales³⁴. In the United States, the flood mechanisms can be summarized as daily precipitation, weekly precipitation, precipitation excess and snowmelt plus rainfall²⁰. In Australia, changes in antecedent soil moisture modulate flood seasonality and the driving factors are different in small and large floods³⁰. Some studies also attribute flood time change to large scale circulations¹²,⁴⁰. Most of these studies focus on changes in flood timing on a continental or basin scales, the research on the spatial pattern of shifts in flood timing under different climatic conditions on the global perspective is still insufficient³⁰,³³.

The objective of this study is to quantify the shifts in flood timing during the past few decades and to identify its influencing factors at a global scale considering the climatic zones. Based on the map of flood timing changes,

<footnote n="1">¹State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, 818 South Beijing Road, Urumqi 830011, China. ²Sino-Belgian Joint Laboratory for Geo-Information, Urumqi 830011, China. ³Xinjiang Key Laboratory of Water Cycle and Utilization in Arid Zone, Urumqi, China. ⁴National Institute of Water and Atmospheric Research, Christchurch 8000, New Zealand. ⁵CNRS, Institut des Géosciences de l'Environnement, University Grenoble Alpes, 38100 Grenoble, France. ⁶Department of Geography, Ghent University, 9000 Ghent, Belgium. ⁷Sino-Belgian Joint Laboratory for Geo-Information, 9000 Ghent, Belgium. ✉email: chenyn@ms.xjb.ac.cn</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 1.** The spatial pattern of flood timing during 1970–2010: (a) mean flood date, displayed as circular for each hotspot; (b) latitude averaged mean flood date; (c) flood concentration index around corresponding mean occurrence date; (d) latitude averaged flood concentration. Both the mean and concentration of flooding date are aggregated in 19 hotspots based on the Köeppen Climate Classification System<sup>43</sup>. The names of these hotspots are listed below: 1. Dfc-North America, 2. Dfb-Mid America, 3. Dfa-Mid America, 4. Cfa-Southeast US, 5. Dfb-Rocky region, 6. Csb-Pacific, 7. Am-Amazon, 8. Aw-Cerrado, 9. Cfa-Coastal, 10. Dfc-North Europe, 11. Dfb-North Eurasia, 12. Cfb-England, 13. Cfb-West Europe, 14. Dfc-Alp, 15. Cwb-South Africa, 16. Aw-India, 17. BWk-Tarim, 18. Cfa-Japan, 19. BSh-Northeast Australia) (Am/Aw is equatorial monsoon/winter dry climate, BWk/BSh is arid desert cold/hot climate, Cfa/Cfb and Csb/Cwb is warm temperate fully humid hot/warm summer and summer/winter dry warm summer climate, Dfa/Dfb/Dfc is snow fully humid hot/warm/cool summer climate).</figcaption></figure>

we quasi-quantitatively analyze the occurrence time of three potential influential factors, i.e. the extreme precipitation, soil moisture excess and snowmelt and infer the main causes leading to the flood timing shifts.

## Results

**Spatial distribution of mean flooding time and its concentration index.** There is a distinct regional pattern in the spatial distribution of the mean date of occurrence and concentration index of global floods pattern (Fig. 1). In the Northern Hemisphere, there is a clear transition of flood timing in the latitudinal direction except for some regions suffering from regional topographical or monsoon effects (Fig. 1b). Winter flooding typically takes place across the western coast of America, large parts of Western Europe and the Mediterranean (with high concentration), while spring flooding generally takes place in eastern and southeastern America and western continental Europe. It is worth noting that the concentration index of floods in southeast America is quite low, meaning that flood timing in this region is widely dispersed throughout the year due to a combination of tropical and extratropical systems controlling the flooding<sup>25</sup>. Summer flooding is found around north America, the Alps, India and Japan. These patterns have been described in detail in regional studies<sup>3,20,22,41</sup>. In southern Hemisphere, the available flood timing data were mainly available in South Africa, Australia and Brazil. Flooding normally occurs in December, January and February in south Africa, north Australian and central Brazil (with concentration index approaching 1.0 in these regions in Fig. 1c). An interesting phenomenon is observed in Brazil: a contrasting flood behavior is observed in the adjacent Cerrado and Coastal regions, as flood
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Global distribution of the trend in observed flood timing during 1970–2010 with latitude average flooding date (gray line; for display purposes the magnitude of the flood date is scaled by a factor of 0.01) and trend in flood timing (blue line) (a). Boxplots b–f demonstrate the relations between the trend in flood timing and mean flood date in the northern hemisphere (NH) (b) and southern hemisphere (SH) (c), latitude (d), outlet elevation (e) and catchment area (f). The interquartile range, equaling the difference between the 25th and 75th percentiles, is used to characterize variability for each bin. The red lines represent a loess-curve fitted to the values. The station density was also displayed in subplot (b–f). Subplot (g) shows the multiple linear regression analysis between flood timing trend with flood mean date, latitude, elevation and area.</figcaption></figure>

g. Multiple linear regression

Coefficients:

<table><thead><tr><th></th><th>Estimate</th><th>P</th></tr></thead><tbody><tr><td>(Intercept)</td><td>2.946e-1</td><td>0.1773</td></tr><tr><td>Mean flood date</td><td>5.454e-3</td><td>4.02e-06***</td></tr><tr><td>Latitude</td><td>-3.471e-2</td><td>&lt;2e-16***</td></tr><tr><td>Elevation</td><td>-3.981e-4</td><td>0.0968</td></tr><tr><td>Catchment area</td><td>6.294e-7</td><td>0.5042</td></tr></tbody></table>

in the coastal regions was jointly influenced by cold fronts, thunderstorms, and tropical cyclones, which make rainfall-induced floods occurring throughout the year, which leads to a high variability in flood timing with a low concentration (R = 0.31)²⁵,³¹,⁴².

Though the global flood timing exhibits strong seasonality at most stations (i.e., 4590 stations, out of a total of 6167 stations, have a concentration index larger than 0.5 and 2764 stations 0.7), low concentration (i.e., without a specific flooding season) was observed in several hotspots, for example, Hotspot 3 and 4 (mid and southeast United States), Hotspot 9 (Costal of South America), and certain hotspots in Europe (11. Dfb-North Eurasia, 12. Cfb-England and 13. Cfb-West Europe) due to the combined system as mentioned above²²,²³ (Fig. 1c). Stations at high latitude north hemisphere and equatorial areas showed high concentration while stations at high latitude south hemisphere demonstrated low concentration (Fig. 1d). We excluded these hotspots with concentration indexes less than 0.7 when studying the influential factors in flood timing shifts.
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Long-term temporal evolutions of flood timings and their influential factors for 13 hotspots with concentrate index > 0.7. Solid lines show the median timing over the entire hotspot: Green: timing of the observed floods; purple: 7-day maximum precipitation; orange: snow melting indicator and blue: timing of the calculated maximum soil moisture excess (SME). All data were subjected to a 10-year weighted moving average filter. The shaded bands indicate the timing variability within the year (±0.5 standard deviations). Vertical axes show month of the year (note different starting months in these panels). The precipitation and SME data were started from 1979 in south Africa as the CPC precipitation data was started from 1979.</figcaption></figure>

**Shift in flood timing.** Figure 2A shows the spatial patterns in the changes of global flood timing during 1970–2010. The trend in flood timing ranges from – 22 days per decade towards earlier floods to 28 days per decade towards delayed floods using the adjusted Theil-Sens slope. Floods tend to occur earlier in north and mid-America, the Mexican Gulf region, Great Lake, north of Australia and large parts of Europe. In contrary, at the Amazon, southern Brazil, India, Japan and northeast Australia, flood timing tends to be delayed (Table S1).

To some extent, the trend in flood timing is significantly related to mean date of floods and catchment latitude (Fig. 2b, c, d, g). Generally, Winter/spring flooding (approx. day 330~61 in the northern hemisphere) tends to occur earlier (Fig. 2b, c) (e.g. Hotspots 1, 2, 5, 6, 10, 11, 12, 13), while a summer or early autumn flooding (e.g. Spots 16, 18) is normally featured with delayed flood peaks. The higher the latitude, the earlier the flood timing in the north hemisphere (Fig. 2d). Catchment elevation and area also have an impact on flood timing but their effects were not statistically significant (Fig. 2e, f, g).

**Influencing factors on the shifting of flood timing.** As floods are the results of the interplay between the precipitation, soil moisture excess and snow processes³,⁴⁴, below we analyzed the shifting in these variables and their relations with flood timings. For the north and mid-America and north Europe (Hotspot 1, 2, 10), the spring snowmelt generally occurs during spring or early summer (Fig. 3). In these regions, the temporal changes of the timing of snowmelt maxima matches well with the timing of the floods, which lines up with previous regional studies²²,⁴⁵. At the west Pacific coast of America (Hotspot 6), Amazon (Hotspot 7), Cerrado Savanna (Hotspot 8), South Africa (Hotspot 15), India (Hotspot 16) and Northeast Australia (Hotspot 19), the precipitation and soil moisture excess could occur during the same period in a year (within 10 days lag), and the extreme precipitation and soil moisture excess (normally driven by a several-day heavy precipitation) jointly contributed to the floods, which could be inferred from the similar changing pattern and contemporaneous occurrence times of the extreme precipitation, soil moisture excess and floods (Fig. 3). Coincidentally, these hotspots with delayed flood timing are located in the monsoon dominated hotspots, including south Africa (Hotspot 15), India (Hotspot 16), Japan (Hotspot 18) and the Amazon areas (Hotspot 7).

**Discussion**

This study explored the spatial pattern of global flood timing, its shifting and reasons. Compared to those studies on the mean flood date and flood-timing shifting in Europe and north America, this study lines up with the existing literature³,²³,²⁴,⁴⁶.
On a global or regional scale, one could analyze the causes behind the observed pattern based on the overlapped time window of floods and other factors. For example, the earlier floods in Scandinavia (in Hotspot 10) are likely to be predominately influenced by snowmelt, of which the timing of the maximum 7-day value is the same with the flood timing (Fig. 3). This method has proven to be effective and relatively simple and practical compared to hydrological modelling or earth-system simulations in interpreting the flood mechanisms³,²⁰.

The flood timing is expected to occur later in the monsoon dominated regions despite few sites showing statistically significant trends (Figure S1a). For South Africa, the spatial pattern of the flood timing changes (based on limited stations) is similar with those based on the reanalysis data³². For India and Japan, the flood is more related to the monsoon indicators and delayed flood peaks have been observed, of which may be related to increased atmospheric effective specific capacity and slower response to seasonal solar forcing⁴⁷. For the south of Brazil (HS 8, Cerrado), the flood timing is getting delayed, which is consistent with the analyzed outcomes based on the local water-year flood¹⁰. In future studies in the monsoon regions, more focus should be on the relation between the floods and large-scale circulation conditions as floods may be more dependent on these circulations³².

The flood timing occurred earlier in the snow melting dominated regions (Figure S1b), such as the Great Lake region (HS 2), Scandinavian region (HS 10), East Germany (HS 13), South English Channel and the United Kingdom region (HS12). The continuous global warming is considered a possible mechanism pushing forward the snow melting floods⁴⁶,⁴⁸. If the snow processes become less important in the future, the flood seasonality patterns in these regions may alter completely as other flood generation processes become dominant, e.g., more rain on snow events²².

The changes in soil moisture excess have also resulted in earlier floods in some tropical or warm temperate areas, such as North Australia (HS 19). Soil moisture excess should be a more realistic factor to predict floods compared to the precipitation in several regions (e.g., South Africa[HS15], northeast Australia [HS19]), as the timing of soil moisture excess coincide better with flood timing.

Despite above promising result, there is still some space for future improvement. For example, only the timing of the annual 7-day maximum flood was studied, which excludes the fact that there might be more than one flood peak for some stations (e.g. two flood seasons). The regional interpretation of flood timing shifts is influenced by the limited and non-uniform spatial coverage, inconsistency of datasets and number of hydrometric stations, which causes high uncertainty, especially for the data-scarce region. Additionally, the flood timing in each calendar year rather than water year is recorded in GISM, which also increased the uncertainty in certain stations with a sequence of dependent flood events being counted in two successive years. Nested catchments were not identified. Furthermore, the combined effect from rainfall and soil moisture excess has not been distinguished yet.

Flood timing shifting will certainly have implications on water management practice and food security. For regions with a strong signal of flood timing changes, local practices of flood management should be adjusted to adapt to the shifted flood peaks, such as preparing for the earlier flood period, discharging reservoirs for later flood⁴⁹. Another practical influence of shifted flood timing is closely related to food security. The water shortage may be prevailing during the water demanding season due to earlier or later flood peaks, resulting in reduction of crop yield and becoming a threat to food security.

## Conclusions

This work provides a study on global flood timing shifts and contributions from three main influential factors (i.e., maximum 7-day precipitation, soil moisture excess, and snowmelt). The shifts in flood timing are calculated with circular index based on the observed flood data from 1970 to 2010. The causality relationship of flood timing and its three influential factors were investigated based on visual interpretation.

Despite the difference in the spatial distribution of the flood timing shifts, there are some spatial patterns: Earlier floods are more prevailing in high-latitude regions, mainly caused by earlier snowmelt due to global warming. Regions with delayed floods were mainly distributed in the monsoon-affected areas (e.g., South Africa, India and Japan), where the delayed floods were mainly the interplay of delayed rainstorm and soil moisture excess.

Future climate change will have a big impact on the global precipitation and snowmelt. This will eventually lead to the shifting on global flood time. This study provides an insight on the shifting mechanism of flood timing, and could be used to support decisions on the global flood mitigation and the impact from future climate change.

## Data and methods

**Data.** In this study, we analyzed a large data set of flood observations to study the shifting in flood timing at 6167 gauging stations from 1970 to 2010. The Global Streamflow Indices and Metadata (GSIM) archive and the Global Runoff Data Centre (GRDC) dataset were jointly used.

For climatic variables, daily precipitation data from the Global Historical Climatology Network (GHCN), the Asian Precipitation Highly Resolved Observational Data Integration Towards Evaluation of Water Resources (APHRODITE), the CPC Global Unified Precipitation data provided by the NOAA/OAR/ESRL PSL, South America Daily Gridded Precipitation (1979~) and the NCEP/NCAR Reanalysis Products 1 were also used. The CRU TS4.05 maximum and minimum temperature dataset was interpolated to daily scale using spline interpolation method in R Package "stats". The daily temperature data were used to calculate the potential evapotranspiration and snowmelt.

We restricted our analysis to sites with no less than 37 years of data (flood) during period 1970–2010 (allowing for a maximum of 4 years’ missing data), resulting in a total of 6167 sites across the world. The study period is selected based on data availability. For Africa and Asia, few gauges were used due to data availability. Note that for the south America, the CPC precipitation only records daily precipitation during 1979 to recent. More details on the datasets and preprocessing are described in Text S1 in Supporting Information.
**Definition of flood and its influential factors.** In this study, floods are represented by the annual maximum 7-day flow. Regarding the flood timing influential factors, we assume floods were only caused by the following climatic factors, i.e. the precipitation, soil moisture excess and snow melt, since most floods are caused by strong rainfall, strong snowmelt or rain on saturated soil$^{3,19,20,44}$. Only natural process factors were being analyzed, omitting human factors such as land use change and engineering of rivers, which are however indirectly reflected in the gauged data.

*Extreme precipitation* in this case, we assume that floods were resulted from the largest 7-day precipitation. Moreover, the occurring times of flood and extreme precipitation should fall into the same period of the year.

*Soil moisture excess* in this case, floods are triggered by high antecedent soil moisture conditions and a strong rainfall during long periods (up to months) rather than short duration rainfall events. The soil moisture excess combines daily precipitation, evapotranspiration, and the antecedent soil moisture and is defined as the daily precipitation amount minus the available soil moisture storage capacity$^{16,19,48}$. Soil moisture excess can also be referred to as surface runoff based on the calculation procedures (Eqs. 1, 2).

$$Su(t) = \max[\min(Su(t-1) + P(t), Su_{max}) - E(t), 0] \quad (1)$$

$$E(t) = \min[0.75 \times ET0, Su(t-1)] \quad (2)$$

where $Su(t)$ and $Su(t-1)$ represent the soil moisture storages on day $t$ and day $t-1$. $Su_{max}$ is the maximum soil moisture storage fixed at 25 mm. Changing this to 75 mm led to the soil moisture excess being 0 for most days as the precipitation will always smaller than evaporation. The changes of the maximum soil moisture did not substantially affect the global results$^{19}$. $P(t)$ and $E(t)$ demonstrate the precipitation and evapotranspiration amount on day $t$, respectively. $ET0$ is the reference evapotranspiration, calculated by the Hargreaves-Samani method in the R evapotranspiration package$^{50,51}$. $E$ is scaled to 75% of its daily value because not all $ET0$ were utilized for evapotranspiration.

*Snowmelt* In snowmelt induced floods, the maximum annual flow is assumed to be caused by the largest 7-day snowmelt. The latter is calculated based on an empirical degree-day model$^{51}$.

$$Ss(t) = \begin{cases} \min(fdd * (TT(t) - t_{crit}), snowpack(t-1)) & TT(t) > t_{crit} \\ 0 & TT(t) \le t_{crit} \end{cases} \quad (3)$$

$$snowpack(t) = \begin{cases} snowpack(t-1) - Ss(t), & TT(t) > t_{crit} \\ snowpack(t-1) + P(t), & TT(t) \le t_{crit} \end{cases} \quad (4)$$

where $Ss(t)$ and $snowpack(t)$ are the snowmelt and snow storage on day $t$ (mm), respectively. $TT(t)$ is the daily maximum temperature on day $t$, $t_{crit}$ the temperature threshold for snowmelt and for the delimiting rainfall and snowfall, which is set at 5 °C. $fdd$ is the snow melting factor, established at 2 mm per day per °C$^{19}$.

**Seasonality characteristics of the flood date and its influencing factors.** In this study, we calculated the mean flooding date and the flood concentration at each station. Then the trend in the flooding timing is estimated using the adjusted Theil-Sen slope method.

The seasonality of flooding and its influencing factors is characterized by circular statistics$^{3,19,20,24}$. When applying this analysis, we firstly excluded the flow stations whose floods are uniformly distributed within the year (through the Kuiper's test with a significance level 0.05) and 5763 stations were used in the following analysis. To calculate the mean flood date and its concentration statistics, the flood occurrence date, $D_i$ (day of the year), was converted into an angular value $\theta_i$. For station $i$,

$$\theta_i = D_i \cdot \frac{2\pi}{m_i} \quad (5)$$

where $D_i=1$ corresponds to January 1st and $D_i=m_i$ to December 31st, and $m_i$ is the Julian days in year $i$ (i.e. 365 in regular years and 366 in leap years). The average date of occurrence $\bar{D}$ of flood (at a station) is defined as:

$$\bar{D} = \begin{cases} \tan^{-1}\left(\frac{\bar{y}}{\bar{x}}\right) \cdot \frac{\bar{m}}{2\pi}, & \bar{x} > 0, \bar{y} > 0 \\ \left[\tan^{-1}\left(\frac{\bar{y}}{\bar{x}}\right) + \pi\right] \cdot \frac{\bar{m}}{2\pi}, & \bar{x} \le 0 \\ \left[\tan^{-1}\left(\frac{\bar{y}}{\bar{x}}\right) + 2\pi\right] \cdot \frac{\bar{m}}{2\pi}, & \bar{x} > 0, \bar{y} < 0 \end{cases} \quad (6)$$

where $\bar{x} = \frac{1}{n} \sum_{i=1}^{n} \cos(\theta_i)$ and $\bar{y} = \frac{1}{n} \sum_{i=1}^{n} \sin(\theta_i)$ ($\theta_i$ is defined in Eq. 5) are the cosine and sine components of the average date, respectively, $\bar{m} = \frac{1}{n} \sum_{i=1}^{n} m_i$ represents the average number of days per year and $n$ the total number of floods at each station$^3$.

To illustrate the variability of the flood occurrence date, a concentration statistic R of the occurrence date around the average date ($\bar{D}$) is defined:

$$R = \sqrt{\bar{x}^2 + \bar{y}^2}, \quad 0 \le R \le 1 \quad (7)$$

where $R=0$ indicates that flood occurrence dates are widely dispersed or sitting opposite on the unit circle throughout the year, and $R=1$ indicates that all flooding events are concentrated on the same ordinal day.
To quantify the trend in the timing of floods and meteorological divers, the adjusted Theil-Sen slope estimator was applied<sup>52</sup>. This non-parametric estimator β was chosen for its robustness and insensitivity to missing values and outliers. The trend estimator β was estimated as:

$$ \beta = \text{median} \left( \frac{D_j - D_i + k}{j - i} \right) \qquad (8) $$

where $k = \begin{cases} -\bar{m}, & D_j - D_i > -\frac{\bar{m}}{2}, \\ \bar{m}, & D_j - D_i < -\frac{\bar{m}}{2}, \end{cases}$ is the adjustment factor and β denotes the slope of the flood timing changes<sup>3</sup>.

**Dominant climatic factors in flood timing changes.** In this study, 19 hotspots are identified according to the Köppen-Geiger classification<sup>43</sup> and the similarity of mean flood dates (Fig. 1). Names of these hotspots consist of the climate zone and geographical location. Note that the names are only indicative for a region and do not exactly correspond to any exactly defined geographic area.

Since the available flood data are limited to the flood dates of each year instead of daily streamflow, we identified the dominant flood-influencing factor by comparing the flood dates with the timing of each candidate mechanism<sup>19</sup>. Our assumption is that two variables are causally related if they are in the same occurrence time window (i.e., occur simultaneously) and have similar trends.

## Data availability

The flood date data used in this study can be downloaded from https://doi.pangaea.de/10.1594/PANGAEA.887477 and https://www.bafg.de/GRDC/EN/Home/homepage_node.html. The precipitation can be downloaded from https://www.ncei.noaa.gov/products/land-based-station/global-historical-climatology-network-daily for America, Europe and Australia, from http://aphrodite.st.hirosaki-u.ac.jp/products.html for Asia region, from https://psl.noaa.gov/data/gridded/data.cpc.globalprecip.html and https://psl.noaa.gov/data/gridded/data.ncep.reanalysis.html for south America. The temperature data can be downloaded from https://cran.dta.uea.ac.uk/cru/data/hrg/cru_ts_4.05/.

Received: 18 June 2022; Accepted: 4 November 2022

Published online: 07 November 2022

## References

1. Hirabayashi, Y. et al. Global flood risk under climate change. Nat. Clim. Change 3, 816–821. https://doi.org/10.1038/nclimate1911 (2013).

2. Dottori, F. et al. Increased human and economic losses from river flooding with anthropogenic warming. Nat. Clim. Change 8, 781–786. https://doi.org/10.1038/s41558-018-0257-z (2018).

3. Blöschl, G. et al. Changing climate shifts timing of European floods. Science 357, 588–590 (2017).

4. Winsemius, H. C. et al. Global drivers of future river flood risk. Nat. Clim. Change 6, 381–385 (2016).

5. Corringham, T. W., Ralph, F. M., Gershunov, A., Cayan, D. R. & Talbot, C. A. Atmospheric rivers drive flood damages in the western United States. Sci. Adv. 5, eaax4631 (2019).

6. Boulange, J., Hanasaki, N., Yamazaki, D. & Pokhrel, Y. Role of dams in reducing global flood exposure under climate change. Nat. Commun. 12, 417. https://doi.org/10.1038/s41467-020-20704-0 (2021).

7. Mallakpour, I. & Villarini, G. The changing nature of flooding across the central United States. Nat. Clim. Change 5, 250–254. https://doi.org/10.1038/nclimate2516 (2015).

8. Yin, J. et al. Large increase in global storm runoff extremes driven by climate and anthropogenic changes. Nat. Commun. 9, 1–10 (2018).

9. Wasko, C., Sharma, A. & Lettenmaier, D. P. Increases in temperature do not translate to increased flooding. Nat. Commun. 10, 1–3 (2019).

10. Wasko, C., Nathan, R. & Peel, M. C. Trends in global flood and streamflow timing based on local water year. Water Resour. Res. 56, e2020WR027233 (2020).

11. Do, H. X. et al. Historical and future changes in global flood magnitude—Evidence from a model-observation investigation. Hydrol. Earth Syst. Sci. 24, 1543–1564. https://doi.org/10.5194/hess-24-1543-2020 (2020).

12. Barichivich, J. et al. Recent intensification of Amazon flooding extremes driven by strengthened Walker circulation. Sci. Adv. 4, eaat8785 (2018).

13. IPCC. Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge University Press, 2021).

14. Blöschl, G. et al. Changing climate both increases and decreases European river floods. Nature 573, 108–111 (2019).

15. Hall, J. et al. Understanding flood regime changes in Europe: A state-of-the-art assessment. Hydrol. Earth Syst. Sci. 18, 2735–2772 (2014).

16. Bertola, M. et al. Do small and large floods have the same drivers of change? A regional attribution analysis in Europe. Hydrol. Earth Syst. Sci. 25, 1347–1364. https://doi.org/10.5194/hess-25-1347-2021 (2021).

17. Do, H. X., Mei, Y. & Gronewold, A. D. To what extent are changes in flood magnitude related to changes in precipitation extremes?. Geophys. Res. Lett. 47, e2020GL088684. https://doi.org/10.1029/2020gl088684 (2020).

18. Hodgkins, G., Dudley, R., Archfield, S. A. & Renard, B. Effects of climate, regulation, and urbanization on historical flood trends in the United States. J. Hydrol. 573, 697–709 (2019).

19. Berghuijs, W. R., Harrigan, S., Molnar, P., Slater, L. J. & Kirchner, J. W. The relative importance of different flood-generating mechanisms across Europe. Water Resour. Res. 55, 4582–4593. https://doi.org/10.1029/2019WR024841 (2019).

20. Berghuijs, W., Woods, R., Hutton, C. & Sivapalan, M. Dominant flood generating mechanisms across the United States. Geophys. Res. Lett. 43, 4382–4390. https://doi.org/10.1002/2016GL068070 (2016).

21. Dettinger, M. D. & Diaz, H. F. Global characteristics of stream flow seasonality and variability. J. Hydrometeorol. 1, 289–310. https://doi.org/10.1175/1525-7541(2000)001%3c0289:gcosfs%3e2.0.co;2 (2000).

22. Hall, J. & Blöschl, G. Spatial patterns and characteristics of flood seasonality in Europe. Hydrol. Earth Syst. Sci. 22, 3883–3901. https://doi.org/10.5194/hess-22-3883-2018 (2018).

23. Villarini, G. On the seasonality of flooding across the continental United States. Adv. Water Resour. 87, 80–91. https://doi.org/10.1016/j.adwatres.2015.11.009 (2016).
24. Ye, S. et al. Understanding flood seasonality and its temporal shifts within the contiguous United States. *J. Hydrometeorol.* **18**, 1997–2009. https://doi.org/10.1175/jhm-d-16-0207.1 (2017).

25. Do, H. X., Westra, S., Leonard, M. & Gudmundsson, L. Global-scale prediction of flood timing using atmospheric reanalysis. *Water Resour. Res.* https://doi.org/10.1029/2019wr024945 (2020).

26. Mediero, L., Santillán, D., Garrote, L. & Granados, A. Detection and attribution of trends in magnitude, frequency and timing of floods in Spain. *J. Hydrol.* **517**, 1072–1088. https://doi.org/10.1016/j.jhydrol.2014.06.040 (2014).

27. Burn, D. H. & Whitfield, P. H. Changes in floods and flood regimes in Canada. *Can. Water Resour. J.* **41**, 139–150. https://doi.org/10.1080/07011784.2015.1026844 (2016).

28. Huang, H. et al. Changes in mechanisms and characteristics of western U.S. floods over the last sixty years. *Geophys. Res. Lett.* **49**, e2021GL097022. https://doi.org/10.1029/2021GL097022 (2022).

29. Collins, M. J., Hodgkins, G. A., Archfield, S. A. & Hirsch, R. M. The occurrence of large floods in the United States in the modern hydroclimate regime: Seasonality, trends, and large-scale climate associations. *Water Resour. Res.* https://doi.org/10.1029/2021W R030480 (2022).

30. Wasko, C., Nathan, R. & Peel, M. C. Changes in antecedent soil moisture modulate flood seasonality in a changing climate. *Water Resour. Res.* https://doi.org/10.1029/2019wr026300 (2020).

31. Cassalho, F., Beskow, S., de Mello, C. R., Oliveira, L. F. & de Aguiar, M. S. Evaluation of flood timing and regularity over hydrological regionalization in southern Brazil. *J. Hydrol. Eng.* https://doi.org/10.1061/(asce)he.1943-5584.0001815 (2019).

32. Ficchi, A. & Stephens, L. Climate variability alters flood timing across Africa. *Geophys. Res. Lett.* **46**, 8809–8819. https://doi.org/10.1029/2019gl081988 (2019).

33. Mao, Y. et al. Flood inundation generation mechanisms and their changes in 1953–2004 in global major river basins. *J. Geophys. Res. Atmos.* **124**, 11672–11692 (2019).

34. Bertola, M., Viglione, A. & Blöschl, G. Informed attribution of flood changes to decadal variation of atmospheric, catchment and river drivers in Upper Austria. *J. Hydrol.* **577**, 123919. https://doi.org/10.1016/j.jhydrol.2019.123919 (2019).

35. Schroter, K., Kunz, M., Elmer, F., Mühr, B. & Merz, B. What made the June 2013 flood in Germany an exceptional event? A hydro-meteorological evaluation. *Hydrol. Earth Syst. Sci.* **19**, 309–327 (2015).

36. Nie, Y. & Sun, J. Moisture sources and transport for extreme precipitation over Henan in July 2021. *Geophys. Res. Lett.* **49**, e2021GL097446. https://doi.org/10.1029/2021GL097446 (2022).

37. Zion, M. S. et al. Investigation and modeling of winter streamflow timing and magnitude under changing climate conditions for the Catskill Mountain region, New York, USA. *Hydrol. Process.* **25**, 3289–3301 (2011).

38. Zhang, Q. et al. Magnitude, frequency and timing of floods in the Tarim River basin, China: Changes, causes and implications. *Glob. Planet Change* **139**, 44–55. https://doi.org/10.1016/j.gloplacha.2015.10.005 (2016).

39. Pall, P. et al. Anthropogenic greenhouse gas contribution to flood risk in England and Wales in autumn 2000. *Nature* **470**, 382–385. https://doi.org/10.1038/nature09762 (2011).

40. Edwards, T. L. et al. Projected land ice contributions to twenty-first-century sea level rise. *Nature* **593**, 74–82. https://doi.org/10.1038/s41586-021-03302-y (2021).

41. Bartiko, D., Oliveira, D. Y., Bonumá, N. B. & Chaffe, P. L. B. Spatial and seasonal patterns of flood change across Brazil. *Hydrol. Sci. J.* **64**, 1071–1079. https://doi.org/10.1080/02626667.2019.1619081 (2019).

42. Teixeira, M. D. S. & Satyamurty, P. Trends in the frequency of intense precipitation events in southern and southeastern Brazil during 1960–2004. *J. Clim.* **24**, 1913–1921 (2011).

43. Beck, H. E. et al. Present and future Köppen-Geiger climate classification maps at 1-km resolution. *Sci. Data* **5**, 180214. https://doi.org/10.1038/sdata.2018.214 (2018).

44. Sivapalan, M., Bloschl, G., Merz, R. & Gutknecht, D. Linking flood frequency to long-term water balance: Incorporating effects of seasonality. *Water Resour. Res.* https://doi.org/10.1029/2004wr003439 (2005).

45. Arheimer, B. & Lindström, G. Climate impact on floods: Changes in high flows in Sweden in the past and the future (1911–2100). *Hydrol. Earth Syst. Sci.* **19**, 771–784. https://doi.org/10.5194/hess-19-771-2015 (2015).

46. Singh, J., Ghosh, S., Simonovic, S. P. & Karmakar, S. Identification of flood seasonality and drivers across Canada. *Hydrol. Process.* https://doi.org/10.1002/hyp.14398 (2021).

47. Song, F. et al. Emergence of seasonal delay of tropical rainfall during 1979–2019. *Nat. Clim. Change* **11**, 605–612 (2021).

48. Stein, L., Pianosi, F. & Woods, R. Event-based classification for global study of river flood generating processes. *Hydrol. Process.* **34**, 1514–1529. https://doi.org/10.1002/hyp.13678 (2020).

49. Yang, Y., Yang, L., Chen, X., Wang, Q. & Tian, F. Climate leads to reversed latitudinal changes in chinese flood peak timing. *Earths Future* **10**, e2022EF002726 (2022).

50. Guo, D., Westra, S. & Maier, H. R. An R package for modelling actual, potential and reference evapotranspiration. *Environ. Model. Softw.* **78**, 216–224. https://doi.org/10.1016/j.envsoft.2015.12.019 (2016).

51. Hargreaves, G. H. & Samani, Z. A. Reference crop evapotranspiration from temperature. *Appl. Eng. Agric.* **1**, 96–99 (1985).

52. Sen, P. K. Estimates of the regression coefficient based on Kendall's tau. *J. Am. Stat. Assoc.* **63**, 1379–1389 (1968).

## Acknowledgements

This research was supported by the National Natural Science Foundation of China (42071046, 42130512), the Open Project of Key Laboratory of Xinjiang, China (2022D04043) and the Western Young Scholar of the Chinese Academy of Sciences. The authors gratefully acknowledge the Youth Innovation Promotion Association of the Chinese Academy of Sciences (No. 2019431). No conflict of interest exits in the submission of this manuscript.

## Author contributions

G.F. and Z.L. conceived the idea, performed the analysis and wrote the draft. J.Y. and Y.C. designed the research.

C.A and W.D. contributed to results analysis and paper writing. P.D. provided comments and revised the manuscript. All authors approved the content of the manuscript and contributed to the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-23748-y.

Correspondence and requests for materials should be addressed to Y.C.

Reprints and permissions information is available at www.nature.com/reprints.
**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**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/.

© The Author(s) 2022