Groundwater Recharge Assessment using Geographic Information System APLIS Method in Donorojo Karst Area, Pacitan

Karst aquifers can be a source of water supply, especially for the community in Donorojo District, which is one of the areas with the worst drought crisis in Pacitan Regency and requires a sustainable solution to the problem of water needs. Therefore, the study and management of karst formations are very important because of their abundance and potential in forming subsurface aquifer karst aquifers. The recharge rate is one of the basic parameters in the management of the consumption and maintenance of this resource. In addition, the distribution of aquifer locations, aquifer characteristics, and the quality of groundwater forming the aquifer need to be known. This study aims to assess aquifers in the Karst area of Donorojo, Pacitan based on recharge rate and spatial distribution. The APLIS method can estimate surface recharge rates and present the results as a map of the spatial distribution of aquifer recharge rates by utilizing a Geographical Information System (GIS). The results of the analysis using the APLIS method, the groundwater recharge rate in the Donorojo Pacitan karst area is divided into 4 classes, namely very low, low, moderate, and high. Almost the entire Donorojo karst area has a high groundwater recharge rate, this means that the area needs to be used as a groundwater protection zone and it is important to carry out good groundwater management, especially to overcome the problem of drought.

Effect of temporal averaging of meteorological data on predictions of groundwater recharge

Accurate estimates of infiltration and groundwater recharge are critical for many hydrologic, agricultural and environmental applications. Anticipated climate change in many regions of the world, especially in tropical areas, is expected to increase the frequency of high-intensity, short-duration precipitation events, which in turn will affect the groundwater recharge rate. Estimates of recharge are often obtained using monthly or even annually averaged meteorological time series data. In this study we employed the HYDRUS-1D software package to assess the sensitivity of groundwater recharge calculations to using meteorological time series of different temporal resolutions (i.e., hourly, daily, weekly, monthly and yearly averaged precipitation and potential evaporation rates). Calculations were applied to three sites in Brazil having different climatological conditions: a tropical savanna (the Cerrado), a humid subtropical area (the temperate southern part of Brazil), and a very wet tropical area (Amazonia). To simplify our current analysis, we did not consider any land use effects by ignoring root water uptake. Temporal averaging of meteorological data was found to lead to significant bias in predictions of groundwater recharge, with much greater estimated recharge rates in case of very uneven temporal rainfall distributions during the year involving distinct wet and dry seasons. For example, at the Cerrado site, using daily averaged data produced recharge rates of up to 9 times greater than using yearly averaged data. In all cases, an increase in the time of averaging of meteorological data led to lower estimates of groundwater recharge, especially at sites having coarse-textured soils. Our results show that temporal averaging limits the ability of simulations to predict deep penetration of moisture in response to precipitation, so that water remains in the upper part of the vadose zone subject to upward flow and evaporation.

Catchment tracers reveal discharge, recharge and sources of groundwater-borne pollutants in a novel lake modelling approach

Groundwater-borne contaminants such as nutrients, dissolved organic carbon (DOC), coloured dissolved organic matter (CDOM) and pesticides can have an impact the biological quality of lakes. The sources of pollutants can, however, be difficult to identify due to high heterogeneity in groundwater flow patterns. This study presents a novel approach for fast hydrological surveys of small groundwater-fed lakes using multiple groundwater-borne tracers. Water samples were collected from the lake and temporary groundwater wells, installed every 50 m within a distance of 5–45 m to the shore, were analysed for tracer concentrations of CDOM, DOC, total dissolved nitrogen (TDN, groundwater only), total nitrogen (TN, lake only), total dissolved phosphorus (TDP, groundwater only), total phosphorus (TP, lake only), δ18O ∕ δ16O isotope ratios and fluorescent dissolved organic matter (FDOM) components derived from parallel factor analysis (PARAFAC). The isolation of groundwater recharge areas was based on δ18O measurements and areas with a high groundwater recharge rate were identified using a microbially influenced FDOM component. Groundwater discharge sites and the fractions of water delivered from the individual sites were isolated with the Community Assembly via Trait Selection model (CATS). The CATS model utilized tracer measurements of TDP, TDN, DOC and CDOM from the groundwater samples and related these to the tracer measurements of TN, TP, DOC and CDOM in the lake. A direct comparison between the lake and the inflowing groundwater was possible as degradation rates of the tracers in the lake were taken into account and related to a range of water retention times (WRTs) of the lake (0.25–3.5 years in 0.25-year increments). These estimations showed that WRTs above 2 years required a higher tracer concentration of inflowing water than found in any of the groundwater wells around the lake. From the estimations of inflowing tracer concentration, the CATS model isolated groundwater discharge sites located mainly in the eastern part of the lake with a single site in the southern part. Observations from the eastern part of the lake revealed an impermeable clay layer that promotes discharge during heavy precipitation events, which would otherwise be difficult to identify using traditional hydrological methods. In comparison to the lake concentrations, high tracer concentrations in the southern part showed that only a smaller fraction of water could originate from this area, thereby confirming the model results. A Euclidean cluster analysis of δ18O isotopes identified recharge sites corresponding to areas adjacent to drainage channels, and a cluster analysis of the microbially influenced FDOM component C4 further identified five sites that showed a tendency towards high groundwater recharge rate. In conclusion, it was found that this methodology can be applied to smaller lakes within a short time frame, providing useful information regarding the WRT of the lake and more importantly the groundwater recharge and discharge sites around the lake. Thus, it is a tool for specific management of the catchment.

Empirical relationship between soil permeability and resistivity, and its application for determining the groundwater gross recharge in Marpoyan Damai, Pekanbaru, Indonesia

In this article, a novel approach to measuring the gross groundwater recharge rate over a district area is presented. It is based on soil resistivity measurements. It is nondestructive, and allows quick, low-cost mapping of infiltration over large areas, which makes it particularly relevant for land planning. An empirical relationship between soil resistivity, ρ, and permeability, k, must be established before the permeability is related to the gross recharge rate theoretically. The method has been applied in Marpoyan Damai, Pekanbaru, Indonesia. The empirical relationship between resistivity and permeability was determined as log k = −0.1588 ρ + 6.628. The administrative districts within Marpoyan Damai with the lowest and highest annual gross groundwater recharge rate were found to be Wonorejo (about 61,000 m3/a for a 0.3 km2 drainage area) and East Sidomulyo (about 881,000 m3/a for 1.44 km2).

Calculation of Agricultural Drain Spacing Taking into Account Regularity of Water Exchange in the Vadose Zone

Conventional analytical formulae for calculation of subsurface drain spacing for maintaining a desired water table depth in agricultural areas, such as Hooghoudt’s formula, are based on using the ratio between the soil saturated hydraulic conductivity Ks and the groundwater recharge rate q. It is well known that selection of the q value as one of the principle drainage criterion is one of the problems of the drain spacing calculation. In this paper, it is illustrated that for steady state conditions and, in case of homogeneous soil profile, the ratio q/Ks can be substituted by an analytical function that takes into account the regularity of infiltration through the vadose zone. This function can be derived from the soil moisture content in the root zone and other well-known hydrodynamic soil parameters. An example of drain spacing calculation is presented.