Copyright © 2019 EarthSoft, Inc • Modified: 31 May 2019
Particle tracking enables users to visualize the movement of conservative particles in groundwater that is consistent with a set of water elevation measurements. Particle tracks are calculated by estimating the water velocity within each cell of a rectangular grid and advecting the particle by a distance equivalent to the product of the velocity magnitude and the user-selected time increment. The new particle location is determined and the particle repeatedly moved for a finite number of steps specified by the users. EnviroInsite interpolates the measured values to points on the rectangular grid. The same interpolation routines are used for particle tracking as are used elsewhere in EnviroInsite for interpolating data prior to contouring.
Prior to generating particle tracks, develop an appropriate potentiometric surface. Use the same grid and interpolation method /parameters for both the groundwater elevation contours and the particle tracking.
View a training video on particle tracking here.
Click Plot> 2D Data from the main menu and select Particle Track. The Particle Tracking dialog box opens. Modify the particle tracking properties on the Query tab, Format tab, Options tab, Grid tab, Interpolation tab, Virtual Points tab, and EQuIS Query tab as desired. Click the OK button to save changes.
This tab allows the user to select the interpolation scheme and the parameters of the interpolation method. The correct selection of interpolation parameters is critical for generation of contours that accurately reflect the field data and our expectations of how the values vary between the measured data points. The default parameters are frequently adequate, although some improvement can be anticipated through trial and error.
Method: Select the contouring method from inverse distance, kriging, and natural neighbor.
Search Limit Radius: The cutoff distance establishes the maximum distance from the interpolated point for which a measured value will be used in the interpolation function. Used by both kriging and inverse distance routines.
The kriging routines are derived from the kt3d routine of the Geostatistical Software Library (GSLIB) authored by Clayton Deutsch and Andre Journel (www.gslib.com).
Variogram Model: Select from the spherical, exponential, gaussian, power, and hole variogram models using conventions of GSLIB.
Sill: Set the variogram model sill or c value.
Nugget: Set local scale variability.
Range a max: set maximum horizontal range.
Range a min: Set minimum horizontal range.
Range a vert: Set vertical range.
Angle: Rotation angle in degrees clockwise about 12 o'clock of axis defining the direction of maximum horizontal range. In some cases it may be advantageous to increase the a max value and to align this axis with the direction of ground-water flow.
Inverse Distance Parameters
Z Scale: The vertical distances are multiplied by this value prior to calculating the interpolation weight assigned to a measured value. Z Scale values greater than one will result in lower weights assigned to measured values that are not at the same elevation as the grid point.
Exponent: The measured values are assigned weights that are equal to one over the distance between the measured value and the grid point raised to some power. This exponent is the power used in that calculation.
Smooth Distance: The inverse distance method tends to cause a kind of bubbly surface with the bubbles coinciding with measurement points. This type of trend can be reduced by introducing a non-zero smoothing distance. The smoothing distance is effected by adding this distance to the separation distance between the measurement point and the grid point prior to calculating the interpolation weight.
Virtual Points Tab
Virtual points are used to control generating contours with sparse data. In those cases, it may be advantageous to create virtual measurement points that will control the resulting contours.