Control Groundwater Flow Research Articles

Characteristics and Mechanism of Hematite Dissolution and Release on Arsenic Migration in Heterogeneous Materials

The migration of arsenic in groundwater is influenced by the heterogeneity of the medium, and the presence of iron minerals adds complexity and uncertainty to this effect. In this study, a stratified heterogeneous sand column with an embedded hematite lens at the coarse-to-medium sand interface was designed. We introduced an arsenic-laden solution and controlled groundwater flow to investigate the spatiotemporal characteristics of arsenic migration and the impact of hematite dissolution. The results showed that the medium structure significantly influenced the arsenic migration and distribution within the lens-containing sand column. The clay layers directed the lateral migration of arsenic, and the arsenic concentrations in deeper layers were up to seven times greater than those on the surface. The extraction experiments of solid-phase arsenic revealed that the main adsorption modes on quartz sand surfaces were the specific adsorption (F2) and adsorption on weakly crystalline iron–aluminum oxides (F3), correlating to the specific and colloidal adsorption modes, respectively. Monitoring the total iron ions (Fe(aq)) revealed rapid increases within the first 14 days, reaching a maximum on day 15, and then gradually declining; these results indicate that hematite did not continuously dissolve. This study can aid in the prevention and control of arsenic contamination in groundwater.

De-risking green energy from mine waters by developing a robust hydrogeological conceptual model of the UK Geoenergy Observatory in Glasgow

Mine-water geothermal resources have potential to provide low-carbon heating and cooling in many areas; however, this potential has not been fully realised due to technical, economic and policy challenges. The UK Geoenergy Observatory (UKGEOS) in Glasgow was developed to provide an at-scale research facility designed to help de-risk mine-water geothermal usage. The limited knowledge of the hydrogeological systems altered by former mining activities is a key determinant of the long-term sustainability of water and heat abstraction/reinjection. This work presents a hydrogeological conceptual model developed using groundwater monitoring data obtained during the construction of the Observatory between 2020 and 2022, results from initial pumping tests performed in 2020, and results of hydrochemistry analysis from 25 sampling rounds collected between 2019 and 2022. The analysis of the data provides evidence of the dominant role of mine workings in controlling groundwater flow, with high intra-mine connectivity; increased fracturing in sandstones above mine workings; and limited inter-mine connectivity. Groundwater recharge is meteoric, mean residence times are >50 years, and there is a general upwards circulation from the deeper mine levels to the superficial deposits and the River Clyde. Faults play a significant role in limiting the extent of the highly transmissive mine workings, but there remains uncertainty surrounding the role of the faults in connecting different mine workings and their hydraulic behaviour in nonmined units. The conceptual model, that will be refined as new data become available, will be used to help guide monitoring and sampling programs and plan research activities in the Observatory.

Engineering properties of microcapsules self-healing geopolymer cutoff wall backfill under dry-wet cycles

Vertical cutoff wall has been employed for decades to control groundwater flow and subsurface contaminant transport. Environmental stress fluctuations induced by alternating dry and wet conditions impair the permeability and durability of the cutoff wall, leading to the dispersion of contaminants. The objective of this study is to incorporate self-healing microcapsules into existing geopolymer cutoff wall backfill (GCWB) to form self-healing geopolymer cutoff wall backfill (SHGCWB) that hold promise for better durability and performance. The in-situ polymerization method was used to develop single-walled and double-walled microcapsules. The microcapsules use sodium silicate as the healing agent encapsulated in single-walled polyurethane (PU) and double-walled polyurethane/melamine-formaldehyde (PU/MF) microcapsules. The effect of microcapsules on the unconfined compressive strength (UCS) and hydraulic conductivity of SHGCWB were elaborated. The durability and hydraulic conductivity variation of SHGCWB in dry-wet cycle was thoroughly investigated by Scanning Electron Microscopy (SEM) test for understanding the self-healing mechanism. The overall performance demonstrated the significant potential of the use microcapsules as a self-healing approach for cutoff walls.

