Hydrogeochemistry of the Goulburn Valley region of the Murray Basin, Australia: implications for flow paths and resource vulnerability

Groundwater in the Goulburn region of the Murray Basin (Australia) contains solutes derived mainly from evapotranspiration of rainfall-derived marine aerosols, silicate dissolution, and ion exchange. ¹⁴C data indicate that groundwater in the Shepparton Formation recharges vertically across the region, whereas groundwater in the Calivil–Renmark Formation shows a greater component of lateral flow. The overall pattern of geochemical and stable isotope variations implies that long-term vertical groundwater flow into the Calivil–Renmark Formation has occurred over thousands of years. Elevated C, N, and F concentrations, together with variable Cl/Br ratios and ¹⁸O values, suggest that short-term (years to decades) flow of surface water into the Calivil–Renmark Formation also occurs locally. The high degree of vertical flow implies that the high-quality groundwater resources of the Calivil–Renmark Formation are vulnerable to surface contamination. Groundwater in both the Shepparton and Calivil–Renmark Formations yields ¹⁴C ages of up to 20 ka that imply that, overall, recharge rates are low and that, consequently, the groundwater resource in both formations could be impacted by over abstraction.

---

Resumen Las aguas subterráneas de la región de Goulburn, en la cuenca del río Murray (Australia), contienen solutos procedentes principalmente de la evaporación de agua de lluvia, disolución de silicatos e intercambio iónico. Aunque la variabilidad espacial de la hidrogeoquímica es reducida, el contenido total de sales disueltas en las aguas subterráneas de la formación no confinada de Shepparton es generalmente mayor que el del acuífero inferior confinado de la Formación Calivil–Renmark. Los datos de ¹⁴C, ³H y geoquímica sugieren que la recarga del acuífero de la Formación Shepparton es por infiltración (vertical) en toda la región, mientras que la Formación Calivil–Renmark tiene una mayor proporción de aportes laterales. Sin embargo, a escala local se produce recarga vertical en la Formación Calivil–Renmark, hecho que indica que este recurso es vulnerable a la contaminación. Las aguas subterráneas de ambas formaciones tienen edades de hasta 20.000 años, según datos de ¹⁴C, lo que implica que las tasas de recarga son bajas y, en consecuencia, que podría haber impactos asociados a la sobreexplotación del recurso. Las aguas subterráneas recientes de la región se caracterizan por relaciones molares muy variables de Cl/Br (de 50 a 1.200), reflejando la química de una precipitación semiárida que se ve modificada posteriormente por disolución de halita, asociada con terrenos secos y riego con aguas salinas. Las aguas subterráneas más profundas presentan relaciones más uniformes de Cl/Br (entre 500 y 800), y probablemente reflejan que la recarga se produjo en condiciones climáticas más húmedas que las actuales.

---

Résumé Les eaux souterraines de la région de Goulburn dans le bassin de Murray (Australie) contiennent des solutés provenant essentiellement de lévaporation des pluies, de la dissolution de silicates et déchanges dions. Alors quil y a peu de variations spatiales du chimisme des eaux souterraines, les concentrations en sels dissous totaux des eaux souterraines de laquifère libre de la formation de Shepparton sont en général plus élevées que celles des eaux souterraines de laquifère captif sous-jacent de la formation de Calivil-Renmark. Des données de carbone-14, de tritium et de géochimie indiquent que laquifère de la formation de Shepparton est rechargé verticalement dans toute la région, tandis que laquifère de la formation de Calivil-Renmark présente une plus forte composante découlement latéral. Toutefois, localement, il se produit une recharge verticale de laquifère de la formation de Calivil-Renmark, ce qui indique que cette ressource est vulnérable à la contamination. Les eaux souterraines des formations de Shepparton et de Calivil-Renmark présentent toutes deux des âges carbone-14 de plus de 20 ka, ce qui implique des taux de recharge très faibles et, par conséquent, que leur ressource peut être affectée par des prélèvements excessifs. Les eaux souterraines récentes dans cette région sont caractérisées par un rapport molaire Cl/Br très variable (de 50 à 1200), qui reflète le chimisme de pluies en région semi-aride modifié en conséquence par la dissolution de la halite associée à la salinité due aux terres sèches et à lirrigation. Les eaux souterraines plus profondes possèdent des rapports Cl/Br plus uniformes (de lordre de 500 à 800) et reflètent probablement des eaux de recharge sous des conditions climatiques plus humides que les conditions actuelles.

     

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.     

