# Self-potential signals related to tree transpiration in a Mediterranean climate

Kaiyan Hu, Bertille Loiseau, Simon D Carrière, Nolwenn Lesparre, Cédric Champollion, Nicolas K Martin-Stpaul, Niklas Linde, Damien Jougnot

► To cite this version:

Kaiyan Hu, Bertille Loiseau, Simon D Carrière, Nolwenn Lesparre, Cédric Champollion, et al. Self-potential signals related to tree transpiration in a Mediterranean climate. Hydrology and Earth System Sciences, 2025, 29, pp.2997-3018. 10.5194/hess-29-2997-2025. hal-04727524v2

HAL Id: hal-04727524

https://hal.sorbonne-universite.fr/hal-04727524v2

Submitted on 26 Sep 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>

# Self-potential signals related to tree transpiration in a Mediterranean climate

Kaiyan Hu<fnref n="1" />, Bertille Loiseau<fnref n="2" />,<fnref n="7" />, Simon D. Carrière<fnref n="2" />,<fnref n="7" />, Nolwenn Lesparre<fnref n="3" />, Cédric Champollion<fnref n="4" />, Nicolas K. Martin-StPaul<fnref n="5" />, Niklas Linde<fnref n="6" />, and Damien Jougnot<fnref n="2" />

**Correspondence:** Damien Jougnot (damien.jougnot@upmc.fr)

Received: 3 August 2024 – Discussion started: 23 September 2024

Revised: 1 March 2025 – Accepted: 1 April 2025 – Published: 16 July 2025

**Abstract.** Plant transpiration is a crucial process in the water cycle, and its quantification is essential for understanding terrestrial ecosystem dynamics. While sap flow measurements offer a direct method for estimating individual tree transpiration, their effectiveness may be limited by the use of point sensors; species-specific calibration requirements; and baseline uncertainties, particularly the assumption of negligible nighttime flow, which may not always hold. Self-potential (SP), a passive geophysical method, holds potential for constraining transpiration rates, though many questions remain regarding the electrophysiological processes occurring within trees. In this study, we continuously measured tree SP and sap velocity on three tree species for 1 year in a Mediterranean climate. Using wavelet coherence analysis and variational mode decomposition, we explored the empirical relationship between tree SP and transpiration. Our analysis revealed strong coherence between SP and sap velocity at diurnal timescales, with coherence weakening and phase shifts increasing on days with higher water supply. We estimated electrokinetic coupling coefficients using a linear regression model between SP and sap velocity variations at the diurnal scale, resulting in values typically found in porous geological media. During dry seasons, the electrokinetic effect emerges as the primary contribution to tree SP, indicating its potential utility in assessing transpiration rates. Our

results emphasize the need for improved electrode configurations and physiochemical modelling to elucidate tree SP in relation to transpiration.

## 1 Introduction

Climate change and anthropogenic activities exert significant influence on terrestrial ecosystems, particularly impacting vegetation (Mottl et al., 2021; Nolan et al., 2018). In anticipation of these growing impacts, understanding the effects of phenomena such as droughts on forests is essential (Anderegg et al., 2016). Evapotranspiration (ET) stands as a critical component of the water cycle, accounting for approximately 65% of precipitated water at the continental scale (Fisher et al., 2017; Oki and Kanae, 2006). Transpiration, constituting 60%-80% of global terrestrial evapotranspiration (Bachofen et al., 2023; Dix and Aubrey, 2021; Jasechko et al., 2013; Kuang et al., 2024; Schlesinger and Jasechko, 2014; Wei et al., 2017), can reach 95% in tropical forests during dry seasons (e.g. Kunert et al., 2017). Monitoring of transpiration is crucial for enhancing our understanding of plant-water relations, ecohydrological processes, and water resource management, with implications for agriculture, ecology, and climate science. Plant transpiration can be measured

<footnote n="1">¹Hubei Subsurface Multiscale Image Key Laboratory, School of Geophysics and Geomatics, China University of Geosciences, Wuhan 430074, China</footnote>

<footnote n="2">²Sorbonne Université, CNRS, EPHE, UMR 7619 METIS, 75005 Paris, France</footnote>

<footnote n="3">³Université de Strasbourg, CNRS, EOST, ENGEES, ITES UMR 7063, 67000 Strasbourg, France</footnote>

<footnote n="4">⁴Université de Montpellier, UMR 5243 GM (CNRS/UM/UA), Montpellier, France</footnote>

<footnote n="5">⁵INRAE, URFM, Domaine Saint Paul, INRAE Centre de Recherche PACA, Domaine Saint-Paul, France</footnote>

<footnote n="6">⁶Institute of Earth Sciences, University of Lausanne, Lausanne, Switzerland</footnote>

<footnote n="7">⁷Université de Montpellier, UMR 5151 HSM (CNRS/UM/IRD), Montpellier, France</footnote>
using various techniques, including sap flow methods that quantify water movement through the xylem (e.g. Goulden and Field, 1994; Granier et al., 1996; Kume et al., 2010), porometry to assess leaf-level transpiration and stomatal conductance (e.g. Damour et al., 2010; Zhang et al., 1997), and flux towers (eddy covariance) that estimate stand-scale evapotranspiration, with tree transpiration inferred by partitioning soil evaporation (e.g. Kurpius et al., 2003; Scanlon and Kustas, 2012). The heat-balance method is commonly used to measure sap flow velocity in individual plants due to its accuracy, simplicity, and ability to provide continuous measurements (Poyatos et al., 2016; Smith and Allen, 1996). Heat dissipation techniques are widely adopted because of their ease of implementation and high precision (Granier, 1987; Flo et al., 2019; Poyatos et al., 2021; Wang et al., 2023). However, these methods represent point measurements, reflecting localized flow rather than providing an integrated assessment across the entire stem (Flo et al., 2019; Oliveras and Llorens, 2001). They also require species-specific empirical calibrations and regular “zeroing” recalibration during the night (Dix and Aubrey, 2021; McCulloh et al., 2007; Moreno et al., 2021; Poyatos et al., 2016).

In recent years, geophysical methods have gained considerable attention for studying ecohydrological systems (e.g. Carrière et al., 2021a, b; Dumont and Singha, 2024; Harmon et al., 2021; Hermans et al., 2023; Jayawickreme et al., 2014; Loiseau et al., 2023; Luo et al., 2020; Voytek et al., 2019). Among geophysical methods, the passive self-potential (SP) method has shown promise for estimating sap velocity (Gibert et al., 2006; Gindl et al., 1999). Measurable SP signals naturally occur between trees and the surrounding soil in response to bioelectrical effects within plants. Tree SP measurements are often positive when measuring the voltage difference of the stem with respect to the surrounding soil, the so-called biopotentials (e.g. Fensom, 1957, 1963; Gibert et al., 2006; Tattar and Blanchard, 1976; Zapata et al., 2020). Tree transpiration processes facilitate the transport of water and solutes within the xylem (e.g. Kim et al., 2014). Additionally, sugar concentration gradients can generate turgor pressure differences, driving flow within the phloem (e.g. van Bel, 2003). These natural processes can trigger electrokinetic effects through the advection of net electrical charges and electro-diffusive effects driven by electrochemical potential gradients, leading to the generation of biopotentials and measurable SP signals. In theory, SP signals provide an integrated response between the two electrodes that is sensitive to the transpiration process.

SP signals are often measured by non-polarizable electrodes and a high-impedance voltmeter to avoid current leakage arising from the measurements. Due to their stability, non-polarizable Ag/AgCl electrodes have been used to measure SP signals within the trunk in laboratory experiments (e.g. Gil and Vargas, 2023; Gindl et al., 1999; Oyarce and Gurovich, 2010). However, this kind of electrode contains gels or electrolytes that dry up after 1 to 2 months for small-

sized electrodes, thereby leading to unreliable measurements (Fensom, 1963). Such electrodes are hence deemed impractical for prolonged, continuously monitored field experiments. Hubbard et al. (2011) developed a miniaturized Petiau-type Pb/PbCl₂ electrode (Petiau, 2000) for measuring SP induced by a biogeobattery within a small laboratory column, potentially offering a suitable solution with an extended lifespan. However, its installation presents challenges, particularly in shallow depths beneath bark, in ensuring effective contact with the sapwood. Instead, in the natural outdoor environment, polarizable stainless metal electrodes are often used for tree SP measurements (e.g. Fensom, 1963; Gibert et al., 2006; Le Mouël et al., 2024; Zapata et al., 2020, 2021). A non-polarizable electrode is one in which the electrode potential is independent of the current passing through the circuit. However, in both polarizable and non-polarizable electrodes, electrode-related effects are superimposed on the primary effects of interest, often leading to observations that are difficult to interpret. The electrode potentials depend on the electrode design, the temperature, and the surrounding and internal hydrochemical environments (Jougnot and Linde, 2013). Due to electrode-related effects, the measured electrical voltage often shows important drifts during long-term monitoring, related, for instance, to electrode ageing (e.g. Perrier et al., 1997; Hu et al., 2020). Hao et al. (2015) conducted experiments to measure plant-related voltages using electrodes made of different metals and found differences exceeding 10 mV when using copper electrodes with 5 mm differences in diameter. When employing metal electrodes to measure SP on trees, it is hence crucial to ensure consistency in terms of electrode design (materials and geometry). Consequently, the installation and choice of electrodes used to measure tree SP should be made very carefully.

Even in the absence of electrode effects, the interpretation of tree SP signals is complicated due to the different mechanisms contributing to the observed data. These include both electrokinetic and electrochemical processes occurring within the sapwood (Fensom, 1957; Fromm and Lautner, 2007, 2012; Gibert et al., 2006; Gil and Vargas, 2023; Love et al., 2008; Tattar and Blanchard, 1976). Moreover, Barlow (2012) and Le Mouël et al. (2024) suggest that luni-solar tides could induce electric fields, possibly further complicating the interpretation of SP measurements on trees. However, some long-term measurements have not observed clear monthly patterns related to lunar cycles (e.g. Fensom, 1963; Zapata et al., 2021). Uncertainties surrounding these mechanisms make quantitative interpretations of SP data in terms of the water dynamic within the soil-plant system very challenging.

Given the considerable ambiguity surrounding the mechanisms governing the coupling magnitude, which determines the mathematical relationship between SP and sap flow observations, further experimentation is essential to clarify their correlation in situ. The existing literature on plant-related SP predominantly focuses on individual trees (Belashev, 2024;
Gibert et al., 2006; Hao et al., 2015; Le Mouël et al., 2010; Voytek et al., 2019) or different tree species, often with non-continuous or short-term measurements (e.g. Gil and Vargas, 2023; Oyarce and Gurovich, 2010; Pozdnyakov, 2013; Zapata et al., 2020). Recently, Le Mouël et al. (2024) analysed SP signals from various trees using the singular spectrum decomposition method, suggesting that sap flow and multi-scale tree SP may be partially or entirely induced by terrestrial gravitational tides, although without providing data on sap flow. The following pertinent questions remain unanswered regarding tree SP observations: (1) do different tree species exhibit distinct time-varying SP patterns under identical environmental conditions? (2) Do trees of the same species display SP characteristics that depend on the weather conditions? (3) Will two trees of the same species show similar SP responses to sap velocity? (4) Is SP a reliable indicator of sap velocity within the trunk? (5) When are electrokinetic effects dominating tree SP responses?

To address these questions, we review tree SP theory in Sect. 2, followed by a description of long-term SP monitoring experiments conducted at three distinct sites in the French Mediterranean region, along with an outline of data analysis methods, in Sect. 3. In Sect. 4.1, we present and analyse the data obtained from measurements on four different trees at three sites. Section 4.2 and 4.3 investigate the short-term characteristics of tree SP in terms of correlations with sap velocity, precipitation, evapotranspiration, and other meteorological data. Finally, Sect. 5 discusses these results by quantitatively relating tree SP signals to various electrical potential sources such as hydraulic properties and electrochemical effects.

## 2 Theoretical background concerning tree SP

Pioneering research of tree SP monitoring dates back to the early 1960s, when Fensom (1963) embarked on extensive studies to monitor SP on different trees. Earlier studies of the bioelectrical potential can be found in Fensom (1957). Fensom (1963) observed that the electrical voltages were the strongest during growing seasons and proposed that voltage measurements of plants might detect cambial growth. Subsequent studies by Gindl et al. (1999) measured both sap velocity and SP on tree stems in a well-controlled laboratory experiment, uncovering a significant linear correlation between the two. They suggested that the transpiration process acts as a driving force in xylem sap, generating electrokinetic effects that explain the relationship between the sap velocity and the SP data, suggesting the potential value of using SP data to estimate plant sap flow. Through continuous observations of an individual poplar tree at various heights and positions, Gibert et al. (2006) found that the relationship between the SP and the sap flow velocity is not entirely linear. They suggested the presence of electrically active structures within the xylem, leading to electro-diffusive effects within plants.

Voytek et al. (2019) simulated the SP generated by root water uptake and reconstructed the observed time-varying patterns of SP in the root zone, confirming that plant root water uptake affects the observed electrical voltages. However, the numerical simulations accounting for the electrokinetic coupling process induced by the soil water movement in the root zone could not predict the amplitude of the observed SP.

