Numerical parametric study on the influence of location and inclination of large-scale asperities on the shear strength of concreterock interfaces of small buttress dams

When assessing the sliding stability of a concrete dam,the influence of large-scale asperities in the sliding plane is often ignored due to limitations of the analytical rigid body assessment methods provided by current dam assessment guidelines.However,these asperities can potentially improve the load capacity of a concrete dam in terms of sliding stability.Although their influence in a sliding plane has been thoroughly studied for direct shear,their influence under eccentric loading,as in the case of dams,is unknown.This paper presents the results of a parametric study that used finite element analysis(FEA)to investigate the influence of large-scale asperities on the load capacity of small buttress dams.By varying the inclination and location of an asperity located in the concrete-rock interface along with the strength of the rock foundation material,transitions between different failure modes and correlations between the load capacity and the varied parameters were observed.The results indicated that the inclination of the asperity had a significant impact on the failure mode.When the inclinationwas 30and greater,interlocking occurred between the dam and foundation and the governing failure modes were either rupture of the dam body or asperity.When the asperity inclination was significant enough to provide interlocking,the load capacity of the dam was impacted by the strength of the rock in the foundation through influencing the load capacity of the asperity.The location of the asperity along the concrete-rock interface did not affect the failure mode,except for when the asperity was located at the toe of the dam,but had an influence on the load capacity when the failure occurred by rupture of the buttress or by sliding.By accounting for a single large-scale asperity in the concrete-rock interface of the analysed dam,a horizontal load capacity increase of 30%e160%was obtained,depending on the inclination and location of the asperity and the strength of the foundation material.

Towards implementation of alloy-specific thermo-fluid modelling for laser powder-bed fusion of Mg alloys

Multi-physics thermo-fluid modeling has been extensively used as an approach to understand melt pool dynamics and defect formation as well as optimizing the process-related parameters of laser powder-bed fusion(L-PBF).However,its capabilities for being implemented as a reliable tool for material design,where minor changes in material-related parameters must be accurately captured,is still in question.In the present research,first,a thermo-fluid computational fluid dynamics(CFD)model is developed and validated against experimental data.Considering the predicted material properties of the pure Mg and commercial ZK60 and WE43 Mg alloys,parametric studies are done attempting to elucidate how the difference in some of the material properties,i.e.,saturated vapor pressure,viscosity,and solidification range,can influence the melt pool dynamics.It is found that a higher saturated vapor pressure,associated with the ZK60 alloy,leads to a deeper unstable keyhole,increasing the keyhole-induced porosity and evaporation mass loss.Higher viscosity and wider solidification range can increase the non-uniformity of temperature and velocity distribution on the keyhole walls,resulting in increased keyhole instability and formation of defects.Finally,the WE43 alloy showed the best behavior in terms of defect formation and evaporation mass loss,providing theoretical support to the extensive use of this alloy in L-PBF.In summary,this study suggests an approach to investigate the effect of materials-related parameters on L-PBF melting and solidification,which can be extremely helpful for future design of new alloys suitable for L-PBF.

Integration of bio-inspired limb-like structure damping into motor suspension of high-speed trains to enhance bogie hunting stability

Hunting stability is an important performance criterion in railway vehicles.This study proposes an incorporation of a bio-inspired limb-like structure(LLS)-based nonlinear damping into the motor suspension system for traction units to improve the nonlinear critical speed and hunting stability of high-speed trains(HSTs).Initially,a vibration transmission analysis is conducted on a HST vehicle and a metro vehicle that suffered from hunting motion to explore the effect of different motor suspension systems from on-track tests.Subsequently,a simplified lateral dynamics model of an HST bogie is established to investigate the influence of the motor suspension on the bogie hunting behavior.The bifurcation analysis is applied to optimize the motor suspension parameters for high critical speed.Then,the nonlinear damping of the bio-inspired LLS,which has a positive correlation with the relative displacement,can further improve the modal damping of hunting motion and nonlinear critical speed compared with the linear motor suspension system.Furthermore,a comprehensive numerical model of a high-speed train,considering all nonlinearities,is established to investigate the influence of different types of motor suspension.The simulation results are well consistent with the theoretical analysis.The benefits of employing nonlinear damping of the bio-inspired LLS into the motor suspension of HSTs to enhance bogie hunting stability are thoroughly validated.

The First Five Years of a Phase Theory for Complex Systems and Networks