Groundwater flow through Pleistocene glacial deposits in the rapidly urbanizing Rouge River–Highland Creek watershed,City of Scarborough,southern Ontario,Canada

The city of Scarborough lies on the eastern margin of the Greater Toronto Area of southern Ontario, Canada, along the northern coastline of Lake Ontario. The City has a population of 500,000 and is presently one of the fastest growing communities in Canada. The City is expanding northwards onto rural land on the south slope of the large Pleistocene glacial Oak Ridges Moraine system. The moraine system is underlain by a thick (150 m) succession of tills, sands and gravels and is a regionally-significant recharge area for three principle aquifer systems that discharge to numerous watercourses that flow to Lake Ontario. Protection of deeper aquifers from surface-generated urban contaminants is a particular concern. A groundwater flow model using Visual MODFLOW was developed for the 350-km² Rouge River–Highland Creek (RRHC) drainage basin using an extensive GIS-based collection of subsurface geological, geophysical and hydrogeological data, maps of land use and surficial geology. The RRHC model was calibrated against point water level data, known potentiometric surfaces of the principal aquifers and baseflow measurements from streamflow gauging stations and determined to be within acceptable limits. Water balance calculations indicate that 70% of the basin recharge (106,000 m³/day) enters the Upper Aquifer along the crest and immediate flanks of the Oak Ridges Moraine. To the south, Upper Aquifer water moving through fractured till aquitards accounts for more than 75% of recharge to deeper aquifers. Water quality data confirm previous observations that urban- and rural-sourced contaminants (chlorides and nitrates) present in Upper Aquifer waters are moving rapidly into deeper aquifers. Some 83% of total RRHC recharge water is ultimately discharged as baseflow to creeks draining to Lake Ontario; the remainder discharges to springs and along eroding lakeshore bluffs. Model results demonstrate that deeper aquifers are poorly protected from urban contaminants and that long-term protection of ground and surface water quality has to be a priority of municipal planners if the resource is not to be severely degraded. Electronic Publication     

Groundwater flow modelling within a coastal alluvial plain setting using a high-resolution hydrofacies approach; Bells Creek plain, Australia

Ground penetrating radar (GPR) has proved to be an extremely useful geophysical tool, in conjunction with direct geological data, to develop a realistic, macroscopic, subjective-based conceptual model of aquifer architecture within a shallow coastal alluvial plain. Subsequent finite-difference groundwater modelling has not only enabled determination of the dominant groundwater flow paths for the plain, but has also quantified the effects of within-facies and between-facies sedimentary heterogeneity on those flow paths. The interconnection of narrow, unconfined alluvial channels and a broad, semi-confined alluvial delta is ensuring that most fresh groundwater that enters the plain in the form of precipitation or recharge from lateral bedrock hills, is discharged into the eastern coastal wetlands via that alluvial delta aquifer.     

Groundwater flow modeling in the Zhangye Basin,Northwestern China

The Zhangye basin is in the middle reaches of the Heihe River, northwestern China. Heavy abstraction of groundwater since the 1970s in the area is for agricultural, industrial and drinking water supplies and has led to a substantial decline in the potentiometric surface. A three-dimensional regional numerical groundwater flow model, calibrated under transient conditions, has been developed and used to predict the drawdown for the period from 2000 to 2030 under two different groundwater management scenarios.     

Constraining flow paths of saline groundwater at basin margins using hydrochemistry and environmental isotopes: Lake Cooper,Murray Basin,Australia

