Semi-arid Catchment Research Articles

A large-scale study of hyporheic nitrate dynamics in a semi-arid catchment, the Tafna River, in Northwest Algeria

Very few studies have evaluated the role of the hyporheic zone (HZ) in aquatic ecosystem functioning or the factors driving hyporheic exchange flows on a large scale, especially in a semi-arid environment such as the Tafna watershed in Algeria. To understand this role through time and space, hydrogeochemical parameters, particularly nitrate concentrations (NO3–N), were measured monthly between February 2013 and April 2014 in surface water (SW), interstitial water (IW) and groundwater (GW) along a 170-km stretch from the river’s source to the outlet into the main stream of the Tafna River and its Isser and Chouly tributaries. The significant longitudinal evolution of NO3–N in the HZ was related to stream order and distance from the source. Moreover, the study indicated a significant difference between nitrate concentrations in the riffles (R) and pools (P) of the Tafna wadi. Principal component analysis (PCA) revealed the considerable impact of agriculture on nitrate concentrations. This study indicated on a large scale that nitrate-rich HZ contributed to increasing surface nitrate concentrations in upwelling sites and could be an important nitrate source for downstream SW, particularly during low-water (LW) periods when the Tafna can run dry on the surface. Thus, these results underline the importance of hyporheic zones functioning to the water quality of the watershed (process of enrichment and retention of nitrogen).

Combining 3D data and traditional soil erosion assessment techniques to study the effect of a vegetation cover gradient on hillslope runoff and soil erosion in a semi-arid catchment

In this paper, we examine the effect of vegetation on soil erosion, runoff generation and sediment transport on saline rangeland hillslopes. Rainfall simulations were conducted at a fixed 114 mm/h intensity on 6 m × 2 m erosion plots with varying degrees of vegetation. Plots were grouped into three categories (L, M and H) based on their canopy cover (L: <5%, M: 5–19%, H: >19%) and selected to limit variations in slope across canopy cover groups. Runoff and sediment samples were combined with three dimensional (3D) reconstruction data used to monitor soil surface microtopographic changes. Runoff initiation was significantly delayed on the L plots but cumulative runoff after 20 min of rainfall simulation indicated a positive effect of vegetation on infiltration processes. Cumulative sediment after 20 min of rainfall was similar across vegetation cover categories. The 3D data suggest that vegetation reduced net sediment delivery from the plots by primarily increasing opportunities for deposition while marginally affecting gross soil erosion. Plots with H vegetation cover experienced lower 3D-estimated erosion volumes but average erosion depth on L plots was lower than that on plots with M and H vegetation covers. Lower runoff volumes on M and H plots may have been compensated by greater runoff erosivity on these plots as runoff was concentrated in a narrower inter-patch space compared to L plots. This study highlights the need for an increased integration between traditional runoff measurement techniques and 3D reconstruction methods.

Characteristics of wet dissolved carbon deposition in a semi-arid catchment at the Loess Plateau, China

Abstract. Wet dissolved carbon deposition is a critical node of the global carbon cycle, but little is known about dissolved organic and inorganic carbon (DOC and DIC) concentrations and fluxes in the semi-arid areas of the Loess Plateau Region (LPR). In this study, we measured variations in DOC and DIC concentrations in rainfalls at Yangjuangou Ecological Restoration and Soil and Water Observatory. Rainwater samples were collected in 16 rainfall events from July to September and the event-based, monthly concentrations and fluxes of DOC and DIC were quantified. The results showed that the event-based concentrations and fluxes of DOC and DIC were highly variable, ranging from 0.56 to 28.71 mg C L−1 and from 3.47 to 17.49 mg C L−1, respectively. The corresponding event-based fluxes ranged from 0.21–258.36 mg C m−2 and from 4.12 to 42.32 mg C m−2. The monthly concentrations of DOC and DIC were 24.62 and 4.30 (July), 3.58 and 10.52 (August), and 1.01 and 5.89 (September) mg C L−1, respectively. Thus, the monthly deposition fluxes of DOC and DIC were 541.64/94.60, 131.03/385.03, and 44.44/259.16 mg C m−2 for July, August, and September, respectively. In addition, the concentrations of DOC and DIC for the concentrated rainfall season (July–September) in the studied catchment were 7.06 and 7.00 mg C L−1, respectively. The estimated annual wet dissolved carbon depositions were 1.91 and 1.89 g C m−2 yr−1 for DOC and DIC, respectively. The results of this study suggest the variation in concentrations and fluxes of DOC and DIC and explore that these variation may be related to the dissolved carbon source and the rainfall characteristics during the concentrated rainfall season in the semi-arid catchment of the LPR. Furthermore, these results also suggest that dissolved carbon may be an important external input of carbon into terrestrial ecosystems.

