# From shifting cultivation to teak plantation: effect on overland flow and sediment yield in a montane tropical catchment

Olivier Ribolzi, O. Evrard, Sylvain Huon, Anneke de Rouw, Norbert Silvera, Keo Oudone Latsachack, Bounsamai Soulileuth, Irène Lefevre, Alain Pierret, Guillaume Lacombe, et al.

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

## To cite this version:

Olivier Ribolzi, O. Evrard, Sylvain Huon, Anneke de Rouw, Norbert Silvera, et al. From shifting cultivation to teak plantation: effect on overland flow and sediment yield in a montane tropical catchment. Scientific Reports, 2017, 7, pp.3987. 10.1038/s41598-017-04385-2. hal-02453748

HAL Id: hal-02453748

https://hal.science/hal-02453748

Submitted on 24 Jan 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>
SCIENTIFIC REPORTS

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

OPEN

# From shifting cultivation to teak plantation: effect on overland flow and sediment yield in a montane tropical catchment

Received: 15 December 2016

Accepted: 15 May 2017

Published online: 21 June 2017

Olivier Ribolzi<fnref n="1" />, Olivier Evrard<fnref n="1" />,², Sylvain Huon³, Anneke de Rouw⁴, Norbert Silvera⁵, Keo Oudone Latsachack⁶, Bounsamai Soulileuth⁶, Irène Lefèvre², Alain Pierret⁵,⁶, Guillaume Lacombe⁷, Oloth Sengtaheuanghoung<fnref n="8" /> & Christian Valentin⁵

Soil erosion supplies large quantities of sediments to rivers of Southeastern Asia. It reduces soil fertility of agro-ecosystems located on hillslopes, and it degrades, downstream, water resource quality and leads to the siltation of reservoirs. An increase in the surface area covered with commercial perennial monocultures such as teak plantations is currently observed at the expanse of traditional slash-and-burn cultivation systems in steep montane environments of these regions. The impacts of land-use change on the hydrological response and sediment yields have been investigated in a representative catchment of Laos monitored for 13 years. After the gradual conversion of rice-based shifting cultivation to teak plantation-based systems, overland flow contribution to stream flow increased from 16 to 31% and sediment yield raised from 98 to 609 Mg km⁻². This result is explained by the higher kinetic energy of raindrops falling from the canopy, the virtual absence of understorey vegetation cover to dissipate drop energy and the formation of an impermeable surface crust accelerating the formation and concentration of overland flow. The 25-to-50% lower ¹³⁷Cs activities measured in soils collected under mature teak plantations compared to soils under other land uses illustrate the severity of soil erosion processes occurring in teak plantations.

Soil erosion is exacerbated by the intensification of agriculture in montane steep cultivated catchments of Southeastern Asia, where it delivers significant quantities of sediment to rivers<fnref n="1" />. This excessive sediment supply leads to water reservoir siltation² and transports particle-bound contaminants towards downstream areas³⁻⁵. Among these pollutants, there are growing concerns regarding the transport of persistent pathogenic bacteria with sediment and their development in the streambed of tropical catchments⁶, where the use of unclean stream water is associated with a high number of fatal and debilitating diseases⁷.

The impact of land-use change on soil erosion has been investigated in various regions of Southeastern Asia<fnref n="8" />, although most studies were restricted to the plot and hillslope scales⁹,¹⁰. The main factors known to control soil erosion at the plot scale<fnref n="1" />, ¹¹ are rainfall, overland flow and surface characteristics (e.g. land use, vegetation cover, slope gradient, slope length, surface crusting). At this scale, inter-rill erosion processes dominate and mainly detach and mobilize soil surface particles. Then, at the hillslope scale, inter-rill erosion is combined with overland flow, which may concentrate and form rills and gullies, thereby removing and transporting additional particles originating from the soil subsurface¹²⁻¹⁴. Finally, at the catchment scale, the dominance of overland flow as the

<footnote n="1">¹Institut de Recherche pour le Développement (IRD), Géosciences Environnement Toulouse (GET), UMR 5563 (CNRS-UPS-IRD), Toulouse, France. ²Laboratoire des Sciences du Climat et de l'Environnement (LSCE/IPSL), Unité Mixte de Recherche 8212 (CEA-CNRS-UVSQ), Université Paris-Saclay, Gif-sur-Yvette, France. ³Université Pierre et Marie Curie (UPMC), UMR 7618 iEES (UPMC-CNRS-IRD-INRA-Université Paris 7-UPEC), Paris, France. ⁴Institut de Recherche pour le Développement (IRD), UMR 7618 iEES (UPMC-CNRS-IRD-INRA-Université Paris 7-UPEC), Paris, France. ⁵Institut de Recherche pour le Développement (IRD), iEES-Paris, Centre IRD France-Nord, Bondy, France. ⁶Institut de Recherche pour le Développement (IRD), Department of Agricultural Land Management (DALaM), Vientiane, Lao PDR. ⁷International Water Management Institute (IWMI), Southeast Asia Regional Office, Vientiane, Lao PDR.</footnote>

<footnote n="8">⁸Department of Agricultural Land Management (DALaM), Ministry of Agriculture and Forestry, Vientiane, Lao PDR. Correspondence and requests for materials should be addressed to O.R. (email: olivier.ribolzi@ird.fr)</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Monitoring devices and sampling sites in the Houay Pano catchment. This original map was created using Surfer 13.2.438 from Golden Software, LLC.</figcaption></figure>

main supply of sediment export at the outlet needs to be more thoroughly investigated, as this phenomenon is in competition with other processes (e.g. reinfiltration and sedimentation, subsurface flow contribution). Although the direct measurement of overland flow is possible at the plot and hillslope scales, it is difficult to carry out similar measurements at the catchment scale. At this larger scale, indirect tracer-based approaches provide an effective method to calculate the respective contributions of subsurface and surface flow. Typically, water mixing models using geochemical properties (e.g. electric conductivity, $\delta^{18}$O) are applied to quantify these contributions$^{15,16}$.

In forested catchments, subsurface contributions were shown to provide the dominant contribution to streamflow$^{17,18}$, whereas surface overland flow was identified as the main supply in cultivated catchments where an almost impermeable crust may be formed at the soil surface$^{15, 19, 20}$. Most land-use change studies investigated the conversion of forests into cropland, showing the associated increase in soil erosion$^{21, 22}$. To our knowledge, only few studies have been conducted on the conversion of annual crops into commercial tree monoculture plantations on the steep hillslopes of humid tropics although this rapid land-use change is affecting extensive areas across this region$^{23}$. The main objective associated with the reforestation of degraded land is often to restore the soil infiltration capacity and to limit soil erosion. However, it was shown that this reforestation may also increase plant water uptake and consequently reduce river baseflow$^{24,25}$.

