# Disentangling the relative impacts of climate change and human activities on fluvial sediment supply to the coast by the world's large rivers: Pearl River Basin, China

Roshanka Ranasinghe, Chuang Shou Wu, John Conallin, Trang Minh Duong, Edward J. Anthony

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

## To cite this version:

Roshanka Ranasinghe, Chuang Shou Wu, John Conallin, Trang Minh Duong, Edward J. Anthony. Disentangling the relative impacts of climate change and human activities on fluvial sediment supply to the coast by the world's large rivers: Pearl River Basin, China. Scientific Reports, 2019, 9 (1), pp.1-10. 10.1038/s41598-019-45442-2. hal-02476670

HAL Id: hal-02476670

https://hal.science/hal-02476670

Submitted on 23 Dec 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

# Disentangling the relative impacts of climate change and human activities on fluvial sediment supply to the coast by the world's large rivers: Pearl River Basin, China

Roshanka Ranasinghe<fnref n="1" />,²,³, Chuang Shou Wu⁴,⁵, John Conallin<fnref n="1" />,⁶, Trang Minh Duong<fnref n="1" />,² & Edward Jamal Anthony⁷,⁸

The world's large rivers are under stress and experiencing unprecedented changes in hydrology, ecosystems, and fluvial sediment loads. Many of these rivers terminate at the great deltas of the world (home to 500 million people), which depend on fluvial sediments for their very existence. While fluvial sediment loads of large rivers have already been greatly modified by human activities, climate change is expected to further exacerbate the situation. But how does the effect of climate change on fluvial sediment loads compare with that of human impacts? Here, we address this question by combining historical observations and 21st century projections for one of the world's largest 25 rivers containing two mega dams; Pearl River, China. Our analysis shows that variations in fluvial sediment supply to the coast from the Pearl river over a ~150 year study period are dominated by human activities. Projected climate change driven 21st century increases in riverflow will only compensate for about 1% of the human induced deficit in sediment load, leading to the coastal zone being starved of about 6000 Mt of sediment over the remainder of this century. A similar dominance of human impacts on fluvial sediment supply is likely at other heavily engineered rivers.

Received: 8 January 2019

Accepted: 6 June 2019

Published online: 25 June 2019

The flood plains and the deltaic coastal zones of the world's large rivers, inhabited by some 2.7 billion people, have historically been the cornerstone of human civilization<fnref n="1" />⁻⁶. However, due to rapidly increasing demands for water, food, land and power, a number of human-induced stressors are resulting in large scale and potentially irrevo-cable changes in these natural systems⁷⁻¹⁰. One such change is the sustained erosion and ensuing land losses in the heavily populated deltas and adjacent coastlines (referred to hereon simply as *coastal zones* for convenience) due to the reduction of fluvial sediment supply to coast. More often than not, upstream dam constructions is the main driver of such reductions in fluvial sediment supply to the coast¹¹,¹². For example, the Three Gorges Dam in China has been retaining about 80% of sediment from upstream since it entered into operation in 2003, resulting in a massive reduction in fluvial sediment supply to the coast by the Yangtze river¹³⁻¹⁵. By 2005, more than 45,000 reservoirs over 15 meters high (i.e. mega dams) had been constructed across the world, which has trapped around 50% of the natural sediment load¹⁶ leading to dramatic decreases in fluvial sediment supply to the coast¹¹,¹².

Potential future land losses in coastal zones due to decreases in fluvial sediment supply is an increasing global concern due to the combination of (a) the significant increase in the construction of/plans for new mega dams

<footnote n="1">¹Department of Water Science and Engineering, IHE Delft Institute for Water Education P.O. Box 3015 2610, DA, Delft, The Netherlands. ²Harbour, Coastal and Offshore Engineering, Deltares, PO Box 177, 2600 MH, Delft, The Netherlands. ³Water Engineering and Management, Faculty of Engineering Technology, University of Twente, PO Box 217, 7500 AE, Enschede, The Netherlands. ⁴Zhejiang Institute of Hydraulics & Estuary, Hangzhou, 310020, China. ⁵Key Laboratory of Estuarine and Coastal Research in Zhejiang Province, Hangzhou, 030020, China. ⁶Institute for Land Water and Society, Charles Sturt University, Albury - Wodonga, Elizabeth Mitchell Dr, Thurgoona, NSW, 2640, Australia. ⁷Aix Marseille Univ, CNRS, IRD, INRA, Coll France, CEREGE, Aix-en-Provence, France. ⁸CNRS, UG, IFREMER, LEEISA USR 3456, Centre de recherche de Montabo, Cayenne, Guyane française, France. Correspondence and requests for materials should be addressed to R.R. (email: r.ranasinghe@un-ihe.org)</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** (a) Map of the Pearl River basin including elevation, tributaries and gauging stations; (b) Location of the Pearl River basin in China.</figcaption></figure>

