Assessing Water Withdrawals in Scarce-Data Transboundary Areas by Use of Dynamic Precipitation–Flow Relationships: The Case of the Hasbani River Basin

: The present study proposes a nonstandard solution to the problem of assessing water withdrawals (AWW) in the scarce-data transboundary basin. The applied AWW method operates with the open-source available data on precipitation and river ﬂow and thereby overcomes the usual restriction due to lack of data on shared water use in the Middle East. Analysis of dynamic precipitation-ﬂow relationships enable to separate the effect of water withdrawals from the total decline of river ﬂow under the decreasing precipitation. This study is the ﬁrst which provides complete information (1972–2020) on water withdrawals from total, surface, and base ﬂow of the Hasbani River (Lebanon). The resulting values that exceed by far earlier published estimates were conﬁrmed by (i) indirect indices (area of irrigated land and population), and (ii) validation of the AWW method based on independent data on water use (Israel Water Authority). The study results are useful for water balance estimations, as well as for management of water resources in the Jordan River headwaters basin and in the entire Lake Kinneret Basin. The AWW method can be applied to other transboundary basins and enables historical and real-time monitoring of water withdrawals as a necessary database for settlement of riparian water relations.


Introduction
In the Middle East, most of the river flow originates outside the countries, and transboundary waters represent over two-thirds of the overall region water resources [1,2]. In a transboundary basin, no riparian state can independently dispose of its water, and all countries are compelled to consider each other's water rights and demands. Under grave shortage of water resources and uncertainty of regular water supply, such situation inevitably causes tension and discord in riparian water relations. It also implies the need for agreements and cooperation on water issues despite seemingly insurmountable geopolitical, socio-economic, and management problems and the negative historical experience of such efforts. It is evident that absence of the complete picture of water-related dynamics in transboundary basins and the lack of a basin-wide vision over the shared waters may result in disputes among the riparian states and ultimately in an unsustainable and biased use of the resources. There is no doubt on the importance and need of reliable data for making water-use decisions, but unfortunately reality is different. The following examples (randomly selected from many similar ones) illustrate the real situation: "The absence of official data (from Syria) on the amount of water diverted from the Yarmouk has left much room for speculation over the years" [1], or "The water balance scheme in the Jordan River Basin is based on the available data to date, which were provided from the different sources and not believed to be so accurate or highly consistent" [3], etc. Such situations are due to strained international relations between the riparian states in the transboundary basin caused by internal economic, political, and national interests. All this may result in secrecy and even distortion of the data on water use [4][5][6][7]. However, in many known cases (see below in the present paper), the situation was aggravated not only and not so much by the availability (or secrecy) of data but because of the real lack or absence of data. Both of these circumstances determine and confirm the relevance of the proposed study topic which is assessing water withdrawals in the transboundary basin by using alternative data on precipitation and river flow.
The studied Hasbani River (named Snir in Israel) is one of the major tributaries of the Upper Jordan River (UJR) which contributes more than a quarter of the total inflow to it ( Figure 1). The catchment area of the transboundary Hasbani basin is shared between Lebanon and Israel at respective rates of 97% and 3%. Hence, most of the flow formation and water use occur in the Lebanese part of the basin. The present study revealed significant decrease (about 50%) of the Hasbani River annual flow during 1940-2017 and only 13-15% decrease of the annual precipitation gauged in this period at the representative stations. Such significant flow decline of the Hasbani River cannot be caused only by the precipitation decrease but seems to be also the result of the developing water use.
Water 2021, 13, x FOR PEER REVIEW 2 of 34 balance scheme in the Jordan River Basin is based on the available data to date, which were provided from the different sources and not believed to be so accurate or highly consistent" [3], etc. Such situations are due to strained international relations between the riparian states in the transboundary basin caused by internal economic, political, and national interests. All this may result in secrecy and even distortion of the data on water use [4][5][6][7]. However, in many known cases (see below in the present paper), the situation was aggravated not only and not so much by the availability (or secrecy) of data but because of the real lack or absence of data. Both of these circumstances determine and confirm the relevance of the proposed study topic which is assessing water withdrawals in the transboundary basin by using alternative data on precipitation and river flow.
The studied Hasbani River (named Snir in Israel) is one of the major tributaries of the Upper Jordan River (UJR) which contributes more than a quarter of the total inflow to it ( Figure 1). The catchment area of the transboundary Hasbani basin is shared between Lebanon and Israel at respective rates of 97% and 3%. Hence, most of the flow formation and water use occur in the Lebanese part of the basin. The present study revealed significant decrease (about 50%) of the Hasbani River annual flow during 1940-2017 and only 13-15% decrease of the annual precipitation gauged in this period at the representative stations. Such significant flow decline of the Hasbani River cannot be caused only by the precipitation decrease but seems to be also the result of the developing water use. A drastic decrease in the flow of rivers in Lebanon has been recorded during the last three decades [8]. Decreases of spring and stream flow were documented both in the Jordan River basin [9][10][11], as well as in the Litani basin in Lebanon [12]. Data on available water Water 2021, 13, 1440 3 of 32 (AW) in Lake Kinneret ) indicated a decreasing trend which was concurrent with the decreasing trend in precipitation and the increasing water use in the entire Lake watershed. According to the Israeli Water Authority, the increasing water use explains at least three quarters of the total decrease in the AW during this period [10,13]. Concurrently, statistical analyses of rainfall time series in Israel have not detected significant decrease in rainfall trends [14,15].
There is no major disagreement that both climate change and increasing human activity impact the water balance of the Jordan River headwaters and of the entire Lake Kinneret Basin. However, the relative contribution of each factor remains to be debated especially considering the unspecified Lebanese water use from the Hasbani River [9,13,[16][17][18][19]. The published values of the Hasbani water use by Lebanon are rather uncertain and insignificant (see Section 6.1 below). The question remains open whether these estimated values are sufficient to explain the significant decrease in the flow of the Hasbani River (1940-2017), even taking into account the observed decrease in precipitation.
Due to insufficient data on water use in the Hasbani River basin, it is evident that (i) any clarification and updating of these data would be an essential contribution to assessing the water balance in the UJR basin, and (ii) the problem requires an alternative approach and a nonstandard solution allowing to differentiate climate and human impacts on the river flow regime. For this purpose the present study applies the method proposed by the authors [20] for assessing water consumption in a transboundary basin under conditions of scarce available data on water use. The method is based on analyzing the alternative data on precipitation and river flow and was efficiently applied to the Yarmouk River  up to the Maqaren Station (Syrian part of the basin). In the present paper, the method was conceptually revised taking into consideration the outcomes of recent widespread discussions of the difference between "water consumption" and "water supply" (use, withdrawal, abstraction). By definitions of the USGS Glossary [21], the object of this study is accepted as assessing water withdrawals, but not water consumption (as in the previous publication [20]). The title of the present article has been updated accordingly. Nevertheless, the calculation technique remained unchanged. Therefore, withdrawals are estimated as the total amount of water removed from the source (or did not reach the downstream hydrometric station and not included in the gauged flow), regardless of how much of that total is consumptively (or non-consumptively) used [22]. The designation AWW (assessing water withdrawals) of the method is hereinafter used.
Assessing the impact of water use on water resources is one of the urgent problems of modern hydrology. The publications offer a huge number of approaches, methods, and models for its solution in general and for particular cases. These studies apply nonparametric methods to analyze the changing point and annual trend in various hydrometeorological time series (e.g., rainfall, temperature, flow and groundwater levels), and afterwards, they employ an integrated model of surface-groundwater flow to simulate effects of alternative water-use scenarios on groundwater levels and streamflow. The differences are in kind and scale of the model, research goal, and results. Some typical examples of similar studies are presented in Appendix A.
In contrast to the modeling flow with consideration of water use as an influencing driver (in the studies presented above), the AWW method offers an inverse solution to the problem. The latter assesses water withdrawal using data on the basin input and output (measured precipitation and flow) within the simple statistical model. This is the novelty of the method and its specificity.
The main objective of the present study can be formulated as (i) Assessing water withdrawals from the Hasbani river flow (Lebanon) by the AWW method; (ii) validation of the AWW method and confirmation of reliability of the obtained results.
The paper is built as following: Section 2 presents a brief description of the study area; Section 3 presents the available data and introduces the AWW method which is adopted for the case study in Section 4; Section 5 presents a wide range of results such as estimated withdrawals from the Hasbani River for the study period, validation of the AWW method, confirmation of estimates by indirect indices, and estimated withdrawals from surface and base flow; in Section 6 results obtained are discussed; conclusions and recommendations are provided in Section 7.
It is important to emphasize that the use of the Hasbani river flow in Lebanon has been for decades a source of tension and of discord between Lebanon and Israel. To the present day there is no agreement on the use of these transboundary waters. The proposed study focuses on hydrological aspects of the problem (specifically, on assessing water withdrawals from the Hasbani River by alternative hydrometric and precipitation data), disregarding geopolitical, socio-economic, and management issues, and without touching upon the issue of transboundary water relations and disagreement between both riparian countries.

Study Area
The Hasbani (Snir) River is one of the major tributaries of the Upper Jordan River (UJR) (Figure 1). It is a significant flow source, which contributes 26% of the total inflow to UJR while the remaining flow is formed by the Dan and Hermon Rivers [10,16,17]. In the present study, contribution of the Hasbani River updated to 30%.
The main sources of the Jordan River are located in the southern foothills of the Hermon Range, which is an elongated anticline, mostly built of >2000 m Jurassic karstic limestone. The Mt. Hermon range is 55 km long and 25 km wide. Its high regions (above 1000 m a.s.l., the summit 2814 m a.s.l.) receive the highest annual precipitation (>1400 mm/year).
The Hasbani River originates from the Hasbani springs (outflowing from karst in the northern part of the basin), flows southwards for about 20 km until reaching the point at which its flow is substantially increased by the Wazzani spring. About 0.5 km downstream the river reaches the Lebanese-Israeli border, flows for about 4 km along the border, and after 7 km joins with Hermon (Banias) and Dan Rivers in northern Israel near Sede Nehemya (Yosef Bridge), to form the Upper Jordan River. The surface catchment area of the transboundary Hasbani basin (614 sq. km) is shared between Lebanon and Israel as 97% and 3%, respectively. Most of the flow is formed and used in the Lebanese part of the basin.
The climate in southern Lebanon is Mediterranean. Precipitation over the Hasbani Basin is characterized by a wet (rainy) season from October to April and by a dry hot season with relatively little to no precipitation from May to September. On Mount Hermon snow usually falls on the elevated areas from December to March, and persists on areas 1400-1900 m a.s.l. and above until March-June (depending on local conditions). The major replenishment of Hasbani River comes from precipitation, as well as from snowmelt and springs which are fed by groundwater from calcareous aquifers recharged by snow and rains. The Hasbani River derives most of its flow from the Wazzani and Hasbani karst springs (Lebanon) with the annual respective contribution of 45 and 30 MCM/y. The rest of the measured discharge is contributed by surface runoff of three tributaries and their seasonal rivulets or flow from other upstream springs. The contribution of the Wazzani spring is very important, since this is the only continuous year-round flow into the river. Water supplies to domestic and agricultural demands in the Hasbani Basin derive from springs, directly from the river and from public and private wells. In this section, origin of data is from [1,16,17,[23][24][25][26][27][28].

Data
Due to lack of sufficient and reliable data on the use of water resources in the transboundary Hasbani river basin, the alternative data on precipitation and river flow was used within the Assessing Water Withdrawal (AWW) method. These data are the annual series of flow of the Hasbani River recorded at the Mayan Barukh and Dan Road Bridge hydrometric stations (1940-2020) and annual series of precipitation gauged at Kfar Giladi (1940-2020) Expanding the study to include a separate assessment of withdrawals from surface and base flow required data on the monthly flow of the Hasbani River (1940-2020) recorded at the Mayan Barukh and Dan Road Bridge hydrometric stations (the Israel Hydrological Service).
It is important to notice that years 1940-2017 are taken as the study period, while series 2018-2020 are considered as independent data.
Validation of the AWW method required the following additional data: (i) annual series  of the Upper Jordan River flow recorded at Sede Nehemya Station (the Israel Hydrological Service) and (ii) data  supplied by the Israel Water Authority on water use in the Upper Catchments of the Jordan River (included the Dan and the Hermon basins in Israel but without the Hasbani basin in Lebanon (Figure 1)).

The AWW Method
The assessing water withdrawal (AWW) method focuses on the detection and attribution of the trends in historical time-series of flow and precipitation at the annual scale. According to [29], the term "Detection" implies that an observed change is significantly different (in a statistical sense) from what can be explained by natural internal variability. Hence, detection is primarily a statistical argument, without explaining the causes for change. The term "Attribution" is the process of establishing the most likely causes (drivers) for the detected change with defined levels of confidence.
At the stage of the trend detection, the AWW method tests the following hypotheses: 1. The significant decrease of linear trend is detected in historical time-series of flow and precipitation while the precipitation trend is less and cannot completely explain the decrease of flow.

2.
Dependence between flow and precipitation (PF relationship) is observed in the study period.

3.
During subsequent sub-periods, the dynamic PF relationship demonstrates the consistent flow decline when for the same precipitation in the basin, less and less runoff was measured at the downstream station.
The next step is the flow trend attribution. If the first and the second hypotheses are confirmed significantly, then a decrease in precipitation (climate signal) can be considered as an obvious (but not the only) driving force behind a decrease in flow. Confirmation of the third hypothesis indicates the influence of additional driver which is similar to the increased water withdrawal. Most likely, this effect is associated with human activities (development of irrigation, economic and domestic needs, etc.), but the influence of increase in total evaporation due to global air warming and other possible reasons is not excluded. The attribution of flow mode changes to specific drivers is discussed in Section 6.2. Until then, we use the term "withdrawals" without detailed interpretation, assuming that withdrawals are the total amount of water removed from the river (or did not reach it) and are not included in the measured flow.
If the third hypothesis is confirmed, the ultimate goal of the study is to assess the influencing factor (water withdrawals) by comparison of the measured annual flow with its "natural" value. The latter is calculated by the PF relationship (natural curve) for the earlier (unregulated) period. Solution of this problem would be trivial if the historical flow data are available (like the Hasbani flow series from 1939/40 to 1970/71 in the present study). Otherwise, the task becomes more difficult. In such case, it is proposed to construct a "natural curve" using the available historical data, which unfortunately may be very limited in transboundary basins [19,20].
On the graph of precipitation-flow relationships (see following figures) the natural curve should be the upper envelope of the field of points. The deviation of the corresponding point (measured yearly flow) downward from the natural PF curve can be regarded as an estimate of the annual volume of withdrawals. However, considering the statistical nature of the PF relationship (due to influence of other unaccounted factors), the corresponding sub-period curve is proposed to be used in the calculation procedure instead of the corresponding point. By this way, the annual volume of withdrawals is estimated as a distance between the natural PF curve and the corresponding sub-period PF curve and as a function of the yearly precipitation: where W is water withdrawals, Q ntr is the natural flow, Q k is flow calculated by the k-th dynamic PF relationship (in accordance with sub-period), X is precipitation, and j is the time step number (year). Thus, the AWW method evaluates water withdrawals as a current shortage (deficit) of flow, which in an earlier "natural" period could have been formed by the same precipitation in the basin, but during the period under consideration it did not reach the gauging station downstream due to upper water use and other potential factors. This formulation is consistent with the interpretation of this term above.
It is important to emphasize that the AWW method enables to separate the influence of climatic precipitation changes on the flow from the anthropogenic impacts (and other potential drivers). In the process of computing water withdrawals (as purpose of the study), the precipitation change is implicitly taken into account by the method definition, because dynamics of the PF relationship (but not declined flow trend) are involved in analysis.
Finally it should be noted that precipitation-flow (PF) relationships are widely used in hydrological research for solving problems of water forecasts and water resources management, for assessing hydrological effects of the climate change, and for other purposes. In each case, the necessary initial step is to determine the PF relationship for its further use with a specific purpose. Wherein input scenarios and output results can be temporally varied, PF relationships are supposed to remain unchanged (e.g., [9,10]). Unlike this approach, the AWW method is based on the retroactive analysis of the dynamic PF relationship which may change primarily as a result of human activity in the basin.
The practical application of the AWW method is described in detail in the following Section 4. Particular attention should be given to Appendix B, which presents estimates of the statistical significance of all calculated trends (Mann-Kendall Test) and correlation coefficients (p-Value Test for Pearson R Score) in this study.

Adaptation to the Case-Study
There is a long and complete series of the Hasbani River flow measured at the Israel hydrometric stations (Mayan Barukh and Dan Road Bridge) since hydrological year 1939/40. The stations are located in close vicinity near the Lebanon-Israel border ( Figure 1). For integrated series 1940-2017 (see Section 4.2), the mean annual flow of the Hasbani River was evaluated as 112 MCM/y but during successive 26-year sub-periods (1940-1965, 1966-1991 and 1992-2017) this value decreased respectively from 128 to 121 and 88 MCM/y. For the last sub-period the mean flow differs significantly from the previous two (p-value is less than 0.05) [30]. Figure 2 demonstrates the significant decrease (55%) of the Hasbani River annual flow during 1940-2017 and only 13% decrease of the annual precipitation recorded for this period at the representative Kfar Giladi station. By the Mann-Kendall Test (Appendix B, Table A1), detected linear trends are statistically significant (with significance level alpha 0.05 and 0.10 for flow and precipitation, respectively). Comparison of trends indicates that the decrease of precipitation is insufficient to explain the considerable flow decline of the Hasbani River, and the flow change seems to be caused by the added drivers such as the developing water use or others. flow during 1940-2017 and only 13% decrease of the annual precipitation recorded for this period at the representative Kfar Giladi station. By the Mann-Kendall Test (Appendix B, Table A1), detected linear trends are statistically significant (with significance level alpha 0.05 and 0.10 for flow and precipitation, respectively). Comparison of trends indicates that the decrease of precipitation is insufficient to explain the considerable flow decline of the Hasbani River, and the flow change seems to be caused by the added drivers such as the developing water use or others. Prior to averaging (AVG Rain), the depths of annual precipitation at Kfar Giladi and Golan Experimental stations are z-score standardized, i.e., transformed to values with a mean of zero and a standard deviation of 1. The z-score standardized variable is the distance between data point and the sample mean in relation to the sample standard deviation. In application of statistics, variables measured on different scales or on a common scale with widely differing ranges are often adjusted by this way for comparison or averaging [31,32].
It is important to notice that according to the regionalization of the rainfall field in the Upper Jordan River (Lake Kinneret Basin) [33], the selected Kfar Giladi and Golan Experimental stations represent different regions: (R1) Upper Galilee and Western Hasbani Basin and (R3) Eastern Hasbani Basin, Northern Golan and Hermon, respectively. Here regions are named as in the referenced paper, in which the reginal division was done considering the specific relationships between the mean annual precipitation depth and elevation. Prior to averaging (AVG Rain), the depths of annual precipitation at Kfar Giladi and Golan Experimental stations are z-score standardized, i.e., transformed to values with a mean of zero and a standard deviation of 1. The z-score standardized variable is the distance between data point and the sample mean in relation to the sample standard deviation. In application of statistics, variables measured on different scales or on a common scale with widely differing ranges are often adjusted by this way for comparison or averaging [31,32].
It is important to notice that according to the regionalization of the rainfall field in the Upper Jordan River (Lake Kinneret Basin) [33], the selected Kfar Giladi and Golan Experimental stations represent different regions: (R1) Upper Galilee and Western Hasbani Basin and (R3) Eastern Hasbani Basin, Northern Golan and Hermon, respectively. Here regions are named as in the referenced paper, in which the reginal division was done considering the specific relationships between the mean annual precipitation depth and elevation.
As a result, the present study operates with the quite long precipitation data series , and AVG Rain is a statistical characteristic of the precipitation field which is more representative and sustained (against random errors) in comparison with the point data. The efficiency of used approach is confirmed by statistics in Table 1. The Hasbani River flow at Mayan Barukh Station in comparison with precipitation (as AVG Rain) during the study period (1940-2017) is presented in Figure 3. The Hasbani River flow at Mayan Barukh Station in comparison with precipitation (as AVG Rain) during the study period (1940-2017) is presented in Figure 3.

Natural Flow Curve
The initial necessary step is to establish the Natural Flow Curve by analyses of PF relationship during the earlier unregulated period  considering the history of the hydrometric stations on the Hasbani River (see map in Figure 1). According to the Hydrological Year-Book of Israel [34], during 1940-1951 the Hasbani (Snir) flow was gauged at the Dan Road Bridge hydrometrical station. In 1952 a part of the Dan River water was diverted through the Qoren Channel into the Snir River upstream of the gauging station. Monthly and yearly flow figures for 1952-1962 represent measured values after deduction of the Qoren input. In 1963 the new Mayan Barukh station (with basin area reduction of less than 1%) started to work upstream of the existing Qoren and Dan Channels (the latter dug in 1962). However, due to technical difficulties in the operation of the new station, the previous downstream Dan Road Bridge Station continued to serve as supporting for flood measurements, and this situation persists to the present date. Since 1963 to the present time the Hydrological Year-Books of Israel [35] include data on the Snir (Hasbani) flow at both stations (with relevant comments), while the present study uses flow measured at the Mayan Barukh station. Figure 4 reveals that during the initial (unregulated) period  there was a tight dependence between Snir (Hasbani) flow and precipitation. The uniformity of the

Natural Flow Curve
The initial necessary step is to establish the Natural Flow Curve by analyses of PF relationship during the earlier unregulated period    Table A2). Points (and curve 2) corresponding to measurements at the Dan Road Bridge (1952)(1953)(1954)(1955)(1956)(1957)(1958)(1959)(1960)(1961)(1962) are placed mostly lower and were not taken into account. In all probability, this was due to above mentioned hydrometrical problems in this period connected with the inaccurate measurement of the water diversion from Dan to Snir.
Water 2021, 13, 1440 9 of 32 period without 1952-1962 years) is accepted as the Natural Flow Curve of the Snir (Hasbani) River with the significant PF correlation coefficient 0.952 (Appendix B, Table  A2). Points (and curve 2) corresponding to measurements at the Dan Road Bridge (1952)(1953)(1954)(1955)(1956)(1957)(1958)(1959)(1960)(1961)(1962) are placed mostly lower and were not taken into account. In all probability, this was due to above mentioned hydrometrical problems in this period connected with the inaccurate measurement of the water diversion from Dan to Snir.

