Jump to content
GMS, SMS, and WMS User Forum

One layer to multilayer

Chaachoua Mohamed

Recommended Posts


I have an hydrodynamic model elaborated with GMS 6.0, whose features are as follows:

1. The aquifer is unconfined

2. The model consists of a single layer

3. the model is simulated in transient of 35 years.

When I subdivided the model into several layers. I kept the first layer unconfined and the lowers confined (i think i should do so). For lowers layers I subdivised the specific yield by the the thickness of each layer to obtain the specific storage). The pumping is applicated through all the layers pro rata their thickness, I have considered the kz = Kh. But I do not get the same results of the original model.

How I could choose the parameters of my model multilayer among the original model to obtain the same results (the specific yield should become the specific storage for lower layers) and how to choose the vertical permeability of the layers?

I am very happy for any help


Link to comment
Share on other sites

It is not possible to consider the lower layers unconfined. When you subdivise the original layer to more than one layer you get for each layer a new top. This top is necessary above the head. So, each layer will be considered confined wich is more logical.

Thanks for your response and I am always waiting for your help

Link to comment
Share on other sites

Hello Sean,

Thanks for your help. Of course I have considered the lower layers as convertible. I work with GMS 6.0 (lpf mode) so i have only two options : unconfined or convertible. The problem is can I choice parametres (S, Kz) from the original model and how can I reaffect pumping data through all layers so that I obtain the same results like my original model. I have explained in my first message how I have done it but unfortunately I could not obtain the same results.

Thanks for more help

Link to comment
Share on other sites

The simpler parameter is Kz - In your single layer model, there is no restriction to water moving vertically. To get similar results when you split the layer, you need to make Kz = Kh for all layers. Otherwise, the preferred flow path is horizontal rather than vertical. Other people might be able to provide more information on specific storage/storativity, but I've found model results to be less sensitive to that parameter. I suggest selecting a ballpark value from literature based on your soil, then try a couple of sensitivity runs to see how much it impacts the results.

Good luck!

Link to comment
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.

Reply to this topic...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.


  • Create New...