Journal of Geographical Sciences >
Quantification of human and climate contributions to multi-dimensional hydrological alterations: A case study in the Upper Minjiang River, China
Zhang Yuhang (1994-), specialized in hydrometeorological ensemble forecast. E-mail: zhangyh19@mail.bnu.edu.cn |
Received date: 2020-08-26
Accepted date: 2021-03-09
Online published: 2021-10-25
Supported by
Natural Science Foundation of China(51879009)
Natural Science Foundation of China(52079143)
Second Tibetan Plateau Scientific Expedition and Research Program(2019QZKK0405)
National Key Research and Development Program of China(2018YFE0196000)
National Key Research and Development Program of China(2017YFC0404405)
Interdisciplinary Research Foundation of Beijing Normal University for the First-Year Doctoral Students(BNUXKJC1905)
Independent Research Projects of POWERCHINA Chengdu Engineering Corporation Limited(P34516)
Dual factors of climate and human on the hydrological process are reflected not only in changes in the spatiotemporal distribution of water resource amounts but also in the various characteristics of river flow regimes. Isolating and quantifying their contributions to these hydrological alterations helps us to comprehensively understand the response mechanism and patterns of hydrological process to the two kinds of factors. Here we develop a general framework using hydrological model and 33 indicators to describe hydrological process and quantify the impact from climate and human. And we select the Upper Minjiang River (UMR) as a case to explore its feasibility. The results indicate that our approach successfully recognizes the characteristics of river flow regimes in different scenarios and quantitatively separates the climate and human contributions to multi-dimensional hydrological alterations. Among these indicators, 26 of 33 indicators decrease over the past half-century (1961-2012) in the UMR, with change rates ranging from 1.3% to 33.2%, and the human impacts are the dominant factor affecting hydrological processes, with an average relative contribution rate of 58.6%. Climate change causes an increase in most indicators, with an average relative contribution rate of 41.4%. Specifically, changes in precipitation and reservoir operation may play a considerable role in inducing these alterations. The findings in this study help us better understand the response mechanism of hydrological process under changing environment and is conducive to climate change adaptation, water resource planning and ecological construction.
ZHANG Yuhang , YE Aizhong , YOU Jinjun , JING Xiangyang . Quantification of human and climate contributions to multi-dimensional hydrological alterations: A case study in the Upper Minjiang River, China[J]. Journal of Geographical Sciences, 2021 , 31(8) : 1102 -1122 . DOI: 10.1007/s11442-021-1887-z
Figure 1 Location and attributes of the Upper Minjiang River |
Figure 2 Framework of this study |
Table 1 Formulas and description of the selected assessment criteria |
Formulas | Description | Perfect/no skill |
---|---|---|
$NSE=1-\frac{\mathop{\sum }_{i=1}^{n}{{\left( x_{sim}^{i}-x_{obs}^{i} \right)}^{2}}}{\mathop{\sum }_{i=1}^{n}{{\left( x_{obs}^{i}-\overline{{{x}_{obs}}} \right)}^{2}}}$ | Predictive skill of hydrological models and accuracy between simulations and observations | 1/≤0 |
$PCC=\frac{\mathop{\sum }_{i=1}^{n}\left[ \left( x_{sim}^{i}-\overline{{{x}_{sim}}} \right)\left( x_{obs}^{i}-\overline{{{x}_{obs}}} \right) \right]}{\sqrt{\mathop{\sum }_{i=1}^{n}{{\left( x_{sim}^{i}-\overline{{{x}_{sim}}} \right)}^{2}}}\sqrt{\mathop{\sum }_{i=1}^{n}{{\left( x_{obs}^{i}-\overline{{{x}_{obs}}} \right)}^{2}}}}$ | Linear correlation between simulations and observations | 1/≤0 |
$PBIAS=\frac{\mathop{\sum }_{i=1}^{n}\left( x_{obs}^{i}-x_{sim}^{i} \right)}{\mathop{\sum }_{i=1}^{n}x_{obs}^{i}}\times 100\text{ }\!\!%\!\!\text{ }$ | The percent difference between simulations and observations; model's performance with regard to its ability to maintain the water balance | 0/∞ |
$RMSE=\sqrt{\frac{1}{n}\underset{i=1}{\overset{n}{\mathop \sum }}\,{{\left( x_{obs}^{i}-x_{sim}^{i} \right)}^{2}}}$ | Association of simulations and observations | 0/∞ |
Table 2 Definition of IHA parameters and abbreviations used in the text (modified from Richter et al., 1996) |
Groups and ID | IHA parameters | Abbreviations |
---|---|---|
G11 | Mean value of annual flow | Annual |
G12-G113 | Mean value of 12 months | Jan-Dec |
G21 | Annual 1-day minima | QN1D |
G22 | Annual 3-day minima | QN3D |
G23 | Annual 7-day minima | QN7D |
G24 | Annual 30-day minima | QN30D |
G25 | Annual 90-day minima | QN90D |
G26 | Annual 1-day maxima | QM1D |
G27 | Annual 3-day maxima | QM3D |
G28 | Annual 7-day maxima | QM7D |
G29 | Annual 30-day maxima | QM30D |
G210 | Annual 90-day maxima | QM90D |
G211 | Base flow index | Base flow |
G31 | Julian date of each annual 1-day maxima | Date max |
G32 | Julian date of each annual 1-day minima | Date min |
G41 | Number of low pulse | Lo pulse # |
G42 | Number of high pulse | Hi pulse # |
G43 | Mean duration of low pulse | Lo pulse L |
G44 | Mean duration of high pulse | Hi pulse L |
G51 | Rising rate | Rise rate |
G52 | Falling rate | Fall rate |
G53 | Number of hydrological reversals | Reversals |
Figure 3 Annual streamflow (a) and Mann-Kendall mutation diagnosis (b) at the ZPP station during 1938-2012. UF and UB in (b) are statistics calculated by sequential and inverse streamflow records, respectively |
Table 3 Model performance for daily discharge simulations at the ZPP station. Detailed description of these criteria can be found in Table 1 |
NSE | PCC | PBIAS | RMSE | |
---|---|---|---|---|
Calibration (1961‒1965) | 0.73 | 0.86 | ‒1.51% | 177.35 |
Verification (1966‒1969) | 0.72 | 0.85 | 5.94% | 183.46 |
Figure 4 Observation and simulation of streamflow at the ZPP station during the calibration period (1961-1965) and verification period (1966-1969). (NSE: Nash-Sutcliffe efficiency; PCC: Pearson correlation coefficient) |
Figure 5 Relative changes (a) and relative contributions (b) for all IHA parameters induced by climate change versus those induced by human activities. The vertical dashed line indicates five groups. |
Figure 6 Joint distribution of the change rate for each group of IHA parameters induced by climate change $(\Delta {{I}_{c}})$ versus those induced by human activities $(\Delta {{I}_{h}}).$ The horizontal axis in the figure represents the change rate of IHA parameters caused by climate, and the vertical axis represents the change rate of IHA parameters caused by human. |
Table 4 IHA results, relative change and contributions calculated at the ZPP station. Bolded numbers indicate significant differences in the IHA indicators between two periods according to M-W U test at a significance level of 0.05. Mean (+) and mean (-) are the average values of increased and decreased parameters, respectively. |
Index | obsbp | obsap | simap | obsnorm,bp | obsnorm,ap | simnorm,ap | $\Delta I$(%) | $\Delta {{I}_{c}}$(%) | $\Delta {{I}_{h}}$(%) | $\Delta {{\eta }_{c}}$(%) | $\Delta {{\eta }_{h}}$(%) |
---|---|---|---|---|---|---|---|---|---|---|---|
Annual | 452.1 | 418.6 | 456.4 | 0.478 | 0.352 | 0.494 | -12.6 | 1.6 | -14.2 | 10.1 | 89.9 |
January | 163.0 | 159.1 | 182.9 | 0.213 | 0.198 | 0.292 | -1.5 | 7.9 | -9.4 | 45.7 | 54.3 |
February | 143.2 | 140.5 | 153.3 | 0.235 | 0.221 | 0.284 | -1.4 | 4.9 | -6.3 | 43.8 | 56.3 |
March | 164.2 | 154.6 | 189.8 | 0.306 | 0.266 | 0.413 | -4.0 | 10.7 | -14.7 | 42.1 | 57.9 |
April | 246.7 | 257.4 | 397.6 | 0.185 | 0.205 | 0.464 | 2.0 | 27.9 | -25.9 | 51.9 | 48.1 |
May | 539.8 | 533.1 | 543.2 | 0.408 | 0.396 | 0.413 | -1.2 | 0.5 | -1.7 | 22.7 | 77.3 |
June | 684.4 | 774.5 | 613.6 | 0.393 | 0.517 | 0.295 | 12.4 | -9.8 | 22.2 | 30.6 | 69.4 |
July | 914.2 | 755.7 | 735.9 | 0.506 | 0.363 | 0.345 | -14.3 | -16.1 | 1.8 | 89.9 | 10.1 |
August | 689.3 | 602.1 | 716.3 | 0.425 | 0.325 | 0.456 | -10.0 | 3.1 | -13.1 | 19.1 | 80.9 |
September | 754.1 | 640.6 | 738.9 | 0.568 | 0.417 | 0.548 | -15.1 | -2.0 | -13.1 | 13.2 | 86.8 |
October | 592.6 | 524.5 | 603.1 | 0.461 | 0.350 | 0.478 | -11.1 | 1.7 | -12.8 | 11.7 | 88.3 |
November | 325.3 | 283.4 | 356.8 | 0.486 | 0.343 | 0.593 | -14.3 | 10.7 | -25.0 | 30.0 | 70.0 |
December | 209.2 | 197.6 | 245.7 | 0.410 | 0.344 | 0.618 | -6.6 | 20.8 | -27.4 | 43.2 | 56.8 |
QN1D | 136.9 | 112.6 | 138.7 | 0.596 | 0.400 | 0.610 | -19.6 | 1.4 | -21.0 | 6.3 | 93.8 |
QN3D | 137.7 | 115.1 | 140.0 | 0.584 | 0.404 | 0.602 | -18.0 | 1.8 | -19.8 | 8.3 | 91.7 |
QN7D | 138.9 | 118.6 | 142.5 | 0.566 | 0.409 | 0.593 | -15.7 | 2.7 | -18.4 | 12.8 | 87.2 |
QN30D | 142.2 | 132.2 | 150.8 | 0.375 | 0.302 | 0.439 | -7.3 | 6.4 | -13.7 | 31.8 | 68.2 |
QN90D | 158.3 | 152.2 | 180.2 | 0.346 | 0.312 | 0.470 | -3.4 | 12.4 | -15.8 | 44.0 | 56.0 |
QX1D | 2014.3 | 1868.7 | 1657.0 | 0.509 | 0.442 | 0.346 | -6.7 | -16.3 | 9.6 | 62.9 | 37.1 |
QX3D | 1711.8 | 1602.3 | 1538.0 | 0.450 | 0.393 | 0.359 | -5.7 | -9.1 | 3.4 | 72.8 | 27.2 |
QX7D | 1411.2 | 1338.8 | 1329.7 | 0.425 | 0.376 | 0.369 | -4.9 | -5.6 | 0.7 | 88.9 | 11.1 |
QX30D | 1091.5 | 1003.2 | 979.7 | 0.410 | 0.322 | 0.298 | -8.8 | -11.2 | 2.4 | 82.4 | 17.6 |
QX90D | 891.4 | 811.7 | 821.20 | 0.558 | 0.424 | 0.440 | -13.4 | -11.8 | -1.6 | 88.1 | 11.9 |
Base flow | 0.29 | 0.27 | 0.3 | 0.531 | 0.460 | 0.570 | -7.1 | 3.9 | -11.0 | 26.2 | 73.8 |
Date min | 56.7 | 72.5 | 56.5 | 0.077 | 0.124 | 0.076 | 4.7 | -0.1 | 4.8 | 2.0 | 98.0 |
Date max | 203.8 | 199.6 | 214.4 | 0.418 | 0.384 | 0.503 | -3.4 | 8.5 | -11.9 | 41.7 | 58.3 |
Lo pulse # | 2.1 | 2.6 | 1.9 | 0.302 | 0.369 | 0.276 | 6.7 | -2.6 | 9.3 | 21.8 | 78.2 |
Lo pulse L | 57.9 | 34.3 | 30.8 | 0.423 | 0.250 | 0.225 | -17.3 | -19.8 | 2.5 | 88.8 | 11.2 |
Hi pulse # | 11.3 | 12.0 | 10.2 | 0.431 | 0.473 | 0.367 | 4.2 | -6.4 | 10.6 | 37.6 | 62.4 |
Hi pulse L | 5.0 | 3.8 | 5.7 | 0.231 | 0.135 | 0.281 | -9.6 | 5.0 | -14.6 | 25.5 | 74.5 |
Rising rate | 36.6 | 23.3 | 26.9 | 0.714 | 0.381 | 0.472 | -33.3 | -24.2 | -9.1 | 72.7 | 27.3 |
Fall rate | -18.1 | -15.7 | -13.2 | 0.415 | 0.545 | 0.676 | 13.0 | 26.1 | -13.1 | 66.6 | 33.4 |
Reversals | 107.8 | 129.5 | 87.6 | 0.289 | 0.463 | 0.126 | 17.4 | -16.3 | 33.7 | 32.6 | 67.4 |
Mean (+) | - | - | - | - | - | - | 8.6 | 8.3 | 9.2 | 41.4 | 58.6 |
Mean (-) | - | - | - | - | - | - | -10.2 | -10.8 | -14.3 |
Figure 7 Comparison of different climatic variables and discharge (a) and Pearson correlation coefficient between different IHA parameters (b). In (a), $\overline{{{P}_{bp}}}$ and $\overline{{{P}_{ap}}}$ are the mean precipitation in the baseline and altered periods, respectively; $\overline{{{T}_{bp}}}$ and $\overline{{{T}_{ap}}}$ are the mean temperature in the baseline and altered periods, respectively; $\overline{E{{T}_{bp}}}$ and $\overline{E{{T}_{ap}}} $are the mean evapotranspiration in the baseline and altered periods, respectively; and $\overline{{{Q}_{bp}}}$ and $\overline{{{Q}_{ap}}}$ are the mean discharge in the baseline and altered periods, respectively. |
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|
[24] |
|
[25] |
|
[26] |
|
[27] |
|
[28] |
|
[29] |
|
[30] |
|
[31] |
|
[32] |
|
[33] |
|
[34] |
|
[35] |
|
[36] |
|
[37] |
|
[38] |
|
[39] |
|
[40] |
|
[41] |
|
[42] |
|
[43] |
|
[44] |
|
[45] |
|
[46] |
|
[47] |
|
[48] |
|
[49] |
|
[50] |
|
[51] |
|
[52] |
|
[53] |
|
[54] |
|
[55] |
|
[56] |
|
[57] |
|
[58] |
|
[59] |
|
[60] |
|
[61] |
|
[62] |
|
[63] |
|
[64] |
|
[65] |
|
[66] |
|
/
〈 | 〉 |