Modelling of channel transmission loss processes in semi-arid catchments of southern Africa using the Pitman Model

Abstract. Water availability is one of the major societal issues facing the world. The ability to understand and quantify the impact of key hydrological processes, on the availability of water resources, is therefore integral to ensuring equitable and sustainable resource management. Channel transmission losses are an “under-researched” hydrological process that affects resource availability in many semi-arid regions such as the Limpopo River Basin in southern Africa, where the loss processes amount to approximately 30 % of the water balance. To improve the understanding of these loss processes and test the capability of modelling routines, three approaches using the Pitman model are applied to selected alluvial aquifer environments. The three approaches are an explicit transmission loss function, the use of a wetland function to represent channel-floodplain storage exchanges and the use of a “dummy” reservoir to represent floodplain storage and evapotranspiration losses. Results indicate that all three approaches are able to simulate channel transmission losses with differing impacts on the regional flows. A determination of which method best represents the channel transmission losses process requires further testing in a study area that has reliable observed historical records.

Spatial and temporal dynamics of deep percolation, lag time and recharge in an irrigated semi-arid region

The time required for deep percolating water to reach the water table can be considerable in areas with a thick vadose zone. Sustainable groundwater management, therefore, has to consider the spatial and temporal dynamics of groundwater recharge. The key parameters that control the lag time have been widely examined in soil physics using small-scale lysimeters and modeling studies. However, only a small number of studies have analyzed how deep-percolation rates affect groundwater recharge dynamics over large spatial scales. This study examined how the parameters influencing lag time affect groundwater recharge in a semi-arid catchment under irrigation (in northeastern Iran) using a numerical modeling approach. Flow simulations were performed by the MODFLOW-NWT code with the Vadose-Zone Flow (UZF) Package. Calibration of the groundwater model was based on data from 48 observation wells. Flow simulations showed that lag times vary from 1 to more than 100 months. A sensitivity analysis demonstrated that during drought conditions, the lag time was highly sensitive to the rate of deep percolation. The study illustrated two critical points: (1) the importance of providing estimates of the lag time as a basis for sustainable groundwater management, and (2) lag time not only depends on factors such as soil hydraulic conductivity or vadose zone depth but also depends on the deep-percolation rates and the antecedent soil-moisture condition. Therefore, estimates of the lag time have to be associated with specific percolation rates, in addition to depth to groundwater and soil properties.

Water quality modelling of an impacted semi-arid catchment using flow data from the WEAP model