Solutes in saline groundwater (total dissolved solids up to 37 000 mg/L) in the Lake Cooper region in the southern margin of the Riverine Province of the Murray Basin are derived by evapotranspiration of rainfall with minor silicate, carbonate and halite dissolution. The distribution of hydraulic heads, salinity, percentage modern carbon (pmc) contents, and Cl/Br ratios imply that the groundwater system is complex with vertical flow superimposed on lateral flow away from the basin margins. Similarities in major ion composition, stable (O, H, and C) isotope, and ⁸⁷Sr/⁸⁶Sr ratios between groundwater from the shallower Shepparton Formation and the deeper Calivil – Renmark aquifer also imply that these aquifers are hydraulically interconnected. Groundwater in the deeper Calivil – Renmark aquifer in the Lake Cooper region has residence times of up to 25 000 years, implying that pre-land-clearing recharge rates were <1 mm/y. As in other regions of the Murray Basin, the low recharge rates account for the occurrence of high-salinity groundwater. Shallow (<20 m) groundwater yields exclusively modern ¹⁴C ages and shows a greater influence of evaporation over transpiration. Both these observations reflect the rise of the regional water-table following land clearing over the last 200 years and a subsequent increase in recharge to 10 – 20 mm/y. The rise of the regional water-table also has increased vertical and horizontal hydraulic gradients that may ultimately lead to the export of salt from the Lake Cooper embayment into the adjacent fresher groundwater resources.     

Impact of rehabilitation of Assiut barrage, Nile River, on groundwater rise in urban areas

To make optimum use of the most vital natural resource of Egypt, the River Nile water, a number of regulating structures (in the form of dams and barrages) for control and diversion of the river flow have been constructed in this river since the start of the 20th century. One of these barrages is the Assiut barrage which will require considerable repairs in the near future. The design of the rehabilitation of the barrage includes a headpond with water levels maintained at a level approximately 0.60 m higher than the highest water level in the headpond of the present barrage. This development will cause an increase of the seepage flow from the river towards the adjacent agricultural lands, Assiut Town and villages. The increased head pond level might cause a rise of the groundwater levels and impedance of drainage outflows. The drainage conditions may therefore be adversely affected in the so-called impacted areas which comprise floodplains on both sides of the Nile for about 70 km upstream of the future barrage. A rise in the groundwater table, particularly when high river levels impede drainage, may result in waterlogging and secondary salinization of the soil profile in agricultural areas and increase of groundwater into cellars beneath buildings in the urban areas. In addition, a rise in the groundwater table could have negative impact on existing sanitation facilities, in particular in the areas which are served with septic tanks. The impacts of increasing the headpond level were assessed using a three-dimensional groundwater model. The mechanisms of interactions between the Nile River and the underlying Quaternary aquifer system as they affect the recharge/discharge processes are comprehensively outlined. The model has been calibrated for steady state and transient conditions against historical data from observation wells. The mitigation measures for the groundwater rise in the urban areas have been tested using the calibrated mode.     

Groundwater flow and geochemistry in the lower reaches of the Yellow River: a case study in Shandang Province, China

Water samples were collected from the Yellow River and from wells for chemical and isotopic measurement in the counties of Yucheng and Qihe, to which 6–9×10⁸ m³ of water is diverted annually from the Yellow River. A zone of high electrical conductivity (EC) in groundwater corresponds well on the regional scale with a ridge in groundwater level, which is the main flow path through the region, but has a low gradient. The zone of highest EC along this ridge occurs at a position with the lowest ground altitude in the study area. The unique characteristic of the groundwater is the linear relationship among the principal anions as the result of mixing. The mixing effect is confirmed by its isotopic signature, which was then used to calculate the contributions from three sources: rainfall, old water, and diverted water with an average mixing rate of 18, 17, and 65%, respectively. As an indicator of water movement, Cl^– content varies across a wide range in the profile from 30–10 m with a maximum concentration at about 1.2 m depth. Concentrations are relatively stable at about 2 m, which is the average boundary of the saturated and unsaturated zone. The water from the Yellow River has proved to be dominant in mixing in the aquifer in terms of groundwater flow and geochemistry. Electronic Publication     

Groundwater recharge and flow in the Lower Cretaceous Nubian Sandstone aquifer in the Sinai Peninsula, using isotopic techniques and hydrochemistry