To investigate the impact of land-use change on sediment yield, the monitoring of paired catchments is often conducted in nearby sites where contrasted land uses prevail while biophysical characteristics are as similar as possible. This synchronic approach suffers from several limitations because of the difficulty to select two perfectly identical catchments$^{26}$. In contrast, this approach is easier to conduct in similar zero-order headwater catchments (i.e. hollow valleys referred to as S7 and S8) with similar soil characteristics although different land uses. In addition to this synchronic approach, we applied a diachronic investigation in the larger catchment affected by land use change and draining these hollow valleys. The 0.6-km$^2$ Houay Pano catchment (Fig. 1), located in Northern Laos, is representative of areas that experienced large-scale afforestation with commercial teak tree (*Tectona grandis* L.f.) plantations at the expense of traditional slash-and-burn (i.e. shifting) cultivation. Discharge and sediment yield were monitored for 13 years in this experimental catchment.

The contributions of surface (i.e. overland flow) and sub-surface flow (i.e. groundwater) to river water at the catchment outlet were assessed by measuring the river water electrical conductivity. This technique is easier to implement than the more conventional method based on $\delta^{18}$O, and both techniques were found to provide similar results$^{15}$. The contributions of surface and sub-surface soil sources to stream sediment exports at the catchment outlet were discriminated by measuring the activity of Caesium-137 ($^{137}$Cs) emitted by thermonuclear bomb tests in the 1960s in the sediment transiting the river$^{27}$. This radioisotope characterised by a half-life of 30 years shows highly contrasting activities between the surface of cultivated soils exposed to atmospheric fallout and gully/ channel bank material sheltered from this fallout$^{28-30}$.
<figure><img src="image_5.png" /></figure>

**Figure 2.** Annual hydro-sedimentary and land-use changes in Houay Pano catchment: (a) Total rainfall depth (Rain), flood-triggering rainfall depth (Rain_F) and total stream flow depth (Runoff); (b) Area percentages of the catchment covered by the main land uses (i.e. annual crops; annual crops associated with one-or-two-years-old teak trees; teak plantations; fallow; secondary forest and other land uses including gardens and banana); (c) Box plots (25th and 75th percentiles and median) of infiltration rates computed for each flood event; (d) Box plots of overland flow contribution to total depth of flood events estimated using tracer-based hydrograph separation; (e) Annual sediment yield including suspended load and bed load contributions. This graph was created using Grapher 4.0 from Golden Software, LLC.

## Results

**Land-use change.** The proportion of secondary forests in the catchment decreased from 16% (2002) to 8% (2014) (Fig. 2b). During this period, the variations in the respective proportions of fallow (29–69%) and annual crops (4–42%) were negatively correlated and they showed large inter-annual variations, reflecting the management strategies of replacing fallows with annual crops through slash-and-burn, followed by the regrowth of vegetation in the temporarily abandoned cultivated fields. This entire cycle may take 2 to 6 years (Fig. 2b). Until 2007, the percentage of catchment area covered with teak plantations (Teak) remained low and did not exceed 4%. From 2008, it continuously increased, almost linearly, to reach a maximum of 36% of the catchment area in 2014 (Fig. 2b), illustrating the transition from an annual crop-based system towards an organisation dominated by tree plantations. Accordingly, the monitoring period was divided into a first sub-period (2002–2007) characterized by the alternation between upland rice (*Oryza sativa* L.), Job's tear (*Coix lacryma-jobi* L.) annual crops and bush fallow, in addition to sparse teak plantations, and a second sub-period (2008–2014) characterised by the concomitant expansion of teak plantations and the reduction of annual crops. During the last five years of this second sub-period, most of the teak plantations were more than 3 years old and characterised by limited understorey vegetation cover, which may enhance soil degradation. In 2014, land use greatly differed between the two hollow valleys investigated in details (Fig. 1b), with a dominant contribution of typical shifting cultivation landscape mosaics at S7 compared to a major contribution of teak plantations at S8.

**Rainfall and runoff.** During the 2002–2014 monitoring period of the Houay Pano catchment, annual rainfall depth (Rain) varied between 978 mm (2012) and 1884 mm (2011), with a mean of 1271 mm (SD: 271 mm, CV: 20%) (Fig. 2a). These values are similar to those recorded at the closest reference weather station located in Luang Prabang (58-year long-term mean: 1302 mm; SD: 364 mm; CV: 28%). Only 32% of the rainfall events generated stormflow (i.e. total stream flow during a flood that is calculated when summing surface and subsurface contributions). The mean annual depth corresponding to these flood-triggering rainfall events (Rain_F) was 404 mm (SD: 118 mm; CV: 28%). Annual runoff showed large inter-annual variations (SD: 52 mm; CV: 67%; mean: 77 mm; range: 13–174 mm) (Fig. 2a). In 2014, runoff from the hollow valleys strongly differed (Friedman non-parametric test for paired samples [α = 0.0001] p < 0.0001) between S7 (32 mm) and S8 (200 mm; Fig. 3b).
<figure><img src="image_6.png" /><figcaption>**Figure 3.** (a) Daily rainfall depth, (b) cumulated suspended sediment yield and, (c) overland flow measured in 2014 at the hydrological station installed at the outlet of two zero-order headwater subcatchments (i.e. ephemeral flow): S7 (dominance of traditional shifting cultivation mosaics); S8 (dominance of teak plantations). This graph was created using Excel 14.0 Software.</figcaption></figure>

**Infiltration rate, overland flow and sediment yield.** The annual median values of infiltration rates calculated for each flood event at the catchment scale (InfRate) strongly varied, between 4 and 25 mm h⁻¹ depending on years (Fig. 2c). It first increased from 9 mm h⁻¹ to 20 mm h⁻¹ during the first sub-period (i.e. 2002–2007), and then decreased gradually from 18 to 4 mm h⁻¹ during the second sub-period (i.e. 2008–2014). Although it is characterised by large inter-annual variations, the annual median of overland flow contribution to the total stream flow during floods (F<sub>OLF</sub>) was significantly higher during the last five years of the second study period (2010–2014; inter-annual median: 31%) when soils of the teak plantations were already degraded and characterised by sparse understorey vegetation cover (Mann-Whitney bilateral non-parametric test [α = 0.01] p = 0.004) compared to the first sub-period (2002–2009; inter-annual median: 16%) (Fig. 2d). Similarly, annual sediment yields (SY) exhibited large inter-annual fluctuations between 10 Mg km⁻² and 1,260 Mg km⁻², with a significant increase (Mann-Whitney bilateral non-parametric test [α = 0.05] p = 0.003) from the first (2002–2009; inter-median: 98 Mg km⁻²) to the second (2010–2014; inter-median: 609 Mg km⁻²) sub-periods (Fig. 2e). Bed load contribution to these total annual losses fluctuated between 9 and 42%, with similar averages during the first and second sub-periods (mean of 16% and 15%, respectively). In the hollow valleys in 2014, cumulated sediment yield remained negligible (Friedman non-parametric test for paired samples [α = 0.0001] p < 0.0001) at S7 (127 Mg km⁻²), whereas it was very high at S8 (2499 Mg km⁻²; Fig. 3c).

A correlation analysis was performed to investigate the relationships between the catchment hydrological behaviour, land uses and sediment yields. The main objective of this analysis was to evaluate how the transition from an annual crop-based system towards a tree plantation-based system affected the hydro-sedimentary response of the catchment. We evaluated the significance of the correlation between several variables reported in Table 1. Among these parameters, the areal percentage of more than 3-year-old teak tree monocultures (Teak₃) was considered. Teak₃ was used in order to exclude the youngest plantations which were not yet submitted to burning (Fig. 4a), and in which soils were not yet degraded.

