A minimum working MODFLOW simulation

Today one of my most important tasks is to establish a MODFLOW simulation that can work.  Considering known figures, the simulation has to be very simple.  I may only use BASIC module.

Now I am working with DISCRETION file, which contains information on the size and shape of the cell, as well as layer properties.  The main problem so far is the DEM.  The original file is too dense (with 0.1m resolution), thus I first extracted a part of data, with 100m resolution.  The result is a data file much lighter.  Though the rationale of such resolution choice still remain debetable.

From DEM I can get surface elevation needed by the model.  However, the elevations of layers underneath are also needed.  Thus I need to know the thickness of accotelm and catotelm.  Different from my previous guess, the catotelm depth in different locations may vary quite significantly.  For this period, I will set the catotelm depth to 2m according to Tallis (1973).

Mineral soil depth also varies.  I used Reeve et al.'s (2000) simplication to set it to 0.5m.

I found a useful tool for parameterising the model: PMWIN 5.3.  It provides a visual means to parameterise MODFLOW.  As a simple model I just inputed the elevation of each layer.  The hydroconductivities are left to default value.  Gamble et al. (2003) gives transmissivity of 10 ft^2/d or 0.2 ft^2/d, using different methods.  I choose 10 ft^2/d in this case.  Peat porosity also seems to vary.  I choose 0.85, which may need verification during fieldwork.