How to create an Elevation Raster from Contour Lines with Python, Geopandas, Numpy and Gdal - Tutorial

/

Spatial analysis is such an interesting discipline because it allows the evaluation of every phenomena related to their location. However, for some parts of the data processing the workflow on a GIS Graphical Computer Interface (GUI) can be repetitive and time consuming. Researchers need better and more efficient tools to process more amount of data in less amount of time and even with less quantity of software tools.

We have create a innovative script to generate an elevation raster file from a contour line with several steps of data processing. The script recognizes invalid geometries, simplify the polylines and extract vertices while creates a point geodataframe that is interpolated and geotransformed as a geospatial raster in .tiff format.

Read More

How to define cell / layer confinement in MODFLOW 6 with Model Muse

/

The process and options to define the cell / layer confinement has changed in MODFLOW 6 with respect to MODFLOW 2005. We have done an explanatory video about the process to setup the confinement for cells in MODFLOW 6 with Model Muse. There is also a discussion on the influence of the confined / convertible option on the aquifer response on an example model.

Read More

Produced Water Chemistry and Isotope Analysis with Python and Pandas - Tutorial

/

On hydrocarbon extraction, water is a byproduct that comes from the geologic formations that have also oil and gas. Produced water can be defined as fossil water, or water with a high residence time on the earth, with a particular water chemistry that is important to analyze when we deal with the treatment methods, disposal techniques or safety to drinking water sources. The USGS has developed a produced water geochemical database with more than 100K points that includes spatial information, well description, rock properties, water chemistry and isotopes.

Read More

Basic Example of Stream Flow Routing (SFR) Package with MODFLOW and Model Muse

/

The disconnexion of surface groundwater flow it is fact very common on highly exploited aquifers or under impacts from climate change. Many rivers are disconnected from the groundwater flow regime decreasing their flow amount due to infiltration losses. For the simulation of this surface flow, unsaturated flow and groundwater flow phenomena the MODFLOW package Streamflow Routing (SFR) was developed.

This tutorial develops an numerical model of disconnected surface groundwater flow with MODFLOW SFR and Model Muse. The study case has regional flow and a river network with monthly flow on transient conditions over 2 years.

Read More

Geochemical Speciation Modeling and Analysis with Phreeqc and Julia - Tutorial

/

The speciation modeling allows to calculate the distribution of aqueous species in a solution. Phreeqc is capable to simulate this speciation calculation and we are going to demonstrate this capability on a study case of aqueous species in seawater. Modeling was done with a Phreeqc executable and results were analysed on a interactive Julia script. Both parts were done on a Jupyter Lab session.

Read More

Websites where you can download complete groundwater flow models in MODFLOW

/

Learning groundwater modeling is a challenge for the complexities of the codes, the lack or scarce piezometric and water quality data, the heterogeneity of the porous media and the few resources to explore complete and applied groundwater models with different boundary conditions and predictive scenarios.

Just few websites had uploaded complete groundwater flow models in MODFLOW and are accessible to the general public. We have to admit that most of the websites come from public institutions in the US. We hope that in a near future we can have more resources to get model reports and the complete MODFLOW model.

Read More

Introducing Hataripy, the unofficial Flopy version that creates awesome VTKs

/

A new unofficial version of the Flopy is available as Hataripy. This “fork” has tools for the representation of many features, boundary conditions and output with options for grid smoothing. Documentation for this library is on progress as we expect to introduce other features and identify some bugs.

Currently the library can:

  • Create the model mesh on the active zone with point and cell values

  • Create a shell of the water table

  • Create the geometry of any boundary condition

Read More

Unofficial MODFLOW flow and transport executables for Linux (Ubuntu tested)

/

