A conjunctive management framework for the optimal design of pumping and injection strategies to mitigate seawater intrusion.

Coastal aquifer management (CAM) considering conjunctive optimization of pumping and injection system for seawater intrusion (SI) mitigation poses significant decision-making challenges. CAM needs to pose multiple objectives and massive decision variables to explore tradeoff strategies between the conflicting resources, economic, and environmental requirements. Here, we investigate a joint artificial injection scheme for ameliorating SI by establishing an evolutionary multi-objective decision-making framework that combines simulation-optimization (S-O) modelling with a cost-benefit analysis, and demonstrate the framework on a large-scale CAM case in Baldwin County, Alabama. First, a SI numerical model, using SEAWAT, was configured to predict the vulnerable region as an SI encroachment area with the scenarios of minimum and maximum pumping capacity. As a result, a smaller number of candidate sites were selected in the SI encroachment area for implementing groundwater injection to avoid the computationally infeasible SI optimization with an inordinate number of injection related decision variables. Second, the effective S-O methodology of niched Pareto tabu search combined with a genetic algorithm (NPTSGA), which considers the moving-well option, was applied to discover optimal pumping/injection (P/I) strategies (including P/I rates and injection well locations) between three conflicting management objectives under complicated SI constraints. Third, for practical operation of the P/I schemes, a cost-benefit analysis provides judgment criteria to allow decision-makers to implement more sustainable P/I strategies to capture the different realistic preferences. The implementation of three extreme optimization solutions for the case study indicates that, compared to the initial unoptimized scheme, a maximum increase of a factor of 3 in groundwater extraction rates, a maximum reduction of 17% in extent of SI, and a maximum 82.3 million US dollars in comprehensive benefits are specifically achieved by conjunctive P/I optimization. The robustness in the decision alternatives attributed to the uncertainty in physical parameters of hydraulic conductivity was discovered through global sensitivity analysis. The proposed framework provides a decision support system for multi-objective CAM with combined pumping control and engineering measures for SI mitigation.

The Autonomy of Science as a Civilian Casualty of Economic Warfare: Inadvertent Censorship of Science Resulting from Unilateral Economic Sanctions.

Unilateral coercive international political, diplomatic, and economic sanctions are regular events of international relations and international law within the landscape of foreign affairs. However, while they may be prescribed by international law, or national legal systems, for peace and security reasons they have also been imposed for political grounds by powerful States such as the United States. The US sanctions are now targeting science, academic and university domains. When applied in this way, these sanctions violate international law, principles of human rights, ethics, the autonomy of scientific institutions, and the norm of universalism in science. All of which protect and promote scientific freedom of expression. It is vital that international and domestic law be correctly applied to uphold proper ethical standards and scientific independence in order to protect the work and the freedom of scholarship. In this way, law is the solution, rather than the problem.

World Map of Scientific Misconduct.

A comparative world map of scientific misconduct reveals that countries with the most rapid growth in scientific publications also have the highest retraction rate. To avoid polluting the scientific record further, these nations must urgently commit to enforcing research integrity among their academic communities.

Recruitment Processes in Academia: Does the Emperor Have Any Clothes?

The final outcome of promotion and recruitment processes in universities should be conventional and plausible by the members of the relevant scientific community, to affirm that the processes have been competitive and fair. The objective of this opinion letter is to make a plea for the importance of the post-auditing and quantitative assessment of the selection criteria. It is shown that for an example case the outcome of the post-audit does not look reasonable from an external point of view, at least regarding the research competency.

The effects of water table fluctuation on LNAPL deposit in highly permeable porous media: A coupled numerical and experimental study.