Abstract. The continuous decline in water quality in many regions is forcing a shift from quantity-based water resources management to a greater emphasis on water quality management. Water quality models can act as invaluable tools as they facilitate a conceptual understanding of processes affecting water quality and can be used to investigate the water quality consequences of management scenarios. In South Africa, the Water Quality Systems Assessment Model (WQSAM) was developed as a management-focussed water quality model that is relatively simple to be able to utilise the small amount of available observed data. Importantly, WQSAM explicitly links to systems (yield) models routinely used in water resources management in South Africa by using their flow output to drive water quality simulations. Although WQSAM has been shown to be able to represent the variability of water quality in South African rivers, its focus on management from a South African perspective limits its use to within southern African regions for which specific systems model setups exist. Facilitating the use of WQSAM within catchments outside of southern Africa and within catchments for which these systems model setups to not exist would require WQSAM to be able to link to a simple-to-use and internationally-applied systems model. One such systems model is the Water Evaluation and Planning (WEAP) model, which incorporates a rainfall-runoff component (natural hydrology), and reservoir storage, return flows and abstractions (systems modelling), but within which water quality modelling facilities are rudimentary. The aims of the current study were therefore to: (1) adapt the WQSAM model to be able to use as input the flow outputs of the WEAP model and; (2) provide an initial assessment of how successful this linkage was by application of the WEAP and WQSAM models to the Buffalo River for historical conditions; a small, semi-arid and impacted catchment in the Eastern Cape of South Africa. The simulations of the two models were compared to the available observed data, with the initial focus within WQSAM on a simulation of instream total dissolved solids (TDS) and nutrient concentrations. The WEAP model was able to adequately simulate flow in the Buffalo River catchment, with consideration of human inputs and outputs. WQSAM was adapted to successfully take as input the flow output of the WEAP model, and the simulations of nutrients by WQSAM provided a good representation of the variability of observed nutrient concentrations in the catchment. This study showed that the WQSAM model is able to accept flow inputs from the WEAP model, and that this approach is able to provide satisfactory estimates of both flow and water quality for a small, semi-arid and impacted catchment. It is hoped that this research will encourage the application of WQSAM to an increased number of catchments within southern Africa and beyond.

Regionalising MUSLE factors for application to a data-scarce catchment

Abstract. The estimation of soil loss and sediment transport is important for effective management of catchments. A model for semi-arid catchments in southern Africa has been developed; however, simplification of the model parameters and further testing are required. Soil loss is calculated through the Modified Universal Soil Loss Equation (MUSLE). The aims of the current study were to: (1) regionalise the MUSLE erodibility factors and; (2) perform a sensitivity analysis and validate the soil loss outputs against independently-estimated measures. The regionalisation was developed using Geographic Information Systems (GIS) coverages. The model was applied to a high erosion semi-arid region in the Eastern Cape, South Africa. Sensitivity analysis indicated model outputs to be more sensitive to the vegetation cover factor. The simulated soil loss estimates of 40 t ha−1 yr−1 were within the range of estimates by previous studies. The outcome of the present research is a framework for parameter estimation for the MUSLE through regionalisation. This is part of the ongoing development of a model which can estimate soil loss and sediment delivery at broad spatial and temporal scales.

Semi-arid catchments under change: Adapted hydrological models to simulate the influence of climate change and human activities on rainfall-runoff processes in southern Africa

Semi-arid catchments under change: Adapted hydrological models to simulate the influence of climate change and human activities on rainfall-runoff processes in southern Africa

Modeling of flood generation in semi-arid catchment using a spatially distributed model: case of study Wadi Mekerra catchment (Northwest Algeria)

This paper aims to apply the MERCEDES distributed hydrological model with an event-based mode to simulate flood generation in the Wadi Mekerra basin situated northwest of Algeria. This catchment is characterized by a semi-arid climate, convective thunderstorms, and ephemeral flow. The input data are mainly the daily rainfall-runoff and numerical maps such as slope, flow direction, and land use. The comparison between calculated and measured runoff revealed that the MERCEDES performance to simulate flood generation in Wadi Mekerra basin is encouraging and satisfying which is justified by the Nash–Sutcliffe and observations standard deviation ratio criteria. However, the proposed model tends to underestimate peak runoff due to convective rainfall that may be much localized in space. Furthermore, the sensitivity analysis shows that the hydrological response in the Wadi Mekerra catchment depends strongly on the potential maximum retention parameter which is related to the land use type and varies significantly between seasons according to the vegetation cover dynamics, rainfall intensity, and drought.

Estimation of groundwater recharge via percolation outputs from a rainfall/runoff model for the Verlorenvlei estuarine system, west coast, South Africa

