Using groundwater levels to estimate recharge

Accurate estimation of groundwater recharge is extremely important for proper management of groundwater systems. Many different approaches exist for estimating recharge. This paper presents a review of methods that are based on groundwater-level data. The water-table fluctuation method may be the most widely used technique for estimating recharge; it requires knowledge of specific yield and changes in water levels over time. Advantages of this approach include its simplicity and an insensitivity to the mechanism by which water moves through the unsaturated zone. Uncertainty in estimates generated by this method relate to the limited accuracy with which specific yield can be determined and to the extent to which assumptions inherent in the method are valid. Other methods that use water levels (mostly based on the Darcy equation) are also described. The theory underlying the methods is explained. Examples from the literature are used to illustrate applications of the different methods. Electronic Publication     

Hydrogeochemistry and isotope analyses used to determine groundwater recharge and flow in low-gradient catchments of the province of Buenos Aires,Argentina

The region known as Pampa Plain, in Argentina, is a vast area characterized by slopes of less than 0.05%. The surface sediment is silty sand, mainly Aeolian, referred to as Pampean Loess, which is a phreatic aquifer unit of utmost importance for the water supply of the region. On account of the slight gradient, hydrogeological analyses using only hydraulic measurements are difficult to perform, often leading to confusing results. Thus, the study presented relies on hydrochemical modeling and isotopic determinations as well. The study area comprises three catchments in the inter-mountainous area corresponding to El Moro, Tamangueyú and Seco creeks, and covering 2,570 km² in the province of Buenos Aires. Measurements of piezometric levels and samples of groundwater, surface water and rainwater were carried out between January and March 2005. Major ion results were analyzed by means of hydrochemical graphs and hydrogeochemical modelling using NETPATH. The resulting data show that it is possible to identify local changes in recharge, flow direction and stream/groundwater relationship by using hydrochemical and isotopic information, which may become a useful and more precise tool for the study of particularly flat landscapes.     

Radon in groundwater as a tracer to assess flow velocities: two test cases from Israel

 Two test cases from Israel are presented herein employing the decay rate of radon along the flow path to assess groundwater flow velocities. Groundwater flow reaching the fault zone emerges in several places along the rift fault zone as thermal springs because of deep water confinement. The high water temperature of the surface is indicative of high vertical flow velocities, which maintains the original high temperatures. Knowing the Rn content at a source point and at a given down-gradient, and assuming no Rn addition from the water itself or along the flow path, one can calculate the flow velocity based on the Rn half-life time. The decay of Rn in western Galilee was found to be ∼570–150 pCi/l, and in the Dead Sea area from 5000–2000 pCi/l along a respective flow path of 1000 and 200 m, Based on the above, the calculated flow velocities were compared with those obtained from pumping tests in the study area. The method is applicable, because of the short Rn half-life, to cases of high Rn contents, short distances and high flow velocities. Received: 18 January 2000 · Accepted: 21 March 2000     

调查地下水流的1m深地温测量反演问题研究

为河流堤坝渗漏、山体滑坡等各种地质灾害的防治提供一个有效、经济、简单的新方法为研究目的,基于浅层地温测量研究地下水流,开发了1m深地温测量调查地下水流的反演解释系统,进一步研究1m深地温图推断地下水流的定量解释方法,该方法能够推断地下水流的规模及埋深,具有理论简单、解释容易等特点,从理论及应用效果上验证了该方法的有效性,若与地温曲线拟合法联合应用,将提高拟合反演速度及精度。     

Post-drilling destabilization of temperature profile in borehole Yaxcopoil-1, Mexico

As part of the International Continental Scientific Drilling Program (ICDP), the 1.5-km-deep borehole Yaxcopoil-1, located in the Chixculub meteor impact structure in Mexico, has undergone further study after drilling operations ceased. Temperature logs were repeated ten times at intervals 0.3–0.8, 15, 24 and 34 months after borehole shut-in. The logs bear a distinct signature of transient heat transfer by groundwater flow manifested by a gradual distortion of the linear temperature profile when a cold wave of 0.8–1.6°C amplitude was detected propagating downward from 145 to 312 m at a rate of 4–6 m/month. To understand the nature of this moving anomaly, a 20-day monitoring of the cold wave was carried out at a depth of 307 m that showed further cooling of 0.6°C during the first 16 days of the passage followed by temperature stabilisation. As an explanation of this unusual phenomenon, a theory is proposed, whereby the drilling mud has accumulated within the overlying and cooler highly porous and permeable karstic rocks during the drilling and migrates downward. The observed migration rate suggests a permeability higher than 10^?11 m². This indicates a high vulnerability to contamination of the only freshwater aquifer in the Yucatan region.     