The present work was conducted in the Sinai Peninsula (1) to identify the recharge and flow characteristics and to evaluate the continuity of the Lower Cretaceous Nubian Sandstone aquifer; and (2) to provide information for the aquifer's rational appraisal. Isotopic and hydrochemical compositions combined with the geological and hydrogeological settings were used for this purpose. A considerable depletion in isotopic content (oxygen-18 and deuterium) and low d-excess values exist in the studied groundwater, reflecting the contribution of old meteoric water that recharged the aquifer in pluvial times. Modern recharge also occurs from precipitation that falls on the aquifer outcrops. The wide scatter of the data points around the two meteoric lines, the global meteoric water line (GMWL) and Mediterranean meteoric water line (MMWL), in the δ¹⁸O–δD diagram indicates considerable variation in recharge conditions (amount, altitude, temperature, air masses, distances from catchment, overland flow, etc.). The isotopic composition in the El-Bruk area is minimum (¹⁸O=–9.53‰), very close to the average value of the Western Desert Nubian Sandstone (¹⁸O=–10‰), where the local structural and lithologic conditions retard groundwater flow and the main bulk of water becomes noncyclic. The continuity of the aquifer in northern and central Sinai is evidenced by the isotopic similarity between samples taken from above and below the central Sinai Ragabet El-Naam fault, the distribution of potentiometric head, and hydrogeological cross sections. The combination of isotopic composition in terms of ¹⁸O and chemical composition in terms of TDS and salt contents is the basis for separating the studied groundwater into groups that reflect the recharge sources and isotopic and chemical modifications during flow. Electronic Publication     

Groundwater flow underneath mound spring tufas from geophysical surveys in the southwestern Great Artesian Basin,Australia

Artesian springs, which form carbonate mounds at the surface, occur in groups along the southwestern edge of the Great Artesian Basin (GAB), in northern South Australia. Their underground structure and relationship to faulting are not well understood. This study investigated four mound spring groups (Beresford Spring, Warburton Spring, the Bubbler Spring group, and Freeling Springs), which have different geological settings, using a range of geophysical techniques: self-potential (SP), magnetotellurics (MT) and time-domain electromagnetics (TEM). The results confirmed that despite generally similar mound morphologies at the surface, spring vents in different groups have different mechanisms for sourcing water from the GAB aquifer. SP measurements effectively showed where most of the vertical flow occurs at all groups, while MT and TEM models identified the location of faults, other geological structures and the thickness of the Bulldog Shale aquitard. Beresford and Warburton springs are similar in that the spring locations at both are directly related to points on a regional-scale fault, which appears to provide a low permeability pathway to the surface. At the Bubbler Springs group, a regionally aligned fault identified in MT data seems to be bounding an area where the aquifer is much closer to the surface (and the overlying aquitard has been thinned) where the spring abundance is highest. SP responses indicate that shallow flow paths feeding individual vents are unrelated to the fault. At Freeling Springs, the SP suggests that groundwater is flowing up within shallow aquifers, and that spring location is offset from this upward flow, instead located on range-front faults that provide the final conduit to the surface.     

Groundwater and climate change: a sensitivity analysis for the Grand Forks aquifer,southern British Columbia,Canada

The Grand Forks aquifer, located in south-central British Columbia, Canada was used as a case study area for modeling the sensitivity of an aquifer to changes in recharge and river stage consistent with projected climate-change scenarios for the region. Results suggest that variations in recharge to the aquifer under the different climate-change scenarios, modeled under steady-state conditions, have a much smaller impact on the groundwater system than changes in river-stage elevation of the Kettle and Granby Rivers, which flow through the valley. All simulations showed relatively small changes in the overall configuration of the water table and general direction of groundwater flow. High-recharge and low-recharge simulations resulted in approximately a +0.05 m increase and a –0.025 m decrease, respectively, in water-table elevations throughout the aquifer. Simulated changes in river-stage elevation, to reflect higher-than-peak-flow levels (by 20 and 50%), resulted in average changes in the water-table elevation of 2.72 and 3.45 m, respectively. Simulated changes in river-stage elevation, to reflect lower-than-baseflow levels (by 20 and 50%), resulted in average changes in the water-table elevation of –0.48 and –2.10 m, respectively. Current observed water-table elevations in the valley are consistent with an average river-stage elevation (between current baseflow and peak-flow stages).

