Conceptual model of leachate migration in a granular aquifer derived from the integration of multi-source characterization data (St-Lambert,Canada)

Numerical models encompassing source zones and receptors, based on representative conceptual models and accounting for aquifer heterogeneity, are needed to understand contaminant migration and fate; however, aquifer characterization seldom provides the necessary data. This study aimed to develop a workflow for field characterization and data integration, which could: (1) be adapted to the definition of subwatershed-scale aquifer heterogeneity (over 10 km²) and (2) adequately support mass transport model development. The study involved the field investigation of a shallow granular aquifer in a 12-km² subwatershed in Saint-Lambert-de-Lauzon, Canada, in which a decommissioned landfill is emitting a leachate plume managed by natural attenuation. Using proven field methods, the characterization sequence was designed to optimize each method in terms of location, scale of acquisition, density and quality. The emphasis was on the acquisition of detailed indirect geophysical data that were integrated with direct hydraulic and geochemical data. This report focuses on the first qualitative and geostatistical data integration steps of the workflow leading to the development of a hydrogeological conceptual model. This is a prerequisite for further integration steps: prediction of hydrofacies and hydraulic conductivity (K), geostatistical simulations of K, studies of geochemical processes and numerical modeling.     

An integrated approach to shallow aquifer characterization: combining geophysics and geostatistics

We present a method of aquifer characterization that is able to utilize multiple sources of conditioning data to build a more realistic model of heterogeneity. This modeling approach (InMod) uses geophysical data to delineate bounding surfaces within sedimentary deposits. The depositional volumes between bounding surfaces are identified automatically from the geophysical data by a region growing algorithm. Simple geometric rules are used to constrain the growth of the regions in 3-D. The nodes within the depositional volume are assigned to categorical lithologies using geostatistical realizations and a dynamic lookup routine that can be conditioned to field data. The realizations created with this method preserve geologically expected features and produces sharp juxtapositions of high and low hydraulic conductivity lithologies along bounding surfaces. The realizations created with InMod also have higher variance than models created only with geostatistics and honor the volumetric distribution of sediments measured from field data.     

Field scale heterogeneity of redox conditions in till-upscaling to a catchment nitrate model

Point scale studies in different settings of glacial geology show a large local variation of redox conditions. There is a need to develop an upscaling methodology for catchment scale models. This paper describes a study of field-scale heterogeneity of redox-interfaces in a till aquitard within an area of 600?×?600 m. The results showed significant variation of the depths to the redox-interface and thickness of the aquitard. Nitrate was present above the redox-interface but reduced to non-detectable levels a few metres below the interface. An upscaling approach for an area of 92 km² is proposed. Two models are proposed to predict the depth to the redox-interface in the aquitard and the resulting nitrate recharge concentrations to an underlying aquifer. The first model assumes that the depth to the redox-interface reflects the hydraulic head in the aquitard, and the second model assumes that the depth of the redox-interface is randomly distributed according to a log-normal probability distribution function. The upscaling approach using the random redox model estimated recharge concentrations comparable to the observed concentration in the underlying aquifer. The presented upscaling approach is applicable in distributed catchment models where sub-grid variability cannot be represented by the large grids.     

利用观测孔中地下水位潮汐效应计算天津滨海新区含水层参数

在沿海地区,尤其是围海造陆工程形成的陆域地区地下水水位受潮汐影响较大,使传统水文地质试验求取含水层参数存在较大误差。因此通过合理概化地下水在潮汐作用下运动规律,建立数学模型,推导解析公式求取沿海含水层参数具有重要意义。分析天津滨海新区两处观测孔地下水位及潮汐波动特征,在滞后时间不明显的情况下,利用观测孔水位变幅数据计算了含水层水头扩散系数,并根据承压含水层储水系数经验值进一步获得含水层渗透系数。通过两个观测孔分别计算,对比计算结果互相验证发现,该方法取得了令人满意的结果。利用地下水潮汐效应计算含水层参数可以广泛应用于沿海地区水文地质工作中。     

利用CFCs定年数据计算静升盆地含水层渗透系数

文章利用CFCs(氟氯碳化合物)定年数据,计算山西静升盆地含水层渗透系数。结果表明:补给条件较好的区域与补给条件差的深层承压区相比,地下水年龄偏年轻;从冲洪积扇顶部到扇缘,地下水径流逐渐变缓,含水层渗透系数变小;盆地南部较北部,地下水流动缓慢,渗透系数整体偏小。该方法获得的渗透系数符合实际地质条件。     