over the last two decades⁶,¹⁷⁻¹⁹ (b) the continued growth of coastal communities (~10% of the global population)²⁰⁻²⁵, and (c) the eroding trend along at least 24% and up to 70% of the world's sandy coastlines²⁶. Moreover, although several studies have shown that fluvial sediment loads in large rivers have been more impacted by human activities than climate change over the 20th century (and early 21st century)¹⁵,²⁷⁻²⁹, projected accelerated climate change during the remainder of the 21st century, and beyond³⁰, is widely expected to affect riverflows and sediment loads in the world's largest rivers in an unprecedented way⁸,³¹⁻³³. Nevertheless, to date there is little quantitative understanding of how the relative contributions of climate change and human activities to fluvial sediment loads may evolve over the coming decades in many of the world's important river systems. Such quantitative knowledge is however crucial to enable effective basin scale catchment and coastal zone management²⁹,³⁴,³⁵. Here, we take a step towards addressing this knowledge gap by combining historical observations and 21st century projections to quantify the relative contributions from both human impacts and climate change to the fluvial sediment supplied to the coast by the Pearl River, China over a ~150 year period (mid-20th century - end of the 21st century).

## Pearl river, china

Pearl river, China, is the 2nd largest river in China in terms of water discharge (after the Yangtze river), and being among the largest 25 rivers in the world, is a globally significant system³⁶⁻³⁹. The Pearl river basin (Fig. 1) is home to 230 million people³⁹ and its catchment area spans 450,000 km² with a catchment-wide average precipitation of 1470 mm/yr³⁸. The average annual water discharge of the Pearl river is 280 million m³/yr and the sediment yield in the upper catchment is 71.5 Mt/yr³⁹. The basin has been subjected to increased anthropogenic interventions since the 1980s due to the implementation of China's open door policy and other reform policies such as the Great Western Development Plan¹⁶. Since the late 1980's, the Pearl River delta has been subjected to massive land reclamations, leading to exponential urban and socio-economic growth⁴⁰, resulting in a current population exceeding 50 million in the delta area alone⁴¹. The Pearl River contains two mega dams in its upper reaches: the Yangtan dam (110 m high with a maximum reservoir capacity of 3,432,000,000 m³; constructed 1985–1995) and the Longtan dam (216 m high with a maximum reservoir capacity of 27,270,000,000 m³; constructed 2001–2009), the latter being the world's tallest gravity dam.

The Pearl River basin is a much studied system, with numerous previous studies focusing on diverse system characteristics such as the formation and evolution of the delta⁴², hydrological changes in response to human activities⁴³, sediment trapping rates in dams³⁸,³⁹, delta shoreline changes⁴¹,⁴⁴⁻⁴⁶, estuarine engineering (e.g. sand excavation, reclamation, and dike building)⁴⁷ and fluvial sediment loads³²,³³,³⁸,³⁹,⁴⁸. However, all previous studies on the latter, which is the focus of this study, have concentrated on changes observed over the second half of the 20th century and the first decade or so of the 21st century. Here, for the first time, we investigate fluvial sediment loads in the Pearl River over the extended period 1954–2100; a period that not only encompasses marked historical changes caused by human activities, but also the future effects of projected accelerated climate change³⁰. Specifically, here we analyse changes in fluvial sediment loads in the historical period 1954–2013, and derive future projections of fluvial sediment loads up to the end of the 21st century under a worst case climate change scenario. The analysis of results provides quantitative insights on the relative contributions of climate change and
human activities to variations in fluvial sediment supply to the coast over this ~150 year period spanning the past and the future.

**The past: 1954–2013.** Monthly riverflow and sediment load data measured at 9 stations from 1954–2013 along the Pearl River (see methods and Fig. S1 in Supplementary Material) were used to gain insights on the spatio-temporal evolution of the fluvial sediment transport along the river. The measurement stations (XLT, TE, DA, QJ, DHJK, WZ, GY, SJ and BL) are shown in Fig. 1. To compute riverflow and sediment load at the Pearl river outlet (virtual station PR), the data at stations GY, SJ, and BL were added together.

Previous studies³³,³⁸,³⁹ have indicated that, in terms of fluvial sediment loads in the Pearl river, the historical period from 1954 to 2013 can be divided into four discrete phases. It should be noted that while the exact timing of the four phases may be somewhat different among the West, North and East rivers of the Pearl river basin, Wu et al.³⁹ have shown that the combined fluvial sediment load in the North and East rivers (SJ and BL) is about 10% of the total sediment supply to the coast (GY, SJ and BL combined). Therefore, as the focus of this study is the fluvial sediment supply to the coast from the Pearl river basin as a whole, here we have considered the same phase timings at all stations to enable the basin aggregated calculations that are necessary for this analysis.

From the 1950s to the late 1970s, fluvial sediment loads in the Pearl river system have mainly been governed by river discharge, followed by an increased impact of deforestation till the end of the 1980s. From the early 1990s, dam construction has dominated the fluvial sediment loads in the Pearl River, with the period from the early 1990s to the mid-2000s being affected by the construction of the Yangtan dam, and the final period from the mid-2000s onward being affected by the construction of the Longtan dam. Following this philosophy, here we have divided the historical period 1954–2013 into the four phases shown in Fig. 2.

