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.

Groundwater modeling with several boundary conditions and complex hydrogeological setups require advanced tools for mesh discretizacion that ensures adequate refinement in the zone of interest while preserving a minimal cell account. Type of mesh has to be engineered in a way to preserve computational resources and represent adequately the groundwater flow regime.

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 tutorial 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.

Groundwater model creation requires a complete set of spatial data for the different hydraulic parameters, boundary conditions and other model items. Vector and raster data need to be preprocessed, converted, reprojected to fit the requirements of Model Muse. 

This tutorial covers an applied case of raster and vector data processing for a basin.

Python has awesome packages for machine learning that can be coupled to groundwater models and perform automatic calibration of hydraulic parameters. This tutorial covers the procedure to implement a neural network based on a set of parameters set with corresponding head values; the study case is on a transient pumping test model with an observation point located around 10 meters away from the pump well. The analysis predicts the calibration parameters from the spatially interpolated observed head values. 

We did this tutorial since the documentation or examples on the topic were not available. Flopy can export the model grid and model attributes to shapefile with a coordinate system of reference for the three types of discretizations of MODFLOW 6. We have done an applied example to export the grid of a discretized by vertices model (DISV) to the ESRI Shapefile model. The tutorial also show options for the final grid representation in geopandas.

Having a geological model can enhance numerical models since it allows to represent higher accuracy the potential distribution of hydraulic parameters in the horizontal and vertical direction. The process to implement/insert a geological model into a Modflow model is a challenge due to restrictions on proprietary software and spatial tools; we have done the whole procedure to insert a geological model into a MODFLOW-NWT groundwater model with scripts in Python, Pyvista and others.

The code extracts the cell centroids of the modflow model and then compares its position with the different geological units exported from Leapfrog to Vtk. Once the corresponding lithology of the cell is identified a hydraulic parameter is assigned. The tutorial works from Vtks of a geological model done in Leapfrog, but it can work with any Vtk.

This tutorial shows the complete procedure to convert a geological unit as a Leapfrog mesh (*.msh) to the VTK (Visualization Toolkit) format using Python and GemGIS. Follow these steps to seamlessly transfer your geological data for advanced visualization, analysis and comparison with other model results.

We have done an applied groundwater model for a pumping test on a confined aquifer. The wellsite is located on the Clairborne aquifer (Georgia, USA) and has two observations on different layers besides the well itself. The aquifer test has two periods of pumping and recovery that last 3 days each one and the numerical model was constructed with MODFLOW-6 based on a Python script with the FloPy package. Comparison among observed and simulated wells were done with plots in Matplotlib.