Dynamic PF Relationship
Figure 5a demonstrates that the established Natural Flow Curve for the Snir (Hasbani) River is the upper envelope of the entire field of points related to the following 1972-2017 period. During this period there is a high accuracy dependence (with correlation coefficient 0.918) between flow and precipitation (flow rises as rainfall increases) while the curve is placed lower than the Natural Flow Curve and confirms the effect of the water withdrawals.
A comprehensive analysis of the relationship between precipitation and flow for 1972-2017 enabled to identify two specific sub-periods within this period ( Figure 5b). The first sub-period of 26 years (1972-1997) stands out by the higher accuracy dependence of flow on precipitation (correlation coefficient 0.969) than during the whole period. In the following sub-period (1998-2017) the curve is placed lower and correlation coefficient (0.915) persists almost at the level of the overall PF relationship. All correlation coefficients are statistically significant (Appendix B, Table A2).
Further analysis showed that a more detailed presentation of the PF relationship is possible, but limited by the available data series and the statistical significance of the estimates.
Generally, during successive sub-periods (1972-1997 and 1998-2017) consistent flow decline is observed. This means that for the following sub-period, the same precip-   . On the other hand, the significant decrease of the river flow use (to the level of the initial unregulated period and of the previous sub-period, respectively) in these rainy years could be caused by the geopolitical and socio-economic (regional or local) factors. Though discussion of possible hypotheses is outside the scope of the present study, the authors consider it important (i) to show a possibility of such situation in transboundary basin and (ii) to raise this issue for a special research.

