Jump to content
GMS, SMS, and WMS User Forum


  • Posts

  • Joined

  • Last visited

Everything posted by Jonas

  1. I use Kriging and variograms to describe spatial variability of pilot points for HK and RCH in a PEST parameter estimation simulation. I use the log values for both HK and RCH. For one pilot point set it works perfectly fine to use a logarithmic semivariogram, which also looks reasonable. These HK values are in the range between 1e-4 to 1e-6 m/sek. For pilot point sets with lower values (in the range between 1e-8 and 1e-10) I get the following message when I try to build a logarithmic semivariogram: Warning. There were no points calculated; bad lag parameters! A semivariogram (that is not logarithmic) does however allow me to build a variogram. I have two questions: 1. What is the reason for this warning? 2. Would it be a significant difference in interpolation between a semivariogram and a logarithmic semivariogram? I am happy to receive any comments and ideas on these issues. Kind regards, Jonas
  2. Jonas

    PEST and Kriging

    Thanks Alan, the nightly build solved the issue.
  3. Hi, When i customize the variogram for the pilot points of the HK-field in the nsmc I and II tutorials, and also in my own models, I encounter the problem showed in the attached figure. It seems like the whole field is generated from only one point and not all available ones. I also edited the already set up variogram in the tutorials, and then changed back to its original form and still got the same problem. Any ideas on how to solve this issue? Kind regards, Jonas
  4. For others interested in this issue, the support team has now fixed the bug that caused the issue of PEST and null space MC with the NWT solver. Download the latest nightly build (Jan 25th) from: http://www.aquaveo.com/downloads-gms Woodward, it was a pleasure to read your article, "Uncertainty in the modelling of spatial and temporal patterns of shallow groundwater flow paths: The role of geological and hydrological site information" in Journal of Hydrology where you have applied this method.
  5. Thank you Alan, I have done so. Looking forward for their response. Woodward, the more simplified converges in almost all cases. The problem is that non-unique grids are created.
  6. Now I have simplified the model a lot to only one layer and use the LPF package. It seems that the only solver that creates unique grid files is 2000. For 2005 and NWT it is the same issue as previous specified. Others with similar experience? Kind regards, Jonas
  7. I use the NWT solver with UPW in my model. When I changed to the same solver in the NSMC I tutorial, the same thing happens. Different values in the bpa files but non-unique grid files with HK values between the different solutions. Ideas on this issue are appreciated!
  8. Hi again! I attempt to use Null Space Monte Carlo in a similar manner as presented in PEST - Null Space Monte Carlo I tutorial (http://gmstutorials-10.1.aquaveo.com/MODFLOW-StochasticModeling-NullSpaceMonteCarloI.pdf) As presented in the attached text file I have some model that converge although most models does not converge for some reason (Ideas on how to adress this are also appreciated). As seen in the text file, there are different calculated values for the pilot points between the four converged solutions. The bpa files for the different solutions are also different. All good so far. My issue is that there is no difference between the grid files of the HK parameter for the four different solutions. Any ideas on how to solve this issue? Could it have something to do with the Kriging process? Best regards, Jonas T40_nsmc.txt
  9. I created the zones based on solids. For some reason the materials got the right values according to the zones but some cells got wrong values of HK. I changed this manually by selecting these cells. This worked. Could this be a bug in GMS maybe?
  10. Hi! I have set up a model for a municipal water supply in an glacifluvial esker. I want to calibrate the hydraulic conductivity (HK) against steady-state heads in the esker material with PEST and pilot points. I also have other materials in the model but with a fixed HK. The recharge is also calibrated with PEST in three different polygons. The picture below shows layer four (of five). Left shows materials (only bedrock and esker in this layer) together with locations of the pilot points. The middle part shows calibrated heads after PEST has been running. Right shows HK after calibration with PEST. My issue is that the area with esker material (left) does not overlap the area that is calibrated with pilot points in the rightmost picture (non red area). Something that looks strange in particular is the straight paths from the middle to the south of the model that get the same HK as the bedrock although the material is set to esker. It is the same issue in the other layers of the model. Ideas on how to solve this issue are appreciated. Kind regards, Jonas
  • Create New...