Light Non-Aqueous Phase Liquid (LNAPL) flow on the water table is highly mobile and is sensitive to the fluctuation of groundwater. This process is highly complex and involves the migration of three immiscible phases (i.e. water, LNAPL and air) which need the explicit definition of multiple parameters. A coupled experimental and numerical simulation methodology is performed by using Time Domain Reflectrometer (TDR) and multiphase simulation of a controlled environment to mimic the water table fluctuation and its effect on the LNAPL residual saturation. TDR probes are installed in different locations of a 2D tank (i.e. a cuboid box with relatively low off-plane thickness) and the bulk permittivity of the phases are measured through artificially imposed boundary conditions. The bulk permittivity is then translated into saturation of the three different phases. The translated residual saturations along with the previously measured porous media properties (e.g. porosity and saturated permeability) are then inserted into the numerical simulator (i.e. COMSOL Multiphysics®) and the migration of the three phase in porous media is simulated. The numerical exponents and entry pressures needed for the simulation of the multiphase flow are estimated using the temporal experimental values. The exponents of water LNAPL relative permeability were estimated to be around 2 while the exponents gas LNAPL relative permeability were estimated to be closer to 3. The results, simulated with the optimized parameters, are then evaluated with pictures taken from the transparent face of the 2D tank different stages of the experiment. The temporal evolution of different phase saturation has been compared and validated between the experimental results obtained and interpreted by the TDR probe measurements and the simulations. The relative error stays in the 5 % confidence level for most reported points and only in the highly dynamic flow time steps the error reaches around 12% which are discussed in the text and is accepted due to the highly nonlinear nature of the problem.

DNAPL flow and complex electrical resistivity evolution in saturated porous media: A coupled numerical simulation.

Induced Polarization (IP) is a non-intrusive geophysical method to monitor Dense Non-Aqueous Phase Liquid (DNAPL) contamination and remediation processes underground. In this study, an advanced numerical code simulating DNAPL flow and complex electrical resistivity is presented. The model was validated against existing IP results and image measurements that were carried out previously in a series of 2D tank experiment. Multiphase flow modeling in porous media is coupled with electrical current modeling to simulate the process of DNAPL migration and the associated IP response. This brings a broader view of the contamination in space and time compared to surface and borehole measurements, especially when the results are supported by field measurements or laboratory experiments. The simulations are developed in 3D and are performed in COMSOL Multiphysics®. The simulations using petrophysical relationships for in-phase and quadrature resistivity and the results of the experiments are in complete accordance with each other in the parts of the tank where the saturation of DNAPL is relatively low (i.e. especially in the cone of depression in the pumping scenario). However, the parts associated with high saturation of DNAPL show high errors between the in-phase resistivity simulations and the results from experiments. The present work can be regarded as a preliminary study toward further applications of coupled IP-multiphase flow for more accurate detection and monitoring of DNAPLs. It is suggested that the choice of tool/approach in this study be extended to larger-scale studies for further investigation.

Numerical simulations of high viscosity DNAPL recovery in highly permeable porous media under isothermal and non-isothermal conditions.

We developed a decimetric size model based on coupling generalized Darcy's law and heat-transfer equations to model viscous dense non-aqueous phase liquid (DNAPL) pumping through highly permeable porous media under non-isothermal conditions. The presence of fingering and non-wetting phase ganglia was modeled through an unsteady capillary diffusion coefficient and an arbitrary heterogeneous permeability field. The model was validated using existing experimental data of a simple case, an oil injection in a 2D tank packed with glass beads. Next, we compared the results of this model against a DNAPL extracting situation in the 2D tank to better understand the two-phase flow behavior in highly permeable porous media. We found that natural convection during heating plays an essential role in heat transfer, especially in the wetting phase zone. By adding the dynamic effect (unsteady conditions) we were better able to describe the presence of the ganglia in porous media. We observed good agreement between modeled and experimental oil saturation curves until the breakthrough point, with a mean relative error of about 10% for low and high flow rates, and 8% and 16% after breakthrough for low and high flow rates, respectively. Extracting viscous oil at low flow rates and high temperature generates less fingering and is well described by the generalized Darcy's law. The remobilization of residual non-wetting ganglia after the breakthrough point at the outlet is, however, difficult to simulate using the generalized Darcy's law. In the end, we treated this issue by using a perturbed permeability field to simulate the observed fingering in the 2D tank.

Land subsidence: A global challenge.