Application of multiple-isotope and groundwater-age data to identify factors affecting the extent of denitrification in a shallow aquifer near a river in South Korea

The extent of denitrification in a small agricultural area near a river in Yangpyeong, South Korea, was determined using multiple isotopes, groundwater age, and physicochemical data for groundwater. The shallow groundwater at one monitoring site had high concentrations of NO₃-N (74–83 mg L^?1). The δ¹⁵N-NO₃ values for groundwater in the study area ranged between +9.1 and +24.6‰ in June 2014 and +12.2 to +21.6‰ in October 2014. High δ¹⁵N-NO₃ values (+10.7 to +12.5‰) in both sampling periods indicated that the high concentrations of nitrate in the groundwater originated from application of organic fertilizers and manure. In the northern part of the study area, some groundwater samples showed elevated δ¹⁵N-NO₃ and δ¹⁸O-NO₃ values, which suggest that nitrate was removed from the groundwater via denitrification, with N isotope enrichment factors ranging between ?4.8 and ?7.9‰ and O isotope enrichment factors varying between ?3.8 and ?4.9‰. Similar δD and δ¹⁸O values of the surface water and groundwater in the south appear to indicate that groundwater in that area was affected by surface-water infiltration. The mean residence times (MRTs) of groundwater showed younger ages in the south (10–20 years) than in the north (20–30 years). Hence, it was concluded that denitrification processes under anaerobic conditions with longer groundwater MRT in the northern part of the study area removed considerable amounts of nitrate. This study demonstrates that multi-isotope data combined with physicochemical data and age-dating information can be effectively applied to characterize nitrate contaminant sources and attenuation processes.     

From well-test data to input to stochastic continuum models: effect of the variable support scale of the hydraulic data

The question of how well the true underlying hydraulic conductivity statistics of heterogeneous media are captured by well tests is addressed. The hydraulic conductivity value and the corresponding support volume associated with a theoretical well are correlated, causing a bias in the statistics derived from well-test analyses. Statistics derived from numerically simulated well tests are compared with the known underlying conductivity statistics and the results indicate an under-prediction by simulations at higher hydraulic conductivities. The deviation starts at about mean conductivity and can be as large as an order of magnitude, with the conductivity in the vicinity of the well defining the upper boundary. In other words, the conductivity value interpreted from the well test cannot be larger than the value that the well test first encounters. Consequently, for data in this simulation exercise, the standard deviation, if only determined for the upper range of the conductivity values, would be underestimated by a factor of 1.6–2. While this specific range is likely to depend on the scale and degree of the underlying heterogeneity as well as the duration of the test, the results should be indicative of a more general behaviour and are likely to occur in other heterogeneous data as well.

---

Resumen Se plantea la pregunta de qué tan bien son representadas en las pruebas de pozo, las estadísticas reales de conductividad hidráulica subyacente de medios heterogéneos. Son correlacionados el valor de conductividad hidráulica y el volumen de apoyo correspondiente asociado con un pozo teórico, causando una distorsión en las estadísticas derivadas del análisis de la prueba de pozo. Las estadísticas derivadas de las pruebas de pozo simuladas numéricamente son comparadas con las estadísticas de conductividad subyacente conocidas, y los resultados indican una sub-predicción por las simulaciones hechas con conductividades hidráulicas más altas. La desviación empieza casi con la conductividad media y puede ser tan grande como un orden de magnitud, con la conductividad en la vecindad del pozo definiendo el límite superior. En otras palabras, el valor de conductividad interpretado a partir de la prueba del pozo no puede ser más grande que el valor que la prueba de pozo encuentre primero. Por consiguiente, para los datos en este ejercicio de simulación, la desviación estándar, si solamente fue determinada para el rango superior de los valores de conductividad, se subestimaría en un factor de 1.6–2. Mientras es probable que este rango específico dependa de la escala y del grado de la heterogeneidad subyacente, así como de la duración de la prueba, los resultados deben ser indicativos de un comportamiento más general y son probables también de ocurrir en otros datos heterogéneos.

---