Previous experimental studies on SP between the living plant and the surrounding soil have shown that the electrical voltages generated by physiological, physical, and chemical processes can range from tens to hundreds of millivolts (e.g. Fensom, 1963; Pozdnyakov, 2013; Zapata et al., 2020). These magnitudes are comparatively high compared to SP signals generated by near-surface water or solute flux (e.g. Jougnot et al., 2015; Voytek et al., 2019; Hu et al., 2020). A schematic diagram of a tree is presented in Fig. 1 to summarize the mechanisms of SP generation. SP differences from root to shoot are induced within the sapwood, including the xylem and phloem (Fig. 1b and c). Below, we describe the possible mechanisms generating tree SP.

### 2.1 Self-potential signals generated by xylem and phloem transport

Observed SP signals arise from the superposition of contributions arising from different coupling sources and from electrode-related effects (e.g. Hu et al., 2020; Jougnot et al., 2015; Jougnot and Linde, 2013). The SP signal due to flow and mass transport process can be described by the following Poisson equation in the absence of external electrical current:

$$ \nabla \cdot (-\sigma E) = \nabla \cdot (J_{\text{ek}} + J_{\text{diff}}), \quad (1) $$

where $\sigma$ (S m⁻¹) is the bulk electrical conductivity; $E = -\nabla V$ (V m⁻¹) is the electric field; $\nabla V$ is the electrical potential gradient; and $J_{\text{ek}}$ (A m⁻²) and $J_{\text{diff}}$ (A m⁻²) are the electrokinetic and electro-diffusive current density, respectively. We detail these two contributions below.

