The release and retention of in-situ colloids in aquifers play an important role in the sustainable operation of managed aquifer recharge (MAR) schemes. The processes of colloid release, retention, and associated permeability changes in consolidated aquifer sediments were studied by displacing native groundwater with reverse osmosis-treated (RO) water at various flow velocities. Significant amounts of colloid release occurred when: (i) the native groundwater was displaced by RO-water with a low ionic strength (IS), and (ii) the flow velocity was increased in a stepwise manner. The amount of colloid release and associated permeability reduction upon RO-water injection depended on the initial clay content of the core. The concentration of released colloids was relatively low and the permeability reduction was negligible for the core sample with a low clay content of about 1.3%. In contrast, core samples with about 6 and 7.5% clay content exhibited: (i) close to two orders of magnitude increase in effluent colloid concentration and (ii) more than 65% permeability reduction. Incremental improvement in the core permeability was achieved when the flow velocity increased, whereas a short flow interruption provided a considerable increase in the core permeability. This dependence of colloid release and permeability changes on flow velocity and colloid concentration was consistent with colloid retention and release at pore constrictions due to the mechanism of hydrodynamic bridging. A mathematical model was formulated to describe the processes of colloid release, transport, retention at pore constrictions, and subsequent permeability changes. Our experimental and modeling results indicated that only a small fraction of the in-situ colloids was released for any given change in the IS or flow velocity. Comparison of the fitted and experimentally measured effluent colloid concentrations and associated changes in the core permeability showed good agreement, indicating that the essential physics were accurately captured by the model.

[1] Conversions of compressional seismic waves to electric fields have been measured in two boreholes drilled in an unconfined sandy aquifer on the Gnangara Mound near Perth, Australia. The seismoelectric conversions at both field sites occurred in the vicinity of the water table at 13-m depth and yielded maximum amplitudes of 1 mV/m using a sledgehammer source on surface. Partially cemented layers, inferred from geological and geophysical logs, straddle the water table and may play a role in generating the conversion and influencing its amplitude distribution. The dense vertical sampling used in these borehole experiments reveals spatial and temporal polarity reversals of the interfacial signal which provide new evidence in support of the conceptual model for seismoelectric conversions at interfaces. We demonstrate that the growth rate of the source zone and its maximum vertical extent below the water table are encoded in the polarity of the interfacial signal. These experiments confirm that vertical seismoelectric profiling can be used to gain further insight into seismoelectric conversions and characteristics of interfaces that makes them amenable to detection.

Changes in subsurface temperature distribution resulting from the injection of fluids into aquifers may impact physiochemical and microbial processes as well as basin resource management strategies. We have completed a 2 year field trial in a hydrogeologically and geochemically heterogeneous aquifer below Perth, Western Australia in which highly treated wastewater was injected for large-scale groundwater replenishment. During the trial, chloride and temperature data were collected from conventional monitoring wells and by time-lapse temperature logging. We used a joint inversion of these solute tracer and temperature data to parameterize a numerical flow and multispecies transport model and to analyze the solute and heat propagation characteristics that prevailed during the trial. The simulation results illustrate that while solute transport is largely confined to the most permeable lithological units, heat transport was also affected by heat exchange with lithological units that have a much lower hydraulic conductivity. Heat transfer by heat conduction was found to significantly influence the complex temporal and spatial temperature distribution, especially with growing radial distance and in aquifer sequences with a heterogeneous hydraulic conductivity distribution. We attempted to estimate spatially varying thermal transport parameters during the data inversion to illustrate the anticipated correlations of these parameters with lithological heterogeneities, but estimates could not be uniquely determined on the basis of the collected data.

Population growth and changing climate continue to impact on the availability of natural resources. Urbanization of vulnerable coastal margins can place serious demands on shallow groundwater. Here, groundwater management requires definition of coastal hydrogeology, particularly the seawater interface. Electrical resistivity imaging (ERI) appears to be ideally suited for this purpose. We investigate challenges and drivers for successful electrical resistivity imaging with field and synthetic experiments. Two decades of seawater intrusion monitoring provide a basis for creating a geo-electrical model suitable for demonstrating the significance of acquisition and inversion parameters on resistivity imaging outcomes. A key observation is that resistivity imaging with combinations of electrode arrays that include dipole-dipole quadrupoles can be configured to illuminate consequential elements of coastal hydrogeology. We extend our analysis of ERI to include a diverse set of hydrogeological settings along more than 100 km of the coastal margin passing the city of Perth, Western Australia. Of particular importance are settings with: (1) a classic seawater wedge in an unconfined aquifer, (2) a shallow unconfined aquifer over an impermeable substrate, and (3) a shallow multi-tiered aquifer system over a conductive impermeable substrate. We also demonstrate a systematic increase in the landward extent of the seawater wedge at sites located progressively closer to the highly urbanized center of Perth. Based on field and synthetic ERI experiments from a broad range of hydrogeological settings, we tabulate current challenges and future directions for this technology. Our research contributes to resolving the globally significant challenge of managing seawater intrusion at vulnerable coastal margins.

High quality coastal aquifer systems provide vast quantities of potable groundwater for millions of people worldwide. Managing this setting has economic and environmental consequences. Specific knowledge of the dynamic relationship between fresh terrestrial groundwater discharging to the ocean and seawater intrusion is necessary. We present multi- disciplinary research that assesses the relationships between groundwater throughflow and seawater intrusion. This combines numerical simulation, geophysics, and analysis of more than 30 years of data from a seawater intrusion monitoring site. The monitoring wells are set in a shallow karstic aquifer system located along the southwest coast of Western Australia, where hundreds of gigalitres of fresh groundwater flow into the ocean annually. There is clear evidence for seawater intrusion along this coastal margin. We demonstrate how hydraulic anisotropy will impact on the landward extent of seawater for a given groundwater throughflow. Our examples show how the distance between the ocean and the seawater interface toe can shrink by over 100% after increasing the rotation angle of hydraulic conductivity anisotropy when compared to a homogeneous aquifer. We observe extreme variability in the properties of the shallow aquifer from ground penetrating radar, hand samples, and hydraulic parameters estimated from field measurements. This motived us to complete numerical experiments with sets of spatially correlated random hydraulic conductivity fields, representative of karstic aquifers. The hydraulic conductivity proximal to the zone of submarine groundwater discharge is shown to be significant in determining the overall geometry and landward extent of the seawater interface. Electrical resistivity imaging (ERI) data was acquired and assessed for its ability to recover the seawater interface. Imaging outcomes from field ERI data are compared with simulated ERI outcomes derived from transport modelling with a range of hydraulic conductivity distributions. This process allows for interpretation of the approximate geometry of the seawater interface, however recovery of an accurate resistivity distribution across the wedge and mixing zone remains challenging. We reveal extremes in groundwater velocity, particularly where fresh terrestrial groundwater discharges to the ocean, and across the seawater recirculation cell. An overarching conclusion is that conventional seawater intrusion monitoring wells may not be suitable to constrain numerical simulation of the seawater intrusion. Based on these lessons, we present future options for groundwater monitoring that are specifically designed to quantify the distribution of; (i) high vertical and horizontal pressure gradients, (ii) sharp variations in subsurface flow velocity, (iii) extremes in hydraulic properties, and (iv) rapid changes in groundwater chemistry. These extremes in parameter distribution are common in karstic aquifer systems at the transition from land to ocean. Our research provides new insights into the behaviour of groundwater in dynamic, densely populated, and ecologically sensitive coastal environments found worldwide.