Simulation-based approach for the optimization of ground freezing in tunneling

Artificial ground freezing is used in tunnelling for temporary ground improvement primarily to control groundwater flow and provide excavation support. The principle of ground freezing is based on freeze pipes bored into the ground where a coolant flows through the freeze pipes. Eventually, the ground freezing process converts pore water into ice by withdrawing the heat from the soil. In tunneling, artificial ground freezing is applied to form a closed arch frozen ground around the tunnel. Under the presence of high seepage flow, the formation time of the frozen arch is delayed. In some cases, it cannot be formed, leading to high construction costs or unsafe temporary frozen ground support. This study proposes a strategy for systematically reducing the freezing time through an optimal design arrangement of the freeze pipes. The strategy is based on the combination of the numerical modeling of the ground freezing process using a multifield finite element model in conjunction with a surrogate model to enable real-time prediction.

Durability of geopolymer cutoff wall backfill incorporating reactive MgO particles under dry-wet cycles

Vertical cutoff wall has been employed for decades to control groundwater flow and subsurface contaminant transport. The barrier performance of cutoff wall material may significantly deteriorate under dry-wet cycles due to precipitation, drought, and groundwater fluctuation. This study adopted active magnesium oxide (MgO) particles to enhance the durability of geopolymer cutoff wall backfill (GCWB) under dry-wet cycles. The influence of particle size and proportion of MgO on the unconfined compressive strength (UCS), hydraulic conductivity, self-healing capacity and durability of GCWB were comprehensively examined. The results indicate that the incorporation of active MgO particles hindered the geopolymerization process of GCWB to a certain extent, leading to an increase in hydraulic conductivity and a decrease in UCS. But the overall performance of GCWB was still favorable. The hydraulic conductivity of samples with MgO gradually increased to more than 1 × 10−8 m/s after at least two dry-wet cycles. In the subsequent cycles, the samples with MgO exhibited a pattern of first gradually declining then ascending in hydraulic conductivity. This is due to hydration and carbonization of the remained MgO in matrix during the dry-wet cycles, which converts into oxygen-rich carbonates, forming well-ramified networks and reducing the proportion of macropores. Thus, the durability of GCWB with MgO is significantly more favorable than that devoid of MgO. In particular, large particle size MgO can exert more persistent self-healing effect, endowing the samples with superior durability under the prolonged dry-wet cycles.

Groundwater Flow Model Part Of Sokoto-Rima Hydrological Basin, Northwestern Nigeria

The aim of this work is to model the groundwater system of Zamfara part of Sokoto-Rima hydrological basin, under three objectives which are to evaluate the water budget, assess the interaction between surface water and groundwater, and to predict the temporal and spatial distribution of groundwater head and groundwater flow. MODFLOW-NWT, which is a Newtonian formulation for MODFLOW-2005, was used in the study. The result of sensitivity analysis revealed groundwater (RCH), evapotranspiration (EVT) and hydraulic conductivity (HKZ) had the highest composite scales. The water budget of the calibrated model showed that groundwater recharge was 5,571,293 m3 d-1, contributing 65.7% of total inflow (8,478,903 m3 d-1) to the aquifer. Inflow from river seepage represented 30.7% of the total inflow to the aquifer by 2,597,995 m3 d-1 while seepage from the general head boundary contributed 3.7% of the inflow to the aquifer by 309,615 m3 d-1. The outflow quantification of the aquifer showed that 75.8% (6,428,824 m3 d-1) of the total outflow was accounted to by evapotranspiration, the remaining 15 % (1,273,747 m3 d-1) outflow represented the river recharge, 8.6% (730,656 m3 d-1) to general head boundary (GBH), and 0.5% (45,678 m3 d-1) by pumping wells from the aquifer. Simulated groundwater level ranges between 202.9 m. asl and 688.5 m. asl, with an average level of 414.14 m. asl. Groundwater Flow model results indicated that the topography and geologic structures control groundwater flow in the study area and that base flow to river is an important factor moderating groundwater movement. This implies that the study area currently has sufficient groundwater resources to meet the demand, despite its fragile climate condition.

The impact of hydraulic conductivity anisotropy on the effectiveness of subsurface dam

Hydraulic conductivity is a key parameter controlling groundwater flow and solute transport, and anisotropy in hydraulic conductivity is common in aquifers. Previous studies have investigated the effectiveness of subsurface dams in combating seawater intrusion (SWI) and associated environmental concerns. However, the quantitative impact of hydraulic conductivity anisotropy on the performance of subsurface dam has not been explored. This study goes beyond previous research by examining, for the first time, the effects of subsurface dams on SWI prevention in anisotropic settings. Furthermore, the study optimized subsurface dam designs in anisotropic coastal aquifers through numerical simulations that take into account the environmental impact. Our findings reveal the mechanism of hydraulic conductivity anisotropy's influence on the discharging zone and demonstrate that as the anisotropy ratio (rk) decreases, the coastal discharging zone expands landward, reducing SWI and ultimately lowering the minimum effective dam height (MEDH). This process effectively prevents SWI while maximizing fresh groundwater discharge. Increasing the rk values simultaneously expands both the MEDH and the peak fresh groundwater discharge (PFGD). Additionally, increasing the horizontal hydraulic conductivity of the aquifer significantly increases the range of PFGD, while the distance between the dam and shoreline emerges as the main factor influencing the extent of the MEDH.

Hydro-mechanical behavior and microstructure evolution of red clay-bentonite backfills

Soil-bentonite (SB) cutoff walls have been extensively used as vertical engineering barriers to control groundwater flow and subsurface contaminant migration. As the pollution containment capability of the SB backfills is intimately related to the hydro-mechanical behavior, consolidation tests and hydraulic conductivity tests were carried out to investigate the effects of bentonite content on the hydro-mechanical behavior of red clay-bentonite backfills, and the microstructure evolution of backfills under vertical loading was evaluated based on mercury intrusion porosimetry (MIP) tests and field emission scanning electron microscope (FESEM) observations. The results of microstructure analysis also evaluated the influence of bentonite on the clay fabric which changed the relationship between hydraulic conductivity and void ratio. The results indicated that the void ratio and hydraulic conductivity decreased linearly with the increase of the effective consolidation stresses, and the compression index decreased with the increase of bentonite content. Meanwhile, the pore size distribution (PSD) curve of the backfill specimen almost turned into unimodal characteristics with only micropores distributed after consolidation. Additionally, the flow paths became more tortuous under compression which decreased the hydraulic conductivity. The montmorillonite aggregates in red clay-bentonite backfill formed into tortuous network structure which bound the clay aggregates together and effectively increased the tortuosity of flow paths.

Towards parameter estimation in global hydrological models

Global hydrological models (GHMs) supply key information for stakeholders and policymakers simulating past, present and future water cycles. Inaccuracy in GHM simulations, i.e. simulation results that poorly match observations, leads to uncertainty that hinders valuable decision support. Improved parameter estimation is one key to more accurate simulations of global models. Here, we introduce an efficient and transparent way to understand the parameter control of GHMs to advance parameter estimation using global sensitivity analysis (GSA). In our analysis, we use the GHM WaterGAP3 and find that the most influential parameters in 50% of 347 basins worldwide are model parameters that have traditionally not been included when calibrating this model. Parameter importance varies in space and between metrics. For example, a parameter that controls groundwater flow velocity is influential on signatures related to the flow duration curve but not on traditional statistical metrics. Parameters linked to evapotranspiration and high flows exhibit unexpected behaviour, i.e. a parameter defining potential evapotranspiration influences high flows more than other parameters we would have expected to be relevant. This unexpected behaviour suggests that the model structure could be improved. We also find that basin attributes explain the spatial variability of parameter importance better than Köppen–Geiger climate zones. Overall, our results demonstrate that GSA can effectively inform parameter estimation in GHMs and guide the improvement of the model structure. Thus, using GSA to advance parameter estimation supports more accurate simulations of the global water cycle and more robust information for stakeholders and policymakers.

Predicting Subsurface Architecture From Surface Channel Networks in the Bengal Delta

AbstractGroundwater is the primary source of water in the Bengal Delta but contamination threatens this vital resource. In deltaic environments, heterogeneous sedimentary architecture controls groundwater flow; therefore, characterizing subsurface structure is a critical step in predicting groundwater contamination. Here, we show that surface information can improve the characterization of the nature and geometry of subsurface features, thus improving the predictions of groundwater flow. We selected three locations in the Bengal Delta with distinct surface river network characteristics—the lower delta with straighter tidal channels, the mid‐delta with meandering and braided channels, and the inactive delta with transitional sinuous channels. We used surface information, including channel widths, depths, and sinuosity, to create models of the subsurface with object‐based geostatistical simulations. We collected an extensive set of lithologic data and filled in gaps with newly drilled boreholes. Our results show that densely distributed lithologic data from active lower and mid‐delta are consistent with the object‐based models generated from surface information. In the inactive delta, metrics from object‐based models derived from surface geometries are not consistent with subsurface data. We further simulated groundwater flow and solute transport through the object‐based models and compared these with simulated flow through lithologic models based only on variograms. Substantial differences in flow and transport through the different geologic models show that geometric structure derived from surface information strongly influences groundwater flow and solute transport. Land surface features in active deltas are therefore a valuable source of information for improving the evaluation of groundwater vulnerability to contamination.

Revisiting MODFLOW's Capability to Model Flow Through Sedimentary Structures.

Sedimentary structures have unique geometries and anisotropic hydraulic conductivity, both of which control groundwater flow. Traditional finite-difference simulators (e.g., MODFLOW) have not been able to correctly represent irregular, dipping and anisotropic structures due their use of a simplified conductivity tensor, causing many modelers to turn toward finite-element codes with their sophisticated meshing capabilities. However, the release of MODFLOW 6 with its flexible discretization and multipoint flux approximation scheme prompts us to revisit its capability to compute flow through complex sedimentary structures. Through the use of a novel benchmark and case study, we show that when versions previous to MODFLOW 6 are applied to dipping structures, modeled fluxes and hence flow through the system, can be significantly over or underestimated. For example, effective conductivity for a 30° dipping layer with a 100:1 conductivity ratio is reduced to only 2% of its inputted value. We show that MODFLOW 6, with its XT3D capability and flexible discretization options is far superior to its predecessors, allowing flow through complex sedimentary structures to be simulated more accurately. However, on vertically offset grids, which have been available in all versions of MODFLOW and are often used in practice, loss of accuracy is still a concern when the vertical offset is large, that is, the dip of the sedimentary layer is steep, particularly if the layer is much more conductive than the surrounding material. The hypothesis that vertically offset grids lack sufficient hydraulic connectivity between adjacent model layers to accurately simulate the steeply dipping, highly heterogeneous case is a topic for further investigation.

A Stepwise Modelling Approach to Identifying Structural Features That Control Groundwater Flow in a Folded Carbonate Aquifer System

This paper concerns a stepwise modelling procedure for groundwater flow simulation in a folded and faulted, multilayer carbonate aquifer, which constitutes a source of good quality water for human consumption in the Apennine Range in Central Italy. A perennial river acts as the main natural drain for groundwater while sustaining valuable water-related ecosystems. The spatial distribution of recharge was estimated using the Thornthwaite–Mather method on 60 years of climate data. The system was conceptualized as three main aquifers separated by two locally discontinuous aquitards. Three numerical models were implemented by gradually adding complexity to the model grid: single layer (2D), three layers (quasi-3D) and five layers (fully 3D), using an equivalent porous medium approach, in order to find the best solution with a parsimonious model setting. To overcome dry-cell problems in the fully 3D model, the Newton–Raphson formulation for MODFLOW-2005 was invoked. The calibration results show that a fully 3D model was required to match the observed distribution of aquifer outflow to the river baseflow. The numerical model demonstrated the major impact of folded and faulted geological structures on controlling the flow dynamics in terms of flow direction, water heads and the spatial distribution of the outflows to the river and springs.

The Worth of Stochastic Inversion for Identifying Connectivity by Means of a Long‐Lasting Large‐Scale Hydraulic Test: The Salar de Atacama Case Study

AbstractUnderstanding groundwater flow involves characterizing the spatial variability of aquifer attributes and possible hydraulic connectivity structures. The latter are of crucial importance because high permeability channels may control groundwater flow and contaminant transport. In evaporitic aquifer systems, these preferential channels can consist of karst conduits, developed at different scales, as well as fault zones. These features condition the economic development of salt flats at the Central Andes. Hydraulic connectivity may affect exploitation efficiency by enhancing the inflow of less evaporated brine thus diminishing its mineral content (i.e., concentration of Li and K). In this context, we investigate if it is possible to use (tomographic) stochastic inversion (regularized pilot point method) in order to characterize the presence of connectivity structures in an evaporitic aquifer of great extension (some 1,500 km2) from head response measured at numerous observation points during a year long sequence of three hydraulic tests. Results show that, even though the solution is nonunique, the main preferential flow zones are identified. Numerous inversions yield similar fits to observed drawdowns with maximum errors of few centimeters. Preferential flow is identified not only by elongated high permeability regions, but also by a marked scale effect (model transmissivities are some 30 times larger than their local tests counterparts). The main high conductivity zones are consistent with independent information based on geophysics, isotopes, mixing ratios, piezometric data, and the expected dissolution processes.

An integrated inversion framework for heterogeneous aquifer structure identification with single-sample generative adversarial network

Generating reasonable heterogeneous aquifer structures is essential for understanding the physicochemical processes controlling groundwater flow and solute transport better. The inversion process of aquifer structure identification is usually time-consuming. This study develops an integrated inversion framework, which combines the geological single-sample generative adversarial network (GeoSinGAN), the deep octave convolution dense residual network (DOCRN), and the iterative local updating ensemble smoother (ILUES), named GeoSinGAN-DOCRN-ILUES, for more efficiently generating heterogeneous aquifer structures. The performance of the integrated framework is illustrated by two synthetic contaminant experiments. We show that GeoSinGAN can generate heterogeneous aquifer structures with geostatistical characteristics similar to those of the training sample, while its training time is at least 10 times faster than that of typical approaches (e.g., multi-sample-based GAN). The octave convolution layer and multi-residual connection enable the DOCRN to map the heterogeneity structures to the state variable fields (e.g., hydraulic head, concentration distributions) while reducing the computational cost. The results show that the integrated inversion framework of GeoSinGAN and DOCRN can effectively and reasonably generate the heterogeneous aquifer structures.

Vertical Barriers for Land Contamination Containment: A Review.

Soil pollution is one of the major threats to the environment and jeopardizes the provision of key soil ecosystem services. Vertical barriers, including slurry trench walls and walls constructed with soil mix technology, have been employed for decades to control groundwater flow and subsurface contaminant transport. This paper comprehensively reviewed and assessed the typical materials and mechanical and permeability properties of soil–bentonite, cement–bentonite and soil mix barriers, with the values of mix design and engineering properties summarized and compared. In addition, the damage and durability of barrier materials under mechanical, chemical, and environmental stresses were discussed. A number of landmark remediation projects were documented to demonstrate the effectiveness of the use of barrier systems. Recent research about crack-resistant and self-healing barrier materials incorporating polymers and minerals at Cambridge University and performance monitoring techniques were analyzed. Future work should focus on two main areas: the use of geophysical methods for non-destructive monitoring and the optimization of resilient barrier materials.

Linking the Surface and Subsurface in River Deltas—Part 2: Relating Subsurface Geometry to Groundwater Flow Behavior

AbstractUnderstanding subsurface structure and groundwater flow in deltaic aquifers is essential for evaluating the vulnerability of groundwater resources in delta systems. Deltaic aquifers contain coarse‐grained paleochannels that preserve a record of former surface river channels as well as fine‐grained floodplain deposits. The distribution of these deposits and how they are interconnected control groundwater flow and contaminant transport. In this work, we link depositional environments of deltaic aquifers to stratigraphic (static) and flow and transport (dynamic) connectivity metrics. Numerical models of deltaic stratigraphy were generated using a reduced‐complexity numerical model (DeltaRCM) with different input sand fractions (ISF) and rates of sea‐level rise (SLR). The groundwater flow and advective transport behavior of these deltas were simulated using MODFLOW and MODPATH. By comparing the static and dynamic metrics calculated from these numerical models, we show that groundwater behavior can be predicted by particular aspects of the subsurface architecture, and that horizontal and vertical connectivity display different characteristics. We also evaluate relationships between connectivity metrics and two environmental controls on delta evolution: ISF and SLR rate. The results show that geologic setting strongly influences both static and dynamic connectivity in different directions. These results provide insights into quantitatively differentiated subsurface hydraulic behavior between deltas formed under different external forcing (ISF and SLR rate) and they are a potential link in using information from delta surface networks and depositional history to predict vulnerability to aquifer contamination.

Fault zone hydrogeology in arid environments: The origin of cold springs in the Wadi Araba Basin, Egypt

AbstractThe role of faults in controlling groundwater flow in the Sahara and most of the hyper‐arid deserts is poorly understood due to scarcity of hydrological data. The Wadi Araba Basin (WAB), in the Eastern Sahara, is highly affected by folds and faults associated with Senonian tectonics and Paleogene rifting. Using the WAB as a test site, satellite imagery, aeromagnetic maps, field observations, isotopic and geochemical data were examined to unravel the structural control on groundwater flow dynamics in the Sahara. Analysis of satellite imagery indicated that springs occur along structurally controlled scarps. Isotopic data suggested that cold springs in the WAB showed a striking similarity with the Sinai Nubian aquifer system (NAS) water and the thermal springs along the Gulf of Suez (e.g., δ18O = −8.01‰ to −5.24‰ and δD = −53.09‰ to −31.12‰) demonstrating similar recharge sources. The findings advocated that cold springs in the WAB represent a natural discharge from a previously undefined aquifer in the Eastern Desert of Egypt rather than infiltrated precipitation over the plateaus surrounding the WAB or through hydrologic windows from deep crystalline basement flow. A complex role of the geological structures was inferred including: (1) channelling of the groundwater flow along low‐angle faults, (2) compartmentalization of the groundwater flow upslope from high‐angle faults, and (3) reduction of the depth to the main aquifer in a breached anticline setting, which resulted in cold spring discharge temperatures (13–22°C). Our findings emphasize on the complex role of faults and folds in controlling groundwater flow, which should be taken into consideration in future examination of aquifer response to climate variability in the Sahara and similar deserts worldwide.

Baseflow Age Distributions and Depth of Active Groundwater Flow in a Snow‐Dominated Mountain Headwater Basin

AbstractDeeper flows through bedrock in mountain watersheds could be important, but lack of data to characterize bedrock properties limits understanding. To address data scarcity, we combine a previously published integrated hydrologic model of a snow‐dominated, headwater basin of the Colorado River with a new method for dating baseflow age using dissolved gas tracers SF6, CFC‐113, N2, and Ar. The original flow model predicts the majority of groundwater flow through shallow alluvium (<8 m) sitting on top of less permeable bedrock. The water moves too quickly and is unable to reproduce observed SF6 concentrations. To match gas data, bedrock permeability is increased to allow a larger fraction of deeper and older groundwater flow (median 112 m). The updated hydrologic model indicates interannual variability in baseflow age (3–12 years) is controlled by the volume of seasonal interflow and tightly coupled to snow accumulation and monsoon rain. Deeper groundwater flow remains stable (11.7 ± 0.7 years) as a function mean historical recharge to bedrock hydraulic conductivity (R/K). A sensitivity analysis suggests that increasing bedrock K effectively moves this alpine basin away from its original conceptualization of hyperlocalized groundwater flow (high R/K) with groundwater age insensitive to changes in water inputs. Instead, this basin is situated close to the precipitation threshold defining recharge controlled groundwater flow conditions (low R/K) in which groundwater age increases with small reductions in precipitation. Work stresses the need to explore alternative methods characterizing bedrock properties in mountain basins to better quantify deeper groundwater flow and predict their hydrologic response to change.

A model ensemble generator to explore structural uncertainty in karst systems with unmapped conduits

Karst aquifers are characterized by high-conductivity conduits embedded in a low-conductivity fractured matrix, resulting in extreme heterogeneity and variable groundwater flow behavior. The conduit network controls groundwater flow, but is often unmapped, making it difficult to apply numerical models to predict system behavior. This paper presents a multi-model ensemble method to represent structural and conceptual uncertainty inherent in simulation of systems with limited spatial information, and to guide data collection. The study tests the new method by applying it to a well-mapped, geologically complex long-term study site: the Gottesacker alpine karst system (Austria/Germany). The ensemble generation process, linking existing tools, consists of three steps: creating 3D geologic models using GemPy (a Python package), generating multiple conduit networks constrained by the geology using the Stochastic Karst Simulator (a MATLAB script), and, finally, running multiple flow simulations through each network using the Storm Water Management Model (C-based software) to reject nonbehavioral models based on the fit of the simulated spring discharge to the observed discharge. This approach captures a diversity of plausible system configurations and behaviors using minimal initial data. The ensemble can then be used to explore the importance of hydraulic flow parameters, and to guide additional data collection. For the ensemble generated in this study, the network structure was more determinant of flow behavior than the hydraulic parameters, but multiple different structures yielded similar fits to the observed flow behavior. This suggests that while modeling multiple network structures is important, additional types of data are needed to discriminate between networks.

Integrated Geoenvironmental and Geotechnical Risk Assessment of East Port Said Region, Egypt for Regional Development

The present paper discusses the use of integrated topographical, hydrogeological, and geotechnical data as a tool to assess the suitability of the East Port Said Region (EPR) for the establishment of new major projects. The method comprises geo-environmental and geotechnical characterization of the surficial soil layer and subsurface rock layers and identification of soil behavior along the proposed area. We analyze the digital elevation model to identify the topographic patterns across the project area. The geology of EPR comprises a normal to slightly consolidated fluvio-marine deltaic deposits of Pleistocene to Holocene age with more than 1000 m thick. We carry out a statistical analysis of geotechnical parameters to group the risky effective parameters for regional development. The study correlated the depth-averaged values of the design geotechnical parameters over the entire EPR. It also integrates the geotechnical properties with the topographical and hydrogeological parameters in the Geographic Information System (GIS) platform aiming to identify the spatial distribution and classes of the integrated risks across the EPR. For each geotechnical–environmental domain, it links the averaged values of geotechnical parameters to the groundwater levels and the topographic gradient that control groundwater flow. The methodological approach enters the hydrogeological model of EPR to the GIS environment and crossing with the results of the geotechnical analysis to identify zones of higher water saturation and possible geotechnical risk. The developed integrated geoenvironmental and geotechnical model assesses the prevailed hazards among the entire EPR. The integrated model uses to identify and delineate geotechnical hazards zones and to shape the spatial distribution of the geotechnical hazards classes, based on the dominant geotechnical and geo-environmental specifications. The achieved model shows that the risky zone spread widely along with the middle and southwestern sectors, while the marginal zones spread along with the southeastern and northwestern sectors. The study recommends applying the prefabricated vertical drains or deep mixing technique to stabilize the soft soil in the middle and northern risky zones of EPR. It also recommends the use of gravel drains and a well-compacted sub-base gravel layer for shallow foundations to enhance the properties of unstable clay formation at the southern risk and marginal zones of EPR.