Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.     

Spatial measurements of stream baseflow,a relevant method for aquifer characterization and permeability evaluation. Application to a hard‐rock aquifer,the Oman ophiolite

The structure, functioning and hydrodynamic properties of aquifers can be determined from an analysis of the spatial variability of baseflow in the streams with which they are associated. Such analyses are based on simple low‐cost measurements. Through interpreting the hydrological profiles (Q = f(A)) it is possible to locate the aquifer(s) linked to the stream network and to determine the type of interrelated flow, i.e. whether the stream drains or feeds the aquifer. Using an analytical solution developed for situations with a positive linear relationship, i.e. where the baseflow increases linearly with increasing catchment size, it is also possible to estimate the permeability of the aquifer(s) concerned at catchment scale. Applied to the hard‐rock aquifers of the Oman ophiolite, this method shows that the ‘gabbro’ aquifer is more permeable than the ‘peridotite’ aquifer. As a consequence the streams drain the peridotites and ‘leak’ into the gabbro. The hydrological profiles within the peridotite are linear and positive, and indicate homogeneity in the hydrodynamic properties of these formations at the kilometre scale. The permeability of the peridotite is estimated at 5 · 10^?7 to 5 · 10^?8 m/s. Copyright © 2004 John Wiley & Sons, Ltd.     

Multiple longshore sand bars in the upper Chesapeake Bay

In the upper Chesapeake Bay (Maryland, U.S.A.) field surveys were conducted at 18 multiple longshore sand bar sites. The multiple bar systems were found in water depths less than approximately 2 m (mean sea level), and exhibited mild bottom slopes of 0·0052 or less. The number of bars composing each system ranged from four to 17 and the spacing between the crests typically increased in the offshore direction, ranging from 12 to 70 m. Bar height also typically increased with distance offshore and ranged from 0·03 to 0·61 m. A grain size analysis of crest and trough sediment did not reveal any significant differences and the sediment was categorized as ‘fine sand’. A review of the literature data indicated that the Chesapeake Bay multiple bars possessed similar characteristics to those found in Gelding Bay (Baltic Sea); similarities in fetch, wave height and tidal range between the two bays may account for this finding. The surf-scaling parameter indicated that the multiple bar systems were extremely dissipative with regard to wave energy, and wave height appeared to be an important factor in controlling bar spacing and bar height. A multiple wave break point hypothesis was discussed as a possible mechanism for the formation of Chesapeake Bay multiple longshore bars, and limited observational evidence appeared to support such a mechanism.     

Particulate organic matter in the coastal and estuarine waters of Goa and its relationship with phytoplankton production

In the coastal and estuarine waters of Goa, particulate organic carbon (POC) varied from 0.52 to 2.51 mg l^?1 and from 0.28 to 5.24 mg l^?1 and particulate phosphorus (PP) varied from 0.71 to 5.18 μg l^?1 and from 0.78 to 20.34 μg l^?1, respectively. The mean values of chlorophyll and primary productivity were 1.94 mg m^?3 and 938.1 mg C m^?2 day^?1 in the coastal waters and 4.3 mg m^?3 and 636.5 mg C m^?1 day^?1 in the estuarine waters, respectively.$\text{POC}\text{chl}$ ratios were low in June and October even when POC values were quite high. The POC in surface waters was linearly correlated with the chlorophyll content. Also PP increased when chlorophyll and primary productivity remained high. The results suggest that the phytoplankton was sharply increasing and contributed to POC and PP content. The percentage of detritus calculated from the intercept values of chlorophyll on POC varied from 46 to 76% depending on season. Results indicate that the major portion of POC and PP during postmonsoon (October–January) is derived from phytoplankton production while the allochthonous matter predominate during monsoon (June–September).     

定水头注水引起的含水层水平运动和应变

基于含水层固体颗粒与孔隙水不可压缩的假设 ,本文导出了单井注水情况下泰斯承压含水层水平运动速度与水头之间的基本关系式。然后利用注水井壁处的应力、应变边界条件 ,进一步导出了单井定水头注水引起的泰斯承压含水层水平运动速度、位移和应变解析表达式。该水平位移与应变由两部分组成 :一部分为由注水压力本身引起的经典弹性力学解项 ,它仅随半径而变化 ,与注水时间无关 ;另一部分为由地下水头变化引起的水动力学位移和应变解项。其中 ,含水层水动力学水平位移随时间加长呈指数增长特征 ,水动力学径向应变则表现为近井处拉张、远井处挤压的分区特征 ,且近井拉张区随时间加长逐渐向外扩展。单井注水含水层水动力学水平位移、应变解的导出 ,完善和发展了单孔内压经典弹性平面力学问题解     

Implications of Kali-Hindon inter-stream aquifer water balance for groundwater management in western Uttar Pradesh

