# Will rivers become more intermittent in France? Learning from an extended set of hydrological projections

Tristan Jaouen, Lionel Benoit, Louis Héraut, Eric Sauquet

► To cite this version:

Tristan Jaouen, Lionel Benoit, Louis Héraut, Eric Sauquet. Will rivers become more intermittent in France? Learning from an extended set of hydrological projections. Hydrology and Earth System Sciences, 2025, 29 (15), pp.3629-3671. 10.5194/hess-29-3629-2025. hal-05326845

HAL Id: hal-05326845

https://hal.inrae.fr/hal-05326845v1

Submitted on 22 Oct 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_1.png" /></figure>

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

# Will rivers become more intermittent in France? Learning from an extended set of hydrological projections

Tristan Jaouen¹, Lionel Benoit², Louis Héraut¹, and Eric Sauquet¹

¹INRAE Lyon-Grenoble Auvergne-Rhône-Alpes, RiverLy, Villeurbanne, France

²INRAE PACA, Biostatistique et processus SPatiaux, Avignon, France

**Correspondence:** Tristan Jaouen (trjaouen@laposte.net) and Eric Sauquet (eric.sauquet@inrae.fr)

Received: 31 August 2024 – Discussion started: 19 September 2024
Revised: 4 March 2025 – Accepted: 16 April 2025 – Published: 8 August 2025

**Abstract.** This study aims to assess the changes in the intermittence of river flows across France in the context of climate change. Projections of flow intermittence are derived from the results of the Explore2 project, which is the latest national study that proposes a wide range of potential hydrological futures for the 21st century. The multi-model approach developed within the Explore2 project enables uncertainties in future flow intermittence to be characterized. Combined with discrete observations of flow states, hydrological projections are post-processed to compute the daily probability of flow intermittence (PFI) on each element of the partition of France in hydro-ecoregions (HERs).

The post-processing consists of calibrating logistic regressions between the historical flow states of the National Low-Flow Observatory (ONDE) network and the flow data simulated by the hydrological models involved in Explore2 run with the SAFRAN atmospheric reanalysis as inputs. After calibration, these regressions are used to project daily PFIs for the entire 21st century, based on flow simulations from five hydrological models driven by up to 17 climate projections under RCP2.6, RCP4.5, and RCP8.5 climate change scenarios.

The results show good agreement among the hydrological models regarding the increase in flow intermittence under RCP4.5 and RCP8.5. The projected increase in mean daily PFI between July and October and the shift of the first and last days when PFI exceeds 20% both suggest a gradual intensification and extension of dry spells throughout the century. The southern regions of France are likely to experience greater increases in runoff intermittence than the northern regions, and mountainous regions such as the Alps and the Pyrenees are likely to experience changes in their dynamics

of intermittence with a reduction in winter intermittence and the apparition of or increase in summer intermittence. The uncertainty of these projected changes is larger in northern France due to greater intermodel variability in this region.

## 1 Introduction

Intermittent rivers and ephemeral streams (IRESs) are watercourses characterized by the absence of continuous year-round water flow (Sefton et al., 2019). While water management practices can exacerbate flow intermittence by altering runoff patterns, climatic factors – particularly aridity – remain the primary determinants of these systems (Addor et al., 2018; Hammond and Fleming, 2021; Sauquet et al., 2021; Zipper et al., 2021). In recent decades, climate change has intensified these dynamics, significantly impacting the global water cycle: summer precipitations tend to decrease, and rising temperatures have increased evapotranspiration and driven more frequent heatwaves (Douville et al., 2021). These changes have led to more severe low-flow periods and more frequent and widespread dry spells in Western Europe (Delso et al., 2017; Tramblay et al., 2020; Vicente-Serrano et al., 2020). Meanwhile, IRESs play a critical role as interfaces between terrestrial and aquatic ecosystems. Altered patterns of flow interruption can therefore compromise water quantity and quality in downstream perennial rivers, increase extinction risks for specialized species, and disrupt organic matter, nutrient, and sediment cycles (Bertassello et al., 2022; Galhart et al., 2012; Giezendanner et al., 2021; Finn et al., 2011; Meyer et al., 2007; Sarremejane et al., 2017; Larned et al., 2010). Hence, a better understanding of IRES behaviours and
their interaction with perennial streams is needed to assess the impacts of ongoing and future changes on river ecosystems (Döll and Schmied, 2012; Jaeger et al., 2014; Pumo et al., 2016; Leigh and Datry, 2017) and to address water management challenges (Acuña et al., 2014).

Globally, IRESs account for nearly two-thirds of river networks (Schneider et al., 2017; Messager et al., 2021), with this proportion reaching 39% in France (Snelder et al., 2013). Their spatial and temporal dynamics can be characterized using conceptual models of network contraction (Shaw et al., 2017; Botter and Durighetto, 2020) or through semi-conceptual and machine learning frameworks that integrate hydroclimatic and physiographic factors (Lapides et al., 2021; Durighetto et al., 2022). The strongest predictors of flow intermittence are aridity indices (Vicente-Serrano et al., 2019; Jaeger et al., 2019). However, local variations in permeability, lithology, and topography also influence intermittence patterns in specific regions, including West Africa (Yu and Duffy, 2018; Belemtougri et al., 2021), northern Spain (González-Ferreras and Barquín, 2017), and north-western Australia (Bourke et al., 2021). Studying intermittence at a large scale therefore requires considering the diversity of hydrological processes at play, particularly when multiple watersheds exhibit diverse sensitivities to hydroclimatic conditions. This variability highlights the need for robust methodologies to assess intermittence at national to (semi-)continental scales (Addor et al., 2018; Hammond and Fleming, 2021; Sando et al., 2022; Döll et al., 2024).

Understanding the drivers of river intermittence is also crucial for analysing historical trends and projecting future changes in IRES behaviour. Climate change has already expanded the spatial and temporal extent of intermittent flows worldwide (Zipper et al., 2021; Jaeger et al., 2019; Sando et al., 2022; Gudmundsson et al., 2021), and the Mediterranean basin is particularly affected (Tramblay et al., 2020; De Grolamo et al., 2022). However, interpreting climate projections and emission scenarios at the scale of headwater basins remains rare (Schneider et al., 2013). Indeed, many climate models have a coarse spatial resolution, and high uncertainties emerge from multi-model ensembles (MMEs), which makes accurate predictions of future IRES changes difficult.