---

Resumen El acuífero de los Grand Forks, situado al sur de la Columbia Británica central (Canadá) ha sido utilizado como lugar de estudio para modelar la sensibilidad de un acuífero a los cambios en la recarga y el caudal de los ríos de acuerdo con escenarios previstos de cambio climático en la región. Los resultados sugieren que las variaciones en la recarga al acuífero bajo los diversos escenarios, que han sido modelados en régimen estacionario, tienen un impacto mucho menor en las aguas subterráneas que los cambios en el caudal de los ríos Kettle y Granby, que discurren por el valle. Todas las simulaciones muestran diferencias relativamente pequeñas en la configuración regional de los niveles freáticos y en la dirección general del flujo subterráneo. Las simulaciones de recarga elevada y baja causan un incremento de 0,05 m y un decremento de 0,025 m, respectivamente, en los niveles del acuífero. Los cambios de la elevación del río, simulados para reflejar niveles de flujo mayores que los valores pico (en un 20% y un 50%) resultan en cambios medios de los niveles del acuífero de 2,72 m y 3,45 m, respectivamente. Los cambios simulados en la elevación del río para flujos inferiores al caudal de base (en un 20% y en un 50%) provocan descensos en los niveles de 0,48 y 2,10 m, respectivamente. Los niveles actuales del acuífero en el valle son coherentes con una elevación media del nivel en el río (entre el caudal de base actual y los picos de caudal).

---

Résumé L'aquifère de Grand Forks, situé en Colombie britannique (Canada), a été utilisé comme zone d'étude pour modéliser la sensibilité d'un aquifère à des modifications de la recharge et du niveau de la rivière, correspondant à des scénarios envisagés de changement climatique dans cette région. Les résultats font apparaître que les variations de recharge de l'aquifère pour différents scénarios de changement climatique, modélisées pour des conditions de régime permanent, ont un impact sur le système aquifère beaucoup plus faible que les changements du niveau des rivières Kettle et Granby, qui coulent dans la vallée. Toutes les simulations ont montré des différences relativement faibles dans la configuration d'ensemble de la nappe et dans la direction générale des écoulements. Des simulations de conditions de recharge forte et de recharge faible produisent respectivement une remontée de 0,05 m et un abaissement de 0,025 m, approximativement, des cotes de la nappe pour l'ensemble de l'aquifère. Des changements simulés de la cote du niveau de la rivière, pour refléter des niveaux plus hauts que ceux des pics de crues (de 20 et de 50%), produisent respectivement des remontées de la nappe de 2,72 et 3,45 m en moyenne. Des changements simulés de l'altitude du niveau de la rivière, pour refléter des niveaux plus bas que ceux de basses eaux (de 20 et de 50%), produisent respectivement des abaissements de la nappe de 0,48 et 2,10 m en moyenne. Les altitudes courantes observées de la nappe dans la vallée sont cohérentes avec une cote moyenne du niveau de la rivière (entre les niveaux courants de basses eaux et de crues).

---

     

Groundwater flow patterns in the San Luis Valley,Colorado, USA revisited: an evaluation of solute and isotopic data

Groundwater systems in the San Luis Valley, Colorado, USA have been re-evaluated by an analysis of solute and isotopic data. Existing stream, spring, and groundwater samples have been augmented with 154 solute and isotopic samples. Based on geochemical stratification, three groundwater regimes have been identified within 1,200 m of the surface: unconfined, upper active confined, and lower active confined with maximum TDS concentrations of 35,000, 3,500 and 600 mg/L, respectively. The elevated TDS of northern valley unconfined and upper active confined systems result from mineral dissolution, ion exchange and methanogenesis of organic and evaporate lake sediments deposited in an ancient lake, herein designated as Lake Sipapu. Chemical evolutions along flow paths were modeled with NETPATH. Groundwater ages, and δ¹³C, δ²H and δ¹⁸O compositions and distributions, suggest that mountain front recharge is the principle recharge mechanism for the upper and lower confined aquifers with travel times in the northern valley of more than 20,000 and 30,000 ¹⁴C years, respectively. Southern valley confined aquifer travel times are 5,000 ¹⁴C years or less. The unconfined aquifer contains appreciable modern recharge water and the contribution of confined aquifer water to the unconfined aquifer does not exceed 20%.     

