Basic tutorial to learn the procedure to build, simulate and represent a MODFLOW 6 model. The tutorial shows a introduction to the model file system on steady state conditions. The model for this tutorial is implemented with the following boundary conditions: Drains, Recharge, Wells, and Constant Head. The grid is regular with a width of 50 meters and it has 30 rows and 24 columns; the model has 4 layers and a total thickness of 130 meters. The model is called "hatari01" and is inspired in the "twri" model from the MODFLOW 2005 documentation adapted to MODFLOW 6. The model defines a constant horizontal hydraulic conductivity as well as vertical conductivity. After the simulation a Python code is run on a Jupyter Notebook to create the Unstructured VTK files for the heads, water table and boundary conditions representation as 3D objects in Paraview.

MODFLOW model output representation is key to understand the groundwater flow regime, to study the interaction with surface water and depending ecosystems and to evaluate the impact to anthropogenic and climate change requirements. Until now, there has been few open source software capable of generating3D representations and those software had limited options for color scales, cross sections and other graphical tools. On the research for more options we found Paraview, a open source software for data representation designed to analyze extremely large datasets using distributed memory computing resources.

In order to represent MODFLOW output into Paraview, a VTK file for unstructured grids is needed, this VTK type is called VTU where the "U" comes from unestructured. The tutorial show the complete procedure to process a MODFLOW model output into a VTU file and the representation in Paraview.

MODFLOW computes the groundwater heads over a porous / fractured media upon a series of boundary conditions as recharge, evapotranspiration, drains, well and others on steady and transient conditions. Free and commercial software is available for the MODFLOW model construction and MODFLOW output representation. Despite the fact the capabilities of these softwares, there are some gaps in data processing and representation; isometric views, animation and custom cross sections are still difficult to achieve under the existing tools, specially on multilayered models with transient conditions over series of time steps and stress periods.

There is a particular open source software for data representation that is of our interest, it is called Paraview (paraview.org). This visual application was designed to analyze extremely large datasets using distributed memory computing resources, in fact the term "para" in Paraview comes from the parallelization of computer cores.
 

Flopy is the Python library to create, run and represent results from MODFLOW models, including support for MODFLOW 6. With Flopy, the model definition and set up of boundary conditions can be done with lines of Python code, and some Python tools of optimization and machine learning can interact with the MODFLOW models. This tutorial show a basic example of model creation, configuration and simulation in Python 3 on a Jupyter Notebook. The tutorial also shows the procedure to import the MODFLOW results done with Flopy in Model Muse.

Due to their extension and climate conditions, Tailings Storage Facilities (TSF) are in interaction with groundwater flow having a dual behaviour of discharge and recharge zones. Long term evaluation of environmental impact from tailings have to consider the potential impact from the leachates into the groundwater flow regime and the impact on water quality on the downstream channels. Numerical modeling is a strong tool for the evaluation of flows and quality on the overall impact zone from tailings dam in current and future conditions.

Simulations on mine closure last for decades, and several remediation techniques as infiltration ponds, pump and test schemes, reactive walls have to be simulated on the long term. This tutorial show the procedure to simulate the particle tracking of seepage flow from tailings and estimate the potential affected stream network. The tutorial also discuss some common problems in the particle tracking simulation. 

Modeling work was done in MODFLOW with Model Muse on steady state conditions over a 5 layer model with 101260 cells. 

MODPATH is a particle tracking model that runs on the groundwater flow models simulated on MODFLOW. Particle tracking are three dimensional and can have a sort of options as forward / backward modeling direction in between others. Model Muse is a pre and post-processor of numerical models in MODFLOW, PHAST and SUTRA, this software can simulate some MODFLOW related models as MODPATH and MT3D (particle tracking and mass transport respectively). This tutorial show the basic procedure to set up a particle tracking model with MODPATH using Model Muse over a basic MODFLOW model.

In case you missed or want to see again our webminar streamed on May 15th, 2017 , here you can find the recorded video.