Use of multiple age tracers to estimate groundwater residence times and long-term recharge rates in arid southern Oman

Multiple age tracers were measured to estimate groundwater residence times in the regional aquifer system underlying southwestern Oman. This area, known as the Najd, is one of the most arid areas in the world and is planned to be the main agricultural center of the Sultanate of Oman in the near future. The three isotopic age tracers ⁴He, ¹⁴C and ³⁶Cl were measured in waters collected from wells along a line that extended roughly from the Dhofar Mountains near the Arabian Sea northward 400 km into the Empty Quarter of the Arabian Peninsula. The wells sampled were mostly open to the Umm Er Radhuma confined aquifer, although, some were completed in the mostly unconfined Rus aquifer. The combined results from the three tracers indicate the age of the confined groundwater is < 40 ka in the recharge area in the Dhofar Mountains, > 100 ka in the central section north of the mountains, and up to and > one Ma in the Empty Quarter. The ¹⁴C data were used to help calibrate the ⁴He and ³⁶Cl data. Mixing models suggest that long open boreholes north of the mountains compromise ¹⁴C-only interpretations there, in contrast to ⁴He and ³⁶Cl calculations that are less sensitive to borehole mixing. Thus, only the latter two tracers from these more distant wells were considered reliable. In addition to the age tracers, δ²H and δ¹⁸O data suggest that seasonal monsoon and infrequent tropical cyclones are both substantial contributors to the recharge. The study highlights the advantages of using multiple chemical and isotopic data when estimating groundwater travel times and recharge rates, and differentiating recharge mechanisms.     

Characterization of the shallow groundwater system in an alpine watershed: Handcart Gulch,Colorado, USA

Water-table elevation measurements and aquifer parameter estimates are rare in alpine settings because few wells exist in these environments. Alpine groundwater systems may be a primary source of recharge to regional groundwater flow systems. Handcart Gulch is an alpine watershed in Colorado, USA comprised of highly fractured Proterozoic metamorphic and igneous rocks with wells completed to various depths. Primary study objectives include determining hydrologic properties of shallow bedrock and surficial materials, developing a watershed water budget, and testing the consistency of measured hydrologic properties and water budget by constructing a simple model incorporating groundwater and surface water for water year 2005. Water enters the study area as precipitation and exits as discharge in the trunk stream or potential recharge for the deeper aquifer. Surficial infiltration rates ranged from 0.1–6.2×10^?5 m/s. Discharge was estimated at 1.28×10^?3 km³. Numerical modeling analysis of single-well aquifer tests predicted lower specific storage in crystalline bedrock than in ferricrete and colluvial material (6.7×10^?5–2.0×10^?3 l/m). Hydraulic conductivity in crystalline bedrock was significantly lower than in colluvial and alluvial material (4.3×10^?9–2.0×10^?4 m/s). Water budget results suggest that during normal precipitation and temperatures water is available to recharge the deeper groundwater flow system.     

Using saline tracers to evaluate preferential recharge in fractured rocks,Floyd County,Virginia, USA

Potassium chloride (KCl) and potassium bromide (KBr) tracers were used to explore the role of geologic structure on groundwater recharge and flow at the Fractured Rock Research Site in Floyd County, Virginia, USA. Tracer migration was monitored through soil, saprolite, and fractured crystalline bedrock for a period of 3 months with chemical, physical, and geophysical techniques. The tracers were applied at specific locations on the ground surface to directly test flow pathways in a shallow saprolite and deep fractured-rock aquifer. Tracer monitoring was accomplished with differential electrical resistivity, chemical sampling, and physical monitoring of water levels and spring discharge. KCl, applied at a concentration of 10,000 mg/L, traveled 160 m downgradient through the thrust fault aquifer to a spring outlet in 24 days. KBr, applied at a concentration of 5,000 mg/L, traveled 90 m downgradient through the saprolite aquifer in 19 days. Tracer breakthrough curves indicate diffuse flow through the saprolite aquifer and fracture flow through the crystalline thrust fault aquifer. Monitoring saline tracer migration through soil, saprolite, and fractured rock provided data on groundwater recharge that would not have been available using other traditional hydrologic methods. Travel times and flowpaths observed during this study support preferential groundwater recharge controlled by geologic structure.     