Dynamic PF Relationship
When the years 1940-1972 are used as the study period, the independent data deserve particular attention. These data are precipitation gauged at Kfar Giladi and Golan Experimental stations in years 2018-2020 and the Hasbani flow recorded at the Mayan Barukh station in these years. The respective points are located near the lower curve as continuation of the 1998-2017 sub-period ( Figure 5b). It means that the independent data support the dynamic character of the PF relationship.

Estimated Withdrawals from the Hasbani River
According to the AWW method, the natural PF curve (Figure 5b) serves as a basis A comprehensive analysis of the relationship between precipitation and flow for 1972-2017 enabled to identify two specific sub-periods within this period (Figure 5b). The first sub-period of 26 years (1972-1997) stands out by the higher accuracy dependence of flow on precipitation (correlation coefficient 0.969) than during the whole period. In the following sub-period (1998-2017) the curve is placed lower and correlation coefficient (0.915) persists almost at the level of the overall PF relationship. All correlation coefficients are statistically significant (Appendix B, Table A2).
Further analysis showed that a more detailed presentation of the PF relationship is possible, but limited by the available data series and the statistical significance of the estimates.
Generally, during successive sub-periods (1972-1997 and 1998-2017) consistent flow decline is observed. This means that for the following sub-period, the same precipitation in the basin formed less flow (measured at the Mayan Barukh station) because of increasing water withdrawals.
The years 2003-2004 and 2012-2013 are exceptional because of the excessive annual flow as compared to the PF relationship for 1998-2017 sub-period. Detailed analysis did not show any natural cause for such sharp rise of the runoff coefficient which was observed both for surface and base flow too (see Section 5.4.2). On the other hand, the significant decrease of the river flow use (to the level of the initial unregulated period and of the previous sub-period, respectively) in these rainy years could be caused by the geopolitical and socio-economic (regional or local) factors. Though discussion of possible hypotheses is outside the scope of the present study, the authors consider it important (i) to show a possibility of such situation in transboundary basin and (ii) to raise this issue for a special research.
When the years 1940-1972 are used as the study period, the independent data deserve particular attention. These data are precipitation gauged at Kfar Giladi and Golan Experimental stations in years 2018-2020 and the Hasbani flow recorded at the Mayan Barukh station in these years. The respective points are located near the lower curve as continuation of the 1998-2017 sub-period (Figure 5b). It means that the independent data support the dynamic character of the PF relationship.