This study presents a comprehensive review of the Land subsidence (LS) cases, as a worldwide environmental, geological, and global geohazard concern. Here, 290 case studies around the world mostly conducted in large metropolitan cities (e.g. Bangkok, Beijing, California, Houston, Mexico City, Shanghai, Jakarta, and Tokyo) in 41 countries were collected. The spatial distribution of LS characteristics (e.g. intensity, magnitude, and affected area), impacts, and influential factors are scrutinized. Worldwide attempts to remedy the crisis of LS were also investigated in this review. It is shown that the coastal plains and river deltaic regions are of high-frequent subsided areas around the world (~47% of 290 study areas). The spaceborne monitoring of LS is the more prevalent technique (~ 38% of total cases) compared to the ground-investigation (e.g. geological surveying, leveling, GPS, and modeling). Human-induced LS cases are 76.92% of all the LS cases around the world and groundwater extraction contributes 59.75% of these cases. Strong direct correlations with the exponential trend are observed between the average LS rate (LSavg) with groundwater withdrawal (R2 = 0.950) and groundwater level decline (R2 = 0.888). To understand the influential factors on LS occurrences, the relationship of LS rate with climate factors, hydrogeological characteristics of the aquifer, human-induced factors are investigated. Finally, we provide future research guidelines and implications that need to be expanded in order to better monitor and reduce the impact of the LS phenomenon. The outcomes of this study can be used to derive a framework helpful for interpreting the observed LS phenomena and for forecasting future situations to mitigate or control this geohazard.

Interaction of lake-groundwater levels using cross-correlation analysis: A case study of Lake Urmia Basin, Iran.

Lake Urmia (LU) is the second largest hypersaline lake in the world. Lake Urmia's water level has dropped drastically from 1277.85 m to 1270.08 m a.s.l (equal to 7.77 m) during the last 20 years, equivalent to a loss of 70% of the lake area. The likelihood of lake-groundwater connection on the basin-scale is uncertain and understudied because of lack of basic data and precise information required for physically-based modeling. In this study, cross-correlation analysis is applied on a various time-frames of water level of the lake and groundwater levels (2001-2018) recorded in 797 observation wells across 17 adjacent aquifers. This provides insightful information on the lake-groundwater interaction. The cross-correlation coefficient between the monthly water level of lake and observations wells (rGW-L) and the difference of these two variables (Hf) was calculated for different time-frames. The values of rGW-L (ranged -0.69 to 0.97) and Hf (ranged -53 m to 293 m) indicated the significant role of time-frames of observed dataset on dynamic behavior of lake-groundwater interaction, and exchange fluxes in the study setting. Results suggested two opposing behaviors in lake-groundwater interaction of the study system mainly arise from anthropogenic activity (overexploitation of groundwater for irrigation) and aquifer type (unconfined/pressurized): three out of 17 adjacent aquifers are feeding by the LU and act as "gaining aquifers" (located in northern half of LU) and others discharging into the LU and act as "losing aquifers". This study aimed to provide easy-to-obtain insights into LGWI in the complex setting of LU Basin. It can be considered a preliminary step towards a deeper understanding of the interaction through physically-based analysis and modeling.

Normalized difference vegetation index as the dominant predicting factor of groundwater recharge in phreatic aquifers: case studies across Iran.

The estimation of long-term groundwater recharge rate ([Formula: see text]) is a pre-requisite for efficient management of groundwater resources, especially for arid and semi-arid regions. Precise estimation of [Formula: see text] is probably the most difficult factor of all measurements in the evaluation of GW resources, particularly in semi-arid regions in which the recharge rate is typically small and/or regions with scarce hydrogeological data. The main objective of this study is to find and assess the predicting factors of [Formula: see text] at an aquifer scale. For this purpose, 325 Iran's phreatic aquifers (61% of Iran's aquifers) were selected based on the data availability and the effect of eight predicting factors were assessed on [Formula: see text] estimation. The predicting factors considered include Normalized Difference Vegetation Index (NDVI), mean annual temperature ([Formula: see text]), the ratio of precipitation to potential evapotranspiration ([Formula: see text]), drainage density ([Formula: see text]), mean annual specific discharge ([Formula: see text]), Mean Slope ([Formula: see text]), Soil Moisture ([Formula: see text]), and population density ([Formula: see text]). The local and global Moran's I index, geographically weighted regression (GWR), and two-step cluster analysis served to support the spatial analysis of the results. The eight predicting factors considered are positively correlated to [Formula: see text] and the NDVI has the greatest influence followed by the [Formula: see text] and [Formula: see text]. In the regression model, NDVI solely explained 71% of the variation in [Formula: see text], while other drivers have only a minor modification (3.6%). The results of this study provide new insight into the complex interrelationship between [Formula: see text] and vegetation density indicated by the NDVI. The findings of this study can help in better estimation of [Formula: see text] especially for the phreatic aquifers that the hydrogeological ground-data requisite for establishing models are scarce.