Analysis of the fluvial sediment loads along the main river reach (XLT to PR) during the four discrete phases shows distinct differences in the characteristics of sediment supply to the coast among the different phases. Over the period 1954–1979 (Phase 1), the sediment load increases along a downstream continuum with the upper reach (XLT to QJ) contributing over half of the total fluvial sediment supply to the coast, with further contributions from tributaries downstream. During the subsequent decade or so (1980–1991; Phase 2), fluvial sediment supply to the coast increases, mostly due to an increase in the contribution from the upper reach. From 1992–2006, the reverse behavior is observed where the fluvial sediment supply to the coast decreases below that during Phase 1. In the final period from 2007–2013 (Phase 4), fluvial sediment supply to the coast decreases even further to less than 30% that of Phase 1.

To understand the physical drivers underlying the above described temporally varying fluvial sediment transport regime of the river, here we made strategic use of sediment rating curves (i.e. empirically determined relationships between sediment load and water discharge at measurement stations) developed for the 9 stations and the ocean outlet (PR) (see methods and Supplementary Material). During Phase 1 (1954–1979), the sediment load predicted by the rating curves (when used together with measured riverflows) follows the observed sediment load data at all 10 locations (Fig. 3a). Therefore, Phase 1 can be considered as a water discharge controlled phase (or a largely “natural” phase) devoid of significant human influences. Having established this “baseline”, it is then possible to use the Phase 1 rating curves to compute what the sediment loads “would have been” during the other 3 Phases in the absence of any significant human influences. It should however be noted that the baseline period, during which no significant human impacts are expected to have occurred, might vary among different reaches of the complex Pearl river system. However, along the main West river reach (XLT – PR), the sediment load is completely governed by river discharge over the 1954–1979 period (as shown by Fig. 3a), and hence this period is suitable as the overall baseline period for the purposes of this study.

Following the above hypothesis, the sediment loads predicted when using the Phase 1 rating curves together with measured monthly riverflows during Phases 2–4 were compared with the observed sediment loads during these Phases (Fig. 3b–d and Tables S1 and S2 in Supplementary material). Essentially, any departure of the computed sediment loads from that observed is indicative of non-natural behavior. Figure 3b shows that during Phase 2 (1980–1991), observed values at most stations exceeded the predicted sediment load along the main reach of the Pearl River basin, with an ultimate observed increase of 7% in fluvial sediment supply to the coast. This could be attributed to land use changes in the upper reaches due to Chinese economic reform policies which resulted in much of the area being cleared and utilized for agricultural production⁴⁰. During Phase 3 (1992–2006), however, the observed sediment load sharply decreases compared to the predicted sediment load (Fig. 3c), with the fluvial sediment supply to the coast decreasing by about 25% compared to Phase 1. The main cause for this decrease in sediment load at the middle and lower reaches of the river is the construction of the Yangtan mega dam in the upper reach (indicated by YD in Fig. 2c). During Phase 4 (2007–2013) the sediment load decreases further compared to predicted sediment load (Fig. 3d), resulting in a massive 71% decrease in fluvial sediment supply to the coast compared to Phase 1. This is attributed to the placement of Longtan mega dam (indicated by LD in Fig. 2d) upstream of the Yangtan dam, trapping approximately 99% of the sediment input from the upper tributaries of the river³⁹.

**Climate change versus Human impacts over 150 years (1954–2100).** For the purposes of this study, it is necessary to determine the relative contributions of climate change and human activities on fluvial sediment supply to the coast in a way that is consistent for both the past and the future (given that measured sediment loads are not available for the future). To that end, sediment loads predicted by the rating curves described above were here used together with reasonable assumptions. The computational procedure followed is illustrated in *Methods*.

First, the total fluvial sediment supply to the coast (i.e. at location PR) was computed under the combined effect of climate change and human activities. For Phases 1 to 4, the mean fluvial sediment supply at PR per Phase was computed using the monthly measured riverflow data and rating curves for PR (see Table S3 in Supplementary Material). Post-2014 riverflows were estimated using the 0–10% increase over the 21st century
<figure><img src="image_5.png" /></figure>

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

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

<figure><img src="image_8.png" /><figcaption>**Figure 2.** The spatio-temporal evolution of fluvial sediment loads along the Pearl river from 1954–2013. For station names, please see Fig. 1. YD and LD indicate Yangtan and Longtan dams respectively.</figcaption></figure>

projected for the Pearl River by Yan *et al.*⁴⁹. Here, assuming a high end scenario, an increase of 10% in the runoff (by 2079–2099, relative to 1979–1999) was used in all calculations. The fluvial sediment supply for the period 2014–2099 was then calculated with the rating curve obtained for Phase 4 at location PR, with the assumption that no more significant human interventions will take place at this already heavily engineered system. The fluvial sediment supply to the coast under the combined effect of climate change and human impacts over the period 1954–2100 thus calculated is shown in Fig. 4 (red line).

Subsequently, the fluvial sediment supply the coast may have received due only to natural forcing (including climate change but excluding all human impacts) was computed. When all human impacts are neglected, ideally, the system can be expected to behave as it did during the “natural” Phase 1. Furthermore, several studies have indicated that climate change driven variations in runoff during the latter half of the 20th century in the Pearl River basin are all but negligible⁵⁰. Therefore, in this “natural forcing only” computation, the PR rating curve for
<figure><img src="image_9.png" /><figcaption>**Figure 3.** Comparison between sediment loads predicted using the Phase 1 rating curves (dashed line) and observed sediment loads along the Pearl River during the 4 different Phases.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 4.** The fluvial sediment supply to the coast over the period 1954–2100 under the combined effect of climate change and human impacts (red line) and under natural forcing only, including climate change (blue circles).</figcaption></figure>