Résumé Nous posons ici la question de savoir dans quelle mesure les statistiques de la conductivité hydraulique des milieux hétérogènes pourrait être révélée par des essais de puits. La valeur de la conductivité hydraulique et le volume capté correspondant sont corrélés, créant un biais dans l’analyse des statistiques dérivées des essais de puits. Les statistiques en provenance de simulations numériques d’essais de puits sont comparées avec les statistiques de conductivités connues et les résultats indiquent une sous-évaluation par les simulations, pour les conductivités hydrauliques les plus élevées: la déviation commence à partir de la valeur moyenne de la conductivité et peut atteindre la magnitude d’un ordre de grandeur en considérant la conductivité mesurée au voisinage du puits. Autrement dit, la valeur de la conductivité interprétée via l’essais de pompage ne peut être plus importante que les premières valeurs rencontrées. Par conséquence, pour les données de cet exercice de simulation, la déviation standard sera sous-estimée d’un facteur compris entre 1.6–2 pour les valeurs les plus élevées. Tandis que l’échelle spécifique de valeurs est dépendante de l’échelle et du degré de l’hétérogénéité souterraine, de même que de la durée du test, les résultats pourraient être indicatifs d’un comportement plus général et seraient sans doute observables dans d’autres cas de données hétérogènes.

     

Impact of model complexity and multi-scale data integration on the estimation of hydrogeological parameters in a dual-porosity aquifer

This study investigates the impact of model complexity and multi-scale prior hydrogeological data on the interpretation of pumping test data in a dual-porosity aquifer (the Chalk aquifer in England, UK). In order to characterize the hydrogeological properties, different approaches ranging from a traditional analytical solution (Theis approach) to more sophisticated numerical models with automatically calibrated input parameters are applied. Comparisons of results from the different approaches show that neither traditional analytical solutions nor a numerical model assuming a homogenous and isotropic aquifer can adequately explain the observed drawdowns. A better reproduction of the observed drawdowns in all seven monitoring locations is instead achieved when medium and local-scale prior information about the vertical hydraulic conductivity (K) distribution is used to constrain the model calibration process. In particular, the integration of medium-scale vertical K variations based on flowmeter measurements lead to an improvement in the goodness-of-fit of the simulated drawdowns of about 30%. Further improvements (up to 70%) were observed when a simple upscaling approach was used to integrate small-scale K data to constrain the automatic calibration process of the numerical model. Although the analysis focuses on a specific case study, these results provide insights about the representativeness of the estimates of hydrogeological properties based on different interpretations of pumping test data, and promote the integration of multi-scale data for the characterization of heterogeneous aquifers in complex hydrogeological settings.     

Physically based modelling of shallow landslide sediment yield at a catchment scale

 A shallow landslide erosion and sediment yield component, applicable at the basin scale, has been incorporated into the physically based, spatially distributed, hydrological and sediment transport modelling system, SHETRAN. The component determines when and where landslides occur in a basin in response to time-varying rainfall and snowmelt, the volume of material eroded and released for onward transport, and the impact on basin sediment yield. Derived relationships are used to link the SHETRAN grid resolution (up to 1 km), at which the basin hydrology and final sediment yield is modelled, to a subgrid resolution (typically around 10–100 m) at which landslide occurrence and erosion is modelled. The subgrid discretization, landslide susceptibility and potential landslide impact are determined in advance using a geographic information system (GIS), with SHETRAN then providing information on temporal variation in the factors controlling landsliding. The ability to simulate landslide sediment yield is demonstrated by a hypothetical application based on a catchment in Scotland. Received: 30 October 1996 · Accepted: 25 June 1997     

Semi-automatic shallow landslide detection by the integration of airborne imagery and laser scanning data

Landslide mapping is essential for effective watershed management. In Taiwan, a typhoon or earthquake event can trigger hundreds, even thousands, of shallow landslides in mountainous watersheds. Thus, improving the efficiency of landslide mapping by means of remote sensing techniques is an important issue. This study proposes a new method that uses concurrent aerial laser scanning (ALS) data and color ortho-imagery as input data: the topographic indices of slope, surface roughness, and object height model can be derived from the ALS data and the Green–Red Vegetation Index from the ortho-images. The method first uses these topographic and spectral indices in a global, semi-automatic algorithm to separate landslide from non-landslide pixels. It then offers a region growing tool and a 3D Eraser/Painter to edit detected landslides locally. These global and local operations are designed with a user interface, which is intuitive and user-friendly. Results from four test sites in a mountainous watershed prove that the method is easy, accurate, and suitable for landslide mapping in Taiwan.     

From spatial-continuous electrical resistivity measurements to the soil hydraulic functioning at the field scale