Assessment of sustainable groundwater resources management using integrated environmental index: Case studies across Iran.

Assessing environmentally sustainable GW management (ESGM) needs a deep knowledge of the present and the projected status of GW (GW) quantity and quality. Translations of these data into policy relevant information are usually done through quantitative indices. Despite the availability of a dozen GW sustainability indicators, defining an integrated index based on internationally accepted scientific standards indicators is required. To fill this gap, an in-depth review on the developed indicators/index for evaluation of GW sustainable management (GWSM) from an environmental viewpoint at aquifer scales is provided in this study. Thirteen environmentally related quantitative indicators are adopted for assessment of GWSM, especially in arid regions, depending upon data availability, and relevance of indicators. An integrated ESGM index (ESGMI) is developed based on weighted aggregation of thirteen adopted indicators through multi criteria decision making (MCDM) methods. ESGMI value ranged between 0 and 100, zero value denotes to the worst state or unsustainable GW management (GWM) and 100 indicates the ideal state or GWM is sustainable. Thirty important aquifers across Iran are chosen to implement the ESGMI at the national scale of a country known to be the fifth largest global GW user. ESGMI values for thirty of Iran's aquifers are obtained in the range 15.40 to 68.50 (on average, 49.96). This reveals the unsustainable status of GWM in this country. The results of this study demonstrate that the ESGMI is a promising tool to determine the current state of GW quantity and quality, reveals the effect of policy actions and plans, and contributes to the development and operation of effective sustainable management policies for GW resources. Due to uncertainties and spatio-temporal variabilities of key controlling variables in GW management, sustainability evaluation should be understood as a dynamic and iterative process, requiring persistent monitoring, analysis, prioritization, and modification.

Non-pumping reactive wells filled with mixing nano and micro zero-valent iron for nitrate removal from groundwater: Vertical, horizontal, and slanted wells.

Non-pumping reactive wells (NPRWs) filled by zero-valent iron (ZVI) can be utilized for the remediation of groundwater contamination of deep aquifers. The efficiency of NPRWs mainly depends on the hydraulic contact time (HCT) of the pollutant with the reactive materials, the extent of the well capture zone (Wcz), and the relative hydraulic conductivity of aquifer and reactive material (Kr). We investigated nitrate removal from groundwater using NPRWs filled by ZVI (in nano and micro scales) and examined the effect of NPRWs orientations (i.e. vertical, slanted, and horizontal) on HCT and Wcz. The dependence of HCT on Wcz for different Kr values was derived theoretically for a homogeneous and isotropic aquifer, and verified using particle tracking simulations performed using the semi-analytical particle tracking and pathlines model (PMPATH). Nine batch experiments were then performed to investigate the impact of mixed nano-ZVI, NZVI (0 to 2 g l-1) and micro-ZVI, MZVI (0 to 4 g l-1) on the nitrate removal rate (with initial [Formula: see text] =132 mg l-1). The NPRWs system was tested in a bench-scale sand medium (60 cm length × 40 cm width × 25 cm height) for three orientations of NPRWs (vertical, horizontal, and slanted with inclination angle of 45°). A mixture of nano/micro ZVI, was used, applying constant conditions of pore water velocity (0.024 mm s-1) and initial nitrate concentration (128 mg l-1) for five pore volumes. The results of the batch tests showed that mixing nano and micro Fe0 outperforms these individual materials in nitrate removal rates. The final products of nitrate degradation in both batch and bench-scale experiments were [Formula: see text] , [Formula: see text] , and N2(gas). The results of sand-box experiments indicated that the slanted NPRWs have a higher nitrate reduction rate (57%) in comparison with vertical (38%) and horizontal (41%) configurations. The results also demonstrated that three factors have pivotal roles in expected HCT and Wcz, namely the contrast between the hydraulic conductivity of aquifer and reactive materials within the wells, the mass of Fe0 in the NPRWs, and the orientation of NPRWs adopted. A trade-off between these factors should be considered to increase the efficiency of remediation using the NPRWs system.