To address these challenges, France has implemented the Explore2 project (https://entrepot.recherche.data.gouv.fr/dataverse/explore2-projections_hydrologiques, Sauquet et al., 2025b) that uses a multi-model, multi-scenario approach to simulate hydrological conditions throughout for the 21st century. It assesses uncertainties at each step of the modelling process. However, the resolution of the climate projections involved in Explore2 is too coarse to accurately capture IRES dynamics. Hence, this study extends the results of Explore2 for purposes of generating flow intermittence projections for French headwater streams. To this end, we follow a statistical approach that estimates the daily probability of flow intermittence (PFI) in small streams based on streamflow data simulated by Explore2 for larger perennial rivers. The estimated

PFI serves as a proxy indicator for the intensity and duration of dry periods. The method is calibrated using field observations carried out regularly at more than 3200 upstream river sites prone to drying.

To achieve this, the present work builds on two key studies: Beaufort et al. (2018), which demonstrated consistent performance of PFI estimation across the heterogeneous climate of France using daily discharge, and Sauquet et al. (2021), which projected PFI under a limited set of climate projections with coarse spatial resolution. In contrast, this study incorporates a finer spatial resolution, a broader range of climate scenarios, and a wider variety of hydrological models thanks to the integration of Explore2 projections. It therefore addresses a critical research gap by generating projections of small-stream intermittence across France and throughout the 21st century while quantifying the associated uncertainties.

The remainder of the paper is organized into five sections. The data used are outlined in the second section. Section 3 describes the statistical method linking PFI to hydrological projections. Section 4 presents the results, which focus on the mean daily PFI between July and October (mPFI<sub>7-10</sub>) and the median of the first and last days (respectively, Tf and Tl) when PFI exceeds 20%. Finally, Sect. 5 discusses these results, and Sect. 6 concludes the study.

## 2 Data

### 2.1 Monitoring river flow intermittence

The French Biodiversity Office (OFB; https://ofb.gouv.fr/, last access: 18 August 2024) initiated in 2012 the National Low-Flow Observatory (ONDE) to gain a better understanding of low flows and intermittent rivers. ONDE is a stream intermittence monitoring network comprising 3248 observation sites strategically distributed across the French river network with the aim of characterizing the occurrence and the intensity of summer (low-)flows. The network has been designed to focus on streams with a Strahler order ranging from 1 to 4, which are prone to natural and/or anthropogenic intermittent flow conditions. Most sites (85%) are located on "small streams" with a drainage area $\leq 100 \text{ km}^2$, 75% are situated on streams with a Strahler index of 1 or 2, and 20% of the sites have a drainage area $\leq 10 \text{ km}^2$ (Appendix Figs. G1 and G3).

A key feature of the ONDE network is that it involves the flow state of watercourses, which results from the visual inspection of streams by OFB observers at each location instead of a sensor-based measurement of flow rates. In the present study, two binary categories are defined for ONDE observations: (1) "visible flows" and (2) "dry conditions", which gather non-visible flows (standing water in isolated pools) and totally dry conditions (dry riverbed at or near the ONDE site). We use only data from regular and structured field campaigns (referred to as "usual campaigns"
in the ONDE nomenclature), during which the ONDE sites
are monitored systematically in mainland France around the
25th day of each month from May to September. Here we
use ONDE data from 2012 to 2022 for model calibration and
evaluation.

Based on ONDE data, the PFI is defined as the propor-
tion of small streams under drying conditions at the regional
scale. It is used as a proxy indicator for the intensity and
duration of dry periods. Long-term analysis of the PFI char-
acterizes the spatial and temporal trends of drought across
regions and provides valuable information about IRES inter-
mittent behaviour.

## 2.2 Delineating areas with homogeneous hydrological behaviour: the hydro-ecoregions

Hydro-ecoregions (HERs) are spatially homogeneous areas
defined over France based on natural drivers involved in river
ecosystem functioning, such as geology, topography, and cli-
mate. There are 85 level-2 HERs (HER2s) across France, and
they are derived from the sub-division of 22 level-1 HERs
(HER1s) (Wasson et al., 2002), which were used in a pre-
vious study about the impact of climate change on PFI in
France (Sauquet et al., 2021). In contrast, the present study
investigates flow intermittence at the scale of HER2s in order
to model PFI at a higher spatial resolution. The statistical ap-
proach developed hereafter requires sufficient observations
of flow states for calibration purposes (see Sect. 3.1), which
led us to perform five groupings of HER2 in order to ensure
that enough observations are available for model calibration.
When merging HER2s we checked that they belonged to the
same level-1 HERs to ensure that they share similar environ-
mental characteristics at the large scale. In the end, a parti-
tion of France into 75 entities (HER2 or pool of HER2s) was
considered for this study (Fig. 1).

These 75 HER2s have a median surface area of 4990 km²,
and more than 20% (17/75) have a surface area greater than
10 000 km². The HER2s were the subject of ONDE field
campaigns from May to September between 2012 and 2022,
thus resulting in 4125 campaigns. A total of 78 campaigns
are excluded from the analysis due to their failure to cover
more than 75% of the monitoring stations within a given
HER2. Thus, the data exclusion rate from the ONDE net-
work is less than 2%, leaving 4047 field campaigns with us-
able data.

## 2.3 Future hydrological projections over France: the
Explore2 project

The assessment of changes in the PFI over the 21st century is
based on the large set of hydrological projections produced
by the Explore2 project (Sauquet et al., 2025a). The daily
discharge projections were obtained at about 4000 simula-
tion points distributed along the French river network, with
a constraint on the minimum drainage area imposed by the

spatial resolution of climate projections (64 km²) and another
constraint of an even coverage across France.

Explore2 is a multi-model simulation experiment involv-
ing nine different hydrological models, but some of them are
restricted to regional applications (on average each simula-
tion point has discharge projections simulated by four hydro-
logical models). Finally, with the objective of predicting flow
intermittence at the national scale, we select the five mod-
els with the largest simulation domain: CTRIP (Decharme et
al., 2019), GRSD (de Lavenne et al., 2016), J2000 (Morel et
al., 2023), ORCHIDEE (Huang et al., 2024), and SMASH
(Jay-Allemand et al., 2020). The Explore2 project considers
only natural flows, although some of the hydrological mod-
els involved in the project have the capability of including
human-induced influences. The evaluation of hydrological
models was performed using an extended dataset of near-
pristine catchments (Strohmenger et al., 2023). The simu-
lations used for model evaluation were carried out for the
period 1976-2022 using the French near-surface SAFRAN
meteorological reanalysis (Vidal et al., 2010) as input. In ad-
dition, the hydrological models included in Explore2 have
been evaluated on the common period of availability 1976-
2019 (Sauquet et al., 2025a). After validation, the hydro-
logical models are forced by an ensemble dataset derived
from 17 pairs of global and regional climate models (GCM-
RCM) corrected by the statistical adjustment method ADA-
MONT (Verfaillie et al., 2018) for both historical (1976-
2004) and future climate scenarios RCP2.6, RCP4.5, and
RCP8.5 (2005-2100; Appendix A). In the end, these bias-
corrected climate projections were used as inputs for the hy-
drological models: all 17 GCM-RCM projections are used
under RCP8.5 scenario, whereas 10 of these projections are
used under RCP2.6 and nine under RCP4.5. The baseline pe-
riod is set to 1976-2005 and used thereafter for analysing
changes in IRES behaviour. Overall, the hydrological projec-
tions used thereafter to predict PFI result from combinations
of RCP-GCM-RCM-hydrological models.

# 3 Modelling framework

## 3.1 Linking intermittence to stream discharge at the HER2 scale

The empirical method suggested by Beaufort et al. (2018; Eq. 4) is applied here to estimate the probability of flow intermittence (PFI) at the scale of HER2s. This method follows a two-step process that links observed flow states at ONDE sites to daily streamflows (either gauged or modelled) using logistic regressions.

First, observations from ONDE sites are used to deter-
mine the percentage of sites with “dry conditions” for each
HER2 (*h*) and each campaign date (*j*). Flow intermittence
evolves over relatively long timescales, with seasonal tran-
sitions between flowing and dry states occurring gradually
<figure><img src="image_3.png" /><figcaption>Figure 1. Delineation of level-2 hydro-ecoregions (HER2s) across France. The four HER2s that are used for illustration in this study are highlighted in purple (see Sect. 3.7).</figcaption></figure>

within the year, except in cases of significant rainfall events that can temporarily trigger flow conditions. Given this inertia, the estimation of flow intermittence remains appropriate even if not all sites within a HER2 are monitored on the exact same day, with a shift of only a few days during a given campaign. To further ensure the temporal consistency across campaigns the monitoring is systematically conducted at fixed sites, and data from a given ONDE campaign and HER2 are only considered if at least 75% of the ONDE sites within the HER2 are monitored during the campaign. This approach minimizes potential biases associated with incomplete observations, and ensures that the percentage of dry ONDE sites referred to hereafter as PFI<sub>h</sub>(j) can be considered representative of the PFI.

Subsequently, a logistic regression model is calibrated for each HER2 in order to link the PFI values for day j to stream-flow conditions (Eq. 1, Fig. 2):

$$ \text{PFI}_h(j) = \frac{e^{\beta_0 \cdot \text{HER2}_h + \beta_1 \cdot \text{HER2}_h \times F_{\text{Q-HER2}_h}(j)}}{1 + e^{\beta_0 \cdot \text{HER2}_h + \beta_1 \cdot \text{HER2}_h \times F_{\text{Q-HER2}_h}(j)}}, \quad (1) $$

where $\beta_0 \cdot \text{HER2}_h$ and $\beta_1 \cdot \text{HER2}_h$ are, respectively, the logistic regression intercept and slope coefficient associated with the predictor $F_{\text{Q-HER2}_h}$.

The model uses the mean non-exceedance frequencies of discharge $F_{\text{Q-HER2}_h}(j)$ as a predictor for PFI<sub>h</sub>(j) (Eq. 2). $F_{\text{Q-HER2}_h}$ is regarded as a proxy characterizing the current wet versus dry hydrological conditions at the regional scale. It is derived from the non-exceedance frequency curves of streams, which describe the probability that a given discharge value is not exceeded. For each specific (j) corresponding to an ONDE field campaign and each HER2 (h), the daily empirical non-exceedance frequencies of discharge are spatially averaged over all available n streams whose drainage areas intersect the HER2 of interest (weighted by their drainage area) and temporally averaged over the period [j - 6; j]. This 7 d window allows the inclusion of non-simultaneous response times caused by heterogeneous propagation times in the underground and the river network (the choice of a 7 d window is the result of an optimization process; see Appendix B). In the end, $F_{\text{Q-HER2}_h}$ is defined by

$$ F_{\text{Q-HER2}_h}(j) = \frac{1}{7 \times n} \sum_{k=j-6}^{j} \sum_{s=1}^{n} F_{\text{Q}}(k, s), \quad (2) $$

where $F_{\text{Q}}(k, s)$ represents the non-exceedance frequency of each stream s among the n streams whose drainage areas intersect the HER2 h. This frequency is weighted by the
drainage area of stream *s*. The summation is computed over each day *k* within the [*j* − 6; *j*] period around campaign *j*.

Once calibrated, the logistic regressions are used to con-
vert the hydrological projections from the Explore2 project
into PFI projections.

## 3.2 Calibrating PFI logistic regression models using observed daily discharge data

Logistic regression models are preliminarily fitted between 2012 and 2022 using observed discharge data. The daily discharge data used in the calculation of $F_{Q-HER2_h}$ were compiled using gauging station records from the French reference hydrologic network HydroPortail (Leleu et al., 2014) (https://www.hydro.eaufrance.fr, last access: 4 September 2023). The set of gauging stations were selected on the basis of three different criteria:

- Daily discharge data must be available during the ONDE field campaigns, i.e. within the [*j* − 10; *j* + 10] interval centred on the 25th of each month from May to September for the years 2012 to 2022.

- The human disturbance to flow must be minimal at the gauging stations.

- All selected gauging stations must have a drainage area of less than 2000 km². Large gauged basins with a high Strahler order have been discarded since they are unlikely to behave in the same way as small streams.

This selection process results in a final dataset including 1008 stations (Fig. 3). The non-exceedance frequency curve of each gauging station is computed using its entire available discharge time series, leading to variations in the time periods covered by the discharge data across stations. Nevertheless, limiting the analysis to the period 2012–2022 would likely exclude critical information on hydrological extremes and fail to fully capture the variability in high- and low-flow events, which are pivotal for understanding a river’s behaviour during extreme conditions.

3.3 Assessing the performance of PFI logistic regression models under current climate

Several cross-validation schemes are then considered to test
the robustness of these models. A leave-one-year-out cross-
validation is first carried out by excluding one year at a time
from the training dataset, learning from the data of the re-
maining years and then making estimations for the excluded
year. The robustness under different climate conditions is
also assessed through a differential split-sample test (DSST;
Dakhlaoui et al., 2017), which is a k-fold cross-validation
performed on three distinct groups of hydrological years
categorized as dry, intermediate, or wet (Appendix C). The
years are discriminated according to the aridity index (Bar-
row, 1992), which is one on the main drivers of flow inter-
mittence at the global scale (Sauquet et al., 2021).

The model performance is assessed using several skill
scores: the bias (to detect underestimation or overestima-
tion), the mean absolute error (MAE, to quantify the mag-
nitude of prediction errors), and the root-mean-square error
(RMSE, which provides a quadratic estimation that assigns
relatively higher weight to large errors). Additionally, the
Nash–Sutcliffe model efficiency (NSE) coefficient is com-
puted to compare the variance of the estimation error to the
variance of the observed time series (Nash and Sutcliffe,
1970). The Kling–Gupta efficiency (KGE) complements the
NSE metric as KGE is the Euclidean distance between ob-
served and estimated PFI, computed using the coordinates of
bias, standard deviation, and correlation (Gupta et al., 2009).
The leave-one-year-out analysis results are obtained by aver-
aging the validation metrics computed for each year.

3.4 Selecting simulation points for PFI projections
under climate change

The modelling framework suggested by Beaufort et
al. (2018) and Sauquet et al. (2021) was applied to the 75
HER2s under regional climate projections to assess the po-
tential impact of climate change on flow intermittence dy-
namics.

Discharge data derived from the Explore2 database are
used as input for the logistic regression models introduced
in Sect. 3.1. However, not all the stations in the HydroPortail
database have been simulated by the hydrological models.
To overcome this co-location problem, the choice is made to
identify the simulation points closest to the 1008 gauging sta-
tions selected for the diagnosis in Sect. 3.2 (Fig. 3). The con-
ditions of the application experiment are thus similar to the
ones of the diagnostic experiment. A constraint is applied to
ensure that no simulation point is assigned multiple times un-
less there is no alternative point within a 50 km distance from
the gauging station (Appendix D). This selection is made in-
dependently for each hydrological model. Only 483 points
are assigned for the J2000 model because this model simu-
lates streamflow only for the Loire and Rhône river basins.

3.5 Constructing synthetic non-exceedance frequency curves for PFI projections

In our framework, non-exceedance frequency curves must be
computed for each simulation point prior to PFI estimation.
To avoid problems of extrapolating the non-exceedance fre-
quency of discharge values outside the range of historical
flows when performing PFI projection under changing cli-
mate, we decided to build synthetic flow time series with a
period length of 169 years that combine the historical dis-
charge data simulated with SAFRAN over the period 1976–
2022 with discharge data simulated by the different RCP–
GCM–RCM–hydrological models modelling chains over the
period 1976–2100. This concatenation is made indepen-
dently for each modelling chain and simulation point. As a
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Schematic view of the approach adapted to derive the regional probability of flow intermittence (PFI) from ONDE sites using daily flows either gauged or modelled within a given HER2. The daily flows simulated by Explore2 at three simulation points generate flow time series (column 1) and their corresponding non-exceedance frequency curves (column 2). A logistic regression is then used to link daily flows to ONDE observations, as illustrated here for the campaign of 25 September 2012. For this campaign, daily discharge values simulated between 18 and 25 September 2012 are extracted for each simulation point whose drainage area intersects the HER2. These seven flow values are matched to their respective non-exceedance frequency using the non-exceedance frequency curve. The mean non-exceedance frequency across all simulation points is then associated with the proportion of ONDE sites exhibiting “dry conditions” during this campaign (observed PFI, column 3). The logistic regression is calibrated using data from all ONDE campaigns. Once calibrated, these models allow for the conversion of daily discharge time series into daily PFI time series for the HER2s (column 4).</figcaption></figure>

consequence, while the assessment of the performance of PFI logistic regression models under the current climate is conducted using observed discharge data that can be unbalanced (i.e. covering different observation periods) due to heterogeneous data availability across stations (see Sect. 3.3), the calibration used for projecting PFI under future climate scenarios relies on standardized and balanced synthetic time series which ensures the consistency of the modelling across simulation points and modelling chains.

Finally, each synthetic discharge time series generated by a specific RCP-GCM-RCM-hydrological models chain results in a non-exceedance frequency curve which is used thereafter to calibrate a logistic regression (Eq. 1).

## 3.6 Recalibrating logistic regressions and projecting PFI under climate scenarios

Before being applied to future climate conditions, the logistic regression model of Eq. (1) needs to be recalibrated with past
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Location of the 3302 ONDE sites, the 1008 gauging stations selected from the HydroPortal database, and the Explore2 simulation points selected for the GRSD and J2000 hydrological models (line 1). The maps on the second line show the ratio between the sum of the catchment areas defined by the gauging stations or the simulation points and the surface area of the HER2 they intersect (line 2). These catchment areas may overlap, and their sums are computed without considering their spatial distribution. The maps on the third line show the density of gauging stations or simulation points (line 3). The same results are available for the other hydrological models in Appendix E.</figcaption></figure>

discharges simulated by each hydrological model (here discharges simulated by the different hydrological models using SAFRAN reanalysis data available during the period 2012–2022 are used). The resulting logistic regressions are then re-evaluated for each hydrological model using the median of the bias, MAE, RMSE, NSE, and KGE skill scores over all RCP-GCM-RCM modelling chains.

Finally, after re-calibration, these logistic regressions are used to project PFI based on discharge data simulated for both the historical period (1976–2004) and the projection period (2005–2100) under various climate scenarios (RCP2.6, RCP4.5, and RCP8.5).

## 3.7 Assessing changes in PFI and associated uncertainties

We compare PFI between a reference period (1976–2005) and medium-term (2041–2070) and long-term (2070–2099) horizons. Results are presented at the HER2 scale across France and discussed in more details for four HER2s (Fig. 1), deliberately selected to represent contrasting hydro-climatic

conditions and behaviours: “Haute Normandie Picardie” (HER2-57, North, oceanic temperate climate), “Plaine de Bourgogne” (HER2-81, Northeast, continental temperate climate), “Massif de l’Oisans” (HER2-12, Alps, mountainous climate), and “Plaine méditerranéenne” (HER2-105, Southeast, Mediterranean climate). In the rest of the article, “dry periods” of HER2s are defined as the period of the year when PFI is greater than 20 %. For each HER2, the changes in intensity and seasonality of flow intermittence for 2041–2070 and 2070–2099 relative to 1976–2005 are characterized using the mean daily PFI between July and October (mPFI$_{7–10}$) and the median of the first and last days (respectively, Tf and Tl) of the year with PFI exceeding 20 %.

The analysis of the uncertainty propagation in the modelling framework is examined using the QUALYPSO method, applied to the national average of mPFI$_{7–10}$ (mPFI$_{7–10}$[France]), weighted by HER2 areas and obtained under RCP2.6, RCP4.5, and RCP8.5 climate change scenarios (Evin et al., 2019, 2021). The QUALYPSO method is used to characterize the changes in mPFI$_{7–10}$[France] by estimating its ensemble mean over all the projections, as well
as the total uncertainty associated with this set of projections, the contribution of the various sources of uncertainty (RCP, GCM, RCM, hydrological models, residual and internal climate variability) to the total uncertainty, and the main effect of each individual model and its contribution to the total uncertainty. The contribution of each component to the total uncertainty is estimated as its fixed effect in a linear regression model, which adjusts for imbalances in the number of simulations across components within modelling chains (Appendix F).

A multi-model index of agreement (MIA) is also computed on mPFI<sub>7–10</sub> time series to highlight convergence in changes over the projections (Tramblay and Somot, 2018):

$$MIA = \frac{1}{n} \sum_{k=1}^{n} i_k, \qquad (3)$$

where *n* is the number of projections, $i_k = 1$ for a significant positive trend of the mPFI<sub>7–10</sub> according to the Mann–Kendall test (with $\alpha = 0.1$), $i_k = -1$ for a significant negative trend, and $i_k = 0$ for no significant trend.

# 4 Results

## 4.1 Data pre-processing

The 3302 ONDE sites are located in small streams, while gauging stations in France are monitoring medium-size catchments (Van Meerveld et al., 2020). In addition, all Explore2 simulation sites have a drainage area larger than 64 km². Unsurprisingly, there is little overlap between the distribution of areas drained by ONDE sites (median: 24 km²; quartiles Q1–Q3: 12–50) and those drained by the two sets of 1008 gauging stations (median: 173 km²; Q1–Q3: 85–396) and 1008 simulation points (median: 178 km²; Q1–Q3: 95–400; Appendix Fig. G1). ONDE sites, gauging stations, and simulation points are located at similar elevations (overall median: 168 m; Q1–Q3: 75–308; Appendix Fig. G2). The drainage areas of ONDE sites, gauging stations, and simulation points represent respectively a median coverage of 21% (Q1–Q3: 15–35), 54% (Q1–Q3: 31–94), and 55% (Q1–Q3: 32–96) of the HER2 areas. This corresponds to a median number of 6.1 sites (Q1–Q3: 5–8), 3.6 stations (Q1–Q3: 1.9–5.7), and 3.7 simulation points (Q1–Q3: 2–6) per 1000 km² of HER2 (Fig. 3). These statistics highlight the question of representativeness of the ONDE sites despite the even coverage of the ONDE network.

## 4.2 Model performance

The logistic regressions calibrated using observed flows accurately predict PFI under current climate conditions (Fig. 4; Appendix Fig. H1). The median performance across HER2s is, for example, 0.85 for KGE (Q1–Q3: 0.77–0.89) and 0.79 for NSE (Q1–Q3: 0.64–0.84) when the leave-one-year-out cross-validation is considered. The KGE skill score exceeds

80% in 50 out of 75 HER2s, and the bias is very low. The best performances are observed in sedimentary plains and in Aquitaine, while weaker performances are obtained in mountainous areas such as the Alps, the Massif Central and the Pyrenees. The mean absolute error and RMSE show that the largest absolute deviations between observations and predictions occur in the southeastern part of France. This result is consistent with previous findings (e.g. Fig. 3 in Sauquet et al., 2021): (1) the performance of the logistic regressions is partly correlated with the level of intermittence, and (2) no-flow conditions which mainly occur in winter due to freezing in the Alps are not captured by the ONDE network. Observed values of mPFI<sub>7–10</sub> range between 3% and 37% (median: 14.3%) across France, while values of mPFI<sub>7–10</sub> obtained with SAFRAN range between 2% and 37% (median: 14.4%). The modelling approach is able to reproduce the spatial pattern of flow intermittence (Fig. 4), although the results obtained with the DSST show less satisfactory performance (Appendix Fig. H2 and Table H1). The two skill scores KGE and NSE are sensitive to wet, intermediate, and dry conditions, while the other skill scores show no major change. The values of RMSE and MAE skill scores do not increase much when the calibration dataset is stratified based on climate conditions, which indicates that the proposed model is quite robust against climate fluctuations. The bias, however, indicates an underestimation of PFI in the southwest and in the north during dry years. This lower performance is partly due to a reduced number of years used for calibration purposes and to the difficulty in extrapolating the logistic regression equations under unobserved climate conditions.

The logistic regressions fitted with discharge data simulated by hydrological models also accurately reconstruct PFI from May to September. Compared to the results obtained with observed discharge data, the performance is slightly lower due to the inherent imperfections of the SAFRAN reanalysis compared to the observed climate (Figs. 5, 6; Table 1; Appendix I). The inter-HER2 and inter-projection median estimated mPFI<sub>7–10</sub> across France is 11.3%–13.5% depending on hydrological models, which is close to the 14.3% value (Q1–Q3: 8.4–20.3) derived from the ONDE network. The KGE is slightly lower but exceeds 80% in more than half of the HER2s using the CTRIP (38/75), SMASH (38/75), and GRSD (43/75) models. The inter-HER2 median of the NSE skill score (calculated for all RCP–GCM–RCM modelling chains and for each of these three hydrological models) indicates that between 71% and 74% of the variance in mPFI<sub>7–10</sub> is explained by the logistic regressions. The other skill scores do not show a loss of performance (KGE between 0.80 and 0.81, MAE between 0.04 and 0.05, and RMSE between 0.06 and 0.07) when comparing the results of the leave-one-year-out cross-validation obtained using observed and simulated discharge data. Lower performance is observed with the ORCHIDEE and J2000 models, as the KGE exceeds 80% in only 7/75 and 5/37 HER2s, respectively.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Leave-one-year-out cross-validation assessing the calibration of the logistic regressions using discharge measurements from gauging stations. From left to right: observed mean PFI measured during the calibration period (2012–2022, a), Kling–Gupta efficiency (KGE, b), and absolute bias of PFI predictions (c).</figcaption></figure>

tively. Their median KGE (0.60 and 0.69, respectively), NSE (0.49 and 0.62), MAE (0.06), and RMSE (0.08 and 0.09) are lower than the values obtained during calibration using observed flow data, indicating that these two hydrological models seem less sensitive.

Regardless of the hydrological model of interest, the maps of skill scores show spatial patterns that are consistent with the results obtained with observed discharge data (Fig. 6; Appendix I). Both NSE and KGE remain high in the sedimentary plains and the southern part of the country, while they are less favourable in mountainous regions (Alps, Massif Central, and Pyrenees). An overestimation persists in the southwestern part of France, while underestimations of PFI values are found in the northwestern, northern, and southeastern parts of France and are especially pronounced during dry years. The models also successfully reproduce the interannual variability in PFI, particularly the alternation of dry (e.g. 2019) and wet (e.g. 2015) years (Fig. 5). Thus, the modelling chains demonstrate their ability to simulate PFI over the period 2012–2022 and are therefore deemed reliable to attempt projecting changes in flow intermittence under modified climate conditions. However, note that PFI values are slightly underestimated when the logistic regressions are calibrated on the driest years and then applied to wetter years. This suggests that our future projections are likely to be biased in the same way, since the calibration is based on current conditions before applying the regression models to a potentially drier future climate.

## 4.3 Flow intermittence projections

The logistic regressions fitted for each HER2 are applied to derive the PFI values from the historical (1976–2004) and future (2005–2100) discharges simulated by the hydrological

models forced by RCP–GCM–RCM climate projections over the period 1976–2100.

Figure 7 shows the evolution of mPFI<sub>7–10</sub> and of the start date (Tf) and end date (Tl) of the dry period for 2041–2070 and 2070–2099 compared to 1976–2005. The results are illustrated in Fig. 7 with the median projection given by the GRSD model (see Appendix Figs. J1–J4 for other models) under RCP2.6, RCP4.5, and RCP8.5 climate scenarios. Table 2 details the same results for the four HER2s under RCP8.5.

Overall, the average PFI shows an increase, indicating a gradual intensification of dry periods over the 21st century under the RCP4.5 and RCP8.5 scenarios. The inter-HER2 and inter-projection median of mPFI<sub>7–10</sub> calculated at the national scale (mPFI<sub>7–10</sub>[France]) for the baseline period 1976–2005 ranges from 10% to 12% depending on the hydrological model. However, the median projection of mPFI<sub>7–10</sub> during 1976–2005 is higher than 20% in 12 to 20 of the 75 HER2s (depending on the hydrological model), although it does not exceed 42% in any HER2. Under RCP8.5, mPFI<sub>7–10</sub>[France] could reach values between 13% and 21% by the middle of the century 2041–2070 and between 16% and 29% by the end of the century 2070–2099. By this time, six HER2s could have inter-projection medians of mPFI<sub>7–10</sub> exceeding 50% according to at least two hydrological models: four in the southeast of France, one in the southwest, and one in the northwest of France. Under RCP4.5, changes are more moderate by the end of the century 2070–2099, with mPFI<sub>7–10</sub>[France] ranging between 14% and 20%. The spatial pattern of the changes is not uniform but looks similar between RCP8.5 and RCP4.5. A strengthening divide between the southwest and the northeast of the country is projected. Regions already prone to intermittence are expected to experience an increase in this
<table><caption>**Table 1.** Evaluation of the logistic regressions calibrated using observed discharge data (results of the cross-validation) and using flows simulated by the CTRIP, GRSD, J2000, ORCHIDEE, and SMASH hydrological models. Statistics of the skill scores are summarized by medians and the first and the third quartiles (Q1–Q3, into brackets) across all HER2s.</caption><thead><tr><th>Hydrological model</th><th>NSE (unitless)</th><th>KGE (unitless)</th><th>MAE (unitless)</th><th>RMSE (unitless)</th></tr></thead><tbody><tr><td>Observed discharge data</td><td>0.79</td><td>0.85</td><td>0.10</td><td>0.07</td></tr><tr><td>Leave-one-year-out (2012–2022)</td><td>(Q1–Q3: 0.64–0.84)</td><td>(Q1–Q3: 0.77–0.89)</td><td>(Q1–Q3: 0.06–0.13)</td><td>(Q1–Q3: 0.05–0.09)</td></tr><tr><td>CTRIP</td><td>0.74</td><td>0.80</td><td>0.05</td><td>0.07</td></tr><tr><td></td><td>(Q1–Q3: 0.61–0.80)</td><td>(Q1–Q3: 0.69–0.85)</td><td>(Q1–Q3: 0.04–0.07)</td><td>(Q1–Q3: 0.05–0.10)</td></tr><tr><td>GRSD</td><td>0.74</td><td>0.81</td><td>0.04</td><td>0.06</td></tr><tr><td></td><td>(Q1–Q3: 0.64–0.80)</td><td>(Q1–Q3: 0.71–0.86)</td><td>(Q1–Q3: 0.03–0.06)</td><td>(Q1–Q3: 0.04–0.08)</td></tr><tr><td>J2000</td><td>0.62</td><td>0.69</td><td>0.06</td><td>0.09</td></tr><tr><td></td><td>(Q1–Q3: 0.52–0.71)</td><td>(Q1–Q3: 0.63–0.79)</td><td>(Q1–Q3: 0.04–0.09)</td><td>(Q1–Q3: 0.07–0.13)</td></tr><tr><td>ORCHIDEE</td><td>0.49</td><td>0.60</td><td>0.06</td><td>0.08</td></tr><tr><td></td><td>(Q1–Q3: 0.39–0.63)</td><td>(Q1–Q3: 0.50–0.72)</td><td>(Q1–Q3: 0.04–0.09)</td><td>(Q1–Q3: 0.06–0.11)</td></tr><tr><td>SMASH</td><td>0.71</td><td>0.80</td><td>0.04</td><td>0.06</td></tr><tr><td></td><td>(Q1–Q3: 0.65–0.80)</td><td>(Q1–Q3: 0.72–0.86)</td><td>(Q1–Q3: 0.03–0.06)</td><td>(Q1–Q3: 0.05–0.08)</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Time series of PFI for the four HER2s selected for illustrative purposes. Dots and lines are observed PFI derived from the ONDE network and estimated PFI computed by the logistic regressions using simulated discharge by each hydrological model, respectively.</figcaption></figure>

phenomenon under both emission scenarios. The ratio of the ensemble median of mPFI<sub>7–10</sub> at the end of the century compared to the baseline period revolves around 1.4 for RCP4.5 and 1.9 for RCP8.5. Mountainous regions see the intensity of summer dry periods increase but remain relatively spared.

Shifts of the start and the end of the dry period are partly correlated with changes in mPFI<sub>7–10</sub> (Fig. 7). It appears that

climate change results in both earlier and later dry periods in a fairly symmetrical manner: dry periods are advanced and extended but with different sensitivities structured along a north–south gradient. Under RCP8.5 climate conditions, the dry periods are projected to get (on median) longer in southern France at the end of the century. According to at least two hydrological models, the highest shifts of Tf or Tl could
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Kling–Gupta efficiency (KGE) assessing the calibration of logistic regression using discharge simulated by each hydrological model with SAFRAN-based simulations available during the 2012–2022 period.</figcaption></figure>

exceed 5 weeks in several mountainous HER2s (four HER2s in the Alps, three in the Pyrenees, two in the Jura). In contrast in the northern part of France, the dry periods could be advanced and extended by only 1 or 2 weeks, as observed in HER2-81 “Plaine de Bourgogne” (Table 2). HER2-28 could be one of the most affected in the north due to its impermeable sandy-clay formations, which differ from the neighbouring sedimentary formations. One can note a change in seasonality for the HER2-12 “Massif de l’Oisans” located in the Alps (Fig. 8, example with the GRSD hydrological model). No-flow events were concentrated in winter (retention of water in the snow cover) during the historical period. Under climate change, temperature will be higher leading to less snowfall and more runoff in the river network during winter. The river flow regime will be more sensitive to losses by evapotranspiration, and finally no-flow conditions will likely occur in summer by the end of the 21st century.

## 4.4 Uncertainties in hydro-climatic projections and their impact on PFI projections

Figure 9 illustrates the results of the uncertainty analysis, showing the uncertainty range of mPFI<sub>7-10</sub>[France] throughout the 21st century. These results confirm that dry conditions may occur more frequently in a changing climate, with a 22% projected increase in mPFI<sub>7-10</sub>[France] by the end of the century. The total uncertainty of mPFI<sub>7-10</sub>[France] over all projections results from the accumulation of uncertainties related to RCP scenarios, GCMs and RCMs, hydrological models, residual variability, and internal variability. The confidence interval, and therefore the extent of change, increases over the course of the 21st century due to the divergent results of the modelling chain. This finding is consistent with studies using a similar methodology to assess the uncertainty of panels of flow projections (Evin et al., 2019; Aitken et al., 2023).

The contribution of each component to the total uncertainty varies over time. The fraction of total variance due to residual variability is small and remains stable over time. The main part of the uncertainty in climate change responses is thus explained by RCP, GCM, RCM, and hydrological models. In addition, the contribution of the internal variability is

highly predominant and represents more than 75% of the total uncertainty until the mid-century. This contribution then decreases with time as the total uncertainty increases. It becomes less than 50% of the total uncertainty by the end of the century but remains the greatest contributor to the total uncertainty over the whole simulation period.

Regarding the spatial distribution of the different contributions to the total uncertainty, it appears that climate models are a predominant source of uncertainty and that they are distributed in a spatially heterogeneous manner over France. The uncertainty related to the RCP scenarios increases over time and becomes predominant in the south, surpassing uncertainties related to other steps of the modelling chain (GCM, RCM, and hydrological models). In the mountainous regions, differences between results obtained with RCP4.5 and RCP8.5 are clearer by the end of the century, leading to a contribution of RCP to the total uncertainty reaching 20% at the end of the century. In northern France, the contribution of RCP to the uncertainty is lower compared to hydrological model, GCM, or RCM contributions. Uncertainty related to hydrological models is also spatially structured, with greater uncertainty in the northwestern part of France.

## 4.5 Agreement between changes

Despite uncertainties about the intensity of change, PFI projections still show a good degree of convergence in the direction of change when considering all climate projections driving each hydrological model individually, as well as in the overall convergence of PFI projections (Fig. 10). MIA values are first calculated independently for each hydrological model before all projections are combined for multi-model assessment.

Figure 10 highlights spatial contrasts of the MIA multi-model approach: unsurprisingly, the projections of mPFI<sub>7-10</sub> are more uncertain in the northern part of France. The climate change signal on the proportion of dry periods is not homogeneous across France, with significant uncertainties remaining for the northern part.

The agreement of projections is high for the GRSD and SMASH models across France (MIA values close to +100% under RCP8.5 in Fig. 10 and +80% under RCP4.5; see Ap-
<table><caption>**Table 2.** Statistics of flow intermittence characteristics for the four illustrative HER2s under RCP8.5. The minimum (min), median (med), and maximum (max) are the ensemble minimum, median, and maximum across all the projections and hydrological models. Nb Days PFI >20 %: number of days in dry periods (period of the year when PFI exceeds 20 %).</caption><thead><tr><th rowspan="2">HER2</th><th rowspan="2">Period</th><th colspan="3">mPFI<sub>7-10</sub> (%)</th><th colspan="3">Tf (date, dd/mm)</th><th colspan="3">Tl (date, dd/mm)</th><th colspan="3">Nb Days PFI >20 %</th></tr><tr><th>Min</th><th>Med</th><th>Max</th><th>Min</th><th>Med</th><th>Max</th><th>Min</th><th>Med</th><th>Max</th><th>Min</th><th>Med</th><th>Max</th></tr></thead><tbody><tr><td rowspan="3">Haute Normandie<br>Picardie (57)</td><td>1976-2005</td><td>8</td><td>13</td><td>15</td><td>30/07</td><td>22/08</td><td>05/10</td><td>22/09</td><td>30/10</td><td>04/12</td><td>3</td><td>33</td><td>40</td></tr><tr><td>2041-2070</td><td>7</td><td>14</td><td>23</td><td>15/07</td><td>24/08</td><td>11/10</td><td>11/10</td><td>10/11</td><td>16/12</td><td>1</td><td>38</td><td>98</td></tr><tr><td>2070-2099</td><td>5</td><td>16</td><td>34</td><td>04/06</td><td>21/08</td><td>23/11</td><td>16/10</td><td>12/11</td><td>23/12</td><td>0</td><td>37</td><td>185</td></tr><tr><td rowspan="3">Plaine de<br>Bourgogne (81)</td><td>1976-2005</td><td>12</td><td>19</td><td>22</td><td>10/07</td><td>31/07</td><td>18/08</td><td>18/09</td><td>27/10</td><td>16/12</td><td>31</td><td>65</td><td>83</td></tr><tr><td>2041-2070</td><td>11</td><td>23</td><td>34</td><td>09/07</td><td>02/08</td><td>01/09</td><td>02/10</td><td>04/11</td><td>29/12</td><td>31</td><td>79</td><td>140</td></tr><tr><td>2070-2099</td><td>10</td><td>27</td><td>47</td><td>21/06</td><td>22/07</td><td>03/09</td><td>04/10</td><td>08/11</td><td>31/12</td><td>31</td><td>93</td><td>189</td></tr><tr><td rowspan="3">Massif de<br>l'Oisans (12)</td><td>1976-2005</td><td>6</td><td>11</td><td>17</td><td>26/07</td><td>23/08</td><td>06/10</td><td>18/09</td><td>19/11</td><td>17/03</td><td>23</td><td>32</td><td>128</td></tr><tr><td>2041-2070</td><td>11</td><td>21</td><td>38</td><td>27/06</td><td>01/08</td><td>01/09</td><td>23/09</td><td>30/10</td><td>03/12</td><td>25</td><td>65</td><td>139</td></tr><tr><td>2070-2099</td><td>12</td><td>29</td><td>52</td><td>29/05</td><td>12/07</td><td>21/08</td><td>28/09</td><td>02/11</td><td>15/12</td><td>38</td><td>92</td><td>197</td></tr><tr><td rowspan="3">Plaine méditerranéenne<br>(105)</td><td>1976-2005</td><td>14</td><td>25</td><td>29</td><td>22/06</td><td>17/07</td><td>16/08</td><td>22/09</td><td>23/10</td><td>17/12</td><td>37</td><td>85</td><td>103</td></tr><tr><td>2041-2070</td><td>9</td><td>31</td><td>46</td><td>29/05</td><td>11/07</td><td>28/08</td><td>25/09</td><td>01/11</td><td>25/12</td><td>28</td><td>107</td><td>168</td></tr><tr><td>2070-2099</td><td>9</td><td>41</td><td>57</td><td>02/04</td><td>16/06</td><td>30/08</td><td>02/10</td><td>13/11</td><td>01/01</td><td>33</td><td>155</td><td>228</td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Ensemble median mPFI<sub>7-10</sub> (columns 1 to 3) and ensemble median shift of the start date Tf (columns 4 and 5) and the end date Tl (columns 6 and 7) of the dry period over the two periods 2041–2070 and 2070–2099 under the three RCPs for the GRSD hydrological model, relative to the baseline period 1976–2005. The shift, expressed in weeks, takes a positive value when the duration of the dry period increases. Grey HER2s had no period with a PFI >20 % during the reference period. The same results are available for the other hydrological models in Appendix Figs. J1–J4.</figcaption></figure>

Appendix K). In contrast, CTRIP, J2000, and ORCHIDEE suggest contrasted regional impacts of climate change on PFI under RCP4.5, with a reduction in the mPFI<sub>7-10</sub> in the northern part of France, indicating a differentiated sensitivity of

hydrological models to climate changes. This trend is also observed for ORCHIDEE under RCP8.5, while the other models agree on an increasing drying.
Average proportion of zero flows smoothed over 5 d
Hydrological model: GRSD - HER2 12

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

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Seasonal pattern of PFI for the periods 1976–2005, 2041–2070, and 2070–2099 under RCP2.6, RCP4.5, and RCP8.5 scenarios for the GRSD hydrological model in the Alps (HER2-12). The storylines highlighted in this figure are introduced in Appendix J. The same results are available for HER2-57, HER2-81, and HER2-105 and for the other hydrological models in Appendix Figs. J5–J21.</figcaption></figure>

# 5 Discussion

## 5.1 Modelling framework and assumptions

This study presents projections of the PFI in small streams by HER2s in France over the 21st century. It extends the previous analysis of Sauquet et al. (2021) on the 2012–2018 period by exploiting an extensive national dataset of stream intermittence observations collected on low-order streams (3302 observation sites), which are the most affected by shifts from perennial to intermittent flows (Reynolds et al., 2015; Dhungel et al., 2016). For the first time, we project the future evolution of the PFI using discharge data obtained from 5 hydrological models with conceptual (GRSD, SMASH), surface (CTRIP, ORCHIDEE), or process-oriented distributed (J2000) structures, under a multitude of possible climate scenarios informed by 17 pairs of CMIP5 GCM–RCM models.

The results of this study should be interpreted within the context of its underlying assumptions. Firstly, the study fo-

cuses on unregulated streams to characterize the “natural” hydrology (i.e. without considering water abstraction by anthropogenic activities or the impact of hydraulic engineering structures). This assumption was made in the Explore2 project, which provides the input streamflow simulations. Nevertheless, global water model simulations including direct human impacts by Döll and Zhang (2010) or Gudmundsson et al. (2021) concluded that ecologically relevant flow characteristics will be more altered by climate change than by withdrawals and dams. However, we believe that quantitatively estimating the extent to which flow intermittence is due to direct anthropogenic stressors is crucial. Such an estimation could improve our projections of flow intermittence and inform decision-making aimed at regulating and preventing water stress situations for populations.

Secondly, groundwater levels are not incorporated into the model, although they could potentially enhance the accuracy of the projections. Indeed, similar logistic regression models were used in Beaufort et al. (2018), and they incorporated
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Decomposition of the effects contributing to the variance of projections for relative changes (unitless) of mPFI<sub>7-10</sub>[France]. In the top-left graph, the grey line and the coloured areas respectively show the climate change response and the contribution of time for each step of the modelling chain to the 90% confidence interval. For each model uncertainty and internal variability component, the vertical extent of the corresponding area is proportional to the fraction of total uncertainty explained by the component.</figcaption></figure>

groundwater levels measured by piezometers. Likewise, in Beaufort et al. (2019), climate, hydrological, and morphological data were supplemented with groundwater levels to predict flow intermittence at a local scale. In these studies incorporating subsurface water in the regression models was thought to improve performance, in particular because the contribution of subsurface processes is known to play a decisive role in maintaining baseflow and mitigating flow intermittence. However, information about groundwater levels has not been included in the present study because hydrological projections incorporating groundwater levels have not yet been conducted across the entirety of France (for instance, within the Explore2 project, projections of groundwater levels have only been produced for a restricted set of areas in France).

Thirdly, with only five annual discrete observations of streamflow intermittence over 11 years for the calibration of the logistic regression models, our ability to capture the full range of extremes and variability in these regression models remains imperfect. Yet, visual monitoring remains the most

common technique for observing non-perennial streams. An alternative approach would be to consider citizen science to augment our database, although concerns about data reliability persist, in particular because past studies have shown participant agreement rates ranging from 46% to 70% (Scheller et al., 2024). Data scarcity necessitated conducting this analysis at the HER2 scale, which nonetheless represents a significant improvement of the spatial resolution by increasing the number of modelling subdomains from 22 to 75 compared to previous studies (Sauquet et al., 2021). In future work, a downscaling process could enhance the usability of PFI projections for local stakeholders in water management.

Fourthly, while projections of no-flow events often exhibit significant uncertainties when derived directly from discharge simulated by hydrological models (Evin et al., 2024; Aitken et al., 2023), this study illustrates how categorical and discrete (in space and time) field observations can be combined with conventional stream gauge data to enhance the understanding of small-stream drying dynamics. In this way, despite the divergences between projections induced by cli-
<figure><img src="image_13.png" /><figcaption>**Figure 10.** Agreement between projections of mPFI<sub>7-10</sub> for each hydrological model and inter-model agreement on the change signal of mPFI<sub>7-10</sub> under the RCP8.5 scenario.</figcaption></figure>

mate data, the agreement within this multi-model approach is relatively strong, indicating a consensus toward increasing PFI throughout the 21st century.

## 5.2 A consistent signal across projections

This study confirms that logistic regressions properly capture the relationships between flows in large watersheds and the PFI of small streams. These regressions are calibrated and validated using discharge measurements from gauging stations and subsequently using discharge simulated from SAFRAN climate data. Following this second calibration, projections based on regressions consistently indicate an increase in average PFI and a shift in the start and end dates of dry periods under RCP4.5 and RCP8.5 climate projections, suggesting a progressive intensification and extension of dry periods over the 21st century.

These results are consistent with previous studies indicating a transition of many streams from perennial to intermittent regimes (Jaeger et al., 2014; Reynolds et al., 2015; Dhungel et al., 2016; Schneider et al., 2013). In regions already affected by intermittence, an increase in the intensity and duration of drying periods is likely, a phenomenon also anticipated in other areas around the world (Jaeger et al., 2014; Dai et al., 2011). Increasing intermittence and decreasing flow rates are observed in catchment-scale studies: the fraction of

time with zero discharge increases from 0.05 % to 4.30 % in a Swiss alpine catchment between 2020–2040 and 2080–2100 (Halloran et al., 2023), while the no-flowing phase could extend by up to 12 d between 1980–2009 and 2030–2059 for a catchment in southern Italy (De Girolamo et al., 2022). This trend is also noticeable on a larger scale, as the monitoring of gauging stations both in five regions with Mediterranean climates around the world between 1980 and 2019 (Carlson et al., 2024) and over 452 rivers on the European continent between 1970 and 2010 (Tramblay et al., 2020) shows that approximately 30 % of them have already experienced drier conditions due to climate change, with modified flow regimes or extended periods of drying.

## 5.3 Uncertainties in northern France

In the northern part of France, discrepancies between GCM-RCM-hydrological models projections result in higher uncertainties in PFI projections and pronounced geographical contrasts on the multi-model-based MIA map (Fig. 10). Part of these uncertainties can be attributed to the underestimation of PFI in this region, particularly during dry years (Sect. 4.2, Appendix Fig. H2, and Table H1). However, the primary source of these uncertainties likely stems from uncertainties in future rainfall patterns in this region where the majority of Explore2 projections indicate an increase in winter rainfall
and winter mean flows (for 8 out of 9 hydrological models) and a decrease in summer precipitation (Evin et al., 2024; Sauquet et al., 2024). The annual precipitation by the end of the century remains uncertain due to the compensatory effect between increased winter recharge and an increased evapotranspiration (Ribes et al., 2022; Douville et al., 2021). As a result, some hydrological models predict a shortening of the dry period for certain northern HER2s, including under RCP8.5. Similar uncertainties are observed along the east coast of the USA, where precipitation changes could turn intermittent streams into perennial ones (Dhungel et al., 2016), while alternative climate change scenarios projected that by the 2040s, approximately half of the streams in Washington State would shift from snow-fed to rain-fed, resulting in reduced annual discharge (Reidy Liermann et al., 2012).

## 5.4 Transformation of the snowmelt regime in the Alps

The mountain ranges could be moderately affected by the increased probability of dry conditions, consistent with our previous modelling efforts (Sauquet et al., 2021), but our projections anticipate two specific phenomena in these regions. First, the Pyrenees are more affected than the other mountain ranges, with significant changes impacting the massif and the dependent basins in southwestern France. Additionally, alpine HER2s will probably undergo hydrological regime changes. Higher winter temperatures will lead to reduced snowfall and snow-related intermittence. This snowpack reduction will also reduce groundwater recharge by spring melt. In addition, since soils tend to retain less water in summer due to more intense and brief rainstorms (Rutkowska et al., 2023), an increased summer intermittence is expected despite increasing summer precipitations. Using 16 hydrographic variables describing the magnitude, frequency, timing, duration, and rate of change of the flow regime at 59 primarily selected sites with a Strahler order of 5, Dhungel et al. (2016) also observed the reduction in the snowmelt regime at two Rocky Mountain sites under climate change. Halloran et al. (2023) also conclude that groundwater will play an increasingly important role in ensuring flow in alpine streams and that the shift from perennial to intermittent could occur for alpine streams over the course of the current century.

Furthermore, the reduction in the snowmelt regime in the Alps indicates that we can still make projections consistent with the literature beyond the calibration period (May to September). More generally, the hypothesis of temporal transferability of the models is a strong assumption in data exploitation, as it assumes that climate models, statistical adjustment methods, and hydrological models can simulate the behaviour of the systems they represent in a future hydroclimatic context that is very different from the one in which they were developed (Evin et al., 2024). In this context, it is important to recall the risk that the projections may be underestimated, as demonstrated by the validation of logistic regressions calibrated on wetter years compared to dry years.

## 6 Conclusion

This study assesses the changes in the intermittency of river flows across France in the context of climate change. For the first time, multi-model and multi-scenario hydro-climatic projections are used as a predictor to explore the possible evolution of the daily probability of flow intermittency at the scale of (level-2) hydro-ecoregions. Leveraging monthly monitoring of small streams over 10 summers, we calibrated logistic regressions to transform hydrological projections of large watersheds into regional proportions of flow intermittence for the 21st century. Under both RCP4.5 and RCP8.5 scenarios, robust signals indicate an intensification of dry events, marked by an increased probability of flow intermittence and longer dry periods throughout the year. These changes are projected to be more pronounced in southern France, with greater uncertainty in the northern half of the country. Mountain areas could remain relatively spared from summer dry periods, but shifts in hydrological regimes are anticipated. By the end of the century under RCP8.5, dry-stream phenomena along the Atlantic coast could surpass those currently observed in the Mediterranean region by the usual monitoring campaigns. The evolution of droughts and reduced water availability suggested by these results could lead to significant ecological impacts, including alterations in the structure and function of freshwater ecosystems (such as changes in microbial activity and habitat loss), shifts in soil chemistry, increased carbon and solute fluxes, and sediment mobilization (Geris et al., 2015).
Appendix A: Explore2 modelling

<table><caption>**Table A1.** Global and regional climate model combinations driving the Explore2 hydrological models selected for PFI simulation in the 21st century.</caption><thead><tr><th>Global climate model</th><th>Regional climate model</th><th>Historical</th><th>RCP2.6</th><th>RCP4.5</th><th>RCP8.5</th></tr></thead><tbody><tr><td>CNRM-CERFACS-CNRM-CM5</td><td>CNRM-ALADIN6</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>CNRM-CERFACS-CNRM-CM5</td><td>MOHC-HadREM3-GA7-05</td><td>×</td><td>—</td><td>—</td><td>×</td></tr><tr><td>ICHEC-EC-EARTH</td><td>KNMI-RACMO22E</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>ICHEC-EC-EARTH</td><td>SMHI-RCA4</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>ICHEC-EC-EARTH</td><td>MOHC-HadREM3-GA7-05</td><td>×</td><td>×</td><td>—</td><td>×</td></tr><tr><td>MOHC-HadGEM2-ES</td><td>CNRM-ALADIN63</td><td>×</td><td>—</td><td>—</td><td>×</td></tr><tr><td>MOHC-HadGEM2-ES</td><td>CLMcom-CCLM4-8-17</td><td>×</td><td>—</td><td>×</td><td>×</td></tr><tr><td>MOHC-HadGEM2-ES</td><td>ICTP-RegCM4-6</td><td>×</td><td>×</td><td>—</td><td>×</td></tr><tr><td>MOHC-HadGEM2-ES</td><td>MOHC-HadREM3-GA7-05</td><td>×</td><td>×</td><td>—</td><td>×</td></tr><tr><td>IPSL-IPSL-CM5A-MR</td><td>DMI-HIRHAM5</td><td>×</td><td>—</td><td>—</td><td>×</td></tr><tr><td>IPSL-IPSL-CM5A-MR</td><td>SMHI-RCA4</td><td>×</td><td>—</td><td>×</td><td>×</td></tr><tr><td>MPI-M-MPI-ESM-LR</td><td>CLMcom-CCLM4-8-17</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>MPI-M-MPI-ESM-LR</td><td>ICTP-RegCM4-6</td><td>×</td><td>×</td><td>—</td><td>×</td></tr><tr><td>MPI-M-MPI-ESM-LR</td><td>MPI-CSC-REMO2009</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>NCC-NorESM1-M</td><td>DMI-HIRHAM5</td><td>×</td><td>—</td><td>×</td><td>×</td></tr><tr><td>NCC-NorESM1-M</td><td>GERICS-REMO2015</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>NCC-NorESM1-M</td><td>IPSL-WRF381P</td><td>×</td><td>—</td><td>—</td><td>×</td></tr></tbody></table>
## Appendix B: Sensitivity analysis of the duration of flow measurement interval used as input for the logistic regression to calculate the PFI

A sensitivity analysis was performed to fix the time window of daily discharge that optimizes the calibration of logistic regressions presented in Sect. 3.1. The calibration and validation were performed using the leave-one-year-out method, detailed in Sect. 3.3. The model performance was assessed using the Kling–Gupta efficiency (KGE) (Gupta et al., 2009). The KGE values for the different HER2s are summarized using boxplots for each tested window size. The [j-6; j] window was selected because it corresponds to the highest median KGE score, the narrowest interquartile range, and the highest minimum KGE.

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

**Figure B1.** Kling–Gupta efficiency (KGE) computed on results obtained by leave-one-year-out validations, testing the sensitivity to the time window of daily discharge used for calibrating the logistic regressions. Boxes represent the quartiles Q1 and Q3, the whiskers extend up to 1.5 times the interquartile range above Q3 and below Q1, and points located beyond the whiskers are displayed individually. The dashed line represents the median KGE of the [j-6; j] window.
## Appendix C: Separation of dataset based on dry, intermediate, and wet years

The data from ONDE sites used for validation were collected between May and September over 11 years from 2012 to 2022. The robustness of the model was assessed through a validation process involving three sets of dry, intermediate, and wet years. For each test, model calibration was performed using the years excluded from the validation set.

The hydrological years are distributed into three equal groups of hydrological years (dry years, intermediate years, and wet years) according to the annual aridity index, calculated at the national scale. The aridity index AI was given by the ratio between the total annual precipitation and potential evapotranspiration from 1 August of the previous year to 31 July of the current year (Barrow, 1992; Fig. C1). A set of dry years was formed using hydrological years where potential evapotranspiration exceeded annual precipitation (AI < 1 for 2017, 2019, and 2022). Two sets were then created: the four years with AI greater than 1.4 were classified as wet years (2012, 2015, 2020 and 2021), and the remaining four years with AI between 1.15 and 1.37 were classified as intermediate years (2013, 2014, 2016, 2018).

<figure><img src="image_15.png" /><figcaption>Figure C1. Aridity index for the hydrological years 2012 to 2022.</figcaption></figure>

## Appendix D: Assignment of simulation points

After establishing the models' reliability through validation at secure monitoring stations, the next step involved extrapolation, which entailed linking these monitoring sites to the nearest simulation points in Explore2 to project future hydrological scenarios while preserving the existing data structure. The proximity between gauging station A and a simulation

point in the Explore2 project is measured by the distance.

$$ \text{Dist}(A, B) = \sqrt{(X(A) - X(B))^2 + (Y(A) - Y(B))^2 + (\alpha \times \Delta S_{\text{rel}})^2} \quad (\text{D1}) $$

Here, $(X(A)Y(A))$ and $(X(B)Y(B))$ are coordinates of A and B (in km), respectively, and $\text{Surf}(A)$ and $\text{Surf}(B)$ are the drainage areas of A and B, respectively, that are used to compute the relative difference between the drainage areas (absolute value).

$$ \Delta S_{\text{rel}} = 2 \times \frac{|\text{Surf}(A) - \text{Surf}(B)|}{|\text{Surf}(A) + \text{Surf}(B)|} \quad (\text{D2}) $$

The coefficient $\alpha$ is used to balance the importance of geographical distance and the relative difference in surface, and it was set to 100.
## Appendix E: Location of the ONDE sites, the gauging stations, and the Explore2 simulation points selected for each hydrological model

<figure><img src="image_16.png" /><figcaption>**Figure E1.** Location of the 3302 ONDE sites, the 1008 gauging stations selected from the HydroPortail database, and the Explore2 simulation points selected for each hydrological model (line 1). Maps on the second and third lines respectively show the ratio between the sum of the catchment areas defined by the gauging stations or the simulation points and the surface area of the HER2 they intersect (line 2), and the density of gauging stations or simulation points (line 3).</figcaption></figure>

## Appendix F: QUALYPSO method

The QUALYPSO method is applied to the multi-scenario multi-model ensembles of projections of the national average of mPFI$_{7-10}$ (mPFI$_{7-10}$[France]), weighted by HER2 areas and obtained under RCP2.6, RCP4.5, and RCP8.5 climate change scenarios. According to Evin et al. (2019), we assume that the variable mPFI$_{7-10}$$_{i,j,k,l}$ (t) characterizes the change in mPFI$_{7-10}$ between a year t and the control year c (e.g. c = 1990 in this study) for a given combination of RCP scenario i, GCM j, RCM k, and hydrological model l:

$$m\text{PFI}_{7-10\text{i,j,k,l}}^{*}(t) = \text{mPFI}_{7-10\text{i,j,k,l}}(t) - \text{mPFI}_{7-10\text{i,j,k,l}}(c).$$

The change in mPFI$_{7-10}$ can be split up into

$$m\text{PFI}_{7-10\text{i,j,k,l}}^{*}(t) = \varphi_{i,j,k,l}^{*}(t) + \eta_{i,j,k,l}^{*}(t),$$

where $\varphi_{i,j,k,l}^{*}(t)$ is the climate change response of the RCP-GCM-RCM-hydrological models combination, and $\eta_{i,j,k,l}^{*}(t)$ is the deviation from the climate change response for this RCP-GCM-RCM-hydrological models combination as a result of internal variability, representing natural and stochastic variability in the climate system.

For each time t, the climate change response characterized by $\varphi_{i,j,k,l}^{*}(t)$ of any RCP-GCM-RCM-hydrological models combination can be expressed as

$$\varphi_{i,j,k,l}^{*}(t) = \mu(t) + \alpha_{i}(t) + \beta_{j}(t) + \gamma_{k}(t) + \theta_{l}(t) + \varepsilon_{i,j,k,l}(t),$$

where $\mu(t)$ is the ensemble mean climate change response, $\alpha_{i}(t)$ is the main effect of emission scenario i (i.e. the mean deviation of RCP scenario i from $\mu(t)$), $\beta_{j}(t)$ is the main effect of GCM j, $\gamma_{k}(t)$ is the main effect of RCM k, $\theta_{l}(t)$ is the main effect of the hydrological model l, and $\varepsilon_{i,j,k,l}(t) = \varphi_{i,j,k,l}^{*}(t) - \mu(t) - \alpha_{i}(t) - \beta_{j}(t) - \gamma_{k}(t) - \theta_{l}(t)$ corresponds to residual terms. For each year t, $\varepsilon_{i,j,k,l}(t)$ are assumed to be independent and identically distributed over all scenarios, GCMs, RCMs, and hydrological models and to follow normal distributions, with mean 0 and variance $\sigma^2(t)$.

If $\varphi_{i,j,k,l}^{*}(t)$ and $\eta_{i,j,k,l}^{*}(t)$ are assumed to be independent, the total variance of the change variable mPFI$_{7-10\text{i,j,k,l}}^{*}(t)$ is given by

$$\text{Var}(\text{mPFI}_{7-10\text{i,j,k,l}}^{*}(t)) = \text{Var}(\varphi_{i,j,k,l}^{*}(t)) + \text{Var}(\eta_{i,j,k,l}^{*}(t)),$$
where $Var(\eta_{i,j,k,l}^*(t))$ is the uncertainty associated with internal variability in mPFI$_{7-10i,j,k,l}^*(t)$.

$Var(\phi_{i,j,k,l}^*(t))$ is the uncertainty in the climate change response, calculated as the sum of different uncertainty components.

$$Var(\phi_{i,j,k,l}^*(t)) = Var(\mu(t)) + Var(\alpha_i(t)) + Var(\beta_j(t)) \\ + Var(\gamma_k(t)) + Var(\theta_l(t)) + Var(\varepsilon_{i,j,k,l}(t))$$

RCP uncertainty $Var(\alpha_i(t))$ quantifies the variability in mPFI$_{7-10i,j,k,l}^*(t)$ values across the RCP scenarios. Specifically, $\alpha_i(t)$ represents the differences between the means of mPFI$_{7-10i,j,k,l}^*(t)$ involving each scenario $i$ and the overall mean of mPFI$_{7-10i,j,k,l}^*(t)$ across all scenarios. This component illustrates how differences between emission scenarios contribute to the total uncertainty in the projected mPFI$_{7-10i,j,k,l}^*(t)$.

Similarly, GCM, RCM, and hydrological model contributions to mPFI$_{7-10i,j,k,l}^*(t)$ uncertainties are characterized by the variances for GCMs (Var($\beta_j(t)$)), RCMs (Var($\gamma_k(t)$)), and hydrological models (Var($\theta_l(t)$)).

Residual variability (Var($\varepsilon_{i,j,k,l}(t)$)) captures interaction effects and unexplained variability not accounted for by the main effects.

## Appendix G: Metadata

<figure><img src="image_17.png" /><figcaption>Figure G1. Frequency distribution of the number of ONDE sites, the 1008 gauging stations selected from the HydroPortal database, and the 1008 simulation points from Explore2 as a function of catchment area. NA: missing values.</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>**Figure G2.** Frequency distribution of the number of ONDE sites, the 1008 gauging stations selected from the HydroPortal database, and the 1008 simulation points from Explore2 as a function of elevation NA: missing value.</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>**Figure G3.** Strahler index of ONDE sites. NA: missing value.</figcaption></figure>
## Appendix H: Model performance using observed discharge data from gauging stations as predictors

Over the observation period (2012–2022), logistic regression models estimate the PFI values at HER2 scale (Fig. 4). The median observed PFI across all HER2s and all campaigns is 14.3% (Q1–Q3: 8.4–20.3), while the logistic regression models yield a median value of 14.4% (Q1–Q3: 8.5–20.4). With results of the leave-one-year-out cross-validation, the models explain 73% of the PFI variability according to the NSE (Q1–Q3: 64%–84%) (Table 1; Fig. H2), the KGE exceeds 80% in 50 out of 75 HER2s (Fig. H2), and the bias is very low.

The models are also able to describe the inter-annual variability with the alternation of dry and wet years (Fig. H1). The median KGE and NSE scores remain above 0.71 during k-fold validations considering dry, intermediate, and wet years for calibration. RMSE and MAE values do not increase much when the calibration dataset is stratified based on climate conditions, which indicates that the proposed model is quite robust against climate variations under current conditions.

<figure><img src="image_20.png" /><figcaption>**Figure H1.** Time series of PFI in HER2-13, HER2-57, HER2-81, and HER2-105. Points and lines respectively represent the PFI derived from the ONDE network and the PFI estimated by the logistic regression models using discharge data from gauging stations during the calibration period (2012–2022).</figcaption></figure>
<figure><img src="image_21.png" /><figcaption>**Figure H2.** Kling–Gupta efficiency (KGE, line 1), Nash–Sutcliffe model efficiency coefficient (NSE, line 2), bias (line 3), mean absolute error (MAE, line 4), and root-mean-square error (RMSE, line 5) over the calibration period (2012–2022).</figcaption></figure>
<table><caption>**Table H1.** Validation results of drying probability predictions at the HER2 scale using observed flows from the 1008 gauging stations from the HYDRO database. The results correspond to the inter-HER medians, quartiles, and minimum and maximum values. The leave-one-year-out analysis results are obtained by averaging the validation metrics computed for each year.</caption><thead><tr><th>Validation years</th><th>PFI observed on ONDE network (%)</th><th>Predicted PFI (%)</th><th>Nash–Sutcliffe efficiency (NSE, unitless)</th><th>Kling–Gupta efficiency (KGE, unitless)</th><th>Bias (unitless)</th><th>Mean absolute error (MAE, unitless)</th><th>Root-mean-square error (RMSE, unitless)</th></tr></thead><tbody><tr><td>Leave-one-year-out (2012–2022)</td><td>14.3<br>Q1–Q3 8; 20<br>min–max 3; 37</td><td>14.4<br>Q1–Q3 9; 20<br>min–max 2; 37</td><td>0.79<br>Q1–Q3 0.6; 0.8<br>min–max 0.2; 0.9</td><td>0.85<br>Q1–Q3 0.8; 0.9<br>min–max 0.3; 1.0</td><td>-0.02<br>Q1–Q3 -0.08; 0.05<br>min–max<br>-0.68; 0.70</td><td>0.10<br>Q1–Q3 0.06; 0.13<br>min–max<br>0.02; 0.32</td><td>0.07<br>Q1–Q3 0.05; 0.09<br>min–max 0.03; 0.15</td></tr><tr><td>Wet years (2013, 2014, 2016, 2018)</td><td>8.3<br>Q1–Q3 5; 15<br>min–max 1; 29</td><td>8.7<br>Q1–Q3 5; 13<br>min–max 0; 35</td><td>0.76<br>Q1–Q3 0.6; 0.9<br>min–max<br>-1.7; 1.0</td><td>0.71<br>Q1–Q3 0.6; 0.8<br>min–max<br>-4.0; 1.0</td><td>0.17<br>Q1–Q3 -1.23; 1.25<br>min–max<br>-7.09; 7.32</td><td>0.05<br>Q1–Q3 0.03; 0.09<br>min–max 0.01; 0.23</td><td>0.05<br>Q1–Q3 0.03; 0.08<br>min–max 0.02; 0.15</td></tr><tr><td>Intermediate years (2012, 2015, 2020, 2021)</td><td>12.9<br>Q1–Q3 7; 19<br>min–max 1; 35</td><td>13.6<br>Q1–Q3 7; 19<br>min–max 2; 31</td><td>0.76<br>Q1–Q3 0.6; 0.8<br>min–max<br>-0.7; 1.0</td><td>0.76<br>Q1–Q3 0.7; 0.9<br>min–max<br>-0.2; 1.0</td><td>0.21<br>Q1–Q3 -1.06; 1.41<br>min–max<br>-8.85; 5.76</td><td>0.07<br>Q1–Q3 0.04; 0.09<br>min–max 0; 0.25</td><td>0.06<br>Q1–Q3 0.04; 0.09<br>min–max 0; 0.18</td></tr><tr><td>Dry years (2017, 2019, 2022)</td><td>21.5<br>Q1–Q3 14; 30<br>min–max 2; 51</td><td>21.4<br>Q1–Q3 14; 32<br>min–max<br>4; 49</td><td>0.7<br>Q1–Q3 0.5; 0.8<br>min–max<br>-0.1; 1.0</td><td>0.75<br>Q1–Q3 0.6; 0.8<br>min–max<br>-0.6; 0.9</td><td>-0.71<br>Q1–Q3 -2.90; 1.49<br>min–max<br>-15.13; 12.78</td><td>0.09<br>Q1–Q3 0.06; 0.12<br>min–max 0.01; 0.43</td><td>0.08<br>Q1–Q3 0.06; 0.12<br>min–max 0.03; 0.24</td></tr></tbody></table>

Drying probability: percentage of ONDE sites in a dry state computed for each HER2, averaged by month; Q1–Q3: first and third quartiles; Min: minimum; Max: maximum.
Appendix I: Model performance using discharge data simulated with SAFRAN

<figure><img src="image_22.png" /><figcaption>**Figure 11.** Nash–Sutcliffe model efficiency (NSE) coefficient (line 1), bias (line 2), mean absolute error (MAE, line 3), and root-mean-square error (RMSE, line 4) assessing the calibration of logistic regression using discharge data simulated with SAFRAN data available during the calibration period (2012–2022).</figcaption></figure>
Appendix J: Application to 21st century PFI modelling

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

**Figure J1.** Ensemble median mPFI<sub>7–10</sub> (columns 1 to 3) and ensemble median shift of the start date Tf (columns 4 and 5) and the end date TI (columns 6 and 7) of the dry period over the two periods 2041–2070 and 2070–2099 under the three RCPs for the CTRIP hydrological model, relative to the baseline period 1976–2005. The shift, expressed in weeks, takes a positive value when the duration of the dry period increases. Grey HER2s had no period with a PFI > 20 % during the reference period.
<figure><img src="image_24.png" /></figure>

**Figure J2.** Ensemble median mPFI$_{7-10}$ (columns 1 to 3) and ensemble median shift of the start date Tf (columns 4 and 5) and the end date TI (columns 6 and 7) of the dry period over the two periods 2041–2070 and 2070–2099 under the three RCPs for the J2000 hydrological model, relative to the baseline period 1976–2005. The shift, expressed in weeks, takes a positive value when the duration of the dry period increases. Grey HER2s had no period with a PFI > 20 % during the reference period.
<figure><img src="image_25.png" /></figure>

**Figure J3.** Ensemble median mPFI<sub>7–10</sub> (columns 1 to 3) and ensemble median shift of the start date Tf (columns 4 and 5) and the end date Tl (columns 6 and 7) of the dry period over the two periods 2041–2070 and 2070–2099 under the three RCPs for the ORCHIDEE hydrological model, relative to the baseline period 1976–2005. The shift, expressed in weeks, takes a positive value when the duration of the dry period increases. Grey HER2s had no period with a PFI > 20 % during the reference period.
<figure><img src="image_26.png" /><figcaption>**Figure J4.** Ensemble median mPFI<sub>7-10</sub> (columns 1 to 3) and ensemble median shift of the start date Tf (columns 4 and 5) and the end date TI (columns 6 and 7) of the dry period over the two periods 2041–2070 and 2070–2099 under the three RCPs for the SMASH hydrological model, relative to the baseline period 1976–2005. The shift, expressed in weeks, takes a positive value when the duration of the dry period increases. Grey HER2s had no period with a PFI > 20 % during the reference period.</figcaption></figure>

## J1 Narratives

The following figures present time series of the evolution of PFI dynamics at different horizons and under different RCP scenarios for each hydrological model in HER2-12, HER2-57, HER2-81, and HER2-105. Four contrasting scenarios were highlighted among the Explore2 climate projections to illustrate a diversity of potential changes under RCP8.5. These story lines range from “Strong warming and strong summer (and annual) drying”, selected as one of the most extreme climate projections, to “Moderate warming and precipitation change” with less pronounced alterations, with two alternative projections named “Dry all year, reduced winter recharge” and “Hot and humid all seasons”. They are also illustrated here to present distinct hydrological nuances.
Average proportion of zero flows smoothed over 5 d
Hydrological model: CTRIP - HER2 12

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

**Figure J5.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the CTRIP hydrological model in HER2-12 (d = day).

Average proportion of zero flows smoothed over 5 d
Hydrological model: GRSD - HER2 12

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

**Figure J6.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the GRSD hydrological model in HER2-12.
Average proportion of zero flows smoothed over 5 d
Hydrological model: ORCHIDEE - HER2 12

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

**Figure J7.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the ORCHIDEE hydrological model in HER2-12.

Average proportion of zero flows smoothed over 5 d
Hydrological model: SMASH - HER2 12

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

**Figure J8.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the SMASH hydrological model in HER2-12.
Average proportion of zero flows smoothed over 5 d
Hydrological model: CTRIP - HER2 57

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

**Figure J9.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the CTRIP hydrological model in HER2-57.

Average proportion of zero flows smoothed over 5 d
Hydrological model: GRSD - HER2 57

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

**Figure J10.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the GRSD hydrological model in HER2-57.
Average proportion of zero flows smoothed over 5 d
Hydrological model: ORCHIDEE - HER2 57

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

**Figure J11.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the ORCHIDEE hydrological model in HER2-57.

Average proportion of zero flows smoothed over 5 d
Hydrological model: SMASH - HER2 57

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

**Figure J12.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the SMASH hydrological model in HER2-57.
Average proportion of zero flows smoothed over 5 d
Hydrological model: CTRIP - HER2 81

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

**Figure J13.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the CTRIP hydrological model in HER2-81.

Average proportion of zero flows smoothed over 5 d
Hydrological model: GRSD - HER2 81

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

**Figure J14.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the GRSD hydrological model in HER2-81.
Average proportion of zero flows smoothed over 5 d
Hydrological model: J2000 - HER2 81

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

**Figure J15.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the J2000 hydrological model in HER2-81.

Average proportion of zero flows smoothed over 5 d
Hydrological model: ORCHIDEE - HER2 81

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

**Figure J16.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the ORCHIDEE hydrological model in HER2-81.
Average proportion of zero flows smoothed over 5 d
Hydrological model: SMASH - HER2 81

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

**Figure J17.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the SMASH hydrological model in HER2-81.

Average proportion of zero flows smoothed over 5 d
Hydrological model: CTRIP - HER2 105

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

**Figure J18.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the CTRIP hydrological model in HER2-105.
Average proportion of zero flows smoothed over 5 d
Hydrological model: GRSD - HER2 105

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

**Figure J19.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the GRSD hydrological model in HER2-105.

Average proportion of zero flows smoothed over 5 d
Hydrological model: ORCHIDEE - HER2 105

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

**Figure J20.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the ORCHIDEE hydrological model in HER2-105.
Average proportion of zero flows smoothed over 5 d
Hydrological model: SMASH - HER2 105

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

**Figure J21.** Time series of the evolution of PFI dynamics at different horizons under different RCP scenarios for the SMASH hydrological model in HER2-105.

**Appendix K: Agreement between changes under RCP4.5**

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

**Figure K1.** Agreement between projections of mPFI<sub>7-10</sub> for each hydrological model and inter-model agreement on the change signal of mPFI<sub>7-10</sub> under the RCP4.5 scenario.
*Code availability.* The codes are available at https://github.com/tjaouen/WillRiversBecomeMoreIntermittentInFrance (last access: 17 October 2024; DOI: https://doi.org/10.5281/zenodo.15974001, Jaouen, 2025).

*Data availability.* The Explore2 streamflow projections are available at https://entrepot.recherche.data.gouv.fr/dataverse/explore2-projections_hydrologiques (Sauquet et al., 2025b) and are described in detail by Sauquet et al. (2024) (https://doi.org/10.57745/J3XIPW). In addition, associated data on projected changes and related uncertainties are available at https://doi.org/10.57745/KWH320 (Evin et al., 2024). The ONDE dataset is available at https://hubeau.eaufrance.fr/page/api-ecoulement (OFB, 2025).

*Author contributions.* TJ, LB, and ES designed the experiments, and TJ carried them out. LH organized the data from the Explore2 project and made them available. TJ and ES developed the model code. TJ performed the simulations. TJ prepared the manuscript with contributions from all co-authors.

*Competing interests.* The contact author has declared that none of the authors has any competing interests.

*Disclaimer.* Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

*Special issue statement.* This article is part of the special issue "Drought, society, and ecosystems (NHESS/BG/GC/HESS inter-journal SI)". It is not associated with a conference.

*Acknowledgements.* Scientific oversight of this research was provided within INRAE by Éric Sauquet and Lionel Benoît. This study builds on exploratory work by Aurélien Beaufort and Éric Sauquet, with contributions from all participants in the Explore2 project. The archiving and management of Explore2 data, expertly handled by Louis Héraut, also played a crucial role in this study. Additionally, this research directly applies the QUALYPSO method, developed and adapted to Explore2 by Guillaume Evin. We extend our deep gratitude to all these individuals for their invaluable contributions to this work.

*Financial support.* This research has been supported by the Rhône-Méditerranée-Corse Water Agency under the supervision of Benoît Terrier.

*Review statement.* This paper was edited by Elena Toth and reviewed by Nicola Durighetto and two anonymous referees.

## References

Acuña, V., Datry, T., Marshall, J., Barceló, D., Dahm, C. N., Giné-breda, A., McGregor, G., Sabater, S., Tockner, K., and Palmer, M. A.: Why should we care about temporary waterways?, Science, 343, 1080–1081, https://doi.org/10.1126/science.1246666, 2014.

Addor, N., Nearing, G., Prieto, C., Newman, A. J., le Vine, N., and Clark, M. P.: A Ranking of Hydrological Signatures Based on Their Predictability in Space, Water Resour. Res., 54, 8792–8812, https://doi.org/10.1029/2018WR022606, 2018.

Aitken, G., Beevers, L., Parry, S., and Facer-Childs, K.: Partitioning model uncertainty in multi-model ensemble river flow projections, Climatic Change, 176, 153, https://doi.org/10.1007/s10584-023-03621-1, 2023.

Barrow, C. J.: World atlas of desertification (United Nations Environment Programme), edited by N. Middleton and D. S. G. Thomas. Edward Arnold, London, 1992, ISBN 0 340 55512 2, £89.50 (hardback), ix + 69 pp., Land Degrad. Dev., 3, 249–249, https://doi.org/10.1002/ldr.3400030407, 1992.

Beaufort, A., Lamouroux, N., Pella, H., Datry, T., and Sauquet, E.: Extrapolating regional probability of drying of headwater streams using discrete observations and gauging networks, Hydrol. Earth Syst. Sci., 22, 3033–3051, https://doi.org/10.5194/hess-22-3033-2018, 2018.

Beaufort, A., Carreau, J., and Sauquet, E.: A classification approach to reconstruct local daily drying dynamics at headwater streams, Hydrol. Process., 33, 1896–1912, https://doi.org/10.1002/hyp.13445, 2019.

Belemtougri, A. P., Ducharne, A., Tazen, F., Oudin, L., and Karambiri, H.: Understanding key factors controlling the duration of river flow intermittency: Case of Burkina Faso in West Africa, Journal of Hydrology: Regional Studies, 37, 100908, https://doi.org/10.1016/j.ejrh.2021.100908, 2021.

Bertassello, L. E., Durighetto, N., and Botter, G.: Eco-hydrological modelling of channel network dynamics - Part 2: Application to metapopulation dynamics, Roy. Soc. Open Sci., 9, 220945, https://doi.org/10.1098/rsos.220945, 2022.

Botter, G. and Durighetto, N.: The Stream Length Duration Curve: A Tool for Characterizing the Time Variability of the Flowing Stream Length, Water Resour. Res., 56, e2020WR027282, https://doi.org/10.1029/2020WR027282, 2020.

Bourke, S. A., Degens, B., Searle, J., de Castro Tayer, T., and Rothery, J.: Geological permeability controls streamflow generation in a remote, ungauged, semi-arid drainage system, Journal of Hydrology: Regional Studies, 38, 100956, https://doi.org/10.1016/j.ejrh.2021.100956, 2021.

Carlson, S. M., Ruhí, A., Bogan, M. T., Hazard, C. W., Ayers, J., Grantham, T. E., Batalla, R. J., and Garcia, C.: Losing flow in free-flowing Mediterranean-climate streams, Front. Ecol. Environment., 22, e2737, https://doi.org/10.1002/fee.2737, 2024.

Dai, Z., Amatya, D. M., Sun, G., Trettin, C. C., Li, C., and Li, H.: Climate variability and its impact on forest hydrology on south Carolina coastal plain, USA, Atmosphere, 2, 330–357, https://doi.org/10.3390/atmos2030330, 2011.
Dakhlaoui, H., Ruelland, D., Tramblay, Y., and Bargaoui, Z.: Evaluating the robustness of conceptual rainfall-runoff models under climate variability in northern Tunisia, J. Hydrol., 550, 201–217, https://doi.org/10.1016/j.jhydrol.2017.04.032, 2017.

Decharme, B., Delire, C., Minvielle, M., Colin, J., Vergnes, J.-P., Alias, A., Saint-Martin, D., Séférion, R., Sénési, S., and Voldoire, A.: Recent changes in the ISBA-CTRIP land surface system for use in the CNRM-CM6 climate model and in global off-line hydrological applications, J. Adv. Model. Earth Sy., 11, 1207–1252, https://doi.org/10.1029/2018MS001545, 2019.

De Girolamo, A. M., Barca, E., Leone, M., and Lo Porto, A.: Impact of long-term climate change on flow regime in a Mediterranean basin, Journal of Hydrology: Regional Studies, 41, 101061, https://doi.org/10.1016/j.jehr.2022.101061, 2022.

de Lavenne, A., Thirel, G., Andréassian, V., Perrin, C., and Ramos, M.-H.: Spatial variability of the parameters of a semi-distributed hydrological model, Proc. IAHS, 373, 87–94, https://doi.org/10.5194/piahs-373-87-2016, 2016.

Delso, J., Magdaleno, F., and Fernández-Yuste, J. A.: Flow patterns in temporary rivers: a methodological approach applied to southern Iberia, Hydrol. Sci. J., 62, 1551–1563, https://doi.org/10.1080/02626667.2017.1346375, 2017.

Dhungel, S., Tarboton, D. G., Jin, J., and Hawkins, C. P.: Potential effects of climate change on ecologically relevant streamflow regimes, River Res. Applic., 32, 1827–1840, https://doi.org/10.1002/rra.3029, 2016.

Döll, P. and Schmied, H. M.: How is the impact of climate change on river flow regimes related to the impact on mean annual runoff? A global-scale analysis, Environ. Res. Lett., 7, 014037, https://doi.org/10.1088/1748-9326/7/1/014037, 2012.

Döll, P. and Zhang, J.: Impact of climate change on freshwater ecosystems: a global-scale analysis of ecologically relevant river flow alterations, Hydrol. Earth Syst. Sci., 14, 783–799, https://doi.org/10.5194/hess-14-783-2010, 2010.

Döll, P., Abbasi, M., Messenger, M. L., Trautmann, T., Lehner, B., Lamouroux, N.: Streamflow Intermittence in Europe: Estimating High-Resolution Monthly Time Series by Downscaling of Simulated Runoff and Random Forest Modeling, Water Resour. Res., 60, e2023WR036900, https://doi.org/10.1029/2023WR036900, 2024.

Douville, H., Raghavan, K., Renwick, J., Allan, R. P., Arias, P. A., Barlow, M., Cerezo-Mota, R., Cherchi, A., Gan, T.Y., Gergis, J., Jiang, D., Khan, A., Pokam Mba, W., Rosenfeld, D., Tierney, J., and Zolina, O.: Water Cycle Changes. In Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change [Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S.L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M.I., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J.B.R., Maycock, T.K., Waterfield, T., Yelekçi, O., Yu, R., and Zhou, B.], Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 1055–1210, https://doi.org/10.1017/9781009157896.010, 2021.

Durighetto, N., Bertassello, L. E., and Botter, G.: Eco-hydrological modelling of channel network dynamics - Part 1: Stochastic simulation of active stream expansion and retraction, Roy. Soc. Open Sci., 9, 220944, https://doi.org/10.1098/rsos.220944, 2022.

Evin, G., Hingray, B., Blanchet, J., Eckert, N., Morin, S., and Verfaillie, D.: Partitioning uncertainty components of an incomplete

ensemble of climate projections using data augmentation, J. Climate, 32, 2423–2440, https://doi.org/10.1175/jcli-d-18-0606.1, 2019.

Evin, G., Somot, S., and Hingray, B.: Balanced estimate and uncertainty assessment of European climate change using the large EURO-CORDEX regional climate model ensemble, Earth Syst. Dynam., 12, 1543–1569, https://doi.org/10.5194/esd-12-1543-2021, 2021.

Evin, G., Hingray, B., Reverdy, A., Ducharne, A., and Sauquet, E.: Ensemble de projections Explore2: Changements moyens et incertitudes associées, Recherche Data Gouv [data set], https://doi.org/10.57745/KWH320, 2024.

Finn, D., Bonada, N., Múrria, C., and Hughes, J.: Small but mighty: headwaters are vital to stream network biodiversity at two levels of organization, J. N. Am. Benthol. Soc., 30, 963–980, https://doi.org/10.1899/11-012.1, 2011.

Gallart, F., Prat, N., García-Roger, E. M., Latron, J., Rieradevall, M., Llorens, P., Barberá, G. G., Brito, D., De Girolamo, A. M., Lo Porto, A., Buffagni, A., Erba, S., Neves, R., Nikolaidis, N. P., Perrin, J. L., Querner, E. P., Quiñonero, J. M., Tournoud, M. G., Tzoraki, O., Skoulikidis, N., Gómez, R., Sánchez-Montoya, M. M., and Froebrich, J.: A novel approach to analysing the regimes of temporary streams in relation to their controls on the composition and structure of aquatic biota, Hydrol. Earth Syst. Sci., 16, 3165–3182, https://doi.org/10.5194/hess-16-3165-2012, 2012.

Geris, J., Tetzlaff, D., and Soulsby, C.: Resistance and resilience to droughts: Hydropedological controls on catchment storage and run-off response, Hydrol. Process., 29, 4579–4593, https://doi.org/10.1002/hyp.10480, 2015.

Giezendanner, J., Benettin, P., Durighetto, N., Botter, G., and Rinaldo, A.: A Note on the Role of Seasonal Expansions and Contractions of the Flowing Fluvial Network on Metapopulation Persistence, Water Resour. Res., 57, e2021WR029813, https://doi.org/10.1029/2021WR029813, 2021.

González-Ferreras, A. M., and Barquín, J.: Mapping the temporary and perennial character of whole river networks, Water Resour. Res., 53, 6709–6724, https://doi.org/10.1002/2017WR020390, 2017.

Gudmundsson, L., Boulange, J., Do, H. X., Gosling, S. N., Grillakis, M. G., Koutroulis, A. G., Leonard, M., Liu, J., Müller Schmied, H., Papadimitriou, L., Pokhrel, Y., Seneviratne, S. I., Satoh, Y., Thiery, W., Westra, S., Zhang, X., and Zhao, F.: Globally observed trends in mean and extreme river flow attributed to climate change, Science, 371, 1159–1162, https://doi.org/10.1126/science.aba3996, 2021.

Gupta, H., Kling, H., Yilmaz, K., and Martinez, G.: Decomposition of the mean squared error and NSE performance criteria: Implications for improving hydrological modelling, J. Hydrol., 377, 80–91, https://doi.org/10.1016/j.jhydrol.2009.08.003, 2009.

Halloran, L., Millwater, J., Hunkeler, D., and Arnoux, M.: Climate change impacts on groundwater discharge-dependent streamflow in an alpine headwater catchment, Sci. Total Environ., 902, 166009, https://doi.org/10.1016/j.scitotenv.2023.166009, 2023.

Hammond, J. C. and Fleming, B. J.: Evaluating low flow patterns, drivers and trends in the Delaware River Basin, J. Hydrol., 598, 126246, https://doi.org/10.1016/j.jhydrol.2021.126246, 2021.

Huang, P., Ducharne, A., Rinchiuso, L., Polcher, J., Baratgin, L., Bastrikov, V., and Sauquet, E.: Multi-objective calibration and evaluation of the ORCHIDEE land surface model over France
at high resolution, Hydrol. Earth Syst. Sci., 28, 4455–4476, https://doi.org/10.5194/hess-28-4455-2024, 2024.

Jaeger, K. L., Olden, J. D., and Pelland, N. A.: Climate change poised to threaten hydrologic connectivity and endemic fishes in dryland streams, P. Natl. Acad. Sci. USA, 111, 13894–13899, https://doi.org/10.1073/pnas.1320890111, 2014.

Jaeger, K. L., Sando, R., McShane, R. R., Dunham, J. B., Hockman- Wert, D. P., Kaiser, K. E., Hafen, K., Risley, J. C., and Blasch, K. W.: Probability of Streamflow Permanence Model (PROSPER): A spatially continuous model of annual streamflow permanence throughout the Pacific Northwest, Journal of Hydrology X, 2, 100005, https://doi.org/10.1016/j.hydroa.2018.100005, 2019.

Jaouen, T.: tjaouen/WillRiversBecomeMoreIntermittentInFrance: WillRiversBecomeMoreIntermittentInFrance_Article, Zenodo [code], https://doi.org/10.5281/zenodo.15974001, 2025.

Jay-Allemand, M., Javelle, P., Gejadze, I., Arnaud, P., Malaterre, P.- O., Fine, J.-A., and Organde, D.: On the potential of variational calibration for a fully distributed hydrological model: applica- tion on a Mediterranean catchment, Hydrol. Earth Syst. Sci., 24, 5519–5538, https://doi.org/10.5194/hess-24-5519-2020, 2020.

Lapides, D. A., Leclerc, C. D., Moidu, H., Dralle, D. N., and Hahm, W. J.: Variability of stream extents controlled by flow regime and network hydraulic scaling, Hydrol. Process., 35, e14079, https://doi.org/10.1002/hyp.14079, 2021.

Larned, S. T., Datry, T., Arscott, D. B., and Tockner, K.: Emerging concepts in temporary-river ecology, Freshwater Biol., 55, 717– 738, https://doi.org/10.1111/j.1365-2427.2009.02322.x, 2010.

Leigh, C. and Datry, T.: Drying as a primary hydrological deter- minant of biodiversity in river systems: A broad-scale analysis, Ecography, 40, 487–499, https://doi.org/10.1111/ecog.02230, 2017.

Leleu, I., Tonnelier, I., Puechberty, R., Gouin, P., Viquendi, I., Cobos, L., Foray, A., Baillon, M., and Ndima, P.-O.: Re- founding the national information system designed to manage and give access to hydrometric data, Houille Blanche, 100, 25– 32, https://doi.org/10.1051/lhb/2014004, 2014.

Message, M. L., Lehner, B., Cockburn, C., Lamouroux, N., Pella, H., Snelder, T., Tockner, K., Trautmann, T., Watt, C., and Datry, T.: Global prevalence of non-perennial rivers and streams, Na- ture, 594, 391–397, https://doi.org/10.1038/s41586-021-03565- 5, 2021.

Meyer, J., Strayer, D., Wallace, J., Eggert, S., Helfman, G., and Leonard, N.: The contribution of headwater streams to biodiver- sity in river networks, J. Am. Water Resour. As., 43, 86–103, https://doi.org/10.1111/j.1752-1688.2007.00008.x, 2007.

Morel, M., Pella, H., Branger, F., Sauquet, E., Grenouillet, G., Côte, J., Braud, I., and Lamouroux, N.: Catchment-scale applications of hydraulic habitat models: Climate change effects on fish, Eco- hydrology, 16, e2513, https://doi.org/10.1002/eco.2513, 2023.

Nash, J. E. and Sutcliffe, J. V.: River flow forecasting through con- ceptual models part I – A discussion of principles, J. Hydrol., 10, 282–290, https://doi.org/10.1016/0022-1694(70)90255-6, 1970.

OFB: Observatoire National Des Etiages (ONDE), Office Français de la Biodiversité [data set], https://hubeau.eaufrance.fr/page/ api-ecoulement, last access: 4 March 2025.

Pumo, D., Caracciolo, D., Viola, F., and Noto, L. V.: Cli- mate change effects on the hydrological regime of small non-perennial river basins, Sci. Total Environ., 542, 76–92, https://doi.org/10.1016/j.scitotenv.2015.10.109, 2016.

Reidy Liermann, C., Olden, J., Beechie, T., Kennard, M., Skidmore, P., Konrad, C., and Imaki, H.: Hydrogeomorphic classification of Washington State rivers to support emerging environmental flow management strategies, River Res. Appl., 28, 1340–1358, https://doi.org/10.1002/rra.1541, 2012.

Reynolds, L. V., Shafroth, P. B., and Poff, L.: Modeled in- termittency risk for small streams in the Upper Colorado River Basin under climate change, J. Hydrol., 523, 768–780, https://doi.org/10.1016/j.jhydrol.2015.02.025, 2015.

Ribes, A., Boé, J., Qasmi, S., Dubuisson, B., Douville, H., and Terray, L.: An updated assessment of past and future warm- ing over France based on a regional observational constraint, Earth Syst. Dynam., 13, 1397–1415, https://doi.org/10.5194/esd- 13-1397-2022, 2022.

Rutkowska, A., Osuch, M., Żelazny, M., Banasik, K., and Klimek, M.: Climatic and anthropogenic drivers of zero-flow events in intermittent rivers in Poland, Journal of Water and Land Devel- opment, 57, 52–61, https://doi.org/10.24425/jwld.2023.145335, 2023.

Sando, R., Jaeger, K. L., Farmer, W. H., Barnhart, T. B., McShane, R. R., Welborn, T. L., Kaiser, K. E., Hafen, K. C., Blasch, K., York, B., and Shallcross, A.: Predictions and drivers of sub- reach-scale annual streamflow permanence for the upper Mis- souri River basin: 1989–2018, Journal of Hydrology X, 17, 100138, https://doi.org/10.1016/j.hydroa.2022.100138, 2022.

Sarremejane, R., Cañedo-Argüelles, M., Prat, N., Mykrä, H., Muotka, T., and Bonada, N.: Do metacommunities vary through time? Intermittent rivers as model systems, J. Biogeogr., 44, 2752–2763, https://doi.org/10.1111/jbi.13077, 2017.

Sauquet, E., Beaufort, A., Sarremejane, R., and Thirel, G.: Predicting flow intermittence in France under climate change, Hydrolog. Sci. J., 66, 2046-2059, https://doi.org/10.1080/02626667.2021.1963444, 2021.

Sauquet, E., Evin, G., Siauve, S., Bornancin-Plantier, A., Jacquin, N., Arnaud, P., Bérel, M., Bernus, S., Bonneau, J., Branger, F., Caballero, Y., Colléoni, F., Collet, L., Corre, L., Drouin, A., Ducharne, A., Fournier, M., Gailhard, J., Habets, F., Hendrickx, F., Héraut, L., Hingray, B., Huang, P., Jaouen, T., Jeantet, A., Lanini, S., Le Lay, M., Loudin, S., Magand, C., Marson, P., Mimeau, L., Monteil, C., Munier, S., Perrin, C., Robin, Y., Rous- set, F., Soubeyroux, J.-M., Strohmenger, L., Thirel, G., Toc- quer, F., Tramblay, Y., Vergnes, J.-P., Vidal, J.-P., and Vrac, M.: Messages et enseignements du projet Explore2, Recherche Data Gouv [data set], https://doi.org/10.57745/J3XIPW, 2024.

Sauquet, E., Evin, G., Siauve, S., Aissat, R., Arnaud, P., Bérel, M., Bonneau, J., Branger, F., Caballero, Y., Colléoni, F., Ducharne, A., Gailhard, J., Habets, F., Hendrickx, F., Héraut, L., Hingray, B., Huang, P., Jaouen, T., Jeantet, A., Lanini, S., Le Lay, M., Magand, C., Mimeau, L., Monteil, C., Munier, S., Perrin, C., Robelin, O., Rousset, F., Soubeyroux, J.-M., Strohmenger, L., Thirel, G., Tocquer, F., Tramblay, Y., Vergnes, J.-P., and Vidal, J.-P.: A large transient multi-scenario multi-model ensemble of future streamflow and groundwater projections in France, EGU- sphere [preprint], https://doi.org/10.5194/egusphere-2025-1788, 2025a.

Sauquet, E., Evin, G., Siauve, S., Aissat, R., Arnaud, P., Bérel, M., Bonneau, J., Branger, F., Caballero, Y., Colléoni, F., Ducharne, A., Gailhard, J., Garambois, P., Habets, F., Hen- drickx, F., Héraut, L., Hingray, B., Huang, P., Jaouen, T.,
Jeantet, A., Lanini, S., Le Lay, M., Magand, C., Mimeau, L., Monteil, C., Munier, S., Perrin, C., Robelin, O., Rousset, F., Soubeyroux, J., Strohmenger, L., Thirel, G., Tocquer, F., Tramblay, Y., Vergnes, J., Vidal, J.: Ensemble des projections hydrologiques Explore2 du débit, INRAE & Entrepôt Recherche Data Gouv [data set], https://entrepot.recherche.data.gouv.fr/dataverse/explore2-projections_hydrologiques, last access: 4 March 2025b.

Scheller, M., Van Meerveld, I., and Seibert, J.: How well can people observe the flow state of temporary streams? Frontiers in Environmental Science, 12, 1352697, https://doi.org/10.3389/fenvs.2024.1352697, 2024.

Schneider, A., Jost, A., Coulon, C., Silvestre, M., Théry, S., and Ducharne, A.: Global-scale river network extraction based on high-resolution topography and constrained by lithology, climate, slope, and observed drainage density, Geophys. Res. Lett., 44, 2773–2781, https://doi.org/10.1002/2016GL071844, 2017.

Schneider, C., Laizé, C. L. R., Acreman, M. C., and Flörke, M.: How will climate change modify river flow regimes in Europe?, Hydrol. Earth Syst. Sci., 17, 325–339, https://doi.org/10.5194/hess-17-325-2013, 2013.

Sefton, C. E. M., Parry, S., England, J., and Angell, G.: Visualising and quantifying the variability of hydrological state in intermittent rivers, Fund. Appl. Limnol., 193, 21–38, https://doi.org/10.1127/fal/2019/1149, 2019.

Shaw, S. B., Bonville, D. B., and Chandler, D. G.: Combining observations of channel network contraction and spatial discharge variation to inform spatial controls on baseflow in Birch Creek, Catskill Mountains, USA, Journal of Hydrology: Regional Studies, 12, 1–12, https://doi.org/10.1016/j.ejrh.2017.03.003, 2017.

Snelder, T. H., Datry, T., Lamouroux, N., Larned, S. T., Sauquet, E., Pella, H., and Catalogne, C.: Regionalization of patterns of flow intermittence from gauging station records, Hydrol. Earth Syst. Sci., 17, 2685–2699, https://doi.org/10.5194/hess-17-2685-2013, 2013.

Strohmenger, L., Sauquet, E., Bernard, C., Bonneau, J., Branger, F., Bresson, A., Brigode, P., Buzier, R., Delaigue, O., Devers, A., Evin, G., Fournier, M., Hsu, S.-C., Lanini, S., de Lavenne, A., Lemaitre-Basset, T., Magand, C., Mendoza Guimarães, G., Mentha, M., Munier, S., Perrin, C., Podechard, T., Rouchy, L., Sadki, M., Soutif-Bellenger, M., Tilmant, F., Tramblay, Y., Véron, A.-L., Vidal, J.-P., and Thirel, G.: On the visual detection of non-natural records in streamflow time series: challenges and impacts, Hydrol. Earth Syst. Sci., 27, 3375–3391, https://doi.org/10.5194/hess-27-3375-2023, 2023.

Tramblay, Y. and Somot, S.: Future evolution of extreme precipitation in the Mediterranean, Climatic Change, 151, 289–302, https://doi.org/10.1007/s10584-018-2300-5, 2018.

Tramblay, Y., Rutkowska, A., Sauquet, E., Sefton, C., Laaha, G., Osuch, M., Albuquerque, M. T. D., Alves, M. H., Banasik, K., Beaufort, A., Brocca, L., Camici, S., Csabai, Z., Dakhlaoui, H., DeGirolamo, A. M., Dörflinger, G., Gallart, F., Gauster, T., Hanich, L., Kohnová, S., Mediero, L., Plamen, N., Parry, S., Quintana-Seguí, P., Tzoraki, O., and Datry, T.: Trends in flow intermittence for European rivers, Hydrolog. Sci. J., 66, 37–49, https://doi.org/10.1080/02626667.2020.1849708, 2020.

Van Meerveld, H. J. I., Sauquet, E., Gallart, F., Sefton, C., Seibert, J., and Bishop, K.: Aqua temporaria incognita, Hydrol. Process., 34, 5704–5711, https://doi.org/10.1002/hyp.13979, 2020.

Verfaillie, D., Lafaysse, M., Déqué, M., Eckert, N., Lejeune, Y., and Morin, S.: Multi-component ensembles of future meteorological and natural snow conditions for 1500 m altitude in the Chartreuse mountain range, Northern French Alps, The Cryosphere, 12, 1249–1271, https://doi.org/10.5194/tc-12-1249-2018, 2018.

Vicente-Serrano, S. M., Peña-Gallardo, M., Hannaford, J., Murphy, C., Lorenzo-Lacruz, J., Dominguez-Castro, F., López-Moreno, J. I., Beguería, S., Noguera, I., Harrigan, S., and Vidal, J. P.: Climate, Irrigation, and Land Cover Change Explain Streamflow Trends in Countries Bordering the Northeast Atlantic, Geophys. Res. Lett., 46, 10821–10833, https://doi.org/10.1029/2019GL084084, 2019.

Vicente-Serrano, S. M., McVicar, T. R., Miralles, D. G., Yang, Y., and Tomas-Burguera, M.: Unravelling the influence of atmospheric evaporative demand on drought and its response to climate change, WIREs Clim Change, 11, e632, https://doi.org/10.1002/wcc.632, 2020.

Vidal, J.-P., Martin, E., Franchistéguy, L., Baillon, M., and Soubeyroux, J.-M.: A 50-year high-resolution atmospheric reanalysis over France with the Safran system, Int. J. Climatol., 30, 1627–1644, https://doi.org/10.1002/joc.2003, 2010.

Wasson, J. G., Chandesris, A., Pella, H., and Blanc, L.: Définition des hydro-écorégions françaises métropolitaines. Approche régionale de la typologie des eaux courantes et éléments pour la définition des peuplements de référence d'invertébrés, Irstea, hal-02580774, 190 pp., https://hal.inrae.fr/hal-02580774v1 (last access: 6 September 2024), 2002.

Yu, X. and Duffy, C. J.: Watershed hydrology: Scientific advances and environmental assessments, Water, 10, 288, https://doi.org/10.3390/w10030288, 2018.

Zipper, S., Hammond, J. C., Shanafield, M., Zimmer, M., Datry, T., Jones, C. N., Kaiser, K. E., Godsey, S. E., Burrows, R. M., Blaszczak, J. R., Busch, M. H., Price, A. N., Boersma, K. S., Ward, A. S., Costigan, K., Allen, G. H., Krabbenhoft, C. A., Dodds, W. K., Mims, M. C., Olden, J. D., Kampf, S. K., Burgin, A. J., and Allen, D. C.: Pervasive changes in stream intermittency across the United States, Environ. Res. Lett., 16, 084033, https://doi.org/10.1088/1748-9326/ac14ec, 2021.