Estimated Withdrawals from the Hasbani River
According to the AWW method, the natural PF curve (Figure 5b) serves as a basis for assessing the water withdrawals  in the Hasbani Basin (Lebanon). The annual volume of withdrawals is estimated by Equation (1) Table A1). Figure 6 demonstrates that against the background of decreased measured flow, increased withdrawals make more and more part of the natural flow. Figure 7 presents the verification of the AWW method for the Hasbani River. This graph shows the relationship between the measured Hasbani flow and precipitation   6a. Natural flow of the Hasbani River as a sum of the measured flow and estimated water withdrawals (1972-2020) is presented in Figure 6b. Figure 6a demonstrates the significant decrease (45%) of the Hasbani River annual flow during 1972-2020 and sharper increase (113%) of the water withdrawals. By the Mann-Kendall Test, both detected linear trends are statistically significant (Appendix B, Table A1). Figure 6 demonstrates that against the background of decreased measured flow, increased withdrawals make more and more part of the natural flow.   Summing up, Table 2 presents the comparison of the mean annual values of the measured Hasbani flow and the calculated water withdrawals during successive 12(13) years sub-periods (1972-1983, 1984-1995, 1996-2007, and 2008-2020). These data indicate the following: -Measured flow was almost unchanged (112-113 MCM) during the two initial sub-periods but dropped sharply to 88 and 66 MCM in the third and the fourth Summing up, Table 2 presents the comparison of the mean annual values of the measured Hasbani flow and the calculated water withdrawals during successive 12 (13) years sub-periods (1972-1983, 1984-1995, 1996-2007, and 2008-2020). These data indicate the following: -Measured flow was almost unchanged (112-113 MCM) during the two initial subperiods but dropped sharply to 88 and 66 MCM in the third and the fourth sub-periods (decrease by 21% and 41%, respectively). -Contrarily, water withdrawals increased from 30-31 MCM during the two initial sub-periods to 38 and 48 MCM in the third and the fourth sub-periods (increase by 27% and 60%, respectively). -As a result, the portion of the water withdrawals within the natural flow doubled and reached 42% in 2008-2020. Note to Table 2: Water withdrawals (%) is related to the natural flow (sum of measured flow and withdrawals).

