June 05, 2025

Representation of Flow Direction Vectors from a Model Muse / MODFLOW6 model with mf6Voronoi - Tutorial

June 05, 2025/ Saul Montoya

This tutorial show the complete workflow to determine flow directions from a MODFLOW model done with Model Muse with the Mf6Voronoi utility FlowVectorGenerator that represent flow direction vectors in a Matplotlib figure. The function works not only with normal Modflow6 Dis models like the generated by Model Muse but also with Voronoi based Modflow6 Disv models. There is a capability to insert background images and other functionalities, however the styling options are closed in order to preserve the artistic style of the plot.

May 23, 2025

Create 3D objects (Vtks) from a groundwater model in Model Muse with mf6Voronoi - Tutorial

May 23, 2025/ Saul Montoya

Model Muse is an awesome tool for groundwater model construction and output visualization, however its graphical capabilities are quite limited. Since we are always on the search of bringing groundwater model ‘close’ to the stakeholder, normal audience and coming generations, there is a special interest in providing complementary tools for the normal workflow with Model Muse and Modflow6 that allows the normal user to create awesome 3D objects representing:

May 16, 2025

Simulation of open pit inflows using MODFLOW 6 and mv6Voronoi - Tutorial

May 16, 2025/ Saul Montoya

This applied case presents a step-by-step methodology to simulate groundwater inflow into an open pit using MODFLOW 6—the groundwater flow modeling engine developed by the USGS—in combination with mv6Voronoi, a tool that generates Voronoi meshes.

May 14, 2025

Paraview and mf6Voronoi for Modflow 6 model 3D visualization - Tutorial

May 14, 2025/ Saul Montoya

Because an image is worth a thousand words, and 3D visualization is key for explaining, sharing, and analyzing groundwater model parameters, boundaries and results. The Python package mf6Voronoi offers a structured workflow to generate 3D spatial files on VTK format for the representation of the model parameters, water table, head distribution and others that can be plotted on Paraview. This is an applied case of 3D representation for a andean basin groundwater flow model where we are going to explore the relevant tools and options of mf6Voronoi and Paraview.

April 24, 2025

MODFLOW 6 Syntax plugin for Visual Studio Code - Tutorial

April 24, 2025/ Saul Montoya

There is a plugin for Visual Studio Code that improves MODFLOW 6 model readability. The installation process is quite easy and works for any model without any additional setting or command. This is an applied case of the plugin installation and exploration of a MODFLOW 6 model for a andean basin.

April 21, 2025

Underground mine inflow modeling with MODFLOW6, Flopy and Mf6Voronoi - Tutorial

April 21, 2025/ Saul Montoya

An applied example of modeling the development of underground mine works with time based on a Voronoi mesh MOFLOW 6 model. The case consists of 5 mine levels developed over 5 years that cross the underground water body at determined lengths. The numerical model is developed with Python, Flopy and Mf6Voronoi for spatial discretization, implementation of the boundary conditions based on geospatial data, model build/simulation and 3d representation on Paraview.

February 04, 2025

Basic example of groundwater modeling with MODFLOW6 and Flopy in Google Colab - Tutorial

February 04, 2025/ Saul Montoya

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 MODFLOW 6 executables. This workflow is intended to give a complete online experience of groundwater modeling without any particular computer/software requirement.

January 24, 2025

Topobathymetric elevation generation for flood modeling with geospatial Python

January 24, 2025/ Saul Montoya

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.

October 24, 2024

4 applied courses in groundwater modeling with MODFLOW - Flopy and more

October 24, 2024/ Saul Montoya

Besides doing amazing tutorials in groundwater modeling with MODFLOW and Flopy we also have the most complete set of Flopy courses applied to different areas of groundwater resources evaluation and impact assessment. But, there is more: we have a whole program that trains you from zero knowlege of Python to develop fully featured groundwater flow models.

September 25, 2024

Geological Modeling of Sedimentary Layers from Geospatial Data with Python and Aquifer App - Tutorial

September 25, 2024/ Saul Montoya

A practical example for constructing a geological model using Python and Aquifer App, based on geospatial data in ESRI Shapefile format and raster data in Tiff format. The tutorial demonstrates the complete procedure of spatial data processing using libraries such as Geopandas, Rasterio, and Pyvista to create surface and orientation files, which are then input into Aquifer App. Finally, the project is run locally in Gempy, generating the geometry of lithology and contacts.

August 15, 2024