Relationship between groundwater flow estimated by soil temperature and slope failures caused by heavy rainfall, Shikoku Island, southwestern Japan

Most of the slope failure disasters in a humid area such as that of Japan are caused by heavy rain. However, even for the case of heavy rainfall that occurs once in every 10 years, total area of slope failures seldom exceeds 10% of a watershed. From this background, we focused on the vein-like groundwater flows that increase pore-water pressure, and clarified the relationship between distributions of slope failures and groundwater veins. In this study, a 1-m-depth ground temperature survey and water-chemistry analyses at springs and boreholes were conducted in Zentoku area of Shikoku Island, southwestern Japan, to grasp the distribution of groundwater veins and their sources. Subsequently, slope-stability was analyzed to investigate the relationship between groundwater veins and slope failures at study sites. These results lead to the following conclusions: The slope failures appear to concentrate around shallow groundwater veins and groundwater veins rising from deep layers. This means that slope failures caused by these groundwater veins in addition to rainfall. Two types of groundwater originate in the deep layers: one has short storage time as indicated by the fact that dissolved substances are low; the other is stored for a lengthy period as noted by a high concentration of dissolved substances. By combining the results of stability analyses and distribution of groundwater veins, it is suggested that prediction of zones with high potential for slope failure can be more accurate.     

Groundwater recharge and discharge processes in the Jakarta groundwater basin, Indonesia

Proper management of groundwater resources requires knowledge of the processes of recharge and discharge associated with a groundwater basin. Such processes have been identified in the Jakarta groundwater basin, Indonesia using a theory that describes the simultaneous transfer of heat and fluid in a porous medium. Temperature-depth profiles in monitoring wells are used to determine the geothermal gradient. To examine the rules of groundwater flow in the distortion of the isotherms in this area, several methods are compared. Subsurface temperature distribution is strongly affected by heat advection due to groundwater flow. Under natural flow conditions, the recharge area is assumed to occur in the hills and uplands, which are located on the periphery of the Jakarta basin, and the discharge area is located in the central and northern part of the Jakarta groundwater basin. A transition area, which could act as local recharge and discharge areas, occupies the middle of the lowland. Subsurface temperatures show good correlation with the groundwater flow conditions, and the data yield important information on the location of recharge and discharge areas.     

Multi-aquifer pumping test to determine cutoff wall length for groundwater flow control during site excavation in Tokyo, Japan

Multi-aquifer pumping tests, using a multi-screen pumping well and multi-level piezometers, were carried out for groundwater flow control in a large-scale excavation site in Tokyo, Japan. The site was underlain by multi-layered confined aquifers. In the tests, pumping was carried out using a multi-aquifer pumping well in which a screen depth was chosen arbitrarily. Changes in groundwater pressure heads in each aquifer were measured at each screen position of the multi-aquifer pumping well. Hydraulic conductivity (K) and specific storage (S _s) of not only aquifers, but also for low permeability layers between the aquifers, were estimated using the Cooper-Jacob method, and calibrated by a finite element method (FEM) groundwater model. Four different cutoff wall lengths were assumed for final excavation depth, and correlations among wall length, pumping discharge and drawdown at the back of the cutoff wall were obtained from simulations using the K and S _s parameters in the FEM model. Then, the most suitable wall length was selected based on the simulated correlations considering environmental condition, construction period and cost of the cutoff wall.     

The use of temperature,stable isotopes and water quality to determine the pattern and spatial extent of groundwater flow: Nagaoka area,Japan

The purpose of this study was to trace the groundwater flow system in the Nagaoka area using subsurface temperature distribution, stable isotopes and water quality. Temperature-depth profiles were measured four times in observation wells during the period 2000–2002. Water was sampled from both observation and pumping wells during the same period. A combination of tracers was used to conceptualize the groundwater flow system. The flow system appears to be divided into shallow, intermediate and deep systems. It is clear that land use, pumping and the infiltration rate are affecting the shallow flow system.An erratum to this article can be found at     

