Next Article in Journal
Aerobiological Monitoring in an Indoor Occupational Setting Using a Real-Time Bioaerosol Sampler
Next Article in Special Issue
Flood Modelling and Risk Analysis of Cinan Feizuo Flood Protection Area, Huaihe River Basin
Previous Article in Journal
Assessment of Climatic Impact on Vegetation Spring Phenology in Northern China
Previous Article in Special Issue
Dry and Wet Changes and Vegetation Time-Delay Responses in Western China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of a Function to Fit the Flow Duration Curve and Parameterization of a Semi-Arid Region in North China

1
State Key Laboratory of Eco-Hydraulics in Northwest Arid Region of China, School of Water Resources and Hydropower, Xi’an University of Technology, Xi’an 710048, China
2
Key Laboratory of Water Cycle and Related Land Surface Processes, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
3
School of Water Conservancy and Hydroelectric Power, Hebei University of Engineering, Handan 056002, China
4
Key Laboratory of Water Management and Water Security for Yellow River Basin (Ministry of Water Resources), Yellow River Engineering Consulting Co., Ltd., Zhengzhou 450003, China
*
Author to whom correspondence should be addressed.
Atmosphere 2023, 14(1), 116; https://doi.org/10.3390/atmos14010116
Submission received: 15 December 2022 / Revised: 2 January 2023 / Accepted: 3 January 2023 / Published: 5 January 2023
(This article belongs to the Special Issue Regional Hydrological Processes in a Changing Climate)

Abstract

:
The discharge process has undergone major changes in many river basins throughout the world as a result of the simultaneous influences of global climate change and human activity. Flow duration curves (FDCs) are crucial indicators of river basins’ hydrological processes. However, it is challenging to compare FDCs in a quantitative way. This study will identify the best function with which to fit the flow duration curve in a semi-arid region of North China, so as to quantify the FDC, and parameterize the function of the FDC of the region in order to describe the FDCs of ungauged basins. In this work, six small- and medium-sized catchments in North China are selected as the study area, and three functions, i.e., log normal, generalized Pareto and H2018 functions, were chosen to fit the FDC at nineteen hydrological stations. The relationship between the parameters of the FDC and the basin characteristics, such as the climatic factors and geographical features, were analyzed. A regression formula of the parameters of the FDC function was established, and its spatial and temporal distributions were examined. Based on the evaluation of four indicators, the Nash–Sutcliffe efficiency, the root mean square relative error, the logarithmic Nash efficiency coefficient and the coefficient of determination, the results demonstrate that the H2018 function can match FDCs the best. Through the annual runoff, annual precipitation, precipitation in summer, potential evapotranspiration, catchment area, mean elevation, length of the main channel and maximum flow frequency, the parameters of a, b, and k in the H2018 function can be formulated. The regression formula constructed in this study can obtain a regional flow duration curve with satisfactory performance, which provides a reference for the validation of remote-sensing-based runoff data in ungauged regions.

1. Introduction

It is still a major challenge to scientifically and efficiently manage the water resources in basins, which is primarily due to a lack of hydrological data. The International Association of Hydrological Sciences (IAHS) put forward hydrological projections for basins without data, which are predictions in ungauged basins (PUBs) [1]. This is one of the key and difficult problems in hydrology research and has received extensive attention from many scholars [2,3,4,5]. In basins with no data or a lack of data, it is difficult to establish a plan for regional water resources.
However, there are only a certain number of gauging stations worldwide, and this number is decreasing as a result of funding reductions (i.e., not all streams are gauged). Generally, there are two methods for runoff process research in ungauged basin [6]: One is to establish a watershed hydrological model that only contains physical parameters. Since the model parameters may be measured, they can be used directly in the areas without data. The alternative option is to build a watershed hydrological model with parameters that will be calibrated, in which the model parameters of regions with data are simulated and applied to regions without data [7,8].
FDCs can capture the properties of precipitation and runoff by building a statistical relationship between them, and they have been widely applied in various fields of hydrology [9,10,11,12]. FDCs are built in two main steps: ranking the streamflow data in descending order and plotting the sorted values against the corresponding frequency of exceedance [13]. As a graphical representation of the relationship between flow frequency and flow rate can be simple and comprehensive, an altered graphic with which to describe the whole study period of runoff, from low water to flood characteristics of the traffic condition, can better reflect the basin rainfall runoff characteristics, and it also can be applied into water resources exploitation and protection.
A method with which to estimate the FDCs of ungauged sites based on distance measures that can be related to the catchment area and the climatic parameters has been established by many authors that addressed the topic of FDC prediction at ungauged or partially gauged locations through regional regression [14,15,16,17] and geostatistical interpolation [18]. Spatial nonlinear interpolation methods were developed by several scholars [19,20,21]. Worland et al. (2019) [22] presented a method involving the use of the copula function. Hughes and Smakhtin (1996) [23] proposed a method with which to extend or fill in daily flow time series at a site by using the monthly FDCs of a target site itself.
The Loess Plateau is located in the middle and upper reaches of the Yellow River, with fragile ecological environment and soil erosion [6,24]. With global climate change, the frequency of extreme climate events increases, and the risk of meteorological disasters intensifies. Furthermore, the rivers of the Loess Plateau are under the dual effects of global climate change and human activities; significant changes in flow processes have occurred, and the smallest watersheds are ungauged basins without observational flow data [25]. The National Development and Reform Commission of the Ministry of Water Resources officially issued and implemented the 14th five-year implementation plan for the construction of warping dams in the Yellow River basin and the comprehensive control of soil and water loss in sloping farmland [26], which intends to, over a period of five years, build 1461 warping dams and 2559 sand-blocking dams in the concentrated source area of coarse sediment. A large number of silt dams for soil and water conservation are constructed in small watersheds with no discharge data. If the FDCs of small watersheds can be constructed by regionalization of parameters, they will provide an important reference for the construction of warping dams in small watersheds.
The primary method used in the current study on FDCs is the construction of a trustworthy fitting function with which to infer the flow processes in an ungauged basin. The shapes of FDCs have strong regional differences. Preliminary studies on the Loess Plateau have shown that the low-flow part of the flow duration curve of the Yellow River basin will rapidly decrease and show an obvious S-shape. Blum et al. (2017) [27] omitted intermittent sites with an average daily flow value of zero from their analysis because such intermittent sites require additional methodological considerations. Therefore, it is necessary to optimize the function form of the optimal FDC on the Loess Plateau and to construct a calculation formula of the parameters. Based on this, this study proposed constructing an FDC, analyzing the influential elements of the curve shape, and studying the variation in the parameters of the FDC on both the temporal and spatial scales. Due to the influence of both global climate change and human activities, important signatures of the hydrological processes of river basins have changed significantly, especially flow duration curves (FDCs). However, FDCs are difficult to quantitatively compare and differ between different basins due to climatic and basin characteristics. The shapes of the curves vary greatly between basins. For instance, the research finds that FDCs display an ‘L-shape’ in the Americas [28]; however, the distribution runoff process in North China exhibits an ’S-shape’. With this being the case, the H2018 function [29] has been put forward to describe the FDC of this study area. Therefore, it is necessary to study the variation of the FDC of this basin and obtain FDCs in ungauged regions via parametric analyses without discharge data. Additionally, the parametric formula of the FDC can be applied to validate the remote-sensing-based runoff data in ungauged basins.
The study is organized as follows. Section 2 presents the study area and data. Section 3 presents the methodology employed to fit FDCs at hydrological stations. Section 4 presents the main results. Section 5 discusses the results in a catchment. Section 6 concludes the study and highlights the outlook for further research based on our finding.