The Kali-Hindon inter-stream region extends over an area of 395 km² within the Ganga-Yamuna interfluve. It is a fertile tract for sugarcane cultivation. Groundwater is a primary resource for irrigation and industrial purposes. In recent years, over-exploitation has resulted in an adverse impact on the groundwater regime. In this study, an attempt has been made to calculate a water balance for the Kali-Hindon inter-stream region. Various inflows and outflows to and from the aquifer have been calculated. The recharge due to rainfall and other recharge parameters such as horizontal inflow, irrigation return flow and canal seepage were also evaluated. Groundwater withdrawals, evaporation from the water table, discharge from the aquifer to rivers and horizontal subsurface outflows were also estimated. The results show that total recharge into the system is 148.72 million cubic metres (Mcum), whereas the total discharge is 161.06 Mcum, leaving a deficit balance of −12.34 Mcum. Similarly, the groundwater balance was evaluated for the successive four years. The result shows that the groundwater balance is highly sensitive to variation in rainfall followed by draft through pumpage. The depths to water level are shallow in the canal-irrigated northern part of the basin and deeper in the southern part. The pre-monsoon and post-monsoon water levels range from 4.6 to 17.7 m below ground level (bgl) and from 3.5 to 16.5 m bgl respectively. It is concluded that the groundwater may be pumped in the canal-irrigated northern part, while withdrawals may be restricted to the southern portion of the basin, where intense abstraction has led to rapidly falling water table levels.     

A comparative study of four vulnerability mapping methods in a detritic aquifer under mediterranean climatic conditions

This paper presents the results of a comparative study relating to the application of four vulnerability mapping methods, GOD, AVI, DRASTIC and SINTACS, in a pilot detritic aquifer situated in NW Morocco, known as the Martil–Alila aquifer. The principal objective of this work is to determine the most suitable such methods for this aquifer type within a Mediterranean context, and to show the effect of the rainfall variations that are characteristic of the Mediterranean climate on the degree of vulnerability. The methods applied distinguish five classes of vulnerability, these being irregularly divided up in space, with the division varying according to the method in question. The vulnerability maps obtained by the different methods strongly suggest that the eastern half of the aquifer is more vulnerable to contamination than the western half, for all hydrological situations. The effect of climatic conditions on the degree of vulnerability is well represented by the DRASTIC, according to which the aquifer is moderately to strongly vulnerable during humid hydrological years and weakly to moderately vulnerable during dry ones. For the other methods, this climatic effect is limited to the area occupied by the two predominant classes (“High” and “Low” for GOD and “High” and “Moderate” for SINTACS) while it is null for AVI. In conclusion, DRASTIC appears the most suitable for mapping the vulnerability to contamination of Mediterranean coastal detritic aquifers such as the Martil–Alila aquifer.     

Deciphering potential groundwater zone in hard rock through the application of GIS

Remote Sensing and Geographic Information System has become one of the leading tools in the field of hydrogeological science, which helps in assessing, monitoring and conserving groundwater resources. It allows manipulation and analysis of individual layer of spatial data. It is used for analysing and modelling the interrelationship between the layers. This paper mainly deals with the integrated approach of Remote Sensing and geographical information system (GIS) to delineate groundwater potential zones in hard rock terrain. The remotely sensed data at the scale of 1:50,000 and topographical information from available maps, have been used for the preparation of ground water prospective map by integrating geology, geomorphology, slope, drainage-density and lineaments map of the study area. Further, the data on yield of aquifer, as observed from existing bore wells in the area, has been used to validate the groundwater potential map. The final result depicts the favourable prospective zones in the study area and can be helpful in better planning and management of groundwater resources especially in hard rock terrains.     

Modeling of water–rock interactions in an aquitard of sandy clay in the coastal area near Beihai,China

Unconsolidated sand, gravel and clay deposits near Beihai and in the Leizhou Peninsula in southern China form an unconfined aquifer, aquitard and a confined aquifer. Water and soil samples were collected from the two aquifers in the coastal Beihai area for the determination of chemical compositions, minerals and soluble ions. Hydrogeochemical modeling of three flow paths through the aquitard are carried out using PHREEQC to determine water–rock interactions along the flow paths. The results indicate that the dissolution of anorthite, fluorite, halite, rhodochrosite and CO₂, and precipitation of potash feldspar and kaolinite may be occurring when groundwater leaks through the aquitard from the unconfined aquifer to the confined aquifer. Cation exchanges between Na and Ca can also happen along the flow paths.     

Groundwater recharge of carbonate aquifers of the Silesian-Cracow Triassic (southern Poland) under human impact

A Triassic carbonate unit has been intensively drained by zinc and lead ore mines and numerous borehole fields since the nineteenth century. Its groundwater recharge has increased due to: pumping of water from boreholes, mining activity, and urbanization. An approach to determine the amounts of the recharge at a variety of spatial scales is presented in the paper. Different methods were used to identify and quantify recharge components on a regional and local scale: mathematical modelling was performed for four aquifers included in an aquifer system, an analytical estimation based on the assumption that an average recharge is equal to the average discharge of the hydrogeological system—for six man-made drainage centres, and the method of water level fluctuation (WLF) was applied in one observation borehole. Results of modelling have been supplemented by observation of environmental tracers (δ¹⁸O, δ²H, ³H), noble gases temperatures, and ⁴He_exc in groundwater. The regional aquifer’s current recharge according to estimations performed by means of modelling varies from 39 to 101 mm/year on average. Depending on the aquifer site the average precipitation ranges from 779 to 864 mm/year. In the confined part of the aquifer average recharge ranges from 26 to 61 mm/year. Within outcrops average recharge varies from 96 to 370 mm/year. Current recharge estimated by the analytical method for man-made drainage centres varies from 158 up to 440 mm/year. High values are caused by different recharge sources like precipitation, induced leakage from shallow aquifers, and water losses from streams, water mains and sewer systems. Pumping of water, mining and municipal activities constitute additional factors accounting for the intensified recharge.