Validation of the AWW Method
Until now, the AWW method has not been validated in any case of its applicationneither for the Yarmouk Basin [20] nor for the Hasbani Basin (in the present paper). In both cases, this is due to a lack of independent and reasonably complete data on water use in the study basin. Such data (1975-2014) on water use in the Upper Catchments of the Jordan River (including the Dan and the Hermon basins in Israel but without the Hasbani basin in Lebanon) were supplied by the Israel Water Authority [36].
To apply this data it is necessary to expand the study area to the basin of the Upper Jordan River up to Sede Nehemya Station. This area (named the Upper Catchments of the Jordan River (UCJR)) includes watersheds of the major headwater sources-Hasbani (Snir), Dan and Hermon Rivers (see map in Figure 1). The area of the transboundary UCJR basin covers in total 860 sq. km which is shared (as 70% and 30%, respectively) by Lebanon (Hasbani basin) and Israel (other basins).
The mean annual flow of the Upper Jordan River at Sede Nehemya Station during 1943-2017 is evaluated (in the present study) as 382 MCM/y. However, for three 25-years sub-periods (1943-1967, 1968-1992, and 1993-2017) this value is reduced from 401 to 394 and 355 MCM/y, respectively. Figure 8 demonstrates the apparent decrease (19%) in flow of the Upper Jordan River measured at the Sede Nehemya Station during hydrological years 1943-2017 and only 12% decrease in precipitation gauged during this period at the representative Kfar Giladi station. These data indicate that the precipitation decrease is insufficient to explain the flow decline of the Upper Jordan River. The previously estimated withdrawals from the Hasbani River in Lebanon and data on water use supplied by the Israel Water Authority corroborate the water use in the UCJR basin. Figure 9 shows the water withdrawals in the Hasbani Basin (Lebanon) calculated by the AWW method vs. water use in the UCJR Basin (Israel) during 1975-2014 years. It is worth recalling that 30% of the total flow of the Upper Jordan River is formed in the Hasbani basin, while the remaining 70% is the contribution of the Dan and Hermon Rivers. The obvious (oppositely directed) linear trends were detected for both variables: the water withdrawals in the Hasbani Basin increased by 46% and water use in the UCJR Basin (Israel) decreased by 39%. Figure 8 demonstrates the apparent decrease (19%) in flow of the Upper Jordan River measured at the Sede Nehemya Station during hydrological years 1943-2017 and only 12% decrease in precipitation gauged during this period at the representative Kfar Giladi station. These data indicate that the precipitation decrease is insufficient to explain the flow decline of the Upper Jordan River. The previously estimated withdrawals from the Hasbani River in Lebanon and data on water use supplied by the Israel Water Authority corroborate the water use in the UCJR basin.    This is the first opportunity to compare water use in the Lebanese and Israeli parts of the UCJR Basin. As seen, until 1997 the water withdrawals in the Hasbani Basin varied from half to three quarters of the water use in the UCJR Basin (Israel). Further the water use in the Hasbani Basin increased and reached (and even exceed) the Israel water use, particularly considering the reduction of the latter following the drought in 1999-2000 (by the operative decision of the Israel Water Authority). In 2012-2014, water withdrawals in Lebanon decreased in contrast to its growth in Israel.
Summary for 10-years sub-periods is presented in Table 3. As seen, during four sub-periods the Hasbani withdrawals increased from 30 to 45 MCM (one-third to half of the total SUM). The opposite situation was observed for water use in UCJR (Israel): it dropped and in 2005-2014 approached the Hasbani withdrawals. Ratio between these components was up from 0.47 to 0.74. The main conclusion from Figure 9 and Table 3 is the availability of data on water use in the UCJR Basin shared between Lebanon and Israel (1975-2014): calculated (by This is the first opportunity to compare water use in the Lebanese and Israeli parts of the UCJR Basin. As seen, until 1997 the water withdrawals in the Hasbani Basin varied from half to three quarters of the water use in the UCJR Basin (Israel). Further the water use in the Hasbani Basin increased and reached (and even exceed) the Israel water use, particularly considering the reduction of the latter following the drought in 1999-2000 (by the operative decision of the Israel Water Authority). In 2012-2014, water withdrawals in Lebanon decreased in contrast to its growth in Israel.
Summary for 10-years sub-periods is presented in Table 3. As seen, during four subperiods the Hasbani withdrawals increased from 30 to 45 MCM (one-third to half of the total SUM). The opposite situation was observed for water use in UCJR (Israel): it dropped and in 2005-2014 approached the Hasbani withdrawals. Ratio between these components was up from 0.47 to 0.74.
The main conclusion from Figure 9 and Table 3 is the availability of data on water use in the UCJR Basin shared between Lebanon andIsrael (1975-2014): calculated (by the AWW method) values for Hasbani Basin and independent data of the Israel Water Authority for the remaining basin area. The total water use (SUM) was compared with measured flow of the Upper Jordan River at Sede Nehemya Station ( Figure 10). Insignificant negative trends are detected for both variables during 1975-2014.  In order to validate the AWW method, the present study proposes to test the following hypothesis (hereinafter called as "hypothesis of validation"): the natural flow of the Upper Jordan River at Sede Nehemya Station (1975-2014) calculated as sum of measured flow, withdrawals in the Hasbani basin and water use in the UCJR (Israel) is consistent with the earlier unregulated flow of this river. In the proposed design scheme, the Hasbani water withdrawals (in Lebanon) are the only estimated variable which has to be validated by applying independent data on measured flow of the Upper Jordan River and water use in the Israeli part of the UCJR basin.
Based on previous experience in this study, the Upper Jordan flow was analyzed by comparing it with precipitation. Despite availability of 75-years series of measured flow and precipitation and PF relationship with correlation coefficient 0.881 (Figure 11), selection of the history period with the unregulated flow of the Upper Jordan River at Sede Nehemya Station turned out to be problematic. The problem is that "from 1951/52 part of the flow was diverted from tributaries upstream of station (through western and eastern Hula channels) and the natural regimen of flow was disturbed" [34]. Due to lack of data on water use in these years, only analyses of the PF relationship made it possible to approximately select the sub-period with the relatively unregulated flow (1943)(1944)(1945)(1946)(1947)(1948)(1949)(1950)(1951)(1952)(1953)(1954)(1955) and to establish Natural Flow Curve (Figure 11). Despite the short 13-years series, the correlation coefficient (0.873) is statistically significant (Appendix B, Table A2) and is at the level of the overall PF relationship (0.881 for 1943-2017).
The established Natural Flow Curve is the upper envelope of the entire field of points related to the following 1956-2017 period. There is a high accuracy dependence of flow on precipitation in this period (flow rises as rainfall increases) with correlation coefficient 0.919 while the trend-line is placed lower than the Natural Flow Curve and confirms the water withdrawals. However, assessing water withdrawals in the UCJR by the AWW method is a special problem and is not a necessary part of this study. The necessary and sufficient result is establishing of Natural Flow Curve for the Upper Jordan River (at Sede Nehemya) during 1943-1955 years ( Figure 11). Figure 12 presents this Natural Flow Curve (1) in comparison with Natural Flow Curve (2) which is trend-line of PF relationship based on the reconstructed natural flow (1975-2014). The latter is cal- In order to validate the AWW method, the present study proposes to test the following hypothesis (hereinafter called as "hypothesis of validation"): the natural flow of the Upper Jordan River at Sede Nehemya Station (1975-2014) calculated as sum of measured flow, withdrawals in the Hasbani basin and water use in the UCJR (Israel) is consistent with the earlier unregulated flow of this river. In the proposed design scheme, the Hasbani water withdrawals (in Lebanon) are the only estimated variable which has to be validated by applying independent data on measured flow of the Upper Jordan River and water use in the Israeli part of the UCJR basin.
Based on previous experience in this study, the Upper Jordan flow was analyzed by comparing it with precipitation. Despite availability of 75-years series of measured flow and precipitation and PF relationship with correlation coefficient 0.881 (Figure 11), selection of the history period with the unregulated flow of the Upper Jordan River at Sede Nehemya Station turned out to be problematic. The problem is that "from 1951/52 part of the flow was diverted from tributaries upstream of station (through western and eastern Hula channels) and the natural regimen of flow was disturbed" [34]. Due to lack of data on water use in these years, only analyses of the PF relationship made it possible to approximately select the sub-period with the relatively unregulated flow (1943)(1944)(1945)(1946)(1947)(1948)(1949)(1950)(1951)(1952)(1953)(1954)(1955) and to establish Natural Flow Curve (Figure 11). Despite the short 13-years series, the correlation coefficient (0.873) is statistically significant (Appendix B, Table A2) and is at the level of the overall PF relationship (0.881 for 1943-2017).
The established Natural Flow Curve is the upper envelope of the entire field of points related to the following 1956-2017 period. There is a high accuracy dependence of flow on precipitation in this period (flow rises as rainfall increases) with correlation coefficient 0.919 while the trend-line is placed lower than the Natural Flow Curve and confirms the water withdrawals. However, assessing water withdrawals in the UCJR by the AWW method is a special problem and is not a necessary part of this study. The necessary and sufficient result is establishing of Natural Flow Curve for the Upper Jordan River (at Sede Nehemya) during 1943-1955 years ( Figure 11). Figure 12 presents this Natural Flow Curve (1) in comparison with Natural Flow Curve (2) which is trend-line of PF relationship based on the reconstructed natural flow . The latter is calculated as the sum of the Upper Jordan measured flow, withdrawals in Hasbani basin (Lebanon) and water use in UCJR (Israel).     Both Natural Flow Curves [NFC] turned out to be of the same order of magnitude. The resemblance between curves (1) and (2) is confirmation of the above formulated "hypothesis of validation". This result is verified numerically by comparison of the natural flow calculated by [NFC-1] and [NFC-2] for 10-year sub-periods (Table 4). The difference between two variances is less than 1% (relatively to the natural flow) and less than 5% (relatively to the total water use SUM in the UCJR). If the discrepancy between [NFC-1] and [NFC-2] is attributed only to the calculated withdrawals in the Hasbani Basin (Table 3), which is an unrealistically "strict" condition, then the difference increases to 8-10%.  (1) and (2) in Figure 12.
Thus, the sharing of independent data on water use in the upper Jordan River catchments in Israel [36] with estimates of withdrawals from the Hasbani River provided a unique opportunity to validate the AWW method. On the other hand, the results of the AWW method validation can be considered as a decisive argument in support of the estimated water withdrawals in the Hasbani Basin.

Confirmation of Results by Indirect Indices
Under conditions of data deficiency, it is all the more important to compare the obtained estimates of water withdrawals with the real water demand and water use in the Hasbani basin based on indirect indices such as area of the irrigated agricultural land and population. By FAO [8], irrigation water in Lebanon is estimated by the rate of about 19,780 cub m/ha/year. This rate (of about 2 MCM/sq. km per year) is confirmed by the currently discussed design of the 50 MCM/y capacities of the Ibl al Saqi Dam on the Hasbani River (Lebanon) with a potential irrigation area of about 2600 ha (26 sq. km) [1,37].
According to the 1999 census, potential agricultural land surface in the Hasbani basin amounted to around 30,000 ha. However, only 15,000 ha were cultivated and 1124 ha (~11 sq. km) were irrigated [1,38]. Another study based on remote sensing (2002) states that an area of 17,600 ha was used for agriculture, of which about 9150 ha (~91 sq. km) received full or supplementary irrigation [39]. Supplementary irrigation is defined as the addition of small amounts of water to essentially rain-fed crops in order to improve and stabilize yields when rainfall fails to provide sufficient moisture for normal plant growth. Based on the first (least) estimate of irrigated land (11 sq. km) and excluding supplementary irrigated area, the accepted irrigation water use in the Hasbani basin is not less than 22 MCM per year.
Data on domestic water use in Lebanon are contradictory. According to FAO [8], it is estimated as 220-250 L/capita/day during the dry period and 200 L/capita/day during the wet period. By an earlier estimate of the Ministry of Environment (Lebanon) this figure corresponds to 140 L/capita/day [40], which fits well with the World Bank estimate of 150 L/capita/day for municipal water use [41]. If the latter value is taken into account, then the domestic water rate can be roughly estimated at 55 cub m/ capita/year.
Estimates of the Lebanese population within the Hasbani basin vary from 95,000 to 170,000 [42,43] with the most reliable estimate of 135,000 [44,45]. According to these sources, the large variability of the estimated population reflects not only poor data quality, but also the high employment-related seasonal migration. Based on the foregoing values of the water use rate and population, the domestic water use in the Hasbani Basin is estimated as not less than 7 MCM per year.
It must be emphasized that the above presented estimates relate to the period previous to year 2011 in which the war in Syria erupted and large numbers of refugees started flowing into Lebanon. According to the UN Refugee Agency (UNHCR), by 2017 about 1 million Syrian refugees spread through the country and over 110,000 settled in Southern Lebanon, almost doubling the resident population within the Hasbani Basin [46]. Among other large economic and social problems, this influx of people had a severe negative impact on the environment and put water resources under severe stress (high demand coupled with increasing pollution). But there are no available estimates of the linkage between water resources use and the migration of refugees in Lebanon [47].
By the most moderate estimates of indirect indices (areas of irrigated land, population, and water consumption rates), the total (irrigation and domestic) water use in the Hasbani Basin (before 2011) should be at least 22 + 7 = 29 MCM/y. This value supports the estimation of 34 MCM/y (calculated water withdrawals for 1972-2010 in the present study) and exceeds the previously published estimations (7-11 MCM/y) of the water use in the Hasbani Basin (see Section 6.1 for details).

Assessing Withdrawals from Surface and Base Flow
This section stands out in the present paper for several reasons: - The main purpose of this study is to evaluate the total water withdrawals from the Hasbani River. The title and the content of the paper correspond to this goal. Nevertheless, the separate assessment of withdrawals from the surface water and groundwater (base) flow is of particular complementary interest, especially in the karst basin. -On the other hand, the object of research is the scarce-data transboundary basin. The limitation also applies to data on precipitation and flow, especially concerning a time scale. This study is only based on annual data, which is the advantage of the AWW method as being widely available. -Expanding the study to include assessing withdrawals from surface and base flow will undoubtedly require additional information which may not occur in transboundary basins. The challenge lies in choosing an appropriate flow separation method and making optimal use of the available data. Below are the results of a compromise solution to this problem for the Hasbani River.

Surface and Base Flow Separation
The hydrograph peaks of the Hasbani River reflect the flow responsiveness to rainfall, while the steady part corresponds to groundwater-fed base flow. The further necessary step would be division of the hydrograph into surface water (SW-quick rainfall floods) and groundwater flows (GW-background of long-term flow or BF-base flow). For this purpose there is a wide range of graphical, analytical, and tracer-based methods which are represented in many classical and modern references. Separating the base flow and quick flow components of hydrographs is a common undertaking in hydrology and is used as a part of rainfall-runoff modelling, assessments of groundwater recharge and discharge and flood estimation.
In the present study, the "recursive digital filter" method originating from signal processing is applied. This method, which is designed to separate high from low-frequency signals, is widely used for the hydrograph division to the surface water and base flow components, e.g., [48][49][50][51][52].
The algorithm for the digital filtering separation of the flow hydrograph is described by: where Q is the total flow, q is the quick flow component (high frequency signal), QB is the base flow component (low frequency signal), i is the time step index, and α is the separation (or filter) parameter (0 < α < 1) that enables the separation to be altered. The base flow component (QB) for each time step is constrained to be never greater than the total flow (Q) (or less than 0). Traditionally, most of the digital filtering separation techniques dealt with daily time series. However, within the framework of this study, a special case is of interest, namely developing a technique which focuses on the monthly flow data [53]. The latter paper proposes to modify the previously used daily separation technique to make it work with monthly flow data. The suggested initial step required is to perform calculations based on daily time series (that is, defining a filter parameter and separating flow components). These results are then adjusted to apply the monthly flow records. Adjustments are required because of the obvious differences in flow variability displayed in daily and monthly flow records. The filter parameter value must be optimized so that the flow components estimated by the monthly separation are consistent with the daily "prototype". By this way, the monthly flow separation technique is verified against the results of the daily one. The proposed method [53] is supported by reference to the widespread use of the daily flow separation in South Africa, as well as to a possible application of regional parameters and the hydrological analogy. The paper presents examples of the successful use of the monthly flow separation and discusses the application of monthly base flow data in the context of ecological reserve issues.
At this stage of the present study, the method of the flow separation (one-parameter recursive digital filter by Equations (2) and (3)) has been selected, the possibility of its application at a monthly time scale is confirmed [53], and problem of the daily "prototype" only is unresolved. To solve the latter problem, we used the results of the daily flow separation which was carried out within the framework of the Hydrological Model for Karst Environment (HYMKE) [16]. The conceptual HYMKE model was applied simultaneously to the major tributaries (Dan, Hermon and Snir) of the Upper Jordan River. Within the model, long-term daily series  of base flow and surface flow were calculated for every tributary using the two-parameter recursive digital filter [54]. Unfortunately, these series are not available to users, and the only data that can be extracted from the publication are examples of hydrographs and the mean annual base and surface flow for each tributary.
Nevertheless, using this limited information it was possible to solve the problem. For the Snir (Hasbani) River the parameter α (in Equation (2)) was calibrated with a process of "trial and error" to obtain a base flow index (BFI) as close as possible to the value derived from [16]. BFI is the ratio of cumulative base flow to cumulative total flow over the time period of analysis. It is dimensionless and varies between 0 and 1. BFI is a hydrogeological parameter which is useful in modeling un-gaged basins [55] and for presentation of the effect of geology on basin low flows [56]. BFI can be used both as a standard for comparison of different basins and for calibrating parameters in analytical methods of the flow separation [57].
As seen in Table 5, results of the parameter calibration are quite satisfactory. Comparison between the calculated flow components (with α equal to 0.890) and the HYMKE results  shows that the surface flow is unchanged, the base flow and BFI decrease by 5% and 2%, respectively. Extension of the calculation period to 1940-2020 does not change the surface flow and decreases the base flow and BFI by 3% and 2%. Respectively, monthly (in this study) and daily [16] hydrographs look similar ( Figure 13).    (2) and (3)) at the monthly time-scale (in the present study) (a). Model HYMKE [16]: draft of separated hydrograph (b) and hydrographs (1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000) presented results of an application of the daily two-parameter recursive digital filter to the Hasbani (Snir) River (c).
Before proceeding to analysis of the flow separation results for the Hasbani River, it is worth making reference to several important remarks: -It is generally acknowledged that majority of the proposed numerical techniques for the flow separation (except for tracer-based methods) are not closely related to the underlying physical processes but offer a way forward for hydrologic practice if they can provide objective and repeatable results [52,58]. -Most hydrograph separations (apart from tracer-based separations) lack a physical basis. Therefore, choosing one method or another, introduces an undesirable element of uncertainty and randomness into the analysis [59,60]. - The time intensive and expensive separation of event hydrographs by natural tracers such as environmental isotopes and geochemical constituents is probably the only method to determine realistically the runoff components [61][62][63]. -Advantages of the filtering methods include use of the only total flow record, easy calculations, objectivity, and replicability of results [54,64].
The above presented comments explain the choice of the one-parameter recursive digital filter as the flow separation method in the present study. It is completely unrealistic to assume that data required for tracer-based flow separations are available in the scarce-data transboundary basin (may be, except in some rare cases).  (2) and (3)) at the monthly time-scale (in the present study) (a). Model HYMKE [16]: draft of separated hydrograph (b) and hydrographs (1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000) presented results of an application of the daily two-parameter recursive digital filter to the Hasbani (Snir) River (c).
Before proceeding to analysis of the flow separation results for the Hasbani River, it is worth making reference to several important remarks: -It is generally acknowledged that majority of the proposed numerical techniques for the flow separation (except for tracer-based methods) are not closely related to the underlying physical processes but offer a way forward for hydrologic practice if they can provide objective and repeatable results [52,58]. -Most hydrograph separations (apart from tracer-based separations) lack a physical basis. Therefore, choosing one method or another, introduces an undesirable element of uncertainty and randomness into the analysis [59,60]. - The time intensive and expensive separation of event hydrographs by natural tracers such as environmental isotopes and geochemical constituents is probably the only method to determine realistically the runoff components [61][62][63]. -Advantages of the filtering methods include use of the only total flow record, easy calculations, objectivity, and replicability of results [54,64].
The above presented comments explain the choice of the one-parameter recursive digital filter as the flow separation method in the present study. It is completely unrealistic to assume that data required for tracer-based flow separations are available in the scarcedata transboundary basin (may be, except in some rare cases). Authors of the model HYMKE [16] (whose results are used here to calibrate the parameter α in Equation (2)) note the physical constraint of the recursive digital filter. This method cannot distinguish between surface flow and quick responses of underground flows, and the term "surface flow" is actually attributed to both types of flow. This remark can also apply to the rather formal division of hydrograph into surface and base flow depending on the separation parameters. Therefore, it is important to verify the base flow separation technique by spring flow data as it is done in the model HYMKE. An argument to support the separating flow module in HYMKE is the satisfactory final testing and verification of this conceptual rainfall-streamflow model. To a certain extent, the latter argument supports also the corresponding results in the present study.
Results of the hydrograph separation into surface flow and base (groundwater) flow for the period 1940-2020 are shown in Figure 14. Attention is drawn to the decrease by more than half of both components over the indicated period, with the base flow falling more than the surface flow. It corresponds to the significant decrease (55%) of the Hasbani River annual flow during 1940-2017 ( Figure 2). Both detected linear trends are statistically significant (Appendix B, Table A1). Authors of the model HYMKE [16] (whose results are used here to calibrate the parameter α in Equation (2)) note the physical constraint of the recursive digital filter. This method cannot distinguish between surface flow and quick responses of underground flows, and the term "surface flow" is actually attributed to both types of flow. This remark can also apply to the rather formal division of hydrograph into surface and base flow depending on the separation parameters. Therefore, it is important to verify the base flow separation technique by spring flow data as it is done in the model HYMKE. An argument to support the separating flow module in HYMKE is the satisfactory final testing and verification of this conceptual rainfall-streamflow model. To a certain extent, the latter argument supports also the corresponding results in the present study.
Results of the hydrograph separation into surface flow and base (groundwater) flow for the period 1940-2020 are shown in Figure 14. Attention is drawn to the decrease by more than half of both components over the indicated period, with the base flow falling more than the surface flow. It corresponds to the significant decrease (55%) of the Hasbani River annual flow during 1940-2017 ( Figure 2). Both detected linear trends are statistically significant (Appendix B, Table A1). Thus, it turned out to be possible to separate the Hasbani flow on the surface and base flow components using data on the monthly flow. The obtained annual series (Figure 14) serves as a base for assessing withdrawals from the surface and groundwater flow of the Hasbani River. Thus, it turned out to be possible to separate the Hasbani flow on the surface and base flow components using data on the monthly flow. The obtained annual series (Figure 14) serves as a base for assessing withdrawals from the surface and groundwater flow of the Hasbani River. The established Natural Flow Curve (1940-1971 without 1952-1962) is the upper envelope of the entire field of points related to the following 1972-2017 period.