Phase 1 was used together with the measured riverflows to compute the fluvial sediment supply to the coast for the period 1954–1999. For the entirety of the 21st century, following the projections of Yan *et al.*⁴⁹, an increase of 10% in the runoff (by 2079–2099, relative to 1979–1999) was adopted. In this 21st century sediment load calculation too, the rating curve computed for the natural Phase (i.e. Phase 1) was used because all human impacts are
<table><caption>**Table 1.** Fitted rating curve parameters (a, b) and goodness of fit measures (r, p) at the 9 Pearl river measurement stations and the virtual station at the Pearl river inlet. Phase 1: 1954–1979; Phase 2:1980–1991; Phase 3: 1992–2006; Phase 4: 2007–2013. Please see Fig. 1 for station locations.</caption><thead><tr><th rowspan="2">Station</th><th colspan="5">Phase 1</th><th colspan="5">Phase 2</th><th colspan="5">Phase 3</th><th colspan="5">Phase 4</th></tr><tr><th>a</th><th>b</th><th>r</th><th>p</th><th>a</th><th>b</th><th>r</th><th>p</th><th>a</th><th>b</th><th>r</th><th>p</th><th>a</th><th>b</th><th>r</th><th>p</th><th>a</th><th>b</th><th>r</th><th>p</th></tr></thead><tbody><tr><td>XLT</td><td>0.63</td><td>1.43</td><td>0.95</td><td>&lt;0.001</td><td>0.90</td><td>1.36</td><td>0.88</td><td>&lt;0.001</td><td>0.84</td><td>1.24</td><td>0.89</td><td>=0.002</td><td>0.31</td><td>2.09</td><td>0.97</td><td>&lt;0.001</td><td>0.31</td><td>2.09</td><td>0.97</td><td>&lt;0.001</td></tr><tr><td>TE</td><td>0.03</td><td>1.66</td><td>0.95</td><td>&lt;0.001</td><td>0.12</td><td>1.58</td><td>0.68</td><td>&lt;0.001</td><td>1.04</td><td>0.86</td><td>0.69</td><td>=0.05</td><td>6*10<sup>-9</sup></td><td>1.56</td><td>0.98</td><td>&gt;0.05</td><td>6*10<sup>-9</sup></td><td>1.56</td><td>0.98</td><td>&gt;0.05</td></tr><tr><td>DA</td><td>0.03</td><td>1.78</td><td>0.94</td><td>&lt;0.001</td><td>0.005</td><td>2.32</td><td>0.76</td><td>&lt;0.001</td><td>1*10<sup>-4</sup></td><td>2.02</td><td>0.70</td><td>=0.05</td><td>2*10<sup>-8</sup></td><td>4.38</td><td>0.85</td><td>&gt;0.05</td><td>2*10<sup>-8</sup></td><td>4.38</td><td>0.85</td><td>&gt;0.05</td></tr><tr><td>QJ</td><td>0.02</td><td>1.88</td><td>0.82</td><td>&lt;0.001</td><td>0.001</td><td>2.59</td><td>0.87</td><td>&lt;0.001</td><td>8*10<sup>-6</sup></td><td>3.41</td><td>0.79</td><td>&lt;0.001</td><td>1*10<sup>-8</sup></td><td>4.60</td><td>0.74</td><td>=0.001</td><td>1*10<sup>-8</sup></td><td>4.60</td><td>0.74</td><td>=0.001</td></tr><tr><td>DHJK</td><td>0.007</td><td>1.74</td><td>0.89</td><td>&lt;0.001</td><td>0.0002</td><td>2.48</td><td>0.84</td><td>&lt;0.001</td><td>4*10<sup>-5</sup></td><td>2.64</td><td>0.80</td><td>&lt;0.001</td><td>2*10<sup>-5</sup></td><td>2.66</td><td>0.92</td><td>=0.005</td><td>2*10<sup>-5</sup></td><td>2.66</td><td>0.92</td><td>=0.005</td></tr><tr><td>WZ</td><td>0.005</td><td>1.78</td><td>0.91</td><td>&lt;0.001</td><td>0.0006</td><td>2.24</td><td>0.81</td><td>=0.002</td><td>2*10<sup>-4</sup></td><td>2.23</td><td>0.75</td><td>=0.002</td><td>2*10<sup>-5</sup></td><td>2.58</td><td>0.75</td><td>=0.05</td><td>2*10<sup>-5</sup></td><td>2.58</td><td>0.75</td><td>=0.05</td></tr><tr><td>GY</td><td>0.008</td><td>1.68</td><td>0.92</td><td>&lt;0.001</td><td>0.005</td><td>1.80</td><td>0.76</td><td>=0.005</td><td>4*10<sup>-5</sup></td><td>2.60</td><td>0.89</td><td>&lt;0.001</td><td>4*10<sup>-4</sup></td><td>2.01</td><td>0.95</td><td>=0.001</td><td>4*10<sup>-4</sup></td><td>2.01</td><td>0.95</td><td>=0.001</td></tr><tr><td>SJ</td><td>0.01</td><td>1.61</td><td>0.95</td><td>&lt;0.001</td><td>0.006</td><td>1.86</td><td>0.85</td><td>&lt;0.001</td><td>0.0002</td><td>2.59</td><td>0.96</td><td>=0.001</td><td>0.001</td><td>2.27</td><td>0.88</td><td>&lt;0.01</td><td>0.0001</td><td>2.27</td><td>0.88</td><td>&lt;0.01</td></tr><tr><td>BL</td><td>0.05</td><td>1.33</td><td>0.98</td><td>&lt;0.001</td><td>0.007</td><td>1.89</td><td>0.87</td><td>&lt;0.001</td><td>0.01</td><td>1.57</td><td>0.91</td><td>&lt;0.001</td><td>0.0002</td><td>2.69</td><td>0.99</td><td>&lt;0.01</td><td>0.0002</td><td>2.69</td><td>0.99</td><td>&lt;0.01</td></tr><tr><td>PR</td><td>0.009</td><td>1.59</td><td>0.86</td><td>&lt;0.001</td><td>0.007</td><td>1.69</td><td>0.75</td><td>=0.005</td><td>2*10<sup>-5</sup></td><td>2.58</td><td>0.91</td><td>&lt;0.001</td><td>0.0009</td><td>1.84</td><td>0.96</td><td>&lt;0.001</td><td>0.0009</td><td>1.84</td><td>0.96</td><td>&lt;0.001</td></tr></tbody></table>