The aim of this article is to present a strategy to interpret the hydraulic functioning of a small field area by using measurements of the soil electrical resistivity. The spatial soil electrical resistivity was recorded at a high resolution on a 2 ha area by the MultiContinous Electrical Profiling (MuCEP) device at two dates. These apparent electrical resistivity measurements were firstly interpreted in terms of local electrical resistivity by 1D inverse modelling to estimate the real resistivity of the soil. These interpreted electrical resistivity data were then transformed into soil water content values and soil water potential values by the use of independent punctual data of water content and the use of the water retention curve determined by laboratory data. Our analysis has permitted us to describe the spatial variability and temporal evolution of the hydraulic functioning at high resolution from electrical resistivity data. The interpretation of the water content and matric potential maps demonstrated that some soil hydraulic processes, such as lateral overland flow, can occur in the studied zone. They would never have been detected by local measurements of soil characteristics or by the use of the soil map.     

Hydrogeological framework and estimation of aquifer hydraulic parameters using geoelectrical data: a case study from West Iran

Detailed local geological, geophysical, and hydrogeological investigations were carried out for the alluvial aquifer in the Kangavar basin, West Iran to delineate the architecture of different subsurface geological horizons using lithologs and generated vertical electrical sounding (VES) data. An attempt has also been made to estimate aquifer transmissivity from resistivity data. Forty VESs were recorded with the Schlumberger electrode configuration in the study area; 28 of these were selected for evaluation. The maximum current electrode spacing was 400–500 m. The data obtained were interpreted by computer iterative modeling with curve matching for calibration purposes. In order to ascertain the subsurface geological framework, the general distribution of resistivity responses of the geological formations was obtained and geoelectrical sections along a number of lines were prepared. Probable aquifer horizons from these sections were identified. The transmissivity of the unconfined aquifer was computed by determining the Dar-Zarrouk parameters (longitudinal unit conductance and transverse unit resistance) and were compared with the actual field transmissivity. The results showed a direct relation between aquifer transmissivity and modified transverse resistance.     

Hydrogeologic and hydraulic characterization of aquifer and nonaquifer layers in a lateritic terrain (West Bengal,India)

The hydrogeologic and hydraulic characteristics of a lateritic terrain in West Bengal, India, were investigated. Test drilling was conducted at ten sites and grain-size distribution curves (GSDCs) were prepared for 275 geologic samples. Performance evaluation of eight grain-size-analysis (GSA) methods was carried out to estimate the hydraulic conductivity (K) of subsurface formations. Finally, the GSA results were validated against pumping-test data. The GSDCs indicated that shallow aquifer layers are coarser than the deeper aquifer layers (uniformity coefficient 0.19–11.4). Stratigraphy analysis revealed that both shallow and deep aquifers of varying thickness exist at depths 9–40 and 40–79 m, respectively. The mean K estimates by the GSA methods are 3.62–292.86 m/day for shallow aquifer layers and 0.97–209.93 m/day for the deeper aquifer layers, suggesting significant aquifer heterogeneity. Pumping-test data indicated that the deeper aquifers are leaky confined with transmissivity 122.69–693.79 m²/day, storage coefficient 1.01?×?10^?7–2.13?×?10^?4 and leakance 2.01?×?10^?7–34.56?×?10^?2 day^?1. Although the K values yielded by the GSA methods are generally larger than those obtained from the pumping tests, the Slichter, Harleman and US Bureau Reclamation (USBR) GSA methods yielded reasonable values at most of the sites (1–3 times higher than K estimates by the pumping-test method). In conclusion, more reliable aquifers exist at deeper depths that can be tapped for dependable water supply. GSA methods such as Slichter, Harleman and USBR can be used for the preliminary assessment of K in lateritic terrains in the absence of reliable field methods.     

Transformations and hydraulic captures of petrochemical contaminants in a karst-fractured aquifer

 The development of a petrochemical industry may be responsible for petrochemical contamination of karst-fractured aquifers in an urban water supply system. In the Dawu Well Field, a karst-fractured aquifer in Zibo City, in the east of China, has been seriously polluted by petrochemicals from the operation of petrochemical plants. More than 60 species of organic contaminants have been detected in the water supply wells of the Liuhang-Hougao zone, in the west part of Dawu Well Field. Investigations indicate that contaminants are transported from the petrochemical plants to the karst-fractured aquifer along karst fractures and the Jinling Fault. In the groundwater, concentrations of pollutants vary with depth. Concentrations are greater with depth because of the transmissibility of the Jinling fault. The local convective flow field has a significant influence on transportation of contaminants. Hydraulic barriers can prevent the transportation of contaminants and they can be removed by extraction. Received: 22 October 1999 · Accepted: 14 February 2000     