新疆孔雀河流域地下水演化及其生态效应

近50 a,西北干旱内陆河流域经历了历史时期的水土资源过度开发和近期的生态保护修复的不同阶段,对区域地下水和生态环境产生了重要影响,但地下水长期演化规律及其生态效应认识尚不清晰。以新疆孔雀河流域为典型研究区,基于最新地下水监测数据与历史数据对比,分析流域地下水流场时空演化规律以及流域地下水对河岸带胡杨林的影响作用,探讨地下水长期演化的生态效应。结果表明：1971—2021年,孔雀河流域因过量集中开采地下水引起区域地下水流动系统发生根本性变化,地下水水位整体呈下降趋势,流域内累计降幅大于40 m的面积达204.69 km2,累计降幅20～40 m的面积为1735.81 km2,累计降幅5～20 m的面积为1018.56 km2,并形成地下水降落漏斗;生态输水前,中下游河道长时间断流,河岸带地下水水位埋深累计降幅超过10 m,导致河岸带胡杨林退化萎缩,2016年实施生态输水工程以后,埋深逐渐上升2～4 m,上升至胡杨林生长临界水位以上,2014—2018年天然胡杨林植被覆盖率增大17%。本研究有助于认识西北干旱内陆河流域地下水和生态环境演化规律,为类似区域地下水可持续开发利用与生态保护提供科学参考。     

Effect of advective mass transfer on field scale fluid and solute movement: Field and modeling studies at a waste disposal site in fractured rock at Oak Ridge National Laboratory, Tennessee, USA

Advective mass transfer is a pore scale mass-transfer process that affects fluid and solute movement between pore domains such as fracture and matrix in a structured porous medium. Mechanistically similar to advection in the advection-dispersion of solutes in non-structured porous medium, it redistributes solutes by moving solute and solvent simultaneously between pore domains. While there is much research on diffusive mass transfer that is often referred to as matrix diffusion, there is a lack of information and study for advective mass transfer in the literature. The objective of this research is to study the effects of advective mass transfer on fluid and solute movement between pore domains. First, field hydraulic measurements at a waste disposal site in fractured rock at Oak Ridge National Laboratory (ORNL), Tennessee, USA, are used to calibrate a fracture-matrix, two-pore-domain groundwater flow model. Latin-hypercube sensitivity analysis suggests that the uncertainty of the calibrated model parameters is small and the calibrated flow model is nearly the optimal. Fracture spacing thus obtained is used to calculate diffusive mass transfer coefficients. The individual effects of advective and diffusive mass transfer on solute movement are then quantitatively evaluated. The calculations indicate that pore structure conceptual models may significantly affect the role of advective mass transfer in field and pore-scale mass transfer. In the particular ORNL field site and with a fracture-matrix pore structure model, contribution of advective mass transfer to solute mass movement is about three to eight orders of magnitude smaller than that of diffusive mass transfer.

---