Wetlands are conservation priorities worldwide, due to their high biodiversity and productivity, but are under threat from agricultural and climate change stresses. To improve the water management practices and resource allocation in these complex systems, a modelling approach has been developed to estimate potential recharge for data poor catchments using rainfall data and basic assumptions regarding soil and aquifer properties. The Verlorenvlei estuarine lake (RAMSAR #525) on the west coast of South Africa is a data poor catchment where rainfall records have been supplemented with farmer’s rainfall records. The catchment has multiple competing users. To determine the ecological reserve for the wetlands, the spatial and temporal distribution of recharge had to be well constrained using the J2000 rainfall/runoff model. The majority of rainfall occurs in the mountains (±650 mm/yr) and considerably less in the valley (±280 mm/yr). Percolation was modelled as ∼3.6% of rainfall in the driest parts of the catchment, ∼10% of rainfall in the moderately wet parts of the catchment and ∼8.4% but up to 28.9% of rainfall in the wettest parts of the catchment. The model results are representative of rainfall and water level measurements in the catchment, and compare well with water table fluctuation technique, although estimates are dissimilar to previous estimates within the catchment. This is most likely due to the daily timestep nature of the model, in comparison to other yearly average methods. These results go some way in understanding the fact that although most semi-arid catchments have very low yearly recharge estimates, they are still capable of sustaining high biodiversity levels. This demonstrates the importance of incorporating shorter term recharge event modeling for improving recharge estimates.

Modelling of runoff response in semi-arid catchment of Paraíba State using SWAT

Modelling of runoff response in semi-arid catchment of Paraíba State using SWAT

Cross-Scale Precipitation Variability in a Semiarid Catchment Area on the Western Slopes of the Central Andes

AbstractSpatiotemporal precipitation patterns were investigated on the western slopes of the central Andes Mountains by applying EOF and cluster analysis as well as the Weather Research and Forecasting (WRF) Model. In the semiarid catchment area in the highlands of Lima, Peru, the precipitation is assumed to be a cross-scale interplay of large-scale dynamics, varying sea surface temperatures (SSTs), and breeze-dominated slope flows. The EOF analysis was used to encompass and elucidate the upper-level circulation patterns dominating the transport of moisture. To delineate local precipitation regimes, a partitioning cluster analysis was carried out, which additionally should illustrate local effects such as the altitudinal gradient of the Andes. The results demonstrated that especially during the transition to the dry season, synoptic-scale circulation aloft controls the precipitation (correlation coefficients between 0.6 and 0.9), whereas in the remaining seasons the slope breezes due to the altitudinal gradient mainly determine the precipitation behavior. Further analysis with regard to the spatiotemporal precipitation variability revealed an inversion of the precipitation distribution along the elevational gradient within the study area, mainly during February (29%) and March (35%), that showed correlations with coastal SST patterns ranging between 0.56 and 0.67. WRF simulations of the underlying mechanisms disclosed that the large-scale circulation influences the thermally induced upslope flows while the strength of southeastern low-level winds related to the coastal SSTs caused a blocking of easterlies in the middle troposphere through a reduced anticyclonic effect. This interplay enables the generation of precipitation in the usually drier environment at lower elevations, which leads to a decrease in rainfall with increasing elevation.

INTEGRATION OF SATELLITE RAINFALL DATA AND CURVE NUMBER METHOD FOR RUNOFF ESTIMATION UNDER SEMI-ARID WADI SYSTEM

Abstract. The arid and semi-arid catchments in dry lands in general require a special effective management as the scarcity of resources and information which is needed to leverage studies and investigations is the common characteristic. Hydrology is one of the most important elements in the management of resources. Deep understanding of hydrological responses is the key towards better planning and land management. Surface runoff quantification of such ungauged semi-arid catchments considered among the important challenges. A 7586 km2 catchment under investigation is located in semi-arid region in central Sudan where mean annual rainfall is around 250 mm and represent the ultimate source for water supply. The objective is to parameterize hydrological characteristics of the catchment and estimate surface runoff using suitable methods and hydrological models that suit the nature of such ungauged catchments with scarce geospatial information. In order to produce spatial runoff estimations, satellite rainfall was used. Remote sensing and GIS were incorporated in the investigations and the generation of landcover and soil information. Five days rainfall event (50.2 mm) was used for the SCS CN model which is considered the suitable for this catchment, as SCS curve number (CN) method is widely used for estimating infiltration characteristics depending on the landcover and soil property. Runoff depths of 3.6, 15.7 and 29.7 mm were estimated for the three different Antecedent Moisture Conditions (AMC-I, AMC-II and AMC-III). The estimated runoff depths of AMCII and AMCIII indicate the possibility of having small artificial surface reservoirs that could provide water for domestic and small household agricultural use.

Hydro-sedimentary flow modelling in some catchments Constantine highlands, case of Wadis Soultez and Reboa (Algeria)

Abstract Erosion is a major phenomenon that causes damage not only to soil and agriculture, but also to the quality of the water amounting to tonnes of matter annually transported on the earth's surface. This fact has attracted the interest of researchers to understand its mechanism and explain its causes and consequences. This work is a comparative study of water erosion in the two semi-arid catchments of Wadi Soultez and Wadi Reboa; located in the North-East of Algeria. The approach adopted for the quantification of sediment transport consists on researching the best regressive model to represent the statistical relation between the sediment yield and the measured water discharge at different scales: annual, seasonal and monthly. The available data cover 27 years from 1985-2012. The results show that the power model has given the best correlation coefficient. Results have indicated that Wadi Reboa transported an average of 14.66 hm3 of water and 0.25 million tonnes of sediments annually. While Wadi Soultez has transported 4.2 hm3 of water and 0.11 million tonnes of sediments annually. At a seasonal scale, sediment amounts have showed significant water erosion in autumn with around 44% and secondarily in the spring with 29% in Wadi Soultez. Unlike Wadi Reboa, sediment transport represents 32% and 46% in autumn and spring respectively. Based on the obtained sediment amounts; it is found that the physical factors: such as steep reliefs, vulnerable lithological nature of rocks and poor vegetal cover, have significantly contributed in accelerating soil erosion.

Climate Change Adaptation in a Mediterranean Semi-Arid Catchment: Testing Managed Aquifer Recharge and Increased Surface Reservoir Capacity

Among different uses of freshwater, irrigation is the most impacting groundwater resource, leading to water table depletion and possible seawater intrusion. The unbalance between the availability of water resources and demand is currently exacerbated and could become worse in the near future in accordance with climate change observations and scenarios provided by Intergovernmental Panel on Climate Change (IPCC). In this context, Increasing Maximum Capacity of the surface reservoir (IMC) and Managed Aquifer Recharge (MAR) are adaptation measures that have the potential to enhance water supply systems resiliency. In this paper, a multiple-users and multiple-resources-Water Supply System (WSS) model is implemented to evaluate the effectiveness of these two adaptation strategies in a context of overexploited groundwater under the RCP 4.5 and the RCP 8.5 IPCC scenarios. The presented a case study that is located in the Puglia, a semi-arid region of South Italy characterized by a conspicuous water demand for irrigation. We observed that, although no significant long-term trend affects the proposed precipitation scenarios, the expected temperature increase highly impacts the WSS resources due to the associated increase of water demand for irrigation purposes. Under the RCP 4.5 the MAR scenario results are more effective than the IMC during long term wet periods (typically 5 years) and successfully compensates the impact on the groundwater resources. Differently, under RCP 8.5, due to more persistent dry periods, both adaptation scenarios fail and groundwater resource become exposed to massive sea water intrusion during the second half of the century. We conclude that the MAR scenario is a suitable adaptation strategy to face the expected future changes in climate, although mitigation actions to reduce green-house gases are strongly required.

Assessing the effect of water harvesting techniques on event-based hydrological responses and sediment yield at a catchment scale in northern Ethiopia using the Limburg Soil Erosion Model (LISEM)

Runoff and sediment yield in semi-arid catchments are highly influenced by infrequent but very heavy rains. These events occur over short temporal scales, so runoff and sediment transport can only be understood using an event-based analysis. We applied a hydrological and soil-erosion model, LISEM, to the Gule catchment (~12km2) in northern Ethiopia. The objectives of the study were: (a) to evaluate the performance of LISEM in describing event-based hydrological processes and sediment yield in a catchment under the influence of different water harvesting techniques (WHTs), and (b) to study the effect of the WHTs on catchment-scale event-based runoff and sediment yield. The model performed satisfactorily (NSE>0.5) for most of the events when discharge was calibrated at the main outlet (Gule) and at a sub-outlet (Misbar). Runoff coefficients for the Gule catchment and Misbar sub-catchment were expectedly low due to the implementation of WHTs, which can store runoff from the rains and increase infiltration into the soil. Simulated and measured sediment yields were of similar orders of magnitude. LISEM generally overestimated sediment yield compared to the measurements. The poor performance of LISEM in predicting sediment yield could be attributed to the uncertainty of several factors controlling soil erosion and the inadequacy of LISEM in describing soil erosion on steep slopes. Catchment-scale model simulations indicated that runoff and sediment yield could be effectively reduced by implementing WHTs. The model estimated 41 and 61% reductions in runoff and sediment yield at the Gule outlet, respectively. Similarly, runoff and sediment yield at the Misbar sub-outlet were reduced by 45 and 48%, respectively. LISEM can thus be used to simulate the effects of different existing or new WHTs on catchment hydrology and sediment yield. The results of scenario predictions could be useful for land-use planners who intend to implement different measures of catchment management.

Assessment of soil particle erodibility and sediment trapping using check dams in small semi-arid catchments

Check dams can be used as a source of information for studies on sediment characteristics and soil particle erodibility. In this study, sediment yield and grain size distribution (GSD) were measured in twenty small catchments draining into a rock check dam in NW Iran for different runoffs during 2010–2011. Significant correlations were found between sediment yield and slope steepness, vegetation cover and soil erodibility factor (K) of the catchments. The erodibility of soil particles was determined using the comparison of GSD between sediment and original soil. Clay was the most erodible soil particle which showed 2.05 times more percentage in sediment than the original soil. The erodibility of soil particles were strongly affected by the rainfall erosivity (EI30). Check dams showed more effectiveness in trapping coarse particles (sand and gravel). The effectiveness of check dams in trapping coarse particles enhanced with increase in the remaining capacity of check dams.

Application of Satellite-Based Precipitation Estimates to Rainfall-Runoff Modelling in a Data-Scarce Semi-Arid Catchment

Rainfall-runoff modelling is a useful tool for water resources management. This study presents a simple daily rainfall-runoff model, based on the water balance equation, which we apply to the 11,630 km2 Lesser Zab catchment in northeast Iraq. The model was forced by either observed daily rain gauge data from four stations in the catchment or satellite-derived rainfall estimates from two TRMM Multi-satellite Precipitation Analysis (TMPA) data products (TMPA-3B42 and 3B42RT) based on the Tropical Rainfall Measuring Mission (TRMM) from 2003 to 2014. As well as using raw TMPA data, we used a bias-correction method to adjust TMPA values based on rain gauge data. The uncorrected TMPA data products underestimated observed mean catchment rainfall by −10.1% and −10.7%. Corrected data also slightly underestimated gauged rainfall by −0.7% and −1.6%, respectively. Nash-Sutcliffe Efficiency (NSE) and Pearson’s Correlation Coefficient (r) for the model fit with the observed hydrograph were 0.75 and 0.87, respectively, for a calibration period (2010–2011) using gauged rainfall data. Model validation performance (2012–2014) was best (highest NSE and r; lowest RMSE and bias) using the corrected 3B42 data product and poorest when driven by uncorrected 3B42RT data. Uncertainty and equifinality were also explored. Our results suggest that TRMM data can be used to drive rainfall-runoff modelling in semi-arid catchments, particularly when corrected using rain gauge data.

A spatio-temporal evaluation of the WRF physical parameterisations for numerical rainfall simulation in semi-humid and semi-arid catchments of Northern China

Mesoscale Numerical Weather Prediction systems can provide rainfall products at high resolutions in space and time, playing an increasingly more important role in water management and flood forecasting. The Weather Research and Forecasting (WRF) model is one of the most popular mesoscale systems and has been extensively used in research and practice. However, for hydrologists, an unsolved question must be addressed before each model application in a different target area. That is, how are the most appropriate combinations of physical parameterisations from the vast WRF library selected to provide the best downscaled rainfall? In this study, the WRF model was applied with 12 designed parameterisation schemes with different combinations of physical parameterisations, including microphysics, radiation, planetary boundary layer (PBL), land-surface model (LSM) and cumulus parameterisations. The selected study areas are two semi-humid and semi-arid catchments located in the Daqinghe River basin, Northern China. The performance of WRF with different parameterisation schemes is tested for simulating eight typical 24-h storm events with different evenness in space and time. In addition to the cumulative rainfall amount, the spatial and temporal patterns of the simulated rainfall are evaluated based on a two-dimensional composed verification statistic. Among the 12 parameterisation schemes, Scheme 4 outperforms the other schemes with the best average performance in simulating rainfall totals and temporal patterns; in contrast, Scheme 6 is generally a good choice for simulations of spatial rainfall distributions. Regarding the individual parameterisations, Single-Moment 6 (WSM6), Yonsei University (YSU), Kain-Fritsch (KF) and Grell-Devenyi (GD) are better choices for microphysics, planetary boundary layers (PBL) and cumulus parameterisations, respectively, in the study area. These findings provide helpful information for WRF rainfall downscaling in semi-humid and semi-arid areas. The methodologies to design and test the combination schemes of parameterisations can also be regarded as a reference for generating ensembles in numerical rainfall predictions using the WRF model.

A method for simultaneous estimation of groundwater evapotranspiration and inflow rates in the discharge area using seasonal water table fluctuations

As an indispensable component of groundwater circulation, groundwater evapotranspiration rate (ETG) estimation using water table fluctuations is a hot research topic in the past decades. However, most existing methods for estimating ETG using either diurnal or seasonal water table fluctuations are based on the assumptions that groundwater recharge equals 0 in the whole period and ETG equals 0 in a specific duration of time. The aim of the current study is to develop a method to estimate ETG when both assumptions do not apply. In the discharge area of a semi-arid catchment in the Ordos Plateau, NW China, it is found out that the water table fluctuations are influenced greatly by recharge, ETG and the vertical inflow rate (qin) in the seasonal scale but are controlled by atmospheric pressure instead of evapotranspiration in the diurnal scale, which implies that the existing methods are not applicable. Therefore, we propose a method to estimate ETG and qin using the seasonal water table fluctuations based on the assumptions that qin is stable throughout the year and ETG has a linear relationship to the reference evapotranspiration (ET0). It is found out that qin is around 1.12mm/d, and the ratio of ETG to ET0 is around 0.4. ETG in April through July with the highest ET0 of the year is found to be around 3mm/d, which is much larger than qin, thus leads to significant decline in water table. ETG in November with the lowest ET0 of the unfrozen period is found to be 1.0mm/d, which is slightly smaller than qin, thus results in a weak increase in water table. The vertical component of regional groundwater inflow is found to exceed precipitation recharge, indicating that regional groundwater inflow plays the major role on maintaining the shallow water table in the discharge area. Groundwater evapotranspiration is also found to exceed precipitation recharge in most time of the unfrozen period, indicating a specified-flux upper boundary with water loss in the discharge area of an arid catchment. The proposed method could be applied to similar study areas and is critical to understand the budget of groundwater.