Modeling of Microplastics Migration in Soil and Groundwater: Insights into Dispersion and Particle Property Effects.

Migration of microplastics (MPs) in soil-groundwater systems plays a pivotal role in determining its concentration in aquifers and future threats to the terrestrial environment, including human health. However, existing models employing an advection-dispersion equation are insufficient to incorporate the holistic mechanism of MP migration. Therefore, to bridge the gap associated with MP migration in soil-groundwater systems, a dispersion-drag force coupled model incorporating a drag force on MPs along with dispersion is developed and validated through existing laboratory and field-scale experiments. The inclusion of the MP dispersion notably increased the global maximum particle velocity (vmaxp) of MPs, resulting in a higher concentration of MPs in the aquifer, which is also established by sensitivity analysis of MP dispersion. Additionally, increasing irrigation flux and irrigation areas significantly accelerates MP migration downward from soil to deep saturated aquifers. Intriguingly, vmaxp of MPs exhibited a nonlinear relationship with MPs' sizes smaller than 20 µm reaching the highest value (=1.64 × 10-5 m/s) at a particle size of 8 µm, while a decreasing trend was identified for particle sizes ranging from 20 to 100 µm because of the hindered effect by porous media and the weaker effect of the drag force. Moreover, distinct behaviors were observed among different plastic types, with poly(vinyl chloride), characterized by the highest density, displaying the lowest vmaxp and minimal flux entering groundwater. Furthermore, the presence of a heterogeneous structure with lower hydraulic conductivity facilitated MP dispersion and promoted their migration in saturated aquifers. The findings shed light on effective strategies to mitigate the impact of MPs in aquifers, contributing valuable insights to the broader scientific fraternity.

Aquifer Characterization Using Fiber Bragg Grating Multi-Level Monitoring System.

Groundwater responses measured from multiple wells at different depths are essential for delineating the aquifer heterogeneity using hydraulic tomography (HT). In general, conducting HT requires many wells because traditional well monitoring is usually partially open at a specific depth interval or is fully penetrating. Accordingly, conducting an HT survey is typically costly and time-consuming. To tackle these issues, a new multi-level monitoring system (MLMS) for the HT survey was developed using the fiber Bragg grating (FBG) technique. This FBG MLMS could collect the depth-discrete groundwater observations from a fully penetrated 2-inch well. Three field campaigns were conducted to validate the capability of the FBG MLMS for HT surveys. The results show that the accuracy and stability of this MLMS are reliable and that FBG MLMS is beneficial for conducting an HT survey. Specifically, compared to the traditional monitoring well in an injection event, this FBG MLMS can concurrently cause an increase in the number of cross-hole tests several times and collect many more head observations than the standard methods, resulting in the observed flow fields efficiently reaching ergodic conditions and effectively improving the accuracy of the estimated hydraulic heterogeneity. Therefore, the FBG MLMS could be an alternative MLMS for efficiently and economically conducting an HT survey.

Redundant and Nonredundant Information for Model Calibration or Hydraulic Tomography.

Drawdown data from independent pumping tests have widely been used to validate the estimated hydraulic parameters from inverse modeling or hydraulic tomography (HT). Yet, the independent pumping test has not been clearly defined. Therefore, the goal of this paper is to define this independent pumping test concept, based on the redundant or nonredundant information about aquifer heterogeneity embedded in the observed heads during cross-hole pumping tests. The definition of complete, moderate redundancy and high nonredundancy of information are stipulated using cross-correlation analysis of the relationship between the head and heterogeneity. Afterward, data from numerical experiments and field sequential pumping test campaigns reinforce the concept and the definition.

Two-Dimensional Probabilistic Infiltration Analysis in a Hillslope Using First-Order Moment Approach.

A first-order moment analysis method is introduced to evaluate the pore-water pressure variability within a hillslope due to spatial variability in saturated hydraulic conductivity (Ks ) during rainfall. The influences of the variance of the natural logarithm of Ks (ln Ks ), spatial structure anisotropy of ln Ks , and normalized vertical infiltration flux (q) on the evaluations of the pore-water pressure uncertainty are investigated. Results indicate different responses of pressure head variability in the unsaturated region and the saturated region. In the unsaturated region, a larger variance of ln Ks , a higher spatial structure anisotropy, and a smaller q lead to a larger variability in pressure head, while in the saturated region, the variability in pressure head increases with the increase of variance of ln Ks , the decrease of spatial structure anisotropy, or the increase of q. These variables have great impacts on the range of fluctuation of the phreatic surface within the hillslope. The influences of these three variables on the variance of pressure head within the saturated region are greater than those within the unsaturated region, and the variance of ln Ks has the greatest impact. These results yield useful insight into the effects of heterogeneity on pressure head and uncertainty associated with predicted flow field.

