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.     

Developing constitutive models from EPR‐based self‐learning finite element analysis

A constitutive model that captures the material behavior under a wide range of loading conditions is essential for simulating complex boundary value problems. In recent years, some attempts have been made to develop constitutive models for finite element analysis using self‐learning simulation (SelfSim). Self‐learning simulation is an inverse analysis technique that extracts material behavior from some boundary measurements (eg, load and displacement). In the heart of the self‐learning framework is a neural network which is used to train and develop a constitutive model that represents the material behavior. It is generally known that neural networks suffer from a number of drawbacks. This paper utilizes evolutionary polynomial regression (EPR) in the framework of SelfSim within an automation process which is coded in Matlab environment. EPR is a hybrid data mining technique that uses a combination of a genetic algorithm and the least square method to search for mathematical equations to represent the behavior of a system. Two strategies of material modeling have been considered in the SelfSim‐based finite element analysis. These include a total stress‐strain strategy applied to analysis of a truss structure using synthetic measurement data and an incremental stress‐strain strategy applied to simulation of triaxial tests using experimental data. The results show that effective and accurate constitutive models can be developed from the proposed EPR‐based self‐learning finite element method. The EPR‐based self‐learning FEM can provide accurate predictions to engineering problems. The main advantages of using EPR over neural network are highlighted.     

Is the age problem resolved?

The cosmological, astrophysical, and nucleocosmochronological methods for estimating the age of the universe and the corresponding uncertainties are comparatively studied in the present paper. We are led to the conclusion that the new measurements of cosmological parameters, and the recent estimates of the age of globular clusters have led to the gradual disappearance of the age problem from the arena of modern cosmology.     

The NAIAD experiment for WIMP searches at Boulby mine and recent results

The NAIAD experiment (NaI Advanced Detector) for weakly interacting massive particle (WIMP) dark matter searches at Boulby mine (UK) is described. The detector consists of an array of encapsulated and unencapsulated NaI(Tl) crystals with high light yield. Six crystals are collecting data at present. Data accumulated by four of them (10.6 kg × year exposure) have been used to set upper limits on the WIMP–nucleon spin-independent and WIMP–proton spin-dependent cross-sections. Pulse shape analysis has been applied to discriminate between nuclear recoils, as may be caused by WIMP interactions, and electron recoils due to gamma background. Various calibrations of crystals are presented.     

Dynamical implications of Jupiter's tropospheric ammonia abundance

Groundbased radio observations indicate that Jupiter's ammonia is globally depleted from 0.6 bars to at least 4-6 bars relative to the deep abundance of ∼3 times solar, a fact that has so far defied explanation. The observations also indicate that (i) the depletion is greater in belts than zones, and (ii) the greatest depletion occurs within Jupiter's local 5-μm hot spots, which have recently been detected at radio wavelengths. Here, we first show that both the global depletion and its belt-zone variation can be explained by a simple model for the interaction of moist convection with Jupiter's cloud-layer circulation. If the global depletion is dynamical in origin, then important endmember models for the belt-zone circulation can be ruled out. Next, we show that the radio observations of Jupiter's 5-μm hot spots imply that the equatorial wave inferred to cause hot spots induces vertical parcel oscillation of a factor of ∼2 in pressure near the 2-bar level, which places important constraints on hot-spot dynamics. Finally, using spatially resolved radio maps, we demonstrate that low-latitude features exceeding ∼4000 km diameter, such as the equatorial plumes and large vortices, are also depleted in ammonia from 0.6 bars to at least 2 bars relative to the deep abundance of 3 times solar. If any low-latitude features exist that contain 3-times-solar ammonia up to the 0.6-bar ammonia condensation level, they must have diameters less than ∼4000 km.     

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.     

SHRIMP zircon dating and Sm/Nd isotopic investigations of Neoproterozoic granitoids, Eastern Desert, Egypt

There is an increasing evidence for the involvement of pre-Neoproterozoic zircons in the Arabian–Nubian Shield, a Neoproterozoic crustal tract that is generally regarded to be juvenile. The source and significance of these xenocrystic zircons are not clear. In an effort to better understand this problem, older and younger granitoids from the Egyptian basement complex were analyzed for chemical composition, SHRIMP U–Pb zircon ages, and Sm–Nd isotopic compositions. Geochemically, the older granitoids are metaluminous and exhibit characteristics of I-type granites and most likely formed in a convergent margin (arc) tectonic environment. On the other hand, the younger granites are peraluminous and exhibit the characteristics of A-type granites; these are post-collisional granites. The U–Pb SHRIMP dating of zircons revealed the ages of magmatic crystallization as well as the presence of slightly older, presumably inherited zircon grains. The age determined for the older granodiorite is 652.5 ± 2.6 Ma, whereas the younger granitoids are 595–605 Ma. Xenocrystic zircons are found in most of the younger granitoid samples; the xenocrystic grains are all Neoproterozoic, but fall into three age ranges that correspond to the ages of other Eastern Desert igneous rocks, viz. 710–690, 675–650 and 635–610 Ma. The analyzed granitoids have (+3.8 to +6.5) and crystallization ages, which confirm previous indications that the Arabian–Nubian Shield is juvenile Neoproterozoic crust. These results nevertheless indicate that older Neoproterozoic crust contributed to the formation of especially the younger granite magmas.     

Travel time based approach for the assessment of vulnerability of karst groundwater: the Transit Time Method

Karst aquifers represent important water resources in many parts of the world. Unfortunately, karst aquifers are characterised by high contamination risks. This paper presents a travel time based method for the estimation of karst groundwater vulnerability. It considers (1) physics-based lateral flow within the uppermost weathered zone (epikarst) in a limestone-dominated region and (2) high velocities of vertical infiltration at discrete infiltration points (e.g. sinkholes) or lines (e.g. dry valleys, faults). Consequently, the Transit Time Method honours the actual flow path within the unsaturated zone of a karst aquifer system. A test site in Northern Jordan was chosen for the demonstration of the assessment technique, i.e. the catchment area of the Qunayyah Spring north of the capital Amman. The results demonstrate that zones of highest vulnerability lie within valleys and nearby main fault zones. It also reveals that regions, categorised as protected areas by other methods due to thick unsaturated zones, contribute to a major degree to the total risk.