2. Study Area and Data

2.1. Study Area

In this study, six small- and medium-sized watersheds of the Loess Plateau in North China were selected as the study watersheds, covering 58,677.54 km2, as shown in Figure 1. They include the Wuding River basin, Jialu River basin, Tuwei River basin, Kuye River basin, Gushan River basin, and Huangfu River basin.
The landform types of the northern Shaanxi region are mainly the aeolian sand transition area and the Loess Plateau area with sparse vegetation, and the region is one of the areas of China with serious levels of soil erosion and desertification [30,31]. It has a typical warm temperate continental monsoon climate, with a cold and dry winter, and a warm and rainy summer. The average annual temperature is 3.6~14.3 °C, and the annual rainfall is 300~800 mm; 65% of the precipitation occurs mainly in the summer (June to September) [32].

2.2. Data

In this study, the daily runoff data of the 19 hydrological stations in the six watersheds, as shown in Figure 1, were obtained from the hydrological data of the Yellow River. Due to the lengths of the observed data being inconsistent, the longest data sequence is from 1955 to 2012, and the shortest is from 1975 to 2012. The length of the data series is sufficient for this research. Furthermore, in order to explore the factors influencing FDCs, regional meteorological characteristics and watershed characteristics data were extracted, as shown in Table 1, including the sub-basin area covered by hydrological stations, A (km2), the average annual precipitation, P (mm), the average annual precipitation in summer, Psummer (mm), the average annual potential evapotranspiration in the sub-basins, Ep (mm), the annual runoff, R (108 m3), the mean elevation of the sub-basins, H (m), the length of the main channel, L, and the elevation difference of the sub-basins, ΔH (m) [33].
The DEM data used in this study were from the official website of the US Geological Survey, and which were extracted from the NASA Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model (ASTER GDEM) with a spatial resolution of 30 m. The precipitation data (1961–2015) were from China’s ground precipitation monthly value 0.5 × 0.5 grid data set (V2.0), from the China Meteorological Data Sharing Service Network (http://data.cma.cn) on 1 October 2022.

3. Methods

Three functions, i.e., log normal function, generalized Pareto function and H2018 function, were selected to simulate FDC through four assessment indices, i.e., the Nash–Sutcliffe efficiency, the root mean square relative error, the logarithmic Nash efficiency coefficient, and the coefficient of determination. Additionally, a regression model was then used to construct the formulae of the parameters.

3.1. Log Normal Function

If the function Y = lnX of the random variable x obeys the normal distribution N(μ, σ2), the probability density function of the lognormal distribution of x obeys the parameters μ and σ2 is said to be as follows:
f ( x ) = 1 x 2 π σ exp [ 1 2 σ 2 ( ln x μ ) 2 ]
where X is the random variable, and x is the independent variable of the probability density function, μ is the logarithmic mean parameter, and σ is the logarithmic variance parameter, as well as the scale parameter.

3.2. Generalized Pareto Function

If X obeys the generalized Pareto function G P D ( μ ,   σ ,   ξ ) , then its cumulative distribution function is as follows:
F ( μ , σ , ξ ) ( X ) = 1 ( 1 + ξ ( X μ ) σ ) 1 / ξ
where ξ 0 .
Its corresponding distribution density function (PDF) is as follows:
f ( μ , σ , ξ ) ( x ) = 1 σ ( 1 + ξ ( x μ ) σ ) ( 1 ξ 1 )

3.3. H2018 Function

Han and Tian [29] proposed a new function of the form below:
y = 1 1 + m x m a x x x x m i n n
where m and n are the parameters, x m i n x x m a x .
In flow duration curve fitting, when the frequency increases from 0 to 1, the flow gradually decreases, and the function in the original form needs to be transformed, such that the modified H2018 function form is as follows:
y = a 1 + b x x m a x x k

3.4. Evaluation Indices

In this study, the Nash–Sutcliffe efficiency (NSE) [17,34,35], the root mean square relative error (RMSRE) [36], the logarithmic Nash–Sutcliffe efficiency coefficient (LNSE) and the coefficient of determination (R2) [37,38,39] were used to evaluate the applicability of different distribution functions at the hydrological station.
The Nash efficiency coefficient (NSE) is an evaluation parameter used to evaluate model quality and is generally used to verify the quality of hydrological model simulation results. The NSE is defined as follows:
N S E = 1 t = 1 T ( Q o t Q m t ) 2 t = 1 T ( Q o t Q ¯ o ) 2
where Q o t refers to the observed value, Q m t refers to the simulated value, Q t represents a certain value at time t, and Q ¯ o represents the total average of the observed value.
The value of the NSE is from negative infinity to 1, and if the NSE is close to 1, this indicates good model quality and high model reliability. If the NSE is close to 0, this indicates that the simulation result is close to the average level of the observation value, that is, the overall result is reliable, but the process simulation error is large. If the NSE is much less than 0, the model is not credible.
The logarithmic Nash efficiency coefficient (LNSE) is an evaluation parameter calculated by taking the logarithms of observed and simulated values. The specific calculation formula is as follows:
L N S E = 1 t = 1 T ( ln Q o t ln Q m t ) 2 t = 1 T ( ln Q o t ln Q ¯ o ) 2
In the formula, ln Q o t refers to the logarithm of the observed value, ln Q m t refers to the logarithm of the simulated value, and ln Q ¯ o refers to the logarithm result of the total average of the observed value.
The root mean square relative error, also known as the standard error, is the square root of the ratio between the sum of squares of the deviation between the observed value and the truth value and the number of observations, m, which measures the deviation between the observed value and the truth value. RMSRE is calculated as follows:
R M S R E = 1 N t = 1 N ( Q o t Q m t Q o t ) 2
where N is the number of samples, Q o t is the observed value and Q m t is the simulated value.
The coefficient of determination (R2) is the square of the Pearson correlation coefficient and is a kind of non-deterministic relation, and is a quantity used to study the degree of linear correlation between variables.

3.5. Multiple Regression

Regressive methods have been used to link the different model parameters to some catchment characteristics such as climatic indices, land coverage, and geological and geomorphological parameters. This analysis was performed in this study by dividing the island into 19 subzones. The regressive formula used in this study has the following structure [40]:
y = k 0 + i = 1 n k i C i
where y represents the parameters used in the function which has been chosen, parameters k0 and ki are determined through a regression analysis.
Multiple regression examines how multiple independent variables are related to one dependent variable. Once each of the independent factors has been determined to predict the dependent variable, the information on the multiple variables can be used to create an accurate prediction on the level of effect that they have on the outcome variable.

3.6. Mann–Kendall Test

The Mann–Kendall test, one of the non-parametric statistical test procedures, is frequently utilized by domestic and international researchers due to its simplicity and effective application; it is particularly useful for the investigation of hydrology, meteorology, and other non-normal distribution data. In this study, the absolute value of U is larger than 1.645, which means that it has passed the significant two-sided test of 0.1 when the Mann–Kendall test is applied for the trend test, and the reliability is 0.1 [41,42].

4. Results

4.1. Best Fitting of the FDC

All of the observed flows at the 19 hydrological stations were used to fit the function for the FDC. Additionally, the best function was chosen from the log-normal distribution, generalized Pareto distribution (GPD), and H2018 distribution. Table 2 and Table 3 show the values of the LNSE, NSE, RMSRE, and R2 at each station, respectively.
For comparison, the LNSE had better performance than the H2018 simulation, with a maximum value of 0.9906 at the Shenmu station and a minimum value of 0.7685 at the Zhaoshiyao station. The highest values of the NSE were mostly in the H2018 simulation (the maximum value reached was 0.9515 at the Baijiachuan station). The RMSRE was lower in the simulation of H2018 and the GPD function compared to the log-normal function, and the lowest value occurred at the Wangdaohengta station. Additionally, the correlation coefficient performed well in both the H2018 and log-normal functions; the best one was 0.9800 at the Baijiachuan station. It is thus proven that the H2018 function was most effective in terms of the best fitting of the FDC through comprehensive judgment. Furthermore, Table 4 shows the best-fitted parameters of the H2018 function.
The stations with optimal parameter values of H2018 function fitting were then selected, and a fitting effect diagram of the three functions on FDCs was drawn, as shown in Figure 2 (the Baijiachuan, Shenmu, Wenjiachuan, Huangfu, and Hengshan stations). The Huangfu station obviously indicates that the maximum frequency was lower than 0.6, which was mainly triggered by a large number of days without runoff production. Overall, the H2018 function can fit the variation in the flow processes better, especially for the heads and tails of the observed curves.

4.2. Characteristics of Fitted Curves with H2018

The Baijiachuan station is taken as an example, and the variation characteristics of the H2018 curves in Figure 2a are analyzed as shown in Figure 3. In Table 5, the flow of ten percentiles, Q 1 , Q 5 , Q 10 , Q 20 , Q 50 , Q 70 , Q 80 , Q 90 , Q 95 , Q 99 represent the flow at the percentiles of 1%, 5%, 10%, 20%, 50%,70%, 80%, 90%, 95%, and 99%, respectively. Each flow percentile represents a different section of the FDC: extreme high flows (Q1 and Q5) to extreme low flows (Q95 and Q99) [43,44]. Low flows are of great importance for the functioning of many processes taking place in river ecosystems [45]. The low-flow thresholds adopted in drought studies usually vary between the 70th and 95th percentiles [46,47].
In order to better recognize the characteristics of fitting curves, the values of specific flow percentiles, and ratios of them, were calculated and are presented in Table 5. Q1/Q99 represents the total change in the curve, Q1/Q50 represents the degree of change in the front, and Q1/Q99 represents the degree of change in the tail.
The results show that the curve has a great change in its entirety. The high-flow part and the low-flow part of the curve both show a descending trend. The FDC presents an ‘S-shape’. Additionally, the results indicated that Q1 is 178.6667 times Q99. The change in the tail is larger than that in the front of the curve. In other words, low flows have greater variation. The above changes in the FDC indicate that the high flow is much higher than the low flow, and that the flow has significantly changed in one year. The extreme value of the flow varies greatly.

4.3. Formulization of Parameters

To investigate the variation in the parameters of the FDC, geomorphological and climatic factors of basins were taken into account. If the FDC can be signified by different climatic and basin characteristics, an inverse calculation will directly obtain the formula of the FDC through related factors. The main influence factors are the average annual P (mm), Psummer (mm), R (108 m3), Ep (mm), A (km2), ΔH (m), H (m), L (km), and x m a x . While coal mining has a great influence on the Kuye River, this study deleted the stations (Shenmu, Wangdaohengta, and Wenjiachuan) of this basin [48]) in order to obtain better function forms of the parameters. The results indicate that the parameters of the FDC have strong relationships with A, H, L, Ep, R, Psummer, and x m a x . The specific expressions calculated by multiple linear regression are as follows:
a = 1.1542 × 10 4 R 181.0585 L 716.4113 P 3.111 × 10 14 R 2 + 0.0025 × 10 4 A 2 + 0.7573 P 2 + 5.1899 L 2 + 0.0347 P A 0.0069 P L 0.3278 A L                                                                    
where the R2 = 0.9624, RMSRE = 0.0236, and P = 0.0058 < 0.01.
b = 45.9423 A 0.0173 R 612.7860 H + 0.0021 A 2 0.0337 R 2 + 4.0180 H 2 2.3569    × 10 7 A R 2.1412 × 10 4 A L 0.8274 A H + 2.3568 × 10 4 R L + 1.3409 × 10 5 R H
where the R2 = 0.9924, RMSRE = 0.0224, and P = 0.0044 < 0.01.
k = 26.981 x m a x 2.098 P s u m m e r + 0.339 P 25.05 x m a x 2 0.008 P s u m m e r 2 0.00014 P 2      + 2.413 x m a x P s u m m e r 0.388 x m a x P + 0.002 P s u m m e r P 3.604
where the R2 = 0.9305, RMSRE = 0.0494, and P = 0.0075 < 0.01.
The values of a, b, and k were calculated through the formulae constructed above (Table 6). Compared with the fitted parameter values, the simulated results are relatively similar (Figure 4). In the scatter diagrams, fitted values and simulated values all have a strong linear relationship, where the R2 values all exceed 0.9 and the RMSE values are all under 0.05 (Table 7). In other words, the simulations have better results.