Solute and isotope constraint of groundwater recharge simulation in an arid environment,Abu Dhabi Emirate,United Arab Emirates

A simple, well-constrained simulation of solute increases in a downgradient direction was used in the shallow unconfined aquifer of eastern Abu Dhabi Emirate (United Arab Emirates). The simulation indicates that the observed exponential increase in solute concentrations results from a combination of upward transport of solutes from underlying mudstones and evaporites, and groundwater losses by evaporation. Groundwater recharge and discharge flux in unconfined regional aquifers in arid regions commonly are difficult to estimate because there are few constraints on the flux of water lost or gained from the system. Total dissolved solids (TDS) and deuterium isotopes (δ²H) in groundwater are used to constrain estimated fluxes to the shallow aquifer of eastern Abu Dhabi Emirate. Vertical upward transport of solutes from underlying mudstones and evaporites accounts for solute increases along approximately the first 80 km of the simulated flow path, but a combination of upward solute transport and evaporation is necessary to explain observed solute concentrations beyond 80 km. Mobilization and transport of solutes in the unsaturated zone by recharging precipitation is not a significant factor.

---

Résumé Une simulation simple et bien contrainte de l’augmentation des solutés dans la direction de l’écoulement a été réalisée sur l’aquifère phréatique et libre de l’Est de l’émirat d’Abu Dhabi (émirats Arabes Unis). La simulation indique que l’augmentation observée de solutés, de forme exponentielle, est le produit de la combinaison du transport ascendant de solutés provenant des mudstones et des évaporites sous-jacents, et de l’évaporation des eaux souterraines. La recharge des eaux souterraines et les flux de vidange dans les aquifères régionaux libres des régions arides sont habituellement difficiles à estimer du fait du peu de contraintes agissant sur la perte ou le gain d’eau dans le système. La charge dissoute totale (TDS, en anglais) et les isotopes deutérium (δ²H) des eaux souterraines sont utilisés pour contraindre l’estimation des flux en direction de l’aquifère phréatique de l’Est de l’émirat d’Abu Dhabi. Le transport vertical ascendant des solutés à partir des mudstones et des évaporites fait parti des processus qui concentrent les solutés approximativement le long des 80 premiers kilomètres du trajet des écoulements simulés, mais une combinaison du transport ascendant des solutés et de l’évaporation est nécessaire pour expliquer les concentrations observées en soluté au delà des 80 km. La mobilisation et le transport des solutés par les précipitations efficaces dans la zone non saturée ne représentent pas des facteurs significatifs.

---

Resumen Se utilizó una simulación simple y muy restringida, de incrementos de soluto en la dirección del flujo, en el acuífero libre poco profundo del Emirato oriental de Abu Dhabi, (Emiratos árabes Unidos). La simulación indica que el aumento exponencial observado en las concentraciones del soluto, resulta de una combinación de transporte ascendente de solutos desde las evaporitas y lodolitas subyacentes, y por pérdidas del agua subterránea por evaporación. Los flujos de recarga y descarga de agua subterránea, en los acuíferos libres regionales en regiones áridas, normalmente son difíciles estimar, porque hay algunas restricciones en el flujo de agua perdido o ganado por el sistema. Los sólidos disueltos totales (TDS) y los isótopos de deuterio (δ²H) en el agua subterránea, se usan para forzar los flujos estimados hacia el acuífero poco profundo del Emirato oriental de Abu Dhabi. El transporte ascendente vertical de solutos desde las lodolitas y evaporitas subyacentes, involucra los aumentos del soluto a lo largo de aproximadamente los primeros 80 Km. de la dirección de flujo simulada, pero es necesaria una combinación de transporte ascendente del soluto y la evaporación para explicar las concentraciones del soluto observadas más allá de 80 Km. La movilización y transporte de solutos en la zona no saturada, por recarga debida a precipitación, no es un factor significante.

     

Evidence for deep groundwater flow and convective heat transport in mountainous terrain, Delta County, Colorado, USA