Conceptualization of Karstic Aquifer with Multiple Outlets Using a Dual Porosity Model.

In this study, two conceptual models, the classic reservoir (CR) model and exchange reservoirs model embedded by dual porosity approach (DPR) are developed for simulation of karst aquifer functioning drained by multiple outlets. The performances of two developed models are demonstrated at a less developed karstic aquifer with three spring outlets located in Zagros Mountain in the south-west of Iran using 22-years of daily data. During the surface recharge, a production function based on water mass balance is implemented for computing the time series of surface recharge to the karst formations. The efficiency of both models has been assessed for simulation of daily spring discharge during the recession and also surface recharge periods. Results indicate that both CR and DPR models are capable of simulating the ordinates of spring hydrographs which drainage less developed karstic aquifer. However, the goodness of fit criteria indicates outperformance of DPR model for simulation of total hydrograph ordinates. In addition, the DPR model is capable of quantifying hydraulic properties of two hydrologically connected overlapping continua conduits network and fissure matrix which lays important foundations for the mining operation and water resource management whereas homogeneous model representations of the karstic subsurface (e.g., the CR) do not work accurately in the karstic environment.

Effects of heterogeneities on capillary pressure-saturation-relative permeability relationships.

In theories of multiphase flow through porous media, capillary pressure-saturation and relative permeability-saturation curves are assumed to be intrinsic properties of the medium. Moreover, relative permeability is assumed to be a scalar property. However, numerous theoretical and experimental works have shown that these basic assumptions may not be valid. For example, relative permeability is known to be affected by the flow velocity (or pressure gradient) at which the measurements are carried out. In this article, it is suggested that the nonuniqueness of capillary pressure-relative permeability-saturation relationships is due to the presence of microheterogeneities within a laboratory sample. In order to investigate this hypothesis, a large number of "numerical experiments" are carried out. A numerical multiphase flow model is used to simulate the procedures that are commonly used in the laboratory for the measurement of capillary pressure and relative permeability curves. The dimensions of the simulation domain are similar to those of a typical laboratory sample (a few centimeters in each direction). Various combinations of boundary conditions and soil heterogeneity are simulated and average capillary pressure, saturation, and relative permeability for the "soil sample" are obtained. It is found that the irreducible water saturation is a function of the capillary number; the smaller the capillary number, the larger the irreducible water saturation. Both drainage and imbibition capillary pressure curves are found to be strongly affected by heterogeneities and boundary conditions. Relative permeability is also found to be affected by the boundary conditions; this is especially true about the nonaqueous phase permeability. Our results reveal that there is much need for laboratory experiments aimed at investigating the interplay of boundary conditions and microheterogeneities and their effect on capillary pressure and relative permeability.

Influence of boundary condition types on unstable density-dependent flow.

Boundary conditions are required to close the mathematical formulation of unstable density-dependent flow systems. Proper implementation of boundary conditions, for both flow and transport equations, in numerical simulation are critical. In this paper, numerical simulations using the FEFLOW model are employed to study the influence of the different boundary conditions for unstable density-dependent flow systems. A similar set up to the Elder problem is studied. It is well known that the numerical simulation results of the standard Elder problem are strongly dependent on spatial discretization. This work shows that for the cases where a solute mass flux boundary condition is employed instead of a specified concentration boundary condition at the solute source, the numerical simulation results do not vary between different convective solution modes (i.e., plume configurations) due to the spatial discretization. Also, the influence of various boundary condition types for nonsource boundaries was studied. It is shown that in addition to other factors such as spatial and temporal discretization, the forms of the solute transport equation such as divergent and convective forms as well as the type of boundary condition employed in the nonsource boundary conditions influence the convective solution mode in coarser meshes. On basis of the numerical experiments performed here, higher sensitivities regarding the numerical solution stability are observed for the Adams-Bashford/Backward Trapezoidal time integration approach in comparison to the Euler-Backward/Euler-Forward time marching approach. The results of this study emphasize the significant consequences of boundary condition choice in the numerical modeling of unstable density-dependent flow.