4.4. The Spatial Distribution of Parameters

Although we calculated the above parameters in this study, how they are distributed in North China has not yet been elucidated. Therefore, the temporal and spatial variations in a, b, and k are analyzed in this part.
The specific spatial variation is simulated in Figure 5. It was found that the values in most stations are concentrated within 29,000–31,000. In addition to the great influence of basins on parameter a, the distribution of a/A was also considered. The value of b had a larger proportion from 2000 to 5000. Additionally, k was higher in the north of the basin by comparison.
In addition to the characteristics of the basin itself, Ep, R, and x m a x had relatively large impacts on the parameters’ distributions; a spatial distribution map of these three factors was drawn as Figure 5. Runoff was lower in the southeast, similar to parameter a. Ep was relatively large in the whole study area, and the majority of x m a x exceeded 0.9.

4.5. The Temporal Distribution of Annual FDC Parameters

As for temporal variation, the Baijiachuan station (which covers the entire Wuding River basin) was chosen to study the changes in parameters a, b, and k of the FDC based on annual flow data (FDCANN) (Figure 6). The trend variation in these three parameters can be found when the test statistic Z, calculated through a Mann–Kendall test, is equal to or greater than 1.96, indicating that it exceeded the significance test with 99% confidence. The results show that a has an obvious descending trend (Z = −4.8146), b has an obvious increasing trend (Z = 4.1371), and k has an unobvious increasing trend (Z = 1.4032).
According to the research findings, altering the parameter size scales and precisions of a will significantly modify the height of FDC, as shown in Figure 7. That is to say, a is mainly influenced by runoff. As can be seen from Figure 7, if only the value of a is changed, the height of the curves will shift significantly. Parameter a has changed from 8000 to 38,000 and 80,000 in Figure 7a–c, which illustrates that when a is higher the simulated curve will move upwards, and the Y-axis represents the volume of runoff. The result then shows that the change in a is greatly influenced by runoff. Based on this, the trend of runoff variation was analyzed. The runoff from 1975 to 2012 is similar to the variation in a, which also exhibits a declining tendency [49,50], and Z = −4.4253 via a Mann–Kendall test, which indicates a significant change.
Table 8 presents the characteristic values of a, b, and k, including the average, max and min of the parameters in FDCANN. As can be seen, the min of a occurred after the year 2000, while the max of b and k appeared. Comparatively speaking, there are only tiny differences between the max and min values of a (0.3%), whereas the values of b and k change significantly (67% and 73%, respectively).