being neglected in this calculation. The potential fluvial sediment supply to the coast under only natural forcing (including climate change) for the entire 1954–2100 period is shown in Fig. 4 (blue circles).

Figure 4 clearly shows the dramatic difference between what could have been (under natural forcing only) and the impact of human activities to date (assuming that no significant human interventions, positive or negative, will take place during the remainder of the 21st century). The percentage reduction of fluvial sediment supply to the coast due to human impacts from 1954 to 2013 is 71%, as shown in Fig. 4. The projected CC driven increase in riverflow during the 21st century will only compensate for 1% of this deficit by 2100.

Although a recent analysis of multi-decadal shoreline mobility in response to changes in fluvial sediment supply shows that the subaerial Pearl delta is still advancing seaward⁵¹ (largely due to coastal land reclamations⁴¹), other recent studies⁵²⁻⁵⁴ that have also examined the evolution of the sub-aqueous bathymetry of the Pearl delta have shown that parts of the delta are in fact presently eroding or accreting at a decreasing rate. Therefore, like several other large deltas in the world that have not yet been significantly or visibly impacted, the decrease in fluvial sediment supply to the coast identified here is likely to change the morphological evolution of the coastal zone adjacent to the Pearl River (inclusive of the Pearl delta) in the years ahead. For example, erosion of Yangtze delta is reported to have commenced when the fluvial sediment supply to the coast decreased by 50%⁵⁵. Given that the computed decrease in fluvial sediment supply from 1954–2013 by the Pearl River is already 71%, it is therefore highly likely that the Pearl delta will also face a similar fate, sooner rather than later. Furthermore, sediment starvation by the river-adjacent coastline could lead to sustained coastline recession as has also been observed at the Yangze⁵⁵,⁵⁶, Ebro⁵⁷ and Mekong rivers²⁶. The projected decrease in fluvial sediment loads in the Pearl River could also have significant environmental implications, due to the associated reduction in the ability of the water column to self-decontaminate³⁸,⁴⁸. A potential positive impact of the projected decrease in fluvial sediment supply is increased nearshore zone primary productivity. Enhanced photosynthesis of planktons due to the potential increase in transparency of seawater associated with a reduction in suspended sediment concentration in the water column could result in such an increase in primary productivity in the nearshore zone³⁸,⁴⁸. This type of environmental impacts are not only relevant for the Pearl Delta, arguably one of the key economic regions in China, but also for nearby Macau and Hong Kong, and thus requires attention in future studies.

**Concluding remarks.** The stability of river-adjacent coastal zones (inclusive of deltas) is strongly dependent on the continued supply of fluvial sediment to the coast. Any significant departure from the natural fluvial sediment supply is very likely to result in sustained erosion leading to large scale coastal zone land losses. With the strategic use of historical data and future projections, this study derived quantitative insights into the relative impacts of human activities and climate change on variations in fluvial sediment supply to the coast at the globally relevant Pearl River, China, over a ~150 year period (1954–2100).

This analysis showed that the fluvial sediment supply to the coast from the Pearl river has decreased (from its natural state) by 71% from 1954–2013 due to human activities, and that projected increases in riverflow due to climate change will compensate for just 1% of this deficit by the end of the 21st century (if the impact of human activities is maintained at current levels). Based on these projections, the coastal zone of the Pearl River will be starved of approximately 6,000 Mt of fluvial sediment during the remainder of the 21st century.

As in the case of the Pearl River, the dominance of human impacts over that of climate change on fluvial sediment supply to the coast is likely at other heavily engineered rivers around the world. If remedial measures are not implemented sooner rather than later, the result will be substantial sediment starvation along river-adjacent coastal zones leading to pervasive coastal erosion and ensuing land losses. This will not only threaten existing coastal zone communities and developments but will also compromise foreshadowed coastal zone developments and expansions.