In this paper,we review the development of a phase theory for systems and networks in its first five years,represented by a trilogy:Matrix phases and their properties;The MIMO LTI system phase response,its physical interpretations,the small phase theorem,and the sectored real lemma;The synchronization of a multi-agent network using phase alignment.Towards the end,we also summarize a list of ongoing research on the phase theory and speculate what will happen in the next five years.

MOF-Derived Iron-Cobalt Phosphide Nanoframe as Bifunctional Electrocatalysts for Overall Water Splitting

Transition metal phosphides(TMPs)have emerged as an alternative to precious metals as efficient and low-cost catalysts for water electrolysis.Elemental doping and morphology control are effective approaches to further improve the performance of TMPs.Herein,Fe-doped CoP nanoframes(Fe-CoP NFs)with specific open cage configuration were designed and synthesized.The unique nano-framework structured Fe-CoP material shows overpotentials of only 255 and 122 mV at 10 mA cm^(−2)for oxygen evolution reaction(OER)and hydrogen evolution reaction(HER),respectively,overwhelming most transition metal phosphides.For overall water splitting,the cell voltage is 1.65 V for Fe-CoP NFs at a current density of 10 mA cm^(−2),much superior to what is observed for the classical nanocubic structures.Fe-CoP NFs show no activity degradation up to 100 h which contrasts sharply with the rapidly decaying performance of noble metal catalyst reference.The superior electrocatalytic performance of Fe-CoP NFs due to abundant accessible active sites,reduced kinetic energy barrier,and preferable*O-containing intermediate adsorption is demonstrated through experimental observations and theoretical calculations.Our findings could provide a potential method for the preparation of multifunctional material with hollow structures and offer more hopeful prospects for obtaining efficient earth-abundant catalysts for water splitting.

Evolution of nonmetallic inclusions in 80-t 9CrMoCoB large-scale ingots during electroslag remelting process