5. Discussion

The FDC is investigated in this section due to the impact of coal mining and additional influences on the flow of the Kuye River. The parameters of FDC are quite different from the other river basins in this region. Figure 8, Figure 9 and Figure 10 indicate the temporal variations in three stations (Wenjiachuan, Wangdaohengta, and Shenmu) in the Kuye River. It can be easily found that both a and b have similar tendencies, and a Mann–Kendall test shows all obvious increases or decreases. However, the variation in parameter k is inconspicuous, and it has little difference with the other two stations in Wangdaohengta. Perhaps this was caused by the inconsistent geomorphology. The Wangdaohengta station, in particular, is situated in the upstream area and belongs to the sandy area, as opposed to the other two stations, which are situated in the Loess area. In addition, the annual runoff generally exhibited a downward tendency [31], the mutational point appeared in 1996, and coal output also started to increase at that time. Additionally, the upstream rainfall erosion is the smallest erosion area.
Based on the general survey data of water conservation of the Yellow River Conservancy Committee in 2011, there were 306 key dams in the Kuye River watershed, with a total storage capacity of 316.64 Mm3. Additionally, the majority of the dams were built in the mid-lower reaches. The regions that were at high risk of soil erosion and sediment yield were mainly concentrated in the middle reaches of the watershed [51]. The influence of human activities will affect the regional runoff variations to a large extent.
Furthermore, the characteristic values of a, b, and k at the Kuye River were also calculated (Table 9). The results indicated that the values of a are all smaller than those at the Baijiachuan station, while those of b and k are both larger. These values may affect the formulation of parameters in Section 4.3.

6. Conclusions

With the impact of human activity and climate change, runoff is continually changing, and the majority of the region’s minor watersheds are unmeasured, with insufficient flow data. Therefore, it is essential to develop formulas for FDCs that are widely applicable in the area and that may be expanded upon and used for research in a vast range of unmeasured areas. The study identifies the best function to fit the flow duration curve of a semi-arid region in North China, which is the H2018 function, and parameterizes the H2018 function of the FDCs of the region in order to estimate the FDCs of ungauged basins.
The generalized Pareto, H2018, and log-normal distribution functions are used to fit the flow duration curves of daily discharge at 19 stations in North China. The H2018 distribution function improves the tail and head fitting of the flow duration curve and has excellent performance to represent the flow duration curve, including zero-value discharge. The specific flows at percentiles were determined, and the ratio of the flow at different frequencies was calculated in order to better comprehend the changes of FDCs. As a result, the FDCs present an ‘S-shape’, in which the former and latter halves of the curve descend with large range changes, while the middle section has a smaller variation.
The parameters of a, b and k in H2018 were formulated with hydrometeorological features and basin characteristics, such as the annual runoff, precipitation, precipitation in summer, potential evapotranspiration, sub-basin area covered by hydrological stations, mean elevation, length of main channel, L, and max flow frequency xmax, by means of a regression model. The distributions of a, b, and k were all analyzed on both the spatial and temporal scales. Parameter a has an obvious descending trend, b has an obvious increasing trend, and k has an unobvious increasing trend. The regression formula constructed in this study can obtain a regional flow duration curve with satisfactory performance, which provides a reference for the validation of remote-sensing-based runoff data in ungauged regions.

Author Contributions