It might be that you are a talented groundwater modeler and very proficient with Python and Flopy that you start to realize that you don´t need to work in Windows anymore and want to migrate your favorite groundwater modeling code to other nice, soft, and warm Linux-based operating systems. Or, it might be the case that you want to configure and run MODFLOW models as a part of a webapp or in a web server. For the reasons mentioned before and other reasons that modern world could bring to us, a compiled version of MODFLOW in Linux is needed.

Read More

How to fill missing Elevations to empty Contour Lines with PyQGIS

/

Elevation contour lines without the elevation attribute is common when we import contour lines from Autocad DXF files, but it has also happened that contour data was stored on a hard disk drive without elevation attribute and years later data is found and there is no one to ask to restore the missing information. This tutorial shows a practical procedure to fill missing elevations on contour lines with the use of PyQGIS on QGIS 3. The procedure uses an intersection line that crosses the contour lines where the base elevation and interval is known. There are some specific instructions to run the script that are well described on the video.

Read More

Temporal discretization impact on MODFLOW Model Heads and Water Budget

/

For a transient model, a modeler needs to conceptualize a coherent and feasible time discretization approach, nevertheless searching for temporal discretization guidelines the Internet gives you digital noise: nothing relevant, nothing precise, nothing useful. These is the type of situation where you ask yourself: how many stress periods do I need, how many time steps do I require, which should be an appropriate size of the model output, or more important, will the size and number of stress periods and time steps have an effect on the model output?

Read More

How to Convert a PDF to ESRI Shapefile with Python, Geopandas and Inkscape - Tutorial

/

In order to use the spatial data provided on a report we need procedures to extract the data on effective way. The amount of tools and techniques are quite advanced, and requires several open source software for specific procedures. We have done a complete tutorial with all the step required to extract the vector spatial data of a map reported as PDF into a ESRI shapefile. For this tutorial we have used Inkscape for the conversion of the PDF to DXF, QGIS to extract some information of the DXF, Python and Geopandas on a Jupyter Lab session for spatial translation and scaling.

Read More

The dark side of Python for Hydrology / Hydrogeology and the nightmare of unmet dependencies

/

The near future of data processing for Hydrology / Hydrogeology is written in Python 3, and many universities and institutions are shifting from teaching C++, Matlab or Fortran to Python. It is surprising the amount of tools, packages, codes and Ipython notebooks available for the data processing and data analysis of water related data, even with high success in data analysis reproducibility.

However, there is not always a sunshine on working with Python, especially if you are working on Windows. After you have done your first steps in hydrological series plotting, statistical analysis or even run a neural network, you might want to do more with Python and there is when the situation can become unbearable.

Read More

How to do more Groundwater Modeling (or any modeling) in less amount of time

/

Because time matters and everybody wants results and solutions now, groundwater modelers (or any water resources modeler) have to think how they can do more with the limited time available. This is a restructuration, a reengineering, a new conceptualization of not only what we do, but how we do, and if in the end we are happy with the main results, and secondary results we get from our modeling work.

Groundwater modelers are water resources specialists with high skills in computing, coding, maths and groundwater flow regime comprehension. The type of work, the complexity of the given areas of study, the limited observation data, the limited budget, the short completion times and the overall low perception and low understanding of groundwater flow and groundwater quality makes the modeling work somehow particular where we have to adopt special techniques to deliver good work while preserving the modeler integrity and wellness.

Read More

The soft and hard courses every Applied Hydrogeology Master Program should have

/

There isn’t as well a clear definition of applied hydrogeology, that following the literal meaning should be the hydrogeology in practice. We like the definition of applied hydrogeology as the real world hydrogeology, as the hydrogeology of non-homogeneous, non-infinite, non-isotropic and even non-porous media; applied hydrogeology is the hydrogeology that goes behind the book, far away from the “aquifer, aquitard, aquiclude” classification. Applied hydrogeology is the part of hydrogeology that solve problems, the science that tells us if the groundwater remediation will work, if a industrial complex will have water on the coming decade, or when the salinitization from the pumping well will overcome the reverse osmosis plant capacity.