Resumen La transferencia de masa advectiva es un proceso de transferencia de masa en escala intersticial que afecta el movimiento de fluido y soluto entre ámbitos porosos tal como fractura y matriz en un medio estructurado poroso. Este proceso, mecánicamente similar a la advección en la dispersión-advección de solutos de medios porosos no estructurados, redistribuye los solutos mediante el movimiento simultáneo de soluto y solvente entre ámbitos porosos. Mientras que existe bastante investigación en transferencia difusiva de masa que frecuentemente se conoce como difusión en matriz, existe falta de información y estudio de transferencia advectiva de masa en la literatura. El objetivo de esta investigación es estudiar los efectos de la transferencia advectiva de masa en el movimiento de fluido y soluto entre ámbitos porosos. Primero se utilizaron mediciones hidráulicas de campo en un sitio de depósito de residuos en roca fracturada en el Laboratorio Nacional Oak Ridge (ORNL), Tennessee, USA, para calibrar un modelo de flujo de agua subterránea de ámbito de dos poros fractura-matriz. Análisis de sensitividad hipercúbico-latino sugieren que la incertidumbre de los parámetros del modelo calibrado es pequeña y que el modelo de flujo calibrado es aproximadamente el óptimo. El espaciamiento de fracturas así obtenido se utiliza para calcular los coeficientes de transferencia de masa difusiva. Luego se evalúa cuantitativamente los efectos individuales de transferencia de masa advectiva y difusiva en el movimiento de soluto. Los cálculos indican que los modelos conceptuales de estructura porosa pueden afectar significativamente el papel de transferencia de masa advectiva en escalas de campo e intersticial de transferencia de masa. En el sitio de campo específico ORNL y con un modelo de estructura porosa de matriz-fractura, la contribución de transferencia de masa advectiva al movimiento de masa soluto es aproximadamente tres a ocho órdenes de magnitud más pequeño que la contribución por transferencia de masa difusiva.

---

Résumé Résumé Le transfert advectif de masse est un processus de transfert de masse à léchelle du pore qui affecte les mouvements du fluide et du soluté entre les différents domaines de pores, tel les fractures et la matrice dans un milieu poreux structuré. Mécaniquement similaire à ladvection dans le concept advection-dispersion de solutés dans les milieux non structurés, ce transfert redistribue les solutés simultanément avec le fluide entre les différents domaines poreux. Alors que de nombreuses recherches portent sur les transferts de masse par diffusion, se référant généralement à une diffusion par la matrice, il y a un grand manquement dinformations et détudes sur les transferts de masse par advection dans la littérature. Lobjectif de cette recherche est détudier leffet du transfert advectif de masse de fluide et de soluté entre les différents domaines poreux. Premièrement, les mesures hydrauliques de terrain sur la décharge en milieu fracturé du laboratoire national dOak Ridge ORNL, Tennessee, USA, sont utilisées pour calibrer un modèle hydrogéologique à double porosité fracture-matrice. Lanalyse de sensibilité latin-hypercube suggère que lincertitude sur les paramètres du modèle est faible et que le calibrage est pratiquement optimal. Lespace de fracture résultant permet de calculer les coefficients de transfert de masse par diffusion. Les effets individuels de ladvection et de la diffusion sur les mouvements de solutés sont dés lors évalués. Les calculs indiquent que le modèle conceptuel de la structure des pores peuvent significativement affectés le rôle du transfert advectif de masse à léchelle du pore et du terrain. Dans le cas du site de lORNL et avec un modèle structuré fracture-matrice, la contribution de ladvection au transfert de masse est de lordre de trois-huitième du transfert de masse par diffusion.

     

Application of the groundwater-balance equation to indicate interbasin and vertical flow in two semi-arid drainage basins, Mexico

An analysis of horizontal inflow and outflow in the groundwater-budget equation and the significance for interbasin flow are presented. Two field cases in Mexico, one in the Baja California peninsula and another in central Mexico, highlight the influence of interbasin flow. A significant proportion (approximately 70%) of the abstracted (thermal) groundwater probably originates outside the drainage basin. A conclusion is that a groundwater-balance study is an unsatisfactory method for determining some parameters, such as storativity (S). Specifically, the groundwater-balance approach provides unreliable results when vertical inflow is ignored or cannot be adequately defined. Vertical flow is indicated by the presence of groundwater temperatures as much as 23 °C higher than ambient temperature. Regional faults could be the pathways for upward flow. When vertical inflow is ignored, uncertainty in the estimation of the storativity through regional groundwater-balance calculation results. On the basis of the groundwater-balance equation, a value of S=0.19 appears to represent the confined condition of the developed part of the aquifer; this result is several orders of magnitude higher than would be reasonable according to the geological conditions. Findings are useful in evaluating whether a groundwater resource is being "overexploited". Conclusions are instructive in the application of transient-flow computer models, in which vertical flow of less dense water from beneath is not included. Electronic Publication