Contrary to the total annual overland flow depth (OLF), which is correlated to Rain (p-value < 0.05), F<sub>OLF</sub> is not. In contrast, F<sub>OLF</sub> is positively correlated to Teak₃ (p-value < 0.05; Table 1) and negatively correlated to the sum of fallow and forest areas (Fw&Fo) (p-value < 0.05), and it provides therefore a better proxy of land use change. This is also confirmed by the regression analysis with a positive correlation with annual crops cultivated alone or in combination with one-or-two-years-old teak (AC) (Fig. 4a). InfRate is very significantly and negatively correlated to Teak and Teak₃ (p-value < 0.01) (Table 1).

Sediment yield was strongly correlated to the annual suspended load (SL) (p-value < 0.001), which is the main export route for soil particles (SL account for 80 and 78% of SY, on average, during the first and the second
<table><caption>**Table 1.** Correlation matrix (Pearson correlation coefficients; *p-value < 0.05, **p-value < 0.01, ***p-value < 0.001) between the following variables: annual rainfall depth (Rain); flood-triggering annual rainfall depth (Rain_F); annual overland flow depth (OLF); annual subsurface flow depth (SSF); annual median of overland flow contribution to total depth of each flood event (F<sub>OLF</sub>); annual median of subsurface flow contribution to total depth of each flood event (F<sub>SSF</sub>); annual median of event infiltration rates at the catchment scale (InfRate); annual suspended load (SL); annual bed load (BL); Total sediment yield (SY, i.e. SL + BL); areal percentage of teak tree monocultures (Teak); areal percentage of teak tree monocultures of more than 3-year-old (Teak<sub>&gt;3</sub>); areal percentage of annual crops cultivated alone or in combination with one-or-two-years-old teak (AC); areal percentage of fallow and secondary forest (Fw&Fo).</caption><thead><tr><th></th><th>Rain mm</th><th>Rain_F mm</th><th>OLF mm</th><th>SSF mm</th><th>F<sub>OLF</sub> %</th><th>F<sub>SSF</sub> %</th><th>InfRate mm h<sup>-1</sup></th><th>SL Mg km<sup>-2</sup></th><th>BL Mg km<sup>-2</sup></th><th>SY Mg km<sup>-2</sup></th><th>Teak %</th><th>Teak<sub>&gt;3</sub> %</th><th>AC %</th></tr></thead><tbody><tr><td>Rain</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rain_F</td><td>0.55</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>OLF</td><td><strong>0.69*</strong></td><td>0.49</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SSF</td><td>0.32</td><td>0.18</td><td>0.94</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>F<sub>OLF</sub></td><td>0.43</td><td>0.39</td><td>0.17</td><td>-0.45</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>F<sub>SSF</sub></td><td>-0.42</td><td>-0.39</td><td>-0.17</td><td>0.46</td><td>-1***</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>InfRate</td><td>-0.25</td><td>-0.55</td><td>0.02</td><td>0.28</td><td>-0.52</td><td>0.52</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SL</td><td>0.52</td><td><strong>0.83***</strong></td><td>0.26</td><td>-0.15</td><td><strong>0.62*</strong></td><td>-0.62*</td><td>-0.44</td><td>1</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>BL</td><td><strong>0.62*</strong></td><td><strong>0.75**</strong></td><td>0.45</td><td>-0.16</td><td><strong>0.81**</strong></td><td>-0.81**</td><td>-0.41</td><td><strong>0.82**</strong></td><td>1</td><td></td><td></td><td></td><td></td></tr><tr><td>SY</td><td>0.55</td><td><strong>0.84***</strong></td><td>0.29</td><td>-0.16</td><td><strong>0.66*</strong></td><td>-0.66*</td><td>-0.44</td><td><strong>0.99***</strong></td><td><strong>0.86***</strong></td><td>1</td><td></td><td></td><td></td></tr><tr><td>Teak</td><td>0.05</td><td>0.28</td><td>-0.37</td><td>-0.69*</td><td>0.55</td><td>0.55</td><td>-0.71**</td><td>0.40</td><td>0.39</td><td>0.40</td><td>1</td><td></td><td></td></tr><tr><td>Teak<sub>&gt;3</sub></td><td>0.30</td><td>0.35</td><td>-0.18</td><td>-0.60*</td><td><strong>0.73**</strong></td><td>-0.73**</td><td>-0.72**</td><td>0.49</td><td>0.55</td><td>0.51</td><td><strong>0.93***</strong></td><td>1</td><td></td></tr><tr><td>AC</td><td>0.14</td><td>0.17</td><td>0.37</td><td>0.07</td><td>0.29</td><td>-0.29</td><td>0.07</td><td>0.19</td><td>0.34</td><td>0.21</td><td>-0.40</td><td>-0.24</td><td>1</td></tr><tr><td>Fw&Fo</td><td>-0.18</td><td>-0.42</td><td>-0.05</td><td>0.49</td><td>-0.75**</td><td>-0.75**</td><td>0.52</td><td>-0.53</td><td>-0.68*</td><td>-0.56</td><td>-0.45</td><td>-0.55</td><td>-0.64*</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Statistical relationships (data collected at the S4 monitoring station) between: (a) median overland flow contribution to stream flow during storm events (F<sub>OLF</sub> %) and the two explanatory variables: areal percentage of teak monocultures during year n-1 (Teak<sub>&gt;3</sub> %) and areal percentage of annual crops cultivated alone or in combination with one-or-two-years-old teak trees (AC, %); (b) median sediment yield (SY, Mg km<sup>-2</sup>) and the two explaining variables F<sub>OLF</sub> and annual flood-triggering rainfall depth (Rain_F, mm). This graph was created using Excel 14.0 Software.</figcaption></figure>

sub-periods, respectively). Annual bed load (BL) is also well correlated to SL (*p-value < 0.01). Interestingly, SY (i.e. BL + SL) is correlated to Rain_F (*p-value < 0.001, < 0.01 and < 0.001 respectively; Table 1) and F<sub>OLF</sub> (*p-value < 0.05, 0.01 and 0.05 respectively), which agrees with the regression equation (Fig. 4b). A negative correlation was observed between the areal percentage of Fw&Fo and soil losses. However, it was significant for BL only (*p-value < 0.05) (Table 1).