Denitrification coupled to pyrite oxidation and changes in groundwater quality in a shallow sandy aquifer

This study focuses on denitrification in a sandy aquifer using geochemical analyses of both sediment and groundwater, combined with groundwater age dating (³H/³He). The study sites are located underneath cultivated fields and an adjacent forested area at Oostrum, The Netherlands. Shallow groundwater in the region has high nitrate concentrations (up to 8 mM) due to intense fertilizer application. Nitrate removal from the groundwater below cultivated fields correlates with sulfate production, and the release of dissolved Fe²⁺ and pyrite-associated trace metals (e.g. As, Ni, Co and Zn). These results, and the presence of pyrite in the sediment matrix within the nitrate removal zone, indicate that denitrification coupled to pyrite oxidation is a major process in the aquifer. Significant nitrate loss coupled to sulfate production is further confirmed by comparing historical estimates of regional sulfate and nitrate loadings to age-dated groundwater sulfate and nitrate concentrations, for the period 1950-2000. However, the observed increases in sulfate concentration are about 50% lower than would be expected from complete oxidation of pyrite to sulfate, possibly due to the accumulation of intermediate oxidation state sulfur compounds, such as elemental sulfur. Pollutant concentrations (NO₃, Cl, As, Co and Ni) measured in the groundwater beneath the agricultural areas in 1996 and 2006 show systematic decreases most likely due to declining fertilizer use.     

Hydrogeology and hydrochemistry of a shallow alluvial aquifer,western Saudi Arabia

A hydrogeological and hydrochemical study was conducted on a shallow alluvial aquifer, Wadi Wajj, in western Saudi Arabia to assess the influence of protection measures on groundwater quality. The hydrochemistry was assessed up-gradient and down-gradient from potential contamination sources in the main city in dry and wet seasons prior to and after the installation of major drainage and wastewater facilities. Wadi Wajj is an unconfined aquifer where water is stored and transmitted through fractured and weathered bedrock and the overlying alluvial sediments. Natural recharge to the aquifer is about 5% of rainfall-runoff. Hydrochemistry of the aquifer shows temporal and seasonal changes as influenced by protection measures and rainfall runoff. Both groundwater and runoff showed similar chemical signature, which is mostly of chloride-sulfate-bicarbonate and sodium-calcium type. Groundwater downstream of the city, though of poorer quality than upstream, showed significant improvement after the installation of a concrete runoff tunnel and a wastewater treatment plant. Concentrations of many of the groundwater quality indicators (e.g., TDS, coliform bacteria, and nitrate) exceed US Environmental Protection Agency drinking-water standards. Heavy metal content is, however, within allowable limits by local and international standards. The chemical analyses also suggest the strong influence of stream runoff and sewage water on the groundwater quality.

---

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

---

Résumé Une étude hydrogéologique et géochimique a été menée sur l’aquifère phréatique alluviale Wadi Wajj dans l’Ouest de l’Arabie Saoudite afin d’évaluer l’influence de mesures de protection sur la qualité de l’eau souterraine. L’hydrogéochimie de l’eau a été étudiée en amont et en aval de sources potentielles de contamination dans la ville principale, pendant les saisons sèches et humides, avant et après l’installation de réseaux majeurs de drainage et d’eaux usées. L’aquifère Wadi Wajj est libre, l’eau est stockée et s’écoule dans les sédiments alluviaux et dans le socle fracturé et altéré sous-jacent. La recharge naturelle de l’aquifère représente 5% des eaux de pluie et de ruissellement. L’hydrogéochimie de l’eau de l’aquifère montre que les changements saisonniers et temporaires sont influencés par les mesures de protection et par le ruissellement des eaux pluviales. L’eau souterraine et l’eau de ruissellement ont présenté le même faciès chimique, de type bicarbonatée-sulfatée-chlorée et calco-sodique. En aval de la ville, l’eau souterraine, bien que de moins bonne qualité qu’en amont, a présenté une nette amélioration après l’installation d’un système de récupération et d’une station de traitement des eaux usées. Les concentrations de plusieurs paramètres indicateurs de la qualité de l’eau (tels que la charge totale dissoute, les coliformes, et les nitrates) dépassent les normes de potabilité de l’eau de consommation de l’agence américaine de la protection de l’environnement. Les teneurs en métaux lourds en revanche n’excèdent pas les normes locales et internationales. Les analyses chimiques indiquent aussi l’influence importante de l’écoulement par ruissellement et des eaux usées sur la qualité de l’eau souterraine.