Read More

6 aquifers that every hydrogeologist should visit on a lifetime

/

Perception and interest on the groundwater resources are of high initial interest but long term low popularity, that is why normal and plain hydrogeologist need to promote an innovative perspective of the relation from human society with the groundwater resources as a strong part for every environmental protection plan and sustainable management policy elaboration.

This article has a list of aquifers around the globe, with special characteristics where a normal audience can see how water can exist below surface, how people have interacted with the aquifer

Read More

How to work with Transmissivities in MODFLOW with Model Muse - Tutorial

/

Hydraulic parameters for an aquifer flow can be expressed on different terms like Hydraulic Conductivity (K) or Transmissivity (T). On the default setup of MODFLOW 2005 with Model Muse there is no option to use directly Transmissivity unless a manual conversion is done to transform into K, this is because the selected flow package is Layer Property Flow (LPF).

When the flow package changes to Block Centered Flow (BCF6) it is possible to insert Transmissivities directly as datasets or as a object property. The BCF package also has more options for the layer type, with unconfined and partial confined aquifer options.

Read More

Performance analysis of Open Drone Map to process hundreds of drone images

/

So far, OpenDroneMap is a strong open source alternative for the processing of drone imagery. We have done a recursive analysis over a set of aerial photos on a standard desktop computer to have a panorama of the rates and trends from OpenDroneMap in the computation of different amount of images. Although the rates and trends apply only for a specific dataset and computer configuration, this study can to be taken as a reference for the computational time and as a perspective of how the software works with great amount of data.

Read More

Tutorial of Unstructured Grid Groundwater Model with MODFLOW 6 and Model Muse 4

/

Of the latest devopments in groundwater modeling there are two softwares: Modflow 6 and Model Muse 4. both developed by the USGS. The first software is the latest version of MODFLOW that allows triangular and unstructured grids, and the second is the latest version (from June 23) of the graphical user interface Model Muse that supports Modflow 6.

Unstructured grid it a type of discretization that allows us to have small cells at certain parts of the model while the rest of the model has bigger cells. This optimization of the model grid and cell number decrease the computing time, the size of the output files and the speed of the visualization tools. With unstructured grids we can model a great extension while preserving the right accuracy of the points of interest, we can even insert regional faulting or complex geological setups.

This tutorial shows the complete procedure to create a geospatial model of a alluvial aquifer with the interaction of regional flow, river and wells. The tutorial creates the unstructured grid, boundary conditions, model geometry, simulate flow with MODFLOW 6 and represent results in Model Muse 4.

Read More

Basic Example of Land Subsidence Modeling with MODFLOW SUB and Model Muse 4

/

Cities under high exploitation of groundwater resources face severe problems of land subsidence due to pumping. This complex problem was addressed from earlier versions of MODFLOW as MODFLOW 2000 and has two packages: MODFLOW SUB and  for its simulation with different options for the conceptualization and numerical modeling of normal aquifers and low conductivity interbeds.

Scientific research and professional consulting on the field of land subsidence modeling is limited; this scarcity might be because the coupled complexity of the appliance of numerical groundwater modeling and consolidation theory.

Read More

How to convert an ArcGIS Style to QGIS 3 - Tutorial

/

Spatial data created / processed by commercial or open source software follows standards by institutions like the Open Geospatial Consortium; this standards allow the interoperability of the vector or raster data among different software however standards apply to the position but not the style. Styles from ArcGIS were not easy to convert to formats compatible with QGIS, specially if you don’t have the commercial software.

Governmental offices release spatial data as land use, land cover, infraestructures, and sometimes release styles in ArcGIS formats posing a great obstacle for the QGIS users. This tutorial shows the complete procedure to convert ArcGIS Style to QGIS as *.xml format with a case study of land cover from Costa Rica. The tutorial is developed in Windows, if you are a Linux and Mac users its necessary to install Mdbtools on your own operating system.

Read More