The Tongue Creek watershed lies on the south flank of Grand Mesa in western Colorado, USA and is a site with 1.5 km of topographic relief, heat flow of 100 mW/m², thermal conductivity of 3.3 W m^–1 °C^–1, hydraulic conductivity of 10^-8 m/s, a water table that closely follows surface topography, and groundwater temperatures 3–15°C above mean surface temperatures. These data suggest that convective heat transport by groundwater flow has modified the thermal regime of the site. Steady state three-dimensional numerical simulations of heat flow, groundwater flow, and convective transport were used to model these thermal and hydrological data. The simulations provided estimates for the scale of hydraulic conductivity and bedrock base flow discharge within the watershed. The numerical models show that (1) complex three-dimensional flow systems develop with a range of scales from tens of meters to tens of kilometers; (2) mapped springs are frequently found at locations where contours of hydraulic head indicate strong vertical flow at the water table, and; (3) the distribution of groundwater temperatures in water wells as a function of surface elevation is predicted by the model.     

Using CRD method for quantification of groundwater recharge in the Gaza Strip, Palestine

Rainfall is the main source of groundwater recharge in the Gaza Strip area in Palestine. The area is located in the semi-arid zone and there is no source of recharge other than rainfall. Estimation of groundwater recharge from rainfall is not an easy task since it depends on many uncertain parameters. The cumulative rainfall departure (CRD) method, which depends on the water balance principle, was used in this study to estimate the net groundwater recharge from rainfall. This method does not require much data as is the case with other classical recharge estimation methods. The CRD method was carried out using optimisation approach to minimise the root mean square error (RMSE) between the measured and the simulated groundwater head. The results of this method were compared with the results of other recharge estimation methods from literature. It was found that the results of the CRD method are very close to the results of the other methods, but with less data requirements and greater ease of application. Based on the CRD method, the annual amount of groundwater recharge from rainfall in the Gaza Strip is about 43 million m³. An erratum to this article can be found at     

Application of environmental groundwater tracers at the Sulphur Bank Mercury Mine, California, USA

Boron, chloride, sulfate, δD, δ¹⁸O, and ³H concentrations in surface water and groundwater samples from the Sulphur Bank Mercury Mine (SBMM), California, USA were used to examine geochemical processes and provide constraints on evaporation and groundwater flow. SBMM is an abandoned sulfur and mercury mine with an underlying hydrothermal system, adjacent to Clear Lake, California. Results for non-³H tracers (i.e., boron, chloride, sulfate, δD, and δ¹⁸O) identify contributions from six water types at SBMM. Processes including evaporation, mixing, hydrothermal water input and possible isotopic exchange with hydrothermal gases are also discerned. Tritium data indicate that hydrothermal waters and other deep groundwaters are likely pre-bomb (before ~1952) in age while most other waters were recharged after ~1990. A boron-based steady-state reservoir model of the Herman Impoundment pit lake indicates that 71–79% of its input is from meteoric water with the remainder from hydrothermal contributions. Results for groundwater samples from six shallow wells over a 6–month period for δD and δ¹⁸O suggests that water from Herman Impoundment is diluted another 3% to more than 40% by infiltrating meteoric water, as it leaves the site. Results for this investigation show that environmental tracers are an effective tool to understand the SBMM hydrogeologic regime.

---

Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.

---

Résumé Les concentrations en bore, chlorure, sulfate, δD, δ¹⁸O, et ³H d’échantillons d’eaux de surface et souterraine prélevés dans le banc de soufre de la mine de mercure (SBMM en anglais) en Californie, USA, ont été utilisées pour étudier les processus géochimiques et pour fournir des contraintes à l’évaporation et à l’écoulement des eaux souterraines. La SBMM est une mine de soufre et de mercure abandonnée, adjacente au lac Clear en Californie et sous laquelle se trouve un système hydrothermal.Les résultats des traceurs autres que le tritium (bore, chlorure, sulfate, δD, et δ¹⁸O) ont permis d’identifier des contributions de six types d’eaux à SBMM. Des processus tels que l’évaporation, le mélange, l’entrée d’eau hydrothermale et de possibles échanges isotopiques avec des gaz hydrothermaux ont également été identifiés. Les données de tritium montrent que les eaux hydrothermales et d’autres eaux profondes sont probablement d’age antérieure à la bombe (avant ~1952), alors que la plupart des autres eaux sont issues de la recharge après ~1990. Un modèle de réservoir représentant le lac situé dans la partie Herman Impoundment de l’ancienne mine, en régime permanent et basé sur le bore, montre que 71–79% de l’eau provient des précipitations, le reste provenant de contributions hydrothermales. Les résultats de δD et δ¹⁸O pour des échantillons d’eau souterraine de six puits peu profonds sur une période de 6 mois suggèrent que l’eau de Hermann Impoundment est encore diluée entre 3% jusqu’à plus de 40% lorsqu’elle quitte le site du fait de l’infiltration d’eau météorique. Les résultats de cette étude montrent que les traceurs environnementaux constituent un outil efficace pour comprendre le régime hydrogéologique de la SBMM.