Surface and Base Flow Withdrawals
-Within the latter period, it is possible to identify two specific sub-periods. For the first sub-period of 26 years (1972-1997) there is the high accuracy dependence of surface and base flow on precipitation (correlation coefficient 0.969 and 0.897, respectively). In the following sub-period of 20 years (1998-2017) curves are placed lower with correlation coefficient (0.913 and 0.758). Values of correlation coefficient for both sub-periods are significant (Appendix B, Table A2). Consistent flow decline is observed during successive sub-periods (1972-1997 and 1998-2017). This means that for following sub-period the same precipitation in the basin formed less surface and base flow because of increasing water withdrawals (by above preliminary interpretation of this term).  Table 6 represents the summary of withdrawals from the surface and base flow as compared with the separated flow of the Hasbani River during successive 12 (13) year sub-periods (1972-1983, 1984-1995, 1996-2007 and 2008-2020). These data indicate the following: By the AWW method, the natural PF curves (Figure 15a,b) serve as a basis for assessing the water withdrawals from surface and base flow . The annual volume of withdrawals is estimated by Equation (1) Table A1). Both figures demonstrate that the increasing component of the base flow constitutes an increasing proportion of the total withdrawal during the study period. -During the first to fourth period, the measured base flow dropped slightly less than by half and BF withdrawals doubled. As a result, the portion of BF withdrawals within the natural base flow increased from 17% to 45% (twice as many as SF withdrawals).
These data confirm the trends in Figure 16a, indicate a significant increase in base flow withdrawals, and corroborate contribution of both components to decrease in the total flow of the Hasbani River (Table 2).   Table 6 represents the summary of withdrawals from the surface and base flow as compared with the separated flow of the Hasbani River during successive 12 (13) year sub-periods (1972-1983, 1984-1995, 1996-2007 and 2008-2020). These data indicate the following: - During the first to fourth period, measured surface flow decreased by almost a third and SF withdrawals increased by about one-third. As a result, the portion of SF withdrawals within the natural surface flow increased from 25% to 39%. -During the first to fourth period, the measured base flow dropped slightly less than by half and BF withdrawals doubled. As a result, the portion of BF withdrawals within the natural base flow increased from 17% to 45% (twice as many as SF withdrawals). Note to Table 6: water withdrawals (%) are related to the natural surface or base flow (sum of corresponding measured flow and withdrawals).
These data confirm the trends in Figure 16a, indicate a significant increase in base flow withdrawals, and corroborate contribution of both components to decrease in the total flow of the Hasbani River (Table 2).
Finally, it should be noted that this section, which at first seemed optional or not necessary within the main research topic, provided ultimately significant results. This concerns new information on the structure of withdrawals (with separation to surface and base flow components) and the identification of the increasing contribution of the base flow withdrawals in use of the Hasbani flow since the beginning of the 2000s. The latter indicates the development of a network of wells and an increase in their total production rate. In addition, the estimated surface flow withdrawals are useful for attribution of drivers in decreasing flow trend of the Hasbani River (Section 6.2).