Multiple log-linear and linear regressions were conducted to identify explanatory variables that best predicted the inter-annual variations of F<sub>OLF</sub> and SY. All the variables listed in Table 1 were considered in this analysis. A best sub-set of two explanatory variables was selected in each regression: Teak<sub>&gt;3</sub> and the areal percentage of annual crops cultivated alone or in combination with one-or-two-years-old teak (AC) for F<sub>OLF</sub> (Fig. 3a); Rain_F and F<sub>OLF</sub> for SY (Fig. 4b). While the regression model predicting SY was linear, a stronger log-linear relationship
<figure><img src="image_8.png" /><figcaption>**Figure 5.** $^{137}$Cs activities measured in (a) annual composite sediment samples in 2002 (SL-2002) and 2012 (SL-2012) at the catchment outlet (S4 monitoring station); (b) potential sources of sediment (surface of cultivated soils sampled between 2000–2005 [Surf-2002]; surface of cultivated soils sampled in 2012 [Surf-2012]; subsurface samples collected in 2012 in gullies and channel banks [Subsurf]; (c) Surface soil under teak plantations of different maturity levels: 5-years old plantation [Teak-5 yr]; 19-years old plantation [Teak-19 yr]; 25-years old plantation [Teak-25 yr]). This original graph was created using SigmaPlot 12 software.</figcaption></figure>

<table><caption>**Table 2.** Results of the Tukey-Kramer Honestly Significant Difference (HSV) test conducted on $^{137}$Cs activities measured in the soil and sediment sample groups shown in Fig. 5. SL = suspended sediment load, Surf = soil surface (<5 cm), Subsurf = soil subsurface (gullies and channel banks).</caption><thead><tr><th></th><th>SL-2002</th><th>Surf-2002</th><th>Surf-2012</th><th>Teak-5 yrs</th><th>Teak-19 yrs</th><th>Teak-25 yrs</th><th>Subsurf</th><th>SL-2012</th></tr></thead><tbody><tr><td>SL-2002</td><td></td><td>insignificant</td><td>insignificant</td><td>insignificant</td><td>*p &lt; 0.05</td><td>**p &lt; 0.01</td><td>**p &lt; 0.01</td><td>insignificant</td></tr><tr><td>Surf-2002</td><td>insignificant</td><td></td><td>*p &lt; 0.05</td><td></td><td>**p &lt; 0.01</td><td>**p &lt; 0.01</td><td>**p &lt; 0.01</td><td>*p &lt; 0.05</td></tr><tr><td>Surf-2012</td><td>insignificant</td><td>*p &lt; 0.05</td><td></td><td>insignificant</td><td>insignificant</td><td>*p &lt; 0.05</td><td>**p &lt; 0.01</td><td>insignificant</td></tr><tr><td>Teak-5 yrs</td><td>insignificant</td><td>**p &lt; 0.01</td><td>insignificant</td><td></td><td>insignificant</td><td>insignificant</td><td>**p &lt; 0.01</td><td>insignificant</td></tr><tr><td>Teak-19 yrs</td><td>*p &lt; 0.05</td><td>**p &lt; 0.01</td><td>insignificant</td><td></td><td>insignificant</td><td>insignificant</td><td>**p &lt; 0.01</td><td>insignificant</td></tr><tr><td>Teak-25 yrs</td><td>**p &lt; 0.01</td><td>**p &lt; 0.01</td><td>*p &lt; 0.05</td><td>insignificant</td><td>insignificant</td><td></td><td>**p &lt; 0.01</td><td>insignificant</td></tr><tr><td>Subsurf</td><td>**p &lt; 0.01</td><td>**p &lt; 0.01</td><td>**p &lt; 0.01</td><td>**p &lt; 0.01</td><td>**p &lt; 0.01</td><td>**p &lt; 0.01</td><td></td><td>insignificant</td></tr><tr><td>SL-2012</td><td>insignificant</td><td>*p &lt; 0.05</td><td>insignificant</td><td>insignificant</td><td>insignificant</td><td>insignificant</td><td>insignificant</td><td></td></tr></tbody></table>

was observed between $F_{OLF}$ and its explanatory variables, reflecting the typically skewed and non-normal statistical distributions of hydrological variables. $F_{OLF}$ and SY were statistically predicted using the equations presented in Fig. 4.

**Soil and sediment signatures.** The $^{137}$Cs signature in the composite sample of suspended sediment representative of the annual sediment export from the catchment decreased by more than 50% during the study period (from $2.5 \pm 0.4 \text{ Bq kg}^{-1}$ in 2002, to $1.2 \pm 0.4 \text{ Bq kg}^{-1}$ in 2012; Fig. 5a). In potential sediment sources, $^{137}$Cs activities varied between subsurface soils ($0.4 \pm 0.3 \text{ Bq kg}^{-1}$), surface soils collected in 2002 ($2.3 \pm 0.8 \text{ Bq kg}^{-1}$) and surface soils collected in 2012 ($1.8 \pm 0.4 \text{ Bq kg}^{-1}$; Fig. 5b). These values were used to discriminate material originating from the soil surface and sub-surface in the collected samples, with a high confidence ($p$-value < 0.01; Table 2). Despite a decrease in $^{137}$Cs activities recorded in cultivated soils between 2002 and 2012 ($p$-value < 0.05), which likely reflects the high intensity of erosion that occurred in cropland during the last decade, $^{137}$Cs activities found in surface soils collected in 2002 and 2012 remained significantly higher than those found in subsurface material ($p$-value < 0.01). Interestingly, surface soils under young (5-years old) teak plantations were characterized by activities ($1.6 \pm 0.3 \text{ Bq kg}^{-1}$) similar to those measured in surface soils in 2012 ($1.8 \pm 0.4 \text{ Bq kg}^{-1}$, Table 2) whereas surface soils collected under the most mature (19–25 years old) teak plantations showed significantly different signatures ($1.1 \pm 1.3 \pm 0.2-0.3 \text{ Bq kg}^{-1}$) compared to 2012 surface measurements ($1.8 \pm 0.4 \text{ Bq kg}^{-1}$, $p$-value < 0.05; Table 2) and subsurface ($0.4 \pm 0.3 \text{ Bq kg}^{-1}$, $p$-value < 0.01; Table 2) material, corresponding to an ‘intermediate’ signature between surface and subsurface sources.

Discussion

**Impact of land-use change on infiltration rates and overland flow.** In contrast to teak plantations⁹, fallows are characterised by higher infiltration rates associated with a higher biological activity maintaining greater porosity and lower surface crusting³¹. Several vegetation strata, including understorey, are observed in fallows. By dissipating the kinetic energy of raindrops, they reduce the soil surface crusting. In contrast, under teak canopies where large leaves concentrate rainfall, heavier rain drops fall from higher interception heights, impacting the soil with a greater kinetic energy³⁰, clogging the soil surface porosity and disturbing the soil surface structure, both processes contributing to crusting and the associated reduction in infiltration rates. Consequently, the production of overland flow increases and detaches more soil particles by sheet flow.

The contribution of surface overland flow to the total stream flow increased during the second sub-period. This evolution reflects a change in the hydrological behaviour of the catchment, which coincides with the increase
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Pictures taken in the Houay Pano catchment: (a) teak understorey in natural conditions (vegetation and litter are present); (b) teak understorey after a fire (litter and vegetation virtually cleared); (c) degraded soil surface with erosion crust and exposed teak roots; (d) sheet overland flow under a teak plantation during a rainfall event; (e) the hydrological station (so-called S4 station) equipped with a sediment trap for the quantification of bed load (larger blue PVC tube), an automatic sampler for the estimation of suspended solid load, a water level recorder and V-nutch weir to monitor the stream discharge, and a probe for measuring the electrical conductivity of stream water (in the smaller blue PVC tube).</figcaption></figure>