How to create a Geological Model from a Point Shapefile with Gempy and Aquifer App - Tutorial

August 15, 2024/ Saul Montoya

This is an applied case were we build a geological model only with lithology information stored as a point shapefile. The tutorial covers all steps from raster (array) generation for all surfaces together with the orientation sampling and format of surfaces/orientations as Gempy input files. The generated data was inserted in the Aquifer App that implements an interface to create Gempy models. Finally the lithology and layer surface geometry was exported as Vtk to be represented on Paraview with the initial data to evaluate the accuracy of the simulation.

August 12, 2024

Geological modeling of a petroleum system with Gempy, Gemgis and Aquifer App - Tutorial

August 12, 2024/ Saul Montoya

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.

July 17, 2024

Export 3d Geological Units as Vtk format with Gempy and Aquifer App - Tutorial

July 17, 2024/ Saul Montoya

Aquifer App offers a friendly, clean and powerful way to create Gempy scripts for geological modeling. The platform also allows users to run the scripts and download the whole modeling project that has the 3D geometry of the geological units and faults in Vtk format. This tutorial shows an applied case of geological data visualization, model creation, model run, export and representation of the Vtk files in Paraview.

July 12, 2024

Regional geological modeling with surface topography with Gempy and Aquifer App - Tutorial

July 12, 2024/ Saul Montoya

The coupling of Python scripts with the Aquifer App platform allows us to generate geological models with a versatility not seen before. This is a basic example of regional geological modeling in Aquifer App with surface topography from a digital elevation model in TIF format. The tutorial also makes a 3D representation of the contacts and orientations of the one layer geology in a Jupyter notebook for better analysis of the geological sequence.

July 03, 2024

Geological modeling of a recumbent fold with Gempy and Aquifer App - Tutorial

July 03, 2024/ Saul Montoya

The coupling of Python scripts with the Aquifer App platform allows us to generate geological models with a versatility not seen before. This is a basic example of geological modeling of an recumbent fold in Aquifer App that can create models using Gempy. The tutorial also makes a 3D representation of the contacts and orientations of the layers in a Jupyter notebook for better analysis of the geological sequence.

July 02, 2024

Flood simulation from direct rainfall with Python and Landlab - Tutorial

July 02, 2024/ Saul Montoya

Simulation of direct precipitation events over a surface is important to assess the impacted areas and manage disaster situations. Response to precipitation is evaluated over cells in a 2D model and runoff is routed to the stream networks. Landlab is a Python library for a variety of surface processes and can model shallow water flow over topography with a variety of algorithms. We have developed an applied case of direct rainfall for flood modeling over a geospatial raster where an initial water height is applied and then flows over the surface. Resulting water heights are exported as geospatial rasters.

June 28, 2024

Geological modeling of an anticline structure with Gempy and Aquifer App - Tutorial

June 28, 2024/ Saul Montoya

The coupling of Python scripts with the Aquifer App platform allows us to generate geological models with a versatility not seen before. This is a basic example of geological modeling of an anticline structure in Aquifer App that can create models using Gempy. The tutorial also makes a 3D representation of the contacts and orientations of the layers in a Jupyter notebook for better analysis of the geological sequence.

June 21, 2024

Missing Crop Vegetation Areas Detection with Python and Scikit Learn

June 21, 2024/ Saul Montoya

An applied case for the recognition of missing crop vegetation areas based on a drone orthophoto. Contours have been identified from an enhanced combination of raster bands with a marching squares method to find constant valued contours and then exported as geospatial polygons.

June 18, 2024

Simplify / densify stream network online with Hatari Utils - Tutorial

June 18, 2024/ Saul Montoya

Watershed delineation and stream network determination are common tasks on the hydrological analysis for any area of interest. There is a parameter on the stream network determination that allows us to have simple stream networks that resemble the main water courses or dense stream networks that map all the permanent / temporary streams. This is an applied example of stream network determination on our online platform Hatari Utils where you can specify the number of upstream cells and interactively review the result.

June 04, 2024

Watershed delineation for basins containing internal lakes / reservoirs with Hatari Utils - Tutorial

June 04, 2024/ Saul Montoya

Hydrological basins could have a diversity of shapes due to the morphology and it´s a challenge to have a process that can delineate basins with different characteristics in a short amount of time. Hatari Utils is an online platform for hydrogeological analysis that has a tool for basin delineation, in this tutorial we have tested the capabilities of Hatari Utils to delineate a basin that has an internal lake.

Blog