• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Exploring groundwater quality in semi-arid areas of Algeria: Impacts on potable water supply and agricultural sustainability

    2024-02-26 09:11:04NouaALLAOUAHindaHAFIDHarounCHENCHOUNI
    Journal of Arid Land 2024年2期

    Noua ALLAOUA, Hinda HAFID, Haroun CHENCHOUNI

    1 Laboratory of Natural Resources and Management of Sensitive Environments, Larbi Ben M'hidi University, Oum-El-Bouaghi 04000, Algeria;

    2 Department of Nature and Life Sciences, Faculty of Exact Sciences and Nature and Life Sciences, Larbi Ben M'hidi University,Oum-El-Bouaghi 04000, Algeria;

    3 Laboratory of Algerian Forests and Climate Change, Higher National School of Forests, Khenchela 40000, Algeria

    Abstract: Groundwater quality assessment is important to assure safe and durable water use.In semi-arid areas of Algeria, groundwater represents the main water resource for drinking water supply of the rural population as well as for irrigation of agricultural lands.Groundwater samples from wells and springs were collected from the Gargaat Tarf and Annk Djemel sub-watersheds of the Oum El Bouaghi, Algeria,and were analyzed and compared with the World Health Organization (WHO) standards.Results showed that most of the measured physical and chemical parameters exceeded the quality limits according to the WHO standards.Groundwater had a slightly alkaline water pH (7.00-7.79), electrical conductivity>1500μS/cm, chloride>500 mg/L, calcium>250 mg/L, and magnesium>155 mg/L.Water quality index (WQI)results showed that 68% of the area had excellent water quality, 24% of the samples fell into good category, and only 8% were of poor quality and unsuitable for human consumption.Six wells in the area showed bacterial contamination.Total coliforms (453.9 (±180.3) CFU (colony-forming units)/100 mL),fecal coliforms (243.2 (±99.2) CFU/100 mL), and fecal streptococci (77.9 (±32.0) CFU/100 mL) loads were above the standard limits set by the WHO.These results confirmed that water resources in the study area were strongly influenced by anthropogenic activities and were not recommended for consumption as drinking water.

    Keywords: bacteriological indicator; groundwater; watershed; physical-chemical parameter; water quality index

    1 Introduction

    Groundwater plays a crucial role in ensuring human survival, fostering industrial development,supporting agriculture, and protecting ecosystem.Due to their importance, 97% of continental water are in a liquid state (Yapo, 2010), and these hydrosystems form a real continental ocean(Ntona et al., 2022).The availability and quality of groundwater, which supplies about half of the world's drinking water and more than 40% of its agricultural use, is critical to industrial and agricultural development, and urbanization, as well as food and energy security, and environmental sustainability, also affects human health (Mukherjee et al., 2020; Misstear et al.,2023).

    Groundwater plays a crucial role in the water cycle, establishing vital connections with surface environments, and significantly influencing the functioning of aquatic ecosystems (Chenchouni et al., 2022).Phreatic waters acquire various soluble chemical elements as they circulate underground and come into contact with different geological formations.Additionally, aquifers,such as karst systems, possess specific physical characteristics that can make them more susceptible to surface contamination or the infiltration of contaminants into the subsurface(Misstear et al., 2023).These elements, originating from geological or anthropogenic sources,have the potential to alter groundwater quality and consequently impact ecosystem characteristics(ElKashouty et al., 2022).The quality of groundwater relies not only on the geological composition of the originating soil but also on reactive substances encountered during its flow,ionic exchange, mineral alteration, and the dissolution of authigenic minerals (Etikala et al., 2019;Chenchouni et al., 2023).

    Groundwater, sourced from both confined and surface aquifers, serves as a vital water resource utilized by humans for a multitude of purposes (Prasad and Narayana, 2004).Conversely, the substantial volume of groundwater and its replenishment contribute to the remarkable stability and the dynamics of subsurface ecosystems when compared with the surface systems (Tsujimura et al., 2007).Particularly in arid and semi-arid areas, communities heavily rely on groundwater as a potable water source (Li et al., 2018).

    Among the most water-scarce areas globally, the semi-arid areas of Algeria face significant challenges in terms of water availability.These areas experience reduced precipitation, often below 500 mm and frequently less than 300 mm per year, coupled with high evapotranspiration rate, resulting in scarce surface water resources (Allaoua et al., 2015).In semi-arid aquifers, the quality of groundwater can deteriorate rapidly due to various factors, including over-pumping,posing additional threats (Lar and Gusikit, 2015; Heggy et al., 2022).

    Intensive use of natural resources and increased human activity have created serious problems with groundwater quality (Lapworth et al., 2022).Groundwater is particularly sensitive to different sources of pollution such as uncontrolled wastewater discharges, solid waste disposal,fertilizer, and pesticide applications.These natural hydrosystems are likely to be contaminated by many pollutant sources (Lukubye et al., 2017), such as the results of extensive agro-industrial activities and urbanization in contaminating aquifers including agricultural fertilizers and pesticides, industrial and domestic wastes, dumping in landfills, and pit latrines (Fida et al., 2023).One of the most important types of pollution in rural areas is animal or human fecal contamination (Barnes and Gordon, 2004).Microbial indicators of fecal contamination are employed as markers to assess water quality.This is because fecal matter may harbor pathogenic organisms that can pose potential health risks or adverse effects on human health.Given the significance of scientific research in achieving effective and efficient groundwater quality management, this study aims to review the current state of groundwater quality in the two sub-watersheds: Annk Djemel and Gargaat Tarf, located in the Oum El Bouaghi in Northeast Algeria and understand its impact on drinking water supply and agriculture.Specific objectives include the detailed evaluation of the water physical-chemical characteristics, identification of sources of contamination, and study of the consequences on agriculture from a perspective of sustainability of water resources.The obtained results will be compared with World Health Organization (WHO) standards for drinking water, using the water quality index (WQI) as a means for reporting the data of water quality measurements and provide them in a consistent manner to authorities' management, decision makers, and the public (Qasemi et al., 2023).This comprehensive analysis seeks to provide valuable insights into the water quality status of the mentioned areas, which are primarily influenced by agricultural activities.In addition, the results of this study may inspire discussion about installation and the importance of installation.We will contribute to sustainable groundwater management by introducing a quality control system and establishing a monitoring network to evaluate aquifer degradation.

    2 Materials and methods

    2.1 Study area

    The study area is located at the south of the Oum El Bouaghi, Northeast Algeria and two sub-watersheds (Gargaat Tarf and Annk Djemel) are selected (35°38′-36°00′N, 06°35′-07°20′E;Fig.1).They are separated by the massif of Djebel Tarf that exclusively constituted by sub-horizontal limestone banks of the Cretaceous Age (upper Aptian Stage).The rest of the surface of sub-watersheds is essentially made up of recent geological formations, i.e., Quaternary,formed by massive limestone crusts, polygenic glacis coating the reliefs, ancient silts, arable land,and recent alluvium.Gargaat Tarf is marked by the presence of dune formation or glass.The hydrographic network represented by a poorly developed hair, is of endoreic type.The main wadi,wadi Boulefreis is born in the Aures Area (Allaoua and Hafid, 2019).

    The climate of the area is continental and belongs to the semi-arid bioclimatic stage, with cold winters and hot, dry summers, characterized by irregular precipitation.The prevailing winds are from the southwest, west, and northwest.The minimum precipitation is 107.7 mm and the maximum is 392.0 mm.Seasonal precipitation is 122.0 mm in winter, 108.0 mm in spring, 50.0 mm in summer, and 132.0 mm in autumn.Monthly average temperatures ranged from 6.1°C in December to 38.3°C in August.The average minimum temperature was 2.1°C in December,whereas the warmest month was August with an average of 22.1°C.

    Fig.1 Location and sampling sites in the Annk Djemel and Gargaat Tarf sub-watersheds, Algeria

    2.2 Water sampling and data collection

    The assessment involved the determination and characterization the physical-chemical and bacteriological properties of 25 groundwater samples, collected from wells and springs (Fig.1;Table 1) in 2023.The physical-chemical analyses were carried out on water from private wells and natural springs (Table 1).Because the collection of a water sample is quite delicate, the greatest care must be brought.It conditions the analytical results and the interpretation that will be given.To ensure the samples' homogeneity and representativeness, we collected and stored the water samples following a strict protocol.Polyethylene bottles with a capacity of 500 mL were used for this purpose.Before filling the bottles with water, we thoroughly washed the bottles themselves with nitric acid and then rinsed with distilled water.Additionally, to further minimize any potential contamination, we washed the bottles three times with the same water that was going to be analyzed.Once the cleaning process was completed, the bottles were fully filled with the water samples to be collected.This meticulous procedure aims to maintain the integrity of the water samples and ensure accurate and reliable analyses for the study.

    Table 1 Morphological characteristics of the sample sites in the Annk Djemel and Gargaat Tarf sub-watersheds,Algeria

    2.3 Water physical-chemical analysis

    After conductingin situmeasurements of electrical conductivity (EC), temperature, pH, and dissolved oxygen (DO) using an oximeter CellOX 325 (Xylem Analytics Germany Sales GmbH& Co.KG WTW, Weilheim, Germany) with a WTW multi-parameter device Multiline P3 PH/LF SET (Wissenschaftlich-Technische-Werkst?tten GmbH, Weilheim, Germany), we collected and pre-filtered water samples.The pre-filtered samples were then placed in 1.5 L polyethylene bottles and stored at a temperature of 4°C to preserve their integrity.Chemical analyses were subsequently performed at the Laboratory of the Biology, Department of the Faculty of Exact and Nature Sciences at the Oum El Bouaghi University, Algeria.The determination of calcium was carried out through complexometry using ethylenediaminetetraacetic acid (EDTA) titration in the presence of a colored indicator (Murexide, also known as ammonium purpurate) and a 2-N sodium hydroxide solution.By calculating the difference between total hardness and calcium hardness,we directly obtained the magnesium hardness of the analyzed water.

    This comprehensive analysis and accurate methodology ensure reliable and precise data for the assessment of water quality in the study area (Rodier, 2009).The determination of bicarbonates is carried out by titrimetry, a solution of known concentration of H2SO4was made to act on a precise volume of samples (100 mL) in the presence of methyl orange as an indicator; the reaction that occurs is due to the neutralization of bicarbonate ions by the H+ions of H2SO4.Chlorides are determined volumetrically by Mohr's method, precipitating silver chloride by reaction of chloride ions with silver nitrates and in the presence of 10% K2CrO4as an indicator.SO42-is determined by nephelometry, using a UV-visible spectrophotometer (WTW) set at a wavelength of 420 nm and calibrated prior to any determination of SO42-concentration that is related to the turbidity of the suspension (BaSO4).Nitrates, nitrites, and ammonium were determined by colorimetry.

    2.4 Water bacteriological analysis

    Water samples from wells and springs were carefully collected in sterile bottles to maintain their purity.After collection, we stored and transported these samples to the laboratory in cooler boxes to preserve their integrity and prevent any degradation of water quality during transportation.The laboratory analysis was promptly carried out within the next 4 h to ensure the freshness and reliability of the results.This strict protocol helps to minimize any potential contamination and provides accurate data for the study.Analyzes were performed at the laboratory.

    Bromocresol purple lactose broth (BCP) tubes fitted with a Durham bell were used to count totalEscherichia coliafter a 24-48 h incubation period at 37°C.For confirmatory testing on Schubert's medium fitted with a Durham bell, aliquot positive tubes (lactose fermentation and gas production) and incubate at 44°C for 24-48 h.Addition of 2-3 drops of Kovacs reagent results in the formation of a red ring on the surface, indicating the formation of indole and thus the presence of faecal coliforms.The search for faecal streptococci is carried out on Rothe medium at 37°C for 24 h.From the positive Rothe tubes, a subculture is performed on Litsky medium for 24 h at 37°C.We calculated the results in number of germs per 100 mL according to the Mac-Grady table.

    2.6 WQI

    WQI is a valuable and efficient tool used to evaluate the overall quality of groundwater (Khan and Jhariya, 2017; Adimalla et al., 2018).It is derived from a comprehensive dataset and is presented in a manner that is easily comprehensible for water resource managers and decision-makers (Rana and Ganguly, 2020).WQI offers a consolidated and user-friendly approach to assess groundwater quality, making it a key component in informed decision-making and effective management of water resources.

    WQI is commonly utilized to assess the suitability of groundwater for drinking purposes in various regions worldwide (Khan and Jhariya, 2017; Wu et al., 2017; Adimalla and Venkatayogi,2018; Verma et al., 2018).Five water classes are listed by the WQI value, i.e., excellent (WQI range: 0-25), good (WQI range: 26-50), moderate (WQI range: 51-75), poor (WQI range:76-100), and very poor (WQI>100).To determine if the groundwater is potable, we calculated the WQI using the weighted arithmetic index method:

    whereqnis the quality rating of thenthwater quality parameter;wnis the unit weight for thenthparameter; andnis the total number of parameters included in the WQI calculation.

    The quality index/sub-index (qnth) is a number that reflects the relative value of a certain parameter in polluted waters compared with its standard admissible value:

    wherevnis the estimated value of thenthparameter;snis the standard permissible value of thenthparameter in most cases; andviois the ideal value of thenthparameter.vio=0.0 except for pH and DO.For pH,vio=7.0 andsn=8.5; for DO,vio=14.6 andsn=5.0 mg/L.According to Adimalla and Venkatayogi (2018), we used WQI to determine whether groundwater is suitable for drinking purposes in different parts of the world.Initially, we considered the maximum values for the physical-chemical parameters studied based on the Algerian standards for surface waters (Table 2).These values were utilized to establish the relative importance of each parameter in the overall assessment of water quality.WQI calculation takes into account these weights and constants to provide a comprehensive and standardized evaluation of the water quality, aligning with the established standards for surface waters in Algeria.

    Table 2 Weights of physical-chemical parameters and Algerian standard of surface water quality (OJRA, 2006)

    2.6 Statistical analysis

    In order to compare the values of different physical-chemical parameters in groundwater across the examined sub-watersheds, we computed the means±standard deviations (SD) and the range(minimum-maximum) based on observed data from wells and springs, considering them as replicates within each sub-watershed.To analyze the variation of physical-chemical parameters between the two sub-watersheds, we conducted Welch's two-samplet-test after confirming the normal distribution of data (using the Shapiro-Wilk test) and unequal variances (using theF-test).

    Using Pearson's correlation, the connections between physical and chemical properties of water were examined.The generated correlation matrix was then displayed in a single plot together with the correspondingP-value and correlation coefficient (r).This visualization was accomplished using the Corrplot package (Wei and Simko, 2021).As changes in water characteristics resulting from the growth of one bacterial group may impact the growth of other bacteria either positively or negatively (Loucif et al., 2020), linear regressions and correlation tests were conducted to investigate the interrelationships between the densities of bacterial group.

    To explore the relationships between previously mentioned water parameters and considered bacterial groups, we conducted redundancy analysis (RDA) using the Vegan package in the R software (Oksanen et al., 2022).The RDA biplot was generated using a correlative scaling method.Furthermore, Pearson's correlation tests were performed between the water parameter values observed at each water sample and RDA site scores (weighted sums of response variable scores) on the one hand, and the actual values of meteorological parameters and RDA constraints(linear combinations of climate variables) on the other hand.

    To assess the impact of measured water physical-chemical parameters, represented by the WQI value, on the variation of bacterial loads for each of the three bacteria groups, we employed a generalized linear model (GLM).Bacterial load data were fitted to a Poisson distribution error with a log link function.All statistical analyses in this study were conducted using the R software(R Core Team, 2023).

    3 Results

    3.1 Water physical-chemical property

    The well water's parameters exhibit minimal seasonal variations throughout the year, making it appropriate to characterize each site by calculating the average of measured values during the entire study period for each considered parameter.By taking the average, the data provides a representative and stable representation of the water quality for each site, which enables a comprehensive assessment of the well water conditions without being influenced by short-term fluctuations or seasonal changes.This approach allows for a more reliable and meaningful characterization of water quality at each site for the duration of the study.The values of the physical-chemical parameters of the sampled water are presented in Table 3.The Annk Djemel sub-watershed was characterized by groundwater with a temperature of 16.8°C (±0.3°C).The pH of the water varied between 7.00 and 7.80 with an average of 7.30; as for EC, it varied between 254 and 4160 μS/cm with an average of 2239 μS/cm.DO averaged 6.76 (±0.45) mg/L, regarding calcium carbonates (CaCO3), their concentrations varied between 200 and 1060 mg/L with an average of 531 mg/L, while bicarbonates (HCO3) were 217 (±52) mg/L.The average values of Ca2+were between 109 and 835 mg/L with an average of 383 mg/L, while Mg2+had an average of 147 (±19) mg/L.

    Sulfates ranged from 80 to 472 mg/L and averaged 252 mg/L, chlorides were 591 (±101) mg/L.As for the four nutrients, the average NO3-was 10.49 (±5.13) mg/L, NO2-averaged 0.08 (±0.03)mg/L, as for NH4+was 0.18 (±0.06) mg/L with a range of 0.02-0.70 mg/L and finally for PO43-the values recorded varied between 0.00 and 0.25 mg/L and an average of 0.15 mg/L.

    The water of the Gargaat-Taref sub-watershed marked the following characteristics: the temperature was 16.9°C (±0.26°C) and the pH varied between 6.68 and 8.04 with an average of 7.38.The EC was 3153 (±608) μS/cm, the DO content was 6.79 (±0.32) mg/L.Calcium carbonates averaged 508 (±56) mg/L and bicarbonates ranged from 210 to 833 mg/L with an average of 478 mg/L.Ca2+averaged 446 (±87) mg/L.Mg2+ranged from 32 to 481 with a mean of 163 mg/L.NO3-averaged 4.81 (±1.12) mg/L, NO2-averaged 0.25 (±0.06) mg/L, and NH4+ranged from 0.01 to 0.14 mg/L with an average of 0.25 mg/L.PO43-ranged from 0.00 to 0.29 mg/L and an average of 0.11 mg/L, Cl-averaged 725 (±188) mg/L and SO42-ranged from 74 to 756 mg/ and an average of 452 (±76) mg/L.Despite the spatial variations observed in the physical-chemical parameters of the sub-watersheds, the difference was statistically insignificant(P>0.050), except for bicarbonates (P<0.003), sulfates (P=0.028), and nitrites (P=0.035) (Table 3).

    3.2 Relationships between water parameters

    Correlation tests between the groundwater physical-chemical parameters of the two subwatersheds are summarized in Figure 2.Correlation analysis in the Annk Djemel and Garrat Taref sub-watersheds resulted in significant correlation in 19 of 98 pairs of water attributes.For the Annk Djemel sub-watershed, the highest negative correlations were obtained between the following pairs: DO and temperature (r= -0.94,P<0.001), HCO3and temperature (r= -0.79,P=0.001), HCO3-and EC (r= -0.63,P=0.031), NH4+and DO (r= -0.79,P=0.003), NO2and temperature (r= -0.66,P=0.014), and NO2-and EC (r= -0.63,P=0.008).The highest positive correlations were obtained for Ca2+and CaCO3(r=0.97,P<0.001), and HCO3-and NO2-(r=0.92,P<0.001).Other positive correlations were recorded between the following pairs: HCO3-and DO(r=0.83,P=0.000), CaCO3and EC (r=0.77,P=0.002), Ca2+and EC (r=0.77,P=0.004), Mg2+and CaCO3(r=0.77,P=0.002), Mg2+and EC (r=0.60,P=0.030), Cl-and CaCO3(r=0.78,P=0.002),Cl-and Ca2+(r=0.77,P=0.002) and NH4+and T (r=0.69,P=0.009).On the other hand, a significant correlation in 25 of 98 pairs of water attributes was found in the Gargaat Taref sub-watershed.There was a strong correlation between EC, CaCO3(r=0.91,P<0.001), sulfates(r=0.89,P<0.001), and Cl-(r=0.85,P=0.001).Other significant correlations were noted between EC and CaCO3(r=0.71,P=0.010).A significant correlation was also found between SO42-and HCO3-(r=0.80,P=0.002), SO42-and NO2-(r=0.76,P=0.004), SO42-and NO3-(r=0.73,P=0.007),SO42-and Cl-(r=0.72,P=0.008), NO3-and NO2-(r=0.86,P<0.001), NO3-and HCO3-(r=0.83,P=0.001), NO2-and HCO3-(r=0.81,P=0.002), and PO43-and NO3-(r=0.60,P=0.041).On the other hand, negative correlations included the following pairs: OD and T (r= -0.91,P<0.001), OD and NH4+(r= -0.71,P=0.010), pH and HCO3-(r= -0.71,P=0.009), pH and NO3-(r= -0.61,P=0.035), pH and PO43-(r=071,P=0.010), and pH and Cl-(r=0.59,P=0.044).

    Table 3 Physical-chemical characteristics of groundwater in the Gargaat Tarf and Annk Djemel sub-watersheds,Algeria

    Fig.2 Correlation between water physical-chemical parameters from wells and groundwater in the Annk Djemel (a) and Gargaat Tarf (b) sub-watersheds, Algeria.***, P<0.001 level; **, P<0.010 level; *, P<0.050 level.T,temperature; EC, electrical conductivity; DO, dissolved oxygen.

    3.3 Bacteriological quality of groundwater

    Variations of bacterial load in the groundwater of the Annk Djemel sub-watershed are depicted in Figure 3.Bacteriological examination of the water revealed that, out of the 13 wells analyzed, 6 wells contain a total heterotrophic bacterial load.The average concentration of total coliforms for all wells during the study period was approximately 453.9 (±180.3) CFU/100 mL, with a maximum value of 1900.0 CFU (colony-forming units)/100 mL observed in the well P9.As for fecal coliforms, their average concentration in all wells was approximately 243.2 (±99.2)CFU/100 mL, with a maximum of 1000.0 CFU/100 mL found in the well P3.As for fecal streptococci, we noted that the average concentration in well water was 77.9 (±32.0) CFU/100 mL with a maximum of 324.0 CFU/100 mL in the well P3.On the relative variation, the result of the coefficient of variation (CV) showed that all the concentrations of bacterial load of the examined surface waters were heterogeneous (CV=143% in total coliforms, CV=147% in fecal coliforms,and CV=151% in fecal streptococci).

    Fig.3 Chord diagram displaying the distribution of densities (CFU/100 mL) of main bacterial groups.TC, total coliforms; FC, fecal coliforms; FS, fecal streptococci.P1-P3 and P7-P9 are sampling sites where bacterial groups are positively detected.

    3.4 Interrelationships between bacterial groups

    Interrelationship between bacterial groups is shown in Figure 4.All bacterial densities are correlated with each other.Positive significant correlations (P<0.001) were revealed between FC and fecal coliforms (r=0.96), fecal coliforms and total coliforms (r=0.92), fecal streptococci and total coliforms (r=0.89).The slopes of linear of regression equations for the correlations ranged from 0.161 to 2.927 being the largest for fecal coliforms.

    Fig.4 Correlative scatter plot among bacterial groups detected at groundwater in semi-arid areas of Algeria.Red lines represent linear regressions with the equations given at top of plots above the diagonal.Blue ellipses represent 50% and 95% concentration levels of observations with the center in solid red circle.(a1 and a2),correlations between FS (fecal streptococci) and FC (fecal coliforms), and between TC (total coliforms) and FS;(b1 and b2), correlations between FC and FS, and between TC and FS; (c1 and c2), correlations between FC and TC, and between FS and TC.

    3.5 Effects of water characteristics on bacterial loads

    RDA was employed to examine the relationships between groundwater physical-chemical parameters and bacterial load.The results indicated that the first and second axes explained 96.9% and 2.9% of the eigenvalues, respectively.According to the first axis of RDA, which represents the majority of information, we found that most of the water parameters exhibited negative correlations with pH, DO, HCO3-, and NO2-.These parameters were positioned on the positive side of the axis.Moreover, the distribution of bacterial groups was also plotted on the positive side of the axis (Fig.5; Table 4).This analysis suggests a clear association between the physical-chemical parameters of groundwater and the presence of bacterial groups, indicating potential inter-relationships between these variables.

    Fig.5 Tri-plot of the RDA (redundancy analysis) examining relationships concentrations of water physical-chemical parameters and loads of bacterial groups in groundwater in semi-arid areas of Algeria.EC, electrical conductivity;DO, dissolved oxygen.

    Table 4 Correlations of groundwater physical-chemical parameters within the first three axes of RDA(redundancy analysis)

    On the second axis of RDA that presents only 2.9% of the variation, the bacterial load of total coliforms was positively associated with the concentrations of NO3-, NH4+, SO42-, and the temperature in the water, but negatively with the pH of the water.DO, Mg2+, HCO3-, Cl-, and EC all along the axis were positively correlated with fecal coliforms and to a lesser degree with fecal streptococci.RDA triplot indicated that all the abiotic parameters and the bacterial load were negatively correlated on the first axis, indicating that bacterial load increased according to the increase of a few parameters such as NO3-, temperature, and NH4+; but were not statistically significant (P>0.050).According to the parameters of water and bacterial abundance, there was a clear spatial separation, which allowed us to classify our sample sites into two groups according to the physical-chemical quality and bacterial load: the first group formed of the wells P9, P1, and P3.The water of these wells showed significant mineralization and pollution by NH4+and NO3-as well as a high bacterial load.The first group was opposed to the second group, which included the rest of the wells (P2, P4, P5, P6, P7, P10, S1, S2, and S3), and the three sources characterized by a relatively average mineralization compared with the other sites and almost zero pathogen contamination.

    3.6 Effect of WQI on bacterial loads of contamination-indicators

    WQI of the Gargaat Tarf varied from 17% to 50%, indicating that the water quality varies from excellent to very poor, while WQI of the Annk Djemel indicated that 85% of all samples analyzed were of excellent quality, 15% were good.For the entire study area, we found that 68% of samples fell into the excellent category and 24% of samples fell into good categories and only 8%of poor quality and unfit for human consumption (Fig.6).

    Throughout the study period, there was a significant difference between average value of WQI(P=0.032) in the two sub-watersheds.On the relative variation, the result of CV observed for the two sub-watersheds was largely higher, indicating that WQI was heterogeneous.In comparison,the basin of the Ank Djemel with the CV of 139% had a more heterogeneous water quality than the Gargaat Tarf (CV=87%).Regarding the results of GLM, the effect of WQI on bacterial load showed a negative effect on the three groups of bacteria indicators: fecal coliforms (t= -0.97),fecal streptococci (t= -0.93), and total coliforms (t= -1.00).Although these relationships were negative (Fig.7), GLM analysis revealed that there were no significant effects observed for the bacterial groups.Similarly, when testing the effects of water physical-chemical parameters on the variation of bacterial group loads, GLM also demonstrated no significant effects (Table 5).These results suggest that the variations in bacterial loads were not significantly influenced by the physical-chemical parameters of the water, indicating a lack of significant correlation between these factors.

    Fig.6 Classes of water quality index of groundwater in the Annk Djemel and Gargaat Tarf sub-watersheds, Algeria

    Fig.7 GLM (generalized linear model) plot showing the effect of WQI (water quality index) on the variation of water bacterial groups for groundwater at the Annk Djemel sub-watershed, Algeria.FC, fecal coliforms; FS, fecal streptococci; TC, total coliforms.Solid lines represent linear regressions obtained by a GLM fit with 95%confidence level in light grey.

    Table 5 GLM (generalized linear model) testing the effects of water physical-chemical parameters on the variation of bacterial groups

    4 Discussion

    The assessment of water quality is generally based on the monitoring of microbial elements, in particular the presence of fecal coliform bacteria, as well as on the analysis of physical-chemical characteristics (Guemmaz et al., 2020; Loucif et al., 2020).Water parameters and quality are influenced by a complex interaction between external and internal factors.External factors include weather conditions, substrate characteristics (such as soil and sediment), and pollution sources, while internal factors are the result of biochemical reactions occurring in the water itself(Hacioglu and Dulger, 2009).

    4.1 Water physical-chemical property

    The results of physical and chemical parameters of groundwater were compared with the standard guideline values recommended by the WHO (WHO, 2017).Temperature is one of the physical-chemical parameters used to assess the quality of drinking water.It plays a role in many phenomena, such as the speed of chemical reactions in the body of water, the decrease in the solubility of gases, as well as the accentuation of the tastes and colors of water (Daghara et al.,2019; Calero Preciado et al., 2021; Moussaoui et al., 2024).

    The average temperature of the water in all seasons was between 15°C and 19°C; the thermal variation of the water is about 0.2°C.Groundwater temperature can vary depending on the depth of aquifer and seasons.For example, the European Parliament has established an indicative value for the temperature of water for human consumption at 25°C, considering this temperature as a limit not to be exceeded (Rodier, 2009).

    In the study area, pH values ranged from 6.68 to 8.04, with an average of 7.36.Most groundwater samples were between weakly alkaline and moderately alkaline levels, this phenomenon being attributed to the flow of water through the surface and subsurface carbonate layers (Rao et al., 2012).Slightly alkaline water has the ability to inhibit heavy metal toxicity in the form of precipitated carbonate or bicarbonate, rendering these heavy metals unavailable(Ahipathy and Puttaiah, 2006).The conductivity of water is an essential parameter because it evaluates the capacity of an aqueous medium to allow the passage of an electric current (Sui et al.,2018).Groundwater EC generally ranges from 254 to 7160 μS/cm, and this is closely related to the concentrations of ions present in these waters.These higher values are likely attributable to a high concentration of ions.Therefore, they indicate higher levels of total dissolved solids in groundwater, which contributes to increased salinity (Abboud, 2018).Furthermore, the large variations in EC values reflect significant differences in the predominant geochemical approaches in the study area.

    DO concentration ranged from 3.99 to 9.53 mg/L, with an average of 6.77 mg/L.According to the standards, the maximum allowable limit for DO in drinking water is 8.00 mg/L (WHO, 2017).Concentration of DO plays a pivotal role in groundwater quality as it regulates the valence state of trace metals and restricts the bacterial metabolism of dissolved organic species.Due to these factors, measuring the DO concentration is considered vital in the majority of water quality studies (Bouaroudj et al., 2019; Cheng et al., 2019; Zhang et al., 2022).

    Bicarbonates trigger the groundwater alkalinity (Adams et al., 2001), these elements are found naturally in waters, and their presence is attributed to the dissolution of carbonate formations such as limestone and dolomites (Drever, 2001).Bicarbonate is the most common anion in groundwater and comes mainly from CO2present in the soil (Ahuja et al., 2008).Its concentration varied from 55 to 833 mg/L in water, and its origin can be attributed to the leaching of carbonate formations.Singh et al.(2012) reported that during irrigation, rainfall infiltration, and groundwater movement, carbonates existing in carbonate rocks may have been dissolved and hence supplied to the groundwater system along with recharge water.These processes could account for the low alkalinity levels seen in some of the examined wells, possibly reflecting the hydrochemical immaturity of the groundwater or the shallowness of aquifer (Demetriades, 2011).

    Cl-is found naturally in all types of water, but their main contributing sources are runoff,inorganic fertilizers used in agricultural fields, and wastewater discharges (?adraku, 2021).In the present study, Cl-varied between 64 and 2482 mg/L, with an average of 655 mg/L, which exceeds the permissible limit (250 mg/L).The presence of high levels of Cl-in groundwater makes it potentially hazardous to human health (Pius et al., 2012; Noori et al., 2014).SO42-come from the dissolution and leaching of rocks containing gypsum, iron sulphides, and other sulfur compounds (Arjum et al., 2021).In 75% of groundwater samples, the concentration of SO42-exceeds the WHO recommended limit, which is set at 200 mg/L for direct drinking and domestic use.

    Based on WHO guidelines (WHO, 2017), the concentration of Ca2+in drinking water should not exceed 200 mg/L.However, this limit is often exceeded in groundwater in many parts of the world, including Europe, Asia, and Africa.Ca2+concentrations varied between 109 and 1239 mg/L, with an average value of 250 mg/L.Ca2+originates from the dissolution of carbonate and evaporitic minerals such as calcite, dolomite, aragonite, gypsum, and anhydrite, as well as from the carbonate cement present in geological formations (Bozda? and G??mez, 2013).Mg2+varied between 32 and 481 mg/L, with an average of 155 mg/L.It is very likely that the presence of Ca2+and Mg2+in groundwater comes from the leaching of carbonate minerals such as calcite and dolomite (Magesh et al., 2013).Geospatial distribution of cations showed that concentrations are higher in the Gagraat Taref than the Annk Djemel.It is very likely that the major source of Mg2+in groundwater is from magnesium-containing minerals such as dolomite and sulfate minerals present in the study area (?ener et al., 2017).

    The presence of NO3-in groundwater can be attributed to fertilizer application, human and animal sewage, plant deposits, and other nitrate-rich wastes (Chukwura et al., 2015; Vincy et al.,2015).The concentration of NO3-in the analyzed water samples ranged from 0.01 to 68.22 mg/L,which was below the WHO permissible consumption limits (50.00 mg/L).NH4+, NO2-, and PO43-were well below the maximum desirable limits of the WHO.Their mean values are 0.22, 0.16,and 0.13 mg/L, respectively.The low concentrations of the three elements are attributed to the low contamination of these waters with pollution sources.

    The use of Pearson's correlation matrix analysis is a valuable tool that can reveal the links and associations between hydrogeochemical parameters, allowing their origins and interrelationships to be identified (Guey-Shin et al., 2011; Mgbenu and Egbueri, 2019).According to their report,when the correlation coefficient (r) is greater than 0.70, it indicates a strong correlation between the two parameters, while if the value ofris between 0.50 and 0.70, it suggests a moderate correlation.A good correlation was observed between EC and HCO3-(r=0.63), EC and CaCO3(r=0.77), Ca2+and CaCO3(r=0.77), Mg2+and CaCO3(r=0.78), Ca2+and SO42-(r=0.91), Ca2+and Cl-(r=0.85), and SO42-and HCO3-(r=0.80) indicating that all of these elements could have the same origin.According to Singh and Mukherjee (2015), the increased presence of these ions can be attributed to the dissolution of evaporitic minerals.In addition, the ionic strength of water increases due to the dissolution of these evaporitic minerals, which also promotes the dissolution of sulphate salts, thus leading to an increase in the concentration of Mg2+and Ca2+in groundwater(Alaya et al., 2014).Negative correlations were observed between temperature and DO (r= -0.94)as well as between DO and ammonium (r= -0.71), illustrating the negative interactions between these parameters.Indeed, temperature influences both the solubility of oxygen in water and biological activities such as photosynthesis and respiration, which impacts the short-term oxygen concentration balance (Butcher and Covington, 1995; Daghara et al., 2019; Calero Preciado et al.,2021).

    4.2 Bacteriological quality of groundwater

    Analyzing the abundance and diversity of microorganisms is a useful approach for assessing the quality of drinking water (Favere et al., 2021).However, it is difficult to test the water for all of these organisms, as the isolation and identification of them is extremely complex (Soni et al.,2023).A commonly used indirect method to assess fecal contamination is to analyze coliform bacteria, which are widely considered reliable indicators of such contamination (Kistemann et al.,2002; Nola et al., 2002; Barnes and Gordon, 2004).The coliform bacteria concentrations were 453.9 CFU/100 mL for TC, 243.2 CFU/100 mL for FC, and 77.9 CFU/100 mL for FS in the wells P9, P1, and P3, respectively, due to the application of animal manure or dung.The coliform bacteria concentration was lower in the deepest wells, thus confirming that an artesian well is safer than a surface well or an outcropping spring.Several studies around the world have corroborated these observations, confirming that runoff and short residence times in permeable soils play a crucial role in the direct infiltration of manure into groundwater (Richardson et al.,2009; Masoud et al., 2016).Thus, the practice of open defecation and the use of private septic tanks lead to sewage entering groundwater sources in rural areas around the world (Schijven et al.,2010; Masoud et al., 2016).In the same context, Joseph et al.(2021) also reported that livestock grazing activity elevated bacteria numbers.Whereas Chadwick and Chen (2002) reported that the highest number of fecal coliforms and fecal streptococci were still present in manure.

    The lack of protection of most wells is also a factor that contributes to the enrichment of water in organic matter due to windblown inputs of leaves and other plant debris that will decompose on site.All indicator bacteria were highly correlated with each other, the correlation coefficientrvaries from 0.89 to 0.96 (P<0.001), these positive correlations that were recorded show that these indicators probably have the same source of contamination.Atherholt et al.(2003) suggest that many total coliforms in the environment can only come from fecal sources.According to Francy et al.(2000), a strong correlation between two indicators may indicate that these two indicators come from identical or similar sources of contamination.Atherholt et al.(2003) reported that the existence of both fecal coliforms and fecal streptococci indicators provide strong evidence of fecal contamination.

    4.3 Effects of water characteristics on bacterial loads

    The growth of microorganisms in water is closely influenced by environmental factors (Fister,2016), which establishes a relationship between physical-chemical properties and the amount of bacteria in this study, bacteriological analysis revealed different groups of bacteria with similar tendencies.Similar results were reported by Guemmaz et al.(2020), showing the influence of environmental factors on fecal bacterial loads measured in urban effluents discharged into the arid wadis of Algeria.Variations in environmental factors such as pH, salinity, metal concentrations,and energy exert selective and direct effects on the composition, abundance, and spatial distribution of microbial communities (Dell'Anno et al., 2003; Franklin et al., 2007; Guo et al.,2015; Jordaan et al., 2016; Liu et al., 2018).RDA analysis shows that NO3-, NH4+, SO42-,temperature, Ca2+, EC, Cl-, CaCO3, and Mg2+are the environmental variables that statistically best explain the variations in bacterial distribution.Total coliforms include fecal coliforms and all bacteria in this group have similar biochemical characteristics and are found in soil or water(Rodier et al., 2009).

    In the present study, total coliforms responded positively to increased nitrate.Vandenberg et al.(2005) reported that coliform bacteria use ammonia and nitrate as nutrients.Fecal coliforms are specifically associated with feces (Rodier et al., 2009), and nitrate ions can be produced by the oxidation of ammonium ions present in feces (Bou Saab et al., 2007).Nitrate ions and fecal coliforms therefore share a common origin.Indeed, water temperature is one of the most significant environmental parameters that influence the concentrations of fecal indicator bacteria.Indeed, some studies have reported an inverse relationship between fecal coliforms and water temperature, while other researchers have found a direct relationship between fecal coliforms and water temperature (Loucif et al., 2020; Valenca et al., 2022).In our study, we found that increasing temperature promotes the development of bacterial load of total coliforms, fecal coliforms, and fecal streptococci.These results are consistent with that reported by Chigbu et al.(2005), who also found a correlation between temperature and fecal coliforms.However, our results differ from those reported by Tek et al.(2001), who did not observe a relationship between temperature and these bacterial indicators.

    Variable responses of bacterial populations were also observed in relation to other physical-chemical parameters of groundwater.For example, high values of EC, Mg2+, Ca2+, and CaCO3have also been associated with an increase in the abundance of these bacteria (Nola et al.,2002).In a bacteriological and chemical investigation of groundwater in Cameroon, Nola et al.(2002) discovered that higher EC and concentrations of oxygen, Cl-, Na+, K+, Ca2+, and Mg2+promoted the abundance of fecal coliforms and fecal streptococci.Although Ondieki et al.(2021)claim that total coliforms rise at acidic pH, Loucif et al.(2020) assert that alkaline pH obviously impairs the survival of fecal coliforms.Depending on the environmental parameters available,such as the ideal pH, temperature, amount of nutrients, and suspended particles, the coliform bacteria's survival can always be prolonged, in some cases, and can even be expanded (Juhna et al.,2007).

    4.4 WQI

    A total of 68% of groundwater samples taken during the research period classed as excellent water,with good water about 24%, whereas poor water was 8%.The latter, which exhibits poor quality,may be a result of agricultural influence, effective ion leaching, excessive groundwater use, direct effluent discharge, and other factors (Sahu and Sikdar, 2008).Indeed, WQI is an essential parameter for evaluating the quality of groundwater and determining its suitability for consumption.In reality, WQI is an essential indicator for determining groundwater's quality and establishing whether it is fit for human use (Mishra and Patel, 2001; Naik and Purohit, 2001;Avvannavar and Shrihari, 2008; Rana and Ganguly, 2020; Ram et al., 2021).

    WQI was tested with bacterial load values, which are chosen as indicators of pollution with the GLM test.Thet-test shows negative correlations between WQI and bacterial load.Thus, low values of WQI were found to be closely related to high load values of total coliforms, fecal coliforms, and fecal streptococci, which indicates poor quality for some groundwater samples in the study area; on the other hand, GLM indicated a non-significant effect of this correlation.On the relative variation, the results of CV showed that all the values of WQI were heterogeneous,thus signifying the spatial diversity of study area as well as various pollutants that threaten the groundwater quality.

    This study offers a robust assessment of groundwater quality in semi-arid areas of Algeria,providing valuable insights into the suitability of water for both drinking and agricultural purposes.The comprehensive analysis includes a range of physical-chemical parameters, and the application of WQI enhances the clarity of the findings.However, it is essential to acknowledge the limitations of this study.The sample size is relatively small, potentially limiting the generalizability of the findings.Furthermore, the focus on specific parameters and the absence of consideration for seasonal variations might not fully encapsulate the dynamic nature of groundwater quality.Future research endeavors with larger and more diverse datasets, along with an exploration of temporal variations, would strengthen the overall reliability and applicability of the results.

    4.5 Implementation and perspective

    In the larger context, this study prompts a reevaluation of water management practices,advocating for policies that prioritize the preservation and enhancement of groundwater quality.The result emphasizes the urgency of collaborative efforts involving stakeholders, policymakers,and local communities to ensure access to safe and potable water-a fundamental right for the present and a necessity for future generation.In the face of bacterial contamination and water-related health risks, diverse solutions can be considered.Firstly, improving sanitation infrastructure is crucial to prevent contaminant spills into groundwater sources.Additionally,reinforcing water point protection through filters and physical barriers can mitigate direct contamination.Raising community awareness about proper hygiene practices and preserving water sources is equally imperative.Implementing regular water quality monitoring is essential to swiftly respond to bacterial contamination.Simultaneously, strict adherence to water quality standards and tailored regulation is necessary to oversee potentially contaminating activities.Lastly, fostering multi-stakeholder collaboration among local authorities, government organizations, communities, and water management experts is paramount to coordinate efforts aimed at sustainably restoring and preserving water quality and this concerted action ensures a safe supply of potable water for local populations.

    5 Conclusions

    The assessment of groundwater quality emerges as a very important aspect in ensuring the safety and sustainability of water usage.In semi-arid areas like Algeria, where groundwater serves as the primary water source for rural communities and agricultural activities, its quality holds immense significance.Our study extensively analyzed various parameters in groundwater samples from the Gargaat Tarf and Annk Djemel sub-watersheds of Algeria.The findings showed that physical and chemical parameters surpassed the recommended limits by the WHO.Notably, pH ranged slightly alkaline (7.00-7.79), with elevated levels of EC, Cl-, Ca2+, Mg2+, and other substances, signaling compromised water quality.Despite these concerns, WQI offered a nuanced perspective.A significant portion of the area, approximately 68%, exhibited excellent water quality, while 24%fell under the good category.However, alarmingly, 8% indicated poor quality, rendering them unsuitable for human consumption, emphasizing the imperative need for remediation and management strategies.The presence of bacterial contamination, particularly in six wells, further underscores the influence of human activities on the water resources.Elevated levels of coliforms and streptococci above the WHO standards reflect the anthropogenic impact on water quality,posing health risks, and deeming certain water points unsuitable for drinking.These revelations call for immediate attention and strategic interventions.The significance of this research extends beyond mere observation; it demands actionable steps to rectify the deteriorating quality of groundwater.Efforts must focus on targeted measures to mitigate anthropogenic influences,implement rigorous monitoring frameworks, and devise interventions that safeguard water resources.

    Conflicts of interestsThe authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgements

    Haroun CHENCHOUNI received funding from the European Union's Horizon 2020 Research & Innovation Programme (2211) under the Partnership for Research and Innovation in the Mediterranean Area (PRIMA) Project''SHARInG-MeD'' and from the Directorate-General for Scientific Research and Technological Development(DGRSDT) under the Projets de Recherche Formation-Universitaire (PRFU) Projects (D00L02UN120120230002,D01N01UN120120230005).We express our gratitude to Dr.Haythem DINAR (University of Oum El Bouaghi,Algeria) for his invaluable assistance in plotting the Figure 1.

    Author contributions

    Conceptualization, methodology, investigation, data curation, resources, and writing - original draft preparation:Noua ALLAOUA; writing - original draft and writing - reviewing and editing: Hinda HAFID; Formal analysis,visualization, writing - original draft, and writing - reviewing and editing: Haroun CHENCHOUNI.All authors approved the manuscript.

    亚洲欧美一区二区三区久久| 啦啦啦啦在线视频资源| 国产一区亚洲一区在线观看| 婷婷色麻豆天堂久久| 精品国产乱码久久久久久男人| 80岁老熟妇乱子伦牲交| 老汉色∧v一级毛片| 欧美日韩精品网址| 午夜91福利影院| 国产成人a∨麻豆精品| 免费观看无遮挡的男女| 一级毛片电影观看| √禁漫天堂资源中文www| 久久久久视频综合| 国产欧美日韩综合在线一区二区| 免费观看a级毛片全部| 亚洲天堂av无毛| 久久久久久久久久久久大奶| 中文字幕人妻丝袜制服| av.在线天堂| 天堂中文最新版在线下载| 99精国产麻豆久久婷婷| 亚洲男人天堂网一区| 秋霞伦理黄片| 精品一区二区三卡| 自线自在国产av| 熟妇人妻不卡中文字幕| 99热全是精品| 日韩中文字幕欧美一区二区 | 欧美日韩综合久久久久久| 久久久久精品人妻al黑| 日本黄色日本黄色录像| 日本黄色日本黄色录像| 亚洲精品美女久久av网站| 国产毛片在线视频| 99re6热这里在线精品视频| 国产成人一区二区在线| 高清视频免费观看一区二区| 亚洲第一青青草原| 欧美日本中文国产一区发布| 久久久久精品人妻al黑| 中文字幕人妻丝袜一区二区 | 18+在线观看网站| 熟女av电影| 亚洲男人天堂网一区| 欧美人与善性xxx| www.av在线官网国产| 精品人妻熟女毛片av久久网站| 国产午夜精品一二区理论片| 久久av网站| 国产探花极品一区二区| av有码第一页| 国产白丝娇喘喷水9色精品| 岛国毛片在线播放| 日本黄色日本黄色录像| 欧美激情高清一区二区三区 | 久久久国产精品麻豆| 国产在线免费精品| 国产不卡av网站在线观看| 国产不卡av网站在线观看| 高清视频免费观看一区二区| 女性生殖器流出的白浆| 肉色欧美久久久久久久蜜桃| 免费女性裸体啪啪无遮挡网站| av.在线天堂| 亚洲av欧美aⅴ国产| 亚洲一码二码三码区别大吗| 高清黄色对白视频在线免费看| 丝袜人妻中文字幕| 69精品国产乱码久久久| 婷婷色综合大香蕉| 超色免费av| 久久久亚洲精品成人影院| 大话2 男鬼变身卡| 欧美精品人与动牲交sv欧美| 午夜日本视频在线| freevideosex欧美| 观看av在线不卡| 一级毛片 在线播放| 一本—道久久a久久精品蜜桃钙片| 一本—道久久a久久精品蜜桃钙片| 日韩av免费高清视频| 欧美精品国产亚洲| 国产男女内射视频| 多毛熟女@视频| 亚洲精品久久午夜乱码| 人妻人人澡人人爽人人| 国产老妇伦熟女老妇高清| 欧美 日韩 精品 国产| 伊人久久大香线蕉亚洲五| 国产深夜福利视频在线观看| 精品国产一区二区三区四区第35| 少妇精品久久久久久久| 秋霞伦理黄片| 人人妻人人爽人人添夜夜欢视频| 亚洲精品美女久久av网站| tube8黄色片| 亚洲伊人色综图| 中文字幕亚洲精品专区| 99re6热这里在线精品视频| 亚洲精品久久午夜乱码| 18在线观看网站| 精品人妻一区二区三区麻豆| 国产极品粉嫩免费观看在线| 亚洲精品美女久久久久99蜜臀 | 久久亚洲国产成人精品v| 日本欧美国产在线视频| 国产欧美日韩综合在线一区二区| 日韩在线高清观看一区二区三区| 99久国产av精品国产电影| 久久久久国产一级毛片高清牌| 成年人免费黄色播放视频| 国产精品蜜桃在线观看| 国产1区2区3区精品| 精品国产一区二区三区四区第35| 国产成人aa在线观看| 多毛熟女@视频| 你懂的网址亚洲精品在线观看| 日日啪夜夜爽| 国产 精品1| 制服人妻中文乱码| 精品一区二区三区四区五区乱码 | 女人久久www免费人成看片| 三上悠亚av全集在线观看| 亚洲欧美一区二区三区久久| 爱豆传媒免费全集在线观看| 男女啪啪激烈高潮av片| 国产探花极品一区二区| 日日啪夜夜爽| 成年av动漫网址| 不卡av一区二区三区| 久久影院123| 精品久久久久久电影网| 最近最新中文字幕大全免费视频 | 只有这里有精品99| 精品国产露脸久久av麻豆| 国产一区二区三区综合在线观看| av视频免费观看在线观看| 91成人精品电影| 欧美成人精品欧美一级黄| 久久久久国产精品人妻一区二区| 一本久久精品| 欧美日韩国产mv在线观看视频| 两个人看的免费小视频| 国产男女内射视频| 日本vs欧美在线观看视频| 2022亚洲国产成人精品| 亚洲国产精品一区二区三区在线| 捣出白浆h1v1| 一区二区日韩欧美中文字幕| 伊人久久大香线蕉亚洲五| 国产精品免费视频内射| 中文字幕精品免费在线观看视频| 久久99热这里只频精品6学生| 搡老乐熟女国产| 99久久中文字幕三级久久日本| 在线观看免费日韩欧美大片| 免费女性裸体啪啪无遮挡网站| 日本黄色日本黄色录像| 成人午夜精彩视频在线观看| 91成人精品电影| 老汉色∧v一级毛片| 国产精品99久久99久久久不卡 | 中文精品一卡2卡3卡4更新| 一区二区av电影网| 人人妻人人澡人人看| 久久99精品国语久久久| 男的添女的下面高潮视频| 美女脱内裤让男人舔精品视频| 国产在线免费精品| 国产毛片在线视频| 欧美激情 高清一区二区三区| 久久精品熟女亚洲av麻豆精品| 三上悠亚av全集在线观看| 久久精品国产a三级三级三级| 人妻少妇偷人精品九色| av不卡在线播放| 2021少妇久久久久久久久久久| 亚洲精品av麻豆狂野| 91久久精品国产一区二区三区| 久久久国产一区二区| 美女国产视频在线观看| 国产一区二区三区av在线| 国产成人精品久久二区二区91 | 欧美成人午夜精品| 日本av手机在线免费观看| 一边亲一边摸免费视频| 啦啦啦在线观看免费高清www| 18禁裸乳无遮挡动漫免费视频| 又大又黄又爽视频免费| 久久精品久久久久久噜噜老黄| 波多野结衣一区麻豆| 欧美日韩精品成人综合77777| 亚洲精品自拍成人| 叶爱在线成人免费视频播放| 黑丝袜美女国产一区| 亚洲经典国产精华液单| 亚洲中文av在线| 性高湖久久久久久久久免费观看| 国产一区二区三区综合在线观看| 久久久久国产网址| 大香蕉久久网| 久久久久久久久久久免费av| 大香蕉久久网| 国产高清不卡午夜福利| av有码第一页| 在线天堂最新版资源| 欧美最新免费一区二区三区| 国产成人精品在线电影| 亚洲,欧美,日韩| 国产探花极品一区二区| 女人久久www免费人成看片| 九草在线视频观看| 国产一区二区激情短视频 | 亚洲三级黄色毛片| 日韩精品免费视频一区二区三区| 精品亚洲成a人片在线观看| 在线观看一区二区三区激情| 少妇精品久久久久久久| 日韩电影二区| 欧美97在线视频| 日本欧美视频一区| 天天影视国产精品| 青春草视频在线免费观看| 午夜福利一区二区在线看| 国产精品香港三级国产av潘金莲 | 男人操女人黄网站| 日本wwww免费看| 亚洲美女视频黄频| 色吧在线观看| 激情五月婷婷亚洲| 久久午夜福利片| 日韩一区二区三区影片| 少妇精品久久久久久久| 国产精品.久久久| 汤姆久久久久久久影院中文字幕| 国产97色在线日韩免费| 黄色怎么调成土黄色| 亚洲三区欧美一区| 一区福利在线观看| 你懂的网址亚洲精品在线观看| 久久精品久久久久久噜噜老黄| 久久人人爽av亚洲精品天堂| 不卡av一区二区三区| av线在线观看网站| 十八禁网站网址无遮挡| 日韩一本色道免费dvd| 亚洲综合色惰| 免费av中文字幕在线| 国语对白做爰xxxⅹ性视频网站| 妹子高潮喷水视频| 久久亚洲国产成人精品v| 男女下面插进去视频免费观看| 久久久久久伊人网av| 亚洲情色 制服丝袜| 下体分泌物呈黄色| 男女午夜视频在线观看| 色94色欧美一区二区| 观看av在线不卡| 久久久久人妻精品一区果冻| 国产精品熟女久久久久浪| 亚洲少妇的诱惑av| 国产精品二区激情视频| 伊人亚洲综合成人网| 一边亲一边摸免费视频| 性色avwww在线观看| 九九爱精品视频在线观看| 久久鲁丝午夜福利片| 两性夫妻黄色片| 9色porny在线观看| 最近的中文字幕免费完整| 国产精品香港三级国产av潘金莲 | 人妻系列 视频| 波多野结衣av一区二区av| www.av在线官网国产| 超色免费av| 久久ye,这里只有精品| 老司机影院成人| 国产乱人偷精品视频| 亚洲国产日韩一区二区| 国产精品成人在线| 国产高清国产精品国产三级| 777久久人妻少妇嫩草av网站| 天天躁狠狠躁夜夜躁狠狠躁| 免费大片黄手机在线观看| 国产成人午夜福利电影在线观看| 国产精品av久久久久免费| 伦理电影大哥的女人| 九九爱精品视频在线观看| 亚洲图色成人| av片东京热男人的天堂| 韩国高清视频一区二区三区| 人人妻人人爽人人添夜夜欢视频| 日韩人妻精品一区2区三区| 成人二区视频| 欧美激情 高清一区二区三区| 免费观看av网站的网址| 18禁裸乳无遮挡动漫免费视频| 毛片一级片免费看久久久久| 国产精品久久久久久久久免| 伦精品一区二区三区| 久久久久久久久久久久大奶| 国产成人免费观看mmmm| 午夜91福利影院| 女人被躁到高潮嗷嗷叫费观| 女性生殖器流出的白浆| 少妇人妻 视频| 久久青草综合色| 精品国产超薄肉色丝袜足j| 人妻 亚洲 视频| 欧美亚洲 丝袜 人妻 在线| 两个人看的免费小视频| 国产欧美日韩一区二区三区在线| 伦理电影大哥的女人| 日本午夜av视频| 久久人人爽av亚洲精品天堂| 韩国高清视频一区二区三区| 国产在线免费精品| 99久久中文字幕三级久久日本| 午夜日韩欧美国产| 国产欧美日韩一区二区三区在线| 成年人免费黄色播放视频| 熟女电影av网| 99香蕉大伊视频| 97在线视频观看| 热re99久久精品国产66热6| 亚洲国产日韩一区二区| 欧美日韩精品成人综合77777| 黑人欧美特级aaaaaa片| 晚上一个人看的免费电影| 亚洲av电影在线进入| 色婷婷av一区二区三区视频| 久久韩国三级中文字幕| 亚洲精品国产av蜜桃| 人体艺术视频欧美日本| av片东京热男人的天堂| 人妻人人澡人人爽人人| 国产一区有黄有色的免费视频| 欧美最新免费一区二区三区| 久久久久久伊人网av| 亚洲精品国产av成人精品| 免费黄频网站在线观看国产| 寂寞人妻少妇视频99o| 制服人妻中文乱码| 激情视频va一区二区三区| 宅男免费午夜| 91aial.com中文字幕在线观看| 久久青草综合色| 三上悠亚av全集在线观看| 亚洲经典国产精华液单| 日韩一卡2卡3卡4卡2021年| 夜夜骑夜夜射夜夜干| 国产精品久久久av美女十八| 中文天堂在线官网| 午夜久久久在线观看| 在线 av 中文字幕| 久热这里只有精品99| 欧美精品一区二区大全| 成人国语在线视频| 在线 av 中文字幕| 秋霞伦理黄片| 2022亚洲国产成人精品| 久久久精品94久久精品| 欧美日韩av久久| 美女午夜性视频免费| 国产精品 国内视频| 欧美中文综合在线视频| 亚洲欧美日韩另类电影网站| 18禁国产床啪视频网站| 国产成人精品久久久久久| 午夜久久久在线观看| 人妻系列 视频| 国产伦理片在线播放av一区| 哪个播放器可以免费观看大片| 两个人看的免费小视频| 一级黄片播放器| 久久av网站| 久久精品亚洲av国产电影网| 久久ye,这里只有精品| 成人二区视频| 亚洲经典国产精华液单| 七月丁香在线播放| 美国免费a级毛片| 亚洲国产看品久久| 男的添女的下面高潮视频| www.熟女人妻精品国产| 少妇熟女欧美另类| 国产精品香港三级国产av潘金莲 | 精品午夜福利在线看| 一二三四在线观看免费中文在| 高清视频免费观看一区二区| 搡女人真爽免费视频火全软件| av卡一久久| 久久久久久人妻| 中文字幕人妻丝袜制服| 男女高潮啪啪啪动态图| 麻豆乱淫一区二区| 免费少妇av软件| 午夜福利在线免费观看网站| 777米奇影视久久| 在线 av 中文字幕| 两个人免费观看高清视频| 制服诱惑二区| 亚洲国产欧美网| 成年女人在线观看亚洲视频| 国产精品国产三级专区第一集| 黄色怎么调成土黄色| 99久久人妻综合| 夫妻性生交免费视频一级片| 亚洲国产欧美在线一区| 午夜激情av网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久99精品国语久久久| 日日摸夜夜添夜夜爱| 男女午夜视频在线观看| 精品国产露脸久久av麻豆| 亚洲成国产人片在线观看| 99久久综合免费| 精品亚洲成国产av| 日韩人妻精品一区2区三区| 亚洲,一卡二卡三卡| 丝袜喷水一区| 日韩人妻精品一区2区三区| 美女xxoo啪啪120秒动态图| 丝袜在线中文字幕| 国产一区二区三区综合在线观看| 国产成人91sexporn| 亚洲综合精品二区| 成年动漫av网址| 日本色播在线视频| 热re99久久国产66热| av在线老鸭窝| 亚洲天堂av无毛| 九色亚洲精品在线播放| 人体艺术视频欧美日本| 色视频在线一区二区三区| 国产乱人偷精品视频| a级毛片在线看网站| 亚洲综合精品二区| 十八禁网站网址无遮挡| 丰满乱子伦码专区| 老女人水多毛片| 91aial.com中文字幕在线观看| 亚洲一区中文字幕在线| 90打野战视频偷拍视频| 9色porny在线观看| 考比视频在线观看| 99热网站在线观看| 婷婷色av中文字幕| 1024视频免费在线观看| av不卡在线播放| 大陆偷拍与自拍| 日本爱情动作片www.在线观看| 成年动漫av网址| 日韩中文字幕视频在线看片| 日韩制服丝袜自拍偷拍| 美女大奶头黄色视频| 国语对白做爰xxxⅹ性视频网站| 在线精品无人区一区二区三| 国产精品三级大全| 黑人巨大精品欧美一区二区蜜桃| 性少妇av在线| 欧美国产精品一级二级三级| 精品人妻一区二区三区麻豆| 免费观看a级毛片全部| 欧美日韩精品网址| 天堂中文最新版在线下载| 啦啦啦中文免费视频观看日本| 另类亚洲欧美激情| 在线观看三级黄色| 自线自在国产av| 中文欧美无线码| a 毛片基地| 精品国产乱码久久久久久男人| 啦啦啦中文免费视频观看日本| 免费看不卡的av| 亚洲内射少妇av| 新久久久久国产一级毛片| 亚洲欧洲国产日韩| 韩国高清视频一区二区三区| 男女啪啪激烈高潮av片| 美女国产高潮福利片在线看| 亚洲av电影在线进入| 一区二区三区四区激情视频| 亚洲视频免费观看视频| 在线观看三级黄色| 日本91视频免费播放| 精品一区二区免费观看| 免费看av在线观看网站| 寂寞人妻少妇视频99o| 在线观看免费高清a一片| 高清视频免费观看一区二区| 国产激情久久老熟女| 在线观看一区二区三区激情| 少妇人妻久久综合中文| 在线亚洲精品国产二区图片欧美| 午夜av观看不卡| 建设人人有责人人尽责人人享有的| 一级片免费观看大全| www日本在线高清视频| 久久久久久久久久久免费av| 综合色丁香网| a级片在线免费高清观看视频| 国产男女内射视频| 亚洲少妇的诱惑av| 制服诱惑二区| 一区二区三区激情视频| 电影成人av| 婷婷色综合www| 色网站视频免费| 美女福利国产在线| 老汉色av国产亚洲站长工具| 亚洲欧洲精品一区二区精品久久久 | 热99久久久久精品小说推荐| av网站在线播放免费| 欧美在线黄色| 一本—道久久a久久精品蜜桃钙片| 丝瓜视频免费看黄片| 欧美bdsm另类| 18禁观看日本| 69精品国产乱码久久久| 激情视频va一区二区三区| 亚洲精品日韩在线中文字幕| 免费黄色在线免费观看| 91午夜精品亚洲一区二区三区| 久久久久网色| 最近最新中文字幕大全免费视频 | 大码成人一级视频| 亚洲图色成人| 日韩制服骚丝袜av| 午夜精品国产一区二区电影| 波多野结衣一区麻豆| 蜜桃国产av成人99| 日本午夜av视频| 精品99又大又爽又粗少妇毛片| 中文字幕人妻熟女乱码| 国产无遮挡羞羞视频在线观看| 国产精品成人在线| 性色av一级| 一区二区三区精品91| 性色av一级| 男人操女人黄网站| 激情视频va一区二区三区| 中文欧美无线码| 亚洲精品中文字幕在线视频| 如日韩欧美国产精品一区二区三区| 国产精品久久久久久精品古装| 欧美日韩视频精品一区| 欧美亚洲日本最大视频资源| 中文字幕人妻熟女乱码| 18禁观看日本| 少妇人妻久久综合中文| 精品亚洲成国产av| 国产一区亚洲一区在线观看| 男女午夜视频在线观看| 久久久久视频综合| 999精品在线视频| 一区在线观看完整版| 国产乱来视频区| 久久99精品国语久久久| 婷婷色av中文字幕| 最黄视频免费看| 精品人妻一区二区三区麻豆| 日本vs欧美在线观看视频| av在线老鸭窝| 2018国产大陆天天弄谢| 一区二区日韩欧美中文字幕| 亚洲欧洲精品一区二区精品久久久 | 男女国产视频网站| 久久久欧美国产精品| 久久久久久久久免费视频了| 深夜精品福利| 色视频在线一区二区三区| 女性生殖器流出的白浆| 天天躁日日躁夜夜躁夜夜| 天堂中文最新版在线下载| 黄色视频在线播放观看不卡| 久久免费观看电影| 麻豆乱淫一区二区| 青草久久国产| 中文字幕亚洲精品专区| 精品一区二区三区四区五区乱码 | 久久99精品国语久久久| 美女视频免费永久观看网站| 最新中文字幕久久久久| av在线app专区| 在线观看三级黄色| h视频一区二区三区| 制服人妻中文乱码| 看免费成人av毛片| 久久国产亚洲av麻豆专区| 最近中文字幕2019免费版| 我要看黄色一级片免费的| 一级毛片黄色毛片免费观看视频| 国产精品久久久久成人av| 精品福利永久在线观看| 亚洲精品第二区| 香蕉国产在线看| 黑人巨大精品欧美一区二区蜜桃| 亚洲一级一片aⅴ在线观看| 久久国产精品男人的天堂亚洲| 亚洲国产欧美在线一区| 国产成人午夜福利电影在线观看| 国产人伦9x9x在线观看 | www.熟女人妻精品国产| 亚洲精品美女久久久久99蜜臀 | 在线看a的网站| 久久久久国产精品人妻一区二区| 日本av免费视频播放| 十八禁高潮呻吟视频| 又粗又硬又长又爽又黄的视频| 亚洲精华国产精华液的使用体验| 亚洲 欧美一区二区三区| 美女脱内裤让男人舔精品视频| 春色校园在线视频观看|