---

Resumen Se han utilizado datos de concentraciones de boro, cloruros, sulfatos, δD, δ¹⁸O, y ³H en muestras de aguas superficiales y subterráneas procedentes de la Mina de Mercurio Sulphur Bank (SBMM), California, USA, para estudiar los procesos geoquímicos y caracterizar la evaporación y el flujo del agua subterránea. SBMM es una mina abandonada de azufre y mercurio con un sistema hidrotermal subyacente, cercano al Lago Clear, California.Los resultados de los trazadores que no son ³H (por ejemplo, boro, cloruros, sulfatos δD, y δ¹⁸O) identifican las contribuciones de seis tipos de agua en la SBMM. Se han identificado diversos procesos, como evaporación, mezcla, entrada de agua hidrotermal y la posibilidad de intercambios isotópicos con gases hidrotermales. Los datos de tritio indican que las aguas hidrotermales y otras aguas subterráneas profundas son probablemente pre-bombas (previas a ~1952) en edad mientras que la mayoría de las otras aguas se han recargado después de ~1990. Un modelo estacionario basado en el boro de la fosa del lago Herman Impoundment indica que el 71–79% de su entrada procede de agua meteórica con restos de contribución hidrotermal. Los resultados para las muestras de aguas subterráneas procedentes de seis pozos superficiales en un periodo de 6 meses para δD y δ¹⁸O sugieren que el agua procedente de Herman Impoundment se diluye en otro 3% hasta más del 40% debido a la infiltración de agua meteórica., que sale del sitio. Los resultados de esta investigación muestran que los trazadores ambientales son una herramienta efectiva para entender el régimen hidrogeológico de la SBMM.

     

Using soil and contaminant properties to assess the potential for groundwater contamination to the lower Great Lakes,USA

Contaminant risk factors in surface soil were evaluated within the urbanized Rouge River watershed in southeastern Michigan, USA, which includes metropolitan Detroit. An analytical risk factor model and Geographic Information Systems overlays were used to quantify and characterize the potential impacts of five categories of contaminants including DNAPLs (dense nonaqueous phase liquids), LNAPLs (light nonaqueous phase liquids), PAHs (polynuclear aromatic hydrocarbons), PCBs (polychlorinated biphenyls), and lead. The results indicate that DNAPL compounds released into sand, moraine, and sandy and silty clay soil types have the greatest potential to affect groundwater, and impact the Great Lakes ecosystem and the public health.     

Evaluation of change in subsurface thermal environment due to groundwater flow in the Tokyo Lowland,Japan

Information on the distribution of subsurface temperature and hydraulic heads at 24 observation wells in and around the Tokyo Lowland, the eastern part of the Tokyo Metropolitan area, were examined to make clear the relationship between groundwater and the subsurface thermal environment in the urban area. Minimums in temperature–depth logs due to subsurface temperature increasing at shallow parts were recognized in 21 wells. This fact shows subsurface temperature is affected by ground surface warming in almost all of this area. Deeper than minimums, where the effects of surface warming became relatively small, regional variation is observed as follows: high temperatures are shown in the central part to the southern part, and low temperatures shown in the inland to eastern part. The high temperature area corresponds to an area where the lower boundary of groundwater flow is relatively shallow. This area corresponds also to an area with severe land subsidence resulting from excessive groundwater pumping. It is considered that this high temperature area is formed by the effects of upward groundwater flow affected by hydrogeological conditions and pumping. On the other hand, a comparison between past data (1956–1967) and present data (2001–2003) revealed widespread decreasing temperature in the inland area. This is explained by downward groundwater flow based on an analysis of temperature–depth logs. This fact suggests that subsurface temperature is not only increasing from the effects of surface warming but also decreasing from the effects of groundwater environment change due to pumping.