of the teak cover at the expense of fallow and annual crops. Until it is two-year old, teak is generally interspersed with annual crops (Fig. 2b). From the third year, if the natural understorey vegetation is maintained, it will protect the soil from the raindrop impact and the subsequent surface crusting⁹. However, in the study area, understo-rey vegetation is systematically removed and/or burnt by the farmers, and soil is left bare and exposed to rain-drop and overland flow (Fig. 6a,b and c). Consistently, F<sub>OLF</sub> increases in correlation with the increase in Teak<sub>&gt;3</sub> (Table 1). Nevertheless, during the four last years of the second sub-period, overland flow slightly decreased with the increase of Fw&Fo (Fig. 2b). These observations are consistent with plot-scale measurements conducted in the same catchment, showing that infiltration rates under fallow are higher than those found under teak plantations⁹.

**Impact of land-use change on sediment yields.** Sediment yields were significantly lower during the first sub-period (2002–2007) when the areal percentage of teak plantations was lower than 4%. However, a high SY value was observed in 2007 (833 Mg km⁻²). It was attributed to the particularly high areal percentage of annual crop observed during that specific year. Indeed, this land use is characterised by a low infiltrability due to the formation of soil surface crusting⁹, and a high sensitivity to overland flow and erosion¹. From 2011 onwards, when Teak<sub>&gt;3</sub> (i.e. teak with degraded soil surface conditions) exceeds 20% each year, high SY and F<sub>OLF</sub> values were observed (up to 1260 Mg km⁻² and 44%, respectively, Fig. 2e). These high values cannot be related to the areal percentage of annual crops, which decreased from 40 to 4% during the second sub-period. Rainfall has long been known as a major control factor for soil erosion at the plot scale³², possibly explaining the high SY values observed after 2011. Indeed, high annual rainfall depths were measured in 2011 (1884 mm) and 2013 (1524 mm). However, the SY value measured in 2004 (232 Mg km⁻²) is lower than that of 2012 (609 Mg km⁻²), despite a 23% higher rainfall depth (978 mm and 1305 mm, respectively) and a similar areal percentage of annual crops (29% and 24% respectively). However, the areal percentage of Teak<sub>&gt;3</sub> increased from 2 to 18% between 2004 and 2012. Therefore, the most likely factor that may explain this difference is the increase in the areal percentage of teak plantations. This is demonstrated by the different behaviour observed in contrasted sub-catchments (S7
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Graphical abstract: bottom left side, two pictures of the Houay Pano catchment in 2001 (almost entirely covered by fallow) and 2007 (after burning and before cropping); bottom right side, a picture of overland flow in a teak tree plantation during a rainfall event; centre, the figure shows plots of measurements conducted at the S4 monitoring station: annual sediment yield (SY) including suspended load and bed load contributions, box plots of overland flow contribution ($F_{OLF}$) during flood events, and areal percentages of the catchment covered by the main land uses; top left and right, conceptual diagrams of the extreme land use situations observed in the catchment; top centre: inter-annual median values of SY and $F_{OLF}$ for the two periods 2002–2009 (left) and 2010–2014 (right).</figcaption></figure>

and S8; Fig. 3). A massive soil loss (2499 Mg km⁻²) occurred in the sub-catchment dominated by degraded teak plantations (without understorey) in 2014, whereas this production remained very low (127 Mg km⁻²) for the sub-catchment dominated by the traditional shifting cultivation (Fig. 3c). The difference is particularly visible during the most intense event of the period (i.e. 16 September 2014; Fig. 3a; 187 mm d⁻¹; $I_{max}$ = 144 mm h⁻¹ based on 6-min measurements).

Since sediment yield is best explained and positively correlated to Rain_F and $F_{OLF}$ (Fig. 4b) and $F_{OLF}$ is best explained and positively correlated to Teak$_{>3}$ and AC (Fig. 4a), it is consistent that the higher the surface covered with fallow and secondary forest, the lower is the sediment yield. In contrast, the higher the surface covered with teak plantations and annual crops, the higher is the soil loss.

Whatever the period considered, subsurface flow was the major contributor to stream discharge during floods. This process can lead to the more frequent collapse of channel banks and the remobilisation of material deposited in the channel³³. Groundwater outflow can also decrease the concentration of suspended sediments in the river¹⁶, ²⁹.

Measurements of ¹³⁷Cs activities showed that soil erosion was so intense under teak plantations that the content of their soil surface in ¹³⁷Cs decreased to reach values measured in subsurface material. This trend is reflected by the decrease in ¹³⁷Cs activities with the increasing age of the teak plantations, revealing the major role of erosion in removing surface soil layers through the increased generation of rills and gullies¹². This result is further supported by the fact that teaks were preferentially planted in the most accessible zones from the village, i.e. in areas with the deepest and most fertile soils.

**Synthesis and recommendations for future land-use management in Southeastern Asia.** In Southeastern Asia, incentives to reduce shifting cultivation have induced rapid land-use changes, leading to a dramatic reduction in the cropland area and to a shortening of fallow periods³⁴, which counter-productively resulted in an increase in land degradation and soil erosion as well as in a decline of crop yields³⁵, ³⁶. The lack of sustainability of these transitional farming systems³⁷ has prompted the need for developing alternative production systems. These include perennial and monoculture tree plantations such as rubber (*Hevea brasiliensis*, *Muell. Arg.*) and teak (*Tectona grandis*, *L.*), which are rapidly expanding throughout montane environments²³, as observed in the Houay Pano catchment (Fig. 7).

Most studies that investigated the impact of afforestation in montane cultivated catchments have demonstrated a reduction of soil losses at the hillslope level³⁸. This trend is attributed to rainfall interception by the canopy, and to a reduction of the throughfall kinetic energy and the resulting splash erosion as the soil is better
protected by various strata of vegetation cover³⁹,⁴⁰. Furthermore, under these conditions, the soil surface hydro-dynamic properties (i.e. hydraulic conductivity) recovered, as a result of the intensification of biological activity⁴¹, improving infiltration and the vertical drainage of rainfall and thereby reducing soil losses due to overland flow.

In the region where the current research was conducted, teak tree afforestation induced an opposite hydro-sedimentary response (Fig. 7). Most farmers intentionally kept the soil bare under mature teak trees through the frequent burning of the understorey. This practice resulted from a mix of local beliefs and practical considerations including the supposed need for a reduced competition for soil water and nutrients, and an improved access to and a better circulation within plantations. The soil surface in these managed teak tree plantations is exposed to raindrops having a kinetic energy exceeding that recorded in non-intercepted rainfall; the reason for this high kinetic energy is that the large leaves of teak intercept and concentrate rain drops which subsequently fall from an effective height exceeding the threshold (8 m) above which velocity no longer increases with height due to air resistance⁴². Besides detachment, exposition of soil to the splash process generates the formation of crusts (Fig. 6c) and the progressive clogging of soil porosity. Because of these two processes, surface soil structure is destroyed and soil surface becomes impermeable, which limits subsurface flow and intensifies overland flow, increasing the generation of rill and gully erosion occurring on these steep slopes. Moreover, reforestation with teak may also decrease baseflow due to the increased root water uptake and evapotranspiration⁴³.