Furthermore, the results of this study underline the importance of taking into account catchment effects, particularly those due to human activities, in coastal zone management and planning efforts which have historically looked at the coast mostly from the sea side (i.e. basing all decisions on characteristics of sediment transport driven by oceanic processes such as waves, tides, extreme sea levels, sea level rise etc.). The non-inclusion of a
<figure><img src="image_11.png" /><figcaption>**Figure 5.** Computational procedure followed to calculate fluvial sediment supply to the coast over 1954–2100 under the effect of (a) both climate change and human activities, and (b) natural forcing (incl. climate change) only.</figcaption></figure>

large human induced deficit in fluvial sediment supply to the coast, such as that projected here for the Pearl River, when considering the coastal zone sediment budget could lead to sub-optimal designs and decisions, which may, under some circumstances, place coastal settlements and communities at risk.

## Methods

**Data.** Long-term water discharge and sediment load collected by the Pearl River Water Resource Commission, published data from the Bulletin of Chinese River Sediment (2003–2013)⁵⁸ and the Bulletin of Guangxi Zhuang Autonomous Region River Sediment (2004, 2006, 2007)⁵⁹ were used in this study. The annual averaged data are graphically shown in Figure S1 in Supplementary material.

**Development of sediment rating curves.** Fitting a power function (Eq. 1) between fluvial sediment load ($Q_d$) and water discharge ($Q_s$) is a commonly used method to develop sediment rating curves for rivers⁶⁰⁻⁶⁴. From the empirical relationship between $Q_s$ (km³ yr⁻¹) and $Q_d$ (Mt/yr) the following correlation can be established:

$$Q_d = a \cdot Q_s^b \qquad (1)$$

where, a is the sediment rating coefficient and b is the rating exponent.
Extensive studies of this type of functions have found that the *a* coefficient is an index of erosion severity, or the source of sediment supply. A high *a* parameter is associated with intensely weathered material which can be easily eroded and transported⁶⁵. The exponent *b* represents the erosive power of the river and the influence on the sediment supply from the entire catchment surface. The value of *b* can also be affected by the grain size distribution of the material available for transport in the river basin. Large values of *b* indicate an increase in sediment transport capacity of the river⁶². A detailed discussion on the spatio-temporal differences in the *a* and *b* coefficients as relevant to the Pearl river is provided by Zhang *et al.*⁴⁸. The rating curve approach was exclusively used in this study to develop sediment rating curves for the 9 measurement stations and the virtual station at the ocean outlet of the river (PR). The *a* and *b* coefficient values thus obtained and the associated goodness of fit statistics are given in Table 1.

**Computational procedure.** The procedures followed in this study to compute the fluvial sediment supply to the coast over the period 1954–2100 under the combined effect of climate change and human activities, and under natural forcing alone (including climate change effects) are shown in Fig. 5a,b respectively.

**Data Availability**

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

**References**

1. Hassan, F. A. The dynamics of a riverine civilization: A geoarchaeological perspective on the Nile Valley, Egypt. *World Archaeol.* **29**, 51–74 (1997).

2. Macklin, M. G. & Lewin, J. The rivers of civilization. *Quat. Sci. Rev.* **114**, 228–244 (2015).

3. Bianchi, T. S. *Deltas and Humans*. (Oxford Univ. Press, 2016).

4. Singh, A. *et al.* Counter-intuitive influence of Himalayan river morphodynamics on Indus Civilisation urban settlements. *Nat. Commun.* **8**, 1617 (2017).

5. Doyle, M. *The Source: How Rivers Made America and America Remade Its Rivers*. (W. W. Norton & Company, 2018).

6. Best, J. Anthropogenic stresses on the world's big rivers. *Nat. Geosci.* **12**, 7–21 (2019).

7. Jackson, R. B. *et al.* Water in a Changing World. *Ecol. Appl.* **11**, 1027–1045 (2001).

8. Palmer, M. A. *et al.* Climate change and the world's river basins: anticipating management options. *Front. Ecol. Environ.* **6**, 81–89 (2008).

9. Fuller, M. R., Doyle, M. W. & Strayer, D. L. Causes and consequences of habitat fragmentation in river networks. *Ann. N. Y. Acad. Sci.* **1355**, 31–51 (2015).

10. Veilleux, J. C. & Anderson, E. P. 2015 Snapshot of Water Security in the Nile, Mekong, and Amazon River Basins. *Limnol. Oceanogr. Bull.* **25**, 8–14 (2016).

11. Vörösmarty, C. J. *et al.* Anthropogenic sediment retention: major global impact from registered river impoundments. *Glob. Planet. Change* **39**, 169–190 (2003).

12. Syvitski, J. P. M. Impact of Humans on the Flux of Terrestrial Sediment to the Global Coastal. *Ocean. Science* (80-.). **308**, 376–380 (2005).

13. Yang, H. F. *et al.* Human impacts on sediment in the Yangtze River: A review and new perspectives. *Glob. Planet. Change* **162**, 8–17 (2018).

14. Yang, S. L. *et al.* Downstream sedimentary and geomorphic impacts of the Three Gorges Dam on the Yangtze River. *Earth-Science Rev.* **138**, 469–486 (2014).