Finite analytic method for modeling variably saturated flows.

This paper develops a finite analytic method (FAM) for solving the two-dimensional Richards' equation. The FAM incorporates the analytic solution in local elements to formulate the algebraic representation of the partial differential equation of unsaturated flow so as to effectively control both numerical oscillation and dispersion. The FAM model is then verified using four examples, in which the numerical solutions are compared with analytical solutions, solutions from VSAFT2, and observational data from a field experiment. These numerical experiments show that the method is not only accurate but also efficient, when compared with other numerical methods.

Comparison of approaches for predicting solute transport: sandbox experiments.

The main purpose of this paper was to compare three approaches for predicting solute transport. The approaches include: (1) an effective parameter/macrodispersion approach (Gelhar and Axness 1983); (2) a heterogeneous approach using ordinary kriging based on core samples; and (3) a heterogeneous approach based on hydraulic tomography. We conducted our comparison in a heterogeneous sandbox aquifer. The aquifer was first characterized by taking 48 core samples to obtain local-scale hydraulic conductivity (K). The spatial statistics of these K values were then used to calculate the effective parameters. These K values and their statistics were also used for kriging to obtain a heterogeneous K field. In parallel, we performed a hydraulic tomography survey using hydraulic tests conducted in a dipole fashion with the drawdown data analyzed using the sequential successive linear estimator code (Yeh and Liu 2000) to obtain a K distribution (or K tomogram). The effective parameters and the heterogeneous K fields from kriging and hydraulic tomography were used in forward simulations of a dipole conservative tracer test. The simulated and observed breakthrough curves and their temporal moments were compared. Results show an improvement in predictions of drawdown behavior and tracer transport when the K tomogram from hydraulic tomography was used. This suggests that the high-resolution prediction of solute transport is possible without collecting a large number of small-scale samples to estimate flow and transport properties that are costly to obtain at the field scale.

Cost-effective hydraulic tomography surveys for predicting flow and transport in heterogeneous aquifers.

This study shows how a cost-effective hydraulic tomography survey (HTS) and the associated data estimator can be used to characterize flow and transport in heterogeneous aquifers. The HTS is an improved field hydraulic test that accounts for responses of hydraulic stresses caused by pumping or injection events at different locations of an aquifer. A sequential data assimilation method based on a cokriging algorithm is then used to map the aquifer hydraulic conductivity (K). This study uses a synthetic two-dimensional aquifer to assess the accuracy of predicted concentration breakthrough curves (BTCs) on the basis of the Kfields estimated by geometric mean, kriging, and HTS. Such Kfields represent different degrees of flow resolutions as compared with the synthetically generated one. Without intensive experimentsto calibrate accurate dispersivities at sites, the flow field based on the HTS Kfield can yield accurate predictions of BTC peaks and phases. On the basis of calculating mean absolute and square errors for estimated K fields, numerical assessments on the HTS operation strategy show that more pumping events will generally lead to more accurate estimations of Kfields, and the pump locations need to be installed in high Kzones to maximize the delivery of head information from pumps to measurement points. Additionally, the appropriate distances of installed wells are suggested to be less than one-third of the ln(K) correlation length in x direction.

Basin-scale transmissivity and storativity estimation using hydraulic tomography.

While tomographic inversion has been successfully applied to laboratory- and field-scale tests, here we address the new issue of scale that arises when extending the method to a basin. Specifically, we apply the hydraulic tomography (HT) concept to jointly interpret four multiwell aquifer tests in a synthetic basin to illustrate the superiority of this approach to a more traditional Theis analysis of the same tests. Transmissivity and storativity are estimated for each element of a regional numerical model using the geostatistically based sequential successive linear estimator (SSLE) inverse solution method. We find that HT inversion is an effective strategy for incorporating data from potentially disparate aquifer tests into a basin-wide aquifer property estimate. The robustness of the SSLE algorithm is investigated by considering the effects of noisy observations, changing the variance of the true aquifer parameters, and supplying incorrect initial and boundary conditions to the inverse model. Ground water flow velocities and total confined storage are used as metrics to compare true and estimated parameter fields; they quantify the effectiveness of HT and SSLE compared to a Theis solution methodology. We discuss alternative software that can be used for implementing tomography inversion.