---

Resumen Un estudio hidrogeológico e hidroquímico fue hecho en un acuífero somero de Wadi Wajj, en Arabia Saudi oeste para evaluar la influecia de medidas de protección en la calidad del agua subterránea. La hidroquímica fue evaluada gradiente-arriba y gradiente-abajo de las fuentes potenciales de contaminación de la ciudad principal durante las estaciones seca y lluviosa, antes y después de la instalación de sistemas principales de drenaje y aguas servidas. Wadi Wajj es un acuífero no-confinado donde el agua es almacenada y transmitida a través de roca fracturada y meteorizada, y los sedimentos aluviales que le sobreyacen. La recarga natural del acuífero es de cerca del 5% de la precipitación-escorrentía. La hidroquímica del acuífero muestra cambios temporales y estacionales influenciados por las medidas de protección y la escorrentía de precipitación. Ambas, agua subterránea y escorrentía mostraron composición química similar, siendo mayoritariamente de tipos cloruro-sulfato-bicarbonato y sodio-calcio. El agua subterránea aguas arriba de la ciudad, aunque de calidad más pobre que aguas abajo, mostró significante mejoría después de la instalación de un tunel de concreto para escorrentía y una planta de tratamiento de aguas servidas. Las concentraciones de muchos de los indicadores de calidad de agua subterránea (e.g., STD, coliformes, y nitrato) exceden los estándares de la Agencia de Protección Ambiental USA para agua potable. El contenido de metales pesados está, sin embargo, dentro de los límites permisibles de los estándares locales e internacionales. Los análisis químicos también sugieren la fuerte influencia de la escorrentía y aguas residuales en la calidad del agua subterránea.

     

Intrinsic groundwater vulnerability determination at the aquifer scale: a methodology coupling travel time estimation and rating methods

Groundwater vulnerability has been subject of much research due to the valuable information it provides concerning groundwater protection and exploitation potential. Up to now, most groundwater vulnerability studies adopt subjective systems of rating the various factors and subsequently, their results are often ambiguous and contradicting. Within the present study a methodology for the estimation of intrinsic groundwater vulnerability at the aquifer scale is presented. The methodology is based on travel time estimation from specified sources of pollution to the aquifer. Besides the deterministic calculation of travel times, the methodology provides a rating system for each pollution source, based on its relative severity and the estimated threat that it poses to the aquifer. Therefore, it can be regarded as a hybrid method that couples the advantages provided by the physically based methods with those of the subjective rating systems. The methodology is applied to the Neon Sidirochorion aquifer, Northeastern Greece, an overexploited aquifer where river waters, sea waters and lake waters interact, causing groundwater quality deterioration to the aquifer. The results indicated that the proposed groundwater vulnerability assessment methodology is well capturing pollution related to saltwater intrusion and agricultural activities, while it is concluded that the conceptual model is significantly affecting the vulnerability assessment results and therefore has to be previously developed.     

Assessment of small-diameter shallow wells for managed aquifer recharge at a site in southern Styria,Austria

An approach to establish the recharge component of managed aquifer recharge (MAR) has recently been proposed that uses small-diameter shallow wells installed using relatively inexpensive drilling methods such as direct push. As part of further development of that approach, a generalized procedure is presented for a technical and economic assessment of the approach’s potential in comparison to other systems. Following this procedure, the use of small-diameter wells was evaluated both experimentally and numerically for a site located in southern Styria, Austria. MAR is currently done at the site using a horizontal pipe infiltration system, and system expansion has been proposed with a target rate of 12 l/s using small-diameter wells as one possible option. A short-duration single-well field recharge experiment (recharge rate 1.3–3.5 l/s) was performed (recharge by gravity only). Numerical modeling of the injection test was used to estimate hydraulic conductivity (K). Quasi-steady-state, single-well recharge simulations for different locations, as well as a long-term transient simulation, were performed using the K value calibrated from the field injection test. Results indicate that a recharge capacity of 4.1 l/s was achievable with a maximum head rise of 0.2 m at the injection well. Finally, simulations were performed for three different well fields (4, 6 and 8 wells, respectively) designed to infiltrate a target rate of 12 l/s. The experimental and numerical assessments, supported by a cost analysis of the small-diameter wells, indicate that the small-diameter wells are a viable, cost-effective recharge approach at this and other similar sites.