15. Yang, S. L., Xu, K. H., Milliman, J. D., Yang, H. F. & Wu, C. S. Decline of Yangtze River water and sediment discharge: Impact from natural and anthropogenic changes. *Sci. Rep.* **5**, 1–14 (2015).

16. Nilsson, C., Reidy, C. A., Dynesius, M. & Revenga, C. Fragmentation and Flow Regulation of the World's Large River Systems. *Science* (80-.). **308**, 405–408 (2005).

17. Gross, M. A global megadam mania. *Curr. Biol.* **26**, R779–R782 (2016).

18. Magilligan, F. J., Snedden, C. S. & Fox, C. A. The era of big dam building: it ain't over till it's over. in *The Politics of Fresh Water Access, conflict and identity* (eds Ashcraft, C. M. & Mayer, T.) 78–97 (2017).

19. Zarfl, C., Lumsdon, A. E., Berlekamp, J., Tydecks, L. & Tockner, K. A global boom in hydropower dam construction. *Aquat. Sci.* **77**, 161–170 (2015).

20. Small, C. & Nicholls, R. J. A Global Analysis of Human Settlement in Coastal Zones. *J. Coast. Res.* **19**(3), 584–599 (2003).

21. Valiela, I. *Global Coastal Change*. (Blackwell Pub., 2006).

22. McGranahan, G., Balk, D. & Anderson, B. The rising tide: assessing the risks of climate change and human settlements in low elevation coastal zones. *Environ. Urban.* **19**, 17–37 (2007).

23. Hallegatte, S., Green, C., Nicholls, R. J. & Corfee-Morlot, J. Future flood losses in major coastal cities. *Nat. Clim. Chang.* **3**, 802–806 (2013).

24. Hinkel, J. *et al.* A global analysis of erosion of sandy beaches and sea-level rise: An application of DIVA. *Glob. Planet. Change* **111**, 150–158 (2013).

25. Ranasinghe, R. Assessing climate change impacts on open sandy coasts: A review. *Earth-Science Rev.* **160**, 320–332 (2016).

26. Luijendijk, A. *et al.* The State of the World's Beaches. *Sci. Rep.* **8**, 6641 (2018).

27. Dai, S. B., Lu, X. X., Yang, S. L. & Cai, A. M. A preliminary estimate of human and natural contributions to the decline in sediment flux from the Yangtze River to the East China Sea. *Quat. Int.* **186**, 43–54 (2008).

28. Miao, C., Ni, J., Borthwick, A. G. L. & Yang, L. A preliminary estimate of human and natural contributions to the changes in water discharge and sediment load in the Yellow River. *Glob. Planet. Change* **76**, 196–205 (2011).

29. Anthony, E. J. *et al.* Linking rapid erosion of the Mekong River delta to human activities. *Sci. Rep.* **5**, 14745 (2015).

30. IPCC. Summary for Policymakers. in *Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change* (eds Masson-Delmotte, V. *et al.*) 3–29 (Cambridge University Press, 2013).

31. Knight, J. & Harrison, S. The impacts of climate change on terrestrial Earth surface systems. *Nat. Clim. Chang.* **3**, 24–29 (2013).

32. Liu, F. *et al.* Impacts of ENSO on multi-scale variations in sediment discharge from the Pearl River to the South China Sea. *Geomorphology* **293**, 24–36 (2017).

33. Tan, C. *et al.* Using the wavelet transform to detect temporal variations in hydrological processes in the Pearl River, China. *Quat. Int.* **440**, 52–63 (2017).

34. Auerbach, L. W. *et al.* Flood risk of natural and embanked landscapes on the Ganges-Brahmaputra tidal delta plain. *Nat. Clim. Chang.* **5**, 153–157 (2015).
35. Wohl, E. et al. The Natural Sediment Regime in Rivers: Broadening the Foundation for Ecosystem Management. *Bioscience* **65**, 358–371 (2015).

36. Eisma, D. *Intertidal deposits: river mouths, tidal flats and coastal lagoons*. (CRC Press, Cop., 1998).

37. Zhao, C. H., Zhu, Z. H. & Zhou, D. Z. *Worldwide Rivers and Dams*. (China Water Conservancy and Hydroelectric Press, 2000).

38. Dai, S. B., Yang, S. L. & Cai, A. M. Impacts of dams on the sediment flux of the Pearl River, southern China. *CATENA* **76**, 36–43 (2008).

39. Wu, C. S., Yang, S. L. & Lei, Y. Quantifying the anthropogenic and climatic impacts on water discharge and sediment load in the Pearl River (Zhujiang), China (1954–2009). *J. Hydrol.* **452–453**, 190–204 (2012).

40. Weng, Q. A historical perspective of river basin management in the Pearl River Delta of China. *J. Environ. Manage.* **85**, 1048–1062 (2007).

41. Wu, Z., Milliman, J. D., Zhao, D., Zhou, J. & Yao, C. Recent geomorphic change in LingDing Bay, China, in response to economic and urban growth on the Pearl River Delta, Southern China. *Glob. Planet. Change* **123**, 1–12 (2014).

42. Jia, L., Pan, S. & Wu, C. Effects of the anthropogenic activities on the morphological evolution of the Modaomen Estuary, Pearl River Delta, China. *China Ocean Eng.* **27**, 795–808 (2013).

