Using a MODFLOW grid, generated with GMS, to solve a transport problem with TOUGH2 in complex geological environments: The intertidal deposits of the Venetian Lagoon

The tides of the Venetian Lagoon generally vary between −0.5 and +0.7 m asl. Occasionally, they may reach a maximum of 1.5 m (acqua alta) and a minimum of −0.8 m asl (acqua bassa). Intertidal areas, called “barene,” exist all along the coast of the Lagoon. These areas are characterized by canals that concentrate the flow of water (and the deposition of sands) during the rising and waning of the tides, and that inundate and drain the vegetated areas found between canals (where organic-rich clays are deposited). Therefore, since the area is subject to subsidence, in time, sand dykes (the original canals) become juxtaposed to clayey dykes (the original vegetated areas). In addition, the sands form a continuous hydrogeologic network within the clays, very similar to that of a vascular system that effectively drains the whole “barena” deposits. In order to be effective, measures for monitoring, confining, or remediating the transport of pollutants through this kind of environment must explicitly take into account the geologic complexity. The same complexity must be included in the numerical models that support remediation efforts. At the moment, there appears to be no off-the-shelf graphical interface that is able to manage such complexity for TOUGH2. To attempt to solve this problem we have used a calibrated USGS-MODFLOW model, of the barena of “Passo a Campalto” in the Venetian Lagoon, developed with the GMS graphical interface. The model is made of 42 layers, which, apart from the first layer, are 0.5 m thick; the first layer has the thickness distribution of a dump found on top of the barena deposit at Passo a Campalto. Each layer consists of 100×60 square cells, for a total of 252,000 cells, only about half of which are active. Using a FORTRAN routine, we translate this grid, with all the hydrogeologic boundary conditions, into a TOUGH2 input file, and we provide the additional necessary information for running a TOUGH2 simulation. The results are promising, in that we were able to produce TOUGH2 grids with very complex geology and to run the models with success. For visualization, the results can be imported back into GMS as 3D scatter point sets, or they can be plotted with any adequate plotting software such as MatLab. Developing conceptual and numerical models with an elaborate graphical interface such as GMS effectively allows setting up complex problems while concentrating on their physics.