Based on the log-linear and linear relationships found in this current research (Fig. 4a,b), overland flow and total sediment yield by 2020 were estimated under the following scenario: (i) an expansion rate of teak tree plantations similar to that observed during the second sub-period (2008–2014) with teak plantations covering 75% of the catchment area by 2020; (ii) an unchanged fractional cover of annual crops and (iii) the stability of triggering-flood annual rainfall depth at a value corresponding to the average measured during the entire study period. Under these conditions, the statistical relationships predicted an increase of overland flow contribution during floods, from 28% in 2014 to approximately 40% in 2020, with a concomitant tremendous increase in annual sediment yield, from 429 Mg km⁻²year⁻¹ in 2014 up to approximately 600 Mg km⁻² year⁻¹ in 2020, which exceeds by far the soil loss rate of 250 Mg km⁻²year⁻¹ considered by most soil scientists as the highest tolerable at the plot scale⁴⁴, and which generally exceeds by far soil loss rates observed at the catchment scale.

As a recommendation, removal of ground vegetation regrowth and burning of the understorey vegetation under teak tree plantations should be strongly discouraged to avoid the degradation of the soil surface (soil surface crusting) and the subsequent increased generation of overland flow and soil loss. Alternatively, planting a shade-tolerant crop as an understorey potentially providing some economic return could be tested; candidate species include ginger (Zingiber officinale), cardamom (Elettaria cardamomum), sweet potato (Ipomoea batatas) or broom grass (Thysanolaena maxima). However, improper harvesting practices of these tuber or rhizome crops may also generate severe soil erosion. One promising avenue to improve the system is to train the farmers to forestry with lower planting densities without necessarily reducing the timber yield per surface unit.

## Methods

**Study site and land uses.** The Houay Pano catchment (60 ha) is located 10 km south of Luang Prabang in northern Laos, in the Mekong River basin (Fig. 1). Mean annual temperature is 25.3°C¹⁰. The tropical monsoon climate of the region is characterized by the succession of dry and wet seasons, with 80% of annual rainfall occurring from May to October⁴⁵. Water and sediment fluxes were monitored at the outlet of the 60-ha catchment (Fig. 1a)¹ since 2001 (at the so-called S4 station). In addition, during the rainy season in 2014, fluxes were monitored at the outlet of two zero-order headwater nested catchments (i.e. hollow valleys, so-called S7 and S8 stations). Soils consist of deep (> 2 m) and moderately deep (> 0.5 m) Alfisols, except along crests and ridges where Inceptisols are found⁴⁶. The main land uses in the area are annual crops (upland rice and Job's tear), fallow, teak plantations and secondary forests⁸,⁴⁷. Their respective dominance in the catchment changed during the last decade: the slash-and-burn cultivation system that dominated between 2002 and 2007 has gradually been replaced with teak tree plantations⁴⁸ (Fig. 2b). Detailed land-use surveys and mapping were conducted each year. The areal percentages of each land use (±10%) in the catchment were calculated using QGis (http://www.qgis.org/). S7 and S8 sub-catchments have similar soil characteristics (Fig. 1a), surface areas (0.60 and 0.57 ha) and mean slope gradients (62 and 54%) although, during the monitoring year 2014, they greatly varied in land use. In S7, the main land uses were teak monoculture with an abundant understorey preserved from burning (22% of surface area), and typical shifting cultivation mosaics (78% of the surface including upland rice 23%, Job's tears 14%, fallow 34%, secondary forest 7%). In contrast, the S8 sub-catchment was covered with teak monoculture with a degraded soil surface in lower parts (42% of the total surface), and secondary forest (58%) in upper parts.