Comparison of Obtained Assessments with Published Values
Published sources offer the following estimates of water use in the Hasbani Basin in Lebanon: 7 MCM/y [1], 8 MCM/y [25,26,28], and 9-11 MCM/y [1,27,65]. These low estimates (generally without mention of specific years) were commented by authors as rough and incomplete, with the evident lack of reliable information, particularly for unmonitored groundwater abstraction. According to the reports of the Israel Hydrological Service [17,66], water use in the Hasbani Basin (Lebanon) was assumed to be higher (to 20 MCM/y), but no arguments or estimates were provided.
Of particular interest is the fact that the total water withdrawals in the Hasbani River Basin (Lebanon) estimated by the AWW method (see Section 5.1) exceed by far the previously published estimates of water. Such a significant discrepancy between the results requires a special discussion.
In the first place, attention should be given to the low accuracy of the published estimates, which are most often accompanied by the following author's comments (presented here not literally but close to the text): - The south of Lebanon in particular lack data on water resources based on accurate, sustained, and reliable observation. There are gaps in data, notably on water use and water quality. Some of the data may be fragmented, outdated, or very rough estimates. Historic water use in the south Lebanon may be more than is normally held, though no estimate based on archival record has been attempted [27]. -Estimates of water use from the Upper Jordan River in Lebanon are greatly compromised by poor data and not solidly supported by metered records [26]. -In the Hasbani Basin, southern Lebanon, there is a lack of authoritative knowledge of regional water resources. At the national level the monitoring of hydrological data is fragmented and unreliable, with no central agency analyzing and disseminating relevant information [28,67].
In the Hasbani Basin, 70% of the regional population works in an agricultural sector dominated by small, family-run landholdings [28]. There are a great number of private wells and most of them are unlicensed and unmonitored. It is noted that "as with other parts of Lebanon, information on abstraction rates from illegal wells is wanting, but any estimate of their actual abstractions cannot be considered as accurate" [26]. Information pertaining to private wells is difficult to gather because the owners are reluctant to provide any information regarding the depths and discharges of the wells [25].
Even the number of wells in the Hasbani Basin is uncertain and disputable. According to [26], the Hasbani catchment is served by 9 public wells (generally for domestic supply to the villages), 27-40 private wells (for irrigation and domestic purposes) are known, of which only 2 wells are monitored. According to [25], 3 main public wells are used for water supply in the Hasbani basin, 12 wells were drilled in the study area (4 wells are under the supervision of the Water Authority and 8 others are unmonitored), and discussions with local farmers revealed the existence of 27 additional private wells, most of which are not registered.
The FAO declared that it is practically impossible to determine the exact figure for groundwater abstraction in the Hasbani Basin [8]. By [26], the rough estimate (11 MCM/y) of the Hasbani water used in Lebanon for local irrigation and drinking water includes about 4.4 MCM/y of groundwater abstracted from licensed and unlicensed wells. Unlike it, in the present paper the total withdrawals from the Hasbani River are estimated (by the AWW method) at 30 MCM/y (197230 MCM/y ( -1983) and 48 MCM/y (2008-2020) including respectively 13 MCM/y and 26 MCM/y as the base flow component (Tables 2 and 6). The latter increase is undoubtedly caused by increased number of private (mainly unregistered) wells. By indirect indices (areas of irrigated land, population, and water consumption rates), the total (irrigation and domestic) water use in the Hasbani Basin (before 2011) is estimated as no less than 29 MCM/y (Section 5.3). This figure supports value of 34 MCM/y which in the present study is the calculated water withdrawals for 1972-2010.
It should be noted that water withdrawals estimated in this study are of the same order as Lebanon's share of Jordan River water use (of 35 MCM/y) under the Johnston Plan Allocation. The Johnston Plan (referring only to surface flow) was proposed to the riparian countries in 1955 but was never ratified. Nevertheless, it became a benchmark for water management and negotiations in the Jordan River basin [1,65,68].
Substantial issues for further consideration are: (i) What water withdrawals means? (ii) how efficiently the abstracted water is used? and (iii) which part of the water use has been accounted for in published estimates (as stated above, incomplete, rough and etc.,)?
The answer to the first question is discussed in the next section. Regarding the second question, a large share of water in the distribution systems is lost through leakages (partly recharging the local aquifer). FAO reported 35-50% seepage from the public water supply networks in Lebanon [8]. One would not expect a better situation in the private sector. The answer to the third question remains uncertain.

Flow Trend Attribution: What Are Water Withdrawals?
In the present study, the statistically significant decreasing trend in precipitation and flow of the Hasbani River, as well as the relationship between them is detected and confirmed. This means that a decrease in precipitation (climate signal) can be considered as an obvious driving force behind a decrease in flow. But this driver is not the only one, as the drop in flow exceeds the decreasing trend in precipitation.
Further analysis revealed the dynamic nature of the PF relationship, which indicated a consistent flow decline in subsequent sub-periods when with the same precipitation in the basin, less and less runoff was measured at the downstream station. Most likely, this effect is associated with some water withdrawals primarily caused by human activities in the basin (development of irrigation, economic needs, etc.,). However the influence of other factors is not excluded. So far, we used the term "withdrawals" without detailed interpretation, assuming (by the AWW method definition) that withdrawals are the total amount of water removed from the river (or did not reach it) and not included in the measured flow. Estimates of indirect indices (irrigated land area, population, and water consumption rates) confirmed that water use is a significant driving force behind a decrease in the Hasbani flow and almost completely explains the amount of water withdrawals. Attributing flow changes to other specific drivers (apart from decrease in precipitation and developing water use) is discussed further in this section.
First of all, this refers to the global air warming hypothesis, according to which increase in air temperature leads to increase in evapotranspiration and decrease in flow. For the Jordan River headwaters, potential evapotranspiration (PET) was estimated from daily minimum and maximum temperature data in the Har Qnaan weather station along with extraterrestrial solar radiation [19,69]. While rising temperature in the basin is statistically significant and the calculated annual PET increases from 1130 mm in 1984-1988 to 1160 mm in the 2007-2016 period, this change (about 2.6%) is too small to explain meaningfully the observed streamflow decreases. This is confirmed by the following simplified calculation. If the change in PET (2.6%) is roughly (as an upper and unrealistically high limit) attributed to the actual annual evapotranspiration (ET) in the Hasbani Basin (from 226 mm [16] to 238 mm [70]), then the increase in ET is estimated as less than 5-6 mm (or 3-4 MCM) for 1984-2016 period. This means that global air warming can be considered as only a secondary driver for reducing the Hasbani flow (in compare with Figure 2 and Table 2).
In the present study, the latter conclusion is confirmed by assessing the possible impact of the global air warming on the Hasbani surface flow decrease. This effect can be estimated as not higher than 1 MCM during 1972-2007, assuming that it only caused a slight increase in surface flow withdrawals during this period (Table 6).
Further, the influence of previous dry years on the Hasbani flow decrease was considered as a possible driver. Such effect of "hydrological memory" (or "drought hypothesis") was found for incoming water to the Lake Kinneret [10,71] and for the Dan spring flow [9,16] (in consideration of one previous year and three previous years, respectively). In the present case, such effect could be real due to the facts that more than half of the Hasbani basin area is represented by karst exposures and the river is fed mainly by springs. The impact of the previous dry year on flow and water withdrawals in the current year was analyzed within the AWW method. As a result, a slight tendency was revealed, which turned out to be statistically insignificant against the influence of water use.
Summing up the above discussion, the following can be argued: -Decrease in precipitation (climate signal) and human activities in the basin (anthropogenic factor) are the main driving forces behind a decrease in the Hasbani flow. -Global air warming can be considered only as a secondary driver for reducing the Hasbani flow. -Water withdrawals from the Hasbani River estimated by the AWW method are generally results of the human water use with a negligible contribution due to global air warming.
Thus, the estimates of water withdrawals and water use are quite comparable. Moreover, it is quite possible to assume that for water resources management or for assessing the water balance in a transboundary basin, it can be more important to assess exactly withdrawals (as a water deficit) than the water consumption.
It is important to reiterate that the AWW method enables to differentiate between climatic and anthropogenic impacts on the river flow decrease. This follows from the definition of the AWW method, since the dynamics of the precipitation-flow relationship is analyzed, and not only the trend to a decrease in the flow. Therefore the present study enabled to estimate the relative contribution of two main drivers to the observed decrease of flow. As it was confirmed, a decrease in precipitation is an obvious driving force behind a decrease in the Hasbani flow. Correlation coefficient (R) between precipitation and measured flow (1972-2020) is 0.918, i.e., this driver describes 84% of variance (or dispersion, D) of flow. For natural flow (sum of measured flow and calculated withdrawals) these values are estimated as R = 0.974 and D = 95%, i.e., 11% more taking into consideration the water withdrawals as a balance component of flow.
On the other hand, regarding results of trend analysis, the statistically significant decrease in precipitation and flow of the Hasbani River (1972-2020) was detected as following: for precipitation-15%; for measured flow-45%; and for natural flow (sum of measured flow and calculated withdrawals-19%. As can be seen, in the absence of water withdrawals, the drop in natural flow would be slowed by 26% compared to the trend in measured flow and would approach the trend in precipitation. The latter (i) provides evidence that in this case water use is a more responsible driver in reducing flow than precipitation, and (ii) serves (together with Figure 7) as verification of the AWW method applied to the Hasbani River.

Conclusions and Recommendations
In the present study, the AWW method was proposed as a nonstandard solution to the problem of assessing water withdrawals in the scarce-data transboundary basin. The method operates with the open-source available data on precipitation and river flow and thereby overcomes the usual restriction due to lack of data on shared water use in the Middle East. Analysis of dynamic precipitation-flow relationships enabled to separate the effect of water withdrawals from the total decline of river flow under the decreasing precipitation.
According to the published sources, the present study is the first which provided the complete and rather detailed estimates of water withdrawals from total, surface, and groundwater flow of the Hasbani River in Lebanon (1972-2020). The obtained estimates were confirmed by indirect indices such as area of the irrigated agricultural land and population in the Hasbani Basin. The strong argument in support of the results was validation of the AWW method based on the independent data on water use in the Israeli part of Upper Catchments of the Jordan River (1975-2014, Israel Water Authority).
Careful study revealed and confirmed that (i) decrease in precipitation (climate signal) and developed human activities in the basin (anthropogenic factor) are the main driving forces behind a decrease in the Hasbani flow, and (ii) water withdrawals estimated by the AWW method are generally the results of water use with a negligible contribution of increase in evapotranspiration due to global air warming.
The results of the study are useful for water balance estimations, as well as for management of water resources in the Jordan River headwaters basin and in the entire Lake Kinneret Basin.
The current research should not be regarded only as a regional study but it offers a general approach to solving similar problems. The AWW method can be applied to other transboundary basins under conditions of uncertain (or limited) data on water use, but subject both to the available precipitation and flow data and to confirmation of main hypotheses formulated in Section 3.2.
The AWW method enables historical and real-time monitoring of water withdrawals in any riparian state and in the entire transboundary basin. Undoubtedly, reliable data on water demand and water use should be the necessary basis for settlement of problematic transboundary water relations, especially within the current critical situation in the Middle East due to expanding geopolitical, socio-economic, and humanitarian problems.
In order to maintain the transparency and reproducibility of the research, "Supplementary File.xlsx" is available online at www.mdpi.com/xxx/s1. Acknowledgments: Thanks are due to the Israel Meteorological Service and the Israel Hydrological Service for providing long data series of monthly and yearly precipitation and flow. The authors gratefully acknowledge the Israel Water Authority and personally Gilad Sapir and Omer Sagi for data on the water use in Upper Catchments of the Jordan River. We are all grateful to Rita Chudinov for the help in analysis of precipitation and flow data. The authors acknowledge the anonymous Reviewers whose comments contributed greatly to improve this paper.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
Some typical examples of studies aimed at assessing impact of water use on water resources by the surface-subsurface flow model:

•
Potential effects of changes in water use (in the Carson Valley, Douglas County, Nevada, and the Alpine County, California) were assessed with the MODFLOW model presenting river-aquifer interaction [72]. The model was calibrated during 1971-2005 period and applied for predicting 55-years simulation (2006-2060) to evaluate the long-term effects of different water-use scenarios on groundwater levels and streamflow. Increasing withdrawals resulted in water-table decline, reduced volume of groundwater storage, and decrease in the downstream flow. The model was proposed as a tool that can assist water managers and planners in making informed decisions.

•
The groundwater policy guidance was proposed for the Indus Basin in Pakistan [73].
Using the same MODFLOW model, simulations were performed to evaluate the groundwater dynamics in the future under different scenarios of groundwater pumping, canal infrastructure improvements, and precipitation changes. The results were highlighted with the mapping of changes in water table, pumping cost, and waterlogged area. The results showed that changes in both groundwater abstraction and seepage from the canal system significantly impacted groundwater heads, whereas the effect of changing precipitation was negligible. Under status quo conditions, the average province-wide pumping cost was projected to increase by 270% in 23 years.

•
The integrated surface-subsurface WASH123D model was applied to simulate groundwater levels under the pumping-free scenario over the Pingtung Alluvial Plain in Southwest Taiwan [74]. The significant discrepancies found between the simulated and observed GLs were indicative of excessive (mostly illegal) pumping which induced serious regional land subsidence and seawater intrusion.

•
The water-withdrawal assessment process and online screening tool (WWAT) were developed in Michigan, USA [75]. The screening tool is designed to evaluate intermittent pumping, to account for the dynamics of stream-aquifer interaction, and to apportion streamflow depletion among neighboring streams. The tool is to be used for an initial screening of a proposed new or increased high-capacity withdrawal in order to identify the withdrawals that may cause adverse impacts on resources. Users enter the location, timing, quantity, and if relevant, the screen depth of their proposed groundwater or surface water withdrawals. Using the hydrologic foundation and groundwater model, WWAT associates the depletion with its risk level. Once reviewed, the withdrawal will be registered, registered with modifications, or rejected. By associating different policy actions with different levels of risk, the system compensates for uncertainties in the models used to quantify the Adverse Resource Impact.   [78].