L.M.: Conceptualization, Methodology, Software, Validation, Visualization, Writing—original draft, Writing—review and editing. D.L.: Methodology, Software, Research Methods Provided, Writing—original draft, Supervision, Writing—review and editing, Funding acquisition, Resources. Q.H.: Supervision, review and modification. F.G.: Analysis of data, Writing—review and editing. X.Z.: Method instruction, Writing—review and editing. J.Z.: Method instruction, Writing—review and editing. J.L.: Writing—review and editing, Resources. J.F.: Writing–review and editing, Funding acquisition. G.M.: Writing—review and editing, Funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (NSFC) (Grant Nos. 42071335, 52279025, 52109031 and 52009102), and State Key Laboratory of Eco-hydraulics in Northwest Arid Region, China (Grant No. 2019KFKT-4).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The DEM data used in this study were from the official website of the US Geological Survey, and which were extracted from the NASA Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model (ASTER GDEM). The precipitation data can be obtained from the China Mete-orological Data Sharing Service Network (http://data.cma.cn, accessed on 14 December 2022). The other data can be obtained from the corresponding author (Dengfeng Liu, [email protected]) upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hrachowitz, M.; Savenije, H.H.G.; Blöschl, G.; McDonnell, J.J.; Sivapalan, M.; Pomeroy, J.W.; Arheimer, B.; Blume, T.; Clark, M.P.; Ehret, U.; et al. A decade of Predictions in Ungauged Basins (PUB)—A review. Hydrol. Sci. J. 2013, 58, 1198–1255. [Google Scholar] [CrossRef]
  2. Young, A.R. Stream flow simulation within UK ungauged catchments using a daily rainfall-runoff model. J. Hydrol. 2006, 320, 155–172. [Google Scholar] [CrossRef]
  3. Teutschbein, C.; Grabs, T.; Laudon, H.; Karlsen, R.H.; Bishop, K. Simulating streamflow in ungauged basins under a changing climate: The importance of landscape characteristics. J. Hydrol. 2018, 561, 160–178. [Google Scholar] [CrossRef]
  4. Wolff, W.; Duarte, S.N. Toward geostatistical unbiased predictions of flow duration curves at ungauged basins. Adv. Water Resour. 2021, 152, 103915. [Google Scholar] [CrossRef]
  5. Grimaldi, S.; Volpi, E.; Langousis, A.; Papalexiou, S.M.; De Luca, D.L.; Piscopia, R.; Nerantzaki, S.D.; Papacharalampous, G.; Petroselli, A. Continuous hydrologic modelling for small and ungauged basins: A comparison of eight rainfall models for sub-daily runoff simulations. J. Hydrol. 2022, 610, 127866. [Google Scholar] [CrossRef]
  6. Yu, P.S.; Yang, T.C. Using synthetic flow duration curves for rainfall-runoff model calibration at ungauged sites. Hydrol. Process 2000, 14, 117–133. [Google Scholar] [CrossRef]
  7. Worland, S.C.; Steinschneider, S.; Asquith, W.; Knight, R.; Wieczorek, M. Prediction and inference of flow duration curves using multioutput neural networks. Water Resour. Res. 2019, 55, 6850–6868. [Google Scholar] [CrossRef] [Green Version]
  8. Kim, D.; Kaluarachchi, J. Predicting streamflows in snowmelt-driven watersheds using the flow duration curve method. Hydrol. Earth Syst. Sci. 2014, 18, 1679–1693. [Google Scholar] [CrossRef] [Green Version]
  9. Mimikou, M.; Kaemaki, S. Regionalization of flow duration characteristics. J. Hydrol. 1985, 82, 77–91. [Google Scholar] [CrossRef]
  10. Shao, Q.; Zhang, L.; Chen, Y.D.; Singh, V.P. A new method for modelling flow duration curves and predicting streamflow regimes under altered land-use conditions/une nouvelle méthode de modélisation des courbes de débits classés et de prévision des régimes d’écoulement sous conditions modifiées d’occupation du sol. Hydrol. Sci. J. 2009, 54, 606–622. [Google Scholar] [CrossRef]
  11. Booker, D.J.; Snelder, T.H. Comparing methods for estimating flow duration curves at ungauged sites. J. Hydrol. 2012, 434–435, 78–94. [Google Scholar] [CrossRef]
  12. Majone, B.; Avesani, D.; Zulian, P.; Fiori, A.; Bellin, A. Analysis of high streamflow extremes in climate change studies: How do we calibrate hydrological models? Hydrol. Earth Syst. Sci. Discuss. 2021, 26, 3863–3883. [Google Scholar] [CrossRef]
  13. Ridolfi, E.; Kumar, H.; Bárdossy, A. A methodology to estimate flow duration curves at partially ungauged basins. Hydrol. Earth Syst. Sci. 2020, 24, 2043–2060. [Google Scholar] [CrossRef] [Green Version]
  14. Fennessey, N.; Vogel, R.M. Regional Flow-Duration Curvesfor Ungauged Sites in Massachusetts. J. Water Resour. Plan. Manag. 1990, 116, 530–549. [Google Scholar] [CrossRef]
  15. Rianna, M.; Efstratiadis, A.; Russo, F.; Napolitano, F.; Kout-soyiannis, D. A stochastic index method for calculating annual flow duration curves in intermittent rivers. Irrig. Draunage 2013, 62, 41–49. [Google Scholar] [CrossRef]
  16. Rianna, M.; Russo, F.; Napolitano, F. Stochastic indexmodel for intermittent regimes: From preliminary analysis toregionalisation. Nat. Hazards Earth Syst. Sci. 2011, 11, 1189–1203. [Google Scholar] [CrossRef]
  17. Pugliese, A.; Farmer, W.H.; Castellarin, A.; Archfield, S.A.; Vogel, R.M. Regional flow duration curves: Geostatistical techniques versus multivariate regression. Adv. Water Resour. 2016, 96, 11–22. [Google Scholar] [CrossRef]
  18. Ganora, D.; Claps, P.; Laio, F.; Viglione, A. An approach to estimate nonparametric flow duration curves in ungauged basins. Water Resour. Res. 2009, 45, W10418. [Google Scholar] [CrossRef]
  19. Archfield, S.A.; Vogel, R.M. Map correlation method: Selection of a reference streamgage to estimate daily streamflow at ungaged catchments. Water Resour. Res. 2010, 46, W10513. [Google Scholar] [CrossRef]
  20. Mohamoud, Y.M. Prediction of daily flow duration curves and streamflow for ungauged catchments using regional flow duration curves. Hydrol. Sci. J. 2008, 53, 706–724. [Google Scholar] [CrossRef]
  21. Farmer, W.H.; Archfield, S.A.; Over, T.M.; Hay, L.E.; LaFontaine, J.H.; Kiang, J.E. A Comparison of Methods to Predict Historical Daily Streamflow Time Series in the Southeastern United States; US Geological Survey: Reston, VA, USA, 2014; 34p. [CrossRef]
  22. Worland, S.C.; Steinschneider, S.; Farmer, W.; Asquith, W.; Knight, R. Copula Theory as a Generalized Framework for Flow-Duration Curve Based Streamflow Estimates in Ungaged and Partially Gaged Catchments. Water Resour. Res. 2019, 55, 9378–9397. [Google Scholar] [CrossRef] [Green Version]
  23. Hughes, D.A.; Smakhtin, V. Daily flow time series patching or extension: A spatial interpolation approach based on flow duration curves. Hydrol. Sci. J. 1996, 41, 851–871. [Google Scholar] [CrossRef] [Green Version]
  24. Fu, B.; Wang, S.; Liu Yu Liu, J.; Liang, W.; Miao, C. Hydrogeomorphic ecosystem responses to natural and anthropogenic changes in the Loess Plateau of China. Annu. Rev. Earth Planet. Sci. 2017, 45, 223–243. [Google Scholar] [CrossRef]
  25. Wang, Z.; Xu, M.; Liu, X.; Singh, D.K.; Fu, X. Quantifying the impact of climate change and anthropogenic activities on runoff and sediment load reduction in a typical Loess Plateau watershed. J. Hydrol. Reg. Stud. 2022, 39, 2214–5818. [Google Scholar] [CrossRef]
  26. The Ministry of Water Resources; the National Development and Reform Commission. Implemented the “14th Five-Year Plan” for the Construction of Silt Dams in the Yellow River Basin and the Comprehensive Control of Soil Erosion on Sloping Farmland; Ministry of Water Resources of the People’s Republic of China: Beijing, China, 2021.
  27. Blum Annalise, G.; Archfield Stacey, A.; Vogel Richard, M. On the probability distribution of daily streamflow in the United States. Hydrol. Earth Syst. Sci. 2017, 21, 3093–3103. [Google Scholar] [CrossRef] [Green Version]
  28. Westerberg, I.K.; Guerrero, J.L.; Younger, P.M.; Beven, K.J.; Seibert, J.; Halldin, S.; Freer, J.E.; Xu, C.-Y. Calibration of hydrological models using flow-duration curves. Hydrol. Earth Syst. Sci. 2011, 15, 2205–2227. [Google Scholar] [CrossRef] [Green Version]
  29. Han, S.; Tian, F. Derivation of a sigmoid generalized complementary function for evaporation with physical constraints. Water Resour. Res. 2018, 54, 1734–1736. [Google Scholar] [CrossRef]
  30. Zhao, G.; Mi, X.; Wen, Z.; Wang, F.; Gao, P. Soil erosion, conservation, and eco-environment changes in the loess plateau of china. Land Degrad. Dev. 2013, 24, 499–510. [Google Scholar] [CrossRef]
  31. Song, J.; Yang, Z.; Xia, J.; Cheng, D. The impact of mining-related human activities on runoff in northern Shaanxi, China. J. Hydrol. 2021, 126235, 598. [Google Scholar] [CrossRef]
  32. Li, C.; Wu, P.T.; Li, X.L.; Zhou, T.W.; Sun, S.K.; Wang, Y.B.; Luan, X.B.; Yu, X. Spatial and temporal evolution of climatic factors and its impacts on potential evapotranspiration in Loess Plateau of Northern Shaanxi, China. Sci. Total Environ. 2017, 589, 165–172. [Google Scholar] [CrossRef]
  33. Luan, J.; Liu, D.; Lin, M.; Huang, Q. The construction of the flow duration curve and the regionalization parameters analysis in the northwest of china. J. Water Clim. Change 2021, 12, 2639–2653. [Google Scholar] [CrossRef]
  34. Nash, I.E.; Sutcliffe, I.V. River flow forecasting through conceptual models part i—A discussion of principles-sciencedirect. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  35. Lane, P.N.J.; Best, A.E.; Hickel, K.; Zhang, L. The response of flow duration curves to afforestation. J. Hydrol. 2005, 310, 253–265. [Google Scholar] [CrossRef]
  36. Nruthya, K.; Srinivas, V.V. Evaluating methods to predict streamflow at ungauged sites using regional flow duration curves: A case study. Aquat. Procedia 2015, 4, 641–648. [Google Scholar] [CrossRef]
  37. Nagelkerke, N.J.D. A note on a general definition of the coefficient of determination. Biometrika 1991, 78, 691–692. [Google Scholar] [CrossRef]
  38. Madsen, H. Automatic calibration of a conceptual rainfall–runoff model using multiple objectives. J. Hydrol. 2000, 235, 276–288. [Google Scholar] [CrossRef]
  39. Shafii, M.; De Smedt, F. Multi-objective calibration of a distributed hydrological model (WetSpa) using a genetic algorithm. Hydrol. Earth Syst. Sci. 2009, 13, 2137–2149. [Google Scholar] [CrossRef] [Green Version]
  40. Viola, F.; Noto, L.V.; Cannarozzo, M.; La Loggia, G. Regional flow duration curves for ungauged sites in sicily. Hydrol. Earth Syst. Sci. 2011, 15, 323–331. [Google Scholar] [CrossRef] [Green Version]
  41. Abebe, S.A.; Qin, T.; Zhang, X.; Yan, D. Wavelet transform-based trend analysis of streamflow and precipitation in Upper Blue Nile River basin. J. Hydrol. Reg. Stud. 2022, 44, 2214–5818. [Google Scholar] [CrossRef]
  42. Güçlü, Y.S. Improved visualization for trend analysis by comparing with classical Mann-Kendall test and ITA. J. Hydrol. 2020, 584, 124674. [Google Scholar] [CrossRef]
  43. Panthi, J.; Talchabhadel, R.; Ghimire, G.R.; Sharma, S.; Dahal, P.; Baniya, R.; Boving, T.; Pradhanang, S.M.; Parajuli, B. Hydrologic regionalization under data scarcity: Implications for streamflow prediction. J. Hydrol. Eng. 2021, 26, 05021022. [Google Scholar] [CrossRef]
  44. Fouad, G.; Loáiciga, H.A. Independent variable selection for regression modeling of the flow duration curve for ungauged basins in the united states. J. Hydrol. 2020, 587, 124975. [Google Scholar] [CrossRef]
  45. Garbe, J.; Beevers, L.; Pender, G. The interaction of low flow conditions and spawning brown trout (Salmo trutta) habitat aviability. Ecol. Eng. 2016, 88, 53–63. [Google Scholar] [CrossRef] [Green Version]
  46. Hasan, H.H.; Razali, S.F.M.; Muhammad, N.S.; Hamzah, F.M. Assessment of probability distributions and analysis of the minimum storage draft rate in the equatorial region. Nat. Hazards Earth Syst. Sci. 2021, 21, 1–19. [Google Scholar] [CrossRef]
  47. Tomaszewski, E.; Kubiak-Wójcicka, K. Low-flows in Polish Rivers. In Management of Water Resources in Poland; Springer: Cham, Switzerland, 2021; pp. 205–228. [Google Scholar] [CrossRef]
  48. Luan, J.; Zhang, Y.; Tian, J.; Meresa, H.; Liu, D. Coal mining impacts on catchment runoff. J. Hydrol. 2020, 589, 125101. [Google Scholar] [CrossRef]
  49. Tian, L.; Jin, J.; Wu, P.; Niu, G.Y. Quantifying the impact of climate change and human activities on streamflow in a semi-arid watershed with the budyko equation incorporating dynamic vegetation information. Water 2018, 10, 1781. [Google Scholar] [CrossRef] [Green Version]
  50. Liang, W.; Bai, D.; Jin, Z.; You, Y.; Li, J.; Yang, Y. A study on the streamflow change and its relationship with climate change and ecological restoration measures in a sediment concentrated region in the loess plateau, China. Water Resour. Manag. 2015, 29, 4045–4060. [Google Scholar] [CrossRef]
  51. Xiang, Z.; Ds, A.; Mh, C.; Gw, D.; Yi, L.A. Understanding the influencing factors (precipitation variation, land use changes and check dams) and mechanisms controlling changes in the sediment load of a typical loess watershed, China. Ecol. Eng. 2021, 163, 106198. [Google Scholar] [CrossRef]
Figure 1. The distribution of the six watersheds and nineteen hydrological stations.
Figure 1. The distribution of the six watersheds and nineteen hydrological stations.
Atmosphere 14 00116 g001
Figure 2. Observed and simulated flow duration curves.
Figure 2. Observed and simulated flow duration curves.
Atmosphere 14 00116 g002aAtmosphere 14 00116 g002b
Figure 3. Flow values corresponding to different frequencies. Red circles represent the point of flow percentile.
Figure 3. Flow values corresponding to different frequencies. Red circles represent the point of flow percentile.
Atmosphere 14 00116 g003
Figure 4. Scatter diagrams comparing fitted values and simulated values of (a) Parameter a, (b) Parameter b and (c) Parameter k. Red circles represent the points of fitted values and simulated values and dotted lines represent linear fitting.
Figure 4. Scatter diagrams comparing fitted values and simulated values of (a) Parameter a, (b) Parameter b and (c) Parameter k. Red circles represent the points of fitted values and simulated values and dotted lines represent linear fitting.
Atmosphere 14 00116 g004aAtmosphere 14 00116 g004b
Figure 5. Spatial variation of (a) Parameter a, (b) a/A, (c) Parameter b, (d) Parameter k, (e) Ep, (f) R and (g) xmax.
Figure 5. Spatial variation of (a) Parameter a, (b) a/A, (c) Parameter b, (d) Parameter k, (e) Ep, (f) R and (g) xmax.
Atmosphere 14 00116 g005aAtmosphere 14 00116 g005b
Figure 6. Temporal variation of parameters of (a) Parameter a, (b) Parameter b and (c) Parameter k at the Baijiachuan station. The solid lines represent the value of parameters and dotted lines represent the linear fitting.
Figure 6. Temporal variation of parameters of (a) Parameter a, (b) Parameter b and (c) Parameter k at the Baijiachuan station. The solid lines represent the value of parameters and dotted lines represent the linear fitting.
Atmosphere 14 00116 g006aAtmosphere 14 00116 g006b
Figure 7. Change in the value of parameter a of FDCANN.
Figure 7. Change in the value of parameter a of FDCANN.
Atmosphere 14 00116 g007
Figure 8. Temporal variation of parameters of (a) Parameter a, (b) Parameter b and (c) Parameter k at the Shenmu station. The solid lines represent the value of parameters and dotted lines represent the linear fitting.
Figure 8. Temporal variation of parameters of (a) Parameter a, (b) Parameter b and (c) Parameter k at the Shenmu station. The solid lines represent the value of parameters and dotted lines represent the linear fitting.
Atmosphere 14 00116 g008aAtmosphere 14 00116 g008b
Figure 9. Temporal variation of parameters of (a) Parameter a, (b) Parameter b and (c) Parameter k at the Wangdaohengta station. The solid lines represent the value of parameters and dotted lines represent the linear fitting.
Figure 9. Temporal variation of parameters of (a) Parameter a, (b) Parameter b and (c) Parameter k at the Wangdaohengta station. The solid lines represent the value of parameters and dotted lines represent the linear fitting.
Atmosphere 14 00116 g009aAtmosphere 14 00116 g009b
Figure 10. Temporal variation of parameters of (a) Parameter a, (b) Parameter b and (c) Parameter k at the Wenjiachuan station. The solid lines represent the value of parameters and dotted lines represent the linear fitting.
Figure 10. Temporal variation of parameters of (a) Parameter a, (b) Parameter b and (c) Parameter k at the Wenjiachuan station. The solid lines represent the value of parameters and dotted lines represent the linear fitting.
Atmosphere 14 00116 g010
Table 1. Geomorphological and climatic factors for the catchments of hydrological stations.
Table 1. Geomorphological and climatic factors for the catchments of hydrological stations.
Hydrological StationA (km2)P (mm)Psummer (mm)Ep (mm)R (108 m3)L (km)ΔH (m)H (m)
Huangfu3279.66405.0369.661028.391.2639115.08036281159.31
Qingyangcha673.27495.0990.801041.210.276926.23075301378.66
Shenjiawan1147.43444.3377.141042.190.606763.07076391125.8
Suide3902.10471.7185.081136.61.4241140.25389011202.32
Dianshi467.09435.2778.251143.290.192332.64633831172.77
Hanjiamao2574.20349.5163.021195.340.8387109.47244281271.81
Dingjiagou30,123.08347.6967.691168.878.6968312.913010581291.17
Gaojiachuan3505.83427.8373.341117.673.3177122.44486381171.09
Lijiahe821.04441.0479.021150.790.339450.41804731183.01
Mahuyu376.91449.8580.011124.490.623232.64463991098.68
Gaojiabu2138.81417.7571.321147.512.556868.58254081223.6
Gaoshiya1274.63426.2772.72961.290.473465.14145801176.28
Hengshan2750.67433.1579.751084.730.6441130.14037991385.72
Baijiachuan36,485.42371.1271.001167.539.4079371.976412061265.55
Caoping205.52450.1480.051186.350.068217.59443081075.25
Zhaoshiyao22,138.10368.0267.881185.084.5628234.65519461329.69
Shenmu6930.95390.0765.301338.474.5651139.39446281300.3
Wenjiachuan12,517.48278.6946.881293.975.1815206.07308181260.1
Wangdaohengta3833.95378.2862.941431.291.8385119.46005571332.49
Table 2. The values of the LNSE and NSE.
Table 2. The values of the LNSE and NSE.
StationLNSENSE
Log NormalGPDH2018Log NormalGPDH2018
Huangfu−0.13540.45160.9657−6.0465 × 10310.19400.7390
Qingyangcha0.02530.13540.8442−9.7243 × 1013−5433.40000.5098
Shenjiawan0.41100.56230.9496−5.1788 × 1010−0.00470.6080
Shenmu0.31010.96300.9906−9.2092 × 10160.53040.6072
Wenjiachuan0.33250.71580.9693−6.9900 × 1012−0.00010.8152
Suide0.30660.91800.7868−4.0192 × 106−4.94910.5575
Wangdaohengta0.42160.90030.9700−3.7119 × 1012−0.03340.6194
Dianshi0.62380.51220.9517−1.5541 × 103−0.00010.0091
Hanjiamao−4.58910.80200.7833−1.1153 × 10180.06720.4122
Dingjiagou0.15580.81620.9075−80.45910.04010.4045
Gaojiachuan−14.77430.13370.9449−8.2734 × 10350.00090.0904
Lijiahe0.27030.77190.9588−2.6045 × 1060.00000.1052
Mahuyu0.46040.60210.9482−6.0553 × 106−0.00010.0060
Gaojiabu0.2834−1.20870.9711−1.3548 × 1030.00330.0311
Gaoshiya−3.69810.77770.9193−4.2702 × 10750.64760.7357
Hengshan−0.03180.71060.9852−1.9970 × 1090.01200.3849
Baijiachuan0.21360.84570.9004−5.5778 × 1060.14650.9515
Caoping0.50300.56820.9647−2.5749 × 108−0.00480.5006
Zhaoshiyao−0.84010.81730.76856.0750 × 10220.01460.1287
Table 3. The values of the RMSRE and R2.
Table 3. The values of the RMSRE and R2.
StationRMSRER2
Log NormalGPDH2018Log NormalGPDH2018
Huangfu1.68 × 10140.01940.01100.53370.45690.9078
Qingyangcha7.9230 × 1040.59230.00560.76190.88960.7640
Shenjiawan3.0760 × 1030.01350.00850.74610.10480.9348
Shenmu3.4500 × 1060.00780.00710.61540.74460.9340
Wenjiachuan4.6400 × 1040.01750.00750.46110.24870.9493
Suide16.31300.01980.00540.75240.90070.7868
Wangdaohengta1.5051 × 1040.00790.00480.76120.92460.7929
Dianshi0.28280.00720.00710.76410.01730.3347
Hanjiamao1.4565 × 1070.01330.01060.72140.27540.6483
Dingjiagou0.08540.00930.00730.94640.21970.6731
Gaojiachuan6.6080 × 10150.00730.00690.78550.06950.3794
Lijiahe12.05380.00750.00710.85570.03430.6102
Mahuyu19.40720.00790.00790.89730.01120.4910
Gaojiabu0.28700.00780.00770.72430.07110.2363
Gaoshiya8.75 × 10350.00800.00690.5539-0.8954
Hengshan3.6975 × 1020.00820.00650.85220.16680.6987
Baijiachuan56.29530.02200.00530.53500.44910.9800
Caoping2.0976 × 1020.01310.00920.64670.13500.9192
Zhaoshiyao3.0736 × 1060.00720.00680.79990.41000.4076
Table 4. The best-fitted parameters of the H2018 function.
Table 4. The best-fitted parameters of the H2018 function.
Stationabk       x m a x
Huangfu2.2906 × 1042.5500 × 1×104 0.85050.5771
Qingyangcha1.2597 × 1041.9400 × 1040.55951.0000
Shenjiawan0.1295 × 1040.1241 × 1040.59540.9856
Shenmu28.5827 × 1044.9721 × 1040.70211.0000
Wenjiachuan26.5870 × 1044.2650 × 1040.62450.9839
Suide2.0009 × 1040.8104 × 1040.62510.9919
Wangdaohengta2.5995 × 1041.2002 × 1040.76990.9885
Dianshi1.0043 × 1043.9879 × 1040.41990.9991
Hanjiamao3.0510 × 1041.2216 × 1040.31991.0000
Dingjiagou3.0000 × 1040.1351 × 1040.44540.9987
Gaojiachuan2.0902 × 1040.2300 × 1040.31120.9576
Lijiahe1.7906 × 1043.9500 × 1040.60450.9950
Mahuyu2.0824 × 1047.0983 × 1040.52800.9645
Gaojiabu2.9873 × 1040.4179 × 1040.23341.0000
Gaoshiya2.7352 × 1044.8779 × 1040.80700.8287
Hengshan2.9181 × 1041.9514 × 1040.44681.0000
Baijiachuan3.8000 × 1040.1691 × 1040.44040.9989
Caoping0.2003 × 1041.9222 × 1040.55540.9576
Zhaoshiyao2.9940 × 1040.2684 × 1040.44990.9697
Table 5. Values of specific flows at the percentiles at Baijiachuan.
Table 5. Values of specific flows at the percentiles at Baijiachuan.
Variables of the PercentilesValues
lnQ15.1401
lnQ54.4085
lnQ104.0764
lnQ203.7158
lnQ503.0986
lnQ702.7208
lnQ802.4800
lnQ902.1163
lnQ951.7783
lnQ990.9998
Q1/Q505.2344
Q1/Q99178.6667
Q50/Q9934.1333
Table 6. Comparison of calculation values.
Table 6. Comparison of calculation values.
StationFitted aSimulated aFitted bSimulated bFitted kSimulated k
Huangfu2.2114 × 1042.2906 × 1042.5500 × 1042.5500 × 1040.85110.8505
Qingyangcha2.6906 × 1042.7727 × 1045.2430 × 1045.1764 × 1040.47380.4367
Shenjiawan2.8116 × 1042.8703 × 1043.7916 × 1043.0109 × 1040.55450.5574
Suide2.0703 × 1042.0009 × 1041.1624 × 1040.8104 × 1040.53350.6251
Dianshi2.4002 ×1042.5243 × 10410.6696 × 10410.6878 × 1040.48960.4180
Hanjiamao3.0545 × 1043.0510 × 1041.2410 × 1041.2216 × 1040.30650.3199
Dingjiagou3.0385 × 1043.0000 × 1040.1692 × 1040.1351 × 1040.44850.4454
Gaojiachuan2.0136 × 1042.0909 × 104−0.1668 × 1040.2230 × 1040.31160.3112
Lijiahe2.6322 × 1042.7192 × 1046.7931 × 1046.3613 × 1040.50940.5625
Mahuyu2.3317 × 1042.0824 × 1045.8293 ×1047.0983 × 1040.52180.5280
Gaojiabu3.0067 × 1042.9873 × 1040.8134 × 1040.4179 × 1040.25130.2334
Gaoshiya2.8487 × 1042.7352 × 1044.7378 × 1044.8779 × 1040.77880.8070
Hengshan2.9428 × 1042.9181 × 1041.5422 × 1041.9514 × 1040.41290.4468
Baijiachuan3.7741 × 1043.8000 × 1040.2347 ×1040.1691 × 1040.44170.4404
Caoping2.1555 × 1042.1393 ×10419.3197 ×10419.2220 × 1040.64430.5554
Zhaoshiyao2.9940 × 1042.9940 × 1040.2082 × 1040.2684 × 1040.45820.4499
Table 7. Comparison of fitted values and simulated values.
Table 7. Comparison of fitted values and simulated values.
Indexabk
R20.96240.99240.9305
RMSRE0.02290.04240.0494
Table 8. Characteristic values of a, b, and k at the Baijiachuan station.
Table 8. Characteristic values of a, b, and k at the Baijiachuan station.
ParameterAverageMaxMin
a2.9997 × 1043.0029 × 1042.9937 × 104
b0.1415 × 1040.2732 × 1040.0901 × 104
k0.40891.00800.2671
Table 9. Characteristic values of a, b, and k at the Kuye River.
Table 9. Characteristic values of a, b, and k at the Kuye River.
StationParameterAverageMaxMin
Wangdaohengtaa2.9414 × 1042.9859 × 1042.8174 × 104
b1.4340 × 1044.1827 × 1040.4485 × 104
k0.58810.85160.3231
Shenmua2.9785 × 1043.0046 × 1042.9322 × 104
b0.6122 × 1041.6347 × 1040.0321 × 104
k0.59620.90390.0031
Wenjiachuana2.9833 × 1042.9982 × 1042.9480 × 104
b0.5049 × 1041.2862 × 1040.1756 × 104
k0.60510.82800.2301
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ma, L.; Liu, D.; Huang, Q.; Guo, F.; Zheng, X.; Zhao, J.; Luan, J.; Fan, J.; Ming, G. Identification of a Function to Fit the Flow Duration Curve and Parameterization of a Semi-Arid Region in North China. Atmosphere 2023, 14, 116. https://doi.org/10.3390/atmos14010116

AMA Style

Ma L, Liu D, Huang Q, Guo F, Zheng X, Zhao J, Luan J, Fan J, Ming G. Identification of a Function to Fit the Flow Duration Curve and Parameterization of a Semi-Arid Region in North China. Atmosphere. 2023; 14(1):116. https://doi.org/10.3390/atmos14010116

Chicago/Turabian Style

Ma, Lan, Dengfeng Liu, Qiang Huang, Fengnian Guo, Xudong Zheng, Jing Zhao, Jinkai Luan, Jingjing Fan, and Guanghui Ming. 2023. "Identification of a Function to Fit the Flow Duration Curve and Parameterization of a Semi-Arid Region in North China" Atmosphere 14, no. 1: 116. https://doi.org/10.3390/atmos14010116

APA Style

Ma, L., Liu, D., Huang, Q., Guo, F., Zheng, X., Zhao, J., Luan, J., Fan, J., & Ming, G. (2023). Identification of a Function to Fit the Flow Duration Curve and Parameterization of a Semi-Arid Region in North China. Atmosphere, 14(1), 116. https://doi.org/10.3390/atmos14010116

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop