This is an applied case of model calibration on a regional scale taking into account 10 piezometers over 6 hydrogeological units and a drain conductance. The model has a Voronoi mesh created from the catchment extension, the river network and piezometer location. Python scripting has been optimized with utils from the mf6Voronoi package. Comparison and statistics among observed and calculated heads are finally generated and discussed.

This example develops a groundwater model in MODFLOW6 DISV that implements a Voronoi mesh generated from the basin boundary and river network with the mf6Voronoi packages that allow high performance meshing fully coupled with geospatial data. On the Voronoi mesh, the refinement levels are defined by a minimum cell size, maximum cell size and a multiplier. The applied case covers all steps on model discretization, construction, simulation and 2D visualization.

Basic example of MODFLOW 6 groundwater flow modeling with Flopy. The model is multilayered and runs on transient conditions with the recharge, well and river boundary conditions implemented. The example also load model results and create head distributions with flow direction plots with Matplotlib. The example runs entirely on the Google Colab platform with the use of a Github repository with the repositories. This workflow is intended to give a complete online experience of groundwater modeling without any particular computer/software requirement.

There is a new dataset from the European Space Agency (ESA) called Copernicous DEM (GLO-30) that has global coverage with raster resolution of 30m available to download from the Copernicous Dataspace. We have prepared an applied case to access this dataset, download and clip for an area of study and process it into Hatari Utils that provides the watershed and river network delineation on spatial format. The example also covers the options to determine the accumulation area and explore other outputs from the platform.

Elevation maps that represent surface and river bathymetry are an essential input for flood modeling in software like HEC-RAS. Even with the latest version of high profile open source GIS software as QGIS the combination of a surface elevation map and a river bottom elevation map is a challenge that requires many turnaround, conversion and manual labour.

We have developed a useful script that works with surface and river bottom elevation in a "smart" way and creates a geospatial raster with the topobathymetric elevation by the use of geospatial Python libraries as Shapely and Rasterio. The script also includes some key steps to identify the river body and treat the missing values on the bathymetry map.

We have developed an applied case of Piper diagram generation over water samples from piezometers and springs from an Environmental Impact Assessment (EIA) on a mining project. The procedure covers all steps to format concentrations and generate the Piper diagram online in Aquifer App (aquifer.hatarilabs.com). A discussion of the point distribution is performed based on the chemistry data, spatial distribution, and geology.

This webinar covers an applied case of reservoir modeling at the discharge point of an andean basin. The example covers the steps to set up the watershed,  define the reservoir network, simulate the hydrological alternative and explore output hydrographs and tables.

Construction of complex and high performance groundwater modes requires an appropriate distribution of the hydrogeological units with depth. We have developed an applied case of 3D modeling of hydrogeological units from data provided as raster (*.tiff) files. The webinar all steps from orientation and elevation sampling from raster data, the conceptualization and setup of the Gempy model in Aquifer App and finally the postprocessing of the 3D lithology based on the surface elevation. The input data comes from the USGS research on the aquifer Systems in the Williston and Powder River Structural Basins, United States and Canada.

We have developed an applied case of groundwater model discretized from ESRI Shapefiles with refinement areas. Boundary conditions are also set up from spatial data with the intersect functionality of Flopy. Model surface and layer bottom are imported / processed from xyz point data. The simulation is run for one steady  and ten transient stress periods and results are plotted for head on aerial view and cross section.

In order to improve the accuracy of a groundwater flow simulation we need a strong conceptual model, a high quality observation dataset and a numerical code that can handle specific characteristics of the groundwater flow we want to evaluate. MODFLOW 6 can deal with variable density flow with the BUY package and also with variable viscosity with the VSC package that is implemented on the latest versions of Model Muse. This webinar covers an applied case of variable viscosity modeling due to the injection of salinity and heat to a model with regional groundwater flow. The main features for the implementation of the VSC package are explained and the resulting concentrations on two observation points are evaluated with Flopy codes.

This webinar covers an applied case of regional groundwater modeling with a group of regional fractures that define areas of small cell size. The applied case covers the construction of a groundwater model with a refinement that ranges from 100m on the outer parts of the model to 15 meters in the area close to the faults.

The application of PEST to real conditions is a challenge since it involves the setup of the automated parameter estimation and the complexities of the groundwater response. We have developed an applied case of automatic parameter estimation on a groundwater flow model related to a site drainage with 6 pumps and steel sheet piles. The calibration is based on observations from a pumping well and two outer monitoring points over a period of 120 days. The case involves as well the sensitivity analysis of hydraulic conductivities from heterogeneous soil and sheet piles together with the conductance of neighboring water channels.

There are many topics to take into account when we create a groundwater flow model. One of the most important is to represent the available geological information into the distribution the hydraulic parameters. We have researched a simple and reproducible workflow to create a geological model from point shapefile and insert the modeled geological units into Model Muse with corresponding K values. This webinar cover the whole procedure to create a geological model with Gempy and Aquifer App and the necessary codes to create a xyz shapefile that is later imported in Model Muse with the Kx value extracted from the attribute table.

Aquifer App offers a friendly, clean and powerful way to create Gempy scripts for geological modeling. We have developed an applied case of a 3D petroleum system model for a part of the Williston Basin, USA that contains the Nesson Anticline. Top layer information was provided as raster format where elevation and orientation for random points were extracted with Python codes using the Gemgis package and exported in Gempy input file format. On Aquifer App the processed CSV files were inserted and the geological model was set up with the corresponding geological sequence. 

Lithium is present in some salt lakes along with other brines in different concentrations. In order to simulate the right groundwater flow regime of Lithium brine along other brines we need to model variable density flow from multiple species. This type of modeling is possible on the latest versions of Modflow 6 with the BUY package and Model Muse. 

We have developed an applied case of variable density flow modeling on a andean salt lake that is pumped from 3 wells in 10 years. The model has a grid refinement close to the lake and works with 6 brines: B(OH)3, NaHCO3, NaCl, KHSO4, LiHSO4 and NaHSO4. The example implements the BUY package with all the required parameters together with the configuration of the boundary conditions and initial conditions. Concentration and water balances are processed with Python to evaluate the brine migration over time and the impact of pumping on the water balance.