43. Liu, F. et al. Hydrological responses to the combined influence of diverse human activities in the Pearl River delta, China. *CATENA* **113**, 41–55 (2014).

44. Wu, C. S., Yang, S., Huang, S. & Mu, J. Delta changes in the Pearl River estuary and its response to human activities (1954–2008). *Quat. Int.* **392**, 147–154 (2016).

45. Wang, X.-Z., Zhang, H.-G., Fu, B. & Shi, A.-Q. Analysis on the coastline change and erosion-accretion evolution of the Pearl River Estuary, China, based on remote-sensing images and nautical charts. *J. Appl. Remote Sens.* **7**, 073519 (2013).

46. Zhang, W. et al. Morphological change in the Pearl River Delta, China. *Mar. Geol.* **363**, 202–219 (2015).

47. Zhang, S. et al. Recent changes of water discharge and sediment load in the Zhujiang (Pearl River) Basin. *China. Glob. Planet. Change* **60**, 365–380 (2008).

48. Zhang, W., Wei, X., Jinhai, Z., Yuliang, Z. & Zhang, Y. Estimating suspended sediment loads in the Pearl River Delta region using sediment rating curves. *Cont. Shelf Res.* **38**, 35–46 (2012).

49. Yan, D., Werners, S. E., Ludwig, F. & Huang, H. Q. Hydrological response to climate change: The Pearl River, China under different RCP scenarios. *J. Hydrol. Reg. Stud.* **4**, 228–245 (2015).

50. Tian, Q., Prange, M. & Merkel, U. Precipitation and temperature changes in the major Chinese river basins during 1957–2013 and links to sea surface temperature. *J. Hydrol.* **536**, 208–221 (2016).

51. Besset, M., Anthony, E. J. & Bouchette, F. Multi-decadal variations in delta shorelines and their relationship to river sediment supply: An assessment and review. *Earth-Science Rev.* in press (2019).

52. Liu, F. et al. Stepwise adjustment of deltaic channels in response to human interventions and its hydrological implications for sustainable water managements in the Pearl River Delta, China. *J. Hydrol.* **573**, 194–206 (2019).

53. Yang, L. et al. Morphological response of Lingding Bay in the Pearl River Estuary to human intervention in recent decades. *Ocean Coast. Manag.* **176**, 1–10 (2019).

54. Wu, Z. et al. Geomorphologic changes in the lower Pearl River Delta, 1850–2015, largely due to human activity. *Geomorphology* **314**, 42–54 (2018).

55. Yang, S. L. Impact of dams on Yangtze River sediment supply to the sea and delta intertidal wetland response. *J. Geophys. Res.* **110**, F03006 (2005).

56. Yang, S. L. et al. Drastic decrease in sediment supply from the Yangtze River and its challenge to coastal wetland management. *Geophys. Res. Lett.* **33**, L06408 (2006).

57. Mikhailova, M. V. Transformation of the Ebro River Delta under the Impact of Intense Human-Induced Reduction of Sediment Runoff. *Water Resour.* **30**, 370–378 (2003).

58. Ministry of Water Resources of the People's Republic of China (2003–2013). *Bulletin of China River Sediment* (2003–2013).

59. Guangxi Zhuang Autonomous Region (2004–2007). *Bulletin of Guangxi Zhuang Autonomous Region River Sediment* (2004–2007).

60. Walling, D. E., Russell, M. A., Hodgkinson, R. A. & Zhang, Y. Establishing sediment budgets for two small lowland agricultural catchments in the UK. *CATENA* **47**, 323–353 (2002).

61. Crawford, C. G. Estimation of suspended-sediment rating curves and mean suspended-sediment loads. *J. Hydrol.* **129**, 331–348 (1991).

62. Asselman, N. E. M. Fitting and interpretation of sediment rating curves. *J. Hydrol.* **234**, 228–248 (2000).

63. Xu, K. et al. Simulated sediment flux during 1998 big-flood of the Yangtze (Changjiang) River, China. *J. Hydrol.* **313**, 221–233 (2005).

64. Warrick, J. A. Trend analyses with river sediment rating curves. *Hydrol. Process.* **29**, 936–949 (2015).

65. Yang, S. L., Zhang, J., Dai, S. B., Li, M. & Xu, X. J. Effect of deposition and erosion within the main river channel and large lakes on sediment delivery to the estuary of the Yangtze River. *J. Geophys. Res.* **112**, F02005 (2007).

## Acknowledgements

R.R. is supported by the AXA Research fund and the Deltares Strategic Research Programme ‘Coastal and Offshore Engineering’. EA was supported by the ANR-Belmont Forum Project ‘BF-Deltas: Catalyzing Action Towards Sustainability of Deltaic Systems with an Integrated Modeling Framework for Risk Assessment’.

## Author Contributions

R.R. and C.S.W. conceived the study and jointly performed the analysis with input from J.C., T.M.D. and E.J.A. R.R. and T.M.D. performed the computations related to future projections. R.R. and J.C. wrote the manuscript, with assistance from E.J.A., T.M.D. and C.S.W.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-45442-2.

**Competing Interests:** The authors declare 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_12.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 Cre-
ative 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 per-
mitted 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) 2019