**Hydro-sedimentary monitoring.** Rainfall was measured every 6-min using 1 automatic tipping-bucket gauge and 6 manual gauges across the catchment (Fig. 1). The maximal spatial variability of annual rainfall is estimated at 20%. Stream discharge was calculated from the continuous monitoring of the water level since 2002. Water level was recorded with an OTT Thalimedes device after every 0.5 cm variation. Uncertainties on river discharge are estimated to be below 10%. Suspended sediment was collected using an automatic sampler (Fig. 6e) during each flood event during the 2002–2014 period (samples were collected after every 2-cm water level variations during the rising stage and after every 5-cm water level variations during the falling stage). Bedload sediment was estimated each month using a sediment trap (Fig. 6e). Suspended and bed load sediments were estimated following the approach described by Valentin *et al.*¹. An average of 20 river water (600 mL) samples was collected ~10 cm below the river surface during each flood event. Overland flow was estimated based on electrical conductivity⁶ measured in these samples. Infiltration rates were estimated by calculating the difference between rainfall and overland flow divided by the duration of the rainfall event. Total sediment fluxes correspond to the
export of both bedload and suspended sediment, and these values are provided with 10% uncertainties¹. All hydro-sedimentary data are available on-line (http://msecobs-mip.fr/).

**Soil and sediment sampling.** Soil surface samples (top 0–5 cm) were collected with 100-cm³ soil density cylinders along toposequences⁸ between 2000 and 2005 (n = 55) and in 2012 (n = 32) to characterize potential sources of sediment (Fig. 1)²⁹. In addition, gully and channel bank samples (n = 23) were collected in 2012 by scraping the material over the entire exposed surface of channel banks (up to 50 cm) and gullies (1–5 cm) to characterize subsurface sources (Fig. 1)²⁷. In teak plantations, soil surface samples (top 0–5 cm, n = 35) were collected along toposequences by scraping soil, 20 cm upstream and downstream of trees, as well as between trees, under teak plantations of different ages (5, 19, 25 years). Samples were then regrouped for each type of plantation. Finally, a composite sediment sample representative of the total annual export at the catchment outlet (S4 station) was prepared for 2002 (dominance of slash-and-burn cultivation) and 2012 (period of expansion of teak plantations).

**Fallout radionuclide activities of the soil and sediment samples.** Caesium-137 ($^{137}$Cs) was measured at the Laboratoire des Sciences du Climat et de l'Environnement (Gif-sur-Yvette, France) following the methods detailed by Gourdin *et al.*²⁹. All results are expressed in Bq kg⁻¹ with all activities decay-corrected to 2012. Uncertainties on radionuclide activities reached up to 30%, depending on the counting time and the quantity of material available for analysis.

**Statistical analyses.** Multiple regressions were computed to determine linear or log-linear relationships⁴⁹,⁵⁰ between the independent variables calculated for each year, namely contribution of overland flow during floods ($F_{OLF}$), annual median percentage of total stormflow) and total annual soil loss (SY), and several candidate explanatory variables calculated for each year (i.e. Rain, Rain_F, OLF, $F_{OLF}$, InfRate, SL, BL, SY, Fw&Fo, AC, Teak, Teak_3). To predict each independent variable, the selection of the best set of explanatory variables was guided by 'best subsets' and 'step-wise' regressions, two selection algorithms –available in the Minitab 17 free trial statistical packages (https://www.minitab.com/en-us/products/minitab/) – that intend to maximize the prediction R-squared ($R^2_{pred}$) calculated by leave-one-out cross-validations. $R^2_{pred}$ reflects the ability of the model to predict observations which were not used in the model calibration. An explanatory variable was considered to be statistically significantly different from zero when its *p*-value, derived from the Student's *t*-test, was lower than 0.05. The required homoscedasticity of the model residuals was verified by visual inspection. Possible multi-collinearity among the explanatory variables was controlled with the variance inflation factor⁵¹.

The Tukey-Kramer Honestly Significant Difference (HSV) test was applied to compare $^{137}$Cs activities measured in various soil and sediment sample groups, after rejection of the ANOVA null hypothesis of equal means (Table 2).

**References**

1. Valentin, C. et al. Runoff and sediment losses from 27 upland catchments in Southeast Asia: Impact of rapid land use changes and conservation practices. *Agriculture, Ecosystems & Environment* **128**, 225–238, doi:10.1016/j.agee.2008.06.004 (2008).

2. Thothong, W. et al. Impact of land use change and rainfall on sediment and carbon accumulation in a water reservoir of North Thailand. *Agriculture, Ecosystems & Environment* **140**, 521–533, doi:10.1016/j.agee.2011.02.006 (2011).

3. Chartin, C. et al. Tracking the early dispersion of contaminated sediment along rivers draining the Fukushima radioactive pollution plume. *Anthropocene* **1**, 23–34, doi:10.1016/j.ancene.2013.07.001 (2013).

4. Costa, D., Burlando, P. & Priadi, C. The importance of integrated solutions to flooding and water quality problems in the tropical megacity of Jakarta. *Sustainable Cities and Society* **20**, 199–209 (2016).

5. Rochelle-Newall, E. J. et al. Effect of land use and hydrological processes on Escherichia coli concentrations in streams of tropical, humid headwater catchments. *Scientific reports* **6**, 32974, doi:10.1038/srep32974 (2016).

6. Ribolzi, O. et al. Use of fallout radionuclides ($^{7}$Be, $^{210}$Pb) to estimate resuspension of Escherichia coli from streambed sediments during floods in a tropical montane catchment. *Environmental science and pollution research international* **23**, 3427–3435, doi:10.1007/s11356-015-5595-z (2016).

7. Rochelle-Newall, E., Nguyen, T. M., Le, T. P., Sengtaheuanghoung, O. & Ribolzi, O. A short review of fecal indicator bacteria in tropical aquatic ecosystems: knowledge gaps and future directions. *Front Microbiol* **6**, 308, doi:10.3389/fmicb.2015.00308 (2015).

8. Huon, S. et al. Long-term soil carbon loss and accumulation in a catchment following the conversion of forest to arable land in northern Laos. *Agriculture, Ecosystems & Environment* **169**, 43–57, doi:10.1016/j.agee.2013.02.007 (2013).

9. Patin, J. et al. Analysis of runoff production at the plot scale during a long-term survey of a small agricultural catchment in Lao PDR. *Journal of Hydrology* **426–427**, 79–92, doi:10.1016/j.jhydrol.2012.01.015 (2012).

10. Ribolzi, O. et al. Impact of slope gradient on soil surface features and infiltration on steep slopes in northern Laos. *Geomorphology* **127**, 53–63, doi:10.1016/j geomorph.2010.12.004 (2011).

11. Janeau, J. L. et al. Soil erosion, dissolved organic carbon and nutrient losses underdifferent land use systems in a small catchment in northern Vietnam. *Agricultural Water Management* **146**, 314–323, doi:10.1016/j.agwat.2014.09.006 (2014).

12. Chaplot, V., Coadou le Brozec, E., Silvera, N. & Valentin, C. Spatial and temporal assessment of linear erosion in catchments under sloping lands of northern Laos. *Catena* **63**, 167–184, doi:10.1016/j.catena.2005.06.003 (2005).

13. Valentin, C., Poesen, J. & Li, Y. Gully erosion: impacts, factors and control. *Catena* **63**, 132–153 (2005).

14. Chaplot, V. & Poesen, J. Sediment, soil organic carbon and runoff delivery at various spatial scales. *Catena* **88**, 46–56, doi:10.1016/j.catena.2011.09.004 (2012).

15. Ribolzi, O. et al. Contribution of groundwater and overland flows to storm flow generation in a cultivated Mediterranean catchment. Quantification by natural chemical tracing. *Journal of Hydrology* **233**, 241–257 (2000).

16. Gourdin, E. et al. Sources and export of particle-borne organic matter during a monsoon flood in a catchment of northern Laos. *Biogeosciences* **12**, 1073–1089, doi:10.5194/bg-12-1073-2015 (2015).

17. Pearce, A., Stewart, M. & Sklash, M. Storm runoff generation in humid headwater catchments: 1. Where does the water come from? *Water Resources Research* **22**, 1263–1272 (1986).

18. Sklash, M., Stewart, M. & Pearce, A. Storm Runoff Generation in Humid Headwater Catchments: 2. A Case Study of Hillslope and Low-Order Stream Response. *Water Resources Research* **22**, 1273–1282 (1986).
19. Evrard, O., Vandaele, K., Bielders, C. & Wesemael, Bv Seasonal evolution of runoff generation on agricultural land in the Belgian loess belt and implications for muddy flood triggering. *Earth Surface Processes and Landforms* **33**, 1285–1301, doi:10.1002/esp.1613 (2008).

20. Ribolzi, O. et al. Land use and water quality along a Mekong tributary in northern Lao P.D.R. *Environmental management* **47**, 291–302, doi:10.1007/s00267-010-9593-0 (2011).

21. Sidle, R. C. et al. Erosion processes in steep terrain—Truths, myths, and uncertainties related to forest management in Southeast Asia. *Forest Ecology and Management* **224**, 199–225, doi:10.1016/j.foreco.2005.12.019 (2006).

22. Häring, V., Fischer, H. & Stahr, K. Erosion of bulk soil and soil organic carbon after land use change in northwest Vietnam. *Catena* **122**, 111–119 (2014).

23. Ziegler, A. D. et al. Environmental Consequences of the Demise in Swidden Cultivation in Montane Mainland Southeast Asia: Hydrology and Geomorphology. *Human Ecology* **37**, 361–373, doi:10.1007/s10745-009-9258-x (2009).

24. Bruijnzeel, L. Forestation and dry season flow in the tropics: A closer look. *Journal of Tropical Forest Science*, 229–243 (1989).

25. Ghimire, C., Bruijnzeel, L., Lubczynski, M. & Bonell, M. Negative trade-off between changes in vegetation water use and infiltration recovery after reforesting degraded pasture land in the Nepalese Lesser Himalaya. *Hydrology and Earth System Sciences* **18**, 4933 (2014).

26. Zégre, N., Skaugset, A. E., Som, N. A., McDonnell, J. J. & Ganio, L. M. In lieu of the paired catchment approach: Hydrologic model change detection at the catchment scale. *Water Resources Research* **46**, n/a–n/a, doi:10.1029/2009wr008601 (2010).

27. Ritchie, J. C. & McHenry, J. R. Application of radioactive fallout cesium-137 for measuring soil erosion and sediment accumulation rates and patterns: a review. *Journal of Environmental Quality* **19**, 215–233 (1990).

28. Olley, J., Burton, J., Smolders, K., Pantus, F. & Pietsch, T. The application of fallout radionuclides to determine the dominant erosion process in water supply catchments of subtropical South-east Queensland, Australia. *Hydrological Processes* **27**, 885–895, doi:10.1002/hyp.9422 (2013).

29. Gourdin, E. et al. Suspended sediment dynamics in a Southeast Asian mountainous catchment: Combining river monitoring and fallout radionuclide tracers. *Journal of Hydrology* **519**, 1811–1823, doi:10.1016/j.jhydrol.2014.09.056 (2014).

30. Evrard, O. et al. Combining multiple fallout radionuclides ($^{137}$Cs, $^{90}$Be, $^{210}$Pb$_{x}$) to investigate temporal sediment source dynamics in tropical, ephemeral riverine systems. *Journal of Soils and Sediments* **16**, 1130–1144, doi:10.1007/s11368-015-1316-y (2016).

31. Valentin, C., Rajot, J.-L. & Mitja, D. Responses of soil crusting, runoff and erosion to fallowing in the sub-humid and semi-arid regions of West Africa. *Agriculture, ecosystems & environment* **104**, 287–302 (2004).

32. Wischmeier, W. H. & Smith, D. D. Predicting rainfall erosion losses-A guide to conservation planning. *Predicting rainfall erosion losses-A guide to conservation planning* (1978).

33. Onda, Y., Tsujimura, M. & Tabuchi, H. The role of subsurface water flow paths on hillslope hydrological processes, landslides and landform development in steep mountains of Japan. *Hydrological Processes* **18**, 637–650 (2004).

34. Linquist, B., Trösch, K., Pandey, S., Phouynyavong, K. & Guenat, D. Montane paddy rice: Development and effects on food security and livelihood activities of highland Lao farmers. *Mountain Research and Development* **27**, 40–47 (2007).

35. Lestrelin, G. & Giordano, M. Upland development policy, livelihood change and land degradation: interactions from a Laotian village. *Land Degradation & Development* **18**, 55–76 (2007).

36. Lestrelin, G., Vigiak, O., Pelletreau, A., Keohavong, B. & Valentin, C. in *Natural Resources Forum*. 63–75 (Wiley Online Library).

37. Saito, K. et al. Response of traditional and improved upland rice cultivars to N and P fertilizer in northern Laos. *Field crops research* **96**, 216–223 (2006).

38. Appanah, S. et al. Forest landscape restoration for Asia-Pacific forests: a synthesis. *Forest landscape restoration for Asia-Pacific forests*, 1 (2016).

39. Goebes, P. et al. Throughfall kinetic energy in young subtropical forests: Investigation on tree species richness effects and spatial variability. *Agricultural and Forest Meteorology* **213**, 148–159 (2015).

40. Rahman, S. A., Rahman, M. F. & Sunderland, T. Increasing tree cover in degrading landscapes: ‘Integration’and ‘Intensification’of smallholder forest culture in the alutilla valley, matiranga, Bangladesh. *Small-scale Forestry* **13**, 237–249 (2014).

41. Ghimire, C. P. et al. The effects of sustained forest use on hillslope soil hydraulic conductivity in the Middle Mountains of Central Nepal. *Ecohydrology* **7**, 478–495 (2014).

42. Vis, M. Interception, drop size distributions and rainfall kinetic energy in four Colombian forest ecosystems. *Earth Surface Processes and Landforms* **11**, 591–603 (1986).

43. Bonell, M. et al. The impact of forest use and reforestation on soil hydraulic conductivity in the Western Ghats of India: implications for surface and sub-surface hydrology. *Journal of Hydrology* **391**, 47–62 (2010).

44. USDA. National Soil Survey Handbook. (United States Department of Agriculture, 1995).

45. Ribolzi, O. et al. Assessment of water quality along a tributary of the Mekong river in a mountainous, mixed land use environment of the Lao P.D.R. *Lao Journal of Agriculture and Forestry* **2008**, 91–111 (2008).

46. Chaplot, V., Podwojewski, P., Phachomphon, K. & Valentin, C. Soil Erosion Impact on Soil Organic Carbon Spatial Variability on Steep Tropical Slopes. *Soil Sci. Soc. Am. J.* **73**, 769–779, doi:10.2136/ssaj2008.0031 (2009).

47. Lacombe, G. et al. Contradictory hydrological impacts of afforestation in the humid tropics evidenced by long-term field monitoring and simulation modelling. *Hydrology and Earth System Sciences* **20**, 2691–2704, doi:10.5194/hess-20-2691-2016 (2016).

48. de Rouw, A., Soulilad, B., Phanthavong, K. & Dupin, B. The adaptation of upland rice cropping to ever-shorter fallow periods and its limit. *Poverty reduction and shifting cultivation stabilisation in the uplands of Lao PDR*. Vientiane: NAFRI, 139–148 (2005).

49. Hrachowitz, M. et al. A decade of Predictions in Ungauged Basins (PUB)—a review. *Hydrological sciences journal* **58**, 1198–1255 (2013).

50. Lacombe, G. et al. Multivariate power-law models for streamflow prediction in the Mekong Basin. *Journal of Hydrology: Regional Studies* **2**, 35–48, doi:10.1016/j.ejrh.2014.08.002 (2014).

51. Helsel, D. & Hirsch, R. Statistical methods in water resources, Chapter A3, Techniques of Water-Resources Investigations of the United State Geological Survey. *US Geological Survey, Reston, Virginia*. (Available at: http://water.usgs.gov/pubs/twri/twri4a3/) (2002).

## Acknowledgements

The authors would like to thank the Institut de Recherche pour le Développement (IRD), the Lao Department of Agriculture Land Management (DALaM), the MSEC project (Multi-Scale Environment Changes), as well as the SOERE RBV (Réseau des Bassins Versants) and the International Water Management Institute (IWMI) for their support of the long-term monitoring network MSEC. For this investigation, the authors received financial support from the French National Research Agency (TecItEasy project; ANR-13-AGRO-0007) and the CGIAR research program on Integrated Systems for the Humid Tropics. The ‘SErvice des DOnnées de l’OMP’ (SEDOO) is thanked for the database web portal management. Fieldwork and laboratory analyses conducted by Florian Bourgeois, and hydrograph separations made by Marie Arnoux and William Limon are greatly acknowledged.
**Author Contributions**

O.R. and C.V. designed research; O.R. and OE wrote the main manuscript; O.E., S.H., N.S. and K.O.L. conducted sampling and field measurements; B.S. and A.d.R. performed land-use mapping; O.E. and I.L. conducted gamma spectrometry analyses; G.L. conducted the statistical analyses on hydro-sedimentary data; A.P., O.S and C.V. participated to the redaction; all authors reviewed the manuscript.

**Additional Information**

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017