In combination with theoretical calculations,experiments were conducted to investigate the evolution behavior of nonmetallic inclusions(NMIs)during the manufacture of large-scale heat-resistant steel ingots using 9CrMoCoB heat-resistant steel and CaF_(2)–CaO–Al_(2)O_(3)–SiO_(2)–B_(2)O_(3)electroslag remelting(ESR)-type slag in an 80-t industrial ESR furnace.The main types of NMI in the consumable electrode comprised pure alumina,a multiphase oxide consisting of an Al_(2)O_(3)core and liquid CaO–Al_(2)O_(3)–SiO_(2)–MnO shell,and M_(23)C_(6)carbides with an MnS core.The Al_(2)O_(3)and MnS inclusions had higher precipitation temperatures than the M_(23)C_(6)-type carbide under equilibrium and nonequilibrium solidification processes.Therefore,inclusions can act as nucleation sites for carbide layer precipitation.The ESR process completely removed the liquid CaO–Al_(2)O_(3)–SiO_(2)–MnO oxide and MnS inclusion with a carbide shell,and only the Al_(2)O_(3)inclusions and Al_(2)O_(3)core with a carbide shell occupied the remelted ingot.The M_(23)C_(6)-type carbides in steel were determined as Cr_(23)C_(6)based on the analysis of transmission electron microscopy results.The substitution of Cr with W,Fe,or/and Mo in the Cr_(23)C_(6)lattice caused slight changes in the lattice parameter of the Cr_(23)C_(6)carbide.Therefore,Cr_(21.34)Fe_(1.66)C_(6),(Cr_(19)W_(4)C_(6),Cr_(18.4)Mo_(4.6)C_(6),and Cr_(16)Fe_(5)Mo_(2)C_(6)can match the fraction pattern of Cr_(23)C_(6)carbide.The Al_(2)O_(3)inclusions in the remelted ingot formed due to the reduction of CaO,SiO_(2),and MnO components in the liquid inclusion.The increased Al content in liquid steel or the higher supersaturation degree of Al_(2)O_(3)precipitation in the remelted ingot than that in the electrode can be attributed to the evaporation of CaF_(2)and the increase in CaO content in the ESR-type slag.

Ab initio density response and local field factor of warm dense hydrogen

We present quasi-exact ab initio path integral Monte Carlo(PIMC)results for the partial static density responses and local field factors of hydrogen in the warm dense matter regime,from solid density conditions to the strongly compressed case.The full dynamic treatment of electrons and protons on the same footing allows us to rigorously quantify both electronic and ionic exchange–correlation effects in the system,and to compare the results with those of earlier incomplete models such as the archetypal uniform electron gas or electrons in a fixed ion snapshot potential that do not take into account the interplay between the two constituents.The full electronic density response is highly sensitive to electronic localization around the ions,and our results constitute unambiguous predictions for upcoming X-ray Thomson scattering experiments with hydrogen jets and fusion plasmas.All PIMC results are made freely available and can be used directly for a gamut of applications,including inertial confinement fusion calculations and the modeling of dense astrophysical objects.Moreover,they constitute invaluable benchmark data for approximate but computationally less demanding approaches such as density functional theory or PIMC within the fixed-node approximation.

Effect of hafnium and molybdenum addition on inclusion characteristics in Co-based dual-phase high-entropy alloys

Specific grades of high-entropy alloys(HEAs)can provide opportunities for optimizing properties toward high-temperature applications.In this work,the Co-based HEA with a chemical composition of Co_(47.5)Cr_(30)Fe_(7.5)Mn_(7.5)Ni_(7.5)(at%)was chosen.The refractory metallic elements hafnium(Hf)and molybdenum(Mo)were added in small amounts(1.5at%)because of their well-known positive effects on high-temperature properties.Inclusion characteristics were comprehensively explored by using a two-dimensional cross-sectional method and extracted by using a three-dimensional electrolytic extraction method.The results revealed that the addition of Hf can reduce Al_(2)O_(3)inclusions and lead to the formation of more stable Hf-rich inclusions as the main phase.Mo addition cannot influence the inclusion type but could influence the inclusion characteristics by affecting the physical parameters of the HEA melt.The calculated coagulation coefficient and collision rate of Al_(2)O_(3)inclusions were higher than those of HfO_(2)inclusions,but the inclusion amount played a larger role in the agglomeration behavior of HfO_(2)and Al_(2)O_(3)inclusions.The impurity level and active elements in HEAs were the crucial factors affecting inclusion formation.

Integrating phosphorus management and cropping technology for sustainable maize production

Achieving high maize yields and efficient phosphorus(P)use with limited environmental impacts is one of the greatest challenges in sustainable maize production.Increasing plant density is considered an effective approach for achieving high maize yields.However,the low mobility of P in soils and the scarcity of natural P resources have hindered the development of methods that can simultaneously optimize P use and mitigate the P-related environmental footprint at high plant densities.In this study,meta-analysis and substance flow analysis were conducted to evaluate the effects of different types of mineral P fertilizer on maize yield at varying plant densities and assess the flow of P from rock phosphate mining to P fertilizer use for maize production in China.A significantly higher yield was obtained at higher plant densities than at lower plant densities.The application of single superphosphate,triple super-phosphate,and calcium magnesium phosphate at high plant densities resulted in higher yields and a smaller environmental footprint than the application of diammonium phosphate and monoammonium phosphate.Our scenario analyses suggest that combining the optimal P type and application rate with a high plant density could increase maize yield by 22%.Further,the P resource use efficiency throughout the P supply chain increased by 39%,whereas the P-related environmental footprint decreased by 33%.Thus,simultaneously optimizing the P type and application rate at high plant densities achieved multiple objectives during maize production,indicating that combining P management with cropping techniques is a practical approach to sustainable maize production.These findings offer strategic,synergistic options for achieving sustainable agricultural development.

Development of track geometry inspection equipment for high-speed comprehensive inspection train in China

Purpose–This study aims to analyze the development direction of track geometry inspection equipment for high-speed comprehensive inspection train in China.Design/methodology/approach–The development of track geometry inspection equipment for highspeed comprehensive inspection train in China in the past 20 years can be divided into 3 stages.Track geometry inspection equipment 1.0 is the stage of analog signal.At the stage 1.0,the first priority is to meet the China’s railways basic needs of pre-operation joint debugging,safety assessment and daily dynamic inspection,maintenance and repair after operation.Track geometry inspection equipment 2.0 is the stage of digital signal.At the stage 2.0,it is important to improve stability and reliability of track geometry inspection equipment by upgrading the hardware sensors and improving software architecture.Track geometry inspection equipment 3.0 is the stage of lightweight.At the stage 3.0,miniaturization,low power consumption,self-running and green economy are co-developing on demand.Findings–The ability of track geometry inspection equipment for high-speed comprehensive inspection train will be expanded.The dynamic inspection of track stiffness changes will be studied under loaded and unloaded conditions in response to the track local settlement,track plate detachment and cushion plate failure.The dynamic measurement method of rail surface slope and vertical curve radius will be proposed,to reveal the changes in railway profile parameters of high-speed railways and the relationship between railway profile,track irregularity and subsidence of subgrade and bridges.The 200 m cut-off wavelength of track regularity will be researched to adapt to the operating speed of 400 km/h.Originality/value–The research can provide new connotations and requirements of track geometry inspection equipment for high-speed comprehensive inspection train in the new railway stage.

Simultaneous determination of typical toxic aldehydes formed during food frying and digestion using isotope dilution UHPLC-QqQ-MS/MS

An isotope dilution ultra-performance liquid chromatography-triple quadrupole mass spectrometry method was developed to simultaneously detect two typical kinds ofα,β-unsaturated aldehydes,namely 4-hydroxy-2-hexenal(4-HHE)and 4-hydroxy-2-nonenal(4-HNE),in foods.The proposed method exhibited a linear range of 10-1000 ng/mL with a limit of detection of 0.1-2.0 ng/g and a limit of quantification of 0.3-5.0 ng/g.The recovery rates of these typical toxic aldehydes(i.e.,4-HHE,4-HNE)and their d3-labeled analogues were 91.54%-105.12%with a low matrix effect.Furthermore,this proposed method was successfully applied to a real frying system and a simulated digestion system,wherein the contents of 4-HHE and 4-HNE were determined for both.Overall,the obtained results provide strong support for further research into the production of 4-HHE and 4-HNE resulting from foods during oil digestion and frying.

Effect of gas blowing nozzle angle on multiphase flow and mass transfer during RH refining process

A three-dimensional mathematical model was developed to investigate the effect of gas blowing nozzle angles on multiphase flow,circulation flow rate,and mixing time during Ruhrstahl-Heraeus(RH) refining process.Also,a water model with a geometric scale of 1:4 from an industrial RH furnace of 260 t was built up,and measurements were carried out to validate the mathematical model.The results show that,with a conventional gas blowing nozzle and the total gas flow rate of 40 L·min^(-1),the mixing time predicted by the mathematical model agrees well with the measured values.The deviations between the model predictions and the measured values are in the range of about 1.3%–7.3% at the selected three monitoring locations,where the mixing time was defined as the required time when the dimensionless concentration is within 3% deviation from the bath averaged value.In addition,the circulation flow rate was 9 kg·s^(-1).When the gas blowing nozzle was horizontally rotated by either 30° or 45°,the circulation flow rate was found to be increased by about 15% compared to a conventional nozzle,due to the rotational flow formed in the up-snorkel.Furthermore,the mixing time at the monitoring point 1,2,and 3 was shortened by around 21.3%,28.2%,and 12.3%,respectively.With the nozzle angle of 30° and 45°,the averaged residence time of 128 bubbles in liquid was increased by around 33.3%.

Two-dimensional plane strain consolidation of unsaturated soils considering the depth-dependent stress

In practical engineering,the total vertical stress in the soil layer is not constant due to stress diffusion,and varies with time and depth.Therefore,the purpose of this paper is to investigate the effect of stress diffusion on the two-dimensional(2D)plane strain consolidation properties of unsaturated soils when the stress varies with time and depth.A series of semi-analytical solutions in terms of excess pore air and water pressures and settlement for 2D plane strain consolidation of unsaturated soils can be derived with the joint use of Laplace transform and Fourier sine series expansion.Then,the inverse Laplace transform of the semi-analytical solution is given in the time domain using a self-programmed code based on Crump’s method.The reliability of the obtained solutions is proved by the degeneration.Finally,the 2D plots of excess pore pressures and the curves of settlement varying with time,considering different physical parameters of unsaturated soil stratum and depth-dependent stress,are depicted and analyzed to study the 2D plane strain consolidation properties of unsaturated soils subjected to the depthdependent stress.

Review of advanced road materials, structures, equipment, and detection technologies

As a vital and integral component of transportation infrastructure,pavement has a direct and tangible impact on socio-economic sustainability.In recent years,an influx of groundbreaking and state-of-the-art materials,structures,equipment,and detection technologies related to road engineering have continually and progressively emerged,reshaping the landscape of pavement systems.There is a pressing and growing need for a timely summarization of the current research status and a clear identification of future research directions in these advanced and evolving technologies.Therefore,Journal of Road Engineering has undertaken the significant initiative of introducing a comprehensive review paper with the overarching theme of“advanced road materials,structures,equipment,and detection technologies”.This extensive and insightful review meticulously gathers and synthesizes research findings from 39 distinguished scholars,all of whom are affiliated with 19 renowned universities or research institutions specializing in the diverse and multidimensional field of highway engineering.It covers the current state and anticipates future development directions in the four major and interconnected domains of road engineering:advanced road materials,advanced road structures and performance evaluation,advanced road construction equipment and technology,and advanced road detection and assessment technologies.

Effect of rock joints on lined rock caverns subjected to high internal gaspressure

The storage of hydrogen gas in lined rock caverns(LRCs)may enable the implementation of the firstlarge-scale fossil-free steelmaking process in Sweden,but filling such storage causes joints in the rockmass to open,concentrating strains in the lining.The structural interaction between the LRC componentsmust be able to reduce the strain concentration in the sealing steel lining;however,this interaction iscomplex and difficult to predict with analytical methods.In this paper,the strain concentration in LRCsfrom the opening of rock joints is studied using finite element(FE)analyses,where the large-and small-scale deformation behaviors of the LRC are coupled.The model also includes concrete crack initiation anddevelopment with increasing gas pressure and rock joint width.The interaction between the jointed rockmass and the reinforced concrete,the sliding layer,and the steel lining is demonstrated.The results showthat the rock mass quality and the spacing of the rock joints have the greatest influence on the straindistributions in the steel lining.The largest effect of rock joints on the maximum strains in the steellining was observed for geological conditions of“good”quality rock masses.

Modelling erosion of a single rock block using a coupled CFD-DEM approach

Rock block removal is the prevalent physical mechanism for rock erosion and could affect the stability of dam foundations and spillways.Despite this,understanding of block removal is still inadequate because of the complex interactions among block characteristics,hydraulic forces,and erosive processes acting on the block.Herein,based on a previously conducted physical experiment of erosion of a single rock block,the removal processes of two different protruding blocks are represented by a coupled computational fluid dynamics-discrete element model(CFD-DEM)approach under varied flow conditions.Additionally,the blocks could be rotated with respect to the flow direction to consider the effect of the discontinuity orientation on the block removal process.Simulation results visualize the entire block removal process.The simulations reproduce the effects of the discontinuity orientation on the critical flow velocity inducing block incipient motion and the trajectory of the block motion observed in the physical experiments.The numerical results present a similar tendency of the critical velocities at different discontinuity orientations but have slightly lower values.The trajectory of the block in the simulations fits well with the experimental measurements.The relationship between the dimensionless critical shear stress and discontinuity orientation observed from the simulations shows that the effect of block protrusion becomes more dominant on the block incipient motion with the increase of relative protrusion height.To our knowledge,this present study is the first attempt to use the coupled finite volume method(FVM)-DEM approach for modelling the interaction behavior between the block and the flowing water so that the block removal process can be reproduced and analyzed.

A novel method for geometric quality assurance of rock joint replicas in direct shear testing-Part 2:Validation and mechanical replicability

Each rock joint is unique by nature which means that utilization of replicas in direct shear tests is required in experimental parameter studies.However,a method to acquire knowledge about the ability of the replicas to imitate the shear mechanical behavior of the rock joint and their dispersion in direct shear testing is lacking.In this study,a novel method is presented for geometric quality assurance of replicas.The aim is to facilitate generation of high-quality direct shear testing data as a prerequisite for reliable subsequent analyses of the results.In Part 1 of this study,two quality assurance parameters,smf and V_(Hp100),are derived and their usefulness for evaluation of geometric deviations,i.e.geometric reproducibility,is shown.In Part 2,the parameters are validated by showing a correlation between the parameters and the shear mechanical behavior,which qualifies the parameters for usage in the quality assurance method.Unique results from direct shear tests presenting comparisons between replicas and the rock joint show that replicas fulfilling proposed threshold values of σ_(mf)<0.06 mm and
|V_(Hp100)|
<0.2 mm have a narrow dispersion and imitate the shear mechanical behavior of the rock joint in all aspects apart from having a slightly lower peak shear strength.The wear in these replicas,which have similar morphology as the rock joint,is in the same areas as in the rock joint.The wear is slightly larger in the rock joint and therefore the discrepancy in peak shear strength derives from differences in material properties,possibly from differences in toughness.It is shown by application of the suggested method that the quality assured replicas manufactured following the process employed in this study phenomenologically capture the shear strength characteristics,which makes them useful in parameter studies.

A novel method for geometric quality assurance of rock joint replicas in direct shear testing-Part 1:Derivation of quality assurance parameters and geometric reproducibility

Since each rock joint is unique by nature,the utilization of replicas in direct shear testing is required to carry out experimental parameter studies.However,information about the ability of the replicas to simulate the shear mechanical behavior of the rock joint and their dispersion in direct shear testing is lacking.With the aim to facilitate generation of high-quality direct shear test data from replicas,a novel component in the testing procedure is introduced by presenting two parameters for geometric quality assurance.The parameters are derived from surface comparisons of three-dimensional(3D)scanning data of the rock joint and its replicas.The first parameter,smf,captures morphological deviations between the replica and the rock joint surfaces.smf is derived as the standard deviation of the deviations between the coordinate points of the replica and the rock joint.Four sources of errors introduced in the replica manufacturing process employed in this study could be identified.These errors could be minimized,yielding replicas with smf
0.06 mm.The second parameter is a vector,VHp100,which describes deviations with respect to the shear direction.It is the projection of the 100 mm long normal vector of the best-fit plane of the replica joint surface to the corresponding plane of the rock joint.

VHp100

was found to be less than or equal to 0.36 mm in this study.Application of these two geometric quality assurance parameters demonstrates that it is possible to manufacture replicas with high geometric similarity to the rock joint.In a subsequent paper(part 2),smf and VHp100 are incorporated in a novel quality assurance method,in which the parameters shall be evaluated prior to direct shear testing.Replicas having parameter values below established thresholds shall have a known and narrow dispersion and imitate the shear mechanical behavior of the rock joint.

A Parallel Hybrid Testing Technique for Tri-Programming Model-Based Software Systems

Recently,researchers have shown increasing interest in combining more than one programming model into systems running on high performance computing systems(HPCs)to achieve exascale by applying parallelism at multiple levels.Combining different programming paradigms,such as Message Passing Interface(MPI),Open Multiple Processing(OpenMP),and Open Accelerators(OpenACC),can increase computation speed and improve performance.During the integration of multiple models,the probability of runtime errors increases,making their detection difficult,especially in the absence of testing techniques that can detect these errors.Numerous studies have been conducted to identify these errors,but no technique exists for detecting errors in three-level programming models.Despite the increasing research that integrates the three programming models,MPI,OpenMP,and OpenACC,a testing technology to detect runtime errors,such as deadlocks and race conditions,which can arise from this integration has not been developed.Therefore,this paper begins with a definition and explanation of runtime errors that result fromintegrating the three programming models that compilers cannot detect.For the first time,this paper presents a classification of operational errors that can result from the integration of the three models.This paper also proposes a parallel hybrid testing technique for detecting runtime errors in systems built in the C++programming language that uses the triple programming models MPI,OpenMP,and OpenACC.This hybrid technology combines static technology and dynamic technology,given that some errors can be detected using static techniques,whereas others can be detected using dynamic technology.The hybrid technique can detect more errors because it combines two distinct technologies.The proposed static technology detects a wide range of error types in less time,whereas a portion of the potential errors that may or may not occur depending on the 4502 CMC,2023,vol.74,no.2 operating environment are left to the dynamic technology,which completes the validation.

Groundwater quality in the vicinity of a dumpsite in Lagos metropolis,Nigeria

Inappropriate management of municipal solid waste dumpsites is a major cause of groundwater contamination in developing countries,but the extent of the problem is not known.This study investigated groundwater quality in the vicinity of Olusosun dumpsite in Lagos,Nigeria,the most populous city in sub-Saharan Africa.During 2020,monthly groundwater samples were collected in 17 wells and boreholes used as drinking water sources,and analysed for 20 physico-chemical parameters.Differences between sites and seasons were statistically assessed,together with changes in water quality index(WQI).The results indicated that heavy metals(Pb^(2+),Ni^(+),Mn^(2+),Fe^(2+),Cr^(6+)),cations(Ca^(2+),Mg^(2+),K^(+)),total hardness and pH were the main parameters impairing water quality.Drinking water quality standards from both the World Health Organization and Nigeria government were exceeded more often in the wet season than in the dry season.Some groundwater properties were negatively correlated with distance to dumpsite(e.g.,Fe^(2+),Pb^(2+),NO_(3)^(−)).Significant differences between sites were identified,but with no clear spatial trend.WQI varied from excellent(6%-24%of the sites over the study period)to unsuitable for drinking water purposes(12%-18%),with good quality prevailing at most sites(35%-47%).Although groundwater quality declined at 24%of the sites over 2020,the results indicated improvements compared with previous decades.Remediation strategies must be implemented to safeguard public health and the sustainability of water resources.