The primary contribution of sap ascent occurs in the xylem, where vessels transport water and dissolved minerals from the roots throughout the plant. This process relies on transpiration pulling, creating a negative pressure gradient that drives upward flow through the xylem (Fig. 1c). The xylem vessels feature a capillary-like structure, as illustrated in Fig. 1d. Although the cellulose walls of the xylem are electrically neutral and insulating, they become negatively charged as they attract water molecules (Fensom, 1957). For negatively charged cell walls, the streaming current, generated by the upward movement of excess charge within the Gouy-Chapman diffuse layer as sap flows, represents the advective transport of electrical charges, resulting in a net source of current density $J_{\text{ek}}$. According to the electrokinetic mechanism in the presence of an electrical double layer, the upward flow of sap induces a natural electric field in the opposite direction, as depicted in Fig. 1e.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** The schematic diagram illustrates the components contributing to tree SP. (a) Illustration of a tree, adapted from “Ash tree” created in BioRender (https://BioRender.com/kixmyvj, Hu, 2025a); (b) horizontal and vertical cross sections of the trunk; (c) vertically extracted xylem and phloem, adapted from “Plant Vessels Layout” created in BioRender (https://BioRender.com/0qwjmt9, Hu, 2025b); (d) cellulose structure, adapted from “Cell wall (cellulose)” created in BioRender (https://BioRender.com/76osafb, Hu, 2025c); and (e) a schematic of the electrical double layer.</figcaption></figure>

The electrokinetic cross-coupling coefficient $C_{ek}$ (V Pa$^{-1}$) can be estimated by the ratio of voltage difference $V_{z1} - V_{z2}$ to pressure difference $p_{z1} - p_{z2}$ with no external currents (i.e. electrical flux equals 0). This describes a condition where only electrokinetic effects are present, with no contribution from diffusive currents in Eq. (1), ensuring that the net current remains zero. The $C_{ek}$ is determined using the Helmholtz–Smoluchowski (HS) equation:

$$C_{ek} = \left. \frac{V_{z1} - V_{z2}}{p_{z1} - p_{z2}} \right|_{J=0} = \frac{\epsilon_f \zeta}{\sigma_f \eta_f}, \qquad (2)$$

where $\epsilon_f$ (F m$^{-1}$) and $\sigma_f$ (S m$^{-1}$) denote the dielectric permittivity and the electrical conductivity of the capillary flowing fluid, respectively. The parameter $\zeta$ (V) is the zeta potential at the shear plane of the electrical double layer (Fig. 1e); it typically has a negative sign. Both $\sigma_f$ and $\zeta$ are influenced by temperature and ion concentration. Consequently, the voltage induced by xylem transport depends not only on sap velocity but also on sap temperature and concentration. Based on laboratory experiments under different pore water ionic concentrations, Linde et al., (2007) proposed an empirical relationship

between the logarithm of $C_{ek}$ and $\sigma_f$ for porous media. A similar model could be devised to estimate $C_{ek}$ within vascular plants.

Alternatively, the electrokinetic current density $J_{ek}$ (A m$^{-2}$) can be expressed by the product of the effective excess charge density $\hat{Q}_v$ (C m$^{-3}$) and the Darcy velocity $u$ (m s$^{-1}$) (e.g. Jougnot et al., 2020):

$$J_{ek} = \hat{Q}_v u. \qquad (3)$$

Here, $\hat{Q}_v$ indicates the density of excess charge in the diffuse layer of the electrical double layer model (Fig. 1e) that is effectively dragged by sap flow. Following Kormiltsev et al. (1998) and Jougnot et al. (2019), $C_{ek}$ can be transformed into $\hat{Q}_v$ by

$$\hat{Q}_v = - \frac{\eta_f \sigma}{k} C_{ek}, \qquad (4)$$

where $k$ (m$^2$) denotes the permeability of a porous medium. Under the assumption of 1-D flow and neglecting $J_{diff}$ in Eq. (1), the current density is solely governed by the electrokinetic effect. The equivalent current density induced by
the sap flow can be simplified to the following:

$$ \sigma E = -\sigma \nabla V = -\hat{Q}_v u. \qquad (5) $$

Thus, the $\hat{Q}_v$ can be expressed as (e.g. Jougnot et al., 2012)

$$ \hat{Q}_v = \frac{\sigma \nabla V}{u}. \qquad (6) $$

During the transpiration process, water loss coupled with osmotic changes leads to an increase in solute concentration towards the crown. This implies the potential generation of electro-diffusive potentials due to electrochemical differences. The xylem sap typically contains approximately 10 mM of inorganic nutrients, along with organic nitrogen compounds metabolically synthesized in the root (Nobel, 2009). Common ion species found in plants include protons ($H^+$), calcium ($Ca^{2+}$), potassium ($K^+$), magnesium ($Mg^{2+}$), ammonium ($NH^{4+}$), chlorine ($Cl^-$), and nitrate ($NO^{3-}$), among others (Davies, 2006; Miller and Wells, 2006; Volkov and Markin, 2012). For instance, increases in $K^+$ of xylem sap occur with increasing transpiration demand and consequently increase hydraulic conductivity and capacity of water movement (e.g. Losso et al., 2023; Nardini et al., 2011). However, the concentration distribution along the longitudinal direction exhibits different trends and patterns among different trees (McDonald et al., 2002). According to McDonald et al. (2002), for a Norway spruce tree in late autumn, the total concentration of amino acids, the minerals magnesium, calcium, and potassium in the xylem sap, ranges from 4 mM at the ground level to 6 mM near the crown base. This concentration range aligns with typical salinity levels in pore water (Linde et al., 2007). The electro-diffusive current density $J_{diff}$ (A m$^{-2}$) arises from the electrochemical potential gradient, which is influenced by the activities of different ion species:

$$ J_{diff} = - \frac{k_B T_K}{e_0} \sigma \sum_{i=1}^{N} \frac{t_i^H}{q_i} \nabla \ln a_i, \qquad (7) $$

where $a_i$, $t_i^H$, and $q_i$ represent the activity (i.e. thermodynamically effective concentration), the microscopic Hittorf transport number, and the valence of the ion species $i$, respectively. $k_B$ (kJ$^{-1}$), $e_0$ (C), $T_K$ (K), and $\sigma$ (S m$^{-1}$) are the Boltzmann constant, the elementary charge, the temperature, and bulk electrical conductivity. $N$ is the total number of ion species in xylem sap. The equivalent electrical potentials induced by $J_{diff}$ are expressed by

$$ \nabla \varphi_{diff} = \frac{J_{diff}}{\sigma} = - \frac{k_B T_K}{e_0} \sum_{i=1}^{N} \frac{t_i^H}{q_i} \nabla \ln a_i. \qquad (8) $$

Compared with the electro-diffusive effect, cell membrane potentials across xylem tracheids constitute another source of electrical potential gradient (Fensom, 1957; Hedrich and Schroeder, 1989). Further details on the membrane potential can be found in Nobel (2009) and Spanswick (2006).

Water movement within the xylem is mainly sustained by hydrostatic pressure gradients driven by transpiration processes. Conversely, in the phloem, water flow is influenced not only by hydrostatic pressure gradients but also by osmotic pressures resulting from translocation processes. Specifically, the loading of sugars into the phloem leads to an increase in solute concentration but a decrease in water potential in source cells, creating a higher water potential in sink cells (Fig. 1c). Consequently, water and solutes move upward in the xylem and downward in the phloem (Hölttä et al., 2006). Moreover, osmotic pressures in the phloem cause water from the xylem to radially flow into the phloem (Fig. 1c). Apart from the passive electrokinetic and electro-diffusive mechanisms, actively transporting charged particles into cells may lead to the accumulation of internal electrical potential and the occurrence of electrogenic pumps across membranes (Nobel, 2009).

## 2.2 Previous studies and discussion of tree SP generation

Multiple studies have explored the relationship between sap flow and electrical potential differences along plant stems, often attributing it to electrokinetic effects (e.g. Gindl et al., 1999; Gil and Vargas, 2023; Gibert et al., 2006; Fensom, 1963; Koppán et al., 2002). However, this explanation has come under scrutiny due to inconsistencies observed in measurements. Love et al. (2008) conducted a laboratory experiment involving *Ficus benjamina* trees and their surrounding soil for different pH levels. They proposed that disparities in electrical potential differences between the plant and its surrounding soil are predominantly influenced by the pH difference between the xylem and the soil. Similarly, Zapata et al. (2020), based on measurements from various trees in a Mediterranean forest region, suggested abandoning the electrokinetic mechanism to explain SP differences along the trunk. However, their measurements were conducted over a limited period. Furthermore, fluctuations in pH alter the concentrations of $H^+$ and $OH^-$ ions in aqueous solutions, consequently affecting the zeta potential of the electrical double layer model (e.g. Al Mahrouqi et al., 2017; Leroy et al., 2013). As a result, pH variations contribute to changes in the electrokinetic coupling coefficients and the magnitudes of corresponding SP.

Aside from the observed longitudinal differences in electrical potential, measurements taken across the stem's cross-sectional profile reveal distinct radial polarization (Islam et al., 2017; Zapata et al., 2020). Islam et al. (2017) conducted measurements of short-circuit currents radially along the trunk's growth direction, noting a peak current near the vascular cambium – a layer between the secondary xylem and phloem (refer to Fig. 1b). They observed that strong electrical current generation may result from the low electrical resistance of living tissues, the influx of water into the phloem, and the diffusion of charges around the cambium re-
region. These findings suggest that the measured tree SP also varies with respect to the depth at which electrodes are inserted into the trunk.

As the measured SP is influenced by multiple mechanisms, its composition changes with water and solute movements, exhibiting spatial and temporal variations. Interpreting tree SP observations becomes complex due to these multifaceted mechanisms, which exhibit different patterns under various seasonal and meteorological conditions. Rather than relying on a single generative mechanism and short-term observations, we present below results obtained from continuous, long-term monitoring of tree SP, sap velocity, and meteorological data to comprehensively analyse their correlations over time.

# 3 Materials and methods

## 3.1 Study sites and trees

Three distinct experimental sites were selected in the French Mediterranean region to monitor sap velocity and SP on four trees from three different species (Fig. 2). The characteristics of each site are summarized in Table 1. All sites are located in karstic hydrosystems. The Larzac and LSBB test sites are part of the H+ National Observatory Service (https://hplus.ore.fr, H+, 2025) and the OZCAR research infrastructure (https://www.ozcar-ri.org, OZCAR, 2025; Gaillardet et al. 2018), while the Font-Blanche site belongs to the ICOS research infrastructure (Gielen et al., 2017).

The Font-Blanche Forest, located 30 km east of Marseille and approximately 7 km from the coastline (Fig. 2a), sits at an elevation of 425 m (Ollivier et al., 2021). The bedrock consists of sub-reefal Cretaceous limestone with Urgonian facies. The forest is dominated by Aleppo pines (*Pinus halepensis* Mill.), averaging 13 m in height, with an understorey of holm oaks (*Quercus ilex* L.) with an average height of 6 m (Girard et al., 2012; Simioni et al., 2020). Electromagnetic induction and electrical resistivity tomography (ERT) characterized the spatial variability as a function of soil depth (Carrière et al., 2021a).

The LSBB site, located in the southern part of the Fontaine-de-Vaucluse karst hydrosystem at an altitude of 422 m, is dominated by holm oak (*Quercus ilex* L.), averaging 13 m in height. The bedrock composition is similar to Font-Blanche. Detailed studies of the site's structure and hydrodynamics were conducted by Carrière et al. (2016) using ERT, ground-penetrating radar, and magnetic resonance sounding. Carrière et al. (2020a, b) examined interindividual drought responses along a transect monitored by ERT. Evapotranspiration at the mountain scale was studied using superconducting gravimetry (Carrière et al., 2021b).

The Larzac observatory, situated at an average altitude of 708 m, features a landscape of grasslands, black pine plantations, and deciduous stands dominated by pubescent oaks

(Quercus pubescens), with monitoring focused solely on the oaks. The observatory is part of the Durzon karstic system, primarily composed of dolomite, and has been explored using gravimetry (Fores et al., 2018; Jacob et al., 2008), ERT, and seismic methods (Valois et al., 2016).

Four trees have been equipped with SP and sap velocity sensors (Fig. 2 and Table 1): (1) *Quercus pubescens* at the Larzac observatory with a trunk circumference of 90 cm; (2) *Quercus ilex* L. at the LSBB site with a circumference of 43 cm; and (3) *Pinus halepensis* Mill. and *Quercus ilex* L., with circumferences of 108.3 and 39.6 cm, respectively, at the Font-Blanche site. These measurements allow us to compare the physiological characteristics of different species under different meteorological and geological conditions.

## 3.2 Instrument and data

### 3.2.1 Sap velocity measurements

Data were collected at the three sites throughout 2023 (Table S1 in the Supplement). The raw data related to this study can be found in Hu et al. (2024). To assess transpiration rates, sap velocity was measured using the heat dissipation technique by Granier (1987) at the four trees. Two thermocouple probes were inserted 2 cm into the sapwood of each tree, positioned 10 cm apart. The upper probe applied a constant heating power (see Fig. 2e). The temperature difference between the upper and lower probes is influenced by sap flow surrounding the heated probe. The Granier-type heat dissipation probes were installed at around 1.5 m above the ground of the trees. Thus, by continuously heating in the upper probe, sap velocity could be estimated by measuring the temperature difference between the probe pairs by copper-constantan thermocouples under the thermal imbalance. The original signal was recorded in millivolts and subsequently transformed into temperature with a calibration coefficient (~0.04 mV °C⁻¹; Do and Rocheteau, 2002). Then, the sap velocity was calculated by Granier's empirical equation (Granier, 1987) with the obtained temperature differences. After calibration and correction, the sap velocity is obtained in µms⁻¹ at 30 min intervals (e.g. Moreno et al., 2021).

### 3.2.2 SP measurements

The electrical voltages were measured using a Campbell Scientific CR1000X datalogger, using the voltmeter function with an input impedance of 20 GΩ. The sampling time intervals for tree SP measurements are 1 min at the Larzac and LSBB sites and 10 min at the Font-Blanche site. Details about the measurement parameters, sampling intervals, and units can be found in Table S1. Each tree is equipped with a pair of stainless-steel screws with 3 mm diameter, separated vertically by 10 cm, serving as trunk electrodes (Fig. 2f). These electrodes were horizontally aligned with the corre-
<table><caption>**Table 1.** Features of experimental sites and trees.</caption><thead><tr><th rowspan="2">Property</th><th colspan="3">Sites</th></tr><tr><th>Larzac</th><th>LSBB</th><th>Font-Blanche</th></tr></thead><tbody><tr><td>Tree number</td><td>C2</td><td>A3</td><td>P327</td><td>C325</td></tr><tr><td>Abbreviations</td><td>LaQp</td><td>LSQi</td><td>FBPh</td><td>FBQi</td></tr><tr><td>Tree species</td><td>pubescent oak<br>(Quercus pubescens)</td><td>holm oak<br>(Quercus ilex L.)</td><td>Aleppo pine<br>(Pinus halepensis Mill.)</td><td>holm oak<br>(Quercus ilex L.)</td></tr><tr><td>Circumference of trunk (m)</td><td>0.900</td><td>0.428</td><td>1.083</td><td>0.396</td></tr><tr><td>Xylem anatomy</td><td>ring-porous</td><td>ring-porous</td><td>diffuse-porous</td><td>ring-porous</td></tr><tr><td>Electrical resistivity of sapwood (Ω m)</td><td>282</td><td>—</td><td>369</td><td>250</td></tr><tr><td>Bed rock</td><td>dolomite</td><td>limestone</td><td></td><td></td></tr><tr><td>Mean altitude (m)</td><td>708</td><td>422</td><td></td><td>425</td></tr><tr><td>Soil depth (cm)</td><td>0–80</td><td>0–70</td><td></td><td>0–60</td></tr><tr><td>Mean annual precipitation (mm)</td><td>870</td><td>900</td><td></td><td>650</td></tr><tr><td>Mean annual actual evapotranspiration (mm)</td><td>600</td><td>580</td><td></td><td>460</td></tr><tr><td>Mean annual temperature (°C)</td><td>10.5</td><td>13</td><td></td><td>14</td></tr></tbody></table>

sponding sap flow sensors on each tree (Fig. 2). After peeling away the bark from the trunk, the electrodes were inserted to a depth of 2 cm to ensure contact with the sapwood. The outside part of the electrodes was coated with black insulating rubber (as depicted in Fig. 2f). All trunk electrodes were identical in material, size, and shape. Following installation, the exposed portions of the electrodes were further shielded with epoxy glue to isolate them from environmental factors and minimize the impact of meteorological conditions such as air temperature, humidity, and radiation.

Polarizable stainless-steel electrodes were used in previous outdoor experiments for tree SP measurements, such as Gibert et al. (2006) and Zapata et al. (2021). Gibert et al. (2006) discussed the temperature effects of steel electrodes and found that artificially heating and cooling a tree electrode did not result in significant variations in electrical potential compared to another electrode that was not thermally disturbed. In this work, we neglect temperature effects on electrode potentials.

Non-polarizable Petiau-type (Pb/PbCl₂) electrodes (Petiau, 2000) were selected as reference electrodes and buried in the soil at a depth of 37 cm (Fig. 2g). The raw signal is a voltage between the trunk and the reference, which includes the electrical potential of the tree in situ and the electrode. In this study, we used the voltage measured at the lower electrode minus the higher electrode to indicate the tree SP.

### 3.2.3 Meteorological data

The physiological conditions of trees respond to climate change. We gathered data on actual evapotranspiration (actual ET), precipitation (Prec.), vapour pressure deficit (VPD), and air temperature (Temp.) from the meteorological station located at the Font-Blanche site. Additionally, a meteorological station near the LSBB site supplied us with daily precipitation, mean air temperature, and global radiation data.

Hourly air temperature and precipitation data were obtained from a meteorological station situated less than 10 km from the Larzac site (see Table S1 for detailed information).

## 3.3 Data analysis

### 3.3.1 Wavelet analysis of sap velocity and tree SP

To identify correlations and phase lags between sap velocity and tree SP at different timescales, we calculated the coherences between the two time series using continuous wavelet transformation coefficients. The wavelet coherence $C_{SV,SP}(a, b)$ between sap velocity (SV) and SP signals is calculated by

$$C_{SV,SP}(a, b) = \frac{|S(W_{SV}(a, b) W_{SP}^*(a, b))|^2}{S(|W_{SV}(a, b)|^2) S(|W_{SP}(a, b)|^2)}, \quad (9)$$

where the superscript * denotes the complex conjugate, S denotes the smoothing operator function (Torrence and Compo, 1998), and $W_{SV}(a, b)$ and $W_{SP}(a, b)$ are the continuous wavelet transforms of sap velocity and SP signals at the wavelet scale a and position (frequency or period) b, using the Morlet wavelet (e.g. Grinsted et al., 2004; Linde et al., 2011).

### 3.3.2 Variational mode decomposition

Due to factors such as temperature variations, trunk wounds, and electrode polarization effects, the recorded SP data may exhibit long-period drifts. Additionally, SP measurements are susceptible to high-frequency electromagnetic noise interference. For example, Fourier spectrum analysis of SV data for FBPh and FBQi through 2023 indicates the presence of 1, 2, 3, and 4 d⁻¹ frequency signals (Fig. S1a and b in the Supplement). However, the SP spectrum (Fig. S1c and d) does not exhibit clear dominant frequencies that can be directly characterized through Fourier analysis. As introduced in Sect. 2.1, the generation mechanisms of tree SP
<figure><img src="image_4.png" /><figcaption>**Figure 2.** The observation system at the three forest sites includes the following components: (a) location map of the sites, indicated by red circles; (b–d) observation setups at the Larzac site, the LSBB site, and the Font-Blanche site, respectively; (e–g) sensors for measuring sap velocity, electrodes for measuring tree SP, and reference electrodes buried underground.</figcaption></figure>

are complex and multiple. This means that the measured SP in our experimental sites contains signal contributions originating from various sources. Similar to an electrocardiogram (ECG), the tree SP is an electrophysiological signal. The variational mode decomposition (VMD) approach has been applied to process ECG signals and extract the intrinsic modes from signals to obtain information-containing spikes and low-frequency patterns (long-term drifts) related to physiological activities.

VMD can be used to decompose signal $f(t)$ into sub-signals $u_k(t)$ with different centre frequencies $\omega_k$ (Dragomiretskiy and Zosso, 2014):

$$f(t) = \sum_{k=1}^{N} u_k(t), \qquad (10)$$

where $k$ and $N$ denote the number and total number of decomposed modes. Assuming each mode $u_k$ to be concentrated around a centre frequency $\omega_k$ in its bandwidth, the following minimization problem is used to determine the de-
composition $u_k$:

$$ \min_{\{\omega_k, u_k\}} \left\{ \sum_{k=1}^{N} \left\| \partial_t \left[ \left( \delta(t) + \frac{j}{\pi t} \right) * u_k(t) \right] e^{-j\omega_k t} \right\|_2^2 \right\}, \quad (11) $$

where $j^2 = -1, \delta(t)$ and $*$ denote the Dirac function and convolution, respectively. More details about the optimization process to solve Eq. (11) involving the updating of the modes and the centre frequencies can be found in Dragomiretskiy and Zosso (2014). In this study, we use the frequencies with the centre of gravity and the largest amplitudes of the mode's power spectrum to represent the centre frequency and the dominant frequency of each corresponding mode, respectively.

# 4 Results

## 4.1 Overview of 1 year of raw data

### 4.1.1 Raw time series data

The tree SP data collected at the Font-Blanche site are displayed in Fig. 3 (refer to Sect. 3.2.2). The time-varying sap velocities and SP obtained on the Aleppo pine trunk (Fig. 3c and e) differ from those of the holm oak trunk (Fig. 3d and f), even when the two trees are in proximity and experience identical weather conditions (Figs. 2d and 3a, b). Both trees show rather large SP signals with a more pronounced response to environmental changes on Aleppo pine compared to holm oak. However, the overall time-varying patterns of the raw SP data from trees do not demonstrate a clear correlation with sap velocity, as evidenced by relatively low Pearson correlation coefficients (see Fig. S2). Additional data from the Larzac and LSBB sites are provided in Fig. S3.

### 4.1.2 Wavelet analysis of sap velocity and tree SP

According to Sect. 3.3.1, we calculated the wavelet coherence between the sap velocity and tree SP data collected throughout 2023. As an example, we present the coherence maps for the two trees at the Font-Blanche site in Fig. 4. Wavelet coherence indicates a localized correlation in time-frequency space (e.g. Grinsted et al., 2004). Figure 4b and c illustrate the temporal correlation between SV and SP across different time periods. The arrows in Fig. 4b and c represent the phase relationship between SV and SP. Particularly, the rightward arrows (0° phase) indicate in-phase behaviours, where SV and SP vary synchronously. Leftward arrows (180° phase) indicate anti-phase behaviours, signifying negative correlation or out-of-phase variations between SV and SP. Conversely, upward or downward arrows (±90° phase) reflect phase shifts, suggesting that one signal leads or lags the other by a quarter cycle. In particular, the 1 d period shows high coherence, with arrows pointing left, indicating that SP

is negatively correlated with SV at this timescale. This indicates that SP and SV exhibit stronger correlation at a daily timescale, compared to other periods. Notably, tree SP exhibits a negative correlation with sap velocity across the two different species of trees. Note also that the coherence diminishes during periods of water surplus (Fig. 4a). In contrast, the holm oak at the LSBB site shows lower coherence and the pubescent oak at the Larzac site lacks coherence between sap velocity and SP (Fig. S4). The wide discrepancies for different trees at different test sites might be attributed to different annual rainfall (Table 1 and Fig. S3a–b) and other factors such as electrochemical and electrode-related effects contributing to the measured SP data.

## 4.2 Decomposition results

In our case, VMD is used to separate the measured data into a finite set of components (modes). We initially assume six principal modes, as all coefficients of determination ($R^2$) between the measured tree SP, sap velocity data, and the sum of the corresponding six intrinsic modes using 1-year data from 2023 for all four trees are 1.00, meaning that the results of VMD could effectively reconstruct the measurements. As an example, we present the frequency spectra of six modes decomposed from the tree SP and sap velocity (SV) obtained from the holm oak tree at the Font-Blanche site throughout 2023 (Fig. 5). The centre and dominant frequencies of each decomposed mode decrease sequentially. Besides the expected diurnal variations driven by the rhythm of sap flow (Fig. 5i and j), other high-frequency content and low-frequency trends are considered noise and long-term drifts in this study.

When compared to the decomposition spectrum of SV data (Fig. 5b and d), noise present in the tree SP data tends to manifest at higher frequencies, as indicated by the prominent spectral amplitudes in the higher-frequency range (Fig. 5a, c, e, and g). The low-frequency components, which reflect baseline trends, dominate the last intrinsic modes and occupy the largest proportion of the signals, as shown by the highest amplitudes in the SP and SV spectra (Fig. 5k and l). Apart from the fifth intrinsic mode, the decompositions of SP and SV exhibit distinct frequency characteristics, highlighting the complexity of SP signals. Notably, the fifth intrinsic modes for both datasets reveal a strong diurnal rhythm, with a dominant frequency of $1 d^{-1}$ (Fig. 5i and j). The spectra of the fifth decomposition also display the second-highest amplitudes, indicating that these diurnal components represent a significant portion of the overall SP and SV signals.

We investigated the impact of increasing the number of modes from 6 to 12 when performing data decomposition. Regardless of the number of modes used, the results consistently indicate a diurnal rhythm in the second-to-last intrinsic mode, with similar amplitudes and patterns (refer to Fig. S5). We also conducted an analysis on the tree SP and sap velocity data collected from the neighbouring Aleppo pine (FBPh)
<figure><img src="image_5.png" /><figcaption>**Figure 3.** A dataset of 1-year data collected at the Font-Blanche site at half-hourly intervals from 1 January 2023 to 1 January 2024. (a) Precipitation and air temperature data; (b) actual evapotranspiration (actual ET) and vapour pressure deficit (VPD) data; (c, d) measured sap velocity in blue lines and daily cumulative sap flow in pink lines, for the Aleppo pine (FBPh) and the holm oak (FBQi), respectively; (e, f) measured SP data in blue lines and daily average of SP data in pink lines, for the Aleppo pine (FBPh) and the holm oak (FBQi), respectively.</figcaption></figure>

using the VMD approach (Fig. S6). The frequency bands of each decomposition derived from FBPh are comparable to those from FBQi. Particularly, the high-frequency oscillating characteristics of SP components on both trees are nearly identical. Since our primary interest lies in understanding the diurnal variations associated with transpiration, we employed the same VMD procedure on the whole datasets to extract the fifth decomposed results based on decompositions with six modes.

## 4.3 The relationship between SP and physiological characteristics of trees

### 4.3.1 Correlation coefficients between data collected at the Font-Blanche site

Tree transpiration is affected by water availability, VPD, net radiation, and atmospheric uptake of water (e.g. Alfieri et al.,

2018). Employing the VMD approach and extracting six intrinsic mode functions, we delineated the diurnal rhythms using the fifth subset component of tree SP, sap velocity, and meteorological data collected during a 2-week period in the growing season of 2023. The 2-week period from 16 to 30 April was selected to capture a phase of active transpiration across different species, during which no precipitation was recorded at the Font-Blanche site (as shown in Fig. 3). The large long-term fluctuations in raw SP data across positive and negative ranges and high-frequency interferences are excluded by extracting the fifth decomposed mode. This allows us to focus on relative variations aligned with the dominant sap flow frequency, providing clearer insight into the relationship between SP and SV.

The fifth VMD component exhibited significant correlations between tree SP and sap velocity, notably at the Aleppo pine (FBPh) with a correlation coefficient of -0.88. The
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Wavelet coherence analysis between sap velocity and tree SP data at the Font-Blanche site in 2023. (a) Water surplus, calculated as the difference between precipitation and actual evapotranspiration. Wavelet coherence maps for the (b) Aleppo pine and (c) holm oak at the Font-Blanche site, with yellow highlighted boxes indicating periods of water surplus. Arrows denote the lag/lead phase between the two time series; dashed white lines (b, c) indicate the cone of influence within which edge artefacts are negligible.</figcaption></figure>

cross-correlation coefficients between decomposed sap ve-
locity and SP for the holm oak (FBQi) trees reached –0.66.
Based on the electrokinetic mechanism generated by sap flow
in Sect. 2.1, a larger negative SP should indicate a higher
transpiration rate, which is indicated by the negative signs
of the Pearson correlation between the sap velocity and tree
SP in this case. Interestingly, the fifth extracted VMD com-
ponents of tree SP data on FBPh displayed relatively high
correlations with other data, including tree SP at FBQi (see
Fig. S7). Below, we further analyse the correlations between
the time-varying tree SP and the sap velocity data on differ-
ent trees and sites.

4.3.2 Relationship between tree SP and sap velocity
across four trees

To analyse the linear relationship between sap velocity and
SP, scatter plots are displayed in Fig. 6 using the fifth VMD
modes. The fifth decomposed modes of SP and SV signals
are denoted by “VMD<sub>5</sub><sup>SP</sup>” and “VMD<sub>5</sub><sup>SV</sup>”, respectively. The
fifth VMD results of tree SP and SV data present similar
frequency content (Fig. 5i and j). Here, we assume that the
fifth VMD result of tree SP originates from diurnal varia-
tions of sap velocity. Please note that the negative values in

the fifth decomposed mode of SV signals represent relative
deviations from the baseline, reflecting diurnal fluctuations
rather than indicating actual reverse sap flow. The baseline
component is captured in the sixth decomposed mode, iso-
lating physical sap flow from the oscillatory variations seen
in the fifth decomposed mode. Upon applying linear regres-
sion of SP versus SV, for two species of trees investigated
at the Font-Blanche Forest, the holm oak (Fig. 6c) exhibits
a steeper slope compared to the Aleppo pine (Fig. 6a). Out-
liers were removed under a 95 % confidence level by the lin-
ear regression to obtain the slope coefficients. In contrast, the
correlation of decompositions between the tree SP and SV is
lower on the holm oak at the LSBB (Fig. 6e) than on the
same species at the Font-Blanche (Fig. 6c). Additionally, no
correlation is found on the pubescent oak at the Larzac site
(Fig. 6g).

Using the parameters obtained from the linear regression
(Fig. 6a, c, e, and g), we predicted the tree SP responses as a
function of the sap velocity (refer to Fig. 6b, d, f, and h). The
predictions effectively replicate the time-varying SP patterns
observed in the Aleppo pine at the Font-Blanche (Fig. 6b).
Rainfall events decreased the coherence and increased the
phase shifts between sap velocity and tree SP (Fig. 4). In this
case, due to the rainfall events during the extracted period at
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Frequency spectra of six decomposed modes of the left column (a, c, e, g, i, k) tree SP and the right column (b, d, f, h, j, l) sap velocity data, obtained on the holm oak (FBQi) at the Font-Blanche site within 2023, using the variational method (VMD). Different rows correspond to different modes, where “CF” and “DF” indicate the centre frequency and dominant frequency of the corresponding mode, respectively.</figcaption></figure>

both LSBB and Larzac sites, the VMD<sub>5</sub><sup>SV</sup> could not restore the VMD<sub>5</sub><sup>SP</sup> obtained from the measured data (Fig. 6f and h). Additionally, despite the absence of rainfall, the modelled tree SP for the holm oak at the Font-Blanche fails to accurately capture the highly varying amplitudes observed in the measured SP data (Fig. 6d). Specifically, elevated SP amplitudes are noted on 21, 22, and 24 April.

Since the amplitude of the calculated SP (derived from the fifth decomposed mode of SV signals, VMD<sub>5</sub><sup>SV</sup>) is lower than that of the fifth decomposed mode of SP signals (VMD<sub>5</sub><sup>SP</sup>) for holm oak at the Font-Blanche site on 21, 22, and 24 April (Fig. 6d), we conducted a test to investigate whether the VMD method reduces the relative amplitude of SV signals during decomposition. To achieve this, we compared the raw SP and SV data to the sum of the diurnal features (VMD<sub>5</sub><sup>SP</sup>, VMD<sub>5</sub><sup>SV</sup>) and the baseline component (VMD<sub>6</sub><sup>SP</sup>, VMD<sub>6</sub><sup>SV</sup>) obtained from the VMD-processed data (Fig. 7a and b). The combination of the fifth and sixth decomposed modes is anal-

ogous to low-pass filtering, retaining lower-frequency signals. The highest diurnal amplitudes of the raw and processed SV data remain relatively stable, whereas SP data show notable variations, particularly on 21, 22, and 24 April.

To further investigate whether the underestimation of SP using VMD<sub>5</sub><sup>SV</sup> resulted from the exclusion of high-frequency components, we examined the sum of the first to fifth decomposed modes of SV signals, representing a detrending process. In Fig. 7c, the fifth decomposed mode of SP (VMD<sub>5</sub><sup>SP</sup>) was compared with the SP calculated from $\sum_{i=1}^{5} \text{VMD}_i^{\text{SV}}$. This comparison demonstrated that even when all decomposed modes except the baseline were used to reconstruct diurnal SP from SV data, the higher SP amplitudes observed on 21, 22, and 24 April could not be reproduced (Fig. 7c). To corroborate this observation, we repeated the analysis using a finite impulse response (FIR) filter, which yielded similar results (Fig. 7d). This suggests that the elevated SP magnitudes
<figure><img src="image_8.png" /><figcaption>**Figure 6.** The relationship between the fifth decomposed modes of sap velocity ($VMD_5^{SV}$) and tree SP ($VMD_5^{SP}$) signals on 16–30 April 2023, at the Font-Blanche, LSBB, and Larzac sites. (a, c, e, g) Scatter plots of the fifth decomposed modes of sap velocity and SP signals. For (a) the Aleppo pine (FBPh), (c) the holm oak (FBQi) at the Font-Blanche site, (e) the holm oak at the LSBB (LSQi), and (g) the pubescent oak at the Larzac (LaQp). Black lines indicate linear regression results and red dots indicate points (outliers) outside the 95% confidence level of the regression. (b, d, f, h) Corresponding decompositions of SP signals (orange dots) and calculated SP (black lines) based on the linear relationship between $VMD_5^{SV}$ and $VMD_5^{SP}$, with red dots as outliers outside the 95% confidence level of the linear regression.</figcaption></figure>

observed on these dates may be influenced by factors other than sap flow.

As detailed in Sect. 3, the height of the lower electrode with respect to the upper electrode on the monitored trunks is –10 cm, allowing us to calculate the electrical potential gradient ($\nabla V$) between the measurement points using $VMD_5^{SP}$. Moreover, electrical resistivity tomography has been conducted on trunks at both the Font-Blanche and Larzac sites. The average sapwood resistivities of the Aleppo Pine and the holm oak at the Font-Blanche are 369 and 250 Ω m, respectively (e.g. Moreno et al., 2021). By converting the localized $\nabla V \approx -\frac{VMD_5^{SP}}{10[\text{cm}]}$ and $v \approx VMD_5^{SV}$ to the International System of Units, we estimated the effective excess charge density ($\hat{Q}_v = \frac{\sigma\nabla V}{v}|_{J=0}$) using linear fitting slope coefficients ($\frac{\nabla V}{v}$). Due to the lack of resistivity measurements at LSBB (LSQi) and the presence of rainfall at the Larzac and LSBB sites during the selected periods, excess charge density calculations

were limited to the two trees at the Font-Blanche site, where sap flow consistently dominated SP generation. The resulting $\hat{Q}_v$ values for FBPh and FBQi in spring are 15.2 and 31.2 Cm⁻³, respectively. Calculating the seasonal variations of $\hat{Q}_v$ yields ranges of 6.8–25.7 Cm⁻³ for FBPh and 31.2–68.0 Cm⁻³ for FBQi (Fig. S8, Table S2). These values are consistent with those typically found in porous media (pH ranging from 6 to 8.5), such as sandstones, glass beads, and limestones (Jardani et al., 2007; Revil and Jardani, 2013). These materials correspond to typical pore sizes of 10–100 µm, which are of the same order of magnitude as the diameter of the cells conducting the sap (e.g. Sperry, 2003).
<figure><img src="image_9.png" /><figcaption>**Figure 7.** A period of 2-week SP data during the 2023 growing season for holm oak (FBQi) at the Font-Blanche site, processed using different approaches. (a) Raw SP data, compared with the sum of the last two decomposed models of SP signals; (b) raw sap velocity data, compared with the sum of the last two decomposed modes of SV signals; (c) fifth decomposed mode of SP signals, alongside SP calculated from the sum of the first to fifth decomposed modes of SV signals ($\sum_{i}^{5} VMD_{i}^{SV}$); and (d) SP data filtered using a finite impulse response (FIR) filter (with a bandpass of 1/48 to 1/1.5 h⁻¹ and a filter order of 100), compared with SP calculated from FIR-filtered SV data using the same filter parameters.</figcaption></figure>

# 5 Discussion

## 5.1 Hydraulic properties and tree SP

Water flow along the soil-plant continuum is governed by hydraulic conductivities and water potential gradients. While the water movement in sapwood differs from that in the vadose zone, certain functional similarities exist. For instance, the flow in xylem and soil can both be approximated by Poiseuille's law (e.g. Jougnot et al., 2019; Nobel, 2009) and exhibit similar responses of hydraulic conductivity to water potential gradients (Cai et al., 2022). Although xylem hydraulic conductivity does not fluctuate with minor variations in soil water potential, significant reductions in soil moisture can lead to declines in xylem conductivity, particularly under drought conditions (Brodribb and Hill, 2000; Carminati and Javaux, 2020). When soil water potential falls within a moderate range, xylem hydraulic conductivity remains relatively stable. However, at lower soil water potentials, xylem conductivity decreases (Cai et al., 2022; Kröber et al., 2014).

For the Font-Blanche site, the selected 2-week periods shown in Figs. 8 and S9 correspond to different seasonal water dynamics (Fig. S10). As indicated in Fig. S10b, soil water potential is lower in spring and summer compared to winter and autumn. Notably, the slopes (ratios) of $VMD_{5}^{SP}$ to $VMD_{5}^{SV}$ peak during summer, with similar but lower values in winter and autumn. An empirical relationship between the hydraulic permeability $k$ and $\hat{Q}_{v}$ for porous materials is given by Jardani et al. (2007):

$$ \log_{10}(\hat{Q}_{v}) = -9.23 - 0.82\log_{10}k. \qquad (12) $$

This relationship is also reproduced by a simple, physically based model proposed by Guarracino and Jougnot (2018), which considers a capillary as a first approximation for a xylem vessel. It is seen that decreasing permeability leads to an increasing effective excess charge density, subsequently reducing electrical streaming current density (Eq. 3) and weakening the electrokinetic effect. Based on Eq. (12) and the $\hat{Q}_{v}$ estimates introduced in Sect. 4.3.2,
<figure><img src="image_10.png" /><figcaption>**Figure 8.** The relationship between the fifth decomposed modes of sap velocity (VMD<sub>5</sub><sup>SV</sup>) and tree SP signals (VMD<sub>5</sub><sup>SP</sup>) in (a, b) winter, (c, d) summer, and (e, f) autumn of 2023 for the Aleppo pine (FBPh) at the Font-Blanche site. (a, c, e) Scatter plots of the decomposed modes of sap velocity and SP signals, with black lines indicating linear regression results and red dots indicating points (outliers) outside the 95 % confidence level of the regression. (b, d, f) Corresponding decompositions of SP signals (orange dots) and calculated SP (black lines), based on the linear relationship between VMD<sub>5</sub><sup>SV</sup> and VMD<sub>5</sub><sup>SP</sup>, with red dots as outliers outside the 95 % confidence level of the linear regression.</figcaption></figure>

we calculated the sapwood permeabilities for Aleppo pine (FBPh) as 4.49 × 10<sup>-13</sup>, 2.01 × 10<sup>-13</sup>, 1.06 × 10<sup>-13</sup>, and 5.38 × 10<sup>-13</sup> m<sup>2</sup> in winter, spring, summer, and autumn, respectively. Detailed seasonal characteristics of 2-week data for the Aleppo pine (FBPh) and holm oak (FBQi) at the Font-Blanche site can be found in Table S2.

The ionic content and pH of xylem sap are reported to con-
tribute to differences in hydraulic properties (López-Portillo
et al., 2005; Losso et al., 2023). Referring to Eqs. (2)–
(4) in Sect. 2.1 and assuming the electrical conductivity of
xylem sap σ<sub>f</sub> is 270 μS cm<sup>-1</sup> (Losso et al., 2023) and the di-
electric permittivity of xylem sap ε<sub>f</sub> is 7.10 × 10<sup>-10</sup> F m<sup>-1</sup>
(equivalent to water at 20 °C), the corresponding zeta po-
tentials for FBPh are estimated to be −49.5, −42.8, −38.8,
and −51.1 mV across the four seasonal periods, respectively.
These zeta potentials are of similar magnitude to typical val-
ues at the surface of silica grains in contact with an elec-
trolyte, with concentration in the range of several millimoles
per litre of NaCl (Revil et al., 1999). Although a lower sap
rate is observed, a stronger electrokinetic coupling capacity

in a tree still induces a larger-amplitude SP. This suggests
that the electrokinetic effects generated by sap flow might
have similar characteristics to porous media. Thus, it has po-
tential for short-term predictions using tree SP to evaluate
sap velocity.

While the electrokinetic properties of plant tissues differ
from those of mineral media due to their organic composi-
tion, similarities in surface charge behaviour and ion adsorp-
tion justify the preliminary use of porous media models to
estimate SP and SV in plants. Plant tissues, such as xylem
cell walls composed of cellulose, hemicellulose, lignin, and
pectin, exhibit negatively charged surfaces of plasma mem-
brane and carboxyl groups in pectin, allowing for cation
binding. This behaviour is comparable to negatively charged
mineral surfaces with electrical double layers. Studies have
reported zeta potentials in plant systems ranging from –48
to 23 mV, depending on pH and ionic conditions (Kinraide
et al., 1998), which, while smaller in magnitude, remain com-
parable to those observed in silica grains. Although direct
measurements of xylem zeta potential are limited, research
on lignocellulosic materials suggests they also exhibit negative zeta potentials under near-neutral pH conditions (Hubbe, 2006). Despite these similarities, plant tissues fundamentally differ from mineral media, with their unique composition influencing water and ion transport. Experimental studies are essential to directly measure xylem zeta potential and refine models for plant electrokinetic processes. This would improve the accuracy of SP and SV estimates and advance our understanding of electrokinetic phenomena in vascular plants.

As discussed in Sect. 2.1, several factors affect or alter tree SP. However, we contend that the gravitational forces driven by lunar-solar tides do not dominate tree SP generation as suggested in previous studies (Barlow, 2012; Le Mouël et al., 2024). First, rainfall events do not change gravitational forces, yet tree SP was clearly affected by rainfall (e.g. Fig. 6f and h). Second, a strong link between diurnal tree SP and sap velocity is observed during the growing season without rainfall (e.g. Figs. 4 and 6a, b). While sap flow may not be the sole mechanism for tree SP generation, our results suggest a close relationship between tree SP and sap velocity. Further studies with monitoring of electrical resistivity, solute concentration, and water content in sapwood are needed to elucidate the mechanisms behind tree SP generation and variation.

## 5.2 Electrochemical effects

The results presented in Sect. 4.3 indicate that electrokinetic mechanisms alone fail to fully explain the amplitudes of tree SP and the phase shifts between the diurnal variations of the measured SP and predicted SP from measured sap velocity data. Researchers propose that pH gradients contribute to SP differences (e.g. Gil and Vargas, 2023; Love et al., 2008). According to the Nernst potential, the following applies:

$$ \varphi_{\text{pH}}^{\text{lower}} - \varphi_{\text{pH}}^{\text{upper}} = -\frac{k_B T_K}{e_0 q_i} (\text{pH}^{\text{lower}} - \text{pH}^{\text{upper}}), \quad (13) $$

where the electrochemical effect could dominate SP generation under significant pH contrasts between the measured electrodes; however, the electrokinetic effect associated with transpiration would still persist.

Schill et al. (1996) observed decreasing pH, ranging from 6.6 to 7.4 and an increase in osmotic potential with height in maple trees, with a pH gradient of approximately 0.8 m⁻¹ in the lower trunk. Applying Eq. (13) to a representative pH difference of ΔpH = pH<sup>lower</sup> - pH<sup>upper</sup> = 0.1 over 0.1 m yields a calculated potential difference of -2.53 mV at 20 °C. Such pH-induced effects likely contribute to the underestimation of SP by sap velocity data (Fig. 6d, f, and h). Furthermore, upward longitudinal concentration gradients of alkaline earth cations (K⁺, Ca²⁺, Mg²⁺) in some tree stems (McDonald et al., 2002) suggest complex sap chemistry, which may partly explain the non-linear relationship between SP and sap velocity observed in our study.

## 5.3 Experimental setup concerning tree SP

In our measurements, stainless-steel electrodes were employed to measure trunk SP. Compared to other materials such as copper or aluminium, stainless-steel electrodes demonstrate relatively stable performance (Hao et al., 2015). However, it is essential to note that metal electrodes are sensitive to environmental changes and are at risk of being polarized. Although the distance between the tree electrodes is small, monitoring the temperature and adding two or three electrodes around each measured point can allow corrections for electrode drift induced by electrode-related effects, enabling a more accurate determination of the tree SP related to physiological changes.

As outlined earlier, all analyses conducted in Sect. 4 were performed without considering electrode effects. In future research, we suggest an improved experimental setup that involves distributing electrodes along the trunk with greater spacing (e.g. Fig. S11) and using duplicated measurements to ensure that the conclusions are robust. By enhancing the experimental setup, we will further strengthen our current understanding of the relationship between tree SP and the transpiration process.

Furthermore, the relationship between SP and upward sap flow does not systematically show a linear pattern. To explore the spatial characteristics associated with physiological processes, further observations and modelling are required. We recommend several directions for future research. First, expanding measurements to include multiple trees within the same site would improve the statistical robustness of findings and help clarify species-specific effects from general physiological principles. Second, utilizing ion-selective electrodes can be used to assess the contributions of different ions to SP observations. Such measurements, combined with detailed analyses of xylem sap chemistry, would enable a more comprehensive understanding of the electrochemical effects driving SP signals. Third, we propose that future studies adopt a coupled experimental and modelling approach. Physiochemical models incorporating ionic transport, pH gradients, and electrokinetic effects could help quantify the relative contributions of these processes to observed SP patterns. Such models would benefit from controlled experiments under varying environmental conditions (e.g. soil moisture, temperature, and nutrient availability) to validate their predictive capabilities. With these advancements, SP monitoring can potentially improve transpiration rate assessments, enable non-invasive studies of water transport dynamics, and contribute to our understanding of tree responses to environmental stressors.

## 6 Conclusions

This study presents a unique long-term dataset integrating SP and sap velocity measurements on four trees over a 1-
year period in a Mediterranean climate. Our findings reveal strong diurnal coherence between SP and sap velocity, which weakens with higher water availability and is accompanied by increased phase shifts. Different trees and the same tree across different seasons exhibit distinct SP responses, with varying diurnal SP: sap velocity ratios.

During dry seasons, the strong agreement between observed SP and SP predicted from sap velocity data underscores the electrokinetic dominance of electrokinetic effects. When electrokinetic mechanisms governed SP generation, we estimated the excess charge densities and zeta potentials of holm oak and Aleppo pine at the same site in different seasons using a linear regression model. The obtained values are comparable to the typical range for porous rock materials. Unlike point measurements of sap velocity, SP signals reflect an integrated bioelectrical process between the electrodes but are also influenced by electrode-related effects and environmental factors. These results highlight SP's potential as a tool for transpiration rate assessment in dry conditions. However, discrepancies during wet periods and non-linear relationships in some trees indicate additional physiological and electrochemical processes influencing SP signals.

While this study advances the understanding of tree SP dynamics, certain limitations – such as short electrode separations and the use of single electrode pairs per tree – highlight key areas for future research. Enhancing electrode materials, optimizing experimental designs, and incorporating physiochemical modelling are important for further deciphering the complex bioelectrical processes governing SP signals. These advancements will improve the reliability of SP measurements and broaden their applications in ecohydrology and plant physiology.

*Data availability.* All data used in this study are freely accessible and available on Zenodo (https://doi.org/10.5281/zenodo.12662288, Hu et al., 2024) and will also be published on the H+ National Observatory Service with free access.

*Supplement.* The supplement related to this article is available online at https://doi.org/10.5194/hess-29-2997-2025-supplement.

*Author contributions.* KH: conceptualization, methodology, formal analysis, data curation, the writing of the original draft, reviewing, editing, and visualization. BL: investigation, data curation, reviewing, and editing. SDC: investigation, the writing of the original draft, reviewing, and editing. NoL: conceptualization, reviewing, and editing. CC: investigation, reviewing, and editing. NKMS: supervision and data curation. NiL: conceptualization, methodology, the writing of the original draft, reviewing, editing, and supervision. DJ: conceptualization, methodology, the writing of the original draft, reviewing, editing, project administration, and supervision.

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

*Acknowledgements.* The authors thank the editor Anke Hildebrandt and the two anonymous reviewers for their insightful comments and constructive feedback, which enhanced the quality of this paper.

*Financial support.* This research has been partly supported by the OZCAR Research Infrastructure, the H+ National Observatory Service, the OSU OREME, the MOOSE project (an EC2CO initiative), the Postdoctoral Fellowship Program of the China Postdoctoral Science Foundation (grant nos. GZC20241598 and 2024M753016), and the “CUG Scholar” Scientific Research Funds at China University of Geosciences (Wuhan) (project no. 2023139).

*Review statement.* This paper was edited by Anke Hildebrandt and reviewed by two anonymous referees.

## References

Al Mahrouqi, D., Vinogradov, J., and Jackson, M. D.: Zeta potential of artificial and natural calcite in aqueous solution, Adv. Colloid Interfac., 240, 60–76, https://doi.org/10.1016/j.cis.2016.12.006, 2017.

Alfieri, J. G., Kustas, W. P., and Anderson, M. C.: A brief overview of approaches for measuring evapotranspiration, Agroclimatology: Linking Agriculture to Climate, 60, 109–127, https://doi.org/10.2134/agronmonogr60.2016.0034, 2018.

Anderegg, W. R., Klein, T., Bartlett, M., Sack, L., Pellegrini, A. F., Choa, B., and Jansen, S.: Meta-analysis reveals that hydraulic traits explain cross-species patterns of drought-induced tree mortality across the globe, P. Natl. Acad. Sci. USA, 113, 5024–5029, https://doi.org/10.1073/pnas.1525678113, 2016.

Bachofen, C., Poyatos, R., Flo, V., Martínez-Vilalta, J., Mencuccini, M., Granda, V., and Grossiord, C.: Stand structure of Central European forests matters more than climate for transpiration sensitivity to VPD, J. Appl. Ecol., 60, 886–897, https://doi.org/10.1111/1365-2664.14383, 2023.

Barlow, P. W.: Moon and cosmos: plant growth and plant bioelectricity, in: Plant Electrophysiology: Signaling and Responses, edited by: Volkov, A. G., Springer, Berlin, Heidelberg, Germany, 249–280, https://doi.org/10.1007/978-3-642-29110-4_10, 249–280, 2012.

Belashev, B.: Monitoring the electrical potential difference of pine tree, BIO Web Conf., 93, 01015, https://doi.org/10.1051/bioconf/20249301015, 2024.
Brodribb, T. J. and Hill, R. S.: Increases in water potential gradient reduce xylem conductivity in whole plants. Evidence from a low-pressure conductivity method, Plant Physiol., 123, 1021–1028, https://doi.org/10.1104/pp.123.3.1021, 2000.

Cai, G., Ahmed, M. A., Abdalla, M., and Carminati, A.: Root hydraulic phenotypes impacting water uptake in drying soils, Plant Cell Environ., 45, 650–663, https://doi.org/10.1111/pce.14259, 2022.

Carminati, A. and Javaux, M.: Soil rather than xylem vulnerability controls stomatal response to drought, Trends Plant Sci., 25, 868–880, https://doi.org/10.1016/j.tplants.2020.04.003, 2020.

Carrière, S., Chalikakis, K., Danquigny, C., Davi, H., Mazzilli, N., Ollivier, C., and Emblanch, C.: The role of porous matrix in water flow regulation within a karst unsaturated zone: an integrated hydrogeophysical approach, Hydrogeol. J., 24, 1905–1918, https://doi.org/10.1007/s10040-016-1425-8, 2016.

Carrière, S. D., Ruffault, J., Pimont, F., Doussan, C., Simioni, G., Chalikakis, K., Limousin, J.-M., Scotti, I., Courdier, F., Cakpo, C.-B., Davi, H., and Martin-StPaul, N. K.: Impact of local soil and subsoil conditions on inter-individual variations in tree responses to drought: insights from Electrical Resistivity Tomography, Sci. Total Environ., 698, 134247, https://doi.org/10.1016/j.scitotenv.2019.134247, 2020a.

Carrière, S. D., Ruffault, J., Cakpo, C. B., Olioso, A., Doussan, C., Simioni, G., Chalikakis, K., Patris, N., Davi, H., and Martin-StPaul, N. K.: Intra-specific variability in deep water extraction between trees growing on a Mediterranean karst, J. Hydrol., 590, 125428, https://doi.org/10.1016/j.jhydrol.2020.125428, 2020b.

Carrière, S. D., Martin-StPaul, N. K., Doussan, C., Courbet, F., Davi, H., and Simioni, G.: Electromagnetic induction is a fast and non-destructive approach to estimate the influence of subsurface heterogeneity on forest canopy structure, Water, 13, 3218, https://doi.org/10.3390/w13223218, 2021a.

Carrière, S. D., Loiseau, B., Champollion, C., Ollivier, C., Martin-StPaul, N. K., Lesparre, N., Olioso, A., Hinderer, J., and Jougnot, D.: First evidence of correlation between evapotranspiration and gravity at a daily time scale from two vertically spaced superconducting gravimeters, Geophys. Res. Lett., 48, e2021GL096579, https://doi.org/10.1029/2021GL096579, 2021b.

Damour, G., Simonneau, T., Cochard, H., and Urban, L.: An overview of models of stomatal conductance at the leaf level, Plant Cell Environ., 33, 1419–1438, https://doi.org/10.1111/j.1365-3040.2010.02181.x, 2010.

Davies, E.: Electrical Signals in Plants: Facts and Hypotheses, in: Plant Electrophysiology: Theory and Methods, edited by: Volkov, A. G., Springer, Berlin, Heidelberg, Germany, https://doi.org/10.1007/978-3-540-37843-3_17, 407–422, 2006.

Dix, M. J. and Aubrey, D. P.: Recalibrating best practices, challenges, and limitations of estimating tree transpiration via sap flow, Current Forestry Reports, 7, 31–37, https://doi.org/10.1007/s40725-021-00134-x, 2021.

Do, F. and Rocheteau, A.: Influence of natural temperature gradients on measurements of xylem sap flow with thermal dissipation probes. 1. Field observations and possible remedies, Tree Physiol., 22, 641–648, https://doi.org/10.1093/treephys/22.9.641, 2002.

Dragomiretskiy, K. and Zosso, D.: Variational mode decomposition, IEEE T. Signal Proces., 62, 531–544, https://doi.org/10.1109/TSP.2013.2288675, 2014.

Dumont, M. and Singha, K.: Geophysics as a hypothesis-testing tool for critical zone hydrogeology, WIREs Water, 11, e1732, https://doi.org/10.1002/wat2.1732, 2024.

Fensom, D. S.: The bio-electric potentials of plants and their functional significance: I. An electrokinetic theory of transport, Can. J. Bot., 35, 573–582, https://doi.org/10.1139/b57-047, 1957.

Fensom, D. S.: The bioelectric potentials of plants and their functional significance: V. Some daily and seasonal changes in the electrical potential and resistance of living trees, Can. J. Bot., 41, 831–851, https://doi.org/10.1139/b63-068, 1963.

Fisher, J. B., Melton, F., Middleton, E., Hain, C., Anderson, M., Allen, R., McCabe, M. F., Baldocchi, D., Townsend, P. A., Kilic, A., Tu, K., Miralles, D. D., Perret, J., Lagouarde, J.-P., Waliser, D., Purdy, A. J., French, A., Schimel, D., Famiglietti, J. S., Stephens, G., and Wood, E. F.: The future of evapotranspiration: Global requirements for ecosystem functioning, carbon and climate feedbacks, agricultural management, and water resources, Water Resour. Res., 53, 2618–2626, https://doi.org/10.1002/2016WR020175, 2017.

Flo, V., Martinez-Vilalta, J., Steppe, K., Schuldt, B., and Poyatos, R.: A synthesis of bias and uncertainty in sap flow methods, Agr. Forest Meteorol., 271, 362–374, https://doi.org/10.1016/j.agrformet.2019.03.012, 2019.

Fores, B., Champollion, C., Mainsant, G., Albaric, J., and Fort, A.: Monitoring saturation changes with ambient seismic noise and gravimetry in a karst environment, Vadose Zone J., 17, 1–12, https://doi.org/10.2136/vzj2017.09.0163, 2018.

Fromm, J. and Lautner, S.: Electrical signals and their physiological significance in plants, Plant Cell Environ., 30, 249–257, https://doi.org/10.1111/j.1365-3040.2006.01614.x, 2007.

Fromm, J. and Lautner, S.: Generation, transmission, and physiological effects of electrical signals in plants, in: Plant Electrophysiology: Signaling and Responses, edited by: Volkov, A. G., Springer, Berlin, Heidelberg, Germany, https://doi.org/10.1007/978-3-642-29110-4_8, 207–232, 2012.

Gaillardet, J., Braud, I., Hankard, F., Anquetin, S., Bour, O., Dorfliger, N., de Dreuzy, J. R., Galle, S., Galy, C., Gogo, S., Gourcy, L., Habets, F., Laggoun, F., Longuevergne, L., Le Borgne, T., Naaim-Bouvet, F., Nord, G., Simonneaux, V., Six, D., Tallec, T., Valentin, C., Abril, G., Allemand, P., Arènes, A., Arfib, B., Arnaud, L., Arnaud, N., Arnaud, P., Audry, S., Comte, V. B., Batiot, C., Battais, A., Bellot, H., Bernard, E., Bertrand, C., Bessière, H., Binet, S., Bodin, J., Bodin, X., Boithias, L., Bouchez, J., Boudevillain, B., Moussa, I. B., Branger, F., Braun, J. J., Brunet, P., Caceres, B., Calmels, D., Cappelaere, B., Celle-Jeanton, H., Chabaux, F., Chalikakis, K., Champollion, C., Copard, Y., Cotel, C., Davy, P., Deline, P., Delrieu, G., Demarty, J., Dessert, C., Dumont, M., Emblanch, C., Ezzahar, J., Estèves, M., Favier, V., Faucheux, M., Filizola, N., Flammarion, P., Floury, P., Fovet, O., Fournier, M., Francez, A. J., Gandois, L., Gascuel, C., Gayer, E., Genthon, C., Gérard, M. F., Gilbert, D., Gouttevin, I., Grippa, M., Gruau, G., Jardani, A., Jeanneau, L., Join, J. L., Jourde, H., Karbou, F., Labat, D., Lagadeuc, Y., Lajeunesse, E., Lastennet, R., Lavado, W., Lawin, E., Lebel, T., Le Bouteiller, C., Legout, C., Lejeune, Y., Le Meur, E., Le Moigne, N., Lions, J., Lucas, A., Malet, J. P., Marais-Sicre, C., Maréchal, J. C., Marlin, C., Martin, P., Martins, J., Martinez, J. M., Massei, N., Mauclerc, A., Mazzilli, N., Molénat, J., Moreira-Turcq, P., Mougin, E., Morin, S., Ngoupayou, J. N., Panthou, G., Peugeot, C.,
Picard, G., Pierret, M. C., Porel, G., Probst, A., Probst, J. L., Rabatel, A., Raclot, D., Ravanel, L., Rejiba, F., René, P., Ribolzi, O., Riotte, J., Rivière, A., Robain, H., Ruiz, L., Sanchez-Perez, J. M., Santini, W., Sauvage, S., Schoeneich, P., Seidel, J. L., Sekhar, M., Sengtaheuanghoung, O., Silvera, N., Steinmann, M., Soruco, A., Tallec, G., Thibert, E., Lao, D. V., Vincent, C., Viville, D., Wagnon, P., and Zitouna, R.: OZCAR: The French network of critical zone observatories, Vadose Zone J., 17, 1–24, https://doi.org/10.2136/vzj2018.04.0067, 2018.

Gibert, D., Le Mouël, J. L., Lambs, L., Nicollin, F., and Perrier, F.: Sap flow and daily electric potential variations in a tree trunk, Plant Sci., 171, 572–584, https://doi.org/10.1016/j.plantsci.2006.06.012, 2006.

Gielen, B., de Beeck, M. O., Loustau, D., Ceulemans, R., Jordan, A., and Papale, D.: Integrated carbon observation system (icos): An infrastructure to monitor the european greenhouse gas balance, in: Terrestrial Ecosystem Research Infrastructures, CRC Press, 505–520, ISBN 9781315368252, 2017.

Gil, P. and Vargas, A. I.: Stem electrical potential variations may aid in the early detection of drought stress in fruit-bearing trees, Int. J. Agric. Nat. Resour., 50, 116–129, https://doi.org/10.7764/ijanr.v50i3.2552, 2023.

Gindl, W., Löppert, H. G., and Wimmer, R.: Relationship between streaming potential and sap velocity in *Salix alba* L., Phyton-Ann. Rei Bot., 39, 217–224, 1999.

Girard, F., Vennetier, M., Guibal, F., Corona, C., Ouarmim, S., and Herrero, A.: *Pinus halepensis* Mill. crown development and fruiting declined with repeated drought in Mediterranean France, Eur. J. For. Res., 131, 919–931, https://doi.org/10.1007/s10342-011-0565-6, 2012.

Goulden, M. L. and Field, C. B.: Three methods for monitoring the gas exchange of individual tree canopies: ventilated-chamber, sap-flow and Penman-Monteith measurements on evergreen oaks, Funct. Ecol., 125–135, https://doi.org/10.2307/2390121, 1994.

Granier, A.: Evaluation of transpiration in a Douglas-fir stand by means of sap flow measurements, Tree Physiol., 3, 309–320, 1987.

Granier, A., Biron, P., Bréda, N., Pontailler, J. Y., and Saugier, B.: Transpiration of trees and forest stands: Short and long-term monitoring using sapflow methods, Glob. Change Biol., 2, 265–274, https://doi.org/10.1111/j.1365-2486.1996.tb00078.x, 1996.

Grinsted, A., Moore, J. C., and Jevrejeva, S.: Application of the cross wavelet transform and wavelet coherence to geophysical time series, Nonlin. Processes Geophys., 11, 561–566, https://doi.org/10.5194/npg-11-561-2004, 2004.

Guarracino, L. and Jougnot, D.: A physically based analytical model to describe effective excess charge for streaming potential generation in water saturated porous media, J. Geophys. Res.-Sol. Ea., 123, 52–65, https://doi.org/10.1002/2017JB014873, 2018.

Hao, Z., Wang, G., Li, W., Zhang, J., and Kan, J.: Effects of electrode material on the voltage of a tree-based energy generator, PLOS ONE, 10, e0136639, https://doi.org/10.1371/journal.pone.0136639, 2015.

Harmon, R. E., Barnard, H. R., Day-Lewis, F. D., Mao, D., and Singha K.: Exploring environmental factors that drive diel variations in tree water storage using wavelet analysis, Front. Water, 3, 682285, https://doi.org/10.3389/frwa.2021.682285, 2021.

Hermans, T., Goderniaux, P., Jougnot, D., Fleckenstein, J. H., Brunner, P., Nguyen, F., Linde, N., Huisman, J. A., Bour, O., Lopez Alvis, J., Hoffmann, R., Palacios, A., Cooke, A.-K., Pardo-Álvarez, Á., Blazevic, L., Pouladi, B., Haruzi, P., Fernandez Visentini, A., Nogueira, G. E. H., Tirado-Conde, J., Looms, M. C., Kenshilikova, M., Davy, P., and Le Borgne, T.: Advancing measurements and representations of subsurface heterogeneity and dynamic processes: towards 4D hydrogeology, Hydrol. Earth Syst. Sci., 27, 255–287, https://doi.org/10.5194/hess-27-255-2023, 2023.

Hedrich, R. and Schroeder, J. I.: The physiology of ion channels and electrogenic pumps in higher plants, Annu. Rev. Plant Biol., 40, 539–569, https://doi.org/10.1146/annurev.pp.40.060189.002543, 1989.

Hölttä, T., Vesala, T., Sevanto, S., Perämäki, M., and Nikinmaa, E.: Modeling xylem and phloem water flows in trees according to cohesion theory and Münch hypothesis, Trees, 20, 67–78, https://doi.org/10.1007/s00468-005-0014-6, 2006.

H+: Network of hydrogeological research sites, https://hplus.ore.fr/, last access: 26 June 2025.

Hu, K., Jougnot, D., Huang, Q., Looms, M. C., and Linde, N.: Advancing quantitative understanding of self-potential signatures in the critical zone through long-term monitoring, J. Hydrol., 585, 124771, https://doi.org/10.1016/j.jhydrol.2020.124771, 2020.

Hu, K., Loiseau, B., Carrière, S., Lesparre, N., Champollion, C., Martin-StPaul, N., Linde, N., and Jougnot, D.: Dataset for “Self-potential signals related to tree transpiration in a Mediterranean climate”, Zenodo [data set], https://doi.org/10.5281/zenodo.12662288, 2024.

Hu, K.: Illustration of a tree, BioRender, https://BioRender.com/kixmyvj (last access: 11 July 2025), 2025a.

Hu, K.: Base layout for vertically extracted xylem and phloem, BioRender, https://BioRender.com/0qwjmt9 (last access: 11 July 2025), 2025b.

Hu, K.: Cellulose structure, https://BioRender.com/76osafb (last access: 11 July 2025), 2025c.

Hubbard, C. G., West, L. J., Morris, K., Kulessa, B., Brookshaw, D., Lloyd, J. R., and Shaw, S.: In search of experimental evidence for the biogeobattery, J. Geophys. Res.-Biogeo., 116, G04018, https://doi.org/10.1029/2011JG001713, 2011.

Hubbe, M. A.: Sensing the electrokinetic potential of cellulosic fiber surfaces, BioResources, 1, 116–149, 2006.

Islam, M., Janssen, D., Chao, D., Gu, J., Eisen, D., and Choa, F.: Electricity derived from plants, Journal of Energy and Power Engineering, 11, 614–619, 2017.

Jacob, T., Bayer, R., Chery, J., Jourde, H., Le Moigne, N., Boy, J. P., Hinderer, J., Luck, B., and Brunet, P.: Absolute gravity monitoring of water storage variation in a karst aquifer on the Larzac plateau (Southern France), J. Hydrol., 359, 105–117, https://doi.org/10.1016/j.jhydrol.2008.06.020, 2008.

Jardani, A., Revil, A., Boleve, A., Crespy, A., Dupont, J. P., Barrash, W., and Malama, B.: Tomography of the Darcy velocity from self-potential measurements, Geophys. Res. Lett., 34, L24403, https://doi.org/10.1029/2007GL031907, 2007.

Jasechko, S., Sharp, Z. D., Gibson, J. J., Birks, S. J., Yi, Y., and Fawcett, P. J.: Terrestrial water fluxes dominated by transpiration, Nature, 496, 347–350, https://doi.org/10.1038/nature11983, 2013.
Jayawickreme, D. H., Jobbágy, E. G., and Jackson, R. B.: Geophysical subsurface imaging for ecological applications, New Phytol., 201, 1170–1175, https://doi.org/10.1111/nph.12619, 2014.

Jougnot, D. and Linde, N.: Self-potentials in partially saturated media: The importance of explicit modeling of electrode effects, Vadose Zone J., 12, 1–21, https://doi.org/10.2136/vzj2012.0169, 2013.

Jougnot, D., Linde, N., Revil, A., and Doussan, C.: Derivation of soil-specific streaming potential electrical parameters from hydrodynamic characteristics of partially saturated soils, Vadose Zone J., 11, https://doi.org/10.2136/vzj2011.0086, 2012.

Jougnot, D., Linde, N., Haarder, E. B., and Looms, M. C.: Monitoring of saline tracer movement at vertically distributed self-potential measurements at the HOBE agricultural test site, Voulund, Denmark, J. Hydrol., 521, 314–327, https://doi.org/10.1016/j.jhydrol.2014.11.041, 2015.

Jougnot, D., Mendieta, A., Leroy, P., and Maineult, A.: Exploring the effect of the pore size distribution on the streaming potential generation in saturated porous media, insight from pore network simulations, J. Geophys. Res.-Sol. Ea., 124, 5315–5335, https://doi.org/10.1029/2018JB017240, 2019.

Jougnot, D., Roubinet, D., Guarracino, L., and Maineult, A.: Modeling Streaming Potential in Porous and Fractured Media, Description and Benefits of the Effective Excess Charge Density Approach, in: Advances in Modeling and Interpretation in Near Surface Geophysics, edited by: Biswas, A. and Sharma, S., Springer, Cham, Switzerland, https://doi.org/10.1007/978-3-030-28909-6_4, 61–96, 2020.

Kim, H. K., Park, J., and Hwang, I.: Investigating water transport through the xylem network in vascular plants, J. Exp. Bot., 65, 1895–1904, https://doi.org/10.1093/jxb/eru075, 2014.

Kinraide, T. B., Yermiyahu, U., and Rytwo, G.: Computation of surface electrical potentials of plant cell membranes: correspondence to published zeta potentials from diverse plant sources, Plant Physiol., 118, 505–512, https://doi.org/10.1104/pp.118.2.505, 1998.

Koppán, A., Fenyvesi, A., Szarka, L., and Westergom, V.: Measurement of electric potential difference on trees, Acta Biol. Szeged., 46, 37–38, https://abs.bibl.u-szeged.hu/index.php/abs/article/view/2232 (last access: 26 June 2025), 2002.

Kormiltsev, V. V., Ratushnyak, A. N., and Shapiro, V. A.: Three-dimensional modeling of electric and magnetic fields induced by the fluid flow movement in porous media, Phys. Earth Planet. In., 105, 109–118, https://doi.org/10.1016/S0031-9201(97)00116-7, 1998.

Kröber, W., Zhang, S., Ehmig, M., and Bruelheide, H.: Linking xylem hydraulic conductivity and vulnerability to the leaf economics spectrum – a cross-species study of 39 evergreen and deciduous broadleaved subtropical tree species, PLOS ONE, 9, e109211, https://doi.org/10.1371/journal.pone.0109211, 2014.

Kuang, X., Liu, J., Scanlon, B. R., Jiao, J. J., Jasechko, S., Lancia, M., Biskaborn, B. K., Wada, Y., Li, H., Zeng, Z., Guo, Z., Yao, Y., Gleeson, T., Nicot, J.-P., Luo, X., Zou, Y., and Zheng, C.: The changing nature of groundwater in the global water cycle, Science, 383, eadf0630, https://doi.org/10.1126/science.adf0630, 2024.

Kume, T., Onozawa, Y., Komatsu, H., Tsuruta, K., Shinohara, Y., Umebayashi, T., and Otsuki, K.: Stand-scale transpiration estimates in a Moso bamboo forest: (I) Applicability of

sap flux measurements, Forest Ecol. Manag., 260, 1287–1294, https://doi.org/10.1016/j.foreco.2010.07.012, 2010.

Kunert, N., Aparecido, L. M. T., Wolff, S., Higuchi, N., dos Santos, J., de Araujo, A. C., and Trumbore, S.: A revised hydrological model for the Central Amazon: The importance of emergent canopy trees in the forest water budget, Agr. Forest Meteorol., 239, 47–57, https://doi.org/10.1016/j.agrformet.2017.03.002, 2017.

Kurpius, M. R., Panek, J. A., Nikolov, N. T., McKay, M., and Goldstein, A. H.: Partitioning of water flux in a Sierra Nevada ponderosa pine plantation, Agr. Forest Meteorol., 117, 173–192, https://doi.org/10.1016/S0168-1923(03)00062-5, 2003.

Le Mouël, J. L., Gibert, D., and Poirier, J. P.: On transient electric potential variations in a standing tree and atmospheric electricity, CR Géosci., 342, 95–99, https://doi.org/10.1016/j.crte.2009.12.001, 2010.

Le Mouël, J. L., Gibert, D., Boulé, J. B., Zuddas, P., Courtillot, V., Lopes, F., Gèze, M., and Maineult, A.: On the effect of the luni-solar gravitational attraction on trees, arXiv [preprint], https://doi.org/10.48550/arXiv.2402.07766, 2024.

Leroy, P., Devau, N., Revil, A., and Bizi, M.: Influence of surface conductivity on the apparent zeta potential of amorphous silica nanoparticles, J. Colloid Interf. Sci., 410, 81–93, https://doi.org/10.1016/j.jcis.2013.08.012, 2013.

Linde, N., Revil, A., Boleve, A., Dagès, C., Castermant, J., Suski, B., and Voltz, M.: Estimation of the water table throughout a catchment using self-potential and piezometric data in a Bayesian framework, J. Hydrol., 334, 88–98. https://doi.org/10.1016/j.jhydrol.2006.09.027, 2007.

Linde, N., Doetsch, J., Jougnot, D., Genoni, O., Dürst, Y., Minsley, B. J., Vogt, T., Pasquale, N., and Luster, J.: Self-potential investigations of a gravel bar in a restored river corridor, Hydrol. Earth Syst. Sci., 15, 729–742, https://doi.org/10.5194/hess-15-729-2011, 2011.

Loiseau, B., Carrière, S. D., Jougnot, D., Singha, K., Mary, B., Delpierre, N., Guérin, R., and Martin-StPaul, N. K.: The geophysical toolbox applied to forest ecosystems–A review, Sci. Total Environ., 899, 165503, https://doi.org/10.1016/j.scitotenv.2023.165503, 2023.

López-Portillo, J., Ewers, F. W., and Angeles, G.: Sap salinity effects on xylem conductivity in two mangrove species, Plant Cell Environ., 28, 1285–1292, https://doi.org/10.1111/j.1365-3040.2005.01366.x, 2005.

Losso, A., Gauthey, A., Choat, B., and Mayr, S.: Seasonal variation in the xylem sap composition of six Australian trees and shrubs, AoB Plants, 15, plad064, https://doi.org/10.1093/aobpla/plad064, 2023.

Love, C. J., Zhang, S., and Mershin, A.: Source of sustained voltage difference between the xylem of a potted Ficus benjamina tree and its soil, PloS One, 3, e2963, https://doi.org/10.1371/journal.pone.0002963, 2008.

Luo, Z., Deng, Z., Singha, K., Zhang, X., Liu, N., Zhou, Y., He, X., and Guan, H.: Temporal and spatial variation in water content within living tree stems determined by electrical resistivity tomography, Agr. Forest Meteorol., 291, 108058, https://doi.org/10.1016/j.agrformet.2020.108058, 2020.

McCulloh, K. A., Winter, K., Meinzer, F. C., Garcia, M., Aranda, J., and Lachenbruch, B.: A comparison of daily water use estimates derived from constant-heat sap-flow probe values and
gravimetric measurements in pot-grown saplings, Tree Physiol., 27, 1355-1360, https://doi.org/10.1093/treephys/27.9.1355, 2007.

McDonald, K. C., Zimmermann, R., and Kimball, J. S.: Diurnal and spatial variation of xylem dielectric constant in Norway spruce (*Picea abies* [L.] Karst.) as related to microclimate, xylem sap flow, and xylem chemistry, IEEE T. Geosci. Remote, 40, 2063–2082, https://doi.org/10.1109/TGRS.2002.803737, 2002.

Miller, A. J. and Wells, D. M.: Electrochemical Methods and Measuring Transmembrane Ion Gradients, in: Plant Electrophysiology: Theory and Methods, edited by: Volkov, A. G., Springer, Berlin, Heidelberg, Germany, https://doi.org/10.1007/978-3-540-37843-3_2, 15-34, 2006.

Moreno, M., Simioni, G., Cailleret, M., Ruffault, J., Badel, E., Carrière, S., Davi, H., Gavinet, J., Huc, R., Limousin, J.-M., Marloie, O., Martin, L., Rodríguez-Calcerrada, J., Vennetier, M., and Martin-StPaul, N.: Consistently lower sap velocity and growth over nine years of rainfall exclusion in a Mediterranean mixed pine-oak forest, Agr. Forest Meteorol., 308, 108472, https://doi.org/10.1016/j.agrformet.2021.108472, 2021.

Mottl, O., Flantua, S. G., Bhatta, K. P., Felde, V. A., Giesecke, T., Goring, S., Grimm, E. C., Haberle, S., Hooghiemstra, H., Ivory, S., Kuneš, P., Wolters, S., Seddon, A. W. R., and Williams, J. W.: Global acceleration in rates of vegetation change over the past 18,000 years, Science, 372, 860–864, https://doi.org/10.1126/science.abg1685, 2021.

Nardini A., Salleo S., and Jansen S.: More than just a vulnerable pipeline: xylem physiology in the light of ion-mediated regulation of plant water transport, J. Exp. Bot., 62, 4701–4718, https://doi.org/10.1093/jxb/err208, 2011.

Nobel, P. S. (Ed.): Physicochemical and Environmental Plant Physiology, Academic Press, ISBN 978-0-12-374143-1, 2009.

Nolan, C., Overpeck, J. T., Allen, J. R., Anderson, P. M., Betancourt, J. L., Binney, H. A., Brewer, S., Bush, M. B., Chase, B. M., Cheddadi, R., Djamali, M., Dodson, J., Edwards, M. E., Gosling, W. D., Haberle, S., Hotchkiss, S. C., Huntley, B., Ivory, S. J., Kershaw, A. P., Kim, S. H., Latorre, C., Leydet, M., Lezine, A. M., Liu, K. B., Liu, Y., Lozhkin, A. V., McGlone, M. S., Marchant, R. A., Momohara, A., Moreno, P. I., Muller, S., Otto-Bliesner, B. L., Shen, C.; Stevenson, J., Takahara, H., Tarasov, P. E., Tipton, J., Vincens, A., Weng, C., Xu, Q., Zheng, Z., and Jackson, S. T.: Past and future global transformation of terrestrial ecosystems under climate change, Science, 361, 920–923, https://doi.org/10.1126/science.aan5360, 2018.

Oki, T. and Kanae, S.: Global hydrological cycles and world water resources, Science, 313, 1068–1072, https://doi.org/10.1126/science.1128845, 2006.

Ollivier, C., Olioso, A., Carrière, S. D., Boulet, G., Chalikakis, K., Chanzy, A., Charlier, J.-B., Combemale, D., Davi, H., Emblanch, C., Marloie, O., Martin-StPaul, N., Mazzilli, N., Simioni, G., and Weiss, M.: An evapotranspiration model driven by remote sensing data for assessing groundwater resource in karst watershed, Sci. Total Environ., 781, 146706, https://doi.org/10.1016/j.scitotenv.2021.146706, 2021.

Oliveras, I. and Llorens, P.: Medium-term sap flux monitoring in a Scots pine stand: analysis of the operability of the heat dissipation method for hydrological purposes, Tree Physiol., 21, 473–480, https://doi.org/10.1093/treephys/21.7.473, 2001.

Oyarce, P. and Gurovich, L.: Electrical signals in avocado trees: Responses to light and water availability conditions, Plant Signal. Behav., 5, 34–41, 2010.

OZCAR: French network of Critical Zone Observatories: Research and Applications, https://www.ozcar-ri.org/, last access: 26 June 2025.

Perrier, F. E., Petiau, G., Clerc, G., Bogorodsky, V., Erkul, E., Jouniaux, L., Lesmes, D., Mancnae, J., Meunier, J., Morgan, D., Nascimento, D., Oettinger, G., Schwarz, G., Toh, H., Valiant, M., Vozoff, K., and Yazici-Cakin, O.: A one-year systematic study of electrodes for long period measurements of the electric field in geophysical environments, J. Geomagn. Geoelectr., 49, 1677–1696, https://doi.org/10.5636/jgg.49.1677, 1997.

Petiau, G.: Second generation of lead-lead chloride electrodes for geophysical applications, Pure Appl. Geophys., 157, 357–382, https://doi.org/10.1007/s000240050004, 2000.

Poyatos, R., Granda, V., Molowny-Horas, R., Mencuccini, M., Steppe, K., and Martínez-Vilalta, J.: SAPFLUXNET: towards a global database of sap flow measurements, Tree Physiol., 36, 1449–1455, https://doi.org/10.1093/treephys/tpw110, 2016.

Poyatos, R., Granda, V., Flo, V., Adams, M. A., Adorján, B., Aguadé, D., Aidar, M. P. M., Allen, S., Alvarado-Barrientos, M. S., Anderson-Teixeira, K. J., Aparecido, L. M., Arain, M. A., Aranda, I., Asbjornsen, H., Baxter, R., Beamesderfer, E., Berry, Z. C., Berveiller, D., Blakely, B., Boggs, J., Bohrer, G., Bolstad, P. V., Bonal, D., Bracho, R., Brito, P., Brodeur, J., Casanoves, F., Chave, J., Chen, H., Cisneros, C., Clark, K., Cremonese, E., Dang, H., David, J. S., David, T. S., Delpierre, N., Desai, A. R., Do, F. C., Dohnal, M., Domec, J.-C., Dzikiti, S., Edgar, C., Eichstaedt, R., El-Madany, T. S., Elbers, J., Eller, C. B., Euskirchen, E. S., Ewers, B., Fonti, P., Forner, A., Forrester, D. I., Freitas, H. C., Galvagno, M., Garcia-Tejera, O., Ghimire, C. P., Gimeno, T. E., Grace, J., Granier, A., Griebel, A., Guangyu, Y., Gush, M. B., Hanson, P. J., Hasselquist, N. J., Heinrich, I., Hernandez-Santana, V., Herrmann, V., Hölttä, T., Holwerda, F., Irvine, J., Isarangkool Na Ayutthaya, S., Jarvis, P. G., Jochheim, H., Joly, C. A., Kaplick, J., Kim, H. S., Klemedtsson, L., Kropp, H., Lagergren, F., Lane, P., Lang, P., Lapenas, A., Lechuga, V., Lee, M., Leuschner, C., Limousin, J.-M., Linares, J. C., Linderson, M.-L., Lindroth, A., Llorens, P., López-Bernal, Á., Loranty, M. M., Lüttschwager, D., Macinnis-Ng, C., Maréchaux, I., Martin, T. A., Matheny, A., McDowell, N., McMahon, S., Meir, P., Mészáros, I., Migliavacca, M., Mitchell, P., Mölder, M., Montagnani, L., Moore, G. W., Nakada, R., Niu, F., Nolan, R. H., Norby, R., Novick, K., Oberhuber, W., Obojes, N., Oishi, A. C., Oliveira, R. S., Oren, R., Ourcival, J.-M., Paljakka, T., Perez-Priego, O., Peri, P. L., Peters, R. L., Pfautsch, S., Pockman, W. T., Preisler, Y., Rascher, K., Robinson, G., Rocha, H., Rocheteau, A., Röll, A., Rosado, B. H. P., Rowland, L., Rubtsov, A. V., Sabaté, S., Salmon, Y., Salomón, R. L., Sánchez-Costa, E., Schäfer, K. V. R., Schuldt, B., Shashkin, A., Stahl, C., Stojanović, M., Suárez, J. C., Sun, G., Szatniewska, J., Tatarinov, F., Tesař, M., Thomas, F. M., Tor-ngern, P., Urban, J., Valladares, F., van der Tol, C., van Meerveld, I., Varlagin, A., Voigt, H., Warren, J., Werner, C., Werner, W., Wieser, G., Wingate, L., Wullschleger, S., Yi, K., Zweifel, R., Steppe, K., Mencuccini, M., and Martínez-Vilalta, J.: Global transpiration data from sap flow measurements: the SAPFLUXNET database, Earth Syst. Sci. Data, 13, 2607–2649, https://doi.org/10.5194/essd-13-2607-2021, 2021.
Pozdnyakov, A. I.: Bioelectric potentials in the soil-plant system, Eurasian Soil Sci., 46, 742–750, https://doi.org/10.1134/S1064229313070089, 2013.

Revil, A. and Jardani, A. (Eds.): The Self-potential Method: Theory and Applications in Environmental Geosciences, Cambridge University Press, ISBN 978-1-107-01927-0, 2013.

Revil, A., Pezard, P. A., and Glover, P. W. J.: Stream- ing potential in porous media: 1. Theory of the zeta potential, J. Geophys. Res.-Sol. Ea., 104, 20021–20031, https://doi.org/10.1029/1999JB900089, 1999.

Scanlon, T. M. and Kustas, W. P.: Partitioning evapotran- spiration using an eddy covariance-based technique: Im- proved assessment of soil moisture and land–atmosphere exchange dynamics, Vadose Zone J., 11, vzj2012.0025, https://doi.org/10.2136/vzj2012.0025, 2012.

Schill, V., Hartung, W., Orthen, B., and Weisenseel, M. H.: The xylem sap of maple (Acer platanoides) trees–sap obtained by a novel method shows changes with season and height, J. Exp. Bot., 47, 123–133, https://doi.org/10.1093/jxb/47.1.123, 1996.

Schlesinger, W. H. and Jasechko, S.: Transpiration in the global water cycle, Agr. Forest Meteorol., 189, 115–117, https://doi.org/10.1016/j.agrformet.2014.01.011, 2014.

Simioni, G., Marie, G., Davi, H., Martin-St Paul, N., and Huc, R.: Natural forest dynamics have more influence than climate change on the net ecosystem production of a mixed Mediterranean forest, Ecol. Modell., 416, 108921, https://doi.org/10.1016/j.ecolmodel.2019.108921, 2020.

Smith, D. M. and Allen, S. J.: Measurement of sap flow in plant stems, J. Exp. Bot., 47, 1833–1844, 1996.

Spanswick, R. M.: Electrogenic pumps, in: Plant Electrophysiolog- ogy: Theory and Methods, edited by: Volkov, A. G., Springer, Berlin, Heidelberg, Germany, https://doi.org/10.1007/978-3-540-37843-3_10, 221–246, 2006.

Sperry, J. S.: Evolution of water transport and xylem structure, Int. J. Plant Sci., 164, S115-S127, https://doi.org/10.1086/368398, 2003.

Tattar, T. A. and Blanchard, R. O.: Electrophysiological research in plant pathology, Annu. Rev. Phytopathol., 14, 309–325, https://doi.org/10.1146/annurev.py.14.090176.001521, 1976.

Torrence, C. and Compo, G. P.: A practical guide to wavelet analysis, B. Am. Meteorol. Soc., 79, 61–78, https://doi.org/10.1175/1520-0477(1998)079%3C0061:APGTWA%3E2.0.CO;2, 1998.

Valois, R., Galibert, P. Y., Guerin, R., and Plagnes, V.: Applica- tion of combined time-lapse seismic refraction and electrical re- sistivity tomography to the analysis of infiltration and dissolu- tion processes in the epikarst of the Causse du Larzac (France), Near Surf. Geophys., 14, 13–22, https://doi.org/10.3997/1873-0604.2015052, 2016.

van Bel, A. J.: Transport phloem: low profile, high impact, Plant Physiol., 131, 1509–1510, 2003.

Volkov, A. G. and Markin, V. S.: Phytosensors and Phytoactuators, in: Plant Electrophysiology: Signaling and Responses, edited by: Volkov, A. G., Springer, Berlin, Heidelberg, Germany, 173–206, https://doi.org/10.1007/978-3-642-29110-4_7, 2012.

Voytek, E. B., Barnard, H. R., Jougnot, D., and Singha, K.: Transpiration-and precipitation-induced subsurface water flow observed using the self-potential method, Hydrol. Process., 33, 1784–1801, https://doi.org/10.1002/hyp.13453, 2019.

Wang, J., Turner, N. C., Feng, H., Dyck, M., and He, H.: Heat tracer- based sap flow methods for tree transpiration measurements: a mini review and bibliometric analysis, J. Exp. Bot., 74, 723–742, https://doi.org/10.1093/jxb/erac424, 2023.

Wei, Z., Yoshimura, K., Wang, L., Miralles, D. G., Jasechko, S., and Lee, X.: Revisiting the contribution of transpiration to global ter- restrial evapotranspiration, Geophys. Res. Lett., 44, 2792–2801, https://doi.org/10.1002/2016GL072235, 2017.

Zapata, R., Oliver-Villanueva, J. V., Lemus-Zúñiga, L. G., Luzuriaga, J. E., Mateo Pla, M. A., and Urchueguía, J. F.: Evaluation of electrical signals in pine trees in a mediter- ranean forest ecosystem, Plant Signal. Behav., 15, 1795580, https://doi.org/10.1080/15592324.2020.1795580, 2020.

Zapata, R., Oliver-Villanueva, J. V., Lemus-Zúñiga, L. G., Fuente, D., Mateo Pla, M. A., Luzuriaga, J. E., and Moreno Es- teve, J. C.: Seasonal variations of electrical signals of Pi- nus halepensis Mill. in Mediterranean forests in dependence on climatic conditions, Plant Signal. Behav., 16, 1948744, https://doi.org/10.1080/15592324.2021.1948744, 2021.

Zhang, H., Simmonds, L. P., Morison, J. I., and Payne, D.: Estima- tion of transpiration by single trees: comparison of sap flow mea- surements with a combination equation, Agr. Forest Meteorol., 87, 155–169, https://doi.org/10.